ansys節點信息函數的案例
從形函數與函數的連續可導性到ansys結果中的節點解與單元解的差異
如題，《從形函數與函數的連續可導性到ansys結果中的節點解與單元解的差異》，形函數對結果的影響大部分人都能聯想到二次單元比線性單元求得的結果更精確，但該文要表達的不僅如此，而是從更一般地討論怎么從單元的形函數來理解節點解與單元解之間的差異。
首先討論單元的階次。作為基礎我們應該明白網格與單元的區別，網格是將幾何體離散化后的結構，即組成幾何體的微元，單元是這些微元的幾何、物理或數學屬性（這里我們并不打算詳細討論單元的這些屬性，但是這些知識會方便對本文的理解）。我們經常在使用ansys或其他CAE軟件時經常會遇到單元的選擇以及單元階次的選擇，一般一種單元包括線性單元和二次單元甚至更高級的單元，比如在ansys中經常被使用的shell181（左）和shell281（右），線性單元使用的形函數是一次的多項式，高次單元使用的形函數是高次的多項式，形函數用于描述相鄰節點之間的位移場，所以高次的單元可以更好的描述形狀復雜的幾何體。
不同于常規材料力學中通過平衡方程求解（首先求得的解是力解），有限元方式求解的特點是首先求解出的結果是節點的位移解，即displacement of nodes，所有的節點位移形成了位移場，在空間上位移場一定是連續的，但是不一定是平滑的。哎哎，是不是特別熟悉的感覺，正是和高數中函數的連續性和可導性兩個性質非常相似，不用奇怪，位移場本來就是用函數描述的，所以自然就存在函數的性質，所以用函數的性質來理解就可以方便解釋一些現象了，下圖分別是用兩種形函數描述的位移場，在有限元求解后得到的首先是節點位移解，即圖中5個節點的位移，假如每個節點的位移用坐標x\y\z的函數來表示，然后通過形函數插值得到相鄰節點之間的位移（也是xyz的函數），上圖是用一次形函數插值，下圖是用二次形函數插值。
展開 ANSYS輸出實體模型表面的節點信息 和單元拓撲關系
ANSYS輸出實體模型表面的節點信息
和單元拓撲關系
遇到一個問題，一個給定的實體模型，劃分了solid185的單元，假如實體模型單元劃分如下。需要提取實體模型外表面節點位置信息和單元拓撲關系（也就是每一個單元是由哪幾個節點組成的），目的是方便做其他分析，比如流體分析，提取外表面的節點可以施加溫度載荷。
圖1
對于此問題，在ansys里面很難直接提取所有外表面的節點和單元信息，因為外表面也是實體單元的一個單元面，不可能剝離出來。
因此，想要提取外表面的單元和節點，最好是需要外表面存在平面單元。
對于此，可以采用ansys里面的特殊單元mesh200，這個單元用于面網格的劃分，而且劃分后的單元不參與實際計算。
于是：
et,2,200 ！定義mesh200單元類型
asel,s,ext ！選擇所有的外表面
aatt,,,2 ! 設置劃分單元為mesh200
KEYOPT, 2, 1, 6 ! 4節點的四邊形單元
amesh,all ! 劃分所以的外表面
此時劃分的面網格和原來的實體網格的節點是一一對應的，這就保證了最后輸出的節點的坐標與原來實體模型的對應節點是一一對應的。
此時可以選擇刪除實體模型和實體單元。
展開 