abaqus無限元劃分的案例
無限元在Abaqus靜力分析中的應用
為了便于數值積分，將物理坐標系(x, y)中的剛度矩陣變換到自然坐標系(ξ, η)中進行計算，即有：
式中：|J2|為2階Jacobian行列式，這里用到了物理坐標與自然坐標之間的變換關系：
這里，與普通等參元不同的地方只是Jacobian矩陣的求解：
同樣基于MATLAB的符號變量求偏導功能，可以得到各個節點的映射函數Mi對自然坐標ξ和η的偏導分別為：
因此，基于Gauss-Legendre積分公式，可以將單元剛度矩陣表示為：
因此單元的剛度方程為：
式中：Pe為等效節點載荷矩陣，可表示為：
式中：Pxi和Pyi分別代表作用在節點i上沿x方向和y方向的節點力。
-02-
Abaqus中的無限元
2.1
單元分類
Abaqus針對平面應力/應變、三維應力以及軸對稱問題均提供了相應的無限元，對于Abaqus/Standard還可以使用帶有減縮積分的無限元。此外Abaqus還提供有聲學無限元。無限元在Abaqus單元庫中遵循如圖4所示的命名慣例。
圖4 Abaqus無限元命名慣例
例如，CIN3D8代表一個8節點線性實體無限元。
2.2
用于靜力分析的無限元的節點定義
由于無限元具有方向性，因此無限元的節點定義方式相比于常規單元有一些特殊的要求。在Abaqus中，無限元的節點編號順序必須使得無限元的第一個單元面（對于平面單元為第一個單元邊）能夠與常規單元的單元面（或單元邊）連接。例如，圖5給出了Abaqus提供的無限元及其節點定義。
展開 關于ABAQUS無限元邊界
大家修改完文件里面的單元類型后，容易出現單元節點順序不對的情況，這里只需要大家在劃分網格的時候，使用掃掠，并且需要定義掃掠的路徑！！！！！路徑非常重要，路徑的方向為朝外，這樣就不會出現單元節點順序錯誤的情況，也不用很麻煩的一個個的節點修改。
我還會做一些盾構開挖，抗震，基坑模擬等等，需要的給我留言。
