abaqus應力方向的案例
有限元理論基礎及Abaqus內部實現方式研究系列14： 殼的應力方向
我們關注CAE中的結構有限元，所以主要選擇了商用結構有限元軟件中文檔相對較完備的Abaqus來研究內部實現方式，同時對某些問題也會涉及其它的Nastran/Ansys等商軟。為了理解方便有很多問題在數學上其實并不嚴謹，同時由于水平有限可能有許多的理論錯誤，歡迎交流討論，也期待有更多的合作機會。
iSolver介紹視頻：
http://www.yqgqt.org.cn/college/video/c12884
==第14篇：殼的應力方向 ==
有限元中，物理量用的最多的是標量、矢量和二階張量。其中位移、坐標等都是矢量，而應變、應力等都是二階張量。矢量很容易理解，體的應力等二階張量直接就采用了全局坐標系的也不會有方向理解問題，但梁殼的應力結果很容易搞錯，后處理結果中的S11、S12等的方向有時會覺得和預想的不一致但又不明所以。同時，這個方向也是單元材料的方向，在自編程序時，如果一開始坐標系的定義就弄錯了，那么將直接導致和材料相關的剛度矩陣的錯誤，所以弄清應力的方向定義對自編程序和理解有限元結果都相當重要。
本章將簡單介紹一下數學上張量和Abaqus中殼的應力方向，并說明Abaqus這么選取的意義，最后通過自編程序iSolver來驗證殼的應力方向的正確性。
具體的驗證詳見下方視頻（帶配音）：
https://www.yqgqt.org.cn/college/video/c12884 20.14 理論系列文章14：殼的應力方向
1.1 數學上的張量方向
矢量的方向是一定的，但它的分量都是基于某個坐標系定義的，坐標系不同，那么分量結果也會不同。
矢量可以表示為：
顯然，分量和坐標系的選取有關。譬如我們一般的直角全局坐標系如下，那么分量就是普通的x、y、z三個分量值。
展開 ABAQUS喵星人教你看懂不同類型單元的應力方向
<p><span style="color: rgba(0, 0, 0, 0.9);">應力為典型的張量，具有明顯的坐標相關性，大家常用查看單元應力方向的方法為直接通過整體坐標系判斷XYZ方向，但這種方法僅適用于實體單元，對于其他類型單元（例如殼單元、Beam單元、Truss單元、Cohesive單元等）或特殊坐標系下的實體單元則不再適用，若仍然采用整體坐標系判定方向則會限制對后處理結果的解讀。今天喵星人就通過一個教程帶大家學習不同類型單元的應力方向應該如何看。</span></p><p class="ql-align-center"><strong>1.實體單元</strong></p><p><br></p><p>默認的實體單元應力方向服從整體坐標系，若想查看其他坐標系下的應力情況則需定義其他坐標系，建立的方式既可在前處理內定義，也可在后處理內完成，前后處理中坐標系的定義位置如下圖所示。
展開 workbench怎么展示多方向的結果應力圖
問題：workbench在出圖的時候，不管是模型示意圖，載荷施加圖還是應力結果圖，有時候一張圖顯示不全，效果不是很直觀，就需要多個方向的去截圖顯示。
如圖所示，點擊紅圈部分，會有one，horizontal，vertical，four四種顯示模式，可以是默認的一個界面，橫向兩個界面，豎向兩個界面，四個界面。這樣可以多方向的展示你的結果，每個圖框可以調整方向，還可以單獨做個截面等。
comsol怎么顯示單獨y方向上的應力分布
大佬們有人知道怎么顯示軸向的應力分布嗎，comsol上結構力學表達式太多，不知道選哪一個
Ls-Dyna復合材料任意主方向定義（類似Abaqus離散化方向定義） ￥9.9
<p>對于擁有復雜曲面結構的復合材料薄板，通常需要定義一個變化的材料主方向，下面介紹在Lspp中如何定義。</p><ul><li>對于任意復雜結構的平面，劃分網格后，每個網格的方向是根據節點坐標得到的，總體上呈現隨機性。</li></ul><div contenteditable="false" width="100%">
<figure class="figure-image" data-img="https://img.jishulink.com/202410/attachment/1c788f57a7554bab9067a3554e8759b0.png" style="text-align: center" data-regular="true">
<img src="https://img.jishulink.com/202410/attachment/1c788f57a7554bab9067a3554e8759b0.png" data-mobile-src="https://img.jishulink.com/202410/attachment/1c788f57a7554bab9067a3554e8759b0.png?image_process=/format,webp" data-pc-src="https://img.jishulink.com/202410/attachment/1c788f57a7554bab9067a3554e8759b0.png?
展開 原創：第一方向殘余應力結果中名義厚度0.969含義。
在Moldflow軟件中，1代表產品表面，0是產品中心層，Normalized thickness =0.969，就是非常靠近表皮的斷面位置，只差0.31%的厚度位置。
不了解的話，看下面的圖片說明。Additional information, Why 0.969, not 1 or 0.938，0.969 is middle between 1 and 0.938.(也有人問過這個問題)
【abaqus】個人筆記—應力奇異&應力平均&應力集中
【abaqus】個人筆記—應力奇異&應力平均&應力集中
ABAQUS熱應力分析 附ABAQUS中初始地應力的施加下載
軋輥與Cu層的熱傳導系數
下載地址：ABAQUS中初始地應力的施加
ABAQUS中殼的材料方向
ABAQUS中殼的材料方向
當結構一個方向的尺度（厚度）遠小于其它方向的尺度，并忽略沿厚度方向的應力時，可以用殼單元進行模擬。在ABAQUS中具有兩種殼單元：常規的殼單元和基于連續體的殼單元。
與實體單元不同，每個殼體單元都使用局部材料方向。
1、默認的局部材料方向
局部材料的1和2方向位于殼面內，默認的局部1方向是整體坐標1軸在殼面上的投影，如果整體1軸垂直于殼面，則將整體3方向投影到殼面形成1方向，殼面的正法線方向為3方向，對于殼面內的2方向，利用3x1=2方向（3方向叉積1方向）確定。即局部的1、2、3方向構成右手坐標系。
然而，在更多的情況下，利用默認的局部材料設置并不能順利完成定義，尤其是對于曲面、圓筒等結構，而此時就需要利用其它方法定義合適的材料方向。
2、可變的材料方向
應用局部的直角、圓柱或者球坐標系，可以代替整體坐標系，如下圖所示。定義局部坐標系（x',y',z'）的方向，并使局部坐標軸的方向與材料方向一致。為此，必須先指定一個最接近垂直于殼體的局部軸，以及繞該軸的旋轉量（如果需要）。ABAQUS按照坐標軸的循環順序（1,2,3）及用戶的選擇將坐標軸投影到殼體上，從而構成材料的1方向。例如，如果選擇了x'軸，ABAQUS將y'軸投影到殼體上而構成材料的1方向。由殼法線和材料1方向的叉積來確定2方向。
如果這些局部坐標軸沒有建立理想的材料方向，就需要用到前面設置的繞軸轉動了。在將軸投影前，先按照該轉動量進行轉動，然后投影得到最終的局部材料方向。
展開 Abaqus帶螺紋螺栓接觸應力分析淺析 Abaqus帶螺紋螺栓接觸應力分析淺析
目前的常規做法通常有兩種：1.簡化，用RBE2和beam梁來代替螺栓，這樣不能反映連接螺栓真實應力，圖1為某結構連接螺栓簡化的beam梁應力云圖，沒有接觸應力：
.直接做出來螺栓螺紋采用接觸分析，雖然得出的結果很精確，但這樣前處理工作量大（螺栓和螺紋用六面體網格建模）、計算量大（接觸收斂困難），如圖為某結構帶螺紋螺栓和連接件模型（圖2）和計算得出的結果（圖3）：
圖3 計算結果
那么，有什么好辦法可以不用簡化帶螺紋螺栓，不用直接做出帶螺紋螺栓，又能得到足夠精確的結果？
運用大型通用非線性有限元分析軟件Abaqus，只需要在接觸定義中設置跟實際螺紋形狀有關聯的參數，如牙角、螺距、螺栓小徑等，就可以模擬真實的連接螺栓接觸狀況。既可以得到足夠精確的分析結果，又節省了時間專注進行其他的分析設置。如圖4，為連接螺栓接觸來定義帶螺紋螺栓：
圖4 連接螺栓接觸來定義帶螺紋螺栓
圖5為某結構直徑10MM的帶螺紋的連接螺栓接觸壓力分布云圖：
圖5 某結構直徑10mm帶螺紋的連接螺栓接觸壓力分部云圖
展開 ABAQUS中殼的材料方向
當結構一個方向的尺度（厚度）遠小于其它方向的尺度，并忽略沿厚度方向的應力時，可以用殼單元進行模擬。在ABAQUS中具有兩種殼單元：常規的殼單元和基于連續體的殼單元。
與實體單元不同，每個殼體單元都使用局部材料方向。
1、默認的局部材料方向
局部材料的1和2方向位于殼面內，默認的局部1方向是整體坐標1軸在殼面上的投影，如果整體1軸垂直于殼面，則將整體3方向投影到殼面形成1方向，殼面的正法線方向為3方向，對于殼面內的2方向，利用3x1=2方向（3方向叉積1方向）確定。即局部的1、2、3方向構成右手坐標系。
然而，在更多的情況下，利用默認的局部材料設置并不能順利完成定義，尤其是對于曲面、圓筒等結構，而此時就需要利用其它方法定義合適的材料方向。
2、可變的材料方向
應用局部的直角、圓柱或者球坐標系，可以代替整體坐標系，如下圖所示。定義局部坐標系（x',y',z'）的方向，并使局部坐標軸的方向與材料方向一致。為此，必須先指定一個最接近垂直于殼體的局部軸，以及繞該軸的旋轉量（如果需要）。ABAQUS按照坐標軸的循環順序（1,2,3）及用戶的選擇將坐標軸投影到殼體上，從而構成材料的1方向。例如，如果選擇了x'軸，ABAQUS將y'軸投影到殼體上而構成材料的1方向。由殼法線和材料1方向的叉積來確定2方向。
如果這些局部坐標軸沒有建立理想的材料方向，就需要用到前面設置的繞軸轉動了。在將軸投影前，先按照該轉動量進行轉動，然后投影得到最終的局部材料方向。
abaqus中殼的局部材料方向.pdf
展開 
abaqus反應譜計算的方向余弦問題
請問一下各位大佬，abaqus用振型分解反應譜法進行動力分析時，方向余弦填1,0,0和-1,0,0為什么得出來的結果是一樣的？不應該相反嗎
ABAQUS-真實應力和名義應力轉化
ABAQUS-真實應力和名義應力轉化.doc
有限元理論基礎及Abaqus內部實現方式研究系列38： 梁單元差異（2）-梁截面方向
梁局部坐標系的x方向永遠都是1->2節點，abaqus中稱為t方向，如下方：
設置梁方向時，輸入n1的一個三維方向矢量，簡單起見，n1和t直接取垂直，取默認的0,0，-1：
Abaqus后臺得到
局部坐標系的z方向（即截面2的方向）Abaqus.2= t×n1
局部坐標系的y方向（即截面1的方向）Abaqus.1=Abaqus.2×t
最終t,S.1,S.2滿足右手定則，得到局部坐標系方向
三維顯示為：
2.2.2 Abaqus梁截面幾何尺寸的設置方向
很簡單，梁截面幾何尺寸的設置方向的1、2就是Abaqus的局部坐標系的y、z軸。
2.3 Nastran的梁截面方向
2.3.1 Nastran梁截面屬性關聯的局部坐標系方向
Nastran的局部坐標系的x方向和Abaqus完全一致，都是節點1->2方向t。和Abaqus一樣的模型如下：
Patran中可以設置Bar Orientation，在Nastran中稱為v方向，同樣設置為0,0，-1：
按照Nastran幫助手冊，Patran設置的局部坐標系和Abaqus完全一致，Nastran由t和v首先確定一個局部坐標系Patran的y，z方向，此時
Patran.z= t×v
Patran.y=Patran.z×t
即下圖表示：
另一種說法是先定義Plane1（即局部xy平面）就是v和t所在平面，Plane2垂直與Plane1。其實就是上面的后臺公式
那既然局部坐標系和Abaqus完全一致，那么Abaqus的L型定義的參數輸入到Patran中是否在三維全局坐標系下也完全一致呢？
展開 有限元理論基礎及Abaqus內部實現方式研究系列15： 殼的剪切應力
我們關注CAE中的結構有限元，所以主要選擇了商用結構有限元軟件中文檔相對較完備的Abaqus來研究內部實現方式，同時對某些問題也會涉及其它的Nastran/Ansys等商軟。為了理解方便有很多問題在數學上其實并不嚴謹，同時由于水平有限可能有許多的理論錯誤，歡迎交流討論，也期待有更多的合作機會。
iSolver介紹視頻：
http://www.yqgqt.org.cn/college/video/c12884
==第15篇：殼的剪切應力 ==
自編有限元應力的校核除了Mises等合力外，也應該校核各個應力分量。材料力學中六個應力分量如下：
其中Tau11，Tau22，Tau33為正應力，Tau12,13,23為三個剪切應力，對殼來說，Tau33=0，Tau12為面內剪應力，Tau13,23即為本文所說的橫向剪切應力。
最近在做iSolver殼的應力分量和Abaqus比對時，發現Abaqus的橫向剪切應力和預想的不一致。iSolver按照常用的殼的理論得到的剪切應力是個與厚度無關的常量，但Abaqus的橫向剪切應力分量TSHR13，TSHR23，在各個截面方向積分點section point不一樣。
花了點時間細致的研究了一下，猜測Abaqus中剪切應力TSHR13、23是真實應力，但有限元理論和iSolver中計算的是板殼近似理論中平均剪切應力。本章將介紹殼單元中實際的和板殼近似理論中的剪切應力，也猜測了Abaqus的內部實現流程，最后通過一個算例來驗算Abaqus中的真實的剪切應力，并通過iSolver來計算板殼理論的平均剪切應力。
1.1 殼的真實的剪切應力
剪應力是材料由于抗拒面之間的滑動而產生的沿表面方向的應力。殼的中間層存在剪切應力，這個可以通過下面簡單的例子驗證。
展開 