基于PATRAN/NASTRAN和ADAMS的聯合仿真分析

田路路 朱東方 孫祿君

摘 要 本文將patran/nastran和ADAM分析軟件結合起來對柔性體進行仿真分析， 論述了多柔體動力學的理論基礎、柔性體模型的定義方法、柔性體動力學特性的仿真流程。以太陽能電池陣等效模擬器為例，詳細說明了聯合仿真過程。

【關鍵詞】patran/nastran adams 聯合仿真 柔性體

ADAMS是主要以剛性體為主要分析對象的運動學仿真分析軟件。在對機械系統進行運動學分析時，在較大的載荷作用下或者作加、減速運動時， 柔性機構就會產生較大的彈性變形。因此要更真實地模擬機構的動態性能， 必須將構件作為柔性體進行分析。Patran/nastran是以柔性體為主要分析對象的有限元分析軟件，因此， 要精確模擬整個系統的運動， 取得動力學仿真結果，可以將Patran/nastran與ADAMS軟件結合使用。

1 柔性體的運動學方程和系統動力學方程

1.1 柔性體的運動學方程

為慣性坐標系，為物體的動坐標系，p為混合坐標系中的任意一點，則，p在慣性坐標系中的位移可表示為，其中，為動坐標系在慣性坐標系中的位置，為點p在中未變形時的位置，為p在中的變形，A表示動坐標系到慣性參考坐標系的轉換矩陣。物體坐標系的模態坐標用表示，表示對應于結點p的移動自由度的模態矩陣子塊，則，代入即有。

1.2 柔性體系統的動力學方程

運用拉格朗日乘子法可得柔性體的動力學微分方程

，其中：M為柔性體的質量，W為勢能；K和D分別為模態剛度矩陣和模態阻尼矩陣；和分別表示物體內部由于彈性變形和阻尼引起的廣義力；ξ為選定的廣義坐標，Q為投影到ξ上的廣義力；為約束方程，λ為對應于的拉格朗日乘子。

2 PATRAN/NASTRAN和ADAMS聯合仿真的步驟

patran/nastran中生成模態中性文件：

以patran/nastran中裝配體生成模態中性文件為例：

（1）在patran菜單欄中選擇File-import，彈出對話框，選擇parasolid文件類型，導入*.x_t的文件，刪除無關的零部件；定義材料屬性，輸入材料屬性，如泊松比、彈性模量、密度等；隨后選擇單元類型。之后定義裝配體中零部件剛/柔性。

（2）網格劃分：選擇Meshing-meshers進行網格劃分。隨后進行創建載荷/邊界條件定義六個自由度。

（3）仿真分析：

設置分析類型：單擊Solution Type按鈕，打開其對應的面板，選擇分析的類型NORMAL MODES，設置ADAMS output 選項為Full Run+MNF，設置單位（此單位需和ADAMS中所用單位一致），單擊OK。

設置裝配體各零部件之間的連接關系：若將零部件i和零部件j之間的連接關系用（）來表示，其中n表示零部件個數，表示T或G（TUCH或GLUE）。單擊Analysis-Subcases，選擇彈出對話框中contact Table按鈕，設置零部件之間關系。

完成以上設置后，開始計算求解，Nastran計算后，調入XDB文件，在右側操作面板區選擇點擊查看模態頻率和每一階模態的動態演示。

（4）將生成的mnf文件導入到ADAMS中進行仿真。

在ADAMS中讀入由Patran生成的NMF文件，查看模態頻率是否一致，在ADAMS中給柔性體設置初始速度，仿真查看柔性體動力學特性。

3 Patran/nastran、Adams聯合仿真實例

下面用Patran/nastran和Adams聯合分析太陽能電池陣等效模擬器的動力學特性。太陽能電池陣等效模擬器所用材料的彈性模量為7..1e10Pa，泊松比為0.31，密度為2770kg/m3。

將模型導入Patran中，在patran中對模型進行有限元劃分，設置太陽能電池陣等效模擬器邊界條件即左端固定，將模型交給Nastran計算求解，查看二階模態響應。將Nastran計算生成的MNF文件導入到ADAMS中。設置一階模態初始位移 0 m，初始速0.7m/s，仿真結果顯示，帆板末端沿Z軸方向速度變化趨勢趨于零。

4 結論

本文將ADAMS和PATRAN/NASTRAN結合應用，克服了兩者單獨使用時的缺陷，完成了對柔性體進行動力學和運動學分析。

參考文獻

[1]李增剛.ADAMS入門詳解與實例[M].北京：國防工業出版社，2006.

[2]解本銘，王新旭.基于ANSYS與ADAMS的并聯機構剛柔耦合體動力學分析[J].機電產品開發與創新，2012，25（1）：103-105.

[3]龍凱，賈長治，李寶峰.Patran2010與Nastran2010有限元分析[M].北京：機械工業出版社，2011.

作者簡介

田路路（1990-），女，安徽省人。碩士學位。現為上海航天控制技術研究所助理工程師。研究方向為動力學建模。

作者單位

上海航天控制技術研究所 上海市 200233