[發明專利]基于CAE分析結果的應力場插值方法

在審

• 申請號: 202011331720.8
• 申請日: 2020-11-24
• 公開（公告）號: CN112417693A
• 公開（公告）日: 2021-02-26
• 發明（設計）人: 徐桑迪;陳睿;黃清華
• 申請（專利權）人: 江蘇華路信息科技有限公司
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06F119/14
• 代理公司: 無錫市匯誠永信專利代理事務所(普通合伙) 32260
• 代理人: 王春麗

1.基于CAE分析結果的應力場插值方法，其特征在于，包括如下步驟：

S1，利用CAE軟件對幾何模型進行第一次網格劃分，并進行第一次應力分析，獲得第一次應力分析結果；

S2，對步驟S1的第一次應力分析結果的收斂性進行判斷，如果所述第一次應力分析結果收斂，導出第一次應力分析結果和第一劃分模型；如果所述第一次應力分析結果不收斂，返回步驟S1；

S3，利用CAE軟件對步驟S2中的第一分析模型進行再次網格劃分，再次網格劃分中網格單元的尺寸小于第一次網格劃分中網格單元的尺寸，導出第二劃分模型；

S4，對步驟S2中的第一劃分模型利用單元形函數計算方法獲得粗網格模型的單元積分點坐標及應力值，對步驟S3中的第二劃分模型利用單元形函數計算方法獲得細網格模型的單元積分點坐標；

S5，通過空間應力場插值計算方法對步驟S4中的粗網格模型的單元積分點坐標及應力值、細網格模型的單元積分點坐標進行插值計算，獲得細網格模型的單元積分點的應力值；

S6，將步驟S3中的第二劃分模型與步驟S4的初始溫度場結果進行溫度場的加載計算，獲得第二劃分模型的應力場映射結果。

2.根據權利要求1所述的基于CAE分析結果的應力場插值方法，其特征在于，在步驟S5中，所述空間溫度場插值計算方法包括如下步驟：

(1)設置插值參數k，k為對步驟S3中第二劃分模型的節點搜索近鄰點的數量；

(2)計算比例系數r，比例系數r＝步驟S3中再次網格劃分的網格單元數量/步驟S1中第一次網格劃分網格單元數量；

(3)當0r4時，t＝0.53，當r≥4時，t＝1，t為步驟S3中再次網格劃分的網格單元平均尺寸s的系數

(4)設置插值參數R，R為以步驟S3中第二劃分模型的節點為球心的球體的半徑；

(5)判斷步驟S3中第二劃分模型的節點的近鄰點與步驟(4)球體的位置關系，當近鄰點與球心重合，即與步驟S2中的第一劃分模型的節點重合時，步驟S3中第二劃分模型的積分點應力值即為該近鄰點的應力值；當近鄰點在球體邊緣或這外部時，剔除該近鄰點；當近鄰點在球體內部時，保留進行后續的空間插值；

(6)判斷步驟(5)中保留的近鄰點是否在步驟S3中第二劃分模型的節點的同一側，所述即步驟S3中第二劃分模型的節點p(x,y,z)和所有近鄰點p(x_i,y_i,z_i)是否滿足(x-x_i)或(y-y_i)或(z-z_i)都為同號，如果為同號，則進行外插，如果為異號則保留進行內插，i＝1,2,…,k；

(7)計算步驟S3中第二劃分模型的節點p(x，y，z)與所有近鄰點p(x_i，y_i，z_i)的歐式距離

(8)計算插值權重比

(9)步驟S3中第二劃分模型的節點應力

3.根據權利要求2所述的基于CAE分析結果的應力場插值方法，其特征在于，步驟(1)中，所述k取值4或5。

4.根據權利要求2所述的基于CAE分析結果的應力場插值方法，其特征在于，步驟(4)中，所述R的值為所述t與所述s的乘積。