找回密码
 注册
Simdroid-非首页
查看: 257|回复: 9

[问题讨论] 请教一个移动荷载作用在三维路基上的问题!

[复制链接]
发表于 2007-12-21 17:07:14 | 显示全部楼层 |阅读模式 来自 湖南长沙
我自己先建立了一个2D路基横断面(plane 42 单元),对其进行划分网格后把它转换成3D模型(solid 45) ,然后沿路纵向从一端到另一端作用移动的面荷载,(DO循环实现),前处理完后利用/POST26分别观察作用路径中间处某个节点(路面表面或垂直下面不同深度处)的位移和应力时程曲线,发现规律和别人做的大致相同,但是别人相似问题得到的曲线很平滑,而我的总是有波动,而且没有平滑过度感。我改了不同参数(包括网格划分加细,改变不同速度或者指定KBC值等)还是一样的,变化不大,请问我的问题到底出在哪里?是不是出现在什么命令上面?请各位指教!谢谢!下面是我的命令流(个人感觉前面建模应该没问题):
finish
/clear
/filname,luji3d
/config,nres,2000000    !设置最大子步数
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!土的参数!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  l=5   !路纵向长
lmk=10   !路面宽度
ljk=20   !路基宽度
ljg=30   !路基高度
  i=1.5  !坡度
mc=0.18  !面层
jc=0.38  !基层
ny=0.4   !碾压砼
sbp=0.25  !手摆片石
ss=0.4   !碎石垫层
snt=0.2   !掺4%水泥土
ldt=0.7   !路堤土
nt=1.8   !粘土
yn=9     !淤泥
yjs=1.8   !亚粘土夹砂
ynt=5.7   !亚粘土
yjs0=ljg-nt-yn-yjs-ynt  !亚粘土夹砂

/PREP7

ET,1,PLANE42
ET,2,solid45

            
MP,EX,1,1e9        !面层材料属性 (材料1  沥青混凝土)
MP,PRXY,1,0.25
MP,DENS,1,2500

MP,EX,2,7e8        !基层材料属性 (材料2  水泥稳定碎石基层)
MP,PRXY,2,0.3
MP,DENS,2,2240

MP,EX,3,3e8        !碾压砼材料属性 (材料3)
MP,PRXY,3,0.3
MP,DENS,3,2100

MP,EX,4,1.5e8      !手摆片石材料属性 (材料4)
MP,PRXY,4,0.3
MP,DENS,4,2020
MP,EX,5,5.9e7      !碎石垫层材料属性 (材料5)
MP,PRXY,5,0.3
MP,DENS,5,1840
MP,EX,6,5.2e7      !掺4%水泥土材料属性 (材料6)
MP,PRXY,6,0.3
MP,DENS,6,1910

MP,EX,7,5e7        !路堤土层材料属性 (材料7)
MP,PRXY,7,0.3
MP,DENS,7,1850
MP,MU,7,0.25       !摩擦系数
TB,DP,7,,,  
TBMODIF,1,1,5000   !粘聚力
TBMODIF,1,2,30     !内摩擦角
TBMODIF,1,3,0      !膨胀角
MP,EX,8,3.8e7      !粘土材料属性 (材料8)
MP,PRXY,8,0.3
MP,DENS,8,1930
MP,MU,8,0.25       !摩擦系数
TB,DP,8,,,  
TBMODIF,1,1,25000  !粘聚力
TBMODIF,1,2,15     !内摩擦角
TBMODIF,1,3,0      !膨胀角

MP,EX,9,3.5e7      !淤泥材料属性 (材料9)
MP,PRXY,9,0.3
MP,DENS,9,1720
MP,MU,9,0.25       !摩擦系数
TB,DP,9,,,  
TBMODIF,1,1,25000  !粘聚力
TBMODIF,1,2,10     !内摩擦角
TBMODIF,1,3,0      !膨胀角
MP,EX,10,4.5e7     !亚粘土夹砂材料属性 (材料10)
MP,PRXY,10,0.3
MP,DENS,10,1750
MP,MU,10,0.25      !摩擦系数
TB,DP,10,,,  
TBMODIF,1,1,10000  !粘聚力
TBMODIF,1,2,25     !内摩擦角
TBMODIF,1,3,0      !膨胀角
MP,EX,11,2e7       !亚粘土材料属性 (材料11)
MP,PRXY,11,0.3
MP,DENS,11,1780
MP,MU,11,0.25      !摩擦系数
TB,DP,11,,,  
TBMODIF,1,1,20000  !粘聚力
TBMODIF,1,2,25     !内摩擦角
TBMODIF,1,3,0      !膨胀角

