對(duì)ABAQUS中的umat和vumat的一些理解【轉(zhuǎn)發(fā)】

2017-06-15 by:CAE仿真在線來(lái)源:互聯(lián)網(wǎng)

1.
If UMAT is used with reduced-integration elements or shear flexible shell or beam elements, the hourglass stiffness and the transverse shear stiffness must be specified with the *HOURGLASS STIFFNESS and *TRANSVERSE SHEAR STIFFNESS options, respectively.

2.
1/The difficulty for passing SDV with UMAT or VUMAT is exactly the SAME ! Only slight differences ! The syntax is always : SDVNew=SDVOld+dSDV or SDVNew=function(blablabla) But you already know that !!!!! 2/There are no conveniences except using your little fingers to transform a UMAT to a VUMAT. 3/VUMAT is easier to code because it doesn't require the Jacobian computation DDSDDE(i,j). BenZ. --- In ABAQUS@xxxxxxxxxxxxxxx, "meyh2002" wrote: > Hi, Benz, > > What you said are hundred percent correct and I agree with them.I > disscuss the viscos progrm with you just wanting to get the hints of > how to pass and update variables in UMAT.I got confused because I used > UMAT instead of VUMAT, which seems easier to pass and update stress > and variable.Maybe I should change my UMAT into VUMAT.I will chek the > manual on the difference of the user subroutines. So far I think at > least I should change the integration scheme and user subroutine > interface. Is there a convenience way to do that? > > By the way, you really have a sharp head and are familar with ABAQUS? > How long have you been working on it, if you don't mind to tell? > Because sometimes I feel frustrated because BAQAUS is so complicated > that you can not understand the whole view of it.This make people > waste time and strength. > > > meyh ------------------------ Yahoo! Groups Sponsor --------------------~--> $9.95 domain names from Yahoo!. Register anything.http://us.click.yahoo.com/J8kdrA/y20IAA/yQLSAA/PMYolB/TM--------------------------------------------------------------------~->
Next Message by Date:
Re: problem during running UMAT

--- In ABAQUS@xxxxxxxxxxxxxxx, "benjamin_hagege" wrote: Thanks.I am working on it now.Hope I can get it done soon. meyh > 1/The difficulty for passing SDV with UMAT or VUMAT is exactly the > SAME ! Only slight differences ! The syntax is always : > SDVNew=SDVOld+dSDV or SDVNew=function(blablabla) > But you already know that !!!!! > 2/There are no conveniences except using your little fingers to > transform a UMAT to a VUMAT. > 3/VUMAT is easier to code because it doesn't require the Jacobian > computation DDSDDE(i,j). > > > BenZ. > > > --- In ABAQUS@xxxxxxxxxxxxxxx, "meyh2002" wrote: > > Hi, Benz, > > > > What you said are hundred percent correct and I agree with them.I > > disscuss the viscos progrm with you just wanting to get the hints of > > how to pass and update variables in UMAT.I got confused because I > used > > UMAT instead of VUMAT, which seems easier to pass and update stress > > and variable.Maybe I should change my UMAT into VUMAT.I will chek > the > > manual on the difference of the user subroutines. So far I think at > > least I should change the integration scheme and user subroutine > > interface. Is there a convenience way to do that? > > > > By the way, you really have a sharp head and are familar with > ABAQUS? > > How long have you been working on it, if you don't mind to tell? > > Because sometimes I feel frustrated

3.stateold是上一步的狀態(tài)變量,如果你沒(méi)有用到狀態(tài)變量,這個(gè)一直就是0,不是應(yīng)變。vumat和umat不一樣。
vumat本構(gòu)更新是這樣的:stressnew=stressold+D*straininc,不會(huì)自動(dòng)保存應(yīng)變
而umat里面是這樣的:stress=C*strain,應(yīng)變是作為一個(gè)變量保存起來(lái)的。
http://forum.simwe.com/forum.php?mod=viewthread&tid=1009357&highlight=VUMAT

