| 24小時(shí)熱門版塊排行榜 |
| 1 | 1/1 | 返回列表 |
| 查看: 933 | 回復(fù): 0 | |||
| 【懸賞金幣】回答本帖問(wèn)題,作者小劉的學(xué)習(xí)將贈(zèng)送您 5 個(gè)金幣 | |||
小劉的學(xué)習(xí)新蟲(chóng) (初入文壇)
|
[求助]
在FLAC3D5.0中,加入自定義本構(gòu)模型后,動(dòng)力計(jì)算出錯(cuò)
|
||
|
在FLAC3D中,樁土模型,土使用的自定義的循環(huán)軟化模型,進(jìn)行位移加載,然后出現(xiàn)“fool==dfool” 錯(cuò)誤,這是模型哪里出錯(cuò)了嗎? 命令流如下: new config dyn set dyn off config cppudm model load UWmodel.dll gen zone radcyl p0 0 0 0 p1 20 0 0 p2 0 0 -6 p3 0 15 0 & p4 20 0 -6 p5 0 15 -6 p6 20 15 0 p7 20 15 -6 & p8 0.15 0 0 p9 0 0.15 0 p10 0.15 0 -6 p11 0 0.15 -6 & size 3 15 15 15 ratio 1 1 1 1.15 group clay1 range z 0 -1.3 group clay2 range z -1.3 -6 gen zone reflect dd 270 dip 90 gen zone reflect dd 0 dip 90 interface 1 face range cylinder end1 (0,0,0) end2 (0 ,0 ,-6.1) radius 0.16 & cylinder end1 (0,0,0) end2 (0,0,-6.1) radius 0.14 not gen zone cyl p0 0 0 7.2 p1 0.15 0 7.2 p2 0 0 1.2 p3 0 0.15 7.2 & p4 0.15 0 1.2 p5 0 0.15 1.2 size 2 15 15 gen zone cyl p0 0 0 7.8 p1 0.15 0 7.8 p2 0 0 7.8 p3 0 0.15 7.8 & p4 0.15 0 7.2 p5 0 0.15 7.2 size 2 10 15 gen zone reflect dd 270 dip 90 range z 1.2 7.8 gen zone reflect dd 0 dip 90 range z 1.2 7.8 group pile range z 1.2 7.8 ini z add -7.2 range group pile range name=clay1 group clay1 range name=clay2 group clay2 range name=pile group pile model mohr range clay1 prop shear 2.04e7 bulk 4.42e7 range clay1 prop coh 28200 fric 2.9 range clay1 model mohr range clay2 prop shear 2.04e7 bulk 4.42e7 range clay2 prop coh 28200 fric 2.9 range clay2 model elastic range pile prop bulk 4.42e7 shear 2.04e7 range pile model mech null range z -0.1 0.61 interface 1 prop kn 1e7 ks 1e7 fric 1.5 coh 14100 fix x range x -20.1 -19.9 fix x range x 19.9 20.1 fix y range y -15.1 -14.9 fix y range y 14.9 15.1 fix x y z range z -6.1 -5.9 set grav 0 0 -10 water density 1000 water table origin 0,0,-1.3 normal 0 0 -1 ini dens 1570 range clay1 ini dens 1570 range clay2 ini dens 1570 range pile solve model elastic range pile prop bulk 1.09e8 shear 8.17e7 range pile ini dens 2500 range pile solve model UWmodel range clay1 prop Gmax 20400000 p0 45510 mm 0.5 PR 0.30 CC 28200 Fai 2.9 bb 432 nn 3.44 Rf 0.9 range clay1 prop TauRef 0.8 A0 0.05 Stpp 1 range clay1 prop dens 1570 range clay1 model UWmodel range clay2 prop Gmax 20400000 p0 45510 mm 0.5 PR 0.30 CC 28200 Fai 2.9 bb 432 nn 3.44 Rf 0.9 range clay2 prop TauRef 0.8 A0 0.05 Stpp 1 range clay2 prop dens 1570 range clay2 solve set dyn on ini state 0 ini xv 0 yv 0 zv 0 xdisp 0 ydisp 0 zdisp 0 sel beam id=1 begin=(0,0,0.6) end=(0,0, 0) nseg=10 sel beam id=1 begin=(0,0,0) end=(0,0,-6) nseg=15 sel beam id=1 prop emod=3.8964e10 nu=0.2 dens=2500 & xcarea=7.1e-2 xcj=3.976e-4 xciy=3.976e-4 xciz=0 & pmoment=51200 sel node fix x y z xr yr zr range id=12 def xf_p1 local ad = top_head local xftot = 0.0 loop while ad # null local gp_pnt = mem(ad+1) local xf = gp_xfunbal(gp_pnt) xftot = xftot + xf ad = mem(ad) endloop xf_p1 = xftot end fix x range z 0.55 0.61 group pile def init_ramp global udapp=0 global udmax=2.5e-2 end @init_ramp def ramp while_stepping if dytime < 1 then udapp=udmax else if dytime>1 if dytime<3 then udapp=-udmax else if dytime>3 if dytime<5 then udapp=udmax else if dytime>5 if dytime<7 then udapp=-udmax else if dytime>7 if dytime<8 then udapp=udmax endif endif endif endif endif endif endif endif endif local ad = top_head loop while ad # null local gp_pnt = mem(ad+1) gp_xvel(gp_pnt) = udapp ad = mem(ad) endloop end def find_add global top_head = null loop foreach local gp_pnt gp_list if gp_zpos(gp_pnt) > 0.59 then local mem_head = get_mem(2) mem(mem_head) = top_head mem(mem_head+1) = gp_pnt top_head = mem_head endif endloop end @find_add hist add fish @xf_p1 hist add gp xdisp 0 0 0.6 hist add gp xvel 0 0 0.6 set mech damp comb set dyn multi on solve age 28 |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[文學(xué)芳草園] 伙伴們,祝我生日快樂(lè)吧 +3 | myrtle 2026-03-10 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 279求調(diào)劑 +3 | 莫xiao 2026-03-10 | 4/200 |
|
|
[考研] 哈工大材料324求調(diào)劑 +6 | 閆旭東 2026-03-10 | 8/400 |
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲(chóng)085600 2026-03-09 | 11/550 |
|
|
[考研] 313分生物學(xué)求調(diào)劑 +5 | Yyt楊1 2026-03-09 | 7/350 |
|
|
[考研] A區(qū)一本交叉課題組,低分調(diào)劑,招收機(jī)械電子信息通信等交叉方向 +56 | lisimayy 2026-03-04 | 73/3650 |
|
|
[考研] 327分求調(diào)劑086 +4 | 西紅柿?小帥 2026-03-09 | 7/350 |
|
|
[考研] 0817化學(xué)工程319求調(diào)劑 +7 | lv945 2026-03-08 | 9/450 |
|
|
[考博] 2026博士申請(qǐng) +6 | 起泡酒 2026-03-08 | 6/300 |
|
|
[考研] 334求調(diào)劑 +5 | Trying] 2026-03-07 | 5/250 |
|
|
[考研] 0856求調(diào)劑 +3 | squirtle11 2026-03-07 | 3/150 |
|
|
[考研] 0817找調(diào)劑 +6 | kk扛 2026-03-08 | 6/300 |
|
|
[考研] 安徽農(nóng)業(yè)大學(xué)材料與化學(xué)學(xué)院0856材料招收調(diào)劑 +3 | akakk47 2026-03-05 | 9/450 |
|
|
[考研] 288求調(diào)劑085600材料與化工 +13 | Daunrin 2026-03-07 | 15/750 |
|
|
[考研] 334求調(diào)劑 +8 | Trying] 2026-03-06 | 8/400 |
|
|
[考研] 085701環(huán)境工程專業(yè),初試305,均過(guò)國(guó)家A區(qū)線 +7 | 卡卡來(lái)了@ 2026-03-07 | 8/400 |
|
|
[考研] 化學(xué)290求調(diào)劑 +3 | 兩顆 西柚 2026-03-07 | 4/200 |
|
|
[考研] 一志愿武漢理工085601,初試301分,請(qǐng)問(wèn)能調(diào)劑到湖北嗎 +3 | 肖yang 2026-03-06 | 3/150 |
|
|
[考研] 2026調(diào)劑】考試A區(qū)0703化學(xué)類323分 誠(chéng)求接收 +3 | 卷柏卷柏 2026-03-05 | 4/200 |
|
|
[考研] 不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過(guò)來(lái) +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|