ABAQUS用戶子程序?qū)W習(xí)小結(jié)

1、1.1 FORTRAN語(yǔ)言不區(qū)分大小寫,文件擴(kuò)展名為for,正文從第7列開始,第15列是標(biāo)號(hào)區(qū),第六列是續(xù)行標(biāo)志區(qū)。1.2 FORTRAN語(yǔ)言中的“I-N規(guī)則”:I、J、K、L、M、N開頭的為整型變量,其他開頭為實(shí)型變量;2 DIMENSION COORDS(3表示聲明一個(gè)含3個(gè)元素的數(shù)組,下標(biāo)分別為1、2、3,訪問(wèn)形式為COORDS(n,n為13;3 子程序(*.for文件中如何輸出調(diào)試信息:WRITE(6,*'COORDS(1',COORDS(1,在*.dat文件中可看到輸出,如果希望WRITE輸出到msg文件中,則寫為WRITE(7,*'COORDS.;4 用戶子

2、程序DLOAD中COORDS數(shù)組的含義:COORDS(1也是一個(gè)數(shù)組,存貯單元集合中所有單元積分點(diǎn)的X坐標(biāo),COORDS(2存貯Y坐標(biāo),相應(yīng)INP文件中的寫法為:*DLOADPY,PYNU其中PY為單元集合名稱,定義方法為:*Elset, elset=BEAM, generate1, 5, 1.*ELSET,ELSET=PYBEAM5 DLOAD中F的定義方法:F只有定義在單元積分點(diǎn)上才有效,例如:F=1.0*COORDS (1附一個(gè)簡(jiǎn)單實(shí)例:beam.inp文件:*Heading* Job name: Job-1 Model name: beam*Preprint, echo=NO, mod

3、el=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: (Se

4、ction-1-BEAM:BEAM, (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, inter

5、nal, 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* Name: CFORCE-1 Type: Concentrated force*DLOADPY,PYN

6、U* 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,SNAMECINCLUDE 'ABA_PARAM.INC'CDIMENSION TIME(2, COORDS (3CHARACTER*80 SNAMEWRITE(6,*'