4.LS-Dyna的本構(gòu)子程序,把輸入,輸出變量的名稱(chēng)改一下就直接可以用在Vumat上,不會(huì)有任何問(wèn)題,甚至你都可以用Matlab測(cè)試你的本構(gòu)子程序。
我對(duì)于您的問(wèn)題是這樣認(rèn)為的
1.首先是stressnew的問(wèn)題,對(duì)于小變形來(lái)說(shuō),基本是stress算到多少,后處理現(xiàn)實(shí)的就是多少,如果存在旋轉(zhuǎn)的話,那么可能兩者不一致,因?yàn)閼?yīng)力=應(yīng)力+應(yīng)力增量,這個(gè)是簡(jiǎn)單的公司、、公式,可是在vumat中你計(jì)算的只是這個(gè)應(yīng)力增量的一部分,還有一部分是共旋所引起的,這個(gè)你可以參考一下應(yīng)力客觀率的介紹,比如Jaumann率,這個(gè)時(shí)候應(yīng)力增量有三部分組成第一部份為變形導(dǎo)致的應(yīng)力增加,后面兩部分為旋轉(zhuǎn)導(dǎo)致的,這兩個(gè)部分應(yīng)力增量abaqus會(huì)在你計(jì)算完stressnew以后給你增加上。
2.一般對(duì)于使用vumat的還都到載荷作用的區(qū)域,我會(huì)選擇將這一位置設(shè)置為彈性體,這是個(gè)比較常用的方式,如果你的計(jì)算感興趣區(qū)域不再載荷作用位置。這個(gè)解釋起來(lái)很簡(jiǎn)單,圣維南原理解釋得很清楚,還有這個(gè)一般不會(huì)影響你的計(jì)算結(jié)果,如果你想問(wèn)為什么,這個(gè)其實(shí)可以開(kāi)個(gè)帖子來(lái)說(shuō),這里不詳細(xì)說(shuō)明了。
3.這個(gè)問(wèn)題就是個(gè)人觀點(diǎn)了,首先對(duì)于下降段出現(xiàn)的問(wèn)題,本身就是控制函數(shù)發(fā)生了變化,從橢圓問(wèn)題退化到雙曲問(wèn)題,我是只這個(gè)函數(shù)在空間上的控制方程的變化,這個(gè)時(shí)候一般會(huì)有震顫(因?yàn)榇嬖陂g斷解),但是存在震顫并不是人家abaqus計(jì)算的不對(duì),而是本來(lái)就應(yīng)該有震顫,好比流體激波行為,差分出來(lái)的解本來(lái)就是有震顫的。但是這個(gè)震顫可以消除就是通過(guò)數(shù)值耗散的方法,但是數(shù)值耗散其實(shí)是一中數(shù)值誤差,他將原本正確的震蕩解變?yōu)榱似交慕狻５沁@樣曲線看起來(lái)好看多了。至于為什么abaqu的只帶模型好使,就是因?yàn)槿思以黾恿藬?shù)值耗散。
這個(gè)問(wèn)題我在這個(gè)帖子中回了,你可以看看http://forum.simwe.com/viewthrea ... F0??&page=2
http://forum.simwe.com/forum.php?mod=viewthread&tid=995145&highlight=VUMAT

5.http://forum.simwe.com/forum.php?mod=viewthread&tid=1112102&highlight=VUMAT
VUMAT中如何調(diào)用或表示初始彈性剛度矩陣

————————————————————————————————————————
2015年1月25日周日

非線性靜態(tài)不穩(wěn)定問(wèn)題:
http://forum.simwe.com/forum.php?mod=viewthread&tid=1083142&highlight=不穩(wěn)定

