李文超， 錢 煒， 閆 強， 李宇昊， 鐘柳春

（1.上海理工大學機械工程學院，上海 200093；2.上海交通大學機械與動力工程學院，上海 200240；3.上海特種數(shù)控裝備及工藝工程技術研究中心，上海 201111）

基于DEFORM-3D的干涉配合鉚接仿真研究

李文超1， 錢 煒1， 閆 強1， 李宇昊2， 鐘柳春3

（1.上海理工大學機械工程學院，上海 200093；2.上海交通大學機械與動力工程學院，上海 200240；3.上海特種數(shù)控裝備及工藝工程技術研究中心，上海 201111）

通過借助有限元仿真軟件DEFORM- 3D建立平錐頭鉚釘?shù)某炼疹^型干涉配合鉚接模型，模擬單個鉚釘?shù)母缮媾浜香T接動態(tài)過程，測量在不同釘孔工藝參數(shù)下沿釘孔軸向方向上各測量點的徑向位移，計算出鉚釘與釘孔之間的相對干涉量.通過分析干涉量對鉚接件疲勞壽命性能的影響，得到合理的干涉配合鉚接釘孔工藝參數(shù)，并得到直徑為5 mm的平錐頭鉚釘干涉配合所需要的精確鉚接力，為實際的工藝試驗驗證和生產(chǎn)工作提供可靠的理論數(shù)據(jù).

干涉配合鉚接；孔徑；DEFORM- 3D；數(shù)值模擬

鉆鉚連接是飛機薄壁件連接的主要方式，鉚接的質(zhì)量直接影響飛機的壽命和可靠性.一般情況下，提升飛機的壽命和可靠性有應用新材料、改變結構和改善工藝三種方式，短時間內(nèi)飛機的材料和結構不會有很大的改變.因此，采用先進的生產(chǎn)工藝來提升飛機的壽命和可靠性是現(xiàn)實有效的途徑［1］.干涉配合鉚接是一種新型的鉚接工藝，其實質(zhì)是在鉚接后釘桿與釘孔之間存在一定的干涉量，該干涉量可以成倍的提高結構的疲勞壽命［2-4］，并在鉚接點處獲得良好的密封性.該技術廣泛地應用在飛行器結構部件的連接上，尤其是對疲勞壽命和氣密性要求高的部位，比如飛機增壓艙、機翼整體油箱等部位構件的連接.

通常情況下，在產(chǎn)品生產(chǎn)前需要做大量的鉚接工藝試片，通過橫銑法逐層切開鉚接試片并測量鉚接后的釘桿或者釘孔的直徑，再計算出測量部位的干涉量.做工藝試片前如果有可靠的工藝參數(shù)數(shù)據(jù)做指導，就能大大地縮短確定工藝參數(shù)的時間，降低生產(chǎn)成本.隨著計算機仿真技術的快速發(fā)展和廣泛應用，計算機科學和有限元技術在航空航天構件中結構材料及成形分析上的應用也日益成熟［5］.應用數(shù)值模擬的方法可為實際的生產(chǎn)工藝制定、產(chǎn)品的設計與優(yōu)化提供科學依據(jù).DEFORM- 3D是一款用于模擬金屬塑性成型的有限元分析軟件，應用DEFORM- 3D模擬鉚釘?shù)母缮媾浜香T接過程，分析釘孔直徑對干涉配合鉚接質(zhì)量的影響，確定合理的釘孔工藝參數(shù)，為干涉配合鉚接工藝參數(shù)的制定進行指導是現(xiàn)實有效的方法.

1 干涉配合鉚接過程仿真

1.1 DEFORM- 3D仿真前處理

仿真前處理過程主要包括幾何模型的創(chuàng)建、材料定義、網(wǎng)格劃分、摩擦與邊界條件確定以及模擬控制條件設定等.

1．1．1 幾何模型建立

