ME! 发表于 2013-4-20 15:28

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

ME! 发表于 2013-5-4 10:22

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

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]
查看完整版本: beam3可以定义为空心变截面梁吗?