ABAQUS用户子程序学习小结.doc
4页
ABAQUS用户子程序学习小结[折叠] 1 FORTRAN语言中的“I-N 规则”：I、J、K、L、M、N 开头的为整型变量，其他开头为实型变量；2 DIMENSION COORDS(3)表示声明一个含 3个元素的数组，下标分别为 1、2、3，访问形式为 COORDS(n)，n 为 1～3；3 子程序（*.for）文件中如何输出调试信息：WRITE(6,*)'COORDS(1)',COORDS(1)，在*.dat文件中可看到输出，如果希望 WRITE输出到 msg文件中，则写为 WRITE(7,*)'COORDS...；4 用户子程序 DLOAD中 COORDS数组的含义：COORDS(1)也是一个数组，存贮单元集合中所有单元积分点的 X坐标，COORDS(2)存贮 Y坐标，相应 INP文件中的写法为：*DLOADPY,PYNU其中 PY为单元集合名称，定义方法为：*Elset, elset=BEAM, generate1, 5, 1...*ELSET,ELSET=PYBEAM5 DLOAD中 F的定义方法：F只有定义在单元积分点上才有效，例如：F=1.0*COORDS (1)附一个简单实例：beam.inp文件：*Heading** Job name: Job-1 Model name: beam*Preprint, echo=NO, model=NO, history=NO, contact=NO**** PARTS***Part, name=PART-1*End Part**** ASSEMBLY***Assembly, name=Assembly** *Instance, name=PART-1-1, part=PART-1*Node1, 0., 0.2, 20., 0.3, 40., 0. 4, 60., 0.5, 80., 0.6, 100., 0.*Element, type=B311, 1, 22, 2, 33, 3, 44, 4, 55, 5, 6*Elset, elset=BEAM, generate1, 5, 1** Region: (Section-1-BEAM:BEAM), (Beam Orientation:BEAM)** Section: Section-1-BEAM Profile: Profile-1*Beam Section, elset=BEAM, material=STEEL, temperature=GRADIENTS, section=RECT0.2, 5.0.,0.,-1.*End Instance*Nset, nset=ENDS, instance=PART-1-11, 6*Nset, nset=_M4, internal, instance=PART-1-16,*Nset, nset=_M5, internal, instance=PART-1-11,*End Assembly** ** MATERIALS** *Material, name=STEEL*Elastic210000., 0.3*ELSET,ELSET=PYBEAM** ** BOUNDARY CONDITIONS** ** Name: Disp-BC-1 Type: Symmetry/Antisymmetry/Encastre*Boundary_M4, ENCASTRE** ----------------------------------------------------------------** ** STEP: Step-1** *Step, name=Step-1*Static** ** LOADS** ** Name: CFORCE-1 Type: Concentrated force*DLOADPY,PYNU** ** OUTPUT REQUESTS** ** ** FIELD OUTPUT: F-Output-1** *Output, field, variable=PRESELECT** ** FIELD OUTPUT: F-Output-2** *Output, field*Element OutputSF, ** ** HISTORY OUTPUT: H-Output-1** *Output, history*Node Output, nset=ENDSCF1, CF2, CF3, CM1, CM2, CM3, RF1, RF2RF3, RM1, RM2, RM3, U1, U2, U3, UR1UR2, UR3*El Print, freq=999999*Node Print, freq=999999*End Stepbbb.for文件SUBROUTINE DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,1 JLTYP,SNAME)CINCLUDE 'ABA_PARAM.INC'CDIMENSION TIME(2), COORDS (3)CHARACTER*80 SNAMEWRITE(6,*)'COORDS(3)',COORDS(3)F=1.0*COORDS (1)RETURN END运行方法：在 Abaqus Command提示符后输入：abaqus job=beam user=bbb interactive。
