一种基于虚拟接地控制及结构协同优化的控制方法与流程

本发明涉及机械结构与控制算法转化领域，特别是涉及一种基于虚拟接地控制及结构协同优化的控制方法。

背景技术：

振动是指物体的全部或一部分沿直线或曲线做有一定时间规律或周期的往返颤动，是自然界普遍存在的现象。然而对大多数机械元件来说，振动会降低元件的精度，缩短元件的使用寿命，严重时会造成螺栓等固结件的松动，对操作人员带来危险，影响生产效率。因此，如何消除振动对元件的影响是大多数机械设备所要面对的问题。如果能够将需要隔振的物体悬挂在如天花板、天空等不受振动影响的稳定物体上，或固定在如地面之类的振动影响对其小到可以忽略不计的庞大物体上，将机械元件的振动转换到这些物体上，则可达到减小振动对机械设备影响的目的。这其中无论是将机械元件悬挂在天空中还是固定在地面上，从现实意义来说，都是将机械元件固定在大地上，以地球作为最大的吸振体弱化振动的影响。

对固定不动的机械设备来说，固定在大地上减振是易于实现的，然而对大多数移动机械设备，如飞机，轮船，火车等，将机械元件一直固定在地面上是不现实的，所以需要通过机械结构和控制算法来协同实现。当前，运用在机械部件中的控制算法，大多数都以固定的结构为基准，优化结构参数从而实现控制目标，或以电机为直接驱动装置，通过优化电机输出力实现机械元件的所需力。在这两种实现方式中，电机作为直接驱动装置的方法又因其较高的成本而成为了第二选择。基于结构并优化结构参数的控制方法虽然有着低成本高性能的优势，但这样的选择常常无法保证作为基础的结构是能够让控制算法充分发挥其效能的最优结构，从而导致控制算法对隔振效果的提升有一定的局限性。

在机电相似性理论基础下，惯性的思想被提出，并在弹簧元件、阻尼元件的基础上提出了惯性元件，使得机电网络能够严格对应起来，同时，许多基于电网络的研究方法得以被运用在机械网络的研究上。为机械结构的研究提供了更多的可能。无论是当前基于结构并优化结构参数的控制方法还是基于控制算法优化结构，都无法使得控制参数和结构在同一优化目标下同时为最优选择，这局限了机械设备性能的提升，也使得隔振研究陷入瓶颈。

技术实现要素：

针对现有技术中存在不足，本发明提供了一种基于虚拟接地控制及结构协同优化的控制方法，在接地隔振这一广泛控制概念下，在同一优化目标下，从优选择弹簧、阻尼、惯性单个或多个元件的组合控制方法，并优化机械结构，实现从控制算法到机械结构的协同优化，选出一种最适控制方法和结构，从而达到真正的控制最优。

本发明是通过以下技术手段实现上述技术目的的。

一种基于虚拟接地控制及结构协同优化的控制方法，包括以下步骤：

步骤1：在不影响分析结果的前提下，将整个设备系统拆分为多个单元设备系统，对单个单元设备系统进行分析；

步骤2：假设单元设备系统通过接地结构固定于地面，接地结构由弹簧、阻尼、惯容器三元件构成；设备通过隔振层与基底相连，振动源通过基底将振动传送到隔振层上；

所述步骤2中的接地结构的速度型阻抗表达式y(s)表示为其中i＝1,2,…，表示第i个机械元件，f为函数符号，cki为第i个接地结构的阻尼系数，kki为第i个接地结构的弹簧刚度，qki为第i个接地结构的惯性元件的系数，s为复变数，s＝σ+jω，σ表示复变数s的实部，ω表示复变数s的虚部，j为虚数单位。

步骤3：用动力学方程表示单元设备系统模拟接地模型；

所述步骤3中的动力学方程为：其中，ms为设备的质量，z0表示基底位移，即为振动输入信号，z1表示隔振层位移，z2表示设备位移，则隔振层位移z1＝z2-z0。

