2019年混凝土損傷塑性模型之小梁算例

1、混凝土損傷塑性模型之小梁算例i=s本帖最后由 blueshell 于 2009-3-2 10:54 編輯/i混凝土損傷塑性模型的收斂性一直是一個(gè)很頭疼的問(wèn)題，論壇中有不少高手提供了不少的思路和方法。本人也嘗試了包括ABA中提供的Riks等方法進(jìn)行了測(cè)試分析，最后采用位移加載法，成功實(shí)現(xiàn)模擬，現(xiàn)整 理供討論。od基本資料：假想一根矩形鋼筋混凝土截面梁，梁長(zhǎng)10m，截面尺寸為0.5 x 0.8m，采用5跟鋼筋，每根鋼筋面積為0.0008mA2 (此處為假想值，截面配筋率為 1%，鋼筋采用truss單元，embedded 到20 節(jié)點(diǎn)的混凝土單元中)，混凝土材料采用 Abaqus Verificat

2、ion Manual 2.2.24 Concrete damagedplasticity 例子中給出的混凝土的參數(shù)。采用兩點(diǎn)加載，加載點(diǎn)距梁端 3.5m。支座中心距梁端 0.15m 。 模型采用位移加載法，控制位移值為 0.05m，但算到一半多的時(shí)候，發(fā)現(xiàn)部分損傷值已超過(guò)0.9，就人為終止了計(jì)算。鋼筋未考慮塑性強(qiáng)化。(壓縮文件為完整inp文件)rocksoul 發(fā)表于 2009-2-19 21:34老兄，弄的不錯(cuò)哦。其實(shí)有限元軟件材料非線性計(jì)算不容易收斂真是一個(gè)令人頭疼的問(wèn)題。不收斂原因涉 及本構(gòu)、網(wǎng)格質(zhì)量、模型幾何形狀等。你這個(gè)位移加載法還不錯(cuò)。學(xué)習(xí)了。blueshell 發(fā)表于 2009-

3、2-19 21:40拜讀了混凝土損傷塑性模型應(yīng)用研究(結(jié)構(gòu)工程師2008.04 ) 一文，文中提到增強(qiáng)收斂性方法：求解混凝土問(wèn)題，一般都會(huì)碰到計(jì)算收斂困難。特別是當(dāng)混凝土開(kāi)裂后，應(yīng)變能突然釋放，使計(jì)算變得極 不穩(wěn)定。在ABAQUS中可用RIKS (弧長(zhǎng))法來(lái)處理這種不穩(wěn)定的非線性問(wèn)題，但作者通過(guò)反復(fù)試算發(fā)現(xiàn)，此計(jì)算結(jié)果并不穩(wěn)定，初始弧長(zhǎng)的設(shè)置對(duì)其計(jì)算結(jié)果影響極大，故而舍棄。為了抓住混凝土材料的主要本 質(zhì)，作者認(rèn)為在材料參數(shù)上力求準(zhǔn)確，而在計(jì)算收斂困難及計(jì)算效率很低的情況下，可采用一定的措施來(lái) 幫助收斂。歸納如下：(1)采用位移加載，初始步長(zhǎng)取0.005 ;并采用自動(dòng)時(shí)間步長(zhǎng)。(2 )通過(guò)引入粘

4、性系數(shù)對(duì)本構(gòu)方程進(jìn)行修正，考慮到過(guò)大的粘性系數(shù)會(huì)使結(jié)構(gòu)變硬”經(jīng)反復(fù)試算，粘性系數(shù)取0.0005 可以獲得較理想的計(jì)算結(jié)果。(3)打開(kāi)ABAQUS的連續(xù)性開(kāi)關(guān)，(*controls,analysis=discoutinuous)關(guān)鍵詞來(lái)幫助收斂。(4 )在混凝土有限元分析中，能夠得到問(wèn)題的解有時(shí)遠(yuǎn)比單純追求計(jì)算精度更為重要，所以必要時(shí)還可以通過(guò)調(diào)整收斂設(shè)置(*controls )關(guān)鍵詞的力和位移收斂準(zhǔn)則來(lái)使計(jì)算順利進(jìn)行?！痹诖?，想探討幾個(gè)問(wèn)題：一1. 很多人都提到采用位移加載的效果更好，那么如何確定加載的位移總量和初始步長(zhǎng)？2. 粘性系數(shù)如何選取更為合理？3. Riks method or Mo

5、dified Riks method是一種計(jì)算非線性很有效的方法，為什么在混凝土開(kāi)裂計(jì)算反而比較困難？4. 有不少人在用 UMAT子程序做損傷塑性分析，其計(jì)算效果跟ABA自帶的concrete damageplasticity 相比怎樣？其收斂性又是如何處理的？rocksoul 發(fā)表于 2009-2-19 21:58i=s 本帖最后由 敦誠(chéng) 于2009-3-2 12:16 編輯/i粘性系數(shù)我以前好像算過(guò)，沒(méi)有記錯(cuò)好像千分之一或二也會(huì)獲得不錯(cuò)效果的。比較關(guān)系UAMT損傷本構(gòu)和aba自帶損傷本構(gòu)收斂性這個(gè)問(wèn)題，有做過(guò)不妨透漏下，郭版主可以談?wù)?color=Red 是 dun敦，不是 guo郭敦誠(chéng)/

6、color敦誠(chéng) 發(fā)表于2009-3-2 11:34談?wù)剛€(gè)人對(duì)于abaqus中混凝土材料的理解吧！其實(shí)準(zhǔn)靜態(tài)時(shí)，混凝土在下降段的的網(wǎng)格依賴(lài)效應(yīng)是眾所 周之的產(chǎn)生的原因是運(yùn)動(dòng)方程的退化?，F(xiàn)階段兩種方式可以較好解決的這一問(wèn)題，一種就是樓主所提供的加入粘性參數(shù)為了減緩在峰值以后的迅 速跌落，混凝土材料在拉伸的過(guò)程中岀現(xiàn)準(zhǔn)脆性斷裂，所以這種方式尤其是在研究混凝土受拉破壞的時(shí)候 長(zhǎng)長(zhǎng)采用。另外一種就是比較熱門(mén)的梯度依賴(lài)?yán)碚摚∵@個(gè)我也就不多說(shuō)了！至于粘性系數(shù)選取，我比較贊同blueshell兄和rocksoul兄的意見(jiàn)，千分之幾得到的效果還是比較合理的。 對(duì)于blueshell兄的幾個(gè)問(wèn)題，我是這樣理解的。

