摆动柔顺式吊钩结构参数的多目标优化

doi:10.16183/j.cnki.jsjtu.2020.283

摆动柔顺式吊钩结构参数的多目标优化

葛世程, 郭着雨, 梁熙, 莫宗来, 李军^,

南京理工大学机械工程学院, 南京 210094

Multi-Objective Optimization for Structural Parameters of Swing-Compliant Hook

GE Shicheng, GUO Zhuoyu, LIANG Xi, MO Zonglai, LI Jun^,

School of Mechanical Engineering, Nanjing University of Science and Technology, Nanjing 210094, China

通讯作者: 李军,男,教授,博士生导师,电话(Tel.): 18936030573;E-mail:leejun_nust@163.com.

责任编辑: 石易文

收稿日期: 2020-09-8

Received: 2020-09-8

作者简介 About authors

葛世程(1996-),男,广西壮族自治区南宁市人,硕士生,主要从事火箭炮发射箱吊装系统设计及结构优化研究.

摘要

针对目前野战火箭炮发射箱装填过程的对准问题,提出一种新型末端执行机构——摆动柔顺式吊钩,并对其结构参数进行优化.首先,基于节点位移法建立描述摆动柔顺式吊钩性能的理论模型,并通过MATLAB分析摆动柔顺式吊钩的静态位移、静态应力和摆动曲线,验证模型工作的合理性.同时,采用实验设计方法筛选影响摆动柔顺式吊钩性能的主要结构参数,并建立表征摆动柔顺式吊钩综合性能的响应面模型.优化结果表明,当柔性机构长度为90 mm,其安装高度为23 mm,吊钩末端高度为110 mm时,摆动柔顺式吊钩在对接、起吊、转运和落装时的性能最优.

关键词： 发射箱; 装填系统; 柔性机构; 节点位移法; 实验设计方法

Abstract

To solve the problems of alignment for launch container loading in field-artillery rocket, a novel end-effector, i.e., a swing-compliant hook, is proposed, whose structural parameters are optimized. First, the theoretical model describing the performance of the swing-compliant hook is established based on the node displacement method. The static displacement, static stress, and swing curve of the swing-compliant hook are analyzed by MATLAB, which verifies the rationality of the model. Then, the main structure parameters on the performance of the swing-compliant hook are obtained by using the experimental design method. Additionally, a response surface model characterizing the comprehensive performance of the swing-compliant hook is established. The optimization results show that when the length and the installation height of the compliant mechanism are 90 mm and 23 mm, and the end height of the lifting hook is 110 mm, the swing-compliant hook has an excellent performance in docking, lifting, transferring, and locating.

Keywords： launch container; loading system; compliant mechanism; node displacement method; experimental design method

