beam3可以定义为空心变截面梁吗？
本帖最后由 ME! 于 2013-4-20 16:59 编辑我用beam3定义了空心变截面梁，不知道下面的命令流对不对？
finish
/clear
/prep7
RO = 7850
PEX = 2.06E11
MP,EX,1,PEX !ex: 弹性模量
MP,DENS,1,RO!密度
ET,1,beam3,,,
*DIM,AERA,ARRAY,16!定义横截面面积数组
AERA(1)= 0.0000636
AERA(2)=0.0001374
AERA(3)=0.0001633
AERA(4)= 0.0001767
AERA(5)=0.0002798
AERA(6)=0.0006103
AERA(7)=0.0003196
AERA (8) = 0.0013188
AERA(9)= 0.0002961
AERA(10)= 0.0001767
AERA(11)= 0.0001767
AERA(12)= 0.0002798
AERA(13) =0.0001767
AERA(14) =0.0001767
AERA (15) =0.0001508
AERA(16)= 0.0000725
*DIM,DIAM,ARRAY,16!定义直径数组
DIAM(1) = 15E-2
DIAM(2) = 16E-2
DIAM(3) = 17E-2
DIAM(4) = 17E-2
DIAM(5) = 20.5E-2
DIAM(6) = 29E-2
DIAM(7) = 21.7E-2
DIAM(8) = 21E-2
DIAM(9) = 20E-2
DIAM(10) = 17E-2
DIAM(11) = 17E-2
DIAM(12) = 20.5E-2
DIAM(13) = 17E-2
DIAM(14) = 17E-2
DIAM(15) = 16E-2
DIAM(16) = 12.5E-2
*DIM,IZZ,ARRAY,16!定义截面惯性矩数组
IZZ(1) = 0.01467E-7
IZZ(2) =0.02895E-7
IZZ(3) =0.03778E-7
IZZ(4) =0.03899E-7
IZZ(5) =0.08468E-7
IZZ(6) =0.3452E-7
IZZ(7) = 0.1068E-7
IZZ(8) =0.09345E-7
IZZ(9) = 0.07653E-7
IZZ(10) =0.03899E-7
IZZ(11) =0.03899E-7
IZZ(12) =0.08468E-7
IZZ(13) = 0.03899E-7
IZZ(14) =0.03899E-7
IZZ(15) =0.03016E-7
IZZ(16) =0.009974E-7
R,1,AERA,IZZ,DIAM
N,1,0
N,2,0.0117
N,3,0.0117*2
N,4,0.0117*3
N,5,0.0117*4
N,6,0.0117*5
N,7,0.0117*6
N,8,0.0117*7
N,9,0.0117*8
N,10,0.0117*9
N,11,0.0117*10
N,12,0.0117*11
N,13,0.0117*12
N,14,0.0117*13
N,15,0.0117*14
N,16,0.0117*15
N,17,0.0117*16
TYPE,1
E,1,2
E,2,3
E,3,4
E,4,5
E,5,6
E,6,7
E,7,8
E,8,9
E,9,10
E,10,11
E,11,12
E,12,13
E,13,14
E,14,15
E,15,16
E,16,17
finish
/SOLU
D,1,ALL
D,17,ALL
ANTYPE,MODAL
MODOPT,lanb,20,,,ON
MXPAND,20
SOLVE
finish
输入ansys计算，出现element1 has a zero aera,check real constant table1这个问题，不知道是什么原因，恳请高手！
下面是我用beam3的实心梁截面，运行没有问题
finish
/clear
/prep7
RO = 7850
PEX = 2.06E11
MP,EX,1,PEX !ex: 弹性模量
MP,DENS,1,RO!密度
ET,1,beam3,,,
R,1,7.85E-5, 4.91e-10,0.01
N,1,0
N,2,0.0117
N,3,0.0117*2
N,4,0.0117*3
N,5,0.0117*4
N,6,0.0117*5
N,7,0.0117*6
N,8,0.0117*7
N,9,0.0117*8
N,10,0.0117*9
N,11,0.0117*10
N,12,0.0117*11
N,13,0.0117*12
N,14,0.0117*13
N,15,0.0117*14
N,16,0.0117*15
N,17,0.0117*16
TYPE,1
E,1,2
E,2,3
E,3,4
E,4,5
E,5,6
E,6,7
E,7,8
E,8,9
E,9,10
E,10,11
E,11,12
E,12,13
E,13,14
E,14,15
E,15,16
E,16,17
finish
/SOLU
D,1,ALL
D,17,ALL
ANTYPE,MODAL
MODOPT,lanb,20,,,ON
MXPAND,20
SOLVE
finish
请问高手beam3可以建模成为空心变截面梁单元吗？我看到beam3的实常数是A,I,h，若是圆形截面h应该怎么设置呢，实心的时候我设置的h=d，但是空心的时候，我不知道怎么设置了，我设置为外径，也提示上面的错误
终于发现错在哪里了，现在把正确的贴上来
FINISH
/CLEAR
/prep7
MP,EX,1,2.06E11!ex: 弹性模量
MP,DENS,1,7850!密度
MP,GXY,1,1E12 !gxy为剪切模量
MP,PRXY,1,0.3!泊松比
ET,1,beam3,,,
*DIM,AERA,ARRAY,16!定义横截面面积数组
AERA(1)= 0.0000636
AERA(2)=0.0001374
AERA(3)=0.0001633
AERA(4)= 0.0001767
AERA(5)=0.0002798
AERA(6)=0.0006103
AERA(7)=0.0003196
AERA (8) = 0.0013188
AERA(9)= 0.0002961
AERA(10)= 0.0001767
AERA(11)= 0.0001767
AERA(12)= 0.0002798
AERA(13) =0.0001767
AERA(14) =0.0001767
AERA (15) =0.0001508
AERA(16)= 0.0000725
*DIM,DIAM,ARRAY,16!定义直径数组
DIAM(1) = 15E-2
DIAM(2) = 16E-2
DIAM(3) = 17E-2
DIAM(4) = 17E-2
DIAM(5) = 20.5E-2
DIAM(6) = 29E-2
DIAM(7) = 21.7E-2
DIAM(8) = 21E-2
DIAM(9) = 20E-2
DIAM(10) = 17E-2
DIAM(11) = 17E-2
DIAM(12) = 20.5E-2
DIAM(13) = 17E-2
DIAM(14) = 17E-2
DIAM(15) = 16E-2
DIAM(16) = 12.5E-2
*DIM,IZZ,ARRAY,16!定义截面惯性矩数组
IZZ(1) = 0.01467E-7
IZZ(2) =0.02895E-7
IZZ(3) =0.03778E-7
IZZ(4) =0.03899E-7
IZZ(5) =0.08468E-7
IZZ(6) =0.3452E-7
IZZ(7) = 0.1068E-7
IZZ(8) =0.09345E-7
IZZ(9) = 0.07653E-7
IZZ(10) =0.03899E-7
IZZ(11) =0.03899E-7
IZZ(12) =0.08468E-7
IZZ(13) = 0.03899E-7
IZZ(14) =0.03899E-7
IZZ(15) =0.03016E-7
IZZ(16) =0.009974E-7
*DO,I,1,16
R,1,AERA(I),IZZ(I),DIAM(I)
*ENDDO
N,1,0
N,2,0.0165
N,3,0.03
N,4,0.034
N,5,0.044
N,6,0.049
N,7,0.0555
N,8,0.0585
N,9,0.065
N,10,0.1122
N,11,0.1172
N,12,0.1377
N,13,0.1427
N,14,0.1527
N,15,0.1575
N,16,0.164
N,17,0.187
TYPE,1
*DO,I,1,16
E,I,I+1
*ENDDO
*DO,I,1,16
SECTYPE,I,BEAM,CTUBE !截面类型ID: 截面号TYPE: BEAM:定义此截面用于梁CSOLID:圆形实心截面CTUBE: 圆管
SECDATA,0.004,DIAM(I)/2,16 !Ri, Ro, N可是输入的时候变为外半径，内半径了？？？？？？？？？
*ENDDO
eplot
finish
/SOLU
D,1,ALL
D,17,ALL
/SOLU
ANTYPE,MODAL
MODOPT,lanb,20,,,ON
MXPAND,20
SOLVE
页:
[1]