7、位移總量是根據(jù)一些試驗(yàn)來(lái)判斷的，一般是先找做過(guò)的試驗(yàn)，然后確定模擬的位移加載量；初始增量一般都要設(shè)的很小1e10，當(dāng)然如果選擇restart就沒(méi)有這個(gè)問(wèn)題了，你可以首先預(yù)測(cè)一下塑性開(kāi)始的位置，然后restat的時(shí)候?qū)⒃隽吭O(shè)的非常小就行了！比較贊同rocksoul兄的觀點(diǎn)，千分之幾比較合理。 .-(3. Riks method or Modified Riks method不能直接就是用，因?yàn)樵谑褂没￠L(zhǎng)法的時(shí)候，會(huì)出現(xiàn)負(fù)加載的情況，而混凝土材料的拉壓區(qū)別十分明顯，在計(jì)算拉伸的時(shí)候，riks發(fā)現(xiàn)施加拉伸載荷過(guò)不去，而壓縮能過(guò)得去，就認(rèn)為那里是分叉點(diǎn)，然后開(kāi)始施加壓縮載荷，這樣就發(fā)生了ls所說(shuō)的情況

8、了，一般riks是加在一個(gè)standard 后面來(lái)計(jì)算的！我做過(guò)的umat與damage plastic 在拉伸載荷與壓縮的時(shí)候比較像，但是在 shear和rotation 的時(shí)候 差別還是很明顯的。umat收斂性就是注意應(yīng)力積分方式，這個(gè)說(shuō)起來(lái)就復(fù)雜了，我有半年的時(shí)間是在對(duì) 于應(yīng)力積分進(jìn)行試驗(yàn)！呵呵個(gè)人觀點(diǎn)，歡迎大家提岀意見(jiàn)和建議！xiaoqg2003 發(fā)表于 2009-3-2 13:24文件好大，我的機(jī)器要算一天一夜才能算完沒(méi)能用對(duì)稱(chēng)性縮減一下規(guī)模就好了blueshell 發(fā)表于 2009-3-2 14:24敦版主可是這一塊兒的高人了。拜讀過(guò)你不少的帖子，受益頗多。在此，也再提一下，后來(lái)才

9、發(fā)現(xiàn)首帖中 的inp文件中粘性參數(shù)取為 0.0005，好像是比較大了，估計(jì)跟計(jì)算能順利收斂也有關(guān)系。不過(guò)，現(xiàn)在沒(méi) 有時(shí)間去做對(duì)比分析了。你提到的在收斂性問(wèn)題嚴(yán)重的步之前考慮采用restart的方法很好，應(yīng)能節(jié)省大量分析時(shí)間。關(guān)于弧長(zhǎng)法，我用上邊的模型也算過(guò)，確實(shí)也獲得了負(fù)加載的情況。Umat就更復(fù)雜了，還得花時(shí)間去學(xué)習(xí)學(xué)習(xí)！回應(yīng)xiaoqg2003 ，這個(gè)模型雖然是根小梁，但涉及開(kāi)裂，網(wǎng)格劃分較密，且增量步一般也比較小，所以, 計(jì)算起來(lái)是很費(fèi)時(shí)的，我的那個(gè)計(jì)算結(jié)果文件都超過(guò)2G !得配個(gè)powerful的計(jì)算機(jī)！呵呵敦誠(chéng) 發(fā)表于 2009-3-3 18:15 re : 10 -blueshel

10、l兄過(guò)獎(jiǎng)了，我那是什么高手啊，就是本科一直到碩士做混凝土方面較多！其實(shí)根據(jù)個(gè)人經(jīng)驗(yàn)來(lái)看，一般選擇standard來(lái)做混凝土損傷的話加了rebar或者鋼筋，計(jì)算會(huì)順利一些，收斂性能也好了很多！但是要是單純做素混凝土的損傷還是應(yīng)該選擇explicit ，但是在加載的時(shí)候速度要足夠??！確保整個(gè)過(guò)程是準(zhǔn)靜態(tài)的！bb20040128 發(fā)表于 2009-3-9 20:47quote精品啊，學(xué)習(xí)了！談?wù)剛€(gè)人對(duì)于abaqus中混凝土材料的理解吧！其實(shí)準(zhǔn)靜態(tài)時(shí)，混凝土在下降段的的網(wǎng)格依賴(lài)效應(yīng)是眾所周之的產(chǎn)生的原因是運(yùn)動(dòng)方程的退化。一現(xiàn)階段兩種方式可以較好解決的這一問(wèn)題，一種就是樓主所。_WJysize=2col

11、or=#999999敦誠(chéng) 發(fā)表于 2009-3-2 11:34/colorurl= g-.color=blacksize=2 敦誠(chéng)版主/size -size=2個(gè)人認(rèn)為版主所說(shuō)第三點(diǎn)的分叉點(diǎn)不應(yīng)該是指計(jì)算這種情況的時(shí)候的分叉點(diǎn)導(dǎo)致弧長(zhǎng)法計(jì)算會(huì)岀現(xiàn)的問(wèn)題的分叉點(diǎn)我以我個(gè)人的理解舉個(gè)例子說(shuō)下軸心受壓的理想桿件承受一定軸力后會(huì)岀現(xiàn)分 叉點(diǎn)穩(wěn)定一條路徑強(qiáng)度一條路徑具有不連續(xù)響應(yīng)趨勢(shì) 這個(gè)是弧長(zhǎng)法計(jì)算無(wú)法解決的對(duì)于混凝土的計(jì)算我不太懂但是我認(rèn)為結(jié)構(gòu)響應(yīng)應(yīng)該是連續(xù)的應(yīng)該可以用弧長(zhǎng)法直接計(jì)算 這是我個(gè)人看法版主長(zhǎng)期做混凝土的 請(qǐng)指教/size/color敦誠(chéng) 發(fā)表于 2009-3-14 15:48quote

12、敦誠(chéng)版主個(gè)人認(rèn)為版主所說(shuō)第三點(diǎn)的分叉點(diǎn)不應(yīng)該是指計(jì)算這種情況的時(shí)候的分叉點(diǎn)導(dǎo)致弧長(zhǎng)法計(jì)算會(huì)岀現(xiàn)的問(wèn)題的分叉點(diǎn)我以我個(gè)人的理解舉個(gè)例子說(shuō)下軸心受壓的理想桿件承受一定軸力后會(huì)岀現(xiàn)分叉點(diǎn)穩(wěn)定 一條.發(fā)表于 2009-3-9 20:47/colorsize=2color=#999999bb20040128 url=gre : bb20040128 兄其實(shí)我所指的分叉并不是分叉理論時(shí)的分叉，我也就是借用了一下這個(gè)術(shù)語(yǔ)，不好意識(shí)??！_其實(shí)riks對(duì)于結(jié)構(gòu)失穩(wěn)的分叉求解是可以的，但是對(duì)于材料失穩(wěn)，尤其是混凝土材料時(shí)候就不好用，其實(shí) 簡(jiǎn)單的材料失穩(wěn)談不上什么分叉理論，但是在abaqus計(jì)算中這個(gè)比較簡(jiǎn)單的下降

