引用回帖: 6樓 : Originally posted by
shgao20 at 2017-08-10 08:00:02
不輸入finish,也沒有執(zhí)行如 /post1 之類的命令轉(zhuǎn)換到其它模塊,就不會離開求解模塊。
要改變單元材料屬性,由非線性材料,變?yōu)榫性材料,有兩個辦法:
1 預(yù)先定義好所需的線性材料,然后在載荷步之間使用 ...
大神你好,我按照你提供的方法,在/solu模塊中,利用mpchg命令,改變了單元的材料屬性,但是在后處理,提取結(jié)果時,提示警告:Cumulative iteration 7 may have been solved using different model or boundary condition data than currently stored.POST1 results may be erroneous unless you RESUME from a Jobname.DB file for this substep. 導(dǎo)致無法提取結(jié)果,請問如何解決呢?
我做的是地鐵隧道開挖模擬,考慮地層損失,在求解時,需要轉(zhuǎn)換單元材料屬性,命令流如下:
FINISH
/CLEAR
/filename,sdkw
/prep7
ET,1,SOLID45 !用于模擬土體的三維實(shí)體單元
ET,2,shell181 !用于模擬襯砌的殼單元
r,1,0.35
!輸入土層信息
MP,EX,1,3e7 !輸入彈性模量(圍巖)
MP,PRXY,1,0.3 !輸入泊松比(圍巖)
MP,DENS,1,1830 !輸入密度(圍巖)
TB,DP,1 !采用DP準(zhǔn)則進(jìn)行彈塑性分析
TBdata,,22e4,12 !輸入凝聚力(圍巖)
!輸入間隙材料信息
MP,EX,2,2e5 !輸入彈性模量(圍巖)
MP,PRXY,2,0.3 !輸入泊松比(圍巖)
MP,DENS,2,1830 !輸入密度(圍巖)
TB,DP,2 !采用DP準(zhǔn)則進(jìn)行彈塑性分析
TBdata,,22e4,12 !輸入凝聚力(圍巖)
!輸入襯砌信息
MP,EX,3,3.45e10 !輸入彈性模量
MP,PRXY,3,0.2 !輸入泊松比
MP,DENS,3,2500 !輸入密度
!輸入隧道開挖土體
MP,EX,4,3e7 !輸入彈性模量(圍巖)
MP,PRXY,4,0.3 !輸入泊松比(圍巖)
MP,DENS,4,1830 !輸入密度(圍巖)
TB,DP,4 !采用DP準(zhǔn)則進(jìn)行彈塑性分析
TBdata,,22e4,12 !輸入凝聚力(圍巖)
!輸入土層損失信息
MP,EX,5,3e7 !輸入彈性模量(圍巖)
MP,PRXY,5,0.3 !輸入泊松比(圍巖)
MP,DENS,5,1830 !輸入密度(圍巖)
TB,DP,5 !采用DP準(zhǔn)則進(jìn)行彈塑性分析
TBdata,,22e4,12 !輸入凝聚力(圍巖)
!!!!幾何建模
csys,0
wpcsys,,0
BLOCK,0,100,0,50,0,50 !建立圍巖
allsel
cyl4,50,35,3,,,,50
cyl4,50,35,3.05,,,,50
vptn,all
wpoff,,,25
vsbw,all
wpcsys,-1
wpoff,50,35
wprota,,,-90
vsbw,all
wpoff,,,6
vsbw,all
wpoff,,,-12
vsbw,all
wprota,,90
vsbw,all
wpoff,,,6
vsbw,all
wpoff,,,-12
vsbw,all
csys,0
wpcsys,,0
allsel
vsel,s,loc,y,33.5,36.5
vsel,r,loc,x,48.5,51.5
cm,dwtt,volu
vatt,4,,1
allsel
csys,0
wpcsys,,0
vsel,s,loc,y,32,38
vsel,r,loc,x,37,53
cmsel,u,dwtt
cm,dcss,volu
vatt,5,,1
csys,0
wpcsys,,0
local,12,1,50,35
csys,12
wpcsys,,12
asel,s,loc,x,3
cm,cq,area
aatt,3,1,2
allsel
csys,0
wpcsys,,0
vsel,all
cmsel,u,dwtt
cmsel,u,dcss
vatt,1,,1
allsel
LSEL,S,LENGTH,,3
LESIZE,ALL,,,4,,,,,1
LSEL,S,LENGTH,,0.05
LESIZE,ALL,,,1,,,,,1
*SET,PI,ACOS(-1)
LSEL,S,LENGTH,,PI*3/2
LESIZE,ALL,,,8,,,,,1
allsel
vsel,all
esize,2
mshape,0
vsweep,all
allsel
mshkey,1
amesh,cq
!!!!邊界條件
allsel
csys,0
wpcsys,,0
ACEL,,9.8
NSEL,S,LOC,Y,0
D,ALL,uy
nsel,s,loc,x,0
d,all,ux
nsel,s,loc,x,100
d,all,ux
nsel,s,loc,z,0
d,all,uz
nsel,s,loc,z,50
d,all,uz
allsel,all
!nummrg,all,,,,low !合并重復(fù)各元素并保留低編號號碼
!numcmp,all !壓縮各元素編號號碼
allsel,all
finish
!!!!進(jìn)入求解器
/solu
antype,static
deltim,0.1,0.1,0.2 !時間步設(shè)置,總長0.1,最小0.1,最大0.2
autots,on !使用自動時間步
pred,on !打開時間步長預(yù)測器
lnsrch,on !打開線性搜索
nlgeom,on !打開大位移效果
nropt,full !設(shè)定牛頓-拉普森選項(xiàng)
cnvtol,f,,0.05,2,0.5 !設(shè)定力收斂條件
Save !保存數(shù)據(jù)庫
esel,s,mat,,3 !選擇襯砌單元
ekill,all !殺死襯砌單元,即在隧道未修建前的自重應(yīng)力場中不存在襯砌
esel,all !選擇所有元素
esel,s,live !選擇活單元,即所有圍巖實(shí)體單元
nsle,s !選擇生單元上的節(jié)點(diǎn)
nsel,invert !反向選擇,即選擇了死單元上的節(jié)點(diǎn)
d,all,all !將死單元上的節(jié)點(diǎn)約束所有位移,使其不參與矩陣運(yùn)算
nsel,all !選擇所有節(jié)點(diǎn)
esel,all !選擇所有單元
Save !保存數(shù)據(jù)庫
iswrite,on !打開初應(yīng)力文件生成開關(guān)
solve !進(jìn)行自重地應(yīng)力場模擬計算,第一荷載步
Save
/post1
/solu
isfile,read,sdkw,ist,,2 !讀入初應(yīng)力文件,作為荷載施加
!isfile,list !查看施加的初應(yīng)力荷載
solve !進(jìn)行地應(yīng)力平衡計算(此時的荷載為重力荷載和初應(yīng)力荷載)
!!!開挖土體
!!!以下進(jìn)行土體的開挖操作,先選擇每次開挖的圍巖單元,然后將其賦予“死屬性”
!第一個荷載開挖步
esel,s,mat,,4
nsle,s
nsel,r,loc,z,0,25
esln,r,1
ekill,all
!以下進(jìn)行支護(hù)結(jié)構(gòu)的施加操作,先選擇支護(hù)結(jié)構(gòu)殼單元,然后將其賦予“生屬性”
esel,s,mat,,3
nsle,s
nsel,r,loc,z,0,25
esln,r,1
ealive,all
nsle,s
ddele,all,all !將定義的約束條件刪除
esel,s,mat,,5
nsle,s
nsel,r,loc,z,0,25
esln,r,1
mpchg,2,all !將材料號為5的單元變?yōu)椴牧咸?
!選擇生單元,即包括支護(hù)結(jié)構(gòu)殼單元和未開挖部分圍巖實(shí)體單元
esel,all
esel,s,live
nsle,s
!反向選擇,并將死單元上的節(jié)點(diǎn)約束所有自由度
nsel,invert
d,all,all
nsel,all !選擇所有節(jié)點(diǎn)
esel,all !選擇所有單元
solve
!第二個荷載開挖步
esel,s,mat,,4
nsle,s
nsel,r,loc,z,25,50
esln,r,1
ekill,all
!以下進(jìn)行支護(hù)結(jié)構(gòu)的施加操作,先選擇支護(hù)結(jié)構(gòu)殼單元,然后將其賦予“生屬性”
esel,s,mat,,3
nsle,s
nsel,r,loc,z,25,50
esln,r,1
ealive,all
nsle,s
ddele,all,all !將定義的約束條件刪除
esel,s,mat,,5
nsle,s
nsel,r,loc,z,25,50
esln,r,1
mpchg,2,all
!選擇生單元,即包括支護(hù)結(jié)構(gòu)殼單元和未開挖部分圍巖實(shí)體單元
esel,all
esel,s,live
nsle,s
!反向選擇,并將死單元上的節(jié)點(diǎn)約束所有自由度
nsel,invert
d,all,all
nsel,all !選擇所有節(jié)點(diǎn)
esel,all !選擇所有單元
solve
/post1
/DEVICE,VECTOR,0
SET,1,LAST,1
path,dbcj,2
ppath,1,,50,50,0
ppath,2,,50,50,50
pdef,dbcj,u,y
prpath,dbcj !提取隧道頂部地表豎向位移