求助-The degree of freedom solution is not available!
求解不收敛，后处理总显示The degree of freedom solution is not available!我是初学者，问题太多，请大家帮忙看一下！
/filname,wlszy_gj,1
/title,daocha !道岔
/prep7
/input,vread,txt
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~钢轨节点~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!钢轨
*do,i,1,645
n,i,gg(1,i),gg(2,i)
*enddo
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~钢轨属性~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
et,1,link1
r,1,77.45e-4
mp,ex,1,2.1e11
mp,alpx,1,11.8e-6 !线膨胀系数
type,1
real,1
mat,1
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~钢轨单元~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
*do,i,1,80
e,i,i+1
e,81+i,82+i
*enddo
*do,i,1,65
e,162+i,163+i
e,228+i,229+i
*enddo
*do,i,1,50
e,339+i,340+i
e,390+i,391+i
e,441+i,442+i
e,492+i,493+i
e,543+i,544+i
e,594+i,595+i
*enddo
*do,i,1,27
e,294+i,295+i
*enddo
*do,i,1,16
e,322+i,323+i
*enddo
e,1,340
e,82,391
e,81,442
e,322,544
e,339,595
e,162,493
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~轨枕节点~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
nsel,s,node,,all
nsel,u,node,,163,180
nsel,u,node,,229,246
nsel,u,node,,295,305
*get,ggnmin,node,,num,min
*get,ggnzs,node,,count
ggnh=ggnmin
*do,i,1,ggnzs
ngen,2,1000,ggnh,,1,,,0.05
ggnh=ndnext(ggnh)
*enddo
ngen,2,1000,1001,1081,1,,-0.5325
ngen,2,1000,1340,1390,1,,-0.5325
ngen,2,1000,1442,1492,1,,-0.5325
ngen,2,1000,1317,1322,1,,0.49
ngen,2,1000,1544,1594,1,,0.5325
ngen,2,1000,1391,1441,1,,0.5325
ngen,2,1000,1082,1156,1,,0.5325
ngen,2,1000,1334,1339,1,0.00697,-0.5001
ngen,2,1000,1157,1162,1,-0.00697,0.5001
ngen,2,1000,1595,1645,1,0.00742,-0.5324
ngen,2,1000,1493,1543,1,-0.00742,0.5324
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~轨枕属性~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
et,2,beam3
r,2,0.0704,0.000284,0.22
mp,ex,2,34.5e9
mp,prxy,2,0.167!泊松比
mp,dens,2,2500 !容重
type,2
real,2
mat,2
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~轨枕单元~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
*do,k,1,126
nsel,s,loc,z,0.02,0.07
nsel,r,loc,x,-30.9+(k-1)*0.6,-30.3+(k-1)*0.6
*get,gznzs,node,,count
*get,gznmin,node,,num,min
gznh=gznmin
*dim,y,array,gznzs
*dim,n,array,gznzs
*do,i,1,gznzs
y(i)=ny(gznh)
gznh=ndnext(gznh)
*enddo
*do,i,1,gznzs-1
*do,j,1,gznzs-i
*if,y(j),gt,y(j+1),then
yy=y(j+1)
y(j+1)=y(j)
y(j)=yy
*endif
*enddo
*enddo
*do,i,1,gznzs
n(i)=node(,y(i),,)
*enddo
*do,i,1,gznzs-1
e,n(i),n(i+1)
*enddo
y=
n=
*enddo
!~~~~~~~~~~~~~
*do,k,1,57
nsel,s,loc,z,0.02,0.07
nsel,r,loc,y,-0.6,1.9676
nsel,u,node,,2334
nsel,r,loc,x,44.7+(k-1)*0.6,45.3+(k-1)*0.6
*get,gznzs,node,,count
*get,gznmin,node,,num,min
gznh=gznmin
*dim,y,array,gznzs
*dim,n,array,gznzs
*do,i,1,gznzs
y(i)=ny(gznh)
gznh=ndnext(gznh)
*enddo
*do,i,1,gznzs-1
*do,j,1,gznzs-i
*if,y(j),gt,y(j+1),then
yy=y(j+1)
y(j+1)=y(j)
y(j)=yy
*endif
*enddo
*enddo
*do,i,1,gznzs
n(i)=node(,y(i),,)
*enddo
*do,i,1,gznzs-1
e,n(i),n(i+1)
*enddo
y=
n=
*enddo
!~~~~~~~~~~~~~
*do,k,1,57
nsel,s,loc,z,0.02,0.07
nsel,r,loc,y,1.97,7.7
nsel,a,node,,2334
nsel,r,loc,x,44.7+(k-1)*0.6,45.3+(k-1)*0.6
*get,gznzs,node,,count
*get,gznmin,node,,num,min
gznh=gznmin
*dim,y,array,gznzs
*dim,n,array,gznzs
*dim,x,array,gznzs
*do,i,1,gznzs
y(i)=ny(gznh)
x(i)=nx(gznh)
gznh=ndnext(gznh)
*enddo
*do,i,1,gznzs-1
*do,j,1,gznzs-i
*if,y(j),gt,y(j+1),then
yy=y(j+1)
y(j+1)=y(j)
y(j)=yy
*endif
*enddo
*enddo
*do,i,1,gznzs-1
*do,j,1,gznzs-i
*if,x(j),lt,x(j+1),then
xx=x(j+1)
x(j+1)=x(j)
x(j)=xx
*endif
*enddo
*enddo
*do,i,1,gznzs
n(i)=node(x(i),y(i),,)
*enddo
*do,i,1,gznzs-1
e,n(i),n(i+1)
*enddo
x=
y=
n=
*enddo
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~扣件属性~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
et,3,combin39
r,3,0,0,0.5e-3,12e3,1e-3,16e3
rmore,5e-3,16.1e3
type,3
real,3
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~扣件单元~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
esel,s,real,,1
nsle,s,all
nsel,u,node,,163,180
nsel,u,node,,229,246
nsel,u,node,,295,305
*get,ggnmin,node,,num,min
*get,ggnzs,node,,count
ggnh=ggnmin
*do,i,1,ggnzs
e,ggnh,1000+ggnh
ggnh=ndnext(ggnh)
*enddo
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~传力部件~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
et,4,combin39 !限位器
r,4,0,0,7e-3,1e-10,7.0001e-3,2000e3
type,4
real,4
e,26,183
e,107,249
et,5,combin39 !间隔铁
r,5,0.5e-3,24e3,1e-3,60e3,1.5e-3,140e3
rmore,2e-3,250e3,2.5e-3,395e3
type,5
real,5
*do,i,1,4
e,307+i,223+i
e,324+i,289+i
*enddo
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~道床阻力&路基约束~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
esel,s,real,,2
nsle,s,all
*get,nmin,node,,num,min
*get,nzs,node,,count
nh=nmin
*do,i,1,nzs
nsel,s,node,,nh
esln,s
*get,enmin,elem,0,num,min
*get,enmax,elem,0,num,max
*if,enmin,eq,enmax,then
*get,ec,elem,enmin,leng
zl5=ec/2*4300.1
zl4=ec/2*4300
zl3=ec/2*3900
zl2=ec/2*3400
zl1=ec/2*2600
ngen,2,2000,nh,,,,,0.1
nsel,s,node,,2000+nh !道床约束
d,2000+nh,all
et,6,combin39
r,6,0,0,0.5e-3,zl1,1e-3,zl2
rmore,1.5e-3,zl3,2e-3,zl4,4e-3,zl5
type,6
real,6
e,nh,2000+nh
*else
*get,ec1,elem,enmin,leng
*get,ec2,elem,enmax,leng
ec=(ec1+ec2)/2
zl5=ec/2*4300.1
zl4=ec/2*4300
zl3=ec/2*3900
zl2=ec/2*3400
zl1=ec/2*2600
ngen,2,2000,nh,,,,,0.1
nsel,s,node,,2000+nh !道床约束
d,2000+nh,all
et,6,combin39
r,6,0,0,0.5e-3,zl1,1e-3,zl2
rmore,1.5e-3,zl3,2e-3,zl4,4e-3,zl5
type,6
real,6
e,nh,2000+nh
*endif
esel,s,real,,2
nsle,s,all
nh=ndnext(nh)
*enddo
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~边界条件&力~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
esel,s,real,,2 !y方向约束
nsle,s,all
nsel,r,loc,y,-0.6,-0.5
nsel,a,node,,2157,2162
nsel,a,node,,2493,2543
d,all,uy
esel,s,real,,1
nsle,s,all
d,all,uy
f,390,fx,948762.5 !固定区温度力
f,441,fx,948762.5
f,492,fx,-948762.5
f,594,fx,-948762.5
f,645,fx,-948854.7
f,543,fx,-948854.7
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~温度荷载~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
/solu
esel,s,real,,1 !钢轨升温50度
bfe,all,temp,,50
ALLSEL,ALL
/replot
solve !求解
我是初学者，问题太多，请大家帮忙看一下。wangbin_jd@163.com 是不是没有约束好
有时候模型太大了也不显示
页:
[1]