13、段就會(huì)遇到麻煩，根本原因不在于求解的技巧，而是在于計(jì)算理論本身就存在問(wèn)題，為了彌補(bǔ)這一問(wèn)題，一下方法被采用了，首 先riks對(duì)于這樣的問(wèn)題本身并不是專(zhuān)屬的，但是它的無(wú)量綱形式以及平衡點(diǎn)的求解方式使得它對(duì)于失穩(wěn)結(jié)構(gòu)方面能夠計(jì)算出下降段，而混凝土結(jié)構(gòu)也會(huì)有下降段，所以在此相當(dāng)于做了個(gè)借用！但是用于這種借用本來(lái)就沒(méi)有什么特別的針對(duì)性（即針對(duì)材料失穩(wěn)實(shí)質(zhì)是運(yùn)動(dòng)，平衡方程的退化），所以 岀了毛病是很正常的，；.-當(dāng)時(shí)的第三點(diǎn)也是我在解決這些問(wèn)題時(shí)采用的一種方法，實(shí)話實(shí)說(shuō)的話效果還是一般！請(qǐng)bb20040128 兄多多指導(dǎo)啊！bb20040128 發(fā)表于 2009-3-16 13:46quotere :

14、bb20040128 兄其實(shí)我所指的分叉并不是分叉理論時(shí)的分叉，我也就是借用了一下這個(gè)術(shù)語(yǔ)，不好意識(shí)啊！其實(shí)riks對(duì)于結(jié)構(gòu)失穩(wěn)的分叉求解是可以的，但是對(duì)于材料失穩(wěn)，尤其是混凝土材料時(shí)候就不好用，其 實(shí)簡(jiǎn)size=2color=#999999敦誠(chéng)發(fā)表于 2009-3-14 15:48/colorurl= g呵呵是俺學(xué)習(xí)啦謝謝dlx850113 發(fā)表于 2009-4-1 15:2315樓主你的文章已經(jīng)拜讀了非常不錯(cuò)但是參數(shù)設(shè)置方面我還是有點(diǎn)困惑比如剪脹角ABA給岀的是還有壓縮時(shí)貌似給出了recovery factor，其值為1.0再提混凝土彈塑性損傷模型的參數(shù)定義在網(wǎng)上查了很多資料，幾乎沒(méi)有一個(gè)

15、對(duì)這里面的參數(shù)進(jìn)行過(guò)解釋或者對(duì)其有根據(jù)的來(lái)源進(jìn)行闡述，所以這 里再討論一下。一我對(duì)混凝土 C50，根據(jù)混凝土規(guī)范附錄 C中的單軸受壓受拉應(yīng)力應(yīng)變關(guān)系，進(jìn)行了求解，詳細(xì)見(jiàn)附件問(wèn)題：1. 受壓真實(shí)應(yīng)力一塑性應(yīng)變起點(diǎn)的確定從附件中可知塑性應(yīng)變由正的轉(zhuǎn)變?yōu)樨?fù)值（因?yàn)槭軌海敲雌瘘c(diǎn) 是不是就是在0的地方？2. 受拉真實(shí)應(yīng)力-塑性應(yīng)變起點(diǎn)的確定從附件中可知塑性應(yīng)變均由正的，起點(diǎn)是不該從那一點(diǎn)開(kāi)始呢？3. 彈性模量的確定3.1從附件受壓曲線應(yīng)力一應(yīng)變求岀的彈性模量3.2從附件受拉曲線應(yīng)力一應(yīng)變求岀的彈性模量3.3由規(guī)范查得混凝土受壓受拉彈性模量-p其中3.1值和3.3值接近，但是3.2得值就差得比較多，此

16、時(shí)該如何取?4. 如何由塑性應(yīng)變求受拉受壓損傷參數(shù)？5. 如何確定在循環(huán)荷載下得Wt和Wc值？ -wei1012 發(fā)表于 2007-6-15 20:31對(duì)于上面我發(fā)的幾個(gè)問(wèn)題我現(xiàn)在就自己回答一下：1. 起點(diǎn)就是混凝土塑性應(yīng)變?yōu)?的點(diǎn)，這個(gè)是和abaqus中的定義要求統(tǒng)一的；2. 根據(jù)abaqus中混凝土的受拉曲線來(lái)看，應(yīng)該就是從受拉開(kāi)裂應(yīng)變那點(diǎn)開(kāi)始的；3. 彈性模量就取受壓的彈性模量。如果是對(duì)結(jié)構(gòu)而言，這個(gè)影響不會(huì)太大；如果是對(duì)試驗(yàn)試件，那么應(yīng)該 用實(shí)驗(yàn)得到這個(gè)值；-m4. 受拉受壓損傷參數(shù)根據(jù)參考資料取獲得的。最好的損傷曲線應(yīng)該是abaqus中幫助文檔里面的verification2.2.2

17、3 了。5. 在循環(huán)荷載下得 Wt和Wc :一般為了簡(jiǎn)化計(jì)算，就直接取Wt = 0，Wc = 1。wsg 發(fā)表于 2007-9-5 10:12請(qǐng)問(wèn)wei1012，你給出的附件里面確定是C50的數(shù)據(jù)嗎？你說(shuō)采用的是混凝土規(guī)范附錄C中的單軸受壓受拉應(yīng)力應(yīng)變關(guān)系，可是我對(duì)照了一下，你的數(shù)據(jù)好像不是 C50的。難道采用的是fcm對(duì)應(yīng)的參數(shù)值？ 另外有問(wèn)題想問(wèn)：采用該應(yīng)力應(yīng)變曲線需要對(duì)著表C.1.2來(lái)選取參數(shù)，公式中的混凝土單軸抗壓強(qiáng)度f(wàn)c*可以是fck、fc或fcm。C20的fck只有13.4，而表C.1.2只給出了 fc* = 15數(shù)值以后的參數(shù)值，這樣 想通過(guò)插值獲得參數(shù)都不能，請(qǐng)問(wèn)該如何解決這個(gè)

18、問(wèn)題？2.彈性模量可否去規(guī)范給岀的受壓受拉彈性模量？shashami 發(fā)表于 2007-10-8 23:23第一個(gè)和第二個(gè)問(wèn)題 wei1012說(shuō)的很對(duì)后面幾個(gè)參數(shù)不是可以根據(jù)規(guī)范和混凝土強(qiáng)度計(jì)算得出的嗎？如果有實(shí)驗(yàn)當(dāng)然更好.淵源 1983 發(fā)表于 2007-10-23 20:56這么好的帖子不能沉哦繼續(xù)頂起來(lái)。另外在提一個(gè)問(wèn)題：A 一我用規(guī)范提供的數(shù)據(jù)算 C30受壓塑性應(yīng)變?nèi)秦?fù)值我的塑性應(yīng)變值取不到 0值那么按照LZ的的觀點(diǎn)這個(gè)起點(diǎn)我取什么呢？一fc取的是標(biāo)準(zhǔn)值 E取的3E4南柯一夢(mèng)發(fā)表于2007-10-31 22:35本人也曾對(duì)此問(wèn)題研究過(guò)好長(zhǎng)時(shí)間，據(jù)文獻(xiàn)知受拉受壓損傷參數(shù)應(yīng)該是根據(jù)混凝土