"如果結(jié)構(gòu)局部變形不穩(wěn)定的話,我看到abaqus help里提到會(huì)施加阻尼"
關(guān)于這個(gè)前段時(shí)間我也碰到過(guò)這樣的問(wèn)題,兩個(gè)模型,有一個(gè)尺寸不同,改變了載荷后其中一個(gè)不收斂了,后來(lái)找別人解決的,就是添加了個(gè)穩(wěn)定因子,contact control ,,,,stabilization=1。就收斂了,結(jié)果還比較符合
后來(lái)查了些資料,說(shuō)是增加了阻尼,讓結(jié)構(gòu)盡快的穩(wěn)定,在這個(gè)過(guò)程中誤差小于默認(rèn)的值,才會(huì)收斂的。

本構(gòu)模型中需計(jì)算各個(gè)時(shí)刻積分點(diǎn)的真實(shí)應(yīng)變,也就是對(duì)數(shù)應(yīng)變 logarithmic strain。樓主說(shuō)的true strain 就是應(yīng)該是 logarithmic strain,所以沒(méi)有任何問(wèn)題的
與這些相對(duì)應(yīng)的應(yīng)該是engineering strain or nominal strain

-----------------------------
2015年1月25日周一陰
1.
靜力問(wèn)題選擇了explicit算法的求解時(shí)間預(yù)估
http://forum.simwe.com/forum.php?mod=viewthread&tid=1106152&highlight=Explicit
由于模型較復(fù)雜,接觸較多,雖然是靜力問(wèn)題,但是我選擇了explicit算法,模型的詳細(xì)參數(shù)如下:單元數(shù)量約30W,最小單元尺寸為1.4mm,為了保證單元質(zhì)量,沒(méi)辦法在增大了,我使用了質(zhì)量縮放,但是時(shí)間步長(zhǎng)也僅僅是2.8e-7,模型的加載為3個(gè)集中力、3個(gè)彎矩、1個(gè)重力載荷。我想問(wèn)問(wèn)有做過(guò)類(lèi)似案例的同仁,這樣的模型計(jì)算時(shí)間需要設(shè)定多久?我說(shuō)的是inp文件里面總的計(jì)算時(shí)間設(shè)定多少秒合適?我第一次采用explicit算法計(jì)算靜力問(wèn)題,沒(méi)經(jīng)驗(yàn),所以咨詢(xún)一下,謝謝!
自問(wèn)自答吧,1.首先求解出系統(tǒng)的第一階頻率(系統(tǒng)最低頻率),根據(jù)頻率算系統(tǒng)的周期,這就是準(zhǔn)靜態(tài)顯式求解的最小求解時(shí)間要求,如果條件允許,將求解時(shí)間放大10倍,求解絕對(duì)精準(zhǔn)。我已經(jīng)實(shí)踐過(guò)了,沒(méi)問(wèn)題。