鉚釘為軸對稱結構，鉚接件相對于鉚釘軸也是完全對稱的，所以在有限元分析時只取其1／4模型進行分析.上夾層需锪鉆82°／30°雙錐頭窩，鉚釘與夾層材料的幾何尺寸見表1.應用SolidWorks軟件創(chuàng)建仿真模型，把創(chuàng)建好的模型以.STL接口導入到DEFORM軟件中.定義鉚釘為塑性體，夾層為彈性體，因上下鉚模剛性要比鉚接件大的多，故視其為剛體.

1．1．2 材料屬性定義

在飛機制造業(yè)中常用的鉚釘與結構件的材料分別為2117- T4和2024- T3鋁合金，在DEFORM- 3D中添加這兩種材料，分別命名為2117- T4和2024-T3.設置材料的詳細性能參數(shù)見表2，鉚釘和鋁合金板件在塑性變形時的應力應變曲線如圖1所示［6］，并設定模擬溫度為20℃.

表1 鉚釘和夾層材料的幾何尺寸Tab.1 Geometrial dimensirns of rivet and sheets_

表2 鉚釘與板件的性能參數(shù)Tab.2 Properties of rivet and sheets

圖1 鉚釘和夾層材料塑性變形應力應變曲線Fig.1 Stress-strain curve of rivet and sheets in the stage of plastic deformation

1．1．3 網(wǎng)格劃分

鉚釘?shù)你T接過程可以看作是釘桿上的金屬沿鉚釘軸向和徑向流動的過程，隨著鉚釘?shù)溺叴忠约芭c夾層材料的接觸，材料的初始網(wǎng)格會產(chǎn)生畸變、退化，網(wǎng)格的嚴重畸變會導致求解精度的降低或者計算的不收斂，因此在有限元仿真過程中必須采用自適應網(wǎng)格重劃分技術.DEFORM- 3D具有強大的網(wǎng)格自動劃分與重劃分功能，采用完全四面體網(wǎng)格劃分中的絕對網(wǎng)格劃分形式對鉚接模型中的變形體進行網(wǎng)格劃分.而且，為提高仿真的精度和運算速度，對鉚釘和釘孔接觸區(qū)域進行局部更細的網(wǎng)格劃分，最終劃分好網(wǎng)格的幾何模型如圖2所示.

圖2 幾何模型和網(wǎng)格化分Fig.2 Geometric models and the meshing

1．1．4 摩擦與邊界條件設定添加各模型之間在鉚接過程中存在的接觸關系.鉚釘被鐓粗，因此選擇仿真的摩擦類型為剪切摩擦，根據(jù)軟件推薦的摩擦系數(shù)，設定鉚模與鋁合金模型之間的摩擦系數(shù)為0.12，鋁合金模型之間的摩擦系數(shù)為0.4，點擊全部生成按鈕，生成各模型間的接觸關系.由于采用1／4模型進行仿真模擬，在模型兩對稱面的法向上沒有金屬的流動，因此需要對鉚釘和夾層添加邊界條件.

1．1．5 仿真控制條件設定根據(jù)靜壓鉚接過程中壓鉚速度的實際應用效果，綜合考慮壓鉚效率，常用的壓鉚速度為3.0～

5.0 mm／s.因為不考慮速度對干涉量的影響，取上鉚模的壓鉚速度為3 mm／s.網(wǎng)格劃分最小單元尺寸約為0.3 mm，為了仿真收斂性及節(jié)約仿真時間，設定步長為0.1 mm，仿真步數(shù)為50步，根據(jù)體積不變原理計算出上鉚模的壓鉚位移作為仿真程序停止的條件.條件設定后檢查并生成DB文件，返回主頁面進行仿真.

1.2 DEFORM- 3D仿真后處理

干涉配合鉚接過程也是鐓粗的釘桿使釘孔脹大的過程，干涉量就是釘孔被脹大的量.通過軟件點追蹤功能可測得鉚接后各測量點Pi（i=1，2，3，4）的徑向位移Si，即可計算出該部位的干涉量，測量點的位置如圖3所示.干涉量的表示又分為兩種，一種是絕對干涉量，如式（2）所示；另一種為相對干涉量，如式（3）所示.