步骤4：以设备在垂直振动方向上的振动行程对基底位移的传递函数gi(s)，以及设备在垂直振动方向上的振动加速度对基底位移的传递函数gii(s)为优化目标，采用智能优化算法分别对接地结构和隔振层的速度型阻抗表达式的系数进行优化，当且时，停止优化，其中m表示接地结构的阶数，n表示隔振层的阶数，m、n均为正整数，且m≤m，n≤n，m为接地结构所允许的最大阶数，n为隔振层所允许的最大阶数；

所述步骤4中设备在垂直振动方向上的振动行程对基底位移的传递函数gi(s)为：设备在垂直振动方向上的振动加速度对基底位移的传递函数gii(s)为：

所述步骤4中的智能优化算法可采用改进鱼群算法。

以改进鱼群算法为例，其步骤如下：

(1)鱼群初始化

鱼群中每条人工鱼为给定范围内随机产生的一组实数，设鱼群的大小为n，待优化的参数为αj，取值范围为0＜αj＜1，产生一个1行n列的初始鱼群，每列数值表示一条人工鱼的参数αj；

(2)觅食行为

假设人工鱼群的当前状态为xp，随机选取感知范围内的一个状态xq，x表示人工鱼个体的状态位置，在求解极大值问题中，如果第p条人工鱼所在位置的食物浓度小于第q条人工鱼所在位置的食物浓度yp<yq，其中y表示人工鱼当前所在位置的食物浓度，则向该方向前进一步：

其中xnext表示鱼群向前进一步的位置，rand()表示取值范围内的随机数，step表示鱼群移动的最大步长。反之，则重新选择随机状态xq，xnext＝xp+rand()·step。

(3)聚群行为

假设人工鱼的当前状态xp，探索当前领域内满足||xq-xp||<visual的伙伴数目nf和中心位置xc，其中visual表示人工鱼的感知距离，假如δ表示拥挤度，则表明在当前领域内，中心位置有较多食物且不太拥挤，则朝中心位置的方向前进一步，反之，再次进行觅食行为。

(4)追尾行为

假设人工鱼当前状态xp，探索当前领域内满足||xq-xp||<visual的伙伴数目nf，这些伙伴中yp为最大的伙伴xp，若满足则表明xq的状态有较高的食物浓度且周围环境不太拥挤，则朝xp的方向进一步，反之，重复执行觅食行为。

(5)随机行为

在视野范围内随机选取一个状态，朝着该方向移动，此过程可认为是觅食过程的一个缺省行为，即：

xp|next＝xp+rand()·visual

步骤5：优化后分别得到含有具体数值的接地结构的速度型阻抗表达式y′(s)和隔振层的速度型阻抗表达式t′(s)，判断y′(s)、t′(s)是否满足约束条件，若满足，则继续步骤6；若不满足约束条件，则返回步骤4，继续优化y(s)、t(s)的系数，直到满足该约束条件；

所述步骤5中的约束条件为：当复变数s的虚部等于0时，即im(s)＝0时，im[y′(s)]＝0且im[t′(s)]＝0；且当复变数s的实部大于等于0时，即re(s)≥0时，re[y′(s)]≥0且re[t′(s)]≥0；其中，im为复数中的虚数部位简称，re为复数中的实数部位简称。

步骤6：将接地结构下放到隔振层，并采用控制结构模拟实现接地结构的效果，确保控制结构产生的作用力与接地结构产生的作用力一致；

所述步骤6中的控制结构的速度型阻抗表达式k(s)与y′(s)之间满足：

步骤7：根据串并联结构的速度型阻抗表达式的基本定律，通过带有具体数值的速度型阻抗表达反推求解得到控制结构和隔振层各元件之间的连接关系，从而得到控制结构和隔振层的具体物理结构。

本发明的有益效果在于：

1.一种基于虚拟接地控制及结构协同优化的控制方法实现了设备系统的虚拟接地，最大程度上减小了振动对设备的影响，提升了设备的工作稳定性，延长了机械元件的寿命。

2.采用控制结构实现模拟接地的方法使得设备无需真正固定于地面，拓展了应有领域，使得该控制方法不仅可以应用在非移动设备上，也可应用在实时移动设备上，实现最优隔振效果。

