xugang2000 发表于 2008-2-20 10:50:23

路面板的分析，考虑传力杆

!此处是做的一个路面板的分析，考虑的是3*3块板，见图
!板间的传力杆的连接见图，考虑ux,uy,uz的耦合
!对于耦合之后再产生弹簧单元，弹簧单元在耦合节点处的弹簧系数取值有点疑惑
!因为在这个地方产生两个弹簧单元，我得弹性系数是用E*arnode（i）计算的
!欢迎指正   QQ 461068444 xug@fhdigz.com
finish
/cle
/nopr
keyw,pr_struc,1

/prep7
! 单位 KN ，M
*set,l1,4 !板长
*set,h1,4!板宽
*set,dis_x,0.002 !板与板之间X向间距，实际是0.2m，此处为后边的cpintf提供方便
*set,dis_y,0.005 !板与板之间Y向间距，实际是0.2m，此处为后边的cpintf提供方便
*set,h,0.4 !板厚
*set,k1,6e4!地基反应模量
*set,net,0.2 !网格0.2
*set,length,6.2 ! 履带车接地长度
*set,width,0.9 !履带车接地宽度
*set,pressure,2500/0.9/6.2

et,1,shell43
et,2,combin14
mp,ex,1,2.8e7!C40 混凝土弹性模量
mp,nuxy,1,0.2!泊松比
mp,dens,1,2.5   !密度，自身重力对后面的剪力，弯距是没有影响的，可能只影响弯沉
r,1,h
/title,heavy duty pavement
/rep
/ANG,1
/REP,FAST
k,1
k,2,l1
k,3,l1,h1
k,4,,h1
a,1,2,3,4
lesize,all,net
type,1
mat,1
real,1
amesh,all
agen,3,all,,,l1+dis_x,,,,0,,
allsel,all
agen,3,all,,,,l1+dis_y,,,0,,
!耦合
nsel,s,loc,x,l1,l1+dis_x
nplot
cpintf,ux,dis_x+0.001! 不加这个0.001，好像耦合的节点没有完全选到
cpintf,uy,dis_x+0.001
cpintf,uz,dis_x+0.001
allsel
nsel,s,loc,x,2*l1+dis_x,2*l1+2*dis_x
nplot
cpintf,ux,dis_x+0.001
cpintf,uy,dis_x+0.001
cpintf,uz,dis_x+0.001
allsel
nsel,s,loc,y,h1,h1+dis_y
nsel,u,loc,x,l1,l1+dis_x
nsel,u,loc,x,2*l1+dis_x,2*l1+2*dis_x
nplot
cpintf,ux,dis_y+0.001
cpintf,uy,dis_y+0.001
cpintf,uz,dis_y+0.001
allsel
nsel,s,loc,y,2*h1+dis_y,2*h1+2*dis_y
nsel,u,loc,x,l1,l1+dis_x
nsel,u,loc,x,2*l1+dis_x,2*l1+2*dis_x
nplot
cpintf,ux,dis_y+0.001
cpintf,uy,dis_y+0.001
cpintf,uz,dis_y+0.001
allsel
!设置弹簧实常数,产生弹簧单元
nsel,s,loc,z,0
*get,n1,node,,count
ngen,2,n1,all,,,,,-1
nsel,s,loc,z,-1
*get,n2,node,,num,min
allsel
*dim,rr,array,n1    !设定地基的弹性模量
*do,i,1,n1,1
rr(i)=arnode(i)
r,i+1,rr(i)*k1
type,2
mat,1
real,i+1
e,i,i+n2-1
*enddo
allsel
/pnum,node,0
/pnum,elem,0
/view,1,1,2,3
eplot    !!模型
save
!加荷载
/solu   !!荷载和约束
outres,all,all
antype,static
/psf,pres,norm,2,1
/rep
nsel,s,loc,z,-1
d,all,all
allsel
acel,,,10!重力加速度
time,1
nsubst,1
nsel,s,loc,x,l1+dis_x,l1+dis_x+width-0.01! 第一种情况
nsel,r,loc,y,0,length
esln,s,
eplot
sfe,all,2,pres,,pressure
allsel
lswrite,1
time,2
nsubst,1
esel,s,type,,1
sfedele,all,all,pres
allsel
nsel,s,loc,x,l1+dis_x+l1/2-width/2+0.01,l1+dis_x+l1/2+width/2-0.01
! 第二种情况
nsel,r,loc,y,0,length
esln,s
eplot
sfe,all,2,pres,,pressure
allsel
lswrite,2
time,3
nsubst,1
esel,s,type,,1
sfedele,all,all,pres
allsel
nsel,s,loc,x,2*l1+2*dis_x+l1/2-width/2+0.01,2*l1+2*dis_x+l1/2+width/2-0.01
! 第三种情况
nsel,r,loc,y,0,length
esln,s
eplot
sfe,all,2,pres,,pressure
allsel
lswrite,3
time,4
nsubst,1
esel,s,type,,1
sfedele,all,all,pres
allsel
nsel,s,loc,x,l1+dis_x,l1+dis_x+width-0.01! 第四种情况
nsel,r,loc,y,1.5*h1+dis_y-length/2+0.01,1.5*h1+dis_y+length/2-0.01
esln,s,
eplot
sfe,all,2,pres,,pressure
allsel
lswrite,4
time,5
nsubst,1
esel,s,type,,1
sfedele,all,all,pres
allsel
nsel,s,loc,x,l1+dis_x+l1/2-width/2+0.01,l1+dis_x+l1/2+width/2-0.01
! 第五种情况
nsel,r,loc,y,1.5*h1+dis_y-length/2+0.01,1.5*h1+dis_y+length/2-0.01
esln,s
eplot
sfe,all,2,pres,,pressure
allsel
lswrite,5
time,6
nsubst,1
esel,s,type,,1
sfedele,all,all,pres
allsel
nsel,s,loc,x,2*l1+2*dis_x+l1/2-width/2+0.01,2*l1+2*dis_x+l1/2+width/2-0.01
! 第六种情况
nsel,r,loc,y,1.5*h1+dis_y-length/2+0.01,1.5*h1+dis_y+length/2-0.01
esln,s
eplot
sfe,all,2,pres,,pressure
allsel
lswrite,6