式中，D0為鉚接前孔徑；D為鉚接后孔徑或釘桿直徑.

圖3 干涉量測量點位置Fig.3 Measuring locations of relative interference

2 仿真結果及分析

根據(jù)航天行業(yè)標準《鉚接通用技術條件》，直徑為5 mm的實心鉚釘孔徑范圍為5.10～5.22 mm.為獲得一定的干涉量，干涉配合鉚接的釘孔直徑要比普通鉚接的釘孔直徑小，而且在干涉配合鉚接工藝過程中需要鉸孔，鉸孔直徑公差可控制在0.01 mm之內(nèi)，所以設置孔徑的尺寸分別為5.01，5.02，5.03，5.04，5.05，5.06，5.07，5.08，5.09，5.10，5.11，5.12，5.13和5.14 mm，對各孔徑參數(shù)下的模型進行仿真模擬.為了分析孔徑對干涉配合鉚接質(zhì)量的影響，除孔徑參數(shù)外其它模擬參數(shù)均相同，圖4為不同釘孔直徑下各測量點的相對干涉量.

圖4 不同孔徑下各測量點的相對干涉量Fig.4 Relative interference amount at each measuring points with different apertures

由仿真結果可以看出，上夾層兩測量點P1與P2的相對干涉量隨著孔徑尺寸的增大而減??；下夾層兩測量點P3與P4處的相對干涉量在5.01～5.04 mm的孔徑范圍內(nèi)緩慢增加，隨后相對干涉量隨著孔徑尺寸的增大而逐漸減小.這是因為在上鉚模與釘桿接觸后釘桿逐漸被鐓粗，墩頭窩和釘孔的過渡區(qū)域首先與被鐓粗的釘桿相接觸，由于釘孔間隙非常小，上夾層產(chǎn)生較大的干涉量阻礙了釘桿向下的金屬流動.隨著釘孔直徑的增大（5.01～5.04 mm之間），在釘桿被鐓粗的過程中金屬向下流動的阻力變小，下夾層所得到的干涉量緩慢增加.隨著孔徑的進一步增加，釘孔間隙加大，下夾層獲得的干涉量又逐步減小.

上夾層獲得的干涉量明顯的比下夾層所獲得的干涉量大，而且越靠近鉚釘頭的部位獲得的干涉量就越小.但P1與P2、P3與P4之間的相對干涉量在數(shù)值上比較接近，說明沿鉚釘軸線指向鉚釘頭的方向上各夾層之間獲得干涉量的難易程度不同，在鉚釘被鐓粗金屬向下流動的過程中上夾層更容易獲得較大的干涉量，下夾層相對較難得到大的干涉量，整個鉚釘軸線上獲得的干涉量不夠均勻.要獲得較為均勻的干涉量，需要金屬從釘桿的兩端向中間流動，采用無頭鉚釘即可實現(xiàn)，該方法需要上下鉚模同時壓鉚，對設備要求較高.

根據(jù)前蘇聯(lián)的研究結果，鋁合金材料的相對干涉量應在1.5%～3.0%之間［7］，該范圍內(nèi)的干涉量有利于鋁合金結構件疲勞壽命的提高.由仿真的結果可知，當孔徑小于5.04 mm時，上夾層的相對干涉量大于3%；當孔徑大于5.12 mm時，下夾層的相對干涉量小于1.5%.因此合理的孔徑范圍為5.04～5.12 mm.由仿真數(shù)據(jù)可以得出，當孔徑為5.08～5.10 mm時，P1與P2，P3與P4的干涉量相差最小，推薦釘孔尺寸范圍為5.08～5.10 mm.

由仿真過程中的載荷行程曲線圖5可知，在填滿墩頭窩之前鉚接力緩慢的增加，待填滿墩頭窩之后鉚接力快速增大，最后反求出鉚接力大小約為21.7 kN.

