我自己先建立了一个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 !路基高度
后面附上我的图!
说明:下面的应力和位移全部是竖向的,即移动荷载作用方向。 |