/units,si
/FILNAME,tajia,0
/TITLE,tower modal
/PREP7 !进入前处理器
*SET,D1,1.5 !定义参变量塔顶直径D1、塔底直径D2、塔高H、塔顶壁厚T1和塔底厚度T2
*SET,D2,3
*SET,H,67.8
*SET,T1,0.015
*SET,T2,0.03
ET,1,SHELL93 !定义材料属性和单元类型
MP,EX,1,2e11
MP,PRXY,1,0.3
MP,DENS,1,7850
R,1,T1,T2,T2,T1,,,
CSYS,1
K,1,D2/2,0,, !创建构成塔架1/4
K,2,D2/2,90,,
K,3,D1/2,90,H,
K,4,D1/2,0,H,
LSTR,1,4
LSTR,2,3
L,1,2
L,3,4
AL,1,3,2,4
/PNUM,AREA,1
APLOT !建立塔架1/4
CSYS,0 !建立整个塔架
ARSYM,X,1, , , ,0,0
ARSYM,Y,1,2, , ,0,0 !
K,1001,0,0,1.5,
KWPAVE,1001
CSYS,4
WPROTA,,,-90.000000
WPROTA,,-45,,
/VIEW,1,WP
CSYS,4
K,2001,0.75,0.25,0,
K,2002,0.75,-0.25,0,
K,2003,-0.75,-0.25,0,
K,2004,-0.75,0.25,0,
/PNUM,KP,1
/PNUM,LINE,1
KPLOT
L,2001,2004
L,2002,2003
L,2001,2002
L,2003,2004
LFILLT,17,20,0.25, ,
LFILLT,20,18,0.25, ,
LFILLT,18,19,0.25, ,
LFILLT,19,17,0.25, ,
LPLOT
AL,17,21,22,18,20,23
VOFFST,5,3,100,
/PNUM,AREA,1
APLOT
ASBV,2,1
NUMMRG,ALL, , , ,LOW
APLOT !
AATT,1,1,1,0
LESIZE,ALL,0.5
MSHAPE,1,2D
AMESH,ALL !
CSYS,0 !
K,9999,0,0,67.8
ET,2,MASS21
R,1,1E-6
KMESH,ALL
NSEL,S,LOC,Z,67.8,67.8
CERIG,NODE(0,0,67.8),ALL,ALL,,,,
/SOL
ANTYPE,0
DL,3, ,ALL,
DL,7, ,ALL,
DL,11, ,ALL,
DL,15, ,ALL,
ACEL,0,0,9.8,
F,9229,FY,385732
F,9229,FZ,-891800
F,9229,MX,-1474276
nsel,all
/STATUS,SOLU
SOLVE
FINISH
/POST1
PLNSOL, S,EQV, 2,1.0
我建立的是圆筒型的塔架,即:用一定厚度的板卷成的。
问题:
我进行的是静力分析,其结果为什么与实际的完全不一样,位移达到了34398,最大应力达到0.867e13,在整个的操作过程中,我都是应用的国际单位,其结果为什么会这样,产生的原因是什么,怎样修改?请指教! |