基于ADAMS的推彈滑座與阻鐵在不同傾角下的碰撞動力學分析

申培剛，戴勁松，王茂森

（南京理工大學機械工程學院，江蘇南京 210094）

0 引言

扣機是火炮自動機的關鍵組成部件，它主要用來控制自動機射擊的啟停。本文采用的是機械式扣機方案，其通常安裝在炮箱下部有限的空間內。由于火炮發射機構結構復雜［1］，零部件受力條件惡劣，設計扣機本體結構時必須準確地計算和評估推彈滑座與阻鐵碰撞后的沖擊力。研究推彈滑座與阻鐵在不同傾角下的動力學問題，在提高零部件剛強度及運動的可靠性方面提供了依據，同時據各零部件運動規律的唯一性，對其大致外廓尺寸設計也有一定的參考。通過虛擬樣機對推彈滑座與阻鐵在不同角度下的碰撞性能進行仿真分析，優化分析出了最佳角度下的模型結構，對火炮自動機的作戰性能及可靠性方面提供了一定的參考。

1 動力學建模

1.1 基本假設及模型簡化

據扣機機構實際的運動過程作如下假設及簡化:

1）忽略該機構液壓系統的建模，但以相應的力作用于該機構的機械部件;

2）忽略對仿真結果影響不大的零部件如銷軸及彈簧等;

3）對各構件做剛性處理，忽略構件的彈性變形;

4）對機構中的零部件進行簡化處理，忽略相關的圓角及倒角。

1.2 模型簡述

利用三維建模軟件建立火炮扣機的三維實體模型后，將其轉化為相應的格式導入ADAMS軟件中建立有效地動力學仿真模型。由于火炮射擊時間較短，推彈滑座與阻鐵在撞擊的過程中持續時間非常短且相互間的作用是極其復雜的，同時力又是巨大的，通常做如下假設［2］:

1）在碰撞過程中，與碰撞力相比較，常規力（如重力、彈性力、物有引力等）遠小于碰撞力的沖量，因此常規力的沖量可以忽略不計。

2）在碰撞時，由于時間非常短促，物體的位移可忽略不計。

為了便于研究自動機各構件的撞擊，通常用反映機構的主要特性，而忽略一些次要特性的結構簡圖來代替實際的機構圖?？蹤C部件中推彈滑座與阻鐵的碰撞簡圖如圖1所示。

在機械系統動力學仿真軟件中設定全局坐標系，該模型簡化為14個剛體，該模型拓撲關系結構如圖2所示。其中滾輪、滾輪座及擊發軸通過（h1）固定鉸連接;傳動杠桿、轉動撥叉及回位杠桿以（h2）旋轉鉸連接;推彈滑座與傳動杠桿間以（h3）滑移鉸連接;阻鐵撥叉與阻鐵間以以（h4）旋轉鉸連接;彈簧導桿與阻鐵間以（h5）滑移鉸連接;解脫器、滑座撥叉通過（h6）旋轉鉸;所有零部件間定義實體與實體接觸來傳遞各受力。

圖1 復進時推彈滑座與阻鐵碰撞受力示意圖

圖2 拓撲關系結構

1.3 多體系統動力學方程

多剛體系統動力學是基于經典力學理論中最簡單的情況即自由質點和一般簡單的少數多個剛體，利用虛擬樣機中帶拉格朗日乘子的第一類拉格朗日方程建立多體系統的動力學方程［3］。它以每個零部件質心在慣性坐標中的笛卡爾坐標與歐拉角為廣義坐標，表達式如下:

由零件質心的參考系與地面坐標系的變換矩陣經轉化所得動力學方程為:

其中:T—系統的動能;

qi—為廣義坐標;

Qi—為廣義力;

動力學模型建立完成后由多體系動力學模型得到動力學數學模型，從而帶入求解器進行求解，其建模與求解流程如圖3所示。

2 ADAMS中碰撞力的定義

圖3 動力學建模與求解流程

ADAMS/Solver中采用兩種方法計算碰撞力:一種是回歸法（Restitution），通過定義恢復因數來計算碰撞力;另一種是沖擊函數法［4］（Impact），通過定義剛度系數和阻尼系數來計算碰撞力。由于回歸法的參數設置比較難，且適用于恢復因數已知的情況，因此多用后者來計算碰撞力。沖擊函數法的計算速度較快，使用的參數（如剛度系數及阻尼系數等）能夠更好地控制碰撞力，在仿真中能得到滿意的效果，因此采用沖擊函數法來計算推彈滑座與阻鐵間的碰撞力，其一般表達式為:

式中:k為剛度系數，e為指數，cmax為阻尼系數，d為切入深度，其決定阻尼力的大小，x為沖擊函數中的距離變量，x0兩物體碰撞的初始距離，dx/dt為距離隨時間的變化率（即速度）。為了防止碰撞過程中阻尼力的不連續，式中采用了step（）子函數可參照ADAMS幫助文件［5］。當x≥x0時表示物體不發生接觸，此時碰撞力為零，當x＜x0時，表示兩物體發生碰撞，其碰撞力大小與式（3）中的參數有關。

