ABAQUS-UMAT彈塑本構(gòu)二次開發(fā)的實(shí)現(xiàn)

/前言有限元法是工程中廣泛使用的一種數(shù)值計(jì)算方法。它是力學(xué)、計(jì)算方法和計(jì)算機(jī)技術(shù)相結(jié)合的產(chǎn)物。在工程應(yīng)用中,有限元法比其它數(shù)值分析方法更流行的一個(gè)重要原因在于：相對與其它數(shù)值分析方法，有限元法對邊界的模擬更靈活，近似程度更高。所以,伴隨著有限元理論以及計(jì)算機(jī)技術(shù)的發(fā)展，大有限元軟件的應(yīng)用證變得越來越普及。ABAQUS軟件一直以非線性有限元分析軟件而聞名，這也是它和ＡNSYS，Ｎasｔran等軟件的區(qū)別所在。非線性有限元分析的用處越來越大，因?yàn)樵谒貌牧戏浅?fù)雜很多情況下，用線性分析來近似已不再有效。比方說,一個(gè)復(fù)合材料就不能用傳統(tǒng)的線性分析軟件包進(jìn)行分析.任何與時(shí)間有關(guān)聯(lián),有較大位移量的情況都不能用線性分析法來處理。多年前,雖然非線性分析能更適合、更準(zhǔn)確的處理問題,但是由于當(dāng)時(shí)計(jì)算設(shè)備的能力不夠強(qiáng)大、非線性分析軟件包線性分析功能不夠健全,所以通常采用線性處理的方法。這種情況已經(jīng)得到了極大的改善，計(jì)算設(shè)備的能力變得更加強(qiáng)大、類似AＢAQUＳ這樣的產(chǎn)品功能日臻完善，應(yīng)用日益廣泛。非線性有限元分析在各個(gè)制造行業(yè)得到了廣泛應(yīng)用，有不少大型用戶.航空航天業(yè)一直是非線性有限元分析的大客戶，一個(gè)重要原因是大量使用復(fù)合材料。新一代波音７8７客機(jī)將全部采用復(fù)合材料。只有像AＢAQＵS這樣的軟件,才能分析包括多個(gè)子系統(tǒng)的產(chǎn)品耐久性能。在汽車業(yè),用線性有限元分析來做四輪耐久性分析不可能得到足夠準(zhǔn)確的結(jié)果.分析汽車的整體和各個(gè)子系統(tǒng)的性能要求（如懸掛系統(tǒng)等)需要進(jìn)行非線性分析。在土木工程業(yè)，ABAQＵS能處理包括混凝土靜動(dòng)力開裂分析以及瀝青混凝土方面的靜動(dòng)力分析，還能處理高度復(fù)雜非線性材料的損傷和斷裂問題，這對于大型橋梁結(jié)構(gòu)，高層建筑的結(jié)構(gòu)分析非常有效。瞬態(tài)、大變形、高級材料的碰撞問題必須用非線性有限元分析來計(jì)算.線性分析在這種情況下是不適用的.以往有一些專門的軟件來分析碰撞問題，但現(xiàn)在ＡＢAQUＳ在通用有限元軟件包就能解決這些問題.所以，AＢＡQUＳ可以在一個(gè)軟件完成線性和非線性分析。AＢAＱUS給用戶提供了強(qiáng)大二次開發(fā)接口，尤其是在材料本構(gòu)方面，給用戶開發(fā)符合實(shí)際工程的材料本構(gòu)模型提供了強(qiáng)大幫助,本文將針對其用戶材料子程序展開研究,總結(jié)常用材料模型的開發(fā)方法。目錄TOＣ\o＂1－3"\h\uHYＰERLINＫ＼ｌ”_Toc2３3０391７5”摘要 PAGEＲEF_Toc２330３9175\hIHＹPERLＩNＫ＼l＂_Ｔｏc2３3039176”ABSTＲACT?PＡGEＲEＦ＿Toc2３3039１７6\hＩIHYPERLＩNK1。1。?課題的研究背景?PAGEＲEＦ_Toc2３3039178\ｈ1HYＰERLINK\l”_Ｔoc23303９１79”1。２。 本文的研究內(nèi)容和方法 PAＧEREF_Ｔoｃ23３0３９１79\h2HYＰERLINK\l”＿Ｔｏc23303９180”2。 基于ABAQUS軟件的二次開發(fā)?PAGERＥＦ_Ｔoｃ２330391８0\h3HYＰERLＩＮK\l”_Toｃ2３303９１8１”２。1. AＢAQUＳ介紹?PＡＧEREF_Toc２3303９181＼h3HYPEＲＬINK\l＂_Toc23３0３９１82＂2。2.?AＢAQUS各模塊簡介 ＰAGERＥF_Toｃ２33039１８2＼h３2。4。?ABAＱUＳ的二次開發(fā)語言 PＡGEREF＿Toc233０３9184\h6ＨＹＰERLINＫ\l＂_Tｏc23３0３9185"3. 用戶材料子程序UMＡT ＰAGEREF_Ｔｏｃ2３30391８5\h８HＹPERＬINK＼l＂＿Tｏc233０3９186＂3.1．?ＵMＡT開發(fā)環(huán)境設(shè)置 PAGEＲEF_Toc２330３918６＼ｈ8HYPERLＩNK\l＂_Toｃ２３30３9187＂3.2。 ＵＭAT注意事項(xiàng) PAGＥＲEＦ＿Tｏc２33039187＼h9HYＰERLINK\l"_Ｔoｃ23３03９18８”3.3． UMAＴ接口的原理 ＰAGEＲEＦ＿Toc2330３9１8８\h10HYPＥRLINK＼l”_Toc2330３9189”３.4． UMＡＴ的使用方法?PAGEＲＥF_Toc２３３039189\ｈ1２HＹPERLＩNK＼l”_Toc２330３9190＂４．?材料非線性問題?ＰAＧEREF_Toc233039１90\h１４HYPERLＩＮK\l＂_Toc233039１９1"4.１．?材料的彈塑性本構(gòu)關(guān)系 ＰAＧＥREF＿Toc2330３919１\h14HYＰEＲLＩNＫ＼l＂＿Toｃ2３３0３９192”4.2。 非線性有限元算法理論 ＰAGＥREF_Tｏc2３3０3９1９2＼h17HＹPERLINK４。4．?增量理論切線剛度法公式推導(dǎo)?３303９194\h21ＨYPＥRＬIＮK5。1. 本構(gòu)關(guān)系描述 ＰＡGEＲEF_Tｏｃ23３03919６\h２5HＹPERLＩNK\ｌ＂＿Ｔoｃ２33０３９１9７"５.2． 常剛度法程序設(shè)計(jì)?PAGＥREF＿Tｏｃ23３039１97\h27HＹPERLＩNＫ＼l"＿Toc2330３9198＂5.3.?常剛度法程序編碼?PAＧEREF＿Toc２33０391９８\h２９ＨYＰEＲLINK\ｌ”_Toc2330391９9"5.4。?切線剛度法程序設(shè)計(jì) PAＧＥＲＥF_Tｏc233０39１99＼h32HYPERLＩNK\l＂_Toc233０3920０”5。5。?切線剛度法程序編碼 PAGEＲEF_Ｔｏc２3303920０＼h35HYPＥRLINK\l"_Ｔoc2３3０39201"5。6。?程序的調(diào)試 PAＧEＲEF_Toc233039２01＼h386．?程序驗(yàn)證?PAGEREF_Toｃ2３3039202\h４０HＹＰERLINK\l”_Ｔoｃ２330３９203”6．1.?問題描述?PAGＥREF_Toc23３0３920３\h40HＹＰEＲＬINK\ｌ”_Toc233０３9２04”6．2。?本構(gòu)關(guān)系 PＡGEREF_Tｏc２３303９２04\h41HYＰEＲLIＮK\l＂＿Toc2３３0３92０5"6。3。?ABＡQＵS自帶材料模型計(jì)算 PＡGEREF＿Toc23３039２05＼h4１ＨYＰERLINK\l"＿Toc２33039206"6。4.?常剛度法的ＵMAＴ驗(yàn)證 ＰＡGＥREＦ_Toｃ2３3039２06\ｈ4３HＹＰERLINK\ｌ”＿Ｔoｃ２３3０392０7”6.５． 切線剛度法的ＵＭAT驗(yàn)證?PAGERＥＦ＿Ｔoｃ233039207\h４5HＹPEＲＬIＮK\ｌ"_Tｏc2330３９208”6。6。?兩種算法的比較分析?PＡGEREＦ_Toc２３303920８\h47HYPＥRLＩNK\ｌ”_Toｃ2３３0３92０9＂7.?結(jié)論與展望?PAGEＲEF_Toｃ２33039２０9＼h51ＨYPＥRＬINK\l”＿Ｔoc2330３92１0"7。1。?結(jié)論?ＰＡＧEREF_Toc233039２10＼h51HYPERLINK\l”_Ｔoc23３0392１1"7．２．?展望 PＡGEREF_Toc２3３03921１\ｈ51ＨYＰERLＩNK參考文獻(xiàn)?PAＧEＲEF_Toｃ23３０39213＼h54ＨYＰERLIＮK＼ｌ"_Toｃ２33０39214”附1：ABＡQUS自帶彈塑性材料驗(yàn)證的INP文件?PＡGＥＲEF_Toc23303９214\h５５ＨYＰＥRＬＩＮK＼l＂_Toc233039215＂附2:用于算法驗(yàn)證的IＮＰ文件?ＰAGEREF_Toc23３03９215\h６1摘要ＡBＡQUS軟件功能強(qiáng)大,特別是能夠模擬復(fù)雜的非線性問題，它包括了多種材料本構(gòu)關(guān)系及失效準(zhǔn)則模型,并具有良好的開放性,提供了若干個(gè)用戶子程序接口，允許用戶以代碼的形式來擴(kuò)展主程序的功能.本文主要研究了ABAQUS用戶子程序UMAT的開發(fā)方法,采用ＦORTRAN語言編制了各向同性硬化材料模型的接口程序,研究該類材料的彈塑性本構(gòu)關(guān)系極其實(shí)現(xiàn)方法。本文緊緊圍繞UMAT的二次開發(fā)技術(shù),首先對其接口原理做了詳細(xì)介紹,然后針對非線性有限元增量理論中的常剛度法和切線剛度法的算法理論做了深入的剖析,推導(dǎo)出了常剛度法和切線剛度法的算法理論的具體表達(dá)式,然后分別編制了兩種算法的UMAT程序,最后建立了一個(gè)具體的驗(yàn)算模型，通過與AＢAＱUＳ自帶彈塑性本構(gòu)關(guān)系的計(jì)算結(jié)果相比較,驗(yàn)證兩者的正確性。本文還對常剛度法和切線剛度法得算法效率做了對比,得出了在非線性程度較高時(shí)切線剛度法效率高于常剛度法的結(jié)論。關(guān)鍵字：ABAQＵＳ、UMAＴ、有限元、材料非線性、FORＴＲAN、切線剛度ABSTＲACTＡＢＡQUＳｓoftwarepowerfuｌ，esｐｅcｉallyｔｏｓｉmｕlateｃomｐlexnon—lineaｒproｂleｍ,wｈｉｃhiｎcludｅｓaｗidｅｒangｅofmaterialcｏnｓtitｕtｉｖｅmodeｌａndfaiｌurecriteriａ，aｎdｈasaｇoｏdｏpen，providiｎganumbeｒｏｆuseｒsubroutineiｎｔｅｒfａｃethatａllｏwｓuserstocｏｄeformtoexｐandtheｆuｎｃｔionsｏｆｔheｍainpｒogram．ThispａperstudiｅstｈｅｕｓersubroｕtineUMＡToｆABAQUSdeveｌｏpmｅnｔmethods,tｈeｕsｅofFＯRＴＲAＮｌａngｕａgeisotrｏpｉcｈａrdeningmaterｉaｌmoｄeloftｈeｉnｔerｆａｃeprograｍ，stuｄiedtheeffecｔsoｆsuchｍaｔerialiseｘtremelｙｅlastｉc—ｐｌasticconstｉtutivｅｒeｌatiｏnmethod．ＴhisarticleUＭATtighｔlyaroundtheseconｄaryｄｅｖeloｐmｅntoｆtechｎoｌｏgy，thefｉrｓｔｐrincipleｏfitsｉnterfａｃedetaｉl，andthenforｔhｅtheｏryofnonｌineａrｆiniｔｅｅlemeｎｔiｎcrｅｍeｎtalｓｔiffneｓsofｔｈｅregｕlａrtａngenｔstiffnessｍetｈｏdaｎdｔhｅthｅoryofaｌgorithmstodｏaniｎ-ｄeptｈａnalyｓｉsoｆｄeｄucｅdaregulartangentsｔiffneｓｓanｄｒiｇidityoｆtｈｅlａwoftｈesｐecｉfｉceｘpｒessionofaｌgorithmtｈeoｒy，andｔｈｅｎthepｒｅpaｒａt(yī)ionofｔhetwoａlgoｒithms，respeｃtively，oftheUMATｐrｏgram，andfｉnａllｙtｈeestａｂlisｈmeｎｔoｆａsｐecificmoｄｅlcheckinｇ，bringingwithＡBAQUSeｌasto—ｐlaｓticconstitutiverｅlaｔｉonoftｈeｃalculatedｒesultscｏmpａredｔoｖerｉfythｅcorrecｔｎｅssoｆtｈｅｔwｏ.Tｈisarticleａｌsｏｏfｔeｎｓtiffnessandｔaｎgｅntｓtiｆfnessmethodwastodoａｃoｍpａrｉsoｎofalｇorｉthｍｅｆficiencｙisobtａinedwhenａｈｉgherｄegreeintｈeｎoｎ-ｌineａrtangｅｎｔstiｆfｎessｍeｔｈｏｄmoｒeｅffｉcｉeｎttｈａnthecoｎclusionsoflawoftensｔifｆｎess.ＫEYWORDＳ:ABAＱUS、UMAT、Ｆiｎiｔｅe(cuò)ｌemｅnｔ、Materiaｌｎoｎlinearitｙ、ＦORTRAN、Tanｇeｎtsｔｉffｎｅsｓ緒論課題的研究背景有限單元法基本思想的提出,可以追溯到克勞夫(R．W。Ｃloｕgh)在１9４3年的工作［1]，他第一次嘗試應(yīng)用定義在三角形區(qū)域上的分片連續(xù)函數(shù)和最小位能原理相結(jié)合，來求解Sｔ。Venant扭轉(zhuǎn)問題。1960年克勞夫進(jìn)一步處理了平面彈性問題，并第一次提出了“有限單元法”的名稱,使人們開始認(rèn)識了有限單元法的功效。四十多年來，隨著電子計(jì)算機(jī)的廣泛應(yīng)用和發(fā)展，有限單元法的理論和應(yīng)用都得到迅速的，持續(xù)不斷的發(fā)展，其應(yīng)用己由彈性力學(xué)平面問題擴(kuò)展到空間問題、板殼問題，由靜力學(xué)問題擴(kuò)展到穩(wěn)定問題、動(dòng)力問題和波動(dòng)問題。分析的對象從彈性材料擴(kuò)展到塑性、粘彈性、粘塑性和復(fù)合材料等，從固體力學(xué)擴(kuò)展到流體力學(xué)、傳熱學(xué)等連續(xù)介質(zhì)力學(xué)領(lǐng)域。在工程分析中的作用已從分析和校核擴(kuò)展到優(yōu)化設(shè)計(jì)并和計(jì)算機(jī)輔助設(shè)計(jì).利用有限元軟件解決工程和科學(xué)問題，是有限元理論應(yīng)用于工程設(shè)計(jì)和科學(xué)研究實(shí)踐的主要形式。由于工程設(shè)計(jì)的巨大市場需要,有限元軟件的發(fā)展是很迅速的，目前常用的大型有限元軟件常見的有Sap2００0,AＤINＡ，ＭSC/NASTRAN，MＳCMarc,ANSＹS，ＡBAQUS等，這些軟件的共同特點(diǎn)是具有豐富的單元庫和求解器，強(qiáng)大而可靠的分析功能,人們利用這些軟件解決了很多工程建設(shè)和工業(yè)產(chǎn)品設(shè)計(jì)中遇到的問題，取得了巨大的經(jīng)濟(jì)技術(shù)效益。由于工程問題的千差萬別，不同的用戶有不同的專業(yè)背景和發(fā)展方向，通用軟件不免在具體的專業(yè)方面有所欠缺，針對這些不足,大部分的通用軟件都提供了二次開發(fā)功能，以幫助用戶減少重復(fù)性的編程工作、提高開發(fā)起點(diǎn)、縮短研發(fā)周期、降低開發(fā)成本,并能簡化后期維護(hù)工作，給用戶帶來很多方便?；谕ㄓ密浖脚_進(jìn)行開發(fā)，是目前研究的一個(gè)重要發(fā)展方向.ＡBＡQＵＳ也提供了若干用戶子程序（ＵserＳubrｏｕtｉneｓ）接口，它是一個(gè)功能非常強(qiáng)大且適用的分析工具，與命令行的程序格式相比，用戶子程序的限制少得多,從而使用更加靈活方便.針對ABAQUS所提供的本構(gòu)關(guān)系模型種類有限，無法滿足工程應(yīng)用需要的問題,用戶子程序中的用戶材料子程序（User—definedMateriaMecｈanｉcalBehavｉoｒ，簡稱UMAT)接口可以幫助用戶定義自己的材料本構(gòu)模型和算法，這是AＢAQUS的獨(dú)到之處。由于其操作方便,能被靈活地應(yīng)用于各個(gè)領(lǐng)域中，尤其受到用戶的青睞。本文的研究內(nèi)容和方法ABAQUS中用戶材料子程序UＭＡＴ的開發(fā)主要解決兩方面的問題：本構(gòu)模型的建立和積分算法的選擇。本文主要研究非線性材料的ＵMＡT實(shí)現(xiàn)方法，并重點(diǎn)研究其迭代算法部分，目前，用戶材料子程序UMＡT的迭代算法主要是常剛度法，常剛度法的優(yōu)點(diǎn)在于算法原理較簡單，程序編寫較方便,缺點(diǎn)是當(dāng)遇到復(fù)雜非線性材料時(shí)，其迭代次數(shù)較多,收斂速度也較慢，在這個(gè)情況下，本文采取的是一種迭代次數(shù)較少且收斂速度較快的切線剛度法，具體就是采用ＦORTRＡN語言編制了基于Ｖoｎ—Miseｓ模型的接口程序,并采用切線剛度算法，通過與AＢAQUS自帶本構(gòu)關(guān)系計(jì)算的結(jié)果相比較，驗(yàn)證其正確性.本文的研究工作緊緊圍繞UMAT的二次開發(fā)技術(shù)，首先根據(jù)有限元方法推導(dǎo)材料非線性問題算法的公式,然后參考UＭAＴ接口規(guī)范設(shè)計(jì)程序的算法流程,繼而編寫出該程序,最后建立一個(gè)具體的本構(gòu)和具體的模型做測試，驗(yàn)證程序的正確性，在這一過程中，調(diào)試是一個(gè)非常重要的過程，占用了大量的時(shí)間，在調(diào)試程序時(shí)采用了將中間變量輸出到文本的方式，這樣能明確跟進(jìn)迭代過程,發(fā)現(xiàn)算法或程序的缺陷。本文采用的本構(gòu)關(guān)系是經(jīng)過歸納和抽象的,也就是說本文的程序并不僅僅是只針對某個(gè)具體模型和問題,而是針對所有符合抽象出的各向同性硬化材料，這樣做的好處是能保證程序的通用性和復(fù)用性,避免以后的重復(fù)勞動(dòng)，當(dāng)然,這也是符合ＡBAＱUS軟件設(shè)計(jì)UＭAT接口的宗旨的.基于ABAQUS軟件的二次開發(fā)ＡBＡＱＵS介紹ABAQUS是一套功能強(qiáng)大的基于有限元法的工程模擬軟件［2］，其解決問題的范圍從相對簡單的線性分析到最富有挑戰(zhàn)性的非線性模擬問題.ABAＱUS具備十分豐富的、可模擬任意實(shí)際形狀的單元庫。并與之對應(yīng)擁有各種類型的材料模型庫，可以模擬大多數(shù)典型工程材料的性能，其中包括金屬、橡膠、高分子材料、復(fù)合材料、鋼筋混凝土、可壓縮彈性的泡沫材料以及巖石和土這樣的地質(zhì)材料。作為通用的模擬分析工具,ABAQUS不僅能解決結(jié)構(gòu)分析中的問題,還能模擬和研究各種領(lǐng)域中的問題，如熱傳導(dǎo)、質(zhì)量擴(kuò)散、電子元器件的熱控制(熱一電耦合分析)、聲學(xué)分析、土壤力學(xué)分析(滲流——應(yīng)力耦合分析）和壓電介質(zhì)力學(xué)分析。ABＡQUS為用戶提供了廣泛的功能，且使用起來又十分簡明。最復(fù)雜的問題也可以很容易地建立模型[3]。例如復(fù)雜的多部件問題可以通過對每個(gè)部件定義材料模型和幾何形狀，然后再把它們組裝起來而構(gòu)成.在大部分模擬分析問題中,甚至在高度非線性問題中，用戶也只需要提供結(jié)構(gòu)的幾何形狀、材料性能、邊界條件和荷載工況這樣的工程數(shù)據(jù)就可以進(jìn)行分析。在非線性分析中，ＡＢＡQUＳ能自動(dòng)選擇合適的荷載增量和收斂精度。不僅能選擇這些參數(shù)值,而且能在分析過程中不斷地調(diào)整參數(shù)來保證有效地得到高精度的解，很少需用戶去定義這些參數(shù).ＡBAQUS各模塊簡介ＡＢAＱUＳ有兩個(gè)主要的分析模塊:ＡBAQUS/Standarｄ和ABAQＵS/Explicit。ＡＢAQＵS/Staｎdard還有兩個(gè)特殊用途的附加分析模塊：ＡＢAQUS/Aqua和ＡＢAQUS/Ｄeｓigｎ。另外，還有ＡBAQUS分別與AＤAMS/Flex，C－MＯLＤ和Moldfｌow的接口模塊:ABＡQＵS／ADAMＳ，ABAQＵS/C—MＯLD和ABAＱＵＳ/MOLDFLOW。ABAQUS/CAE是完全的ABＡQUＳ工作環(huán)境模塊,它包括了ABAＱUS模型的構(gòu)造，交互式提交作業(yè)、監(jiān)控作業(yè)過程以及評價(jià)結(jié)果的能力.ABAＱＵS/Vｉewer是ABAQUS/CAE的子集，它具有后處理功能，這些模塊之間的關(guān)系見圖２－SEQ圖2－\＊ＡRABＩC1圖2－１AＢAQUＳ/StandａrdAＢAQUS/Standａrｄ是一個(gè)通用分析模塊,在數(shù)值方法上采用有限元方法常用的隱式積分.它能夠求解廣泛的線性和非線性問題，包括結(jié)構(gòu)的靜態(tài)、動(dòng)態(tài)問題、熱力學(xué)場和電磁場問題等.對于通常同時(shí)發(fā)生作用的幾何、材料和接觸非線性可以采用自動(dòng)控制技術(shù)處理，也可以由用戶自己控制。ABAQUS/ExpｌｉciｔABAQUS/Explicit是一個(gè)在數(shù)值方法上采用有限元顯式積分的特殊模塊，它利用對時(shí)間的顯式積分求解動(dòng)態(tài)有限元方程.它適合于分析諸如沖擊和爆炸這樣短暫、瞬時(shí)的動(dòng)態(tài)問題，同時(shí)對高度非線性問題如模擬加工成型過程中接觸條件的改變等也非常有效。ABAQＵS/ＣAＥABAＱＵS／CＡＥ是ABAQUS進(jìn)行有限元分析的前后處理模塊，也是建模、分析和后處理的人機(jī)交互平臺.該模塊根據(jù)結(jié)構(gòu)的幾何圖形生成網(wǎng)格,將材料和截面的特性分配到網(wǎng)格上，并施加載荷和邊界條件。該模塊可以進(jìn)一步將生成的模型投入到分析模塊中進(jìn)行高效率的后臺運(yùn)行，并對運(yùn)行情況進(jìn)行監(jiān)測，對計(jì)算結(jié)果進(jìn)行后處理。AＢＡＱUS/CＡE的后處理支持ABAＱＵS分析模塊的所有功能,并且對計(jì)算結(jié)果的描述和解釋提供了范圍很廣的選擇,除了通常的云圖，等值線和動(dòng)畫顯示之外，還可以用列表，曲線（包括部分常用運(yùn)算）等其他常用工具來完成對結(jié)果數(shù)據(jù)的處理.該模塊的許多獨(dú)特功能與特點(diǎn)，例如CＡD特征化建模、參數(shù)化建模、適應(yīng)設(shè)計(jì)者要求的數(shù)據(jù)管理系統(tǒng)等極大的方便了AＢAQUS的使用者。ＡBAQUS/AquaABAQＵS／Ａqｕa的一系列功能可以附加在AＢAQＵS/Stanｄard中應(yīng)用.它偏向于模擬海上結(jié)構(gòu)，如海洋石油平臺。它的功能包括模擬波浪，風(fēng)載荷及浮力的影響。在本指南中不討論ABAＱUS/Ａquａ。ＡBAQＵS／ADAMSABAQUS/ADAMS允許ABAQＵS有限元模型作為柔性部件進(jìn)入到MDIＡDＡMS產(chǎn)品族中去進(jìn)行分析。AＢAQUS／Ｃ-MOLDAＢAＱUS/Ｃ－MOLD把注模分析軟件C－ＭＯＬD中有限元網(wǎng)格、材料性質(zhì)和初始應(yīng)力數(shù)據(jù)轉(zhuǎn)換成為ABAQUS輸入文件。ABＡＱUS/DｅｓｉｇｎA(yù)BAQＵS/Deｓiｇn的一系列功能可附加在AＢAＱＵS／Sｔandaｒｄ中進(jìn)行設(shè)計(jì)敏度計(jì)算。ABAＱUＳ/ＭＯLDＦLOＷＡBAQUS/MＯLDＦＬOＷ模塊把MOＬDFLOW分析軟件中的有限元模型信息轉(zhuǎn)換成ＡＢＡQUVS輸入文件的一部分.ABAＱUS的二次開發(fā)平臺ABAQＵＳ的腳本語言接口非常友好，其自嵌的腳本語言是Pｙthｏn[4］,系國際上廣泛使用、功能強(qiáng)大、具有良好開放性的一種面向?qū)ο蟪绦蛟O(shè)計(jì)語言。所以,應(yīng)用Pyｔhｏn在AＢＡQUS中進(jìn)行二次開發(fā)也比較方便,且可移植性強(qiáng)。ＡＢAQUS以基于Pytｈon的語法規(guī)則向二次開發(fā)者提供了許多庫函數(shù)，這些庫函數(shù)主要是用來增強(qiáng)ABAQＵS的交互式(GUI）操作功能。用戶可以通過AＢAＱＵS的交互式（GUI)界面實(shí)現(xiàn)分析對象的特征造型、指定材料屬性、完成網(wǎng)格剖分和控制、提交并監(jiān)控分析作業(yè)，也可以使用ABAＱUＳ腳本語言越過ABＡQUS的交互式(GUI）界面直接高效地向ＡＢAQUS內(nèi)核提交任務(wù)。使用Python可以進(jìn)行參數(shù)化建模,修改交互式建立的模型，還可以一次提交多個(gè)作業(yè)。出了腳本語言接口,AＢAQUS還為用戶提供了功能強(qiáng)大的用戶子程序接口(AbaqｕsＵserSubroutines），以幫助用戶開發(fā)基于ABAQUＳ內(nèi)核的程序，常用的用戶子程序包括ＵEL(Ｕsersubrｏｕtineｔodefiｎｅanｅlemｅnt,用戶單元子程序），UＭAT（Uｓｅｒsubroutiｎetｏｄｅｆｉｎeａmａt(yī)ｅｒｉａl'sｍechanicaｌｂehavioｒ,用戶材料子程序）［5]，其中UMAT的使用最為廣泛，它主要用于用戶開發(fā)自己的材料模型，以彌補(bǔ)ABAQＵS自帶材料模型的不足，幫助用戶完成各種材料分析，功能極為強(qiáng)大.在國外,眾多的有限元分析和研究者熱衷于使用ＡBAQＵS，一個(gè)很重要的原因就在于AＢAQUS給用戶提供了功能強(qiáng)大，使用方便的二次開發(fā)工具和接口,使得用戶可以方便的進(jìn)行富含個(gè)性化的有限元建模、分析和后處理，滿足特定工程問題的需要.通過用戶材料子程序接口，用戶可定義任何補(bǔ)充的材料模型,不但任意數(shù)量的材料常數(shù)都可以作為資料被讀取，而且AＢAQUS對于任何數(shù)量的與解相關(guān)的狀態(tài)變量在每一材料計(jì)數(shù)點(diǎn)都提供了存儲功能，以便在這些子程序中應(yīng)用.ＡBＡQＵＳ的二次開發(fā)語言ABＡＱＵS的二次開發(fā)語言主要有３種：Pｙthoｎ，ＦORTRAＮ，C＋+Pythｏn語言主要用于ＧＵI開發(fā),FＯRTＲAN語言主要用于用戶子程序開發(fā)，而c++語言主要專注于其他高級開發(fā)部分.本文主要是針對用戶子程序的開發(fā),所以采用FORＴRAN語言,下面簡要介紹一下該語言極其特點(diǎn)：FORＴRAＮ語言是世界上第一個(gè)被正式推廣使用的高級語言［6].它是１95４年被提出來的，１956年開始正式使用,至今已有三十多年的歷史，但仍歷久不衰,它始終是數(shù)值計(jì)算領(lǐng)域所使用的主要語言。ＦORＴRAN語言是FormulaTrａｎｓｌaｔion的縮寫，意為“公式翻譯”.它是為科學(xué)、工程問題或企事業(yè)管理中的那些能夠用數(shù)學(xué)公式表達(dá)的問題而設(shè)計(jì)的，其數(shù)值計(jì)算的功能較強(qiáng)。FORTRＡN語言問世以來,根據(jù)需要幾經(jīng)發(fā)展，先后推出了不同的版本,主要版本有FORTRAN77,FOＲＴRAN9０,FOＲＴRAN9５，ABAQUＳ采用FOＲTＲAＮ77,通常用固定格式編寫代碼.ＦOＲTＲAＮ7７語言同C語言一樣，是一種結(jié)構(gòu)化編程語言結(jié)構(gòu)化程序設(shè)計(jì)方法規(guī)定，在結(jié)構(gòu)化的程序中，只能有三種基本結(jié)構(gòu):（1)順序結(jié)構(gòu)這是一種最簡單的基本結(jié)構(gòu)形式，它的特點(diǎn)是,在這個(gè)結(jié)構(gòu)內(nèi)的各個(gè)功能模塊或語句序列，是按其出現(xiàn)的先后順序執(zhí)行的，如賦值語句、輸入／輸出語句等。它有一個(gè)入口和一個(gè)出口,并在入口和出口之間包含著若干個(gè)功能塊，其中每一個(gè)功能塊可以是一個(gè)非轉(zhuǎn)移語句.因此，順序基本結(jié)構(gòu)塊是由一系列的順序執(zhí)行語句組成的。（2）分支選擇結(jié)構(gòu)在給定的條件下,分支選擇結(jié)構(gòu)判斷選擇哪一條路徑執(zhí)行，不同路徑完成的功能是不同的。實(shí)現(xiàn)分支選擇結(jié)構(gòu)主要由塊ＩF語句、ELSE語句、ENＤIＦ語句以及ELSEIF語句組成的IF－ＴＨEN—ELSE結(jié)構(gòu)。（3）循環(huán)結(jié)構(gòu)循環(huán)結(jié)構(gòu)也稱重復(fù)處理結(jié)構(gòu)，即重復(fù)執(zhí)行某一功能塊,直到滿足（或不滿足)某一條件為止。實(shí)現(xiàn)循環(huán)結(jié)構(gòu)的FORTＲAN90語句主要是DO語句、塊IF語句和邏輯IF語句的結(jié)合。以上三種基本結(jié)構(gòu)，是組成結(jié)構(gòu)化程序的基本結(jié)構(gòu)形式.這里有兩層意思：一是結(jié)構(gòu)化的程序中,各個(gè)模塊均由這三種基本結(jié)構(gòu)組成；二是結(jié)構(gòu)化程序本身，從宏觀上也是這三種基本結(jié)構(gòu)形式之一.用戶材料子程序UMＡTUMAT開發(fā)環(huán)境設(shè)置由于UMＡT是采用FORTRＡＮ語言編寫,那么要運(yùn)行UMAT就需要安裝FＯRＴRAN的開發(fā)環(huán)境，同時(shí)還需要ABＡQUＳ的支持，本文采用的ABAQUS版本為6。81，支持ＩNTELFｏｒtｒaｎ9．１-１0。1,IｎtelＦｏｒtran安裝時(shí)又需要安裝MicrosofｔViｓualＳtudio的相應(yīng)版本,經(jīng)過比較,本文選用ABAQUS６。８1+InteｌFoｒtran10。1+MicｒoｓoftＶisuａlＣ++２０05,相對于ABAQUS來說,UＭAＴ開發(fā)環(huán)境的設(shè)置較為繁瑣，這給子程序的使用帶來諸多不便，為了解決這一問題,我用C＃語言編制了AＢAＱＵS子程序編譯環(huán)境設(shè)置工具，只需要將安裝文件解壓到ABAQＵＳ的安裝目錄，運(yùn)行安裝程序就可以了，整個(gè)過程不需要人工干預(yù)，也不需要安裝龐大的ＶisuａlC++20０５，如圖3—1所示圖3—1UMＡT注意事項(xiàng)ＡBAＱUS的用戶子程序是根據(jù)ＡBAQUS提供的相應(yīng)接口，按照Ｆortraｎ語法，用戶自己編寫的代碼。它是一個(gè)獨(dú)立的程序單元,可以獨(dú)立的被存儲和編譯，也能被其它程序單元引用,因此，利用它可帶回大量數(shù)據(jù)供引用程序使用,也可以用它來完成各種特殊的功能。它的一般結(jié)構(gòu)形式是:SUBROＵTINES(x1，x2,……，xｎ)INCＬＵDＥ‘ABA_PＡRAM.INC’（用于ABAQUS／Ｓtaｎｄaｒｄ用戶子程序中)OＲＩNCLUＤE‘VABＡ_ＰARＡＭ.INC'）（用于ABＡQUS／Expｌiｃit用戶子程序中）……REＴURNEＮＤx１,ｘ２，……，xn是ABAQＵＳ提供的用戶子程序的接口參數(shù)，有些參數(shù)是ABAＱＵＳ傳到用戶子程序中的，例如SUBRＯUTIＮEDLOAD中的ＫSTEP、KINC、COORＤＳ,有些是需要用戶自己定義的，例如F,文件abａ＿paｒａｍ.inc和vａba_param。inｃ隨著ABAQUS軟件的安裝而包含在操作系統(tǒng)中，它們含有重要的參數(shù)，幫助AＢAQUS主求解程序?qū)τ脩糇映绦蜻M(jìn)行編譯和鏈接。當(dāng)控制遇到RETUＲN語句時(shí)便返回到引用程序單元中去，ENＤ語句是用戶子程序結(jié)束的標(biāo)志.在一個(gè)算例中，用戶可以用到多個(gè)用戶子程序，但必須把它們放在一個(gè)以．for為擴(kuò)展名的文件中。運(yùn)行帶有用戶子程序的算例同時(shí)有兩種方法：一是在CAE中運(yùn)行,在EDIＴJＯB菜單中的GENERＡL子菜單的UＳERSUＢROUTＩNEFILE對話框中選擇用戶子程序所在的文件即可；另外是在MAND中運(yùn)行,語法如下：aｂaquｓｊob＝ｊoｂ-naｍeuseｒ=｛ｓｏurｃe—｝編制用戶子程序時(shí)應(yīng)注意(１）用戶子程序相互之間不能調(diào)用,但可以調(diào)用用戶自己編寫的Fortｒan子程序和ABＡQUＳ應(yīng)用程序.AＢAＱＵS應(yīng)用程序必須由用戶子程序調(diào)用。當(dāng)用戶編寫Fｏｒｔrａn子程序時(shí)，建議子程序名以Ｋ開頭，以免和AＢAQUS內(nèi)部程序沖突。(2)當(dāng)用戶在用戶子程序中利用OPEN打開外部文件時(shí)，要注意以下兩點(diǎn):一是設(shè)備號的選擇是有限制的，只能取15~１8和大于100的設(shè)備號，其余的都已被ABAＱUS占用;二是用戶需提供外部文件的絕對路徑而不是相對路徑。(3）對于不同的用戶子程序ABAQUＳ調(diào)用的時(shí)間是不同的，有的是在每個(gè)ＳＴEＰ的開始,有的是STEＰ的結(jié)尾，有的是在每個(gè)ＩNCREＭＥＮT的開始等等.當(dāng)ABAＱＵS調(diào)用用戶子程序時(shí),都會把當(dāng)前的STEＰ和IＮＣＲＥＭENＴ利用用戶子程序的兩個(gè)實(shí)參ＫＳTEＰ和KＩNC傳給用戶子程序，用戶可把它們輸出到外部文件中，這樣就可清楚的知道ＡBAQUＳ何時(shí)調(diào)用該用戶子程序。為保證用戶子程序的正確執(zhí)行,子程序的書寫必須遵循ＡBＡQUS的相關(guān)規(guī)定，下面以用戶材料子程序?yàn)槔敿?xì)說明。UMAT接口的原理用戶材料子程序（Usｅr-dｅfｉnedMaｔeriａlMechaniｃaｌBｅhａｖiｏｒ,簡稱UMAT）是ABAＱUS提供給用戶定義自己的材料屬性的Ｆoｒtｒan程序接口［７］［８],使用戶能使用AＢAQＵS材料庫中沒有定義的材料模型。用戶材料子程序ＵMＡT通過與ABAQUS主求解程序的接口實(shí)現(xiàn)與ＡBＡQUＳ的資料交流.在輸入文件中,使用關(guān)鍵詞“＊ＵSＥRMATEＲIAL"表示定義用戶材料屬性。UMAT子程序具有強(qiáng)大的功能，使用UMAＴ子程序：(1）可以定義材料的本構(gòu)關(guān)系，使用ABAQUS材料庫中沒有包含的材料進(jìn)行計(jì)算，擴(kuò)充程序功能。(2）幾乎可以用于力學(xué)行為分析的任何分析過程，幾乎可以把用戶材料屬性賦予ABAＱUS中的任何單元。（3）必須在ＵMAT中提供材料本構(gòu)的雅可比(Jａcｏbian)矩陣，即應(yīng)力增量對應(yīng)變增量的變化率。由于主程序與UMＡＴ之間存在數(shù)據(jù)傳遞,甚至共享一些變量,因此必須遵守有關(guān)UMAT的書寫格式，UMＡT中常用的變量在文件開頭予以定義，通常格式SUBROUＴINEUMAT（ＳTREＳS，ＳTＡTEV，DＤSＤDE,SSＥ,SPD,SCD,1ＲPL，ＤDＳＤDT，DＲPLDE，DRＰLDT，2ＳTRAN，DSTＲAN，TＩME，DＴIME,TEMP，DTEMＰ,ＰRＥDEF,DPRＥD，CMNAＭＥ,３NDI，NSＨR,NTENS，ＮＳTATV,PROPS，ＮＰROPS,CＯORDＳ,DROT，PNEWDT，４CELENＴ，ＤFＧＲＤ0,DFGRD1,ＮＯEＬ,NPT，ＬAYER,KＳPＴ，ＫSTEP，KINＣ）INCLUDＥ‘ABA_PＡＲＡM．INC’CHARＡCＴER＊8０CMNAMＥDIMENSIONSＴＲESS(NTENS),ＳＴＡＴEＶ（NSTＡTV)，１DDSDDＥ（NTENS，NＴENＳ)，DDＳDDＴ(NＴENＳ），ＤRPＬDＥ(NTEＮS）,２STRＡN（NTENS)，ＤSTRAN（NTＥNS)，TIME（２),PREＤEF(1)，DPRＥD（1)，３PRＯPＳ（ＮPROPＳ），COＯRDＳ（3)，DROＴ（３，3），DFGRD0（３,3)，DＦＧRD1（３，３）ｕｓeｒcｏdingｔodefineＤDＳDＤE，STＲESS，STAＴEＶ，ＳＳE,SPD，SCDand，ｉfｎecesｓaｒｙ,RPL，ＤDSDDT,ＤRＰLＤE，DRPLDT，PNEWDTRETUＲNENDUMAＴ中的應(yīng)力矩陣、應(yīng)變矩陣以及矩陣DDＳDDE、DＤSDDT、DRPＬDE等，都是直接分量存儲在前，剪切分量存儲在后.直接分量有ＮＤＩ個(gè)，剪切分量有NSHR個(gè)。各分量之間的順序根據(jù)單元自由度的不同有一些差異，所以編寫ＵMAT時(shí)要考慮到所使用單元的類別.下面對ＵMAT中用到的一些變量進(jìn)行說明：DDSDDE（NTENSNTENS)：一個(gè)NTENS×NTＥＮS的矩陣,稱作Jacoｂian矩陣是應(yīng)力的增量，是應(yīng)變的增量,DDSDDＥ(i，j）表示增量步結(jié)束時(shí)第ｊ個(gè)應(yīng)變分量的改變引起的第i個(gè)應(yīng)力分量的變化。通常Jaｃobiａｎ矩陣是一個(gè)對稱矩陣，除非在“*USERMAＴＥRIAＬ”語句中加入了“ＵNSYMM”參數(shù).STＲESS（ＮTENS）:應(yīng)力張量數(shù)組，對應(yīng)ＮDＩ個(gè)直接分量和ＮSHR個(gè)剪切分量。在增量步的開始,應(yīng)力張量矩陣中的數(shù)值通過ＵMAT和主程序之間的接口傳遞到UMAＴ中，在增量步的結(jié)束UMＡT將對應(yīng)力張量矩陣更新.對于包含剛體轉(zhuǎn)動(dòng)的有限應(yīng)變問題，一個(gè)增量步調(diào)用UＭＡT之前就已經(jīng)對應(yīng)力張量進(jìn)行了剛體轉(zhuǎn)動(dòng)，因此ＵMAT中只需處理應(yīng)力張量的共旋部分。UMAT中應(yīng)力張量的度量為柯西應(yīng)力.STATEV（NSTATＥＶ)：用于存儲與解有關(guān)的狀態(tài)變量的數(shù)組，在增量步開始時(shí)將數(shù)值傳遞到UＭＡT中，也可在子程序USDFLＤ或ＵＥXＰAN中先更新數(shù)據(jù),然后增量步開始時(shí)將更新后的資料傳遞到UMAＴ中。在增量步的結(jié)束必須更新狀態(tài)變量矩陣中的數(shù)據(jù)。和應(yīng)力張量矩陣不同的是:對于有限應(yīng)變問題,除了材料本構(gòu)行為引起的資料更新以外，與解有關(guān)的狀態(tài)變量矩陣中的任何向量或者張量都必須通過旋轉(zhuǎn)來考慮材料的剛體運(yùn)動(dòng)。狀態(tài)變量矩陣的維數(shù)通過ABAＱUS輸入文件中的關(guān)鍵詞“＊ＤＥＰVＡR”定義,關(guān)鍵詞下面數(shù)據(jù)行的數(shù)值即為狀態(tài)變量矩陣的維數(shù)。PＲＯPＳ（ＮPROPS）：材料常數(shù)數(shù)組。材料常數(shù)的個(gè)數(shù)，等于關(guān)鍵詞“*UＳEＲMATERIAL”中“COＮＳTANTS”常數(shù)設(shè)定的值.矩陣中元素的數(shù)值對應(yīng)于關(guān)鍵詞“ＵＳERMATERIＡＬ"下面的數(shù)據(jù)行。SSE,SＰＤ,SCD:分別定義每一增量步的彈性應(yīng)變能，塑性耗散和蠕變耗散。它們對計(jì)算結(jié)果沒有影響,僅僅作為能量輸出。ＳTＲAN(NTＥNS）：應(yīng)變數(shù)組。DSTRＡN(NTEＮS)：應(yīng)變增量數(shù)組。DＴIME：增量步的時(shí)間增量。NDI：直接應(yīng)力分量的個(gè)數(shù)。ＮＳHR：剪切應(yīng)力分量的個(gè)數(shù)。NTENＳ:總應(yīng)力分量的個(gè)數(shù)，NＴEＮS＝NＤＩ+ＮSHR。由于UMAT子程序在單元的積分點(diǎn)上調(diào)用，增量步開始時(shí)，主程序路徑將通過UMAＴ的接口進(jìn)入ＵMAT,單元當(dāng)前積分點(diǎn)必要變量的初始值將隨之傳遞給UMAT的相應(yīng)變量。在UMAT結(jié)束時(shí)，變量的更新值將通過接口返回主程序。ＵＭＡT的使用方法我們知道，有限元計(jì)算(增量方法)的基本問題[７］是：已知第ｎ步的結(jié)果（應(yīng)力，應(yīng)變等）,，然后給出一個(gè)應(yīng)變增量，計(jì)算,UMAT要完成這一計(jì)算，并要計(jì)算DDSDＤE（I,Ｊ)=。是應(yīng)力增量矩陣，是應(yīng)變增量矩陣，DDＳDDＥ(I,J）定義了第J個(gè)應(yīng)變分量的微小變化對第Ｉ個(gè)應(yīng)力分量帶來的變化。該矩陣只影響收斂速度，不影響計(jì)算結(jié)果的準(zhǔn)確性(當(dāng)然，不收斂自然得不到結(jié)果)。有限元計(jì)算的中心問題就是求得節(jié)點(diǎn)的位移(進(jìn)而應(yīng)變、應(yīng)力）,以使內(nèi)力和外力達(dá)到平衡：（３-1）d是節(jié)點(diǎn)位移矩陣，黑體字表示矩陣或矢量。除了小變形、線彈性問題,方程2—1是非性的,要用迭代的方法解出:（3-2）(３—3）i表示一個(gè)增量步內(nèi)的第i次迭代，ｎ表示第n個(gè)增量步.KT是切線剛度,由材料的Jacoｂｉan矩陣結(jié)合單元計(jì)算組裝而得。剛度矩陣其實(shí)就是力對位移的梯度。要想快速收斂,位移增量應(yīng)沿該梯度方向變化,也就是說,如果Jacobiaｎ矩陣不是那么準(zhǔn)確，自然KT也不怎么準(zhǔn)確,那么滿足3-１式的位移被找到的速度也就變慢，甚至發(fā)散,根本找不到。但收斂速度無論慢快,3—1式才是判斷結(jié)果準(zhǔn)確與否的唯一標(biāo)準(zhǔn)。所以Ｊａcobian矩陣不影響結(jié)果的準(zhǔn)確性,只影響收斂速度的快慢。以大變形、非性材料為例,整個(gè)計(jì)算步驟是這樣的：整個(gè)外力不是一次加上，而是一點(diǎn)點(diǎn)加上的,不然會發(fā)散得不到結(jié)果的。所以,每一個(gè)增量步開始時(shí)就是在原來的外力上加上一點(diǎn)點(diǎn),得到.根據(jù)3-２得到位移增量，此時(shí)要知道力對位移的梯度KＴ，以盡快找到滿足平衡條件的位移，由材料的Jacobian矩陣和單元結(jié)合起來組裝得到（此處使用UMAT提供的Ｊａcobiaｎ矩陣）。然后可計(jì)算應(yīng)變增量，調(diào)用UMAT，得到新的應(yīng)力,進(jìn)而得到新的內(nèi)力,所以，程序不在乎新的應(yīng)力是由增量方法得到，還是全量方法得到，而只在乎新應(yīng)力是否準(zhǔn)確.然后回到3—2，如此循環(huán),直至３-2右端為0，也即滿足3-1。這樣第n+１步就完成了,然后開始第n+2步,即外力加上一點(diǎn)點(diǎn)，按同樣的方法求解新的位移。直至整個(gè)外力全部施加并得到滿足３—１的位移。材料非線性問題彈性力學(xué)作為精確理論,從本質(zhì)上都是非線性的,早期Ｃａuchy,Greｅn，Kircｈhｏｆf和Kｅｌｖiｎ在這些方面都作出了重要貢獻(xiàn)。后來又提出超彈性(即具有彈性勢的)有限變形理論,由于理論方程的冗長而復(fù)雜，且工程應(yīng)用也沒有提出這方面要求而被擱置.20世紀(jì)4０年代以后，由于橡膠材料、高分子合成材料的迅速發(fā)展和工業(yè)領(lǐng)域的大量應(yīng)用,非線性彈性與超彈性的研究再次引起科學(xué)和工程界的重視,除了一般理論研究有了新的發(fā)展以外，工程應(yīng)用計(jì)算方法也得到長足的發(fā)展。材料的彈塑性本構(gòu)關(guān)系彈塑性材料進(jìn)入塑性的特征是當(dāng)荷載卸去后存在不可恢復(fù)的永久變形.所以，在卸載情況下，應(yīng)力應(yīng)變之間不再是唯一的對應(yīng)關(guān)系。這是區(qū)別于非線性彈性材料的基本屬性。只以加載時(shí)應(yīng)力應(yīng)變關(guān)系成非線性,還不足以判斷材料是非線性彈性還是彈塑性。但是一經(jīng)卸載就可以看出兩者的區(qū)別。非線性彈性材料沿原路徑返回，而彈塑性材料將依據(jù)不同的加載歷史卸載后產(chǎn)生不同的永久變形。對大多數(shù)材料來說，在單調(diào)加載的情況下，存在一個(gè)明顯的極限應(yīng)力，當(dāng)應(yīng)力低于時(shí),材料保持線彈性。而當(dāng)應(yīng)力達(dá)到以后，則材料開始進(jìn)入彈塑性狀態(tài)。如繼續(xù)加載，然后在卸載，材料始終保持永久的塑性變形.如果應(yīng)力達(dá)到后，應(yīng)力不再增加,而材料變形可以繼續(xù)增加，及變形處于不定的流動(dòng)狀態(tài)，則稱材料為理想彈塑性的。反之如果應(yīng)力達(dá)到后,再增加變形,應(yīng)力也必須增加，則材料是應(yīng)變硬化的，這時(shí)應(yīng)力是塑性應(yīng)變的函數(shù)，可解析為:（4—1)本構(gòu)關(guān)系反應(yīng)著應(yīng)力應(yīng)變之間的關(guān)系。對于彈性材料變形是可以恢復(fù)的；而塑性材料變形是不可以恢復(fù)的。典型的彈塑性應(yīng)變在卸載后要保持一個(gè)永久的變形。如圖3－2圖4-１塑性應(yīng)變有下列特性:（1)總應(yīng)變分為彈性和塑性兩部分,即（4-2）或者：(4－3）（2）塑性變形取決于加載路徑，而應(yīng)力應(yīng)變之間沒有一一對應(yīng)的關(guān)系。所以必須確定二則之間的本構(gòu)關(guān)系，這種本構(gòu)關(guān)系可以用偏微分方程或者增量形式來描述?？傊?，彈塑性理論主要包括以下幾個(gè)方面：（1）應(yīng)變張量的分解；（2）應(yīng)力空間的屈服條件；（3）流動(dòng)法則;（4）強(qiáng)化法則；（5）協(xié)調(diào)性條件。1：本構(gòu)模型塑性力學(xué)的應(yīng)力-應(yīng)變曲線通常有５種簡化模型［8］：(1）理想彈塑性模型，用于低碳鋼或強(qiáng)化性質(zhì)不明顯的材料。（２)線性強(qiáng)化彈塑性模型,用于有顯著強(qiáng)化性質(zhì)的材料。（3）理想剛塑性模型，用于彈性應(yīng)變比塑性應(yīng)變小得多且強(qiáng)化性質(zhì)不明顯的材料。（4）線性強(qiáng)化剛塑性模型，用于彈性應(yīng)變比塑性應(yīng)變小得多且強(qiáng)化性質(zhì)明顯的材料。（5）冪強(qiáng)化模型，為簡化計(jì)算中的解析式，可將應(yīng)力-應(yīng)變關(guān)系的解析式寫為σ=σｙ（ε/εｙ）n,式中σｙ為屈服應(yīng)力，εy為與σy相對應(yīng)的應(yīng)變,ｎ為材料常數(shù)。圖4-2２：屈服條件在復(fù)雜應(yīng)力狀態(tài)下,判斷物體屈服狀態(tài)的準(zhǔn)則稱為屈服條件[9]。屈服條件是各應(yīng)力分量組合應(yīng)滿足的條件。對于金屬材料，最常用的屈服條件為最大剪應(yīng)力屈服條件(又稱Ｔｒesca屈服條件）和彈性形變比能屈服條件(又稱VonＭises條件）.對于巖土材料則常用Treｓca屈服條件、Dｒuｃker—Prager屈服條件和Mohr—Couloｍb屈服條件.對于強(qiáng)化或軟化材料，屈服條件將隨塑性變形的增長而變化,改變后的屈服條件稱為后繼屈服條件.當(dāng)已知主應(yīng)力的大小次序時(shí)，使用Tｒesca屈服條件較為方便;若不知道主應(yīng)力的大小次序，則使用VonＭises屈服條件較為方便。對于韌性較好的材料,VｏｎMises屈服條件與試驗(yàn)數(shù)據(jù)符合較好。VonＭises屈服準(zhǔn)則具體形式是,對于各項(xiàng)同性材料，應(yīng)力偏量第二不變量等于某一定值時(shí),材料開始進(jìn)入了塑性狀態(tài)。（４—4）3:強(qiáng)化法則對理想的彈塑性材料而言，因無強(qiáng)化作用，所以，整個(gè)塑性變形過程中,屈服函數(shù)值保持一個(gè)常量,強(qiáng)化定義了屈服面在應(yīng)力空間的演化準(zhǔn)則。（４-５）其中,是強(qiáng)化參數(shù)。通常采用的強(qiáng)化法則有以下幾種:（１)各向同性強(qiáng)化此法則規(guī)定材料進(jìn)入塑性變形以后，加載曲面在各方向均勻的向外擴(kuò)張，沒有畸變。而其形狀、中心及其在應(yīng)力空間的方位均保持不變［１0］。需要指出的是：各向同性強(qiáng)化法則主要適用于單調(diào)加載情況。如果用于卸載情況,它只適合反向屈服應(yīng)力等于應(yīng)力反轉(zhuǎn)點(diǎn)的材料,而通常材料不具備這種性質(zhì)，因此在塑性力學(xué)中還發(fā)展了其它強(qiáng)化準(zhǔn)則.（2)隨動(dòng)強(qiáng)化此法則規(guī)定材料進(jìn)入塑性狀態(tài)以后，加載曲面在應(yīng)力空間作剛體移動(dòng)而沒有轉(zhuǎn)動(dòng),因此初始屈服面的形狀、大小和方向仍然保持不變。（３）混合強(qiáng)化把各向同性強(qiáng)化模型和隨動(dòng)強(qiáng)化模型加以組合，得到混合強(qiáng)化模型。它假定在塑性變形過程中，加載曲面不但作剛性平移，還同時(shí)在各個(gè)方向作均勻擴(kuò)大.在以上幾種強(qiáng)化模型中，各向同性強(qiáng)化模型應(yīng)用最為廣泛.本文也是采用該硬化法則，這一方面是由于它便于進(jìn)行數(shù)學(xué)處理;另一方面，如果在加載過程中應(yīng)力方向（或各個(gè)應(yīng)力分量的比值）變化不大，采用各向同性強(qiáng)化模型的計(jì)算結(jié)果與實(shí)際情況也比要符合.隨動(dòng)強(qiáng)化模型可以考慮材料的包興格(Bauschiｎgｅｒ)效應(yīng),在循環(huán)加載或可能出現(xiàn)反向屈服的問題中，需要采用這種模型。由于塑性變形與變形歷史有關(guān),因此反映塑性應(yīng)力-應(yīng)變關(guān)系的本構(gòu)關(guān)系用應(yīng)變增量形式給出比較方便。用應(yīng)變增量形式表示塑性本構(gòu)關(guān)系的理論稱為塑性增量理論。增量理論的本構(gòu)關(guān)系在理論上是合理的，但應(yīng)用比較麻煩，因?yàn)橐e分整個(gè)變形路徑才能得到最后結(jié)果.因此，又發(fā)展出塑性全量理論，即采用全量應(yīng)力和全量應(yīng)變表示塑性本構(gòu)關(guān)系的理論。在比例變形的條件下，可通過積分增量理論的本構(gòu)關(guān)系獲得全量理論的本構(gòu)關(guān)系。當(dāng)偏離比例變形條件不多時(shí),全量理論的計(jì)算結(jié)果和實(shí)險(xiǎn)結(jié)果比較接近。本文的程序都是基于增量理論。非線性有限元算法理論對于非線性問題，在有限元求解該問題時(shí)，對一個(gè)自由度總可以表達(dá)成,式中,為基本未矢量。如果是線性問題,與無關(guān)，而是一次項(xiàng),顯然這是一個(gè)線性方程。如果與相關(guān),則方程的出現(xiàn)非一次項(xiàng)，變成非線性問題，在實(shí)際工程中，特別是塑性成型問題,材料的幾何方程，本構(gòu)方程以及邊界條件往往是非線性的也體現(xiàn)在中出現(xiàn)了,所以變?yōu)榱朔蔷€性問題，要得到最基本的未知量，就必須求解非線性方程組１：直接迭代法又稱常剛度法［1１］，這是種最簡.單的求解方法，在每次求解前，利用上次的解來求出這一次的值，然后利用和的倒數(shù)的乘積求出的當(dāng)前值（4—6)表達(dá)為迭代形式（4-７)上式可以看出，這種方法首先需要有一個(gè)初始的值，以便開始迭代。另外，每一次求解都需要對求倒數(shù)，如果求解方程組，就是對剛度矩陣求逆，這種方法在求解中控制兩次求解之差,當(dāng)其值很小時(shí)，就認(rèn)為接近真實(shí)值了,迭代結(jié)束圖4－3２：Ｎeｗton-Raｐhson方法Ｎewｔon-Rａphson方法的算法與常剛度法不同［12],如果得近似表達(dá)式是不成立的，存在著殘余值,即,此式也可以作為近似值與真實(shí)值的差值量度，實(shí)際上在具體計(jì)算時(shí)，也可以控制其值，當(dāng)極小時(shí)，就認(rèn)為接近真實(shí)值了，當(dāng)?shù)诖蔚闹凳钦鎸?shí)解，則可以按照Tａyloｒ級數(shù)展開得到圖4－43：切線剛度法在復(fù)雜非線性問題求解中,剛度與的大小是有一定關(guān)系的，在用增量法來求解這種問題時(shí)，就等于結(jié)構(gòu)任一點(diǎn)處力與位移的曲線的局部梯度，稱為切線剛度［１３],剛度矩陣的倒數(shù)很難用自變量顯示表達(dá),通過增量方式求解,在每一步荷載增量范圍內(nèi)把問題線性化，求解方法與Nｅwｔｏｎ-Rａpｈsoｎ方法相同總結(jié)以上可以得到:以上幾種算法中,通過比較，不難發(fā)現(xiàn)，直接迭代法采用了固定的剛度,適合解決非線性程度不高的本構(gòu)關(guān)系,而切線剛度法采用了變化的剛度，在每一步上都做了實(shí)時(shí)的修正，對非線性程度較高本構(gòu)關(guān)系任然有效，在效率和迭代精度方面,切線剛度法采用的修正更符合非線性材料的應(yīng)力應(yīng)變關(guān)系,具有較大的優(yōu)勢，這也是本文采用切線剛度法計(jì)算的原因,當(dāng)然，非線性有限元算法還有很多,切線剛度法也不見得就是最好的能解決所有問題的算法，但是它是在程序開發(fā)難度不高和精度方面較高的條件下相對來說最好的本文采用的本構(gòu)關(guān)系是同性硬化彈塑性模型［１4]，采用Mｉｓes屈服準(zhǔn)則，下面將根據(jù)J2理論[１5］,分別推導(dǎo)常剛度法和切線剛度法計(jì)算該問題的的算法公式增量理論常剛度法公式推導(dǎo)由應(yīng)力應(yīng)變關(guān)系得:(4-8）（4-9）其中是彈性矩陣，它的表達(dá)式為是塑性應(yīng)變增量，它的表達(dá)式為其中為等效塑性應(yīng)變增量，它的表達(dá)式為為切線模量對于３維空間問題，流動(dòng)方向：等效應(yīng)力應(yīng)力偏量增量理論切線剛度法公式推導(dǎo)本文采用的是一種切線剛度法，其應(yīng)力應(yīng)變關(guān)系為(4-10）(4－1１)３－14兩端左右同乘，得到:對于強(qiáng)化材料，Ｍises準(zhǔn)則代入可得其中，為等效塑性應(yīng)變增量,它的表達(dá)式為由于：由流動(dòng)法則可知應(yīng)用上式得所以得到(4－12)這就是切線剛度法的矩陣表達(dá)式其中為彈塑性矩陣,在ＡBAQUＳ里面稱為雅可比矩陣,它的表達(dá)式為其中是彈性矩陣,它的表達(dá)式為為塑性矩陣，它的表達(dá)式為對于３維空間問題流動(dòng)方向所以最后推導(dǎo)可得，彈塑性矩陣的表達(dá)式其中為切線模量，對本構(gòu)關(guān)系求導(dǎo)得到總結(jié)推導(dǎo)過程：上面的推導(dǎo)過程看似復(fù)雜，其實(shí)核心的問題只有一個(gè),即兩者對于塑性階段應(yīng)力更新的算法不同。常剛度法采用的是：而切線剛度法采用的是：把握好了兩者的本質(zhì)上的區(qū)別，對于兩者的算法設(shè)計(jì)和程序開發(fā)問題便迎刃而解UＭAT程序設(shè)計(jì)和編碼本章將嚴(yán)格按照前一章推導(dǎo)的公式展開程序設(shè)計(jì)和編碼，為了便于編程,本文將本構(gòu)關(guān)系做了抽象化處理，即將其描述成一個(gè)含參數(shù)的表達(dá)式，改變參數(shù)即可應(yīng)用于不同的模型，這樣做的好處是能保證程序的復(fù)用性,這也是本文反復(fù)強(qiáng)調(diào)的使用UMＡT的原則。本構(gòu)關(guān)系描述本文采用各向同性硬化彈塑性材料,材料參數(shù)如下:(5－1)彈性部分:,彈性模量E=2００000Ｍpa,泊松比Mu=0。3圖５-1彈性部分本構(gòu)關(guān)系塑性部分:，為了研究方便,取A=700,B=０.５，C＝400圖5-２塑性部分本構(gòu)關(guān)系將2個(gè)曲線統(tǒng)一到同一個(gè)坐標(biāo)系(為方便顯示，x軸標(biāo)注時(shí)擴(kuò)大了100０倍)圖5－3本構(gòu)關(guān)系由此可以求出兩條線的交點(diǎn)即初始屈服點(diǎn)的應(yīng)力Yield0＝400Mpａ（注意兩條曲線相差了一個(gè)屈服應(yīng)變,因?yàn)閮烧咂鋵?shí)不是一個(gè)坐標(biāo)系）綜上定義的材料常數(shù)見表５—1：表５—1彈性模量Ｅ200000泊松比Mu０。3屈服應(yīng)力Yiｅld0400A７00Ｂ0．5C400注意：上面Ａ,B，C的取值只是為了便于理解和分析本文的材料模型,為了保證程序的通用性，本文的參數(shù)在程序中的A，B，C一律用變量表示。常剛度法程序設(shè)計(jì)算法設(shè)計(jì)1：定義程序需要用到的常數(shù)和變量2:讀取AＢAＱUＳ定義的材料常數(shù)和狀態(tài)變量（這里只定義了一個(gè)狀態(tài)變量），材料常數(shù)為,彈性模量E，泊松比Ｍu,屈服應(yīng)力Ｙｉeld0，參數(shù)A,B，C，并且計(jì)算出剪切模量G，狀態(tài)變量為等效塑性應(yīng)變ＥQPLAS3：讀取應(yīng)力分量，計(jì)算平均應(yīng)力,應(yīng)力偏量以及Mises等效應(yīng)力平均應(yīng)力:應(yīng)力偏量：Misｅs等效應(yīng)力:４:根據(jù)３計(jì)算的Mｉses等效應(yīng)力和2讀取的屈服應(yīng)力Ｙield0比較,如果Ｍｉｓeｓ等效應(yīng)力小于屈服應(yīng)力，表明此時(shí)材料未屈服,那么轉(zhuǎn)到５，否則轉(zhuǎn)到６５:雅可比矩陣，初始化為0,計(jì)算彈性矩陣,按照彈性理論更新應(yīng)力6:雅可比矩陣,初始化為0：計(jì)算切線模量H’,注意到當(dāng)?shù)刃苄詰?yīng)時(shí)對應(yīng)于本構(gòu)關(guān)系的屈服點(diǎn),此時(shí)的Ｈ不能通過上式計(jì)算，可以取此時(shí)的H為彈性模量：計(jì)算等效塑性應(yīng)變增量并更新EQＰLAＳ=DEQPLＡS+ＥＱＰＬＡＳ更新狀態(tài)變量3)計(jì)算流動(dòng)方向?計(jì)算塑性應(yīng)變增量５）更新應(yīng)力算法流程圖圖5-4常剛度法算法流程圖常剛度法程序編碼根據(jù)算法流程，用FOＲＴRAN77固定格式編制了常剛度法的計(jì)算程序如下:SUBROUTINEＵMAT(ＳTRＥSS,STATEV,DDSDDE，SSＥ，SＰD,SＣD，RPＬ，DDSＤDT， 1DRPLDＥ,DRPLDT,STＲＡＮ,DSＴRAN,TIMＥ，ＤＴIME,ＴEMP，DTＥMＰ,PＲEDＥF，DＰRED，2CＭNAME，ＮDI，ＮＳHR,NTＥＮＳ，NＳTAＴV,PＲOPS,ＮPROＰS,ＣＯOＲＤＳ,ＤROＴ，3PNEWDＴ,ＣELENT，DFGRD0,DFＧＲD１，NOEL，ＮPT，LAYER，ＫＳPＴ，KＳTＥP,ＫINＣ）?INCＬUＤE’ABＡ_PARAM.INC＇CHＡRＡCTER＊8CMNAMEＤＩMENSIONＳTRＥSS(NTENS),ＳTＡTEV(NSＴAＴV),DＤＳＤDE（ＮＴＥNＳ，ＮTＥＮS)，１DDSDDT(ＮTENS)，DRＰLDE（ＮTENS),ＳTRAN(NＴENS)，DＳTＲAN(NTENS），2ＰRＥDEＦ(1)，DPRED（1），PＲOPS（ＮＰＲOPS）,COＯRDS（3），ＤROT(３,3)，3DＦＧRＤ０（3,3)，DFGRD1（3，3)，ＤPLAS（6）C———－-—-－————－--—-－－—— 常剛度法CCＣCCＣCCCＣCＣＣCCCＣＣＣCＣCCＣＣC1：定義變量ｃ定義常數(shù) PＡRＡＭETEＲ（ZERＯ＝0．D０，OＮE=1。D０，ＴWO=２．D０,TＨREE=3。D0）C定義材料常數(shù) DOUBLEPRECＩSIONＥ,Mｕ，Yｉeld0，A,Ｂ,CC定義中間變量 ＤＯUBLEPRECＩSIＯＮＨ,Ｗ，CEＧＭA_EQ,CEGＭA＿X，ＣEＧMA_Ｙ，CEGＭA_Z，TAO_XＹ，?1TAO_ＹZ，TAO_ＺX,CEＧＭＡ＿CＰ，CEGMA_ＳX,CEＧＭA＿SＹ，ＣＥGMA＿SZ，TＡO＿SＸY，TAO_ＳYZ2,ＴAO_SZXC定義狀態(tài)變量 DOUBLEＰＲECISIONEQPLAＳC 定義更新變量ＤOＵBＬEPRECＩSIＯNDEQPLASCCCCCCCCCCCCCCＣCCCCCＣCCＣCC２:讀取變量C 彈性模量E=PＲOPS(1)C 泊松比 Mｕ＝PＲOPS（2)C剪切模量?G＝E/ＴWＯ/(ＯNＥ+Mｕ）C 屈服應(yīng)力?Ｙiｅlｄ0＝PROPS（３) C參數(shù)?Ａ=PROPS(４) B=PROPＳ（5） Ｃ=PRＯＰS（6)Ｃ 讀取狀態(tài)參數(shù) EＱPLAＳ=SＴＡTEV（1） Ｙield=A*（ＥQＰＬAＳ+Yielｄ0/E）＊*B+CＣCCCCCＣCCCCＣＣCCCCCCCCCCCCＣ3：讀取應(yīng)力分量，計(jì)算平均應(yīng)力，應(yīng)力偏量以及Miｓeｓ等效應(yīng)力Ｃ6個(gè)應(yīng)力分量 CEGＭＡ_Ｘ=SＴRESS（1） CEGMＡ_Y=STRESS(2）CEGMA_Z＝STRESS(3)?TAO_ＸY＝ＳTRESS（4) TAO_ＺX＝STRESS（5)?ＴＡO_YZ=STＲESS（6）C平均應(yīng)力?CEＧMA＿ＣＰ=(ＣEGMＡ_X＋CＥＧMA＿Y＋ＣEGMA_Ｚ）/TＨREEC６個(gè)應(yīng)力偏量?CEGMA_SX＝CEGMＡ_Ｘ－CEGMA_CP CＥＧMＡ＿ＳY=CEGMA_Y-CEＧMＡ_ＣP CＥＧMＡ_ＳZ=CEGMA_Ｚ—ＣＥGMA＿CP ＴAO_SXY=TAO＿XY TAO_ＳZX＝ＴAＯ_ZX?TAO_SＹＺ=TAO＿YZcMises等效應(yīng)力?CEGMA＿EＱ=SQＲＴ（THRＥE/ＴWO＊（CEGＭＡ_SX＊＊2+CEGMＡ_ＳY*＊2＋CEGMA_SZ＊*2+ 1TWO＊(TAO_SXY*＊2＋ＴAＯ_SYZ*＊2+ＴAO_SZX＊＊2）))?C雅可比矩陣，初始化默認(rèn)為 DＯk1=1，ＮTＥNS DOｋ２=1,NＴENS DDSＤＤＥ（Ｋ1,K2)=ZERO ENDDO ＥNＤDOc計(jì)算彈性矩陣DOＫ1=1,NDI DＯK2=1，ＮDI DDＳDＤE(Ｋ１，k2)＝E＊（ONE—Mu)／(ONE＋Mu)/（ONＥ—TＷO＊Mｕ)*(Mu/（ONＥ—Mu)）?EＮDDＯ DDSDDE（K１，K１）＝E＊（ONE－Mu）/（ONE＋Ｍu）/（OＮE-TＷO*Mu） ENDDO DOＫ1＝NＤI+1，NTＥNＳ?DDSＤDＥ（K1，K1)=E＊(ONE－Ｍu)／（OＮE+Ｍu）/（ＯＮE—TＷO*Mu）*(ONＥ－TWO＊Mu）/TＷＯ/1（OＮE-Ｍu) ENDDＯ ?CＣCCCCCＣCCCCCＣＣCCCCCCCCＣCC4：根據(jù)計(jì)算的Misｅs等效應(yīng)力和讀取的屈服應(yīng)力Yield0比較,CCCCCＣCCCCＣCCCCCCCCCＣCCCCC如果Misｅs等效應(yīng)力小于屈服應(yīng)力,表明此時(shí)材料未屈服,那么轉(zhuǎn)到,否則轉(zhuǎn)到?ＩF(ＣEGＭA_ＥQ．ＬT.Yiｅlｄ0）THENCＣCCCＣCCCCCCＣCCCCCCCCCCCCC5：沒有屈服，按照彈性理論計(jì)算Ｃ按照彈性理論更新應(yīng)力 DOK1＝1，ＮTEＮSDOK2＝1，NＴENSSTREＳS（Ｋ２）＝STREＳS(Ｋ2）＋DDSDＤE(K２，K1)*ＤＳTRAN（K1)ENＤDO?ENDDO ELSEＣCCCＣＣCCＣＣCCCCCＣCCＣCCCCCＣC6：屈服發(fā)生,按照塑性理論計(jì)算C切線模量H,根據(jù)本構(gòu)關(guān)系求導(dǎo)IF（EＱPLAS。ＥＱ。0）?THEN Ｈ=E EＬSE H=Ａ＊B＊EQPLAS*＊(B-ONE） ＥＮDIFC求WW＝9．D0*G/TＷO/ＣEＧMA_EQ*＊2/（H+THREE＊G)Ｃ等效塑性應(yīng)變增量?DEQPLAS=(6．Ｄ0＊G*CEGＭA_EQ/（TWO＊ＣEＧMA＿EQ*＊2＊H+9.Ｄ0*Ｇ＊（CEＧMＡ_SＸ１**２+CEGＭＡ_SY**2+CEGMA_ＳZ＊＊2+TＷＯ*TＡＯ_XY*＊2＋TＷＯ＊TAO_ZX*＊2+ＴＷＯ*TAO_ＹZ2＊*2)）)＊(CEGMA_SX＊DSTRＡN(1）+ＣEGMA_SY＊DＳＴRＡN(2)+CEＧMA_SＺ* 3DＳTＲＡN(3)+TAO_ＸY＊DSTRAＮ（4）+TＡO＿ＺＸ*DSＴRAN(5)+TAＯ＿YZ*ＤSTRAN(6)）C?wｒｉte（１０,＊)”DEＱPＬＡS",ＤＥQＰLAＳ,EＱPLＡS,ＨC更新狀態(tài)變量STATＥV(1）=EQPLAS+ＤＥQＰLＡSＣ計(jì)算塑性應(yīng)變增量? DPLAS（1）=DEQＰＬAS*THRＥE＊ＣEＧMA_ＳX/TWO/ＣEGMA_EQ ＤPLAS（２)＝ＤＥQPＬAＳ*THREE＊ＣEGMA_SY/ＴWO/CEGMA_ＥQ ＤPLＡS（3)=DEQＰLAS＊THREE＊CＥGMＡ＿ＳＺ/ＴＷO/CEGMＡ＿EQ DPLAS(4）＝DEQPLAS＊THREE＊TＡＯ_XＹ/CEＧMA＿EQ DPＬAS(5）＝ＤＥQＰＬAＳ＊TＨＲEE＊TＡＯ_ZX/ＣEGＭＡ_EQ?DPLAS（６）＝DＥＱＰLＡＳ＊THRＥE＊TAO_ZY／CEGMA_EQ?C按照塑性理論更新應(yīng)力 DOＫ1=1,ＮＴEＮSDＯK2=1，ＮTENSSTRESS(Ｋ2)＝STＲEＳS（K2)＋DＤＳＤDE(K2,K1）*(DSTRＡN（K1）—DPＬAS（K1)）ＥNDＤO EＮDＤO?ENDIF?? ??ＣCCCCCＣCCCCCCCCCCCCCＣCCＣCC計(jì)算完成 REＴＵRN?END切線剛度法程序設(shè)計(jì)算法設(shè)計(jì)1：定義程序需要用到的常數(shù)和變量２:讀?。罛AQUS定義的材料常數(shù)和狀態(tài)變量(這里只定義了一個(gè)狀態(tài)變量）,材料常數(shù)為，彈性模量E，泊松比Mu，屈服應(yīng)力Yielｄ0,參數(shù)Ａ,B，C,并且計(jì)算出剪切模量Ｇ，狀態(tài)變量為等效塑性應(yīng)變EＱPＬＡS３:讀取應(yīng)力分量，計(jì)算平均應(yīng)力,應(yīng)力偏量以及Misｅｓ等效應(yīng)力平均應(yīng)力：應(yīng)力偏量：Mises等效應(yīng)力：4：根據(jù)3計(jì)算的Mises等效應(yīng)力和2讀取的屈服應(yīng)力Ｙielｄ0比較，如果Misｅs等效應(yīng)力小于屈服應(yīng)力，表明此時(shí)材料未屈服,那么轉(zhuǎn)到５,否則轉(zhuǎn)到６5:雅可比矩陣,初始化為0，然后計(jì)算彈性矩陣,按照彈性理論更新應(yīng)力6：雅可比矩陣,初始化為0計(jì)算切線模量H注意到當(dāng)?shù)刃苄詰?yīng)時(shí)對應(yīng)于本構(gòu)關(guān)系的屈服點(diǎn)，此時(shí)的Ｈ不能通過上式計(jì)算,可以取此時(shí)的H為彈性模量計(jì)算w根據(jù)前一章推導(dǎo)的公式計(jì)算等效塑性應(yīng)變增量DEＱＰLAS，并更新狀態(tài)變量根據(jù)前一章推導(dǎo)的公式：帶入后可得等效塑性應(yīng)變增量DEQPLAS然后ＥQＰＬＡS=ＤEQＰLＡＳ+EＱPLAS更新狀態(tài)變量計(jì)算雅可比矩陣首先初始化默認(rèn)為０，然后用下式計(jì)算雅可比矩陣注意:aｂaquｓ的剪應(yīng)力方向跟彈性力學(xué)規(guī)定的方向不一致，所以上式的最后兩行的應(yīng)交換更新應(yīng)力算法流程圖圖5－5切線剛度法算法流程圖切線剛度法程序編碼根據(jù)算法流程,用ＦＯRTRＡＮ77固定格式編制了切線剛度法的計(jì)算程序如下:SＵBRＯＵTINEUMAT(ＳＴＲＥSS,STATEV，ＤＤSＤDE，SSE，SPD，SCD,RPL,DDＳDDT,1DRPLＤE，ＤＲPLDＴ，ＳTRAN,DＳTRAＮ,TIME，ＤTIME，ＴEMP，ＤTＥＭP,PREＤEF，DPRED,２CＭNAME,ＮDI，NＳHR,NTEＮS，NSＴATＶ,PＲOPS，NPROPS，COORDS，DＲＯＴ，３PNＥWDT，ＣＥLENT，DFGＲＤ０，DFGRD１,NOEＬ,ＮPT,ＬＡYEＲ,ＫSPT,ＫSTEP，KＩNC）INCＬUDE＇ABＡ_PＡRAM．INC＇CHAＲACTER*8ＣMNAMEDIMENSIＯNSＴRＥSS（NTENS），STAＴEV（ＮSＴATV），DDＳDDE（NTEＮS，NTＥNＳ），1DDＳDDT(NTENＳ），DRＰLDE（NTENS)，ＳTＲAＮ(ＮTEＮＳ),DSTＲAN（ＮＴＥNS)，2PRＥＤEF(1），DPRED（1），ＰROＰS（NPROＰS),ＣＯORＤS（3），DＲOT（3,3)，3DＦGRD0(３，3）,ＤFGRD１（３，3）CCCCCCCCCCCＣCCＣＣCCCCＣCＣCＣC1:定義變量ｃ定義常數(shù)PAＲＡMＥＴEＲ(ＺERＯ＝0.D0,ONＥ=１.D0,TWＯ=2。D0,ＴHREE=３.D0)C定義材料常數(shù)ＤOUBＬEPRECISIONE，Mｕ,Yieｌd０,A,B，CＣ定義中間變量DOＵＢLＥPＲECＩＳIＯNH,Ｗ，CＥＧMA_EQ，CEＧMA＿X，CEGMA_Y,CEGＭＡ_Z，TAＯ_XY，1TAO＿YZ，TAＯ_ＺX,ＣＥGMA_CP，ＣＥGMA_ＳX,CＥGＭA_SY，CEGMＡ＿SZ，TＡO_SXY，ＴAO_SYＺ２,ＴAＯ_SZXC定義狀態(tài)變量DOUＢLEPＲECISＩＯNEQPLASＣ?定義更新變量DOUBLEＰRＥCISIONＤＥQPLASCCCCCCCCCCCＣＣCCCCCCＣCCCCCＣ2：讀取變量C 彈性模量E=PROＰＳ(1）C 泊松比Mu=PROPＳ（2)Ｃ剪切模量G=E/TWO/（ＯNE+Ｍu）Ｃ 屈服應(yīng)力 Yiｅld0=ＰＲＯＰＳ（3)?C參數(shù) Ａ=PRＯPS(４） ?B=PROPS（5）? Ｃ=ＰROPS(6）Ｃ?讀取狀態(tài)參數(shù)??EＱPLAS=SＴATEＶ(1） ?Yield=A*(ＥQＰLＡＳ+Yｉeld0/Ｅ）＊*B＋ＣＣＣCCCＣCCCＣCCＣCCCCＣCCＣCCCＣ