finish
/clear
/prep7
LC1=-7.5
K,1,-0.250,1.125,0 !K,ID,x,y,z
K,2,0,1.125,0
K,3,0.250,1.125,0
K,4,-0.325,-1.125,0
K,5,0,-1.125,0
K,6,0.325,-1.125,0
K,7,,1.125,LC1
L,1,2
L,1,3
L,4,5
L,5,6
L,2,5
L,2,7
ADRAG,1, , , , , ,6
ADRAG,2, , , , , ,6
ADRAG,3, , , , , ,6
ADRAG,4, , , , , ,6
ADRAG,5, , , , , ,6
LC2=-28.750
WPOFF,,,-7.5 !偏移工作平面
csys,4 !激活工作平面
*GET,QSD,KP,,NUM,MAX !获取最大KP号并赋予变量28750QSD
*GET,QSL,LINE,,NUM,MAX
*GET,QSA,AREA,,NUM,MAX
NUMSTR,KP,KPNEXT(QSD) !为自动编号的项目建立起始编号 !函数KPNEXT(QSD)返回关键点编号大于QSD的下一个关键点编号
NUMSTR,LINE,LSNEXT(QSL)
NUMSTR,AREA,ARNEXT(QSA)
K,KPNEXT(QSD),-0.25,1.125,0 !K,ID,x,y,z
K,KPNEXT(QSD)+1,0,1.125,0
K,KPNEXT(QSD)+2,0.25,1.125,0
K,KPNEXT(QSD)+3,-0.4,-1.125,0
K,KPNEXT(QSD)+4,0,-1.125,0
K,KPNEXT(QSD)+5,0.4,-1.125,0
K,KPNEXT(QSD)+6,,-1.125,LC2
L,KPNEXT(QSD),KPNEXT(QSD)+1
L,KPNEXT(QSD)+1,KPNEXT(QSD)+2
L,KPNEXT(QSD)+3,KPNEXT(QSD)+4
L,KPNEXT(QSD)+4,KPNEXT(QSD)+5
L,KPNEXT(QSD)+1,KPNEXT(QSD)+4
L,KPNEXT(QSD)+4,KPNEXT(QSD)+6
ADRAG,LSNEXT(QSL), , , , , ,LSNEXT(QSL)+5
ADRAG,LSNEXT(QSL)+1, , , , , ,LSNEXT(QSL)+5
ADRAG,LSNEXT(QSL)+2, , , , , ,LSNEXT(QSL)+5
ADRAG,LSNEXT(QSL)+3, , , , , ,LSNEXT(QSL)+5
ADRAG,LSNEXT(QSL)+4, , , , , ,LSNEXT(QSL)+5
csys,0 !激活总体直角坐标系
wpcsys,-1,0 !将既有坐标系的XY平面定义为工作平面
LC3=-7.25
WPOFF,,,-36.25 !偏移工作平面
csys,4 !激活工作平面
*GET,QSD,KP,,NUM,MAX !获取最大KP号并赋予变量QSD
*GET,QSL,LINE,,NUM,MAX
*GET,QSA,AREA,,NUM,MAX
NUMSTR,KP,KPNEXT(QSD) !为自动编号的项目建立起始编号 !函数KPNEXT(QSD)返回关键点编号大于QSD的下一个关键点编号
NUMSTR,LINE,LSNEXT(QSL)
NUMSTR,AREA,ARNEXT(QSA)
K,KPNEXT(QSD),-0.25,1.125,0 !K,ID,x,y,z
K,KPNEXT(QSD)+1,0,1.125,0
K,KPNEXT(QSD)+2,0.25,1.125,0
K,KPNEXT(QSD)+3,-0.4,-1.125,0
K,KPNEXT(QSD)+4,0,-1.125,0
K,KPNEXT(QSD)+5,0.4,-1.125,0
K,KPNEXT(QSD)+6,,-1.125,LC3
L,KPNEXT(QSD),KPNEXT(QSD)+1
L,KPNEXT(QSD)+1,KPNEXT(QSD)+2
L,KPNEXT(QSD)+3,KPNEXT(QSD)+4
L,KPNEXT(QSD)+4,KPNEXT(QSD)+5
L,KPNEXT(QSD)+1,KPNEXT(QSD)+4
L,KPNEXT(QSD)+4,KPNEXT(QSD)+6
ADRAG,LSNEXT(QSL), , , , , ,LSNEXT(QSL)+5
ADRAG,LSNEXT(QSL)+1, , , , , ,LSNEXT(QSL)+5
ADRAG,LSNEXT(QSL)+2, , , , , ,LSNEXT(QSL)+5
ADRAG,LSNEXT(QSL)+3, , , , , ,LSNEXT(QSL)+5
ADRAG,LSNEXT(QSL)+4, , , , , ,LSNEXT(QSL)+5
NUMMRG,ALL !重新编号
NUMCMP,ALL
csys,0 !激活总体直角坐标系
wpcsys,-1,0 !将既有坐标系的XY平面定义为工作平面
ASEL,S,AREA, ,1,5 !纵梁1
CM,ZL1,AREA
ASEL,S,AREA, ,6,10 !纵梁2
CM,ZL2,AREA
ASEL,S,AREA, ,11,15 !纵梁3
CM,ZL3,AREA
KGEN,2,2,,,,,-0.62 !纵梁竖肋
KGEN,2,2,,,0.18,,-0.62
*GET,QSD1,KP,,NUM,MAX !获取最大KP号并赋予变量QSD
L,QSD1,QSD1-1
*GET,QSL1,LINE,,NUM,MAX
ADRAG,QSL1, , , , , ,5
KGEN,2,2,,,,-0.4,-5 !纵梁纵肋1
KGEN,2,2,,,0.18,-0.4,-5
KGEN,2,2,,,,-0.4,-37.4005
*GET,QSD2,KP,,NUM,MAX
L,QSD2-1,QSD2-2
L,QSD2,QSD2-2
*GET,QSL2,LINE,,NUM,MAX
ADRAG,QSL2-1, , , , , ,QSL2
NUMMRG,ALL !重新编号
NUMCMP,ALL
KGEN,2,5,,,,0.4,-34.2505 !纵梁纵肋2
KGEN,2,5,,,0.18,0.4,-34.2505
KGEN,2,5,,,,0.4,-43.5
*GET,QSD2,KP,,NUM,MAX
L,QSD2-1,QSD2-2
L,QSD2,QSD2-2
*GET,QSL2,LINE,,NUM,MAX
ADRAG,QSL2-1, , , , , ,QSL2
NUMMRG,ALL !重新编号
NUMCMP,ALL
ALLSEL !显示全部
APLOT
KGEN,2,5,,,,1.15,-4.825 !横梁
*GET,QSD1,KP,,NUM,MAX
KGEN,2,QSD1,,,,,-0.175
KGEN,2,QSD1,,,,,-0.35
L,QSD1,QSD1+1
L,QSD1+1,QSD1+2
KGEN,2,5,,,,0.7,-4.825
*GET,QSD2,KP,,NUM,MAX
KGEN,2,QSD2,,,,,-0.175
KGEN,2,QSD2,,,,,-0.35
L,QSD2,QSD2+1
L,QSD2+1,QSD2+2
KGEN,2,QSD2+1,,,3.294
L,QSD1+1,QSD2+1
L,QSD2+1,QSD2+3
*GET,QSL1,LINE,,NUM,MAX
ADRAG,QSL1-1, , , , , ,QSL1
ADRAG,QSL1-2, , , , , ,QSL1
ADRAG,QSL1-3, , , , , ,QSL1
ADRAG,QSL1-4, , , , , ,QSL1
ADRAG,QSL1-5, , , , , ,QSL1
KGEN,2,QSD2+1,,,1.647 !横梁加劲肋
*GET,QSD3,KP,,NUM,MAX
KGEN,2,QSD3,,,,,-0.120
KGEN,2,QSD3,,,,,0.120
L,QSD3,QSD3+1
L,QSD3,QSD3+2
NUMMRG,ALL !重新编号
NUMCMP,ALL
*GET,QSL2,LINE,,NUM,MAX
ADRAG,QSL2, , , , , ,QSL1-1
ADRAG,QSL2-1, , , , , ,QSL1-1
NUMMRG,ALL !重新编号
NUMCMP,ALL
AGEN,2,16,,,,,-1.38
*GET,QSD1,AREA,,NUM,MAX
AGEN,5,QSD1,,,,,-1.5
*GET,QSD2,AREA,,NUM,MAX
AGEN,2,QSD2,,,,,-1.25
*GET,QSD3,AREA,,NUM,MAX
AGEN,16,QSD3,,,,,-1.6667
*GET,QSD4,AREA,,NUM,MAX
AGEN,2,QSD4,,,,,-1.55
AGEN,3,QSD4+1,,,,,-1.6
*GET,QSD5,AREA,,NUM,MAX
AGEN,3,QSD5,,,,,-1.5
*GET,QSD6,AREA,,NUM,MAX
AGEN,2,QSD6,,,,,-1.5+0.0005
ASEL,S,AREA, ,24,25 !横梁加劲肋
CM,HLJJL,AREA
ASEL,S,AREA, ,19,23 !横梁没有加劲肋
CM,HLMYJJL,AREA
ASEL,S,AREA, ,HLJJL !横梁
ASEL,A,AREA, ,HLMYJJL
CM,HL1,AREA
ASEL,S,AREA, ,16 !纵梁竖肋
ASEL,A,AREA, ,26,52
CM,ZLSL,AREA
ASEL,S,AREA, ,17 !纵梁纵肋
ASEL,A,AREA, ,18
CM,ZLZL,AREA
AGEN,2,HL1,,,,,4.38
AGEN,2,HL1,,,,,-4.25
AGEN,2,HL1,,,,,-4.25-3*1.6667
AGEN,2,HL1,,,,,-4.25-6*1.6667
AGEN,2,HL1,,,,,-4.25-9*1.6667
AGEN,2,HL1,,,,,-4.25-12*1.6667
AGEN,2,HL1,,,,,-4.25-15*1.6667
AGEN,2,HL1,,,,,-4.25-15*1.6667-4.75
AGEN,2,HL1,,,,,-4.25-15*1.6667-4.75-4.5
CMSEL,S,HL1 !两片总横梁
ASEL,A,AREA, ,53,115
CM,HLZ,AREA
CMSEL,S,ZL1
CMSEL,A,ZL2
CMSEL,A,ZL3
CM,ZLL1,AREA
AGEN,2,ZLL1,,,3.294
AGEN,2,ZLSL,,,3.294-0.180 !两片纵梁
AGEN,2,ZLZL,,,3.294-0.180
ASEL,A,AREA, ,ALL
CM,LPZL,AREA
AGEN,2,LPZL,,,3.294*2
AGEN,2,LPZL,,,3.294*4
AGEN,2,LPZL,,,3.294*6
AGEN,2,HLZ,,,3.294 !简支纵梁
AGEN,2,HLZ,,,3.294*3
AGEN,2,HLZ,,,3.294*5
NUMMRG,ALL
NUMCMP,ALL
K,,-0.2,-1.125 !支座面
K,,0.2,-1.125
K,,0.2,-1.125,-0.8
K,,-0.2,-1.125,-0.8
*GET,ZZD,KP,,NUM,MAX
A,ZZD,ZZD-1,ZZD-2,ZZD-3
*GET,ZZM,AREA,,NUM,MAX
AGEN,8,ZZM,,,3.294
ASEL,S,AREA, ,ZZM,ZZM+7
CM,ZZZBM,AREA
AGEN,2,ZZZBM,,,,,-43.5+0.8
ASEL,S,AREA, ,ZZM+8,ZZM+15
CM,ZZYBM,AREA
NUMMRG,ALL
NUMCMP,ALL
ALLSEL !显示全部
APLOT
ET,1,SHELL181 !1=用户定义的单元类型参考号
MP,EX,1,2.06e11 !1=MAT
MP,PRXY,1,0.3
MP,dens,1,7698
MP,alpx,1,1.2e-5
sectype,1,shell,,A1 !1=SECID
secdata,0.02,1,0.0,3 !壳单元厚0.02m ,1=MAT
AATT,1,,1,,1
AESIZE,ALL,0.24
MSHAPE,0,2D !网格形状控制
MSHKEY,1 !映射
AMESH,all
K,,-0.2,-1.125 !支座面
K,,0.2,-1.125
K,,0.2,-1.125,-0.4
K,,-0.2,-1.125,-0.4
*GET,ZZD,KP,,NUM,MAX
A,ZZD,ZZD-1,ZZD-2,ZZD-3
*GET,ZZM,AREA,,NUM,MAX
AGEN,8,ZZM,,,3.294
ASEL,S,AREA, ,ZZM,ZZM+7
CM,ZZZBM,AREA
AGEN,2,ZZZBM,,,,,-43.5+0.4
NUMMRG,ALL
NUMCMP,ALL
ASEL,S,AREA, ,ZZM+8,ZZM+15
CM,ZZYBM,AREA
DA,ZZZBM,UX
DA,ZZZBM,UY
DA,ZZZBM,UZ
DA,ZZYBM,UY
DA,ZZYBM,UZ
DTRAN !将几何模型上的约束传到有限元模型
FINISH
/SOLU
!施加重力加速度
allsel
acel,,9.8
allsel,ALL
SOLVE
FINISH
/post1
|
楼主