2.
[高級(jí)討論區(qū)] 采用能量判斷Explicit擬靜力分析結(jié)果的方法總結(jié)及疑問(wèn)
http://forum.simwe.com/forum.php?mod=viewthread&tid=1117632&highlight=Explicit
在用Explicit開(kāi)展一個(gè)具有多個(gè)接觸關(guān)系模型的擬靜力分析,理想的情況當(dāng)然是使用implicit,但是如果接觸很多且接觸關(guān)系在不斷變化的彈塑性屈曲分析,implicit收斂就很困難,因此選用implicit,畢竟每一步不是力的平衡計(jì)算,總感覺(jué)需要對(duì)結(jié)果進(jìn)行嚴(yán)格的判別后才能使用,因此在論壇和Manual上找了一下相關(guān)準(zhǔn)則,現(xiàn)總結(jié)如下:(聲明:本帖內(nèi)容不屬于原創(chuàng),僅是總結(jié),論壇里已有很多大俠對(duì)這一問(wèn)題進(jìn)行過(guò)論述,對(duì)我?guī)椭艽?但由于看了很多帖子,就不一一注明出處了,在此一并謝過(guò)!)
分析類(lèi)型:具有多個(gè)接觸的彈塑性擬靜力分析,采用Abaqus/explicit模塊
1、Etotal應(yīng)保持常量,或近似常量(1%變化幅度);
問(wèn)題:我的模型始終沒(méi)有達(dá)到,一直為負(fù)且變化幅度大于1%,且在計(jì)算的某一時(shí)刻會(huì)突然變正,且持續(xù)增大到一個(gè)很大的量級(jí)(伴隨著ALLVD的突然增大),增加Linear Bulk Viscosity后,情況會(huì)好點(diǎn),能夠算更多的步數(shù),但是增加這個(gè)參數(shù)后對(duì)結(jié)果的影響還有疑問(wèn),是否應(yīng)該判斷ALLVD/ALLIE的比值,比如小于5%?
2、如果采用縮減積分單元,應(yīng)控制沙漏能,ALLAE/ALLIE應(yīng)控制在一定范圍(5%);
問(wèn)題:Manual里講ALLAE應(yīng)該和ALLSE或ALLKE相比是可以忽略的,那么這個(gè)可以忽略的應(yīng)該如何控制呢?我認(rèn)為應(yīng)該比5%有所放寬吧,畢竟ALLIE與ALLSE相比要大很多,另外當(dāng)擬靜力分析時(shí),ALLKE本來(lái)控制的就很小,此時(shí)不應(yīng)該再用ALLKE來(lái)控制了。
3、如是擬靜力分析,應(yīng)控制ALLKE/ALLIE小于1%,我覺(jué)得這一條還是相對(duì)容易達(dá)到的。
4、Manual上提到,ALLCW和ALLPW應(yīng)該接近于零;
問(wèn)題:我的模型中ALLCW可以保持始終為零,但是ALLPW卻為負(fù)值,且絕對(duì)值在增大,我認(rèn)為我的Etotal為負(fù)也可能與這一原因有關(guān)(論壇上說(shuō),Etotal為負(fù)的一個(gè)可能原因是沙漏能導(dǎo)致模型能量損失),有什么辦法可以減少ALLPW呢?
5、關(guān)于Manual里提到的Etotal應(yīng)當(dāng)恒定為一個(gè)常數(shù),我的理解是:
1)對(duì)于初始就有應(yīng)變能,動(dòng)能的系統(tǒng),Etotal的值應(yīng)該從分析開(kāi)始時(shí)就等于這一應(yīng)變能/動(dòng)能的值(從Manual的Benchmark模型的結(jié)果可以看出);
2)如果系統(tǒng)初始沒(méi)有能量存在,那么Etotal在分析過(guò)程中應(yīng)該接近于0;
3)關(guān)于實(shí)際計(jì)算時(shí),Etotal不能保持恒定,而是始終在變化,我猜想可能是數(shù)值計(jì)算引起的誤差(誤差產(chǎn)生后,其值可正可負(fù));
4)具體到我的模型,我認(rèn)為由于模型在分析過(guò)程中出現(xiàn)多次的屈曲,每一次屈曲可能帶來(lái)計(jì)算誤差的累積,造成Etotal的值在變化,每次屈曲后,Etotal都會(huì)在新的平衡值附近震蕩;
以上是我對(duì)采用能量輸出判斷Explicit進(jìn)行擬靜力分析結(jié)果的方法總結(jié)以及一些還未解決的問(wèn)題,尤其是第5條,屬于個(gè)人猜想,還沒(méi)找到理論依據(jù),希望各位高手指點(diǎn)迷津!

3.[土木] 素砼梁位移加載下降段-Stabilization/Riks/Explicit(呼喚大神@敦誠(chéng))
http://forum.simwe.com/forum.php?mod=viewthread&tid=1069815&highlight=Explicit

---------------------------------------
2015年1月30日

