ABAQUS的約定

2021-01-16 技術鄰CAE學院

來源:CAE愛聯盟

嗯,沒錯,ABAQUS和我們有個約定,約定了什麼呢?比如自由度、單位制、應力應變等符號的含義。其實對應於ABAQUS幫助文檔Abaqus Analysis User's Guide 1.2.2 Conventions,如果沒有遵守這些約定,冥冥之中我們就種下錯誤的小種子!約定的內容總共包含:自由度、坐標系統、單位、時間尺度、曲面方向、應力與應變、旋轉、自由度。

除了軸對稱單元、流體單元和電磁單元以外的其他單元自由度進行如下的約定:

其中,x、y、z默認情況下是分別與系統的整體坐標系X、Y、Z相一致的,但如果使用*Transform對結點進行局部坐標系轉化的話,它們將與局部坐標系中的相關坐標軸一致。Abaqus/Standard中對於殼或梁單元最多可以20個溫度自由度(自由度從11到30)。



對軸對稱單元的平動與旋轉自由度如下規定:

用*transform進行結點坐標系轉換的自由度改變同上。



上述所列自由度並不是同時都能用在某一單元結點上的,不同的分析、不同的單元將擁有適合其分析的自由度,而其他自由度在某些情況下是失效的。


同一般規定一樣,ABAQUS採用的基本坐標系(系統整體坐標系)是直角坐標系,方向遵循右手法則。為便於各種分析,用戶可以自行定義局部坐標系以便於建模、定義材料、定義載荷以及變量輸出等。


建模中,常在node/ngen中加入*system;


材料中,常以*orientation進行定義局部坐標系(尤其對於殼、梁單元);


載荷中,用*transform可以定義局部坐標系下的載荷;


輸出中,先前定義的*transform(用於結點變量)與*orientation可以發揮作用。


實際上,ABAQUS沒有單位的概念,它僅是通過有限元方法對矩陣進行數學運算得到結果,理論上沒有什麼物理意義,但各種變量從人為地角度賦予物理意義以及物理定理的數學表示,從而使得ABAQUS求解出「有意義的」的結果來,由此看來,結果是否有效,人對各種數據變量的主觀把握是重要的。單位的一致性可以保證結果運算不會產生與之相關的問題。一般,ABAQUS建議用一套認可的單位制進行單位定義,比方說,國際單位制。當然,如果你願意通過一系列的轉化(轉化因子可能複雜)的話,可以不遵守單位一致性的約定。


ABAQUS中,旋轉自由度 (4、5、6) 以弧度來表示,而其他角度相關的一般都以真實角度表示(如相角、*ncopy,shift中的旋轉角度),其實,便於好記的說法就是:與長度或三角運算相關的用弧度表示,與旋轉相關的用角度表示。


國際單位制是我們最常用的一套符合一致性要求的單位系統 (SI)。總共有七個基本單位,這裡主要用到其中的五個:長度 (m),質量 (kg)、時間 (s)、溫度 (K)、電流 (A),其他相關單位均在此單位基礎上組合或推導出來。比如力,其單位為牛頓 (N=kg·m/s²);焦耳 (J=N·m=kg·m²/s²)。


然而,有時國際標準單位在分析中使用並不方便。比如說,結構分析中常用到的楊氏模量,其單位常用為MPa(=N/mm²),此時,為單位一致,與之相對應的其他基本單位應為(噸、毫米、秒)。


對我們而言,總是不太習慣使用美式或英式單位制,這是因為其單位命名規則不像國際單位制中的表示的那樣清楚。比方說,1磅力等於1磅質量 (lbm) 乘以重力加速度。如果以磅力、英尺 (ft) 和秒作為基本單位的話,則質量lbm 可以表示為:ft/sec²。但在一般情況下,密度的單位為lbm/in³,所以若要以上述基本單位表示的話,就必須轉化為:sec²/ft4。是不是感覺非常不方便?


而且,從手冊中也容易讓我們在lb 到底代表的是lbm 還是lbf 而感到模糊。我們必須通過查看其由哪些單位組成推導出才能確認是什麼。


