在弹塑性分析中需要设置实效应变值,看看弹塑性理论的应力应变关系说明中有这一概念。这是别人做的例子
用Q235材料的分段线性塑性模型做了一个简易的碰撞分析,
当我改变失效应变的值(0.19到0.75),时,计算结果中,第一方案有
失效单元出现,第二方案中无失效单元,但碰撞结果中两方案的能量吸收,
最大界面力及位移结果一样,想请教一下这是什么原因,我在真实应力-应变曲线中
定义的最大应变为0.19,是不是失效都以此为限,当最大应变到0.19后,即使定义的
失效应变为更高值,计算就到此停止。
附上模型文件,请帮助查一下有没有需要改进的地方,谢谢!
finish
/clear,start
/filname,D50PL8-019
/title,analysis of pipe-buckling
/com,施加位移(长100mm,直径50mm,施加70mm位移)
/prep7
et,1,shell163
keyopt,1,1,10
r,1,5/6,5,1.5e-3,,,,!!!!管壁厚1.5mm
r,2,,3,0.01,,,,
MP,DENS,1,7860
MP,EX,1,2.1E+011
MP,NUXY,1,0.31
TB,PLAW,,,,8,
TBDAT,1,
TBDAT,2,
TBDAT,3,0.19
TBDAT,4,40
TBDAT,5,5
TBDAT,6,1
*DIM,TruStran,,7
*DIM,TruStres,,7
TruStran(1)=0,0.0237,0.0477,0.0711,0.0941,0.1380,0.1900
TruStres(1)=285e6,339e6,378e6,405e6,428e6,461e6,491e6
EDCURVE,ADD,1,TruStran(1),TruStres(1)
k,1,0,0,0
k,2,0,0,0.100
lstr,1,2
circle,1,0.025
adrag,2,3,4,5,,,1
aatt,1,1,1,0,
esize,0.0035
amesh,all
nsel,s,loc,z,-0.001,0.001
d,all,,,,,,all,,,,,
allsel,all
allsel,all
LDELE, 1, , ,1
CM,tube,NODE
allsel,all
WPCSYS,-1,0
wpoffs,,,.105
CYL4, , ,0.035
EDMP,RIGI,2,4,7
MP,DENS,2,3000000
MP,EX,2,2e11
MP,NUXY,2,0.3
TYPE, 1
MAT, 2
REAL, 2
ESYS, 0
SECNUM,
MSHAPE,0,2D
AMESH,5,
ASEL,S, , , 5
ALLSEL,BELOW,AREA
NSLA,S,1
CM,wall,NODE
allsel,all
EDPART,CREATE
EDPV,VELO, 2,0,0,-8,0,0,0,
allsel,all
EDMP,HGLS,1,3,0.1,1.5,0.05, , ,
EDMP,HGLS,2,3,0.1,1.5,0.05, , ,
EDCGEN,ANTS,1,2,0.4,0.3,0,0,0, , , , ,0,10000000,0,0
EDCGEN,ASSC, , ,0.4,0.3,0,0,0, , , , ,0,10000000,0,0
/solu
time,0.020,
EDOPT,ADD,blank,ANSYS
edrst,50,
edhtime,500,
EDHIST,tube
EDHIST,wall
EDOUT,GLSTAT
EDOUT,RWFORCE
EDOUT,MATSUM
EDOUT,RCFORC
EDOUT,RBDOUT
edint,5,0,
EDENERGY,1,1,1,1
SAVE
SOLVE
finish |