beam3可以定义为空心变截面梁吗？

ME!

我用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，但是空心的时候，我不知道怎么设置了，我设置为外径，也提示上面的错误

ME!

终于发现错在哪里了，现在把正确的贴上来

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