19、的滯洄曲線得岀的，估計(jì)沒(méi)有一個(gè)統(tǒng)一的公式計(jì)算。xj983 發(fā)表于 2008-1-20 20:34幫定一下，盼望牛人將混凝土規(guī)范，C30-C60的都計(jì)算一遍就好了斑竹，用于計(jì)算收斂性如何呢kawakuki 發(fā)表于 2008-1-28 09:36樓主采用的參數(shù)如何定出？C50的sc應(yīng)該取值為1920X 10-6，為何取1781X 10-6 ? aa應(yīng)該取1.78為何取1.906 ?有什么說(shuō)法么？向 ABA轉(zhuǎn)換的公式出處在哪里？請(qǐng)指教，謝謝。/size/fontalbertyu 發(fā)表于 2008-5-28 21:45目前對(duì),損傷參數(shù)的確定還是一個(gè)難題，大多數(shù)只能通過(guò)實(shí)驗(yàn)獲取.然而在我過(guò)去的研究中發(fā)現(xiàn)

20、，損傷參數(shù)的 引入對(duì)單向加載的模擬影響并不大.caliber 發(fā)表于 2008-11-23 10:58其實(shí)關(guān)于混凝土損傷參數(shù)的確定是難以有統(tǒng)一的結(jié)果，最合適的方法是試驗(yàn)數(shù)據(jù)。但這個(gè)數(shù)據(jù)僅對(duì)非常接 近試驗(yàn)的各種狀態(tài)下才是有效而真實(shí)的。在工程實(shí)踐計(jì)算中你可以發(fā)現(xiàn)，其實(shí)其影響是非常有限的。通過(guò) 取不同的損傷參數(shù)對(duì)宏觀的結(jié)果的影響基本可以忽略，當(dāng)然做試驗(yàn)對(duì)照是另當(dāng)別論了。彈性模量可以取進(jìn) 入塑性臨界狀態(tài)的割線模量。ABAQUS Concrete damaged plasticity中混凝土的損傷變量 D的計(jì)算方法大家都知道ABAQUS中認(rèn)為損傷變量是塑性應(yīng)變、溫度、場(chǎng)變量的函數(shù)即血=如護(hù)4齊）；0&l

21、t; <1.dc =血（聖血九）;0 < rfc <打損傷變量取值可從 o （代表無(wú)損傷）到1 （完全失去強(qiáng)度）但是，ABAQUS幫助文檔中并未給出損傷變量 Dt,Dc的相應(yīng)計(jì)算公式，如果不輸入損傷變 量，那就不是塑性損傷模型了，而稱(chēng)之為塑性模型。那損傷變量如何計(jì)算呢？首先我們了解下?lián)p傷變量，損傷是指在單調(diào)加載或重復(fù)加載下，材料性質(zhì)所產(chǎn)生的一種劣化現(xiàn)象，材料的損傷可以用損傷變量來(lái)描述，損傷變量是力學(xué)中描述材料內(nèi)部性質(zhì)劣化的宏觀變量，它的定義有很多種，有各向同性彈性損傷模型，各向異性彈性損傷模型，彈塑性損傷 模型和隨即損傷模型等等，因此損傷變量的定義也有很多種。這里個(gè)人推薦兩種

22、比較經(jīng)典的損傷模型的損傷變量計(jì)算方法，可近似用與ABAQUS計(jì)算。1:經(jīng)典損傷理論的混凝土損傷變量計(jì)算方法：D=1-Es/E0其中Es為應(yīng)力應(yīng)變曲線上任一點(diǎn)的割線模量。（經(jīng)過(guò)分析，該公式為彈性損傷模型， 計(jì)算的損傷變量偏大， 不適合ABAQU塑性損傷計(jì)算,輸入后會(huì)報(bào)錯(cuò)）2： Mazars模型混凝土損傷變量計(jì)算方法：Q I f（l A）cl 7<cxp（ - 一 £f）|式中:A,B為材料常數(shù)，由實(shí)驗(yàn)確定，對(duì)一般混凝土材料0.7<A<1, 10 4<B<105, £ f為對(duì)應(yīng)于（T f的應(yīng)變。（計(jì)算公式較復(fù)雜，但有實(shí)驗(yàn)數(shù)據(jù)說(shuō)話，可信度高，可以參考

23、）3:適合鋼纖維混凝土損傷變量計(jì)算方法：根據(jù)材料的全應(yīng)力-應(yīng)變曲線特征，選擇 Weibull分布的密度函數(shù)模擬全應(yīng)力應(yīng)變曲線，又因?yàn)椴牧系膹?qiáng)度服從Weibull統(tǒng)計(jì)分布，可認(rèn)為材料的損傷變量D也服從該分布：式中E pk, S pk為峰荷應(yīng)變值、峰值應(yīng)力。 還有其他的一些模型以后再補(bǔ)充 以上為個(gè)人愚見(jiàn)，若有不妥，歡迎拍磚!1經(jīng)典損傷理論的混凝土損傷變量計(jì)算方法：D=1-Es/E0其中Es為應(yīng)力應(yīng)變曲線上任一點(diǎn)的割線模量。（直觀并且計(jì)算起來(lái)非常簡(jiǎn)單，個(gè)人強(qiáng)烈推薦）請(qǐng)問(wèn)，這里計(jì)算 Es時(shí)，應(yīng)變不是INELASTIC PLASTIC 吧？為何我計(jì)算后帶入abaqus中時(shí)出現(xiàn)錯(cuò)誤呢？” During t

24、he conversion from crushing to plastic strain abaqus found negative and/or decreasing values of plastic strain. Verify that the degradati on data un der *con crete compressi on da mage is correct1的公式是有點(diǎn)問(wèn)題這個(gè)公式采用的Es斜率起點(diǎn)在零點(diǎn)處而損傷公式里面采用的 E起點(diǎn)實(shí)在塑性應(yīng)變處但是ABA定義塑性應(yīng)變又是基于D的感覺(jué)ABA對(duì)于塑性損傷模型定義的一套公式是一種循環(huán)公式 就是這個(gè)變量基于那個(gè)變量那

25、個(gè)變量又基于這個(gè)變量象這里的塑性應(yīng)變和損傷變量就是互相定義計(jì)算的時(shí)候不知道怎么輸入2007- 10-31 13:14 一 一 一上-="zyy1983 wrote:1的公式是有點(diǎn)問(wèn)題這個(gè)公式采用的Es斜率起點(diǎn)在零點(diǎn)處而損傷公式里面采用的E起點(diǎn)實(shí)在塑性應(yīng)變處但是ABA定義塑性應(yīng)變又是基于D的感覺(jué)ABA對(duì)于塑性損傷模型定義的一套公式是一種循環(huán)公式 就是這個(gè)變量基于那個(gè)變量對(duì)，經(jīng)過(guò)分析可知：經(jīng)典損傷理論的混凝土損傷變量計(jì)算方法并不適用與ABAQUS計(jì)算按該公式輸入后確實(shí)會(huì)出現(xiàn)如下錯(cuò)誤：” During the conversion from crushing to plastic stra