PDF (1931KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

葛世程, 郭着雨, 梁熙, 莫宗来, 李军. 摆动柔顺式吊钩结构参数的多目标优化[J]. 上海交通大学学报（自然版）, 2021, 55(11): 1467-1475 doi:10.16183/j.cnki.jsjtu.2020.283

GE Shicheng, GUO Zhuoyu, LIANG Xi, MO Zonglai, LI Jun. Multi-Objective Optimization for Structural Parameters of Swing-Compliant Hook[J]. Journal of shanghai Jiaotong University, 2021, 55(11): 1467-1475 doi:10.16183/j.cnki.jsjtu.2020.283

野战火箭炮的再装填过程需要将发射箱在装填车、发射车与地面之间进行多次对接、起吊、转运和落装,这些操作通常依靠装填车上的装填系统完成^[1].为保证射击密集度,火箭炮利用发射箱底部定位孔与发射车的定位销进行约束并锁定发射箱,对准时需要保持较高的位置精度.而现有绳索吊装技术由于绳索的变形特性常常造成发射箱的转运以及对准难以控制^[2,3].为根本解决绳索吊装问题,刚性装填技术利用吸盘、电磁铁或抓取机构等较复杂的末端执行器起吊工件,能够较好地控制执行器与发射箱的对接以及箱体的转运^[4,5,6].但在对准时需大幅降低运动速度,同时难以利用导向结构降低对接、起吊和落装过程中执行器所需的工作精度,如S400防空导弹的装填.因此,为了实现高效装填,需要末端执行机构有柔性变形能力以补偿误差带来的执行器偏移^[7].

在航空航天、船舶和自动化系统中,柔性对接机构已经获得了较为广泛的应用.文献[8]设计的一种2自由度柔顺移动并联机构,在有较大的径向平移偏差和角度偏差情况下仍能完成柔顺对接.文献[9]针对多组销孔对接问题,提出一种力柔顺控制算法降低销孔间的接触力.文献[10]设计了一种柔顺对接与撤收机构,提高机器人与火箭定位、对接、加注和撤收的效率.由于野外地形复杂,发射箱的对准要求执行机构有更大的误差容许范围,同时还能承载发射箱的质量,然而以上要求对接机构均难以满足.在中小型航天器对接机构中,锥-杆式对接机构因其结构简单、误差容许范围大获得了较好的应用.文献[11]设计了一种卫星柔性对接机构,能在20 mm横向对接容差及±5° 角度容差范围内实现可靠对接.

由于航天器工作环境的特殊性,对于锥-杆式对接机构的研究以分析碰撞过程为主,鲜有对对接机构在重力作用下的工作特性进行分析^[12].此外,航天器允许的对接时间较长、对接次数较少,因而对于对接机构的运动特性和使用寿命的研究也较少.综上所述,有必要研究锥-杆式对接机构的应用性能,并合理设计结构尺寸以满足综合性能的要求.

本文针对野战火箭炮发射箱的装填要求,设计了一种适用于发射箱柔性装填系统的摆动柔顺式吊钩末端执行机构,较好地控制发射箱的转运和对准,同时降低执行器所需的精度.为获得摆动柔顺式吊钩的最佳结构尺寸,首先基于节点位移法建立摆动柔顺式吊钩的参数化理论模型,其次对摆动柔顺式吊钩的结构参数进行筛选,最后采用主要参数建立响应面模型并对其进行优化.优化结果可为进一步的吊钩结构优化和工程实际提供理论依据.

1 摆动柔顺式吊钩工作原理

摆动柔顺式吊钩是发射箱柔性装填系统的主要功能部件,如图1所示.其固定于吊装平台,工作时与发射箱上的横梁对接,实现对发射箱的装填.具体工作过程如下,吊臂及吊装平台根据指令运动至发射箱上方,使吊钩与发射箱上的横梁对准,进一步做U型运动完成与横梁的对接并顺势起吊发射箱.吊臂及吊装平台将发射箱转运至发射车上方,使发射箱底部的定位销孔与发射车上的定位销对准,同时向下运动完成落装.落装完成后,吊钩作U型运动完成脱钩.

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 柔性装填系统示意图

Fig.1 Schematic diagram of compliant loading system

摆动柔顺式吊钩几何模型如图2所示.摆动柔顺式吊钩为被动柔顺式机构,可以自然顺从以补偿吊钩与目标位置的偏差.其中,柔性机构是实现装填系统柔性功能的部件,通过防松螺栓分别与固定支座和转动支座联接,用于限制吊钩摆动.固定支座、转动支座和摆臂互相配合形成虎克铰链,使吊钩可以在空间进行2自由度摆动.

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 摆动柔顺式吊钩几何模型

Fig.2 Geometric model of swing-compliant hook

摆动柔顺式吊钩的工作原理如图3所示.空载状态下(见图3(a)),在柔性机构的限制下,吊钩快速稳定且能保持与固定支座的位置关系,即摆动偏角最小,允许的最大角度偏差为2°.负载状态下(见图3(b)),吊钩在发射箱的重力作用下被动柔顺至稳定状态,使落装前摆臂轴线均与重力方向接近平行,即摆动偏角最大,允许的最大角度偏差为2°,以利于根据吊装平台位姿确定发射箱位置姿态.

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 摆动柔顺式吊钩工作原理

Fig.3 Working principle of swing-compliant hook

2 摆动柔顺式吊钩理论模型

2.1 参数及坐标系定义

在研究柔性机构限制下吊钩的摆动特性时,为方便后续的参数化建模,仅保留摆动柔顺式吊钩的约束情况,摆动柔顺式吊钩的初步简化图如图4所示.其中:l_c、h₂和h₁分别为柔性机构的长度、高度和安装高度;b₁和b₂为柔性机构边缘宽度和中部宽度;h₃为吊钩末端高度;A₁、A₂、B₁、B₂、C₁和C₂为螺栓联接点;D为接触点;E₁为吊钩等效质量点.

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 摆动柔顺式吊钩简化图

Fig.4 Simplification of swing-compliant hook

在结构参数的研究中,节点位移法将优化对象离散为若干柔性梁,以梁的联接位置为研究对象,并在给定条件下计算满足设计要求的最优结构^[13].坐标系定义及柔性机构离散如图5所示.其中:C₁₁、C₂₁、D₁和E₂依次为点C₁、C₂、D和E₁在吊钩摆动后的新位置;α和β分别为吊钩绕固定支座和转动支座的角度;O'-x₁y₁z₁为吊钩运动坐标系.由于柔性机构的具体形状尚未确定,仅将柔性机构简化成A₁D、A₂D、B₁D、B₂D、C₁D、C₂D、A₁C₁、A₂C₂、B₁C₁和B₂C₂柔性梁.

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 坐标系定义和柔性机构离散

Fig.5 Definition of coordinate system and discretization of compliant mechanism

本文关注的是柔性机构在给定变形下产生的抵抗力,而不是柔性机构的变形情况,因此将柔性梁作为一个梁单元处理.设简化后的梁单元固定点为完全约束,受力点为自由端,截面为圆形,且只考虑拉压和弯曲两种变形情况,不计自身重力.

考虑到吊钩的运动包括摆动中心(如图2虚线所示的两回转轴交点)的直线运动和吊钩的摆动,为便于分析,建立如下坐标系.

(1) 吊钩运动坐标系采用直角坐标表示吊钩的运动特性和结构,见图5中的O'-x₁y₁z₁,坐标原点位于吊钩摆动中心.以柔性机构未变形时的摆臂轴向为z₁轴正向,转动支座轴线为y₁轴,固定支座轴线为x₁轴.结合图4中的6个结构参数及偏摆角可以获得各点在O'-x₁y₁z₁中的坐标.

(2) 装填车坐标系(见图5中的O-XYZ),用于描述吊钩相对装填车的位置.为便于计算,将O'-x₁y₁z₁的初始位置与O-XYZ设在同一位置.为通过修改O-XYZ中的重力方向模拟吊装平台的倾斜状态,将O-XYZ与大地固定,且Z轴竖直向下.

由于柔性机构的性能随倾斜角度的减小而提高,根据发射车与装填车间的倾角要求,为保证设计结果的可靠性,O-XYZ中的重力方向均取(0.2,0.2,1) mm.

2.2 吊钩及柔性机构静力学模型

2.2.1 吊钩位移模块该方程用于计算吊装平台倾斜一定角度时,吊钩在柔性机构弹力和重力作用下相对于固定支座的偏角.

(1) 梁单元弹力

A₁D单元 $L_{A_{1} D}$ 抵抗力分解图如图6所示.在O'-x₁y₁z₁中,将接触点D设为受力点,D₁为变形后的受力点,螺栓联接点A₁为固定点,将梁单元变形分解为拉伸变形和弯曲变形.拉伸抵抗力F可表示为

F=ES

$\frac{Δ l}{l}$

$\frac{ES (|L_{A_{1} D_{1}}| - |L_{A_{1} D}|)}{|L_{A_{1} D}|}$

(1)

式中:Δl为伸长量;l为梁原长;E为材料弹性模量;S为梁单元的原截面面积.变形后梁单元对受力点的拉伸抵抗力方向为D₁A₁.

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 A₁D单元抵抗力分解图

Fig.6 Resistance decomposition of A₁D element

梁单元的弯曲视为由梁单元转动平面内、垂直于 $L_{A_{1} D}$ 且作用在受力点上的力推动.因此,梁单元产生的对应弯曲抵抗力F'为

F'=

$\frac{3 w_{D} EI}{l'^{3}}$

$\frac{3 EI |L_{A_{1} D_{1}}| \sin θ}{{|L_{A_{1} D_{1}}|}^{3}}$

(2)

式中:w_D为弯曲挠度;I为轴变形后梁截面对其中性轴的惯性矩;l'为梁变形后的总长度;θ为 $L_{A_{1} D_{1}}$ 与 $L_{A_{1} D}$ 的夹角.泊松比取为0.45.变形后梁单元对受力点的弯曲抵抗力方向为 $L_{A_{1} D}$ ×( $L_{A_{1} D}$ × $L_{A_{1} D_{1}}$ ).

其中,C₁D和C₂D两个梁单元仅在吊钩摆角β不等于0 rad的情况下产生弯曲抵抗力.由于C₁D和C₂D的联接点受力均会影响吊钩摆动,采用固定点与受力点轮换的方式求作用在C₁、C₂、D上的梁单元弹力.

(2) 吊钩合力矩

根据O'-x₁y₁z₁在O-XYZ中的位置获得梁单元的拉伸抵抗力和弯曲抵抗力在O-XYZ各轴方向上的投影 $F_{X_{i}}$ 、 $F_{Y_{i}}$ 和 $F_{Z_{i}}$ ,其中i=1,2,…,12为受力点序号.设在O-XYZ中,O'-x₁y₁z₁原点的坐标为O'(O'_X, O'_Y, O'_Z),相应的受力点坐标为(x_i, y_i, z_i).则弹力对于x₁轴的合力矩为

$M_{x_{1}}$ =∑

$[- F_{Y_{i}} (z_{i} - O'_{Z}) + F_{Z_{i}} (y_{i} - O'_{Y})]$

(3)

转动支座轴线C与 $L_{C_{21} C_{11}}$ 平行,又与摆臂轴线 $L_{O' D_{1}}$ 垂直,且联接点C₁₁和C₂₁的受力不计入转动支座轴的合力矩方程.因此,弹力对于转动支座轴的合力矩为

M₀=∑

$[\frac{L_{C_{21} C_{11}} \times L_{O' D_{1}}}{|L_{C_{21} C_{11}} \times L_{O' D_{1}}|} [F_{X_{i}} F_{Y_{i}} F_{Z_{i}}]]$

(h₁+h₂)

(4)

根据摆臂以及钩的质量可获得作用在E₂点的等效重力G_m,考虑到吊钩的受载情况,取摆臂线密度ρ=9.6 g/mm,同时为了扩大各因素差异,取钩的质量m=10 kg.同理,结合给定重力方向以及式(3)可获得重力对x₁轴的力矩M $″_{x_{1}}$ ,在式(4)中添加吊钩末端长度h₃可获得重力对转动支座轴的力矩M″₀.

负载工况下,满载发射箱质量为7 t,取发射箱作用在吊钩的等效重力G_M=17.5 kN.在式(3)和 (4)中增加G_M可以分析负载状态下的吊钩摆角.

(3) 吊钩平衡位置

在同一算例中,柔顺式吊钩的参数和重力方向一定,因此吊钩拥有唯一空间位姿使其关于x₁轴和转动支座轴的合力矩为0.使用黄金分割法查找预估范围内α和β的解,设空载工况下的解为α₁和β₁,负载工况下的解为α₂和β₂,求解模型如下:

$\begin{array}{l} \min : M (α, β) = |M_{x_{1}} + M ″_{x_{1}}| + |M_{0} + M ″_{0}| \\ subjectto : - 0.2 rad < α < 0.2 rad \\ - 0.2 rad < β < 0.2 rad \\ M (α, β) \geq 10 (N \cdot mm) \end{array}\}$

(5)

2.2.2 柔性机构应力模块柔性机构应力模块用于计算吊钩与固定支座成一定摆角时,柔性机构简化后各个梁单元的应力及梁单元的总应力.

(1) 单个梁单元拉应力

在吊钩位移模块基础上,获得梁单元拉应力为

σ_j=E

$\frac{Δ l}{l}$

(6)

式中:j=1,2,…,N为梁单元序号.

(2) 单个梁单元最大弯曲应力

以A₁D梁为例,最大弯曲应力为

σ'_j=

$\frac{F' |L_{A_{1} D_{1}}| d' \cos θ}{2 I}$

(7)

式中:d'为变形后的截面直径.

由于一般材料的拉伸屈服极限和压缩屈服极限大致相同,不再区分拉应力和压应力,所以单个梁单元的最大应力为 $|σ_{j}|$ +σ'_j,并以此为该梁单元应力,可获得全部梁单元总应力为

σ_sum=∑(

$|σ_{j}|$

+σ'_j)

(8)

全部梁单元应力中的最大值为

σ_max=max(

$|σ_{j}|$

+σ'_j)

(9)

2.2.3 静态检验为验证吊钩位移模块和柔性机构应力模块的正确性,在Abaqus软件中使用B33梁单元建立柔性机构并创建对应约束和载荷模拟柔性机构.柔性机构应力模块中通过设定吊钩摆角与重力方向平行实现负载工况.通过静力分析,得到的柔性机构空载位移s和负载应力σ如图7所示.

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 Abaqus静态检验

Fig.7 Static verification by Abaqus

由图7可知,动点D₁的空载位移较大,柔性机构负载应力较小.相同条件下,理论模型分析结果显示,动点D₁的空载位移为2.8 mm,柔性机构最大负载应力为5.37 MPa.分析结果差距较小,认为吊钩位移模块和柔性机构应力模块可以描述摆动柔顺式吊钩性能.

2.3 吊钩动力学模型

该模型综合考虑了柔性机构作用在吊钩上的弹力、阻尼力和吊钩自身重力,并结合吊装平台的运动和第二类拉格朗日方程,通过数值迭代的方式获得吊装平台运动状态或重力改变后吊钩的摆动曲线.

2.3.1 梁单元阻尼力柔性机构阻尼为结构阻尼,与结构、材料和载荷相关,因此采用等效黏性阻尼系数近似计算,为更好地体现阻尼特性,材料的阻尼系数δ取为0.2 kg/s.

以A₁D梁单元为例,在O-XYZ中,设螺栓联接点A₁为t₁时刻固定点,接触点D为运动点;在t₂时刻分别为A'₁和D₁,t₂大于t₁.可获得梁单元作用在动点D₁的拉伸阻尼力为

F₂=δ

$\frac{(|L_{A'_{1} D_{1}}| - |L_{A_{1} D_{1}}|)}{t_{2} - t_{1}}$

(10)

同理,梁单元作用在动点D₁的弯曲阻尼力为

F'₂=δ

$\frac{θ}{(t_{2} - t_{1}) |L_{A'_{1} D_{1}}|}$

(11)

对于C₁D和C₂D两个梁单元,采用固定点和运动点轮换的方式求其对于C₁₁、D₁和C₂₁的阻尼力,且只考虑β变化产生的阻尼力.

拉伸阻尼力F₂方向为 $L_{D_{1} A'_{1}}$ ,弯曲阻尼力F'₂的方向为 $L_{A'_{1} D_{1}}$ ×( $L_{A_{1} D}$ × $L_{A'_{1} D_{1}}$ ),不难得出梁单元阻尼力在各轴上的投影为F $'_{X_{i}}$ 、F $'_{Y_{i}}$ 和F $'_{Z_{i}}$ ,其中i=1,2,…,12为运动点序号.

为方便计算,当t₁和t₂间隔较小时,认为所求阻尼力为t₂时刻阻尼力,且t=0时,阻尼力为0.最后,由式(3)和(4)易得阻尼力对于x₁轴的合力矩M $'_{x_{1}}$ 及对于转动支座轴的合力矩.

2.3.2 吊钩动力学方程设转动支座绕x₁轴的转动惯量为J,吊钩质心与回转中心距离H=h₁+h₂+h₃.在O-XYZ中,设运动坐标系原点(回转中心)坐标为(x,y,z),吊钩等效质心坐标为(x₂,y₂,z₂).为保证转动支座强度,取J=0.5 kg·m²,为保证柔性机构安装,取摆臂长l_a=h₁+h₂+10 mm,则有:

x₂=x+Hsinβy2₌y-Hcosβsinαz2=_z+Hcosβcosα

(12)

采用第二类拉格朗日方程构建吊钩动力学方程,则有:

$\frac{d}{d t} (\frac{\partial T}{\partial {\overset{\cdot}{q}}_{k}})$

$\frac{\partial T}{\partial q_{k}}$

=Q_k

(13)

k=1,2,…,N

式中:q_k为该质点系的独立广义坐标;Q_k为主动力对应于q_k的广义力;T为系统动能.本文只对吊钩相对吊装平台的转动进行分析,涉及的动能只包括摆臂、吊钩动能及转动支座转动动能.

由式(13)可得吊钩摆动加速度为

$\ddot{α}$

=[Q_α+(2/3)ρ

$l_{a}^{3} \overset{\cdot}{α} \overset{\cdot}{β}$

cosβsinβ+(1/2)ρ

$\ddot{y} l_{a}^{2}$

cosβcosα+(1/2)ρ

$\ddot{z} l_{a}^{2}$

cosβsinα+2mH²

$\overset{\cdot}{α} \overset{\cdot}{β}$

cosβsinβ+

$\ddot{my}$

Hcosβcosα+

$\ddot{mz}$

Hcosβsinα]/[(1/3)ρ

$l_{a}^{3}$

cos²β+mH²cos²β+J]

(14)

$\ddot{β}$

=[Q_β-(1/2)ρ

$\ddot{x} l_{a}^{2}$

cosβ-(1/3)ρ

$l_{a}^{3} {\overset{\cdot}{α}}^{2}$

cosβsinβ-(1/2)ρ

$\ddot{y} l_{a}^{2}$

sinβsinα+(1/2)ρ

$\ddot{z} l_{a}^{2}$

sinβcosα-

$\ddot{mx}$

Hcosβ-mH²

${\overset{\cdot}{α}}^{2}$

cosβsinβ-

$\ddot{my}$

Hsinβsinα+

$\ddot{mz}$

Hsinβcosα]/[(1/3)ρ

$L_{a}^{3}$

+mH²]

(15)

式中:Q_α= $M_{x_{1}}$ +M $″_{x_{1}}$ +M $'_{x_{1}}$ 和Q_β=M₀+M″₀+M'₀分别为弹力、重力和阻尼力对于x₁轴和转动支座轴的广义力,其数值等于合力矩.

2.3.3 吊钩运动的数值迭代在描述吊钩的非线性运动时,结合直线插补的思想,将吊钩的非线性连续运动拆分为吊装平台的平动和吊钩的摆动,如图8所示^[14-15].

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 吊钩运动分解图

Fig.8 Decomposition of motion of lifting hook

为完成分析,吊装平台的运动和吊钩的初始状态需提前给出.而吊钩的受力与吊装平台的运动、柔性机构和吊钩的结构尺寸有关,且随时间变化,难以通过显函数表示.为获得吊钩的摆动曲线,需要进行数值迭代^[16^-^17].具体迭代步骤如下.

(1) 设吊装平台位移函数为x=f(t),则有:t₁时刻到t₁+Δt时刻吊装平台位移dx=f(t₁+Δt)-f(t₁);t₁时刻吊装平台的速度为 ${\overset{\cdot}{x}}_{t_{1}}$ =f'(t₁)和加速度为 ${\ddot{x}}_{t_{1}}$ =f″(t₁).同理可得,y和z方向安装平台的位移和速度.

(2) 根据吊钩t₁时刻的运动状态及受力,由吊钩动力学方程可得:t₁时刻吊钩摆动加速度 ${\ddot{α}}_{t 1}$ ;t₁+Δt时刻的角速度 ${\overset{\cdot}{α}}_{t_{1} + Δ t}$ = ${\overset{\cdot}{α}}_{t_{1}}$ + ${\ddot{α}}_{t_{1}}$ Δt.为补偿速度损失造成的位移损失,t₁时刻到t₁+Δt时刻的角位移dα= ${\overset{\cdot}{α}}_{t_{1} + Δ t}$ Δt,β也由此得到.

Δt取为0.001 s.考虑到吊钩的稳定时间可以反映柔顺式吊钩阻尼性能,设吊钩的稳定判据为当 $\overset{\cdot}{α}$ <0.005 rad/s、Q_α<10 N·mm时,认为吊钩的α摆角稳定,稳定时间为t_α;当 $\overset{\cdot}{β}$ <0.005 rad/s、Q_β<10 N·mm时,认为吊钩的β摆角稳定,稳定时间为t_β.

2.3.4 动力学检验空载状态下,吊钩动力学模型中的吊装平台位移被设为0,初始状态的吊钩摆角和角速度被设为0,既吊装平台倾斜后吊钩在重力作用下的摆动.为验证吊钩动力学模型的正确性,对E=0、δ=0情况下的摆角进行分析并与Adams动力学软件分析结果对比,对比结果如图9所示,其中:t_si为仿真时间.由图9可知,二者拟合程度高,由于模型间仍有细微差别,误差会随时间积累,但在一定程度上仍可以证明吊钩动力学模型的正确性.

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9 吊钩动力学模型检验

Fig.9 Verification of dynamics model of lifting hook

3 摆动柔顺式吊钩参数优化

3.1 优化参数及目标

柔性机构性能主要涉及4个优化目标,如表1所示.除了前文提到的偏角要求,柔性机构在最大应力满足材料强度要求的同时,平均应力应尽量小以延长使用寿命^[18].吊钩摆动时间也应尽可能短,以缩短对接前的准备时间.

表1 摆动柔顺式吊钩优化目标

Tab.1 Optimization objectives of swing-compliant hook

工况	目标
空载	α₁、β₁最小; t_α、t_β最短
负载	α₂、β₂最大;σ_sum、σ_max最小

新窗口打开| 下载CSV

考虑到摆动柔顺式吊钩的安装尺寸和受载情况,各参数的预设值及其优化范围如表2所示.

表2 摆动柔顺式吊钩参数设置

Tab.2 Parameter settings of swing-compliant hook

参数	最小值	步距	最大值	预设值
l_c/mm	70	5	130	100
b₁/mm	10	5	70	40
b₂/mm	30	5	90	60
h₁/mm	10	5	70	40
h₂/mm	40	5	100	70
h₃/mm	110	5	170	140

新窗口打开| 下载CSV

3.2 灵敏度分析

为分析各参数对优化目标的影响,进而筛选出对摆动柔顺式吊钩性能影响较大的参数,基于柔性机构理论模型分别使用单因素分析法和二水平析因设计进行灵敏度分析.

3.2.1 单因素分析因梁单元直径d与柔性机构最终结构形式和尺寸相关,弹性模量E与材料的特性有关,本文不再对二者做讨论,取d=15 mm,E=24 MPa.以表2数据进行单因素分析,选取对性能影响较大的参数,如表3所示.从表3中可以看出,以预设值设计的摆动柔顺式吊钩的吊钩空载偏角和负载偏角接近极限值,需要对结构参数进行优化.

表3 单因素分析结果

Tab.3 Results of single-factor analysis

性能参数	初始性能参考取值	正相关参数/mm	负相关参数/mm
α₁/rad	-0.025	l_c	b₁
β₁/rad	0.005	b₁	h₁
α₂/rad	-0.189	l_c	b₁、h₁
β₂/rad	0.158	b₁	h₁
σ_sum/MPa	25.21	h₁	l_c、b₁
σ_max/MPa	5.37	b₁、 h₁	l_c
t_α/s	10.35	l_c、h₃	h₁、b₂
t_β/s	6.51	l_c、b₂、h₃	h₁

新窗口打开| 下载CSV

3.2.2 二水平析因分析采用实验设计方法对柔性机构进行二水平析因设计^[19].设计因子的选择及取值见表2.分别取吊钩空载的偏角和θ₁= $|α_{1}|$ +β₁、吊钩负载的偏角和θ₂= $|α_{2}|$ +β₂、柔性机构负载应力响应σ_r=σ_sum+σ_max以及吊钩空载稳定时间响应t_r=t_α+t_β作为响应量.利用Design-Expert软件进行实验设计,求得各参数贡献度如表4所示.

表4 参数对响应量的贡献度

Tab.4 Contribution of parameters to responses

参数	θ₁贡献度/%	θ₂贡献度/%	σ_r贡献度/%	t_r贡献度/%
l_c	7.822	2.046	-21.556	22.798
b₁	-37.772	-0.713	-2.420	-0.104
b₂	-2.454	-0.169	0.403	-8.854
h₁	-31.504	-76.900	62.038	-32.084
h₂	-0.279	-0.310	-1.241	-4.160
h₃	0.619	4.425	0	6.454

新窗口打开| 下载CSV

由表4可以看出,参数l_c和h₁对不同目标均有较大影响,参数b₂和h₂贡献度相对较小.参数b₁只对θ₁的影响较大,对其他目标影响范围较小,而参数h₃整体影响较好.因此,结合表3内容,选取l_c、h₁和h₃进行优化.

3.3 响应面模型分析

采用含交叉项的二次型响应面函数建立响应面模型^[20,21]:

$\tilde{y}$

=u₀+

$\overset{n}{\sum_{j = 1}}$

u_jx_j+

$\overset{n}{\sum_{i = 1}} \overset{n}{\sum_{j = 1}}$

u_ijx_ix_j

(16)

式中: $\tilde{y}$ 为响应面模型预测值;x_i为n维自变量的第i个分量;u₀、u_j和u_ij分别为常数项、一次项及二次项待定系数,可由最小二乘法确定.l_c、h₁和h₃作为实验因子,因子高低水平与二水平析因设计的上下限一致.

为建立表征吊钩空载时偏角、稳定时间和负载时偏角、柔性机构应力的响应面模型,定义输出响应为R_i(i=1,2,…,8),分别为 $|α_{1}|$ 、β₁、 $|α_{2}|$ 、β₂、σ_sum、σ_max、t_α和t_β.采用Box-Behnken方法布置实验点位置及拟合二阶响应面三水平设计,利用已有理论模型进行实验设计,并采用最小二乘法对参数进行拟合,可得响应面近似模型为

R=uχ+u₀

(17)

R= ${[\begin{array}{l} R_{1} & R_{2} & R_{3} & R_{4} & R_{5} & R_{6} & R_{7} & R_{8} \end{array}]}^{T}$

χ= ${[\begin{array}{l} a & h_{1} & h_{3} & a h_{1} & a h_{3} & h_{1} h_{3} & a^{2} & h_{1}^{2} & h_{3}^{2} \end{array}]}^{T}$

u₀= ${[\begin{array}{l} 0.002 & 0.010 & 0.160 & 0.186 & 37.111 & 6.712 & - 1.777 & - 2.013 \end{array}]}^{T}$

u= $[\begin{array}{l} 8.47 \times 10^{- 6} & - 2 \times 10^{- 5} & 4.02 \times 10^{- 5} & - 4.5 \times 10^{- 6} & 1.17 \times 10^{- 6} & - 1.4 \times 10^{- 6} & 2.15 \times 10^{- 6} & 3.73 \times 10^{- 6} & - 2.6 \times 10^{- 8} \\ - 1.1 \times 10^{- 4} & - 1.2 \times 10^{- 4} & 5.07 \times 10^{- 5} & 1.91 \times 10^{- 7} & - 4.8 \times 10^{- 8} & - 3.7 \times 10^{- 7} & 4.71 \times 10^{- 7} & 9.12 \times 10^{- 7} & - 3.5 \times 10^{- 8} \\ 4.67 \times 10^{- 4} & - 3.8 \times 10^{- 4} & 9.65 \times 10^{- 5} & 2.44 \times 10^{- 6} & - 5.6 \times 10^{- 7} & 3.13 \times 10^{- 7} & - 1.6 \times 10^{- 6} & - 4.6 \times 10^{- 7} & - 7.6 \times 10^{- 8} \\ - 5.4 \times 10^{- 4} & - 4.3 \times 10^{- 4} & 1.6 \times 10^{- 4} & - 6.7 \times 10^{- 7} & 1.83 \times 10^{- 7} & 5.86 \times 10^{- 7} & 2.35 \times 10^{- 6} & 2.37 \times 10^{- 7} & - 3.2 \times 10^{- 7} \\ - 0.3169 & 0.3589 & 2.71 \times 10^{- 3} & - 1.63 \times 10^{- 3} & 0 & 0 & 1.12 \times 10^{- 3} & 3.64 \times 10^{- 4} & - 9.7 \times 10^{- 6} \\ - 0.0442 & 0.0843 & - 4.7 \times 10^{- 4} & - 3.5 \times 10^{- 4} & 0 & 0 & 1.12 \times 10^{- 4} & 7.86 \times 10^{- 6} & 1.66 \times 10^{- 5} \\ - 0.0027 & 0.1275 & 1.31 \times 10^{- 3} & - 2.45 \times 10^{- 3} & 7.81 \times 10^{- 4} & - 7.8 \times 10^{- 4} & 8.47 \times 10^{- 4} & 1 \times 10^{- 3} & 1.61 \times 10^{- 5} \\ 0.0375 & 0.006 & 0.018 & - 1.6 \times 10^{- 5} & 3.59 \times 10^{- 4} & - 5.8 \times 10^{- 4} & - 2 \times 10^{- 4} & 4.13 \times 10^{- 4} & 7.76 \times 10^{- 5} \end{array}]$

式中: R为响应列阵; χ为实验因子列阵; u₀为常数项列阵;u为总系数矩阵.

同时,求得复相关系数R²和修正的复相关系数 $R_{adj}^{2}$ ,用于检验响应面近似模型对实验数据的拟合度.R²趋于1,表示响应面近似模型与原模型相似度高;R²与 $R_{adj}^{2}$ 越接近,表示模型不含冗余参数,分析结果如表5所示.由表5可知,该模型可以准确描述设计参数与响应之间的关系^[20].

表5 响应面模型误差分析

Tab.5 Error analysis of response surface model

输出响应	R²	$R_{adj}^{2}$
R₁	0.9991	0.9980
R₂	0.9984	0.9963
R₃	0.9969	0.9930
R₄	0.9994	0.9987
R₅	0.9996	0.9991
R₆	0.9999	0.9999
R₇	0.9985	0.9966
R₈	0.9981	0.9956

新窗口打开| 下载CSV

3.4 优化分析

以表1的性能为目标进行参数优化.通过Design-Expert软件中的响应面模型求解实验因子的优化值,部分实验因子组合如表6所示.h₃的值是优化范围的最小值,因此不是最优值,考虑到吊钩尺寸h₃=110 mm.对接性能中θ₁对于性能的影响较大且工程中θ₂更容易满足,参数l_c对响应θ₁为正影响,参数h₁对响应θ₁为负影响.为了使θ₁趋于最小值,选取l_c=90 mm、h₁= 23 mm.参数b₁和b₂对响应θ₁为负影响,考虑到结构尺寸和加工制作,二者均取为70 mm.参数h₂对响应θ₁和t_r为负影响,取其最大值为100 mm.

表6 实验因子的部分优化解集

Tab.6 Partial optimal solutions of the experimental factors

序号	l_c/mm	h₁/mm	h₃/mm
1	90.61	22.39	110
2	94.87	23.32	110
3	106.05	28.49	110

新窗口打开| 下载CSV

经由理论模型分析,优化后的摆动柔顺式吊钩性能如表7所示.优化后的摆动柔顺式吊钩满足吊钩空载和负载的偏角要求,同时空载时吊钩稳定时间和负载时柔性机构应力较小.与初始性能对比可以看出,吊钩空载和负载偏角得到了一定程度的优化,同时柔性机构应力和吊钩稳定时间均得到了改善.

表7 优化结果比较

Tab.7 Comparison of optimization results

性能参数	初始性能参数取值	优化后性能参数取值
α₁/rad	-0.025	-0.012
β₁/rad	0.005	0.007
α₂/rad	-0.189	-0.180
β₂/rad	0.158	0.169
σ_sum/MPa	25.21	20.33
σ_max/MPa	5.37	4.91
t_α/s	10.35	8.38
t_β/s	6.51	5.73

新窗口打开| 下载CSV

4 结论

经过上述研究,本文主要结论如下.

(1) 对现有的发射箱装填方式进行分析,提出一种用于发射箱装填的新型柔性装填方案,并在此基础上设计了摆动柔顺式吊钩.

(2) 为描述摆动柔顺式吊钩的性能,在节点位移法基础上建立了以结构参数为变量的参数化理论模型,如吊钩位移模块、柔性机构应力模块和吊钩动力学模型.

(3) 为了综合优化摆动柔顺式吊钩性能,根据筛选出对性能影响较大的关键参数构建表征摆动柔顺式吊钩性能的响应面模型.优化结果表明,当柔性机构长度为90 mm,安装高度为23 mm,吊钩末端高度为110 mm时,该摆动柔顺式吊钩的综合性能最优.

(4) 本文提出的节点位移法建模、筛选关键结构参数、响应面模型优化研究方法适用于机械系统关键部分的结构参数优化,可为进一步的优化(如拓扑优化)提供理论依据.

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

李军, 马大为, 曹听荣, 等. 火箭发射系统设计[M]. 北京: 国防工业出版社, 2008.