为什么我用ansys命令流建立beam3模型，定义了实常数，去得面积为
本帖最后由 ME! 于 2013-7-7 21:46 编辑finish
/clear
/prep7
RO = 7850
PEX = 2.06E11
MP,EX,1,PEX !ex: 弹性模量
MP,DENS,1,RO!密度
MP,PRXY,1,0.3
ET,1,beam3,,,
ET,32,mass21 !mass ux uy uz
keyopt,32,3,4
R,32,0.0248,0.0248,0.0248
ET,33,mass21
keyopt,33,3,4
R,33,0.0577,0.0577,0.0577
ET,34,mass21
keyopt,34,3,4
R,34,0.0577,0.0577,0.0577,
ET,38,mass21
keyopt,38,3,4
R,38,0.134,0.134,0.134
ET,313,mass21
keyopt,313,3,4
R,313,0.0577,0.0577,0.0577
ET,314,mass21
keyopt,314,3,4
R,314,0.0577,0.0577,0.0577
ET,315,mass21
keyopt,315,3,4
R,315,0.02051,0.02051,0.02051
ET,316,mass21
keyopt,316,3,4
R,316,0.06647,0.06647,0.06647
ET,22,COMBIN14
keyopt,22,2,0
R,22,164e6
ET,23,COMBIN14
keyopt,23,2,0
R,23,164e6
*DIM,DIAM,ARRAY,16!定义直径数组
DIAM(1) = 15E-3
DIAM(2) = 16E-3
DIAM(3) = 17E-3
DIAM(4) = 17E-3
DIAM(5) = 20.5E-3
DIAM(6) = 29E-3
DIAM(7) = 21.7E-3
DIAM(8) = 21E-3
DIAM(9) = 20E-3
DIAM(10) = 17E-3
DIAM(11) = 17E-3
DIAM(12) = 20.5E-3
DIAM(13) = 17E-3
DIAM(14) = 17E-3
DIAM(15) = 16E-3
DIAM(16) = 12.5E-3
*DIM,AERA,ARRAY,16!定义横截面面积数组
AERA(1)= 0.000126449104
AERA(2)=0.000150796447
AERA(3)=0.000176714587
AERA(4)= 0.000176714587
AERA(5)=0.0002797980
AERA(6)=0.000610254372
AERA(7)=0.0003195707
AERA (8) = 0.0002960951
AERA(9)= 0.0002638938
AERA(10)= 0.00017671459
AERA(11)= 0.00017671459
AERA(12)= 0.00027979808
AERA(13) =0.00017671459
AERA(14) =0.00017671459
AERA (15) =0.00015079645
AERA(16)= 0.00007252981
*DIM,IZZ,ARRAY,16!定义截面惯性矩数组
IZZ(1) = 0.022839869465450E-7
IZZ(2) = 0.030159289474462E-7
IZZ(3) =0.038987655704901E-7
IZZ(4) =0.038987655704901E-7
IZZ(5) = 0.084682642404833E-7
IZZ(6) = 0.084682642404833E-7
IZZ(7) = 0.106834468521627E-7
IZZ(8) = 0.093455018336514E-7
IZZ(9) =0.076529197041447E-7
IZZ(10) =0.038987655704901E-7
IZZ(11) = 0.038987655704901E-7
IZZ(12) = 0.084682642404833E-7
IZZ(13) =0.038987655704901E-7
IZZ(14) =0.038987655704901E-7
IZZ(15) = 0.030159289474462E-7
IZZ(16) = 0.009973605607059E-7
*DO,I,1,16
R,I,AERA(I),IZZ(I),DIAM(I)
*ENDDO
N,1,0
N,2,0.0164
N,3,0.029
N,4,0.039
N,5,0.049
N,6,0.0555
N,7,0.0585
N,8,0.065
N,26,0.09
N,9,0.1122
N,10,0.1172
N,11,0.1222
N,12,0.1272
N,13,0.1377
N,14,0.1477
N,15,0.1577
N,16,0.164
N,17,0.187
TYPE,1
*DO,I,1,16
E,I,I+1
R,I
*ENDDO
TYPE,32
real,32
e,2
TYPE,33
real,33
e,3
TYPE,34
real,34
e,4
TYPE,38
real,38
e,26
TYPE,313
real,313
e,13
TYPE,314
real,314
e,14
TYPE,315
real,315
e,15
TYPE,316
real,316
e,16
BRG =-0.010
N,18,0.039,BRG
N,20,0.049,BRG
N,22,0.1477,BRG
N,24,0.1577,BRG
TYPE,22 ! CREATE BEARING
REAL,22
E,18,4
E,20,5
TYPE,23 ! CREATE BEARING
REAL,23
E,22,14
E,24,15
finish
/SOLU
D,18,ALL
D,20,ALL
D,22,ALL
D,24,ALL
/SOLU
ANTYPE,MODAL
MODOPT,lanb,20,,,ON
MXPAND,20
SOLVE
finish
TYPE,1
*DO,I,1,16
E,I,I+1
R,I
*ENDDO
是不是红色部分字体有问题啊，不能这么定义吗？老是弹出警告
*** WARNING *** CP = 324.547 TIME= 21:12:49
Zero thickness element 1 may result in wrong thermal or bending stress.
Check real constant table 1.
*** ERROR *** CP = 324.547 TIME= 21:12:49
Element 1 has a zero area.Check real constant table 1.
*** WARNING *** CP = 324.703 TIME= 21:12:53
Zero moment of inertia for element 1.Check real constant table 1.
页:
[1]