26、in abaqus found negative and/or decreasing values of plastic strain. Verify that the degradati on data un der *con crete compressi on damage is correct也就是說(shuō)按（1）式計(jì)算結(jié)果損傷變量偏大，導(dǎo)致計(jì)算得到的塑性應(yīng)變值 隨著非彈性應(yīng)變的增加變?yōu)樨?fù)值或減小，ABAQUS認(rèn)為損傷變量定義有錯(cuò)誤。感覺(jué)ABAQUS理論手冊(cè)不太全面了。herojoe wrote:那怎么辦呀？要是用 2 公式的話，沒(méi)有試驗(yàn)呀難道用嘗試法，呵呵不要著急，我在測(cè)試當(dāng)中，在近期給與

27、調(diào)試結(jié)果。論壇以前都沒(méi)有專(zhuān)門(mén)討 論損傷變量的計(jì)算方法。并且ABAQUS理論手冊(cè)中說(shuō)的也很含糊，要怪就要ABAQUS公司不給計(jì)算公式了。但是ABAQUS中混凝土塑性損傷模型的損傷思路確實(shí)很明確： 即描述損傷問(wèn)題時(shí)用損傷因子D=1- 3 =(AA')/A其中A為無(wú)損傷狀態(tài)的橫截面面積，A'為損傷后的有效承載截面積。然而在受損材料中，從細(xì)觀上對(duì)缺陷形式和損傷確定有效承載面積A'是很困難的，因此為了間接測(cè)定損傷，通過(guò)應(yīng)變等價(jià)原理認(rèn)為D=1-E'/E其中E'為受損材料彈性模量，也稱(chēng)為有效彈性模量。E的計(jì)算原理為卸載線的斜率。但是關(guān)鍵問(wèn)題是：ABAQUS并未給出卸載

28、線的斜率怎么計(jì)算，近似用割線模量Es代替的話計(jì)算的損傷變量又偏大，最后導(dǎo)致submit job后出現(xiàn)報(bào)警，只能通過(guò)其他損傷變量計(jì)算方法來(lái)計(jì)算D 了。這個(gè)問(wèn)題還有待進(jìn)一步商榷。2007-11-1 10:13 一 一上-=''''cyasj_whu 老大，有辦怙法了嗎？還有，就是想請(qǐng)教。我在做構(gòu)件的骨架曲線時(shí)也要考慮損傷嗎？用塑性模型可以嗎？我之前用了塑性模型來(lái)模擬過(guò)一個(gè)試驗(yàn)結(jié)果，在軸壓比較小的時(shí)候吻合的還是蠻好的，可軸壓比較大時(shí)就有差距了，我認(rèn)為是混凝土本構(gòu)有問(wèn)題，因?yàn)槲也捎玫膯屋S應(yīng)力應(yīng)變，可在鋼管混凝 土中，混凝土應(yīng)該是三向受力狀態(tài)，所以在軸壓比較大時(shí)就出現(xiàn)了偏

29、差。不知道我的分析對(duì)不對(duì)，請(qǐng)大俠指教。謝謝！2007-11-1 15:29 一 一上-=''''根據(jù)測(cè)試I養(yǎng)第一種方法為彈性損傷 模型，計(jì)算出的損傷變 量必會(huì)出現(xiàn)如下錯(cuò)誤"Duri ng the con vers ion from crush ing to plastic stra in abaqus found negative and/or decreasing values of plastic strain. Verify that the degradati on data un der *con crete compressi on dam

30、age is correct 不適合ABAQUS計(jì)算按照第二種方法，根據(jù)著名Mazars損傷演化得的損傷變量通過(guò)計(jì)算，適合 ABAQUS對(duì)混凝土計(jì)算，第三種方法適合對(duì)鋼纖維混凝土計(jì)算。問(wèn)題解決。cyasj_whu 修改于 2007-11-01 15:32herojoe wrote:我在做構(gòu)件的骨架曲線時(shí)也要考慮損傷嗎？用塑性模型可以嗎？我之前用了塑性模型來(lái)模擬過(guò)一個(gè)試驗(yàn)結(jié)果，在軸壓比較小的時(shí)候吻合的還是蠻好的，可軸壓比較大時(shí)就有差距了，我認(rèn)為是混凝土本構(gòu)有問(wèn)題，因?yàn)槲也捎玫膯屋S應(yīng)力應(yīng)變，可在鋼管混凝 土中，混凝土應(yīng)該是三向受力狀態(tài)，所以在軸壓比較大時(shí)就出現(xiàn)了偏差。 不知道我的分析對(duì)不對(duì)1：不輸

31、入損傷變量，就是塑性模型了，原理基于彈塑性力學(xué)，當(dāng)然可以 用。2:輸入的是單軸應(yīng)力-應(yīng)變關(guān)系沒(méi)錯(cuò)，約束混凝土是三向受力，軟件根據(jù) 你輸入的單軸應(yīng)力-應(yīng)變關(guān)系運(yùn)用一定的關(guān)系（比如，破壞準(zhǔn)則、非線性指 標(biāo)、等效應(yīng)力-應(yīng)變等）將單軸泛化到多軸。其實(shí)有些東西還是需要力學(xué)專(zhuān)業(yè)的專(zhuān)家們來(lái)解釋?zhuān)瑢?duì)于非力學(xué)專(zhuān)業(yè)，比如說(shuō)我，除了學(xué)習(xí)三大力學(xué)外，光個(gè)損傷力學(xué)就把我整損傷了。更何況還有其他的比如塑性力學(xué)、固體力學(xué)、斷裂力學(xué)、連續(xù)介質(zhì)力學(xué)等等等等。2007-11-1 16:36 一 一上-=''''謝謝，cyasj_whu !挺 幽默的哈！1、那我就用第2種算法先嘗試一下，那這樣的話，

32、我的A.B值就要自己嘗試了。能否給個(gè)推薦值呢2、 你說(shuō)根據(jù)一定關(guān)系（比如，破壞準(zhǔn)則，非線性指標(biāo)、等效應(yīng)力-應(yīng)變等） 將單軸泛化到多軸， 那我能不能在 ABAQUS中修改這些準(zhǔn)則呢，因?yàn)槲矣X(jué) 得ABAQUS在將單軸轉(zhuǎn)化為三軸的時(shí)候，混凝土的峰值提高的不多。呵呵，不好意思，老麻煩你我用方法2做了損傷變 量，結(jié)果還是顯示During the con vers ionfrom crush ing to plastic stra in abaqus found n egative an d/or decreas ing values of plastic stra in. Verify that the

33、 degradati on data un der *con crete compressi on damage is correct“錯(cuò)誤，怎么辦呢？謝謝2007-11-27 16:52 一 一 一上- _ :怎么會(huì)呢？滋我經(jīng)過(guò)測(cè)試了的，用該公式不會(huì)出現(xiàn)那種錯(cuò)誤，而且有些人用過(guò)了沒(méi)出錯(cuò)，咋到你那就有問(wèn)題了呢你要注意A,B的取值。我這列的是拉損傷，你查詢(xún)下參考文獻(xiàn)，拉損傷與 壓損傷的A,B取值范圍不一樣。cyasj_whu 修改于 2007-11-27 17:102007-11-28 11:33 一 一 一上-_ ="我受拉的時(shí)候A取得 0.8 , B取得50000.受壓的時(shí)候A=1

34、.3,B=1500 ，不知道有沒(méi)有問(wèn)題，而且我算受壓的時(shí)候到了最后De值超過(guò)1 了。到底是哪出錯(cuò)了呢？這個(gè)公式是不是適用于所有的本構(gòu)關(guān)系？推薦一篇論文，基于能量損失的混凝土損傷模型，余志武#182008- 1-15 19:39 一 一 一上-幾個(gè)月前參考樓上說(shuō)的那篇文獻(xiàn)計(jì)算過(guò)損傷變 量，效果還行，就是很 難收斂，我也不知道到底要怎么調(diào)。樓上的你如果要設(shè)定受拉損傷的時(shí)候 把塑性應(yīng)變稍微調(diào)大一些再 試試（對(duì)比單軸受拉曲線）2008-1-17 10:26 一 一上-='''' chenjua n62 :你說(shuō)效 駅果還行是什么意思，收 斂嗎？不收斂你怎么看 出效果還行的

