| 1 | 1/1 | 返回列表 |
| 查看: 928 | 回復(fù): 0 | ||
| 【懸賞金幣】回答本帖問題,作者小劉的學(xué)習(xí)將贈送您 5 個金幣 | ||
小劉的學(xué)習(xí)新蟲 (初入文壇)
|
[求助]
在FLAC3D5.0中,加入自定義本構(gòu)模型后,動力計算出錯
|
|
|
在FLAC3D中,樁土模型,土使用的自定義的循環(huán)軟化模型,進行位移加載,然后出現(xiàn)“fool==dfool” 錯誤,這是模型哪里出錯了嗎? 命令流如下: 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ā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料與化工,270求調(diào)劑 +16 | YXCT 2026-03-01 | 18/900 |
|
|---|---|---|---|---|
|
[考研] 材料與化工,291分,求調(diào)劑 +8 | 咕嚕咕嚕123123 2026-03-04 | 10/500 |
|
|
[考研] 307求調(diào)劑 +4 | xuyuan111 2026-03-04 | 5/250 |
|
|
[考研] 0855機械工程324求調(diào)劑 +3 | 鄭SK 2026-03-02 | 7/350 |
|
|
[考研] 材料專碩346求調(diào)劑 +3 | 旺一下 2026-03-04 | 3/150 |
|
|
[考研] 一志愿314求調(diào)劑 +7 | 202111120625 2026-03-03 | 7/350 |
|
|
[考研] 347求調(diào)劑 +5 | 啊歐歐歐 2026-03-03 | 5/250 |
|
|
[考研] 281求調(diào)劑 +3 | Y?l?h 2026-03-04 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | G212 2026-03-03 | 5/250 |
|
|
[考研] 材料工程269求調(diào)劑 +7 | 白刺玫 2026-03-02 | 7/350 |
|
|
[考研] 一志愿華中科技大學(xué),化學(xué)專業(yè)344分,求調(diào)劑 +6 | 邢xing1 2026-03-02 | 6/300 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +15 | February_Feb 2026-03-01 | 17/850 |
|
|
[考研] 284求調(diào)劑 +6 | 天下熯 2026-03-02 | 6/300 |
|
|
[考研]
材料工程專碩283求調(diào)劑
5+8
|
,。 2026-03-02 | 10/500 |
|
|
[考研] 298求調(diào)劑 +10 | 人間唯你是清歡 2026-02-28 | 14/700 |
|
|
[考研] 求調(diào)劑 +11 | yunziaaaaa 2026-03-01 | 13/650 |
|
|
[考研] 材料工程274求調(diào)劑 +5 | Lilithan 2026-03-01 | 5/250 |
|
|
[考研] 275求調(diào)劑 +3 | L-xin? 2026-03-01 | 6/300 |
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-02 | 4/200 |
|
|
[考研] 317一志愿華南理工電氣工程求調(diào)劑 +6 | Soliloquy_Q 2026-02-28 | 11/550 |
|