MP,EX,12,3.5e7     !亚粘土夹砂材料属性 (材料12)
MP,PRXY,12,0.3
MP,DENS,12,1750
MP,MU,12,0.25      !摩擦系数
TB,DP,12,,,  
TBMODIF,1,1,10000  !粘聚力
TBMODIF,1,2,25     !内摩擦角
TBMODIF,1,3,0      !膨胀角

x1=lmk+mc*i+jc*i
x2=x1+ny*i+sbp*i
x3=x2+ss*i+snt*i
y1=yjs0+ynt+yjs
y2=y1+yn+nt
y3=y2+ldt+snt
y4=y3+ss+sbp
y5=y4+ny+jc
k,1,0,0,0              !!!!!!!!!!建立模型!!!!!!!!!!!
k,2,x3+ldt*i,0,0
k,3,ljk,0,0
k,4,ljk,yjs0,0
k,5,x3+ldt*i,yjs0,0
k,6,0,yjs0,0
k,7,0,yjs0+ynt,0
k,8,x3+ldt*i,yjs0+ynt,0
k,9,ljk,yjs0+ynt,0
k,10,ljk,y1,0
k,11,x3+ldt*i,y1,0
k,12,0,y1,0
k,13,0,y1+yn,0
k,14,x3+ldt*i,y1+yn,0
k,15,ljk,y1+yn,0
k,16,ljk,y2,0
k,17,x3+ldt*i,y2,0   !边坡点
k,18,0,y2,0
k,19,0,y2+ldt,0
k,20,x3,y2+ldt,0  !边坡点
k,21,x2+ss*i,y3,0  !边坡点
k,22,0,y3,0
k,23,0,y3+ss,0
k,24,x2,y3+ss,0  !边坡点
k,25,x1+ny*i,y4,0 !边坡点
k,26,0,y4,0
k,27,0,y4+ny,0
k,28,x1,y4+ny,0  !边坡点
k,29,lmk+mc*i,y5,0    !边坡点
k,30,0,y5,0
k,31,0,y5+mc,0
k,32,lmk,y5+mc,0   !边坡点
k,33,0,0,l
k=1
*do,j,1,31
k=j+1
l,j,k
*enddo
l,30,27
l,26,23
l,22,19
l,18,13
l,12,7
l,6,1
l,32,29
l,28,25
l,24,21
l,20,17
l,17,14
l,14,11
l,11,8
l,8,5
l,5,2
l,15,10
l,9,4
l,1,33