本構(gòu)模型有較為強(qiáng)烈的下降段時(shí),計(jì)算容易表現(xiàn)出部化效應(yīng)(損傷、塑性變形、或者應(yīng)變?cè)诰植繂卧?導(dǎo)致計(jì)算不穩(wěn)定,出現(xiàn)所謂網(wǎng)格依賴(lài)特性,就是網(wǎng)格劃分少時(shí)計(jì)算結(jié)果好,劃分多時(shí)表現(xiàn)出局部化效應(yīng)計(jì)算不穩(wěn)定,打個(gè)比方來(lái)說(shuō),對(duì)于一個(gè)truss單元,承受受壓荷載,如果其材料本構(gòu)有明顯下降段,在explicit分析中如果劃分1-3個(gè)單元,計(jì)算能穩(wěn)定性好,能計(jì)算出比較好的承載力曲線,如果劃分10個(gè)以上單元,在計(jì)算時(shí)構(gòu)件中的一個(gè)單元會(huì)率先進(jìn)入下降段,其壓縮應(yīng)變迅速增長(zhǎng),其它單元會(huì)出現(xiàn)拉伸(雖然此時(shí)構(gòu)件總體還是處于壓縮狀態(tài)),計(jì)算根本得不到下降段曲線,一旦越過(guò)最高點(diǎn),構(gòu)件的承載力就迅速降為0,理論上要通過(guò)所謂的梯度理論才能解決這個(gè)問(wèn)題,但梯度理論目前很不成熟,需要顛覆傳統(tǒng)有限元的單元構(gòu)建方式,使用困難,目前能較好解決方法就是如敦誠(chéng)版主版主所說(shuō)的引入粘滯。對(duì)于你的情況,建議你將模型中的損傷變量粘滯化使之能考慮應(yīng)變率效應(yīng),看效果是不是好一些。給一篇文獻(xiàn)希望對(duì)你有幫助。
[子程序] 損傷在vumat中的實(shí)現(xiàn)
損傷在vumat中的實(shí)現(xiàn)
http://forum.simwe.com/thread-966598-1-1.html
(出處: Simwe仿真論壇(forum.simwe.com),CAE/CAD/CAM/,FEA/FEM/有限元分析論壇---(手機(jī)驗(yàn)證注冊(cè)))

http://forum.simwe.com/forum.php?mod=viewthread&tid=966598&highlight=VUMAT

關(guān)于step,increment,attempt,iteration的關(guān)系的認(rèn)識(shí)
http://forum.simwe.com/thread-795512-1-1.html
(出處: Simwe仿真論壇(forum.simwe.com),CAE/CAD/CAM/,FEA/FEM/有限元分析論壇---(手機(jī)驗(yàn)證注冊(cè)))

step,increment,attempt,iteration,的關(guān)系
1)step 分析步
2)increment 時(shí)間增量步
3)attempt 減小增量步的嘗試,即“cutback”
4)iteration 迭代
在一個(gè)計(jì)算中有可能用到多步分析,比如建一個(gè)土石壩,每激活(add)一個(gè)填筑層就

是一個(gè)分析步step;
在每個(gè)step中,如果考慮非線性,step就會(huì)分成幾個(gè)增量步(increment)進(jìn)行計(jì)算;
在每個(gè)increment中,會(huì)有減小增量步的嘗試(attempt),在每個(gè)attemp中,要進(jìn)行

迭代計(jì)算(iteration)
如果迭代收斂,則在下一個(gè)increment中會(huì)增大時(shí)間增量步(比如第一個(gè)

increment=0.2,則下一個(gè)會(huì)增大為0.3)
如果迭代無(wú)法達(dá)到收斂,則ABAQUS會(huì)自動(dòng)減小時(shí)間增量步(減小increment),即

所謂的“cutback”,如果仍然不能收斂,則會(huì)繼續(xù)減小時(shí)間增量步,默認(rèn)的cutback

最大次數(shù)為5次,也就是attempt最大=5,如果5次之后仍不能收斂則ABAQUS會(huì)停止分