另外,難理解的還有兩個單位,一個是斯(勒格)-slug(=lbf sec²/ft),其被定義為在受到一磅重的力作用時產生每秒鐘一英尺的加速度的質量單位。另一個是磅達-poundal,被定義為使質量1磅的質點發生1尺/秒加速度的力。以下兩個有用的轉化為slug=g·lbm 和lbf=g·poundals,這裡g 為重力加速度值。


ABAQUS共有兩種時間計法,一種是分析步時間 (step time),另一種是總分析時間 (total time)。除了線性擾動分析(它不考慮時間,頻域分析),分析步時間是從每一分析步開始計算,而總分析時間則是從第一個分析步開始計算的所有分析步的時間積累(包括*restart步)。


完全的空間曲面定義需要有局部面方向的定義,這樣在諸如基於單元接觸曲面的切向滑移方向或是殼單元的應力應變方向均能得到一致的表示。為此,ABAQUS作如下規定:






對於幾何線性分析而言,以默認的材料方向(初始參考中定義)就可以將應力應變表示出來。


對於幾何非線性分析,在ABAQUS/Standard中的小應變殼單元 (S4R5, S8R, S8R5,S8RT, S9R5, STRI3, 和STRI65) 使用總體拉格朗日應變算法,應力應變可以相對於參考構型的材料方向絕定。墊片單元是小應變小位移單元,默認情況下其應力應變值也是以初始參考構型定義的行為方向輸出。


對於有限膜應變單元(所有的膜單元以及S3/S3R, S4, S4R, SAX, 和 SAXA單元)和在ABAQUS/Explicit中的小應變單元,其材料方向是隨著曲面的平均剛性旋轉運動,而變形成當前構型的材料方向。此時,這些單元的應力應變則是根據當前的參考構型中的材料方向給出的(更詳細地說明可以參考ABAQUS相關手冊)。用戶可以決定與*section print和*section file相關的局部坐標系統是固定不動還是隨著曲面的平均剛性運動而旋轉。


在定義材料特性時,在ABAQUS中應力應變的各個分量定義如下順序:

比如,一個完全各向異性、線彈性矩陣為:

1-、2-和3-方向的確定依靠所選擇的單元類型,對於實體單元就是整體坐標系的空間方向。對於殼單元和膜單元,1-、2-方向默認上是殼或膜曲面的局部方向,可以通過*orientation對三個方向進行適合自己分析的修改。


對於實體單元的幾何非線性分析,默認的方向並不會隨著材料的旋轉而變化。然而,如果通過*orientation定義的方向則會旋轉。


ABAQUS/Explicit在分析中,在內部以不同的順序完成對各應力的存儲。對於幾何非線性分析,無論是否使用*orientation,這些被存儲的變量總是隨著材料的方向旋轉的,這一點在子程序VUMAT被使用時尤顯重要。


在連續體單元中定義各向異性材料行為時,用*orientation所定義的材料方向是很重要且必要的。


在分析中,如果應力值一直為零,在矩陣存儲時將被忽略。例如,對於平面應力分析,ABAQUS僅存儲兩個面內的正應力和一個面內的剪應力。


ABAQUS輸出工程剪應變:

在ABAQUS中,應力使用的是柯西應力或真實應力值,即每單位面積(變形後的面積)上的力大小。詳細地請看手冊"Stress measures" Section 1.5.2 of theABAQUS Theory Manual。


對於幾何非線性分析,存在多種不同的應變計量方法。不像真應力,它沒有很清楚的真應變的計算方法。對同一種物理變形,在大應變分析中不同的應變測量方法會給出不同的應變值,當然其值所能反應的實際性也就不同。如何選擇最好的應變測量方法,這依靠分析類型、材料行為以及在一定程度上也依靠個人喜好!詳見"Strain measures" Section 1.4.2 of the ABAQUS Theory Manual。


默認情況下,在ABAQUS/Standard中應變以「綜合」總應變 (E) 輸出,對於大應變的殼、膜和實體單元,還有兩種其他的總應變計法可以輸出:自然應變 (LE) 和公稱應變 (NE)。


在ABAQUS/Explicit中,自然應變 (LE) 是默認的應變輸出,也可以要求公稱應變輸出。而「綜合」應變在ABAQUS/Explicit中是無法得到的。


