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, generate 1, 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*Node 1, 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, generate 1, 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-1 1, 6*Nset, nset=_M4, internal, instance=PART-1-1 6,*Nset, nset=_M5, internal, instance=PART-1-1 1,*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)C INCLUDE 'ABA_PARAM.INC'C DIMENSION TIME(2), COORDS (3) CHARACTER*80 SNAME WRITE(6,*)'COORDS(3)',COORDS(3) F=1.0*COORDS (1) RETURN END运行方法：在Abaqus Command提示符后输入：abaqus job=beam user=bbb interactive。