3.一种基于虚拟接地控制及结构协同优化的控制方法突破了原有“先结构后优化”的模式，在优化参数的同时优化结构，使得最终优化获得的结构与参数同时最优，也最为协同。

附图说明

图1为整个设备系统模拟接地控制模型

图2为单元设备系统接地模型

图3(a)为y′(s)降阶后的结构框图，图3(b)为t′(s)降阶后的结构框图

图4为单元设备系统模拟接地控制模型

具体实施方式

下面结合附图以及具体实施例对本发明作进一步的说明，但本发明的保护范围并不限于此。

如图1所示，整个设备系统由许多个如图2所示的单元设备系统组合而成的，在分析整个设备系统时，在不影响分析结果的前提下，可将整个设备拆分为多个单元设备系统进行分析。

如图2所示，假设单元设备系统通过接地结构固定于地面，接地结构由弹簧、阻尼、惯性三元件构成；设备通过隔振层与基底相连，振动源通过基底将振动传送到隔振层上。接地结构的速度型阻抗表达式y(s)可用弹簧元件、阻尼元件、惯性元件综合表示，其速度型阻抗表达式可表示为：其中i＝1,2,…，表示第i个机械元件，f为函数符号，cki为第i个接地结构的阻尼系数，kki为第i个接地结构的弹簧刚度，qki为第i个接地结构的惯性元件的系数，s为复变数，s＝σ+jω，σ表示复变数s的实部，ω表示复变数s的虚部，j为虚数单位。同样的，隔振层的速度型阻抗表达式为：其中kgi表示隔振层的第i个弹簧刚度，cgi隔振层的第i个阻尼系数，qgi隔振层的第i个惯性元件的系数。

可根据结构的速度型阻抗表达式s的最高阶数定义结构的阶数，假设接地结构为m阶结构，隔振层是n阶结构，其中m、n均为正整数，且m≤m，n≤n，m为接地结构所允许的最大阶数，n为隔振层所允许的最大阶数，则接地结构的速度型阻抗表达式y(s)可表示为：

其中β、α均为系数；

隔振层的速度型阻抗表达式t(s)可表示为：

其中λ、τ均为系数。

单元设备系统模拟接地模型的动力学方程可表示为：

其中ms为设备的质量，z0表示基底位移，即为振动输入信号，z1表示隔振层位移，z2表示设备位移，则隔振层位移z1＝z2-z0。

单元设备系统接地时，设备在垂直振动方向上的振动行程对基底位移，即振动输入信号的传递函数gi(s)为：

单元设备系统模拟接地时，设备在垂直振动方向上的振动加速度对基底位移，即振动输入信号的传递函数gii(s)为：

gi(s)、gii(s)的值越接近于0，设备受到振动的影响越小。

以gi(s)、gii(s)为优化目标，可选择鱼群算法、遗传算法、神经网络等智能优化算法，对y(s)、t(s)速度型阻抗表达式的系数β、α、λ、τ进行优化，当且时，停止优化，将β、α、λ、τ、m、n的数值代回y(s)和t(s)，则接地结构的速度型阻抗表达式可通过优化后的参数具体表示为速度型阻抗表达式y′(s)，则隔振层的速度型阻抗表达式可通过优化后的参数具体表示为速度型阻抗表达式t′(s)。

判断速度型阻抗表达式y′(s)、t′(s)是否满足约束条件：(1)当复变数s的虚部等于0时，即im(s)＝0时，im[y′(s)]＝0且im[t′(s)]＝0；(2)当复变数s的实部大于等于0时，即re(s)≥0时，re[y′(s)]≥0且re[t′(s)]≥0。im为复数中的虚数部位简称，re为复数中的实数部位简称。若y′(s)、t′(s)能够满足上述2个约束条件，则继续下一步，若不满足约束条件，则返回上述优化步骤，继续优化y(s)、t(s)的系数，直到满足该判断条件。

