[發明專利]一種Nurbs基有限元法中多重邊界節點的配點處理方法

有效

• 申請號: 201410230786.6
• 申請日: 2014-05-28
• 公開（公告）號: CN104008242B
• 公開（公告）日: 2017-01-04
• 發明（設計）人: 莫蓉;王文慧;張洪海;高瑾宇
• 申請（專利權）人: 西北工業大學
• 主分類號: G06F17/50
• 分類號: G06F17/50
• 代理公司: 西北工業大學專利中心61204
• 代理人: 陳星

技術領域

本發明涉及有限元技術領域，具體為一種Nurbs基有限元法中多重邊界節點的配點處理方法。

背景技術

在構建零部件三維模型，進行有限元分析時，常常用到等幾何分析技術，該技術是將物理域中的問題轉換為參數域進行求解，如果被分析對象有多個零件則零件的內外邊界可以通過配點法進行處理。用這種方法在處理邊界條件的過程中，可能形成一些其上有多個控制方程的節點，我們稱之為多重邊界節點，比如既是內邊界節點同時也是外邊界節點的節點，以及即屬于一個面也屬于其他多個面的節點，如圖1中所示的節點A和B。

如果所選的配點位于多重邊界節點上，此時待求的離散變量與一個多重邊界節點對應，比如采用開節點矢量表達的兩條NURBS曲線的公共端點。對應于多重邊界節點的離散形式要同時考慮所需滿足的多個方程，比如即需要滿足區域內離散的控制方程，也需要滿足一個或多個邊界區域上的方程。

有兩種傳統辦法來處理這種多重邊界節點。概述如下：一是為每個多重節點變量，根據所需滿足的多個方程推導出唯一一個方程，以該方程應用于多重邊界節點的離散形式，其缺點需要額外的數學推導過程，不能使同一個CAD模型應用于不同的PDE方程。另一種方法是將多重邊界節點近似為多個距離相近的節點，使每個節點對應一個控制方程。然而在等幾何分析中，多重邊界節點通常不宜做這樣的近似，這主要是由于等幾何法中的多重邊界節點有以下特性：1)精確邊界特性。通常由NURBS技術建模的邊界在離散時，可以保持幾何邊界與仿真邊界精確相同，而在采用多個節點逼近多重邊界節點時，幾何邊界的精確性會有損失。2)多重邊界節點通常是參數空間控制矢量端點處的控制變量，如果采用多個節點近似的方法，將會產生尺度非常小的小單元。

發明內容

要解決的技術問題

為適應等幾何分析中多重邊界節點的精確邊界特性和盡量避免出現尺度過小的小單元，便于以統一的方式生成有限元代數方程，本發明提出了一種Nurbs基有限元法中多重邊界節點的配點處理方法，將多重邊界節點上的控制變量確定問題轉換為優化問題，使離散后生成的矩陣方程為過約束矩陣方程，使得多重邊界節點處所對應的多個方程可以在與其他節點上的離散方程一起求解，在簡化了多重邊界節點的離散處理的同時，保證了幾何邊界的精確性。

技術方案

本發明的技術方案為：

所述一種Nurbs基有限元法中多重邊界節點的配點處理方法，其特征在于：采用以下步驟：

步驟1：對零部件中存在多重邊界節點的求解區域進行離散，并進行升階及節點插入，以滿足求解精度的要求；

步驟2：生成求解區域離散后的控制節點，將控制節點按照所對應的基的順序依次連接生成求解區域的控制框架；

步驟3：為每個控制節點配置一個待求的控制變量，并根據控制節點的點、線、面關系，判斷控制變量所對應的空間位置：

如果控制節點位于控制框架的內部，則所對應的控制變量位于求解區域內部；

如果控制節點位于控制框架的某個表面上，則所對應的控制變量位于求解區域的一個表面上；

如果控制節點位于控制框架的多個表面的交集上，則所對應的控制變量位于求解區域的多個表面的交集上；

如果多個控制框架有交集，則其交集中的控制節點所對應的控制變量位于求解區域的內表面上；

步驟4：對于已經分類的控制變量中的內部變量，根據所采取的有限元法，生成所有內部變量所對應的有限元代數方程；

步驟5：采用配點法為所有不屬于多重邊界節點的內表面節點所對應的控制變量生成有限元代數方程；采用配點法為所有不屬于多重邊界節點的外邊界節點所對應的控制變量生成有限元代數方程；

步驟6：根據多重邊界節點所對應的控制變量所屬每個外表面所需滿足的方程及所屬每個內表面所需滿足的連接關系，采用配點法為多重邊界節點所對應的控制變量生成多個所對應的有限元代數方程；

步驟7：將步驟4、步驟5和步驟6中生成的有限元代數方程合并，生成待求解的過約束代數方程組，應用奇異值分解法求解該過約束代數方程組，最終獲得所有控制節點對應的控制變量。

有益效果

本發明的有益效果是，1)實現簡單；2)適用于不同類型的多重邊界節點；3)易于在同一個CAD模型上切換不同的PDE方程，而求解器不需改變；4)所獲得的解，包括與多重邊界節點上所對應的變量的值，是符合所有求解條件的某種最優解。

附圖說明

圖1多重邊界節點示例

圖2問題求解區域及兩種外邊界區域