35、？2008-1-17 21:42 一 一 一上-我說(shuō)的是計(jì)算的損傷變量，提交任務(wù)也沒(méi)有報(bào) 錯(cuò)，只是算不過(guò)去，用vertificatio n2.2.23 里的那段本構(gòu)也算不過(guò)去。提交任務(wù)得時(shí)候都沒(méi)有報(bào)錯(cuò)已經(jīng)不同意了，我算得提交任務(wù)得時(shí)候就 不行了。為什么我用上面得方法算得就出錯(cuò)那#232008-1-18 16:58 一 一 一上-如果損傷太大，提交文件得時(shí)候就直接出錯(cuò) 了，那么怎樣調(diào)節(jié)一下 它才能不出錯(cuò)那？2008-3-4 11:44.L亠-t自己調(diào)整損傷量值 D 廠廠可以通過(guò)文件的提交可是在答辯時(shí)應(yīng)該有依據(jù)拜托樓上的兄弟們Mazars模型混凝土損傷變量計(jì)算方法的相關(guān)文獻(xiàn)和書(shū)籍介紹一下謝了還有一

36、個(gè)問(wèn)題就是鋼板的材料為什么你們都沒(méi)有考慮其循環(huán)荷載下的 強(qiáng)化準(zhǔn)則呢我定義的是*harde nin g=COMBINED,DATA TYPE=HALF CYCLE4.20E+08,04.70E+08,0.02365.00E+08,0.26.00E+08,0.5*Cyclic Harde ning3.5e+07, 0.014615.4e+07, 0.07037e+07, 0.131497.8e+07, 0.189932008-3-5 15:53.L亠- t我計(jì)算的時(shí)候鋼管用的線性隨動(dòng)強(qiáng)化模型，因 為非線性的不會(huì)定義，請(qǐng)教liqian25420怎么定義阿，把你的定義的鋼材弄進(jìn)去，提交任務(wù)時(shí)提 示出錯(cuò)

37、。我用線性隨動(dòng)強(qiáng)化算得骨架曲線數(shù)值比試驗(yàn)值低，不知道用非線 性隨動(dòng)強(qiáng)化怎樣。在ANSYS上能不能模 擬損傷，如何在圖上顯示出來(lái)？高手請(qǐng)指點(diǎn)一下，謝謝!誰(shuí)有Mazars模型混凝0分土損傷變量計(jì)算方法的相關(guān)文獻(xiàn)和書(shū)籍介紹一下謝了#282009- 3-18 10:12 一 - 一上-Iily8232 wrote:0 分我用方法2做了損傷變量，結(jié)果還是顯示 During the conversion from crushing to plastic strain abaqus found neg ative and/or decreasing values of plastic strain. Ver

38、ify that the degr adation data under *concret .我計(jì)算時(shí)也出現(xiàn)了此現(xiàn)象，肯定是損傷因子計(jì)算有問(wèn)題。不知道怎么算 ?。磕奈桓呤种更c(diǎn)下啊2007-5-25 18:22 一 一 一上-以前沒(méi)用過(guò)弧長(zhǎng)法，現(xiàn)在需用它做一個(gè)靜力彈塑性分析，在一結(jié)構(gòu)中施加水平力（是通過(guò)重力形式施加的，只是重力加速度方向改為水平方向），觀察力與位移的關(guān)系， 但是出現(xiàn)一個(gè)問(wèn)題就是當(dāng)計(jì)算到一定程度以后，增量步開(kāi)始出現(xiàn)負(fù)值，也許正常計(jì)算也是這樣，可以解釋為超過(guò)結(jié)構(gòu)的極限承載力，出現(xiàn)了下降段。但是提取結(jié)果一看卻發(fā)現(xiàn)后期頂層位移開(kāi)始逐漸減小。不知是何原因？還有一個(gè)問(wèn)題，我試圖采用位移方式施

39、加，但是算到一定程度，開(kāi)始退 出，錯(cuò)誤如下：A zero displacement solution was found in the first iteration of a riks step. The an alysis cannot proceed. The most likely cause is that all displaceme nt degrees-of-freedom are con stra ined to be zero.請(qǐng)問(wèn)是材料極限強(qiáng)度達(dá)到了極限強(qiáng)度，但是還沒(méi)有產(chǎn)生大位移即極限應(yīng)變?cè)O(shè)置太小的原因造成的嗎？ 請(qǐng)高手指點(diǎn)！謝謝不論是位移加載還是FORCE加載，都應(yīng)該是

40、加在邊界條件上，因此無(wú)論是什么模型和什么 iteration方法，位移加載時(shí)，位移是持續(xù)穩(wěn)定地增長(zhǎng)的.FORCE加載也是這 樣,只是當(dāng)結(jié)構(gòu)達(dá)到極限荷載時(shí)，結(jié)構(gòu)將不能繼續(xù)承載，程序?qū)?huì)提示不能 收斂.理論上講，應(yīng)該是這樣在靜力荷載作用，結(jié)構(gòu)位移應(yīng)該是有增無(wú)減的。但是我在用riks法進(jìn)行靜 力分析試圖得到其結(jié)構(gòu)從到屈曲后的性能。但是缺算的結(jié)果確實(shí)是在位移增大到一定程度后，出現(xiàn)了減小。這是和我以前做過(guò)靜力所不同的，也是 以前所沒(méi)見(jiàn)過(guò)的問(wèn)題。不過(guò)以前再結(jié)構(gòu)靜力分析時(shí)沒(méi)有用過(guò)riks方法。因?yàn)閞iks在做彈塑性靜力分析更容易收斂，并且可以做出力-位移曲線。但是出現(xiàn)了上述問(wèn)題。如果兄臺(tái)有時(shí)間倒是可以做做試

