《ADINA例题手册-命令流文件.docx》由会员分享，可在线阅读，更多相关《ADINA例题手册-命令流文件.docx（166页珍藏版）》请在课桌文档上搜索。
1、*本例提供一种利用rigidlink实现beam与solid连接的方法*FeprogramPROGRAM=ADlNA*HEADING,Examplem1:Beamconnectwithasolid*建模*BODYB1.OCKNAME=IOPTloN=CENTEREDPOS1TION=VECTORORIENTAT=SYSTEMCXl=OCX2=0CX3=0SYSTEM=ODXl=IDX2=1DX3=1*COORDINATESPOINTSYSTEM=O9000.50100020*1.INESTRAIGHTNAME=IPl=9P2=10*一定义材料及单元组*MATERIA1.E1.ASTICNAME
2、=IE=200000NU=O.3DENSTTY=7800*egroupThreedsolidnme=imateril=i*Cross-SECTIOpipename=idiameter=o.2thicknes=o.05sc=otc=otorfac=iSSHEARf=OTSHEArf=OSO1.ID=NO*EGROUPBEAMNAME=2SUBTYPE=THREE-DMATERIA1.=1SECTION=1*定义单元尺寸及划分网格*SUBDIVIDEBODYNAME=IMODE=1.ENGTHSIZE=O.21*SUBDIVIDE1.INENAME=IMODE=1.ENGTHSIZE=O.21GB
3、ODYN0DES=27GROUP=I1O*G1.INEN0DES=2AUXP0INT=5GR0UP=21*定义刚性连接*RIGID1.INK-NONAme=ISlavenod=OMASTERNO=ODlSP1.ACE=DEFAu1.T2661096DEFAU1.T12701096DEFAU1.T2721096,DEFAU1.T,2761096DEFAU1.T,*定义约束及荷载*FIXITYNAME=ZF,Z-TRNS1.TION,*FIXBOUNDARYFACESFIXITY=A1.1.BODY=I5ZF*1.OADMOMENTNAME=IMAGNITUDESMX=OMY=OMZ=I*APP1
4、.Y-1.OADBODY=I1MOMENT1,POINT,10O1OO-1OOO,N0,OOlO*设置动力分析及大位移*masterAnalysis=Dynamic-Direct-Integration*ANA1.YSISDYNAMIC-DIRECT-INTEGRATIONMETHOD=NEWMARKDE1.TA=O.5A1.PHA=O.25*KINEMATICSDISP1.CE=1.RGE*TIMESTEPNAME=DEFAU1.TC1.EAR501* 这是一个shell与solid连接的小例子，利用meshingcreatemesh-shelltransition来完成。* 可以通过dis
5、playgeometrymeshpolt-modify,点击elementdepiction,再在弹出窗口中把SheIlelementattributes由mid-SUrfaCe改为TOP/Bottom。*FEPROGRAMPROGRAM=ADlNA*建立body并对其分网*BODYB1.OCKNAME=IOPTION=CENTEREDPOSITION=VECTORORIENTAT=SYSTEMCXl=OCX2=0CX3=0SYSTEM=ODXl=IDX2=1DX3=0.2*SUBDIVIDEBODYNAME=IMODE=1.ENGTHSIZE=O.2*MATERIA1.E1.ASTICNAM
6、E=IE=2E+11NUR.3DENSITY=7800*EGROUPThreedsolidname=imaterial=i*GBODYN0DES=8NCOTNCID=BOUNDARIESGROUP=I10*-建立Shell的几何模型并指定厚度0.1*COORDINATESPOINTSYSTEM=O90.50.500100.5-0.500*1.INESTRAIGHTNAME=IPl=9P2=10*SURFACEEXTRUDEDNAME=I1.INE=IDX=2DY=ODZ=OSYSTEM=OPCOINCID=YESptoleran=ie-05ndiv=ioption=vector*Sfthick
7、nesso.*一一对Shell单元组指定厚度0.2（Shell单元厚度以GeOnIetry中定义的为准）*SUBDIVIDESURFACENAME=IMODE=DIVISIONSNDIV1=5NDIV2=5*GSURFACEN0DES=4NColNCID=BOUNDARlES1* 用Shel1transition功能连接Shel1和SOIid单元，自动生成过渡Shell单元*Transttion-Ebodyi=ofacei=iedge=igroup1=2body2=iface2=6group2=i* 施加边界条件及荷载*FIXBOUNDARY1.INESFIX1TY=1.1.4,A1.1.,*
8、1.OADMASS-PROPORTIONA1.NAME=IMGN1TUD=9.8AX=OAY=OAZ=TINTERPRE=BODY-FORCE*APP1.Y-1.OADBODY=I1MASS-PROPORTIONA1.1,MODE1.,00100-1000,N0,OOlO*Ayia6eIAA1113I-OIoEpAyEO2aI1.OODAIo,IA6Aa1EAAUEYEAaA1.-.0u6uAyia,E66A9OOUaiNSoaeEDDiFeprogramPRoGRAM=AD工NA*A1.*COORDINATESPOINTSYSTEM=Oentriesnameyz1 -80-402 -8060
9、3 80-404 80605 -10-106 -10107 10-108 10109 -9-910 -9911 9-912 9913 00SURFACEVERTEXNAME=IPl=4P2=2P3=6P4=8SURFACEVERTEXNAME=2Pl=6P2=2P3=lP4=5SURFACEVERTEXNAME=3Pl=7P2=5P3=lP4=3SURFACE*VERTEXNAME=5Pl=8P2=6P3=10P4=12SURFACE*VERTEXNAME=6Pl=IOP2=6P3=5P4=9SURFACE*VERTEXNAME=7Pl=IlP2=9P3=5P4=7SURFACE*VERTEX
10、NAME=8Pl=8P2=12P3=llP4=7SURFACEVERTEXNAME=9Pl=13P2=8P3=6P4=13P2=8P3=7P4=3SURFACEVERTEXNAME=4Pl=4SURFACEVERTEXNAME=IOPl=13P2=6P3=5P4=13SURFACEVERTEXNAME=IlPl=13P2=5P3=7P4=13SURFACEVERTEXNAME=12Pl=8P2=13P3=7P4=8*afi0o060*FIXITYNAME=YF,y-translation,0FIXITYNAME=ZF,z-translation,FIXBOUNDARY1.INESFIXITY=
11、ZF5,YF,9,ZF,12,YF,*1.OADMASS-PROPORTIONA1.NAME=IMAGNITUD=9.8AX=OAY=OAZ=-IInterpre=Body-ForceAPP1.Y-1.OADBODY=O1,MASS-PROPORTIONA1.,1MODE1.00100-1000,NO,OOlOSTRAIN-FIE1.DNAME=IA=-1630000B=27000C=O.3AXESCONSTANTNAME=IAX=IAY=OAZ=OBX=OBY=IBZ=OSET-AXES-STRSURFACES1123to121235I02AAioYa6MATERIA1.MOHR-COU1.
12、OMBNAME=IE=6E+08NU=O.2PHI=25PSI=OCOH=8000000TCUT=2000000DENSITY=2700DI1.ATION=NOMDESCRIP=,NONE,MATERIA1.E1.ASTICNAME=2E=1E+1ONU=O.2DENSITY=2200A1.PHA=Omdescrip=none,*egroupTwodsolidname=iSubtype=Strainmaterial=iInitials=Elementstrain-f=iegroupTwodsolidname=2subtype=strainmaterial=2*SIYa3,oOOl0jftSUB
13、DIVIDE1.INENAME=IMODE=DIVISIONSNDIV=IO1to12/14/17/18/20SUBDIVIDE1.INENAME=21MODE=DIVISIONSNDIV=821to24*Gsurfacenodes=4NC。工NC工D=BoUNDAR工ESgroup=i1to4/9to12Gsurfacenodes=4NCOlNC工D=BoUNDArj:ESgroup=25to8*SI,aYOaEuEASurf-E1.EMdattwodsolidentriessurfacetbirthtdeath540.10660.10780.10820.109030.110050.1110
14、70.112010.1EeOAEa2o-OioAaDi-masterAnalysis=STATicinitials=defKINEMATICSDISP1.ACE=1.ARGE*TIMESTEPNAME=DEFAU1.TC1.EAR1010HEADINGSTRING=,Example4:Thesideslopewithanchoringcalculation,*Feprogramprogram=adina*XA1.DtCOORDINATESPOINTentriesnameyz1002 0103 10104 1185 1276 1357 1528 1719 20010 141011 15912 1
15、68.513 18414 22020310.9021020411.902820512.902720613.902520715.902220817.902120920.90201.INESTRAIGHTNAME=I*Pl=2P2=l1.INESTRAIGHTNAME=2Pl=IP2=9*1.INESTRAIGHTNAME=3*Pl=2P2=31.INESTRAIGHTNAME=4Pl=203P2=10*1.INESTRAIGHTNAME=5Pl=209P2=14*1.INEPO1.Y1.INENAME=6TYPE=SP1.INE3/4/5/6/7/8/9*1.INEPO1.Y1.INENAME=
16、7TYPE=SP1.INE203/204/205/206/207/208/2091.INEPO1.Y1.INENAME=8TYPE=SP1.INE10/11/12/13/14SURFACEPATCHNAME=IEDGE1=3EDGE2=1EDGE3=2EDGE4=6SURFACEPATCHNAME=2EDGE1=4EDGE2=7EDGE3=5EDGE4=8*SIQIGEpX。EO0FIXITYNAME=Y,y-translation,*FIXITYNAME=Z,z-translation,FIXBOUNDARY1.INES1,Y2,Z,5,Z,1.OADPRESSURENAME=IMAGNIT
17、UD=-1500000*APP1.Y-1.OADBODY=O1,PRESSURE11,1.INE,30100-1000,N0,OO1O2,PRESSURE11,1.INE,40100-1000,N0,OO1O*o,*CGROUPC0NTACT2NAME=ISUBTYPE=STRAINContactsurfaname=isolid=nobody=oorientat=automatic610*Contactsurfaname=2solid=nobody=oorientat=automatic710*CONTACTPAIRNAME=ITARGET=ICONTACTO=2* OApAADfaDEJiA
18、0COORDINATESPOINTentriesnameyz20310.0021020411.002820512.002720613.002520715.002220817.002120920.00200ANiA2AAioYaeMATERIA1.E1.ASTICNAME=IE=2E09NU=O.2DENSITY=1800*MATERIA1.E1.ASTICNAME=2E=1E+O8NU=O.2DENSITY=2000*egroupTwodsolidname=iSubtype=Strainmaterial=iegroupTwodsolidname=2SUBTYPE=STRAINmaterial=
19、2* l0jnz6Di2l0ifiSUBDIVIDESURFACENAME=IMODE=1.ENGTHSIZE=O.51/2Gsurfacenodes=4ncoincid=boundariesgroup=i*Gsurfacenodes=4ncoincid=nogroup=22* 0AajEACOORDINATESPOINTentriesnameVZ100l13910027810031471004861005155100694*1.INESTRAIGHTNAME=9Pl=1002P2=10011.INESTRAIGHTNAME=IOPl=1004P2=10031.INESTRAIGHTNAME=
20、IlPl=1006P2=1005*1.INE-E1.EMDATTRUSS900.01to1100.01*SIOAajEA2AAioYOaeMATERIA1.E1.ASTICNAME=3E=2E+11NU=O.3DENSITY=7800EGROUPTRUSSNAME=3SUBTYPE=GENERA1.MATERIA1.=30PTI0N=REBARRB-1.INE=IREBAR-1.INENAME=I91011*EeOA：0OiiimasterAnalysis=Staticidof=100iiiTIMEFUNCTIONNAME=IIF1.IB=IC1.EAR001001TIMESTEPNAME=D
21、EFAU1.TC1.EAR1010KINEMATICSDISP1.ACE=1.ARGE*Feprogramprogram=adina*读入模型*READSlope2d.in*-删除有限元模型及约束和荷载条件*DE1.ETE-FE-MOPRoGRAM=ADlNA*FIXBOUNDARY1.INESFIX1TY=1.1.C1.EAR*DE1.ETEPP1.Y-1.OADFIRST=I1.AST=I*DE1.ETEAPP1.Y-1.OADFIRST=21.AST=2*-修改锚杆坐标*COORDINATESPOINT10010.513910020.57810030.514710040.58610050
22、.515510060.594*拉伸成体*VO1.UMEEXTRUDEDNAME:ISURFACE=ODX=IDY=ODZ=ONDIV=21/2*定义边界条件及荷载*FIXBOUNDARYSURFACES4Y/5,Z,10Z*PP1.Y-1.OD1PRESSURE1SURFACE30100-1000NOOOlO2PRESSURE1,SURFACE*80100-1000,N0,OOlO*定义接触*CGROUPC0NTACT3NAME=I*Contactsurfaname=isolid=nobody=o61O*ContactsurfaNAME=2solid=nobody=o91O*CONTACTPA
23、IRNAME=ITARGET=IC0NTACT0=2*定义岩体的单元组并划分网格*egroupThreedsolidname=imaterial=i*egroupTiireedsolidNAME=2material=2*EGROUPTRUSSNAME=3SUBTYPE=GENERA1.MATERIA1.=3OPTlON=REBARRB-1.INE=I*GVO1.umenodes=8Ncoincid=BoundariesGROUP=IMESHING=MAPPED1*GVO1.UMEN0DES=8NcoINCID=NOGROUP=2MESHlNG=UAPPED2*士Ay06afaamEljA0*A
24、6Af1i-eAy6Ap%6pe/E2i_1QfdifFeprogramprogram=adinaalCOORDINATESPOINTentriesnameyz1234O11OOO0.30.51.INESTRAIGHTNAME=IPl=4P2=3SP1.IT-1.INENAME=IUSP1.IT=O.2SP1.IT-1.INENAME=3*USP1.IT=O.1SURFACEEXTRUDEDNAME=I1.INE=4DX=ODY=O.OlDZ=O.05SYSTEM=OPcoincid=Yesndiv=ioption=vectorSURFACEVERTEXNAME=2Pl=3P2=4P3=lP4
25、=2O2AAioYaeMATERIA1.E1.ASTICNAME=IE=2E+11NU=O.3DENSITY=7800egroupTwodsolidNAME=ISubtype=Strainmaterial=iegroupTwodsolidNAME=2Subtype=Strainmaterial=iafi0oE0FIXBOUNDARY1.INESFIXITY=A1.1.10,A1.1.,AY=OAZ=-IN0,00101.OADMASS-PROPORTIONA1.NAME=IMAGNITUD=9.8AX=OInterpre=Body-ForceAPP1.Y-1.OADBODY=O1MASS-PR
26、OPORTIONA1.1MODE1.,00100-10000A6zYCGROUPC0NTACT2NAME=ISUBTYPE=STRAINContactsurfaname=isolid=nobody=oIlOContactsurfaname=2solid=nobody=o410*CONTACTPAIRNAME=ITARGET=ICONTACTO=2FRICTION=0.Ol*,aYa3zoXg)l04nSUBDIVIDEMODE1.MODE=DIVISIONSNDIV=IOSUBDIVIDE1.INENAME=IMODE=DIVISIONSNDIV=2010*Gsurfacenodes=4NeC
27、)工NC工D=BoUNDAR工ESgroup=i2*Gsurfacenodes=4ncoincid=nogroup=2*fieAffiEafiiiimasterAnalysis=Dynamic-Direct-Integration*ANA1.YSISDYNAMIC-DIRECT-INTEGRATIONMETHOD=NEWMARKDE1.TA=O.5A1.PHA=O.25*KINEMATICSDISP1.ACE=1.ARGE*TIMESTEPNAME=DEFAU1.TC1.EAR1000.01*本例为了说明线与体的接触方法，此方法同样适用于边与体的接触。利用菜单modelcontact-cont
28、actbody定义线的接触边界。*FeprogramPROGRAM=ADINA*建模*BODYB1.OCKNAME=IOPTION=CENTEREDPOSITION=VECTORORIENTAT=SYSTEMCXl=OCX2=0CX3=0SYSTEM=ODX1=5DX2=5DX3=1*COORDINATESPOINTSYSTEM=O90030104000*1.INESTRAIGHTNAME=IPl=IOP2=9*-定义材料及单元组*MATERIA1.E1.ASTICNAME=IE=2E+11NU=O.3DENSITY=7800*MATERIA1.NON1.INEAR-E1.ASTICNAUE=
29、2DENSITY-7800-100012E+6*egroupThreedsolidname=imaterial=i*1.ine-E1.EMDATtruss100.0001*EGROUPTRUSSNAME=2MATERIA1.=2*通过ContactBody来定义线与体的接触*CGROUPCONTACTSNAME=I*CONTACTBODYNAME=ISO1.ID=NOBODY=IDD,1.INE,1*ContactsurfaNAME=2solid=yesbody=i11051O61O*CONTACTPAIRNAME=ITARGET=2CONTACTO=1* -施加边界条件及荷载*Ftxboun
30、daryfacesfixity=allbody=i5,A1.1.,*1.OADMASS-PROPORTIONA1.NAME=IMAGNITUD=9.80000000000000,Ax=o.ooooooooooooooAY=o.ooooooooooooooaz=-i.oooooooooooooo,Interpre=Body-Force*apply-loadbody=i1,MASS-PROPORTIONA1.,1,MODE1.*00100-1000,NO,OOlO* 划分网格*SUBDIVIDE1.INENAME=IMODE=DIVISIONSNDIV=20*GBODYN0DES=8NCOTNCI
31、D=BOUNDARIESGROUP=I1*G1.INEN0DES=2NCOINCID=ENDSNCENDS=12NCTO1.ERA=1E-O5GR0UP=21* 设置时间步及动力分析选项*masterAnalysis=Dynamic-Direct-Integration*analysisDynamic-Dtrect-IntegrationMETHoD=NEwMARKdelta=o.5alpha=o.25*kinemticsDisplace=1.arge*RAY1.EIGH-DAMA1.PHA=O.05BETA=O.05*TIMESTEPNAME=DEFAU1.T觥1.EAR600.005400
32、0.001*本例采用二维轴对称单元计算过廉装配问题*HEADINGSTRING=*Example7:Assemblyanalysiswithover-contact*Feprogramprogram=adina*建模*COORDINATESPOINT1000200.203006400.26500.192.5600.193.5700.32.5800.33.5*SURFACEVERTEXNAME=IPl=4P2=3P3=lP4=2*SURFACEVERTEXNAME=2Pl=8P2=6P3=5P4=7*-定义约束*FIXITYNAME=ZFz-translation,*FIXBOUNDARYPOI
33、NTS1 ,A1.1.,/3A1.1.,/7ZF,*定义材料及单元组*MATERIA1.E1.ASTICNAME=IE=2E+10NU=O.3DENSITY=7800*egroupTwodsolidname=iSubtype=Axisymmetricmaterial=i*EGRO1.TTwodsolidNAME=2Subtype=Axisymmetricmaterial=i*定义接触组,指定一个时间单位把初始穿透消除*CGROUPC0NTACT2NAME=ISUBTYPE=AXISYMMETRICTlME-PEN=I*Contactsurfaname=isolid=no4*Contactsur
34、faNAME=2solid=no6*CONTACTPAIRNAME=ITARGET=IC0NTCT0=2*一定义单元尺寸*SUBDIVIDESURFACENAME=IMODE=DIVISIONSNDIVl=IONDIV2=30*SUBDIVIDESURFACENAME=2MODE=DIVISIONSNDIV1=5NDIV2=10* 划分网格，在划分第二格网格时，要指定nochecking*Gsurfacenodes=9ncoincid=boundariesGRoUP=I1*Gsurfacenodes=9ncoincid=noGRoUP=22* 设置控制选项*masteranlysis=stat
35、ic*iterationMethod=Full-Newtonline-sea=yesmax-iter=15*TIMESTEPNAME=DEFU1.TC1.EAR200.05*Feprogramprogram=adina*读入模型*READOver-Contact2D.in*-删除有限元模型及约束和荷载条件*DE1.ETE-FE-MOPRoGRAM=ADlNA*FIXBOUNDARYPOINTSC1.EAR1 A1.1.3,A1.1.,* 旋转成体*VO1.UMEREVO1.VEDNAME=IMODE=AXISSURFACE=OANG1.E=180AXIS=Z1.NDIV=181/2* 定义边界
36、条件及荷载*FIXITYNAME=XFX-TRNS1.TION,*FIXBOUNDARY1.INES17ZF*FIXBOUNDARYSURFACES1XF2XF6,XF,11XF* 定义接触*CGROUPC0NTACT3NAME=ITTME-PEN=I*Contactsurfaname=isolid=no510*Contactsurfaname=2solid=no8o*CONTACTPAIRNAME=ITARGET=IC0NTACT0=2*-定义单元组并划分网格*egroupThreedsolidname=imaterial=i*egroupThreedsolidNAME=2mteril=i*
37、Gvolumenodes=8Ncoincid=Boundariesgroup=1meshing=mapped1*GVO1.UMEN0DES=8NcoINCID=NOGROUP=2MESHING=MAPPED2* 本例主要说明两个问题：1、各向异性材料的用法；2、体单元外贴一层壳单元的做法* 各向异性材料使用时与一般材料的主要不同之处在于需要定义材料轴* 体单元外贴一层壳单元的方法是给体划分完体单元后，再给它的某一个面划分壳单元*FEPROGRAMPROGRAM=ADlNA*建模*BODYB1.OCKNAME=IOPTION=CENTEREDPOSITION=VECTORORIENTAT=SYS
38、TEMCX1=2.5CX2=0CX3=0SYSTEM=ODX1=5DX2=0.2DX3=0.4*一定义约束及荷载*Ftxboundaryfacesbody=i3A1.1.*1.OADFORCENAME=IMAGNITUD=10000FZ=-I*APP1.Y-1.OADBODY=I1 ,FORCE,1,POINT,60100-1000,N0,00102 ,FORCE,1,POINT*50100-1000,N0,0010* 一定义材料及单元组*MATERIA1.E1.ASTICNAME=IE=1E+O9NU=O.167DENSITY=2400*MATERIA1.ORTHOTROPICNAME=2E
39、A=5E+10EB=1E+1OEC=1E+1ONUB=O.3NUAC=O.1NUBC=O.1GAB=IE+10GAC=1E+1OGBC=IE+10DENSITY=1200WRINK1.E=NOW-TIME=O*egroupThreedsolidnme=imteril=i*EGROUPSHE1.1.NAME=2MATERIA1.=2THICKNES=O.01* -定义材料轴AXESCONSTANTNAME=IAX=IAY=OAZ=OBX=OBY=IBZ=O*SET-AXES-MATFACESBODY=I1112* 一定义单元尺寸及划分网格*SUBDIVIDEBODYNAME=IMODE=1.EN
40、GTHSIZE=O.1*SUBDIVIDEEDGENAME=IOBODY=IMODE=1.ENGTHSIZE=O.32/4/12*GBODYNoDES=8GROUP=IMESHTNG=MAPPED1*GFACEN0DES=4NCOINCID=A1.1.GR0UP=21*Feprogramprogram=adina* EeOA-OIoAaDiHEADINGSTRING=,Example9:ADINA-MConstructPressureVesseland1.inearAnalysis,*masterAnalysis=Staticidof=iii* A1.-BODYSPHERENAME=IDIMENSIO=RADIUSCXl=OCX2=0CX3=0RADIUS=ITransformatiscalename=imode=systemsx=isy=iSZ=Cl.4BODYTRANSFORMEDNAME=IOPTION=MOVETRANSFOR=I*BODY
