[发明专利]一种飞机薄壁件自动钻铆多姿态变形建模方法

摘要

公开了一种飞机薄壁件自动钻铆多姿态变形建模方法。该方法包括步骤将飞机薄壁件和托板之间的定位夹紧关系简化为连续梁受多个支承座支承的模型；计算各托板的支持力，在所有中间托板位置将薄壁件截开并将托板支承简化为简支支承；求解出各个托板的弯矩后，进一步利用两个托板的受力平衡模型，计算自动钻铆过程中承载飞机薄壁件的各个托板的支持力；将薄壁件按照托板位置划分成相应的小区域，针对每个小区域求解薄壁件变形，其中，每个小区域薄壁件受力变形模型简化为托板位置薄壁件简支支承，薄壁件表面仅受均布载荷的变形计算；最后计算出薄壁件中面沿方向的变形量。

主权项

一种飞机薄壁件自动钻铆多姿态变形建模方法，其包括如下步骤：(1)将飞机薄壁件和托板之间的定位夹紧关系简化为连续梁受多个支承座支承的模型，其中，简化后各跨的薄壁件所受重力用均布力qw表示qw＝Gi/Li，其中，Gi表示第i跨的薄壁件重力，Li表示第i跨的长度；(2)计算各托板的支持力，在所有中间托板位置将薄壁件截开并将托板支承简化为简支支承，其中，利用下述公式1‑1求解各个托板位置的弯矩：式中：li表示i‑1到i托板间的跨距m；Mi表示i托板截面上的弯矩KN·m；ai表示ωi的形心ci到i支承的距离m；ωi表示i‑1到i支承间载荷弯矩图的面积m2；bi+1表示ωi+1的形心ci+1到i+1支承的距离m；Ei表示i‑1到i托板间薄壁件材料的弹性模量GPa；Ii表示i‑1到i托板间薄壁件简化为连续梁的截面惯性矩；当薄壁件的材料相同时，各个跨距段薄壁件的弹性模量Ei相同，则公式1‑1变形为：式中，li、Ii、ωi、ai、bi+1(i＝0，1，2，3，…，n)为已知条件，n+1个托板中有n‑1个中间托板，列出n‑1个方程，加上方程M0＝0和Mn＝0，求解出n‑1个中间托板所承受的弯矩Mi；(3)求解出各个托板的弯矩后，进一步利用两个托板的受力平衡模型，计算自动钻铆过程中承载飞机薄壁件的各个托板的支持力：(最前端托板i＝0)   (公式1‑3)(中间托板i＝1，2，…，n‑1)  (公式1‑4)(最后端托板i＝n)   (公式1‑5)(4)将薄壁件按照托板位置划分成相应的小区域，针对每个小区域求解薄壁件变形，其中，每个小区域薄壁件受力变形模型简化为托板位置薄壁件简支支承，薄壁件表面仅受均布载荷的变形计算，其中：薄壁件简支支承边界条件满足：T11＝0、S＝γ12Eh/2(1+μ)＝0、ω＝0、M11＝0、M12＝M21＝0，其中：T11表示中面内力；其中，S表示剪切力；γ12表示平面剪切应变；E表示弹性模量；h表示板厚；μ表示泊松比；ω表示法向位移；M11表示壁板受到的弯矩；M12和M21表示壁板受到的扭矩；薄壁件自由边界条件满足：μ1＝0，μ2＝0，ω＝0，M11＝M011＝0，φ2＝0，其中：μ1表示沿x方向位移；μ2表示沿y方向位移；M011表示表示壁板各个区域左边支撑处的弯矩；φ2表示沿y轴的转角；建立满足以上边界条件的薄壁件变形计算，设其中：F和f为求解过程假设的函数；L表示托板均匀排布时，相邻两个托板间的距离；将公式1‑6和公式1‑7代入薄壁件的静力方程并将所受载荷展开成傅里叶级数得到公式1‑8和公式1‑9：其中：R表示壁板半径；μ表示壁板材料泊松比；k0表示h表示壁板厚度；k表示和壁板尺寸及材料相关的参数，等于K表示和壁板厚度及壁板材料相关的参数，等于α表示自动钻铆机工作过程中绕全局坐标系X轴转动的角度；β表示自动钻铆机工作过程中绕全局坐标系Y轴转动的角度；公式1‑8的特解为：求解公式1‑8的齐次解，设：代入得：上式简化为：式中：根据方程求出特征值λn(n＝1，2，…8)，则方程的齐次解为：Fnh=A1eλ1y+A2eλ2y+A3eλ3y+A4eλ4y+A5eλ5y+A6eλ6y+A7eλ7y+A8eλ8y]]>全解为：Fn=Fnh+2qn[1-cos(nπ)]Kk(1-μ)nπ[(nπL)8-1-μ2kR2(nπL)6+1-μ2kR2k0(nπL)4]]]>对于公式1‑9，设：将上式代入公式1‑9得到特征方程：求得特征根为：λm=±(nπL)2+2k01-μ(m=9,10)]]>于是：由Fn、fn可以得出F、f，又由公式：其中：表示微分算子计算出F1、F2，同时，M11＝D(κ1+uκ2)；M22＝D(κ2+uκ1)其中：τ12表示沿壳中面1方向的扭曲变形；τ21表示沿壳中面2方向的扭曲变形；D表示和壁板材料参数相关的常数；κ1表示沿壳中面1方向的弯曲变形；κ2表示沿壳中面2方向的弯曲变形φ1表示沿x轴的转角；φ2表示沿y轴的转角；薄壁件中面位移u1、u2、ω、φ1、φ2：将已经计算出的F1、F2、f分别代入公式1‑17到公式1‑21并结合边界条件计算出薄壁件中面沿方向的变形量，式中：