ABAQUS用户子程序学习小结

1、1.1 FORTRAN语言不区分大小写,文件扩展名为for,正文从第7列开始,第15列是标号区,第六列是续行标志区。1.2 FORTRAN语言中的“I-N规则”:I、J、K、L、M、N开头的为整型变量,其他开头为实型变量;2 DIMENSION COORDS(3表示声明一个含3个元素的数组,下标分别为1、2、3,访问形式为COORDS(n,n为13;3 子程序(*.for文件中如何输出调试信息:WRITE(6,*COORDS(1,COORDS(1,在*.dat文件中可看到输出,如果希望WRITE输出到msg文件中,则写为WRITE(7,*COORDS.;4 用户子程序DLOAD中COORDS数

2、组的含义: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=

3、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:BE

4、AM, (Beam Orientation:BEAM* Section: Section-1-BEAM Profile: Profile-1*Beam Section, elset=BEAM, material=STEEL, temperature=GRADIENTS, section=RECT 0.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=P

5、ART-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* Name: CFORCE-1 Type: Concentrated force*DLOADPY,PYNU* OUTPUT REQUE

6、STS* 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,SNAMECINCLUDE ABA_PARAM.INCCDIMENSION TIME(2, COORDS (3CHARACTER*80 SNAMEWRITE