al,1,46,5,37
al,2,3,4,46
al,5,45,7,6
al,4,48,8,45
al,7,44,11,36
al,8,9,10,44
al,11,43,13,12
al,10,47,14,43
al,13,42,17,35
al,14,15,16,42
al,17,41,19,18
al,19,20,21,34
al,21,40,23,22
al,23,24,25,33
al,25,39,27,26
al,27,28,29,32
al,29,38,31,30
allsel
lsel,s,,,37
lsel,a,,,46
lsel,a,,,3
lesize,all,,,5 !5等分
allsel
lsel,s,,,6
lsel,a,,,45
lsel,a,,,48
lesize,all,,,3 !3等分
allsel
lsel,s,,,36
lsel,a,,,44
lsel,a,,,9
lesize,all,,,2 !2等分
allsel
lsel,s,,,12
lsel,a,,,43
lsel,a,,,47
lesize,all,,,5 !5等分
allsel
lsel,s,,,35
lsel,a,,,42
lsel,a,,,15
lesize,all,,,2 !2等分
allsel
lsel,s,,,2
lsel,a,,,4
lsel,a,,,8
lsel,a,,,10
lsel,a,,,14
lsel,a,,,16
lesize,all,,,5!5等分
allsel
lsel,s,,,1
lsel,a,,,5
lsel,a,,,7
lsel,a,,,11
lsel,a,,,13
lsel,a,,,17
lsel,a,,,19
lsel,a,,,21
lsel,a,,,23
lsel,a,,,25
lsel,a,,,27
lsel,a,,,29
lsel,a,,,31
lesize,all,,,20 !20等分
allsel
lsel,s,,,18
lsel,a,,,41
lesize,all,,,1 !1等分
allsel
lsel,s,,,34
lsel,a,,,20
lesize,all,,,1 !1等分
allsel
lsel,s,,,22
lsel,a,,,40
lesize,all,,,1 !1等分
allsel
lsel,s,,,33
lsel,a,,,24
lesize,all,,,1 !1等分
allsel
lsel,s,,,26
lsel,a,,,39
lesize,all,,,1 !1等分
allsel
lsel,s,,,32
lsel,a,,,28
lesize,all,,,1 !1等分
allsel
lsel,s,,,30
lsel,a,,,38
lesize,all,,,1 !1等分
allsel

mshape,0,2d              !划分网格              
mshkey,1                 !映射网格
type,1                  
mat,1                    !按材料属性划分网格,从上到下依次从1-12种材料属性,顺序和前面一致               
amesh,17
mat,2                     
amesh,16
mat,3                     
amesh,15                  
mat,4                     
amesh,14                  
mat,5                     
amesh,13
mat,6                     
amesh,12
mat,7                     
amesh,11
mat,8                     
amesh,9,10,1
mat,9                     
amesh,7,8,1
mat,10                     
amesh,5,6,1
mat,11                     
amesh,3,4,1        
mat,12                     
amesh,1,2,1

ARSYM,X,all, , , ,0,0      !对称
EPLOT
nummrg,all                 !压缩和合并
numcmp,all
lsel,s,,,49
lesize,all,,,10
asel,s,,,1,34,1
VDRAG,all, , , , , , 49    !映射到3D模型
aclear,1,34,1              !删除原来的平面
allsel
vsel,s,,,34
vsel,a,,,17
eslv,r                     !使单元与体相关
EMODIF,all,mat,1           !修改单元属性(由2D映射到3D单元)1为材料编号
allsel

vsel,s,,,33
vsel,a,,,16
eslv,r
EMODIF,all,mat,2
allsel
vsel,s,,,32
vsel,a,,,15
eslv,r
EMODIF,all,mat,3
allsel
vsel,s,,,31
vsel,a,,,14
eslv,r
EMODIF,all,mat,4
allsel
vsel,s,,,30
vsel,a,,,13
eslv,r
EMODIF,all,mat,5
allsel
vsel,s,,,29
vsel,a,,,12
eslv,r
EMODIF,all,mat,6
allsel
vsel,s,,,28
vsel,a,,,11
eslv,r
EMODIF,all,mat,7
allsel
vsel,s,,,27
vsel,a,,,26
vsel,a,,,9
vsel,a,,,10
eslv,r
EMODIF,all,mat,8
allsel