41、試。我上述題目所顯示的錯(cuò)誤信息有意味什么呢。懇求答復(fù)！謝謝！另外我還要補(bǔ)充一句，當(dāng)結(jié)構(gòu)超過(guò)了其極限承載力，用弧長(zhǎng)法一般依然可以收斂的，只是荷載會(huì)調(diào)整變小而已。所以，事實(shí)上弧長(zhǎng)法是可以做出結(jié)構(gòu)的下降段的。有時(shí)間 咱們可以共同研究一下弧長(zhǎng)法的用處及使用方法。還望不吝賜教。2007-5-26 22:05 一 一 一上-我在做考慮結(jié)構(gòu)的壓柑 失穩(wěn)影響的時(shí)候， 也出現(xiàn)過(guò)類(lèi)似情況； 當(dāng)改變一些 脆弱”桿件的本構(gòu)時(shí)，會(huì)呈現(xiàn)不同的XY值曲線，當(dāng)本構(gòu)較弱時(shí)，有一些情況是結(jié)構(gòu)的力，到了最大值后下降，然后按彈性 的E為斜率，位移不斷減小，這就是LZ說(shuō)的情況；事實(shí)是計(jì)算的問(wèn)題，當(dāng)所加荷載結(jié)構(gòu)不能承受時(shí)， 有可能會(huì)出現(xiàn)

42、兩種情況，一，結(jié)構(gòu)被壓屈，承載力下降，但位移仍然增大；二，結(jié)構(gòu)被壓屈，出現(xiàn)塑性的情況下卸載”即位移減小，位移-荷載曲線很類(lèi)似滯回曲線的形狀；這是我做算例所發(fā)現(xiàn)的情況，由于我對(duì)具體的算法及有限元理論了解有限， 希望能有人解答。我想請(qǐng)教的是怎么樣才可以時(shí)結(jié)構(gòu)采用力加載的方式，能夠出現(xiàn)樓上所說(shuō)結(jié)構(gòu)進(jìn)入塑性，而位移繼續(xù)增大，同時(shí)可以出現(xiàn)卸載的下降段呢。我在稍有變動(dòng)的結(jié)構(gòu)模型中計(jì) 算時(shí)描繪出了下降段，但是同樣的設(shè)置，同樣的材料本構(gòu)，另一個(gè)模型卻 不能分析出下降段，這是什么原因呢？正如我上面提到的情況一樣，應(yīng)用力加載的方式，結(jié)構(gòu)是不穩(wěn)定的，不論你采用弧長(zhǎng)法或任何其他 方法，力是不斷增加的，而結(jié)構(gòu)已經(jīng)進(jìn)入軟

43、化段，無(wú)論理論上或數(shù)值上都不 可能收斂如果改為位移加載，即使結(jié)構(gòu)進(jìn)入軟化段，但弧長(zhǎng)法能夠收斂 到相應(yīng)的荷載加載模擬中如果出現(xiàn)荷載和位移都減小段，它只可能出現(xiàn)在結(jié)構(gòu)某一構(gòu)件中或者你控制結(jié)構(gòu)開(kāi)始卸載，而不是在加栽.UNL_Structuralwrote:不論你采用弧長(zhǎng)法或任何其他方法，力是不斷增加的至少對(duì)于弧長(zhǎng)法，力是可以增加或者減少的。UNL_Structural wrote:如果改為位移加載，即使結(jié)構(gòu)進(jìn)入軟化段，但弧長(zhǎng)法能夠收斂到相應(yīng)的荷載在abaqus中，弧長(zhǎng)法和位移加載不可能在同一步內(nèi)共存。meditera n wrote:我在做考慮結(jié)構(gòu)的壓柑失穩(wěn)影響的時(shí)候，也出現(xiàn)過(guò)類(lèi)似情況；當(dāng)改變一些

44、脆弱”桿件的本構(gòu)時(shí)，會(huì)呈現(xiàn)不同的XY值曲線，當(dāng)本構(gòu)較弱時(shí)，有一些情況是結(jié)構(gòu)的力，到了最大值后下降，然后按彈性 的E為斜率，位移不斷減小，這就是LZ說(shuō)的情況；事實(shí)是計(jì)算的問(wèn)題，當(dāng)所加荷載結(jié)構(gòu)不能承受時(shí)，有可能會(huì)出現(xiàn)兩種情況，一，結(jié)構(gòu)被壓屈，承載力下降，但位移仍然增大；二，結(jié)構(gòu)被壓屈，出現(xiàn)塑性的情況下卸載”即位移減小，位移-荷載曲線很類(lèi)似滯回曲線的形狀；這是我做算例所發(fā)現(xiàn)的情況，由于我對(duì)具體的算法及有限元理論了解有限，希望能有人解答。本貼與樓主提出的問(wèn)題無(wú)關(guān)，請(qǐng)開(kāi)辟新話題！參考論壇發(fā)帖規(guī)則第 13條,謝謝合作！編輯： walter2003顯然本貼和樓主提出的問(wèn)題是有關(guān)系的，mediteran的這段

45、話事實(shí)是計(jì)算的問(wèn)題，當(dāng)所加荷載結(jié)構(gòu)不能承受時(shí)，有可能會(huì)出現(xiàn)兩種情況，一，結(jié)構(gòu)被壓屈，承載力下降，但位移仍然增大；二，結(jié)構(gòu)被壓屈，出現(xiàn)塑性的情況下卸載”即位移減小，位移-荷載曲線很類(lèi)似滯回曲線的形狀”從某個(gè)角度回答了樓主的問(wèn)題，在abaqus里，用riks法分析彈塑性結(jié)構(gòu)，如果塑性加載但控制荷載減少的話，可能會(huì)得到彈性卸載的解。對(duì)我來(lái)說(shuō)，walter2003的判斷和處理難以理解。接受wo ng兄的批評(píng)，是我思考問(wèn)題不夠仔細(xì)。該貼打分已更改，歡迎wo ng兄多點(diǎn)評(píng)。謝謝！walter2003首先對(duì)wo ng兄的關(guān)注表示感謝，同時(shí)也感覺(jué)仁兄對(duì)riks有很深的認(rèn)識(shí)。我仍然有一個(gè)疑問(wèn)，實(shí)際加載并未達(dá)到我