3 結 論

應用有限元分析軟件DEFORM- 3D模擬了平錐頭鉚釘?shù)母缮媾浜香T接過程，提高釘孔工藝參數(shù)對平錐頭鉚釘干涉配合鉚接質(zhì)量影響的認識，增加鉚接成型的預見性.通過仿真分析得到合理的釘孔工藝參數(shù)，推薦材料為2117- T4，直徑為Φ5 mm的平錐頭鉚釘干涉配合鉚接的孔徑范圍為5.08～5.10 mm，鉚接變形需要的壓鉚力為21.7 KN.應用有限元分析軟件研究工藝參數(shù)對鉚接質(zhì)量的影響，獲得基本能反映實際情況的試驗數(shù)據(jù)，避免一些不合理的參數(shù)設置，從而可提高生產(chǎn)效率，降低成本.

圖5 載荷-行程曲線Fig.5 Riveting power-stroke curve

［1］ 曹增強，馮東格.理想干涉配合鉚接［J］.航空科學技術，2012，23（4）：84-86.

［2］ 郝小忠，王珉，陳文亮.自動鉆鉚最佳干涉量試驗研究［J］.航空制造技術，2012（12）：77-79.

［3］ Mu W，Li Y，Zhang K，et al.An effective method of studying interference-fit riveting for 2117-T4 aluminum slug rivet［C］／／2010 International Conference on Computer and Communication Technologies in Agriculture Engineering，Chengdu：IEEE，2010：303-307.

［4］ Cao Z，Cardew-Hall M.Interference-fit riveting technique in fiber composite laminates［J］.Aerospace Science and Technology，2006，10（4）：327-330.

［5］ 徐學春，胡廣洪，董萬鵬.Deform在航天航空工業(yè)中的應用［J］.精密成形工程，2012，4（6）：53-55.

［6］ 宋丹龍，李原，駱彬，等.CFRP／Al復合構件無頭鉚釘壓鉚力建模與仿真分析［J］.西北工業(yè)大學學報，2012，30（4）：558-564.

［7］ 張昕.復合材料干涉配合連接有限元分析［D］.西安：西安電子科技大學，2011.

（編輯：董 偉）

Simulation of Interference-fit Riveting by Using the Software DEFORM-3D

LIWen-chao1， QIANWei1， YANQiang1， LIYu-hao2， ZHONGLiu-chun3
（1.School of Mechanical Engineering，University of Shanghai for Science and Technology，Shanghai 200093，China；2.School of Mechanical Engineering，Shanghai Jiaotong University，Shanghai 200240，China；3.Shanghai Special Numerical Control Equipment and Technology Engineering Research Center，Shanghai 201111，China）

By using the finite element analysis software DEFORM- 3D，the models for flat cone head rivets were established under the technological condition of countersunk head interference-fit riveting.The dynamic process of interference-fit riveting of a single rivet was simulated.The radial displacement at a series of measuring points along the axial direction of nail holes under different technological parameters were measured，and the relative amount of interference between the rivets and rivet holes were calculated.The effect of interference amount on the fatigue life of structure was analysed.The reasonable parameters for aperture processing were achieved and the necessang riveting force for 5 mm diameter flat cone head rivets was accurately estimated.The study provides reliable dates for actual process test and production.

interference-fit；aperture；DEFORM- 3D；numerical simulation

TH 131.1；TP 391.9文獻標示碼：A

1007-6735（2014）05-0479-04

10.13255／j.cnki.jusst.2014.05.013

2013-08-28

上海市科技創(chuàng)新行動計劃資助項目（12DZ2250500）

李文超（1987-），男，碩士研究生.研究方向：航空航天自動化裝配與制造.E-mail：liwenchao100＠126.com

錢 煒（1964-），男，副教授.研究方向：機械結構設計與理論.E-mail：qyc0510＠163.com