析,顯示錯(cuò)誤:too many attempts made for this increment:analysis

terminated.
increment時(shí)間增量步有最小值,默認(rèn)的是1e-5,如果increment減小到比這還小

,ABAQUS就會(huì)停止分析,出現(xiàn)錯(cuò)誤:time increment required is less than the

minimum specified.
increment的值可以在關(guān)鍵字*static中修改:
*static
1., 1., 1e-05, 1.
分別為初始增量步,分析時(shí)間步,最小增量步,最大增量步
可以用關(guān)鍵字*Step設(shè)定一個(gè)分析步中increment的最大步數(shù),如:
*Step,INC=600 (the maximum number of increments in a step,默認(rèn)的是100 )
*static和*Step中的increment是相同的,*Step,INC默認(rèn)為100,而*static中默認(rèn)為

1e-5,并不是100*(1e-5)=1,這兩個(gè)數(shù)都是限值,即number of increments最大為100

,而increment最小為1e-5。

關(guān)于explicit質(zhì)量放大的一點(diǎn)小心得
http://forum.simwe.com/thread-828207-1-1.html
(出處: Simwe仿真論壇(forum.simwe.com),CAE/CAD/CAM/,FEA/FEM/有限元分析論壇---(手機(jī)驗(yàn)證注冊(cè)))

本構(gòu)模型有較為強(qiáng)烈的下降段時(shí),計(jì)算容易表現(xiàn)出部化效應(yīng)(損傷、塑性變形、或者應(yīng)變?cè)诰植繂卧?導(dǎo)致計(jì)算不穩(wěn)定,出現(xiàn)所謂網(wǎng)格依賴(lài)特性,就是網(wǎng)格劃分少時(shí)計(jì)算結(jié)果好,劃分多時(shí)表現(xiàn)出局部化效應(yīng)計(jì)算不穩(wěn)定,打個(gè)比方來(lái)說(shuō),對(duì)于一個(gè)truss單元,承受受壓荷載,如果其材料本構(gòu)有明顯下降段,在explicit分析中如果劃分1-3個(gè)單元,計(jì)算能穩(wěn)定性好,能計(jì)算出比較好的承載力曲線,如果劃分10個(gè)以上單元,在計(jì)算時(shí)構(gòu)件中的一個(gè)單元會(huì)率先進(jìn)入下降段,其壓縮應(yīng)變迅速增長(zhǎng),其它單元會(huì)出現(xiàn)拉伸(雖然此時(shí)構(gòu)件總體還是處于壓縮狀態(tài)),計(jì)算根本得不到下降段曲線,一旦越過(guò)最高點(diǎn),構(gòu)件的承載力就迅速降為0,理論上要通過(guò)所謂的梯度理論才能解決這個(gè)問(wèn)題,但梯度理論目前很不成熟,需要顛覆傳統(tǒng)有限元的單元構(gòu)建方式,使用困難,目前能較好解決方法就是如敦誠(chéng)版主版主所說(shuō)的引入粘滯。對(duì)于你的情況,建議你將模型中的損傷變量粘滯化使之能考慮應(yīng)變率效應(yīng),看效果是不是好一些。給一篇文獻(xiàn)希望對(duì)你有幫助。
粘滯化不僅僅是為了考慮應(yīng)變加載速率效應(yīng)對(duì)本構(gòu)行為的影響,它能提高模型的穩(wěn)定性,在一定程度上消除網(wǎng)格依賴(lài),在考慮應(yīng)變率的本構(gòu)模型中如果將加載時(shí)間取的足夠長(zhǎng),則可認(rèn)為沒(méi)有考慮應(yīng)變率效應(yīng),但這樣處理卻能提高模型的穩(wěn)定性,部分的消除網(wǎng)格依賴(lài)和局部化效應(yīng)。

