abaqus平面問題計算的案例
Abaqus+PyQt+Python平面變形歐拉角計算
5 相對歐拉角計算
利用計算絕對歐拉角時得到的坐標系文件，計算平面變形相對歐拉角，如下圖所示，計算平面2相當于平面1、平面3相對與平面1的相對歐拉角。
計算結果如下圖所示。
6 小結
上述軟件用的算法申請了發明專利，軟件申請了軟著。CAE工程師，也可以自制軟件工具，解決重復性、復雜性數據處理等工作痛點。
2021年8月24日于西昌衛星發射中心
基于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計算內存的問題：aba_sint_cap
因此，查閱了網絡上相關帖子和幫助文檔，加上親身經歷，整理如下：
方法一：可以直接在abaqus_v6.env中設置memory='80%'或其他方式，詳見幫助文檔
方法二：可通過高級系統設置--->環境變量--->變量名：ABA_SINT_CAP 值：16384
（說明：16384這個值可以更大，但似乎沒有意義！）
基于上述方法設置了足夠大的內存后，但仍然出現該錯誤。具體原因：
Abaqus中每8GB的物理內存，使用的內存將不會超過2GB
在滿足上述條件的基礎上，最多使用內存限制為16GB。
遇到這個問題后，建議重新調整模型，減小inp文件的體量。目前尚未找到特別好的方法！
展開 
ABAQUS建模計算的常見問題（二）
1.DAT文件的錯誤與警告信息說明
當ABAQUS無法順利完成模擬分析任務時，系統會在MSG文件中會出現對應的錯誤位置說明，以供操作者查看檢測。
具體的錯誤類型有：零主元、負特征值、過多迭代次數和數值奇異等。
2.錯誤舉例
零主元
解決方法：出現零主元（zero pivot）的最可能原因是模型中存在過約束。這時就要調整模型，可以從邊界條件、接觸關系設置等著手。
（2）負特征值
解決方法：負特征值（Negative Eigenvalue）的出現，可能是由于單元發生了過度的變形、剛體位移沒有消除，或者應力應變關系中存在負斜率等。需要注意：負特征值出現并不一定意味著模型不對，只要該提示不是出現在最后一次迭代，計算分析就沒錯；否則，需要調整模型。
（3）過多迭代次數
解決措施：過多次迭代出現的原因是，軟件在計算時，會按照設定的增量進行迭代，達到一定次數如果不能收斂時，增量就會自動減小；隨后重復迭代；多次“減小”之后還不能收斂，計算就會停止。這種情況比較常見，往往是由于同時還存在這其他的警告或錯誤信息。可以先嘗試解決其他的問題。
（4）數值奇異
解決措施：數值奇異（Numerical Singularity）時，一般是因為模型中存在剛體位移。所以需要檢查模型的約束，包括邊界約束和各個部件的約束。.
（5）局部塑性變形過大
解決措施：此錯誤的出現并不一定意味著模型不對。若是出現收斂問題，應該查看是否局部有過大的變形。在一些分析中，比如開孔的鋼板，孔的邊緣容易出現過大的塑性變形，如果這部分不屬于分析的關鍵，可以局部提高材料性質。
來源：土木CAE筆記
展開 ABAQUS建模計算的常見問題（一）
若無，則說明已經消除了過約束，計算正常；否則，則要調整模型。
參考文獻：基于ABAQUS的有限元分析和應用 莊茁等。
來源：土木CAE筆記
abaqus反應譜計算的方向余弦問題
請問一下各位大佬，abaqus用振型分解反應譜法進行動力分析時，方向余弦填1,0,0和-1,0,0為什么得出來的結果是一樣的？不應該相反嗎
ABAQUS橡膠工業中有限元計算問題過盈配合
摘要：本文分析了過盈配合的有限元計算時用到超彈性本構時可壓縮性對計算結果的影響情況，得到在過盈配合中必須考慮這種可壓縮性的結論并分析考慮可壓縮性的原因。
1、引言
過盈配合是橡膠工業中的一種常見的配合方式。橡膠為超彈性材料，有限元計算中通常假定為不可壓或者幾乎不可壓。本文首先給出一種不可壓橡膠模型過盈配合的理論解，并與ABAQUS 計算解進行比較。進一步本文探討過盈配合中假定橡膠不可壓時遇到的問題， 提出處理過盈配合中橡膠計算的方法。
2、可壓模型理論解與 ABAQUS 數值解的比較
2.1、理論解
理論解模型如圖 1，內層為鋼，中間不可壓橡膠，最外層為鋼給出橡膠和橡膠之間的過盈量求整個結構的應力應變狀態假設平面應變狀態。
本構方程：
對于鋼：
對于橡膠：
2.1 材料性質：
鋼：E=210000v=0.3橡膠：C10=0.461312, C20=0.01752, C30=8.8e-05 ，其余為 0，（三次多項式模型，材料不可壓縮）
2.2.2 幾何特性
如圖 2 所示， R59.50 為內層鋼的半徑和中間層橡膠的內徑， R73.00 為中間層橡膠的外徑，R71.10 為外層鋼的內徑， R80.00 為外層鋼的外徑。
展開 