求大神,关于模拟砖块受压不收敛的问题!!
浏览:10984
本人用ansys10.0模拟砖块受压,砖块用solid65单元①在求解控制阶段,本人用弧长法,能顺利跑出结果。鉴于本人以后打算用接触单元模拟砖墙,故不打算用弧长法求解。
②在求解阶段重新编写命令流后却无法跑出结果。求大神指导~!
finish
/clear
/prep7
l=240
b=115
h=53
et,1,solid65 !!!!定义单元类型
keyopt,1,1,1 !不计形函数附加项
keyopt,1,7,1 !开裂后拉应力释放系数,慢慢释放拉应力
keyopt,1,3,2
mp,ex,1,13429 !砖
mp,prxy,1,0.17
fc=10.83 !砖块抗压强度
ft=1.0 !砖块抗拉强度
tb,concr,1
tbdata,,0.3,0.85,1.0,-1 !定义砖相关参数,关闭压碎(fc=-1)
tb,miso,1,,12 !砖墙应力应变曲线,采用miso模型
tbpt,,0.00015,fc*0.186
tbpt,,0.0003,fc*0.354
tbpt,,0.00045,fc*0.502
tbpt,,0.0006,fc*0.63
tbpt,,0.00075,fc*0.74
tbpt,,0.0009,fc*0.83
tbpt,,0.00105,fc*0.902
tbpt,,0.0012,fc*0.954
tbpt,,0.00135,fc*0.986
tbpt,,0.0015,fc
tbpt,,0.00165,fc*0.98
tbpt,,0.0018,fc*0.96
et,2,solid45
mp,ex,2,2.1e5
mp,prxy,2,0.3
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!建模
block,0,240,0,10,0,115
block,0,240,10,63,0,115
block,0,240,63,73,0,115
vglue,1,2,3
esize,20
vsel,s,loc,y,10,63
vatt,1,,1
vsel,inve
vatt,2,,2
vsel,all
vmesh,all
!---------------------------------
asel,s,loc,y,0
da,all,all
asel,s,loc,y,73
da,all,ux
da,all,uz
sfa,all,,pres,30
allsel,all
!-----------------------------------
/solu
antype,0 !定义分析类型
NLGEOM,Off !定义小变形
NROPT,FULL,,ON !定义牛顿-拉弗森选项
cnvtol,u,,0.05 !定义收敛准则为位移
cnvtol,f,,0.05
solcontrol,on !打开非线性求解设置与算法控制
outres,all,all
autots,off !关闭自动时间步控制
lnsrch,1 !打开线性搜索
ncnv,2 !如果不收敛时结束而不退出
time,100 !定义第1载荷步
nsub,150 !指定子步数
neqit,30 !每一子步中方程的迭代次数限值
allsel,all
solve
②在求解阶段重新编写命令流后却无法跑出结果。求大神指导~!
finish
/clear
/prep7
l=240
b=115
h=53
et,1,solid65 !!!!定义单元类型
keyopt,1,1,1 !不计形函数附加项
keyopt,1,7,1 !开裂后拉应力释放系数,慢慢释放拉应力
keyopt,1,3,2
mp,ex,1,13429 !砖
mp,prxy,1,0.17
fc=10.83 !砖块抗压强度
ft=1.0 !砖块抗拉强度
tb,concr,1
tbdata,,0.3,0.85,1.0,-1 !定义砖相关参数,关闭压碎(fc=-1)
tb,miso,1,,12 !砖墙应力应变曲线,采用miso模型
tbpt,,0.00015,fc*0.186
tbpt,,0.0003,fc*0.354
tbpt,,0.00045,fc*0.502
tbpt,,0.0006,fc*0.63
tbpt,,0.00075,fc*0.74
tbpt,,0.0009,fc*0.83
tbpt,,0.00105,fc*0.902
tbpt,,0.0012,fc*0.954
tbpt,,0.00135,fc*0.986
tbpt,,0.0015,fc
tbpt,,0.00165,fc*0.98
tbpt,,0.0018,fc*0.96
et,2,solid45
mp,ex,2,2.1e5
mp,prxy,2,0.3
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!建模
block,0,240,0,10,0,115
block,0,240,10,63,0,115
block,0,240,63,73,0,115
vglue,1,2,3
esize,20
vsel,s,loc,y,10,63
vatt,1,,1
vsel,inve
vatt,2,,2
vsel,all
vmesh,all
!---------------------------------
asel,s,loc,y,0
da,all,all
asel,s,loc,y,73
da,all,ux
da,all,uz
sfa,all,,pres,30
allsel,all
!-----------------------------------
/solu
antype,0 !定义分析类型
NLGEOM,Off !定义小变形
NROPT,FULL,,ON !定义牛顿-拉弗森选项
cnvtol,u,,0.05 !定义收敛准则为位移
cnvtol,f,,0.05
solcontrol,on !打开非线性求解设置与算法控制
outres,all,all
autots,off !关闭自动时间步控制
lnsrch,1 !打开线性搜索
ncnv,2 !如果不收敛时结束而不退出
time,100 !定义第1载荷步
nsub,150 !指定子步数
neqit,30 !每一子步中方程的迭代次数限值
allsel,all
solve