abaqus平面的案例
Abaqus中平面應力單元高斯積分點的順序
可以輸出umat接口中的變量coords進行查看
write(*,"(A,I4)") "npt = ", npt
write(*,"(A,3ES16.8)") "coords = ", coords
結果為：
npt = 1
coords = -5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 2
coords = 5.77350269E-01 -5.77350269E-01 1.00000000E-02
npt = 3
coords = -5.77350269E-01 5.77350269E-01 1.00000000E-02
npt = 4
coords = 5.77350269E-01 5.77350269E-01 1.00000000E-02
因此Abaqus中平面應力單元高斯積分點的順序為：
展開 Abaqus+PyQt+Python平面變形歐拉角計算
5 相對歐拉角計算
利用計算絕對歐拉角時得到的坐標系文件，計算平面變形相對歐拉角，如下圖所示，計算平面2相當于平面1、平面3相對與平面1的相對歐拉角。
計算結果如下圖所示。
6 小結
上述軟件用的算法申請了發明專利，軟件申請了軟著。CAE工程師，也可以自制軟件工具，解決重復性、復雜性數據處理等工作痛點。
2021年8月24日于西昌衛星發射中心
SolidWorks平面模型導入ABAQUS建立軸對稱模型
（3）有了螺栓截面的草圖，接下來應用曲面工具中的平面工具按鈕，為螺栓零件區域建立截面模型。
圖5
如圖5所示，利用平面工具，根據草圖2生成了螺栓零件的截面模型，這時，在曲面實體下有了相應的截面實體列表。
圖6
重復新建草圖→平面工具過程完成全部6個零件截面建模，結果如圖6。最后可以隱藏草圖1，使得圖形區的圖面顯得較為干凈。
然后可以另存為Parasolid格式的文件，以供ABAQUS導入使用。
（4）如圖7所示，在ABAQUS中作為裝配導入Parasolid文件。在ABAQUS中自動創建了6個零件實例，這樣就可以為每個零件實例劃分網格和賦予材料、建立零件之間的接觸關系，然后加載分析。
ABAQUS導入的面模型默認是在三維空間中，為了分析軸對稱模型，需要回到部件位置對每一個零件編輯，改為軸對稱模型。如圖7所示。
圖7
為了在螺栓上施加預緊力，需要在螺栓桿部適當位置進行一次切分。如圖8所示。
圖8
后面在ABAQUS中的操作都是ABAQUS使用者所熟悉的（賦予材料、建立接觸、添加約束、添加螺栓預緊力等），完善模型后進行分析，結果如圖9所示。
圖9
上述過程還是比較簡單的，ABAQUS使用者有的可能不熟悉SolidWorks的草圖繪制和特征工具的操作，SolidWorks是公認學習曲線非常平緩的軟件，簡單的摸索就能用起來。需要注意的是：要找到SolidWorks里的曲面工具欄；在ABAQUS中導入時注意，要進行接觸分析需要從裝配位置右鍵導入；還需注意默認導入時三維空間（的曲面），要進行平面或者軸對稱分析，需要回到部件位置對每一個部件修改為二維平面或軸對稱，以使得模型的空間維度是正確的。
展開 ABAQUS UEL - 損傷材料本構簡單應用于4節點平面單元 ￥300
利用ABAQUS自定義單元子程序，既可以開發新的單元，同時也可以定義新的材料本構模型。本文以損傷模型簡單應用于4節點平面單元為案例，介紹ABAQUS UEL的開發和使用。
如上圖所示，該單元包含4個節點，每個節點有兩個自由度，分別在水平（X）和垂直（Y）方向運動。節點1的兩個自由度被固定，節點4的水平自由度被固定，節點2的垂直自由度被固定。節點3和節點4在垂直方向上向上運動，位移為0.1mm。該正方形單元的邊長為100mm。在input文件里，坐標表示為，
定義節點組合與邊界條件為，
為了讓模型收斂性更好，采用quasi-newton 求解器。時間步設置為，
在文件夾中通過Powershell提交job和子程序，
單個單元的變形為，
采用不同的 ??
，在后處理中得到損傷因子的變化，
相對應的力-時間關系為，
對于多個單元的情況，比如9單元組成的模型，
具體介紹見知乎：ABAQUS UEL - 損傷材料本構簡單應用于4節點平面單元 - 知乎 (zhihu.com)
相對應的UEL代碼和input文件在付費內容中，
展開 
ABAQUS UEL-梯度損傷模型應用于4節點平面單元 ￥600
本文詳細介紹了如何將梯度損傷模型應用于4節點平面單元，并在有限元模型中進行模擬。
ABAQUS提供了UEL（user defined element）給使用者進行開發。筆者利用UEL開發4節點平面單元，其邊界條件如下圖所示。其中，節點1的X、Y方向被限制住，節點2的Y方向被限制，節點4的X方向被限制，節點3、4的Y方向有豎向位移0.1mm。單元為100*100mm的二維正方形。
每個節點除了X和Y方向的位移，還帶有非局部應變（nonlocal strain）。
單個單元模型，
多個單元模型，
具體內容可參見知乎文章：
ABAQUS UEL-梯度損傷模型應用于4節點平面單元 - 知乎 (zhihu.com)
相應的input文件和uel代碼付費可見，
展開 abaqus中平面應力應變厚度對切削力的影響 ￥5
在鋁合金的二位正交切削仿真中，不同的平面應力應變厚度的對切削力的影響結果
以上為設定值為1的情況
基于abaqus的鋼筋混凝土平面框架倒塌性能分析 ￥100
<p>結構在遭遇偶然突發事件后， 不可避免的會導致結構局部破壞或者損傷， 如果剩余結構不能有效的承擔結構初始破壞和損傷造成的內力變化， 剩余結構就會發生進一步破壞， 造成多米諾骨牌式的連鎖反應，從而造成大范圍嚴重破壞乃至倒塌，這就是通常所說的連續倒塌。附件中只有一個cae有限元模型。</p><p><img src="https://img.jishulink.com/upload/201905/597ca43812cb414e98ab1fd96e276a82.jpg" alt="2019-05-07_105121.jpg"></p><p><img src="https://img.jishulink.com/upload/201905/77bbfce9cdc84985b02d22088e6933bc.jpg" alt="2019-05-07_105131.jpg"></p><p><img src="https://img.jishulink.com/upload/201905/5dd8c1f0916b47cab1cbf1df2992706f.jpg" alt="2019-05-07_105149.jpg"></p><p><img src="https://img.jishulink.com/upload/201905/984ad184d1ba41209ee701b4d0aec1de.jpg" alt="2019-05-07_105208.jpg"></p><p><br></p><p><img src="https://img.jishulink.com/upload/201905/48d8ca94aa6e42768f67ad19803e150b.jpg" alt="2019-05-07_105234.jpg"></p><p><img src="https://img.jishulink.com
展開 基于ABAQUS的UEL子程序定義4節點平面應變等參單元的剛度問題
摘要：
采用基于ABAQUS的UEL子程序開發4節點平面應變等參單元，采用雙線性形函數，4點高斯積分，本構關系為線彈性各向同性材料，得到的單元剛度矩陣和ABABUS自帶的CPE4單元的單元剛度矩陣（剛度矩陣輸出方式為*element matrix output, elset= ALLE, stiffness=yes, OUTPUT FILE=USER DEFINED）不同；對比ANSYS的單元剛度矩陣，結果顯示兩者也不相同。問題出在哪里呢？本文檔將對此問題進行回答。
本文可以作為ABAQUS高級子程序UEL的入門級教程，做UEL的應該關注下！
基于ABAQUS的UEL子程序定義4節點平面應變等參單元的剛度問題（技術鄰 藍牙）.pdf
展開 針對平面應力問題的YLD2000-2D屈服準則及其在ABAQUS中UMAT子程序的實現
Barlat在2003年提出了專門針對平面應力問題的各向異性屈服準則，該屈服準則對于各向異性材料具有很高的精度，得到了廣泛的應用。
YLD2000-2D屈服面示意圖
Yld2000-2d屈服準則由下式給出
其中
矩陣X′和X″的元素分別由柯西應力的下列線性變換獲得
L′和L″的分量由下式求得
積分算法采用徑向返回算法，該方法是穩健和精確的。
當彈性試算超出屈服面時，則需要進行塑性修正
使其滿足
公式9可以通過牛頓法進行迭代求解。
計算的應力應變曲線如下圖所示
B, F. Barlat A , et al. "Plane stress yield function for aluminum alloy sheets—part 1: theory." International Journal of Plasticity 19. 9(2003):1297-1319.
王海波, 萬敏, 閻昱,等. 屈服準則在有限元軟件中實現的正確性驗證[J]. 固體力學學報, 2010, 031(002):173-180.
最后，有需要歡迎通過微信公眾號聯系我們。
微信公眾號：320科技工作室。
展開 誰大牛知道如何讓后半段下降嗎
abaqus平面內受彎實驗。上下固定。左右加相同力 ，abaqus模擬的力-位移曲線一直上升不下降。有人知道怎么調嗎
Abaqus 使用中的一些經驗
ABAQUS對于平面應變和軸對稱的處理，好像不是一般有限元書上講的，對于平面應變，ABAQUS 還是按照3個正應力1個剪應力來做的，當然軸對稱是4個應力分量（不過z方向變成了2方向，而不是三維問題中的3方向），個人認為對于軸對稱不需作特殊處理，只需記住2方向的不同即可；但是平面應變問題，按照ABAQUS ，3方向是z方向，在這個方向上應變始終為0，那么在umat中怎么應力更新，得到非0的應力？在一般有限元書上，E,V都要作一下轉換的（對于平面應變）。
無限元在Abaqus靜力分析中的應用
-02-
Abaqus中的無限元
2.1
單元分類
Abaqus針對平面應力/應變、三維應力以及軸對稱問題均提供了相應的無限元，對于Abaqus/Standard還可以使用帶有減縮積分的無限元。此外Abaqus還提供有聲學無限元。無限元在Abaqus單元庫中遵循如圖4所示的命名慣例。
圖4 Abaqus無限元命名慣例
例如，CIN3D8代表一個8節點線性實體無限元。
2.2
用于靜力分析的無限元的節點定義
由于無限元具有方向性，因此無限元的節點定義方式相比于常規單元有一些特殊的要求。在Abaqus中，無限元的節點編號順序必須使得無限元的第一個單元面（對于平面單元為第一個單元邊）能夠與常規單元的單元面（或單元邊）連接。例如，圖5給出了Abaqus提供的無限元及其節點定義。
（a）平面應力/應變無限元
（b）軸對稱無限元
（c）三維無限元
圖5 Abaqus提供的無限元
例如，對于圖5（a）中的4節點和5節點平面應力/應變無限元，只有4節點無限元的節點1和節點2構成的單元邊才能與常規平面應力/應變單元的單元邊連接；而在5節點無限元中，只有節點1、節點2以及中間節點5構成的單元邊才能與常規的二階平面應力/應變單元的單元邊連接。
此外，在顯式動力學分析中，無限元節點中不屬于第一個面的區域的處理方式與其他分析有所不同。在無限域的方向這些節點會遠離有限元網格。對于顯式分析，這些節點的位置是沒有意義的，并且在顯式動力學分析中不能將載荷和邊界條件賦予到這些節點上。但在通用靜力分析中，這些處于遠離有限元網格區域一側的節點在單元定義中同樣重要，并且能夠賦予載荷和邊界條件。因此如果無限域存在某些邊界條件，如對稱約束，則必須在無限元上也賦予相應的邊界條件。
展開 abaqus形狀記憶聚合物結構的熱-力學有限元模擬
為了研究形狀記憶聚合物相關結構的形狀記憶過程，以往常常需要使用Fortran語言去編寫復雜繁瑣UMAT(用戶材料子程序)，現在本人采用了一種適合對SMP復雜結構進行有限元模擬的方法，該方法不需要寫umat子程序，分別利用有限元仿真軟件ABAQUS中內置的廣義Maxwell模型和Neo-Hookean模型來描述材料的粘彈性行為和超彈性行為。然后針對SMP的板結構，通過ABAQUS軟件對它們的形狀記憶過程進行了有限元模擬分析，得到了應力-應變-溫度三者間的關系。模擬結果表明：本文介紹的這種新方法能夠準確地模擬SMP的形狀記憶過程。
一、SMP熱粘彈性本構模型
根據Tobushi等人的研究，得到了用應力率表示應變率的微分形式的SMP力學一維本構方程：
二、SMP板結構的有限元模擬
1、有限元模型建立
在ABAQUS中建立SMP平面板模型如圖1所示，尺寸為100mm×40mm，選擇Shell進行建模，指定厚度為5mm。網格劃分一共有160個單元，從計算效率考慮，每一個單元尺寸設置為2mm，采用S4R殼單元，即為四節點減縮積分殼單元，計算方式采用Full-Newton求解法。
圖1 有限元模型
在相互作用模塊，需要將板的兩個短邊分別耦合到兩個控制點，控制點與邊之間設置MPC-beam耦合，圖1中的RP-1和RP-2分別為兩邊的控制點。材料屬性設置用到了SMP本構模型。
2、分析過程設置
分為四個步驟：高溫變形、應力凍結、低溫卸載和升溫恢復。具體步驟如下所述：
(1) 初始階段：將RP-1上的U2、U3、UR1和UR3自由度約束住，將RP-2上的U1、U2、U3、UR1和UR3自由度約束住，設置預定義溫度場。
展開 有限元理論基礎及Abaqus內部實現方式研究系列19： Abaqus幾何非線性的設置和后臺
介紹Abaqus的S4單元如何來消除剪切自鎖以及S4R如何來抑制沙漏的。
http://www.yqgqt.org.cn/content/post/350865
1.2 Abaqus的幾何非線性單元設置
為了有限元編程的方便，Abaqus和ansys等商用有限元程序將幾何非線性的設置和小應變/大應變設置完全分開的。
（1）幾何非線性開關的設置：在Abaqus中只要簡單的在Step中勾選NL Geom這個開關就行。
（2）單元的設置：在Mesh的Assign Element Type單元中設置，但有時Abaqus的設置完全看不出到底采用哪種幾何非線性，譬如下面梁的設置，大家能找到哪地方設置小應變和大應變嗎？
而有些時候不得不說Abaqus界面上和后臺計算方法之間的關聯還是非常精準的，譬如殼單元只有membrane strain才有small（小應變）和finite strain（大應變）之分，為什么不直接來一個：Strain選擇small和finite呢？因為Abaqus把平面應變的貢獻分為了兩部分：membrane和bending部分，而后臺計算時只有membrane才是采用的真實應變的方式計算，也就是類似對數方式，而bending部分的應變后臺并不是采用對數方式計算的，這也是為何Abaqus的S4的設置中只有membrane strain的設置原因。
在Abaqus中第一個幾何非線性開關的設置將控制后面的單元設置采用線性和非線性。相當于幾何非線性的開關就是總閘，如果總閘關閉，那么后面的單元設置無論是小應變還是大應變單元都只會采用工程應變來計算，開啟后，才會根據你選擇的小應變還是大應變單元來決定采用Green應變還是真實應變。
展開 基于Matlab的有限元網格自動生成算法 | Q4、Q8、Abaqus單元網格
網格繪制
Q4單元網格生成（橫向排序）
繪圖修飾
眾所周知，Matlab的可視化能力強的一批，接下來木木稍微修改一下patch函數里面的參數，即可更改填充面的顏色以及標記的形狀：
修改繪圖細節
Abaqus-Q4單元
在這個小節，木木想要依據自己的理解，解釋一下:
Abaqus-Q4單元節點編碼順序為什么是1-2-4-3，而不是按照順時針或逆時針的順序編排。
如下圖所示，是在Abaqus中建立一個CPS4(Q4)單元的節點順序：
Abaqus-CPS4單元節點順序
模型節點讀取
在Abaqus建立平面板模型，劃分為
個單元，節點、單元編碼如下圖所示：
Abaqus劃分網格的節點單元編碼
由上圖可知，Abaqus在進行CPS4單元節點編碼時，時一排一排地排序，所以說，在一個單元中節點編碼先是同一方向上的1-2，再是第二行的3-4。
展開 