lssolve,1,6,1
finish
/post1
!shell,mid
/VIEW,1,,,1
/RGB,INDEX,100,100,100, 0    !图像背景白色显示
/RGB,INDEX, 80, 80, 80,13   
/RGB,INDEX, 60, 60, 60,14   
/RGB,INDEX, 0, 0, 0,15
/plopts,info,1
/plopts,date,0
/rep
set,1
/title,mx
etable,mx_1,smisc,4
pletab,mx_1,avg
/image,save,mx_1,jpeg
/title,nx
etable,nx_1,smisc,7
pletab,nx_1,avg
/image,save,nx_1,jpeg

set,2
/title,mx
etable,mx_2,smisc,4
pletab,mx_2,avg
/image,save,mx_2,jpeg
/title,nx
etable,nx_2,smisc,7
pletab,nx_2,avg
/image,save,nx_2,jpeg
set,3
/title,mx
etable,mx_3,smisc,4
pletab,mx_3,avg
/image,save,mx_3,jpeg
/title,nx
etable,nx_3,smisc,7
pletab,nx_3,avg
/image,save,nx_3,jpeg
set,4
/title,mx
etable,mx_4,smisc,4
pletab,mx_4,avg
/image,save,mx_4,jpeg
/title,nx
etable,nx_4,smisc,7
pletab,nx_4,avg
/image,save,nx_4,jpeg

set,5
/title,mx
etable,mx_5,smisc,4
pletab,mx_5,avg
/image,save,mx_5,jpeg
/title,nx
etable,nx_5,smisc,7
pletab,nx_5,avg
/image,save,nx_5,jpeg

set,6
/title,mx
etable,mx_6,smisc,4
pletab,mx_6,avg
/image,save,mx_6,jpeg
/title,nx
etable,nx_6,smisc,7
pletab,nx_6,avg
/image,save,nx_6,jpeg

xugang2000 发表于 2008-2-20 10:52:44

这个是计算图示

xugang2000 发表于 2008-2-20 10:55:06

细部图

xugang2000 发表于 2008-2-20 11:06:15

付上一篇论文 ，大家多看看，我得思路和论文是一样的

xugang2000 发表于 2008-2-20 11:08:22

程序有乱码的话打开下边这个，pspad写的

yuefl 发表于 2008-2-25 17:25:32

下载了 谢谢分享 好好学习

mylang888 发表于 2008-4-12 15:25:08

楼主是个好人，我先下下来看看

xsrose376 发表于 2008-9-18 15:32:43

谢谢分享。

yequnshan 发表于 2009-7-2 15:04:40

学习学习谢谢
页: [1]
查看完整版本: 路面板的分析，考虑传力杆