默認情況下,ABAQUS/Standard輸出到.dat 文件和.fil 文件的應變為綜合應變,這是對於所有可以在材料體下將應變率數值積分獲得有限應變的單元都適用的。

其中,ε[n+1]和εn 分別是分析中第n+1次和第n 次的總應變,ΔR 為旋轉張量的增量,Δε 為從第n 次增量到第n+1次增量的總應變增量。對於使用正轉坐標系統(右手法則)的單元(使用*orientation有限應變殼單元、膜單元和實體單元)上式可以簡化為:

應變增量可以通過對變形率D 在整個時間增量上積分得到:

這種應變計法對於彈塑性(粘彈性材料)或彈性-蠕變材料都是適合的,這是因為塑性應變和蠕變應變都是通過相同的積分方式得到的。在這樣的材料裡,彈性應變是很小的(因為屈服應力相對於彈模來說是較小的),此時總(綜合)應變可以直接地與塑性應變和蠕變應變相對比。


如果應變的主方向關於材料方向的變化而旋轉變化的,那麼最後所得到的應變是不能和總變形相聯繫的。此時,無論採用的是何種坐標系統。如果主應變方向保持固定,那麼應變就是變形率的積分:

這與稍後將要討論的自然應變等效。


在ABAQUS/Standard中,對於小應變殼和梁單元,默認的應變計法E 為格林應變:

此時,F 為變形適量梯度,而I 為特徵張量。這種應變計法適合於在小應變、大旋轉分析中的單元使用。分量εG 代表沿原定義方向的應變,不能在有彈塑性或超彈性材料行為的有限應變分析中使用小應變殼和梁單元,因為可能會導致不正確的分析結果!


公稱應變NE 被定義為:

其中,V=√F·FT 為左拉伸張量,λi 為主拉伸,ni  為在當前材料參數下的主拉伸方向。因此,公稱應變也就是在主材料參數方向下長度變化對原長度的比值,從而直觀地解釋了變形。


自然應變LE 被定義為:

式中,變量同上。


在ABAQUS中許多結構模型都是根據應力不變量進行闡述的,不變量的概念是指:數值不會隨著坐標系的轉換而變化,這些應力不變量為:


等效壓應力:

Mises等效應力:

第三偏應力不變量:

這中間,S 為偏應力,被定義成:


空間中,對於有限旋轉作如下規定:定義Φx、Φy、Φz 為關於整體X、Y和Z軸的旋轉量(也就是指在一個節點的4、5、6自由度)。


定義,px=Φx/Φ,py=Φy/Φ,pz=Φz/Φ,方向P 為旋轉同方向,其中Φ 為右手法則得出的旋轉角度(弧度表示)。值不是唯一確定的,在大旋轉問題中,如果總旋轉量超過2π,可以加或減任意倍的2π(以得到在-2π到2π內的值),而這可能會導致旋轉分量的不連續輸出。


這個規則提供了在大多數分析中對於運動邊條和彎矩的簡便輸入以及輸出的簡單解釋。由ABAQUS產生的對旋轉的輸出,代表的是關於一固定軸從原參考構型到當前構型的旋轉。此輸出不同於對結點的旋轉歷史輸出,而且,這個規則成為了小旋轉分析甚至是施加到初始有限旋轉上的小旋轉情況(這種情況可在關於一個預變形狀態下小振動研究考慮)最通常的規則。


因為有限旋轉的方向性,它不是累積的,它們被指定的方式與其它種邊條被指定的方式有所不同。在整個分析步上,指定的旋轉增量必須是從分析步開始的構型到該步末所需要的構型旋轉結點所需要的旋轉量。對於總旋轉量,在分析中它將結點旋轉到最終的位置,但如果它是在其他的幾個初始參考構型中施加的,那麼在一個分析步中就使結點旋轉到總量來說是不夠的(也可以說是不正確的)。如果對於某結點施加的一個旋轉增量 (ΔΦx , ΔΦy , ΔΦz )是在該分析步開始(也就是上一分析步結束時)時使該結點從其旋轉邊界條件 (Φx , Φy , Φz )t 旋轉並最終在該步結束時到其最終位置的話,那麼,在該步末的模型此結點的邊界條件諸如旋轉向量就應該是 (Φxt+ΔΦx , Φyt+ΔΦy , Φzt+ΔΦz )。如果旋轉向量的方向保持不變,那麼定義旋轉邊條和總旋轉向量的方法是相同的。









