24小时热门版块排行榜

返回列表

王禄友

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 6.5
帖子: 5
在线: 1.2小时
虫号: 2696818
注册: 2013-10-03
专业: 机械动力学

[求助] ansys分析不收敛

/prep7
!定义线的半径、线的长度、线间距、计算线根数
Bro=0.1
Bl=20
Bs=0.5
Nbeam=bl/bs
!定义单元类型，beam188定义两种单元类型为了便于选择操作，而接触单元的keyopt不同
!定义单元实常数和材料常数、线线摩擦系数、线截面
Et,1,beam188
Et,2,beam188
Et,3,conta176
Et,4,targe170
Keyopt,3,3,1
Keyopt,3,12,2
R,1
R,2,bro,bro,0.001
Mp,ex,1,1.84e5
Mp,prxy,1,0.29
Mp,dens,1,7.85e-6
Mp,mu,2,0.1
Sectype,1,beam,csolid
Secdata,bro
!创建梁模型：此种方法可保证线方向一致，而布尔运算则不一定
!圆截面可不设方位点，但设置方位点可保证变形后的单元显示正确
K,1,,,bs/2
*do,i,2,nbeam+1
K,i,bs/2+(i-2)*bs,,bs/2
*enddo
K,nbeam+2,bl,,bs/2
*do,i,1,nbeam+1
L,i,i+1
*enddo
Lgen,nbeam,all,,,,,bs
Cm,cml1,line
Lsel,none
Kmax=kpinqr(0,14)
K,kmax+1,bs/2
*do,i,2,nbeam+1
K,kmax+i,bs/2,,bs/2+(i-2)*bs
*enddo
K,kmax+nbeam+2,bs/2,,bl
*do,i,1,nbeam+1
L,kmax+i,kmax+i+1
*enddo
Lgen,nbeam,all,,,bs
Cm,cml2,line
Lsel,all
Kmax=kpinqr(0,14)+1
K,kmax,0,1000*b
Lsel,s,tan1,x
Latt,1,1,1,,,kmax,1
Lsel,s,tan1,z
Latt,1,1,2,,,kmax,1
Lsel,all
Esize,bs/8
Lmesh,all
!创建纵线与横线之间的接触对，通过实常数2识别
Type,3
Real,2
Mat,2
Esel,s,type,,1
Nsle
Esurf
Esel,s,type,,2
Nsle
Type,4
Esurf
!施加边界条件
Nsel,s,loc,x,0
Nsel,a,loc,x,bl
Nsel,a,loc,z,0
Nsel,a,loc,z,bl
D,all,all
Cm,dnode1,node
Allsel,all
!在交叉处的重合节点施加反向力对，使其发生上下穿越变形，此力可略大些，可先脱空穿越
*do,i,1,nbeam
*do,j,1,nbeam
*if,mod(i+j,2),eq,0,then
Ftemp=125
*else
Ftemp=-125
*endif
Esel,s,type,,1
Nsle,s
F,node(bs/2+(i-1)*bs,0,bs/2+(j-1)*bs),fy,ftemp
Esel,s,type,,2
Nsle,s
F,node(bs/2+(i-1)*bs,0,bs/2+(j-1)*bs),fy,-ftemp
*enddo
*enddo
!进入求解层，设置求解选项；第一载荷步杀死全部接触对并求解
/solu
antype,static
rescontrol,define,all,1
Nlgeom,on
Outres,all,all
Nsubst,20,,5
Time,1
Esel,s,type,,3,4
Ekill,all
Allsel,all
solve
!第二载荷步定义收敛准则及载荷子步，激活纵横梁之间的接触对，将力对设置为0值，求解
time,2
cnvtol,f,,0.005
nsubst,100,,20
esel,s,type,,3,4
ealive,all
esel,s,type,,1,2
nsle,s
cmsel,u,dnode1
f,all,fy,0
allsel,all
solve
finish
如果把线的长度bl改为10就可以收敛了，但是真实情况是45呀，新人，只有这么多金币，求大神帮忙！！！

回复此楼