abaqus計算砂土的案例
PFC休止角（堆積角）計算的滾筒方法——砂土顆粒
休止角的計算常用的有自然下落，然是這個方法受限制于顆粒數(shù)目，而且隨機性較大。
而滾筒法在低顆粒數(shù)目情況下也能有比較好的模擬效果。
這里給出砂土顆粒的滾筒模擬結(jié)果。
首先在生成圓形容器并且生成顆粒，加重力平衡
new
domain extent -1 1
wall generate circle position 0 0 radius 1
ball distribute porosity 0.28 radius 0.006 0.009 range annulus center 0 0 radius 0 [1-0.009]...
plane o 0 0 dip 0 belowcmat default model linear method deformability emod 100e6 kratio 1.74 property fric 0.5
ball attribute density 2.7e3 damp 0.7
cycle 2000 calm 50
set gravity 9.8
solve
save sample
以上為得到初始狀態(tài)所需的代碼。
然后簡單的給圓形容器一個滾動速度便可以了，注意墻和球之間的摩擦效率要很大，這樣模擬效率高。
展開 ABAQUS CEL (例7) 3D模擬巖土靜力觸探(CPT)貫入砂土 ￥66.67
ABAQUS CEL (例7) 3D模擬巖土靜力觸探(CPT)貫入砂土
一、模型背景
1) 3D模擬CPT貫入砂土全過程，為土體大變形模擬；
2) CPT錐角為60度，貫入速度為0.02m/s，貫入深度為10米;
3) 土本構(gòu)采用摩爾庫倫模型，以模擬CPT貫入砂土時土的應(yīng)力應(yīng)變行為。
二、模型的建立
利用軸對稱性，土層和CPT都只建立90度模型：
模型全貌
CPT的錐尖
CPT與土層初始位置
三、模擬的效果
初始狀態(tài)
CPT貫入土層中
CPT貫入時的土體應(yīng)變分布圖
CPT貫入時的土體應(yīng)力分布圖
CPT貫入時的砂土速度場流動分布圖
展開 ABAQUS CEL (例2) 海上開口樁貫入砂土 (附完整的input file) ￥66.67
ABAQUS CEL (例2) 海上開口樁貫入砂土 (附完整的input file)
背景
該input file模擬的是開口鋼樁在砂土中的貫入過程，樁相較于土的剛度大很多，因而采用解析剛體來模擬開口鋼管樁，土的本構(gòu)采用摩爾庫倫本構(gòu)模型。
所建模型
模擬結(jié)果
樁貫入土中時，土體水平應(yīng)力分布圖
樁貫入土中時，土體豎向應(yīng)力分布圖
樁貫入土中時，土體應(yīng)變分布圖
樁貫入土中時，樁內(nèi)外壁土體的流動矢量圖（土體速度分布圖）
在Abaqus軟件中計算礫石床對混凝土擋土墻的側(cè)壓力
在Abaqus軟件中計算礫石床對混凝土擋土墻的側(cè)壓力
在下圖中，您可以看到一個高度為3米的沙丘，由一個混凝土擋土墻固定。在這個例子中，挖掘面上施加了7千帕的頂部壓力。
為了解決這個問題，使用了兩個分析步，在第一步施加了載荷的重量和壓力，第二步是混凝土擋土墻的位移。
這是從Sample Glenn的第7章的例子3中選擇并用Abaqus軟件編寫的。在這個例子中，我們將Abacus軟件的結(jié)果與Sam Helleny的書的結(jié)果進行了比較。
在下圖中，您可以看到執(zhí)行器故障期間凹痕處塑性應(yīng)變的分布。
該任務(wù)的目的是在三種條件下獲得側(cè)壓力
第一模式：
靜止土壓力
當(dāng)擋土墻無法向后移動或朝向溝渠移動時，土壤上的壓力仍然存在。
在這個問題的模擬的第一部分中，我們使用Abaqus軟件計算下降狀態(tài)下土壤的相鄰側(cè)向載荷，并將其與Sam 與Helewani的書的答案進行比較。
在下圖中，您可以看到樣本簿中找到的停滯狀態(tài)的并排副作用圖
在下圖中，您可以看到從Abaqus軟件獲得的基臺側(cè)波動圖
第二種模式
土壓力或主動土壓力
當(dāng)擋土墻具有向后移動并遠(yuǎn)離土壤背面的能力時，就會發(fā)生土壤壓力。
在這個問題的模擬的第二部分中，我們使用Abacus軟件計算土壤抑制器的超導(dǎo)性，并將其與Sam Helleny的書的結(jié)果進行比較。
展開 
ABAQUS混凝土損傷塑性模型-C30EXCEL計算表格（含計算結(jié)果）
這是根據(jù)GB50010-2010中混凝土結(jié)構(gòu)設(shè)計規(guī)范中的混凝土本構(gòu)模型，結(jié)合文獻所述的損傷因子定義，編制的計算C30混凝土非彈性應(yīng)變和損傷因子的EXCEL表格。也是邊學(xué)變做，希望能和大家多交流。
C30砼本構(gòu)（損傷塑性模型）.rar
ABAQUS混凝土損傷塑性損傷因子計算依據(jù).rar
ABAQUS混凝土損傷塑性模型-C30EXCEL計算表格（含計算結(jié)果） ￥3.7
計算表格（如下）中標(biāo)黃部分的參數(shù)可自行設(shè)定后，EXCEL程序會自動計算“抗拉強度、非彈性應(yīng)變、受拉損傷因子”。
ABAQUS混凝土損傷塑性模型-2010混凝土結(jié)構(gòu)設(shè)計規(guī)范中C50混凝土-彈模34400Mpa-損傷因子計算及EXCEL
這是我自己計算的2010規(guī)范中ABAQUS混凝土損傷塑性模型-2010混凝土結(jié)構(gòu)設(shè)計規(guī)范中C50混凝土-彈模34400Mpa-損傷因子計算及EXCEL
首先用自己的數(shù)據(jù)計算2010規(guī)范中規(guī)定的混凝土本構(gòu)關(guān)系
然后借助文件夾中02版規(guī)范的方法，計算損傷因子。
以后還會有詳細(xì)計算方法，此數(shù)據(jù)僅供參考。
2010規(guī)范用C50混凝土損傷塑性本構(gòu)關(guān)系數(shù)據(jù)-彈模34400MPa-帶損傷因子-自己數(shù)據(jù)計算得出.rar
展開 ABAQUS中混凝土CDP模型的參數(shù)計算(下附python程序代碼) ￥19
<p>因為要仿真混凝土破壞實驗，考慮用abaqus里面的CDP模型，查閱了相關(guān)資料進行了理論總結(jié)，并根據(jù)理論編寫計算程序。</p><p>ABAQUS中CDP 模型中采用的是混凝土在單軸受力狀態(tài)下的應(yīng)力和非彈性應(yīng)變，非彈性應(yīng)變根據(jù)混凝土的單軸應(yīng)力-應(yīng)變曲線換算。
ABAQUS UEL 二次開發(fā)（Koyna混凝土壩地震響應(yīng)計算）
概述：以Koyna混凝土壩為對象進行地震響應(yīng)計算。將自編的八節(jié)點UEL和二十節(jié)點UEL應(yīng)用到計算中。分別進行了混凝土壩模態(tài)計算和地震時程計算。
其中，在模態(tài)計算中共設(shè)置四種計算工況，分別為：ABAQUS-C3D8、UEL-C3D8、ABAQUS-C3D20、UEL-C3D20。
在地震時程計算中設(shè)置兩種計算工況，分別為：ABAQUS-C3D8、UEL-C3D8。
計算結(jié)果表明，自編UEL與ABAQUS自帶單元結(jié)果一致。
（）模型信息
Koyna混凝土重力壩位于印度孟買東南200 km處，1967年12月11日，Koyna混凝土重力壩遭受里氏6.5級的地震（Koyna地震），該地震給大壩和水電站等水工建筑帶來了巨大的損壞，給下游數(shù)十萬居民的人身生命財產(chǎn)安全造成了巨大損失，該大壩地震案例事后成為諸多學(xué)者進行地震作用下壩基動力相互作用、混凝土材料的動態(tài)力學(xué)性能等領(lǐng)域研究的對象。
Koyna混凝土重力壩的橫斷面尺寸如下圖：
該壩壩高113 m，壩頂寬度14.8 m，壩底寬度70 m，正常蓄水位91.75 m，壩基尺寸分別取向上游、下游延伸兩倍壩高，向地基方向同樣延伸兩倍壩高，向橫河向延伸20 m，最終的三維幾何模型示意圖如下圖：
采用六面體單元離散，有限元計算模型如下圖：
該模型的壩體和壩基共計用17950個六面體單元，其中壩體5000個，壩基12950個，壩基的網(wǎng)格采用疏密漸進過渡的方式避免計算結(jié)果在網(wǎng)格突變處不連續(xù)。
展開 