接地结构所允许的最大阶数m与隔振层所允许的最大阶数n可根据设备工作环境、结构设计、结构所允许的最大元件数量、安装空间等实际工程运用情况来决定。虽然结构的速度型阻抗表达式的阶数越高，隔振性能越好，但随着阶数的升高，与之对应的机械结构也会随之变得复杂，用于实现结构的元件数量也随之增长。为了确保结构的可实施性和结构的合理性，将结构的速度型阻抗表达式的阶数控制在3以下，因此单元设备系统接地模型在m＝3且n＝3时的情况最复杂。

当m＝3且n＝3时，a、b、c、d、e、f、g均为及优化后的接地结构的具体的速度型阻抗表达式的系数，a’、b’、c’、d’、e’、f’、g’均为及优化后的隔振层结构的具体的速度型阻抗表达式的系数，为了易于通过速度型阻抗表达式求解元件结构，分别对y′(s)、t′(s)采用以下降阶步骤：

(1)分别串联或并联地提取出一个阻尼，y′(s)中串联的阻尼系数记作cuy，y′(s)中并联的阻尼系数记作cby；t′(s)中串联的阻尼系数记作cut，t′(s)中并联的阻尼系数记作cbt；其中cby＝0或cuy为无穷大，cbt＝0或cut为无穷大。

(2)并联地提取出一个惯性元件，y′(s)中并联的惯性元件的系数记作quy；t′(s)中并联的惯性元件的系数记作qut。

通过上述步骤分别得到y′(s)、t′(s)降阶后的结构框图，分别如图3(a)(b)所示。

则y′(s)降阶后的速度型阻抗表达式可表达为y″(s)：

则t′(s)降阶后的速度型阻抗表达式可表达为t″(s)：

其中，速度型阻抗表达式y1(s)、t1(s)均为双二次型速度型阻抗表达式，它们的标准速度型阻抗表达式分别为：

其中ay，by，cy，dy，ey，f’y是y1(s)的系数，且它们与quy，cby，cuy满足以下条件：

a＝quycuydy；b＝cuy(quyey+cbydy+ay)；

c＝cuy(quyfy+cbyey+by)；d＝cuy(cbyfy+cy)；

1＝quydy；e＝cuydy+quyey+cbydy+ay；

f＝cuyey+quyfy+cbyey+by；g＝cuyfy+cbyfy+cy；

at，bt，ct，dt，et，ft是t1(s)的系数，且它们与qut，cbt，cut满足以下条件：

a′＝qutcutdt；b′＝cut(qutet+cbtdt+at)；

c′＝cut(qutft+cbtet+bt)；d′＝cut(cbtft+ct)；

1＝qutdt；e′＝cutdt+qutet+cbtdt+at；

f′＝cutet+qutft+cbtet+bt；g′＝cutft+cbtft+ct；

且满足：

求解得到quy，cby，cuy，qut，cbt，cut以及各系数ay，by，cy，dy，ey，f’y，at，bt，ct，dt，et，ft的数值，将这些数值带回表达式y1(s)、t1(s)、y″(s)、t″(s)，从而得到带有具体数值的速度型阻抗表达式y″(s)、t″(s)，分别记为表达式y″′(s)、t″′(s)。

将接地结构下放到隔振层，并采用控制结构模拟实现如图4所示的接地结构的效果，采用主动或半主动的控制方式使得控制结构产生的作用力与接地结构产生的作用力一致，则控制结构的速度型阻抗表达式k(s)与y′(s)之间满足：

其中z0、z2分别为z0、z2的拉氏变换。

根据串并联结构的速度型阻抗表达式的基本定律，即若弹簧元件、阻尼元件、惯性元件，用c表示阻尼的系数，k表示弹簧的系数，m表示惯性元件的系数，若三元件并联，则三元件并联结构的速度型阻抗表达式ybing可表示为：若三元件串联，则三元件串联结构的速度型阻抗表达式ychuan可表示为：可通过带有具体数值的速度型阻抗表达式k(s)、t″′(s)反推求解得到控制结构和隔振层各元件之间的连接关系，从而得到控制结构和隔振层的具体物理结构。

在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示意性实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。

尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。