ANSYS輸出實體模型表面的節點信息

和單元拓撲關系

遇到一個問題，一個給定的實體模型，劃分了solid185的單元，假如實體模型單元劃分如下。需要提取實體模型外表面節點位置信息和單元拓撲關系（也就是每一個單元是由哪幾個節點組成的），目的是方便做其他分析，比如流體分析，提取外表面的節點可以施加溫度載荷。

                                               

圖1

對于此問題，在ansys里面很難直接提取所有外表面的節點和單元信息，因為外表面也是實體單元的一個單元面，不可能剝離出來。

因此，想要提取外表面的單元和節點，最好是需要外表面存在平面單元。

對于此，可以采用ansys里面的特殊單元mesh200，這個單元用于面網格的劃分，而且劃分后的單元不參與實際計算。

于是：

et,2,200                   ！定義mesh200單元類型

asel,s,ext                   ！選擇所有的外表面

aatt,,,2                             ! 設置劃分單元為mesh200

KEYOPT, 2, 1, 6      ! 4節點的四邊形單元

amesh,all                   ! 劃分所以的外表面

 

此時劃分的面網格和原來的實體網格的節點是一一對應的，這就保證了最后輸出的節點的坐標與原來實體模型的對應節點是一一對應的。

此時可以選擇刪除實體模型和實體單元。

刪除后，顯示的就只有外表面的單元和節點，通過菜單欄List下面的Nodes，可選擇列表輸出節點坐標，如下所示：

圖2

同樣，可通過List輸出單元拓撲信息，如下：

先用命令：

numcmp,all

壓縮所有的單元和節點的編號，也就是重新編號，否則可能不是從1開始。

最后顯示如下：

圖3

這兩個文件都可以保存，之后需要可以再寫程序提取。

這里有一個問題，現在得到的表面的節點號和原來實體模型對于位置的節點號不是對應的，處理這個問題需要重新寫程序，用什么語言都可以，Python，C++等等，目的是讀取ansys輸出的節點信息文件，讀出固定坐標處對應的節點號，通過對比可以找到所以的節點對應關系。