[發明專利]一種基于多場耦合的加速度計振動時效仿真方法有效
| 申請號: | 201210004049.5 | 申請日: | 2012-01-06 |
| 公開(公告)號: | CN102542117A | 公開(公告)日: | 2012-07-04 |
| 發明(設計)人: | 陳穎;劉丙棟;侯澤兵;康銳 | 申請(專利權)人: | 北京航空航天大學 |
| 主分類號: | G06F17/50 | 分類號: | G06F17/50 |
| 代理公司: | 北京慧泉知識產權代理有限公司 11232 | 代理人: | 王順榮;唐愛華 |
| 地址: | 100191*** | 國省代碼: | 北京;11 |
| 權利要求書: | 查看更多 | 說明書: | 查看更多 |
| 摘要: | |||
| 搜索關鍵詞: | 一種 基于 耦合 加速度計 振動 時效 仿真 方法 | ||
1.一種基于多場耦合的加速度計振動時效仿真方法,其特征在于:該方法具體步驟如下:
步驟一:激光焊接過程仿真;包括:
a.結合焊接殘余應力產生和消除的機理來選擇有限元單元;激光焊接過程中,熱流密度很高的激光束作用于面積很小的焊接部位,使材料局部熔化、汽化,造成被焊部位受熱不均,冷卻過程中的溫度變化產生殘余應力;而VSR消除殘余應力的原理是產生了塑性變形,因此選擇熱結構耦合單元,且需要同時滿足以下條件:(1)三維六面體耦合單元;(2)具有溫度自由度;(3)具有結構自由度;(4)能夠進行瞬態動力學仿真;(5)能夠產生塑性變形;
b.輸入各部件的材料屬性:結合仿真類型輸入仿真需要用到的各部件材料屬性,仿真過程中的材料為石英、低膨脹合金、不銹鋼合金;材料參數包括熱仿真中需要的導熱率,結構和振動仿真中需要的密度、彈性模量、泊松比、熱膨脹系數,為了描述材料的塑性屬性,選用雙線性隨動強化模型,需要輸入最多五個溫度點的材料的屈服極限和屈服后的彈性模量;
c.建立加速度計幾何模型;首先對加速度計的各結構部件進行相應的簡化,只建立被焊接的軟磁體、預負載環,忽略其他的部件;其次,根據結構的對稱性,只建立軟磁體、預負載環的二分之一模型,為了方便將整體劃分為質量較高的六面體單元,將軟磁體切分成幾部分,分別進行建模;
d.劃分網格,形成有限元模型;包括給各幾何部件賦予材料屬性以及利用掃略的方式生成網格;為了生成均勻的六面體單元,將構成結構的特征線劃分成長度大致相同的等份,例如,對于圓柱體,特征線為上下圓面的外環線以及徑向線;
e.施加約束及載荷,包括對稱約束、零位移約束、熱邊界條件、預緊力、熱流密度載荷;在結構對稱面上施加法向對稱約束;在焊接過程中,加速度計的固支部位施加各向零位移約束和預緊力;在所有與空氣相接觸的面上施加對流換熱系數邊界條件;在焊接部位施加熱流密度載荷;
f.進行熱結構耦合仿真,選擇仿真類型為瞬態,設置結果輸出條件,保存數據文件并進行求解;
步驟二:自然冷卻,保留殘余應力,保持模型不變,刪除熱流密度載荷,仿真自然冷卻的過程,仍然進行熱結構耦合瞬態仿真;
步驟三:施加振動應力,進行時域動態仿真;包括:
a.施加VSR過程的約束及載荷;包括刪除激光焊接仿真過程中的零位移約束和預緊力,根據VSR試驗中加速度計的固定方式施加零位移約束,并在固支點上施加隨時間變化的正弦振動載荷;
b.進行振動瞬態動力學仿真,選擇仿真類型為瞬態仿真,打開瞬態效應選項,設定載荷為階躍加載,設置載荷時間和載荷子步數,設置結果輸出條件,保存數據文件并進行仿真計算,求得在施加振動后焊接部位應力變化。
2.根據權利要求1所述的一種基于多場耦合的加速度計振動時效仿真方法,其特征在于:在步驟一e中所述的在焊接部位施加熱流密度載荷,其具體做法的步驟如下:
1)定義熱源中心的x、y、z坐標分別為a、b、c;定義焊接電壓為Q,有效熱半徑為r,則熱流密度為q?max=3*Q/(π*r**2);
2)通過*get命令讀取預負載環外圓周面上單元編號的最小值ne?min和最大值ne?max;
3)循環變量為i,i值從ne?min開始增大,步長為1,讀取每個單元的中心坐標分別為xsy、ysy、zsy,則該單元中心點距離熱源中心的距離為rr=sqrt((xsy-a)**2+(ysy-b)**2+(zsy-c)**2);
4)如果rr≤r,則對該單元施加熱流密度;如果rr>r,則對該單元不施加熱流密度;
5)當i值增大至ne?max時,則結束循環。
3.根據權利要求1所述的一種基于多場耦合的加速度計振動時效仿真方法,其特征在于:在步驟三a中所述的在固支點上施加隨時間變化的正弦振動載荷,其具體做法的步驟如下:
1)定義正弦振動最大振幅為DDAMP,振動頻率為FREQ,振動周期數為NC,每個振動周期上取NN個點,則每個載荷步的時間間隔為dt=1/FREQ/NN,總的載荷步數為NT=NC*NN;
2)設置振動開始時間為tt,循環變量為i,i值從1開始增大,步長為1,則每經過1次循環時間為tt=tt+dt,該時刻振動的幅值為dd=DDAMP*sin(FREQ*2*π*tt);
3)選擇需要施加載荷的節點,并以斜坡方式施加振動應力dd,求解;
4)當i值增大至NT時,則結束循環。
該專利技術資料僅供研究查看技術是否侵權等信息,商用須獲得專利權人授權。該專利全部權利屬于北京航空航天大學,未經北京航空航天大學許可,擅自商用是侵權行為。如果您想購買此專利、獲得商業授權和技術合作,請聯系【客服】
本文鏈接:http://www.szxzyx.cn/pat/books/201210004049.5/1.html,轉載請聲明來源鉆瓜專利網。