vsel,s,,,25
vsel,a,,,24
vsel,a,,,7
vsel,a,,,8
eslv,r
EMODIF,all,mat,9
allsel
vsel,s,,,23
vsel,a,,,22
vsel,a,,,5
vsel,a,,,6
eslv,r
EMODIF,all,mat,10
allsel
vsel,s,,,21
vsel,a,,,20
vsel,a,,,3
vsel,a,,,4
eslv,r
EMODIF,all,mat,11
allsel
vsel,s,,,19
vsel,a,,,18
vsel,a,,,1
vsel,a,,,2
eslv,r
EMODIF,all,mat,12
allsel
asel,s,loc,y,0             !边界条件
da,all,all
allsel
asel,s,loc,x,ljk
asel,a,loc,x,-ljk
da,all,ux
da,all,uz
allsel
asel,s,loc,z,0
asel,a,loc,z,l
da,all,uz
da,all,ux
allsel
save
allsel
/solu
antype,trans
!autots,on
deltim=0.0015
tim=0
timint,on
*do,i,0,5,0.5(移动次数)
tim=tim+0.03
time,tim
!nsel,s,loc,y,32.51
!nsel,r,loc,x,0.5,1   !大车后轮1-1
!nsel,r,loc,z,i-4.5,i-4
!cm,hl1,nodes
!sf,hl1,pres,150000
!allsel
!nsel,s,loc,y,32.51
!nsel,r,loc,x,2,2.5   !大车后轮1-2
!nsel,r,loc,z,i-4.5,i-4
!cm,hl2,nodes
!sf,hl2,pres,150000
!allsel
!nsel,s,loc,y,32.51
!nsel,r,loc,x,0.5,1   !大车后轮2-1
!nsel,r,loc,z,i-3,i-2.5
!cm,hl3,nodes
!sf,hl3,pres,150000
!allsel
!nsel,s,loc,y,32.51
!nsel,r,loc,x,2,2.5   !大车后轮2-2
!nsel,r,loc,z,i-3,i-2.5
!cm,hl4,nodes
!sf,hl4,pres,150000
!allsel

nsel,s,loc,y,32.51
nsel,r,loc,x,0,0.5   !大车前轮1-1
nsel,r,loc,z,i-0.5,i
cm,ql1,nodes
sf,ql1,pres,100000
allsel
!nsel,s,loc,y,32.51
!nsel,r,loc,x,2,3   !大车前轮1-2
!nsel,r,loc,z,i-0.25,i
!cm,ql2,nodes
!sf,ql2,pres,10000
!allsel
nsubst,5
outres,all,all
solve
SFDELE,all,all
*enddo
finish
!注记:网格细分为0.5m*0.5m一个单元,施加大车面荷载200KN,前轮100000Pa,后轮150000Pa
!轮间距2m  轴距1.5m
!速度作用时间为0.03s(速度为60m/s),
!测点取路段中部,荷载作用正中间下面不同深度,
!l=5 路纵向长  lmk=10   !路面宽度  ljk=20   !路基宽度  ljg=30   !路基高度
后面附上我的图!
说明:下面的应力和位移全部是竖向的,即移动荷载作用方向。

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有账号?注册

×
 楼主| 发表于 2007-12-21 17:09:15 | 显示全部楼层 来自 湖南长沙
Simdroid开发平台
请教各位:为什么我的时程曲线波动那么大?
回复 不支持

使用道具 举报

 楼主| 发表于 2007-12-26 21:43:52 | 显示全部楼层 来自 湖南长沙
怎么没有人回答我的问题啊?
版主能帮帮忙么?
回复 不支持

使用道具 举报

 楼主| 发表于 2007-12-27 11:28:37 | 显示全部楼层 来自 湖南长沙
自己顶一下!望各位大人帮帮忙!谢谢
回复 不支持

使用道具 举报

 楼主| 发表于 2007-12-28 17:03:51 | 显示全部楼层 来自 湖南长沙
我再顶一下
回复 不支持

使用道具 举报

 楼主| 发表于 2008-1-5 19:23:17 | 显示全部楼层 来自 湖南长沙
有没有人能帮帮忙
回复 不支持

使用道具 举报

 楼主| 发表于 2008-1-6 20:45:38 | 显示全部楼层 来自 湖南长沙
^^^^^^^^^^^^
回复 不支持

使用道具 举报

发表于 2008-1-6 21:21:39 | 显示全部楼层 来自 云南昆明
可能是边界约束不足吧!!!
回复 不支持

使用道具 举报

发表于 2009-4-22 17:59:34 | 显示全部楼层 来自 陕西西安
我也想知道。。。
回复 不支持

使用道具 举报

发表于 2016-6-3 19:48:19 | 显示全部楼层 来自 河北石家庄
不知道问题有没有解决?
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-29 06:41 , Processed in 0.055527 second(s), 13 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表