單元為1的時(shí)候邊界都是施加在這個(gè)單元上的,這樣就相當(dāng)于增加了單元的約束,提高了剛度,所以衰減行為可以比較緩慢的釋放,而多個(gè)單元模型中很多單元沒(méi)有這個(gè)優(yōu)勢(shì),所以還是快速衰減,并發(fā)生扭曲和失效,建議你定義一個(gè)單元?jiǎng)h除變量, 讓單元在damage達(dá)到一定數(shù)值的時(shí)候立刻刪除。同時(shí)在單元類(lèi)型里面設(shè)置一下比較大的體積粘滯系數(shù),看看還有沒(méi)有這種情況。
最后減小增量步!做兩個(gè)分析部一個(gè)讓材料基本接近損傷的起始,第二個(gè)分析部設(shè)置步長(zhǎng)很小,比如1e-9,然后讓他每一步一存儲(chǔ)。來(lái)看看損傷行為的發(fā)展是不是被explicit跳躍了。這樣做可能odb文件會(huì)很大,但是值得,我以前做過(guò),odb文件甚至達(dá)到過(guò)20G。
損傷在vumat中的實(shí)現(xiàn)
http://forum.simwe.com/thread-966598-1-1.html
(出處: Simwe仿真論壇(forum.simwe.com),CAE/CAD/CAM/,FEA/FEM/有限元分析論壇---(手機(jī)驗(yàn)證注冊(cè)))

-----------------------------------------------

2015年1月31日
對(duì)本構(gòu)方程積分的方法有exlicit and implicit,前者不需要迭代,運(yùn)算較快,但是它是有條件穩(wěn)定算法,而且精確度不高,還可能引起不正常卸載的情況發(fā)生,因此多半人們采用后者積分.對(duì)implicit算法的本構(gòu)方程是一個(gè)非線性方程組,用Newton-Raphson方法求解.這個(gè)就是叫做return mapping的算法,那個(gè)radial return method只是它的一個(gè)特列.
當(dāng)要求解軟化材料的問(wèn)題時(shí),剛度矩陣會(huì)變負(fù),現(xiàn)在的解決方法有引入附加條件,比如sphercial path control and arc-length control.參考:
1.Simo, J.C., Hughes, T.J.R., 1998. Computational Inelasticity. Springer-Verlag, New York.
2.O. C. Zienkiewicz and R. L. Taylor, The Finite Element Method
3.Crisfiled MA. Non-linear Finite Element Analysis of Solids and Structures. Volume 1: Essentials. Wiley:. Chichester, 1997
http://www.baisi.net/thread-680016-1-1.html

-------------------------------------------------

2015年3月10日
用顯示算法進(jìn)行動(dòng)力學(xué)分析時(shí),monitor里面的total energy一會(huì)負(fù)值一會(huì)正值,這個(gè)是為什么啊。還有那個(gè)kinetic energy和total energy 是什么意義,有什么關(guān)系?
kinetic energy是動(dòng)能;total energy 是總能量,
在進(jìn)行動(dòng)力學(xué)分析時(shí),一般來(lái)說(shuō)在整個(gè)計(jì)算過(guò)程中能量是守恒的,即計(jì)算過(guò)程中,total energy 應(yīng)該不變的,但是kinetic energy應(yīng)該是變化的,除非速度和質(zhì)量不變,同時(shí)可能在求解過(guò)程中存在摩擦能、沙漏能、外力做功、內(nèi)能與動(dòng)能的相互轉(zhuǎn)換等等,所以求解過(guò)程中可能會(huì)出現(xiàn)各種能量之間的互相轉(zhuǎn)換,你要全權(quán)監(jiān)控,total energy的變化是由什么引起的?可以肯定的是如果是有沙漏能引起的,那結(jié)果就是不可靠的。

轉(zhuǎn)自http://blog.sina.com.cn/s/blog_64fe8b9d0102vipx.html

開(kāi)放分享：優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才