3 推彈滑座與阻鐵碰撞的動力學仿真

3.1 仿真參數設置

采用三維建模軟件來完成對火炮扣機三維實體模型的建立，把模型導入ADAMS后對各個零件定義相關的約束、接觸和施加外力［6］，在接觸命令中涉及到剛度、阻尼等各系數查閱碰撞參數表，其接觸參數的選取如表1所示。

表1 推彈滑座與阻鐵碰撞參數設置

為了提高仿真的速度，縮短仿真時間，設定仿真的時間為0.1 s，仿真的步數為100。

3.2 動力學仿真及分析

上述設置完成后，開始進行仿真。隨著傾角φ的增大，為保證推彈滑座與阻鐵能夠在不同傾角φ下很好的解脫，在滾輪座上施加的外力也會隨之增加?？紤]到應在施加外力相同的情況下比較推彈滑座與阻鐵的碰撞力，由于在大傾角時施加的初始外力會對小傾角的仿真曲線有很大，鑒于其對仿真結果準確性的影響，在φ≤50°時，分析的是在推彈滑座與阻鐵后坐到位后，在液壓系統不工作阻鐵已抬起時推彈滑座的水平方向的位移大小，如圖4—圖6所示。

從圖4—圖6中知，在液壓系統不工作時推彈滑座在外力的驅動下復進，在最大位移處時阻鐵的抬起阻止了推彈滑座的繼續移動，經幾次碰撞后停下。推彈滑座在最低點處位移大小依次為 26.150 8 mm，25.805 mm 及25.588 8 mm。隨著傾角的逐漸增大，推彈滑座在后坐到位后運動的最大位移是逐漸減小的。在φ=47°時為26.150 8 mm，當 φ =50°時為25.588 8 mm，這是因為在液壓系統不工作下阻鐵已抬起，阻鐵能夠阻擋住推彈滑座的運動。在停射時出于安全考慮推彈滑座的位移越小越好。

在φ≥51°時分析推彈滑座與阻鐵在解脫模型下的碰撞。由于該碰撞為斜碰撞，即在x（水平）方向及y（豎直）方向推彈滑座與阻鐵的碰撞力都存在，因此每個傾角對應的碰撞力仿真曲線有兩條，如圖7—圖10所示。

從圖7—圖10的曲線看出水平、豎直方向的碰撞曲線在大約2.5 ms時發生突變，經過一段時間后逐漸減小，出現這種現象的原因是推彈滑座與阻鐵兩物體在瞬間突然接觸發生撞擊，物體的動量在極短的時間內發生明顯變化造成的。隨著角度φ的增大，阻鐵與推彈滑座在水平方向和豎直方向的最大碰撞力呈現越來越大的趨勢，且隨著傾角φ的增大碰撞作用的時間也逐漸增長。在φ=51°時，水平方向的最大碰撞力為4 569.553 7 N，豎直方向的最大碰撞力5 810.763 7 N;當φ=54°時，水平方向的最大碰撞力為13 475.093 1 N，豎直方向的最大碰撞力為16651.1337 N，水平、豎直方向的力增大了接近三倍，且φ=51°時，其碰撞力較小且作用時間相對其它來說也較短。由圖4—圖6知φ=51°時推彈滑座的最大位移比φ=50°時的要小，因此在φ=51°時構件的剛強度及可靠性方面結構最優。

4 結語

本文利用三維建模軟件建立了機械式的火炮扣機模型，利用機械系統動力學仿真分析軟件ADAMS對火炮扣機的推彈滑座與阻鐵構件進行了碰撞動力學分析，通過在小傾角下比較推彈滑座的位移的大小，在大傾角時比較推彈滑座與阻鐵的碰撞力的大小，優化分析出了最佳角度下的模型結構，在提高零部件剛強度及運動的可靠性方面為其提供了依據。

［1］郭保全，彭峰生，潘玉田，等.扣機電磁鐵優化設計［J］.火炮發射與控制學報，2005（3）:39-41.

［2］赫雄.ADAMS動力學仿真算法及參數設置分析［J］.傳動技術，2005，19（3）:27-30.

［3］孫世賢，黃圳圭，唐乾剛，等.理論力學教程［M］.長沙:國防科技大學出版社，2009.

［4］范成建，雄光明，周明飛.虛擬樣機軟件MSC.ADAMS應用與提高［M］.北京:機械工業出版社，2006.

［5］陳德民，槐創峰，張克濤，等.精通ADAMS2005/2007虛擬樣機技術［M］.北京:化學工業出版社，2010.

［6］魏衍俠.基于ADAMS的手拋式機器人碰撞動力學分析［J］.機械工程與自動化，2011（2）:82-84.