46、所設(shè)置的上限。我瀏覽了一下幫助文件，發(fā)現(xiàn)適用的荷 載沒(méi)有g(shù)rav，所以我懷疑是我施加荷載的方式不對(duì)。但是現(xiàn)在還有一個(gè)問(wèn)題，我改用了集中荷載的方式加載，最后都不能收斂，錯(cuò)誤信息是：A zerodisplacement solution was found in the first iteration of a riks step. The an alysis cannot proceed. The most likely cause is that all displaceme nt degrees-of-freedom are con stra ined to be zero. 不知道是什么原

47、因造成 的。請(qǐng)wong兄以及各位大俠賜教！小弟謝謝了！感謝wo ng哥們兒替我出頭，呵呵，我想，斑竹每天需要看的帖子很多， 偶爾有點(diǎn)誤會(huì)也是可以理解的， 我就當(dāng)是對(duì)我的嚴(yán)格要求吧；建議Isk同學(xué)上傳inp文件，我想實(shí)際運(yùn)行一下， 這樣更穩(wěn)妥，也直觀。不知方便否；另外，合適的荷載沒(méi)有g(shù)rav ”不知道是什么意思， 是說(shuō)grav加載方式不適合做弧長(zhǎng)法嗎？ 我做過(guò)很多用grav的弧長(zhǎng)法的算例，都可以的； 不知道是我理解有誤還是其他原因；再說(shuō)句題外話，非線形有限元計(jì)算中， 計(jì)算結(jié)果不收斂真的是很讓人頭疼的事， 想成功地解決它，還需要大家多討論， 自己多感悟。謝謝理解，不過(guò)你可以自己PM給我，呵呵，我會(huì)

48、虛心接受批評(píng)的。再次感謝wo ng兄的指正。walter2003Won g:在abaqus里，用riks法分析彈塑性結(jié)構(gòu)，如果塑性加載但控制荷載減少 的話，可能會(huì)得到彈性卸載的解。這句話明顯有矛盾，塑性加載能得到彈性卸載的解，那么結(jié)構(gòu)到底是加載還 是卸載？#132007-5-28 23:52 一 一 一上-UNL_Structuralwrote:塑性加載能得到彈性卸載的解，那么結(jié)構(gòu)到底是加載還是卸載 我可沒(méi)說(shuō)過(guò)塑性加載能得到彈性卸載的解。塑性加載荷載減少是一個(gè)可能的解 彈性卸載荷載減少是另外一個(gè)可能的解 兩個(gè)解都可能符合 Riks迭代的約束條件。2007-5-29 02:29 一 一 一上-那

49、么你能否定義塑性加載荷載減少這到底是什么狀態(tài)呢？有算例或試驗(yàn)實(shí)現(xiàn)過(guò)這種狀態(tài)嗎#152007-5-29 02:37一 一 一上-UNL_Structuralwrote:那么你能否定義塑性加載荷載減少這到底是什么狀態(tài)呢？有算例或試驗(yàn)實(shí)現(xiàn)過(guò)這種狀態(tài)嗎？用你的話來(lái)回答UNL_Structural wrote:結(jié)構(gòu)已經(jīng)進(jìn)入軟化段#162007-5-29 02:57OK,那么在結(jié)構(gòu)進(jìn)入軟化段，它實(shí)際上存在兩個(gè)路經(jīng)（PATH）. 個(gè)是彈性卸載，一個(gè)是塑性加載. 如果是力加栽控制，實(shí)際上這兩個(gè)路經(jīng)無(wú)法區(qū)分 （所謂分叉）,理論上和程序 上都無(wú)法實(shí)現(xiàn).RIKS也不可能實(shí)現(xiàn)這一功能.2007-5-29 03:35

50、一 一 一上-UNL_Structural wrote:理論上和程序上都無(wú)法實(shí)現(xiàn).RIKS也不可能實(shí)現(xiàn)這一功能如果有解的話，檢測(cè)控制位移就可以區(qū)分了。#182007-5-29 03:59 一 一 一上-現(xiàn)有的有限元不可能是如此SMART吧.按照你的意思，在有限元計(jì)算，對(duì)應(yīng)與荷載加載位置，相應(yīng) 的位移要記錄下來(lái)，并和上一步比較況且就在當(dāng)前步中，非線性迭代也不 可能一步搜索到正確解，所以位移要上下調(diào)整，你如何確定它是加載還是 卸載其實(shí)在有限元計(jì)算中，應(yīng)該先確定是加載還是卸載，然后才能確定是 用相應(yīng)的彈性還是塑性剛度距陣，最后才能計(jì)算出相應(yīng)的位移#192007-5-29 04:06 一 一 一上-它

51、要不是如此SMART, 你又怎么會(huì)掏錢(qián)買(mǎi)呢:D2007-5-29 07:54mediteran 兄：關(guān)于 riks法中的荷載引用，幫助文件是這樣說(shuō)的：The followi ng loads can beprescribed in a Riks an alysis:Concen trated no dal forces can be applied to the displaceme nt degrees of freedom (1這個(gè)模型之前我做過(guò)動(dòng)力分析，所以里面有動(dòng)力計(jì)算的痕跡，另外，我在計(jì)算時(shí)支撐施加了初始 缺陷，為方便，在這里我把初始缺陷的設(shè)置處理掉了。我的支撐兩端設(shè)置 為鉸接，采用

52、了先前仁兄在另外一個(gè)帖子中提到的release的方法，確實(shí)很好使。我之前用的建立多個(gè)part，形成多個(gè)實(shí)體，之后在設(shè)置鉸的方式，不理想，因?yàn)楫?dāng)一個(gè)節(jié)點(diǎn)連接多個(gè)桿件時(shí)很容易過(guò)約束。我用現(xiàn)在的這個(gè)模型用riks計(jì)算時(shí)，還沒(méi)有算出下降段，請(qǐng)老兄指點(diǎn)！ ！718898-f24-j-3ji ngli.i np ( 206.5 K)下載次數(shù) 25我今天又用riks計(jì)算了一下，發(fā)現(xiàn)一個(gè)問(wèn)題，當(dāng)弧長(zhǎng)增量步取太大時(shí)不能收斂（正常），而弧長(zhǎng) 增量步太小時(shí)也不能收斂，錯(cuò)誤信息就是如同我前面發(fā)帖所說(shuō)。而增量步 適中時(shí)，可以收斂，但是有一點(diǎn)，我計(jì)算的這個(gè)模型算到一定程度又往回 算了。也許正如mediteran兄所提到的結(jié)構(gòu)剛度較小時(shí)可能出現(xiàn)塑性情況 下卸載”同時(shí)位移減小。請(qǐng)問(wèn)哪位仁兄可以指點(diǎn)一下，如何可以控制riks的計(jì)算使不出現(xiàn)往回算的情況呢？UNL_Structural wrote:其實(shí)在有限元計(jì)算中，應(yīng)該先確定是加載還是卸載，然后才能確定是用相應(yīng)的彈性還是塑