相關焦點

  • ABAQUS與你我的約定
    嗯,沒錯,ABAQUS和我們有個約定,約定了什麼呢?比如自由度、單位制、應力應變等符號的含義。
  • 詳述使用ABAQUS如何計算動剛度--實例解析
    原標題:詳述使用ABAQUS如何計算動剛度--實例解析 下面實例講述在abaqus中計算動剛度的詳細步驟: 如下圖所示,一根細長立方體,一端完全固定,在另一端的一個節點施加單位簡諧激振力。
  • 情人節歌曲:約定歌詞
    >情人節歌曲:約定歌詞2015-07-28 10:02:26出處:其他作者:佚名  約定歌詞  約定  詞:姚若龍 曲:陳小霞  遠處的鐘聲迴蕩在雨裡  我們在屋簷底下牽手聽  幻想教堂裡頭那場婚禮  是為祝褔我倆而舉行  一路從泥濘走到了美景  習慣在彼此眼中找勇氣  累到無力總會想吻你  才能忘了情路艱辛  你我約定難過的往事不許提
  • 海賊王:布魯克和拉布的約定是什麼,路飛是否會履行約定?
    索隆的目標是成為世界第一大劍豪,而布魯克加入路飛是想要完成和拉布的約定,那這個約定是什麼呢?拉布是一條鯨魚,當年倫巴海賊團想要進入偉大航路冒險的時候,因為太過於危險就將拉布放在了顛倒山之前,而布魯克等人和拉布的約定就是等他們繞世界一圈後就會回來,那樣就可以見面了。
  • 《崩壞3》紺海約定活動詳情 紺海約定活動內容獎勵
    導 讀 崩壞3紺海約定活動有什麼獎勵?值得抽嗎?這是崩壞3本周新上線的補給活動之一,小編這裡帶來了崩壞3紺海約定活動介紹,感興趣的玩家們就一起來看看吧!
  • 歌曲《約定》歌詞賞析
    我們約定,難過的往事不許再提,從今以後,我們應該珍惜現在的美好。也答應對方永遠都不讓對方擔心,要做快樂的自己照顧好自己,就算某一天不在對方身邊了。我們約定,一旦爭吵很快要喊停,既然相愛,為何要用爭吵來傷害彼此。
  • 一條道走到黑……跟羅冠聰隔空「約定」後,林夕又跟黎智英約定...
    【環球網報導 記者 趙友平】跟正被通緝的「亂港頭目」羅冠聰隔空「約定」後,香港作詞人林夕又跟亂港頭目黎智英約定周日「深情對談」。林夕12日晚在個人臉書帳號宣稱:約定未改,跟黎智英對談如常舉行。他放出的海報中有「不改約定」「深情對談」字樣,並預告兩人對談的時間為8月16日(星期日)下午,到時會網上直播。香港警方8月10日早上以涉嫌違反香港國安法拘捕黎智英等人,11日凌晨1點多,林夕就跳出來,在臉書上發文支持黎智英,並稱兩人本來在本周六有一個對談節目,「現在看來,希望是個不遠的約定」。
  • 《約定》,天地間最美麗的歌詞
    然而,林夕最優美歌詞並不是以上歌詞,而是為王菲出道時所作的《約定》。忘掉天地仿佛也想不起自己仍未忘相約看漫天黃葉遠飛  就算會與你分離 悽絕的戲  要決心忘記我便記不起  明日天地只恐怕認不出自己  仍未忘跟你約定假如沒有死
  • 《奇蹟暖暖》湖畔的約定高分搭配詳解
    導 讀 奇蹟暖暖湖畔的約定怎麼搭配 第二卷第二章湖畔的約定高分搭配。下面一起來看看吧。
  • SpringBoot約定大於配置到底是什麼意思?
    spring Boot採用約定大約配置的方式,大量的減少了配置文件的使用。使用SpringBoot的同學可能經常看到這幾個字「約定大於配置」,那麼它到底是什麼意思呢?百度一下意思如下:約定優於配置(convention over configuration),也稱作按約定編程,是一種軟體設計範式,旨在減少軟體開發人員需做決定的數量,獲得簡單的好處,而又不失靈活性。
  • 為什麼合同明明約定了,卻依然是套路?
    · 這是優選保公眾號第510篇原創首發文章 ·最近保險市場上出現了一個怪相:各種增額終身壽險泛濫式的接連上市,紛紛都在合同中明確約定了「基本保額每年以xx%遞增」,大部分增額終身壽險的合同中是約定每年3.5%遞增,最高的甚至達到了4%,但是最終給到消費者的收益卻沒有一個能超過3.5%
  • 買賣合同,貨物的交貨時間沒有約定,怎麼辦?
    法律知識要點:買賣合同中的交貨日期是非常重要的合同內容,但是有的當事人在訂立合同時,由於各種因素,沒有約定交貨的日期,或者約定了,但是不能明確交貨的具體日期。例如,合同僅約定三十天內交貨,但是沒有說明是從合同籤訂之日起三十日,還是付款之日起三十日,三十日的起算日期約定不明,容易導致爭議。
  • CF11月靈狐的約定活動網址介紹
    CF11月靈狐的約定活動網址介紹 發布時間:2020-11-18 15:41 來源:52pk 作者:
  • 林夕和被通緝「港獨」羅冠聰,還「約定」上了!
    正被通緝的「亂港頭目」羅冠聰,近日多次在社交媒體上宣稱自己「人在英國,心繫香港」,這一說法還得到了香港作詞人林夕的回應,7日,林夕在個人社交媒體上更改歌曲《約定》的歌詞,稱「還記得當天抗爭多難捱」、「仍未忘跟你約定心念沒有死」,文末還帶上了羅冠聰的話題,並配圖「約定」。
  • 《海賊王》布魯克和拉布的約定,執著的等待感動了許多人
    如果說《海賊王》中哪裡最讓人感動,那麼布魯克和拉布的約定,一定會讓你淚目。50年前,倫巴海賊團與鯨魚拉布相互約定,在環繞偉大航路一圈後,再來到雙子岬與拉布見面。如今整整50年過去了,拉布依然在等待著這個約定,而約定的人也遲遲未到。拉布見不到昔日的夥伴們,便想獨自前往偉大航路,尋找當年約定的夥伴們。
  • 《與神的約定》韓彩英,演技打動了觀眾的心
    昨天(24日)首播的MBC新周末特別企劃電視劇《與神的約定》(編劇洪英熙,導演尹宰文)在製作藝E&M),刻畫了曾經熱愛雙方的韓彩英和裴秀彬離婚,以及因為兒子的病而重逢的故事。人情愛情劇《與神的約定》講述了為挽救垂死的子女而選擇超越世界倫理和道德的兩對夫婦的故事。《與神的約定》每周六晚9點播出。
  • 【最高院•裁判文書】當事人約定對罰息計收複利的,人民法院是否...
    1年期貸款基礎利率加25個基點,按季結息,每季度末月的20日為結息日,21日為付息日;對逾期貸款從逾期之日起就逾期部分按本款約定的罰息利率計收罰息,直至清償本息為止;對借款人不能按期支付的利息以及罰息,以本條第3款約定的結息方式按本款約定的罰息利率計收複利;逾期貸款罰息利率為該浮動利率水平上加收30%;第十二條第一款約定:下列事項之一即構成或視為借款人在本合同項下的違約事件:1.借款人未按本合同的約定履行對貸款人的支付和清償義務
  • 勞務分包合同中,能否約定由勞務承包方提供輔料?
    所分享內容為作者個人觀點,僅供讀者學習參考,不代表本號觀點勞務分包合同中,能否約定由勞務承包方提供輔料?閱讀提示:在勞務分包關係中,勞務承包方通常僅承擔勞務,並計取人工費,即包工不包料,俗稱「清包」或「包清工」。但如果在勞務分包合同中,約定由勞務承包方負責提供輔料,是否會影響合同性質的認定呢?