CN108897213B - 一种基于面内调控的褶皱薄膜主动振动控制方法 - Google Patents
一种基于面内调控的褶皱薄膜主动振动控制方法 Download PDFInfo
- Publication number
- CN108897213B CN108897213B CN201810523112.3A CN201810523112A CN108897213B CN 108897213 B CN108897213 B CN 108897213B CN 201810523112 A CN201810523112 A CN 201810523112A CN 108897213 B CN108897213 B CN 108897213B
- Authority
- CN
- China
- Prior art keywords
- matrix
- vibration
- model
- film
- order
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明涉及一种基于面内调控的褶皱薄膜主动振动控制方法,主要是通过分析褶皱形态对薄膜结构动态特性的影响,建立了褶皱薄膜的动态模型;利用模态叠加法得到以固有振型为基础的空间向量,建立了不同载荷下褶皱薄膜的动力学模型;根据各阶模态对系统振动变形的贡献,对动力学模型进行缩聚,从而建立了面向控制的薄膜结构动力学模型;采用模型参考自适应控制器,结合PID控制器,提出了一种基于面内调控多模型切换的PID‑MRAC控制方法,本发明的方法无需在薄膜表面粘贴作动器,可有效地抑制薄膜面外振动,不会影响薄膜本身的特性,保证薄膜结构的稳定性,比传统的PID控制方法下的振动时间缩短了约39.7%。
Description
技术领域
本发明属于控制工程技术领域,具体涉及一种基于面内调控的褶皱薄膜主动振动控制方法,用于更简单且更有效地抑制薄膜振动。
背景技术
褶皱的出现显著降低了薄膜的表面精度,对薄膜结构的动态特性也存在一定影响,在无空气阻尼的太空环境中,薄膜结构一旦受到外界扰动,膜面振动将很难停止,加之薄膜本身的柔性特性对保证在轨稳定性和可控性带来一定难度,因此对褶皱薄膜的调控面临着严峻的挑战。国内外学者已开展的研究,目前,已有的动态调控方法主要有两类。
一类是通过在薄膜表面粘贴作动器来抑制薄膜的振动。Renno针对薄膜条结构,提出了一种非线性控制方法。Ruggiero在薄膜镜表面分布多个压电作动器,设计了LQR(linear quadratic regulator)控制系统,Ferhat在此基础上提出了基于变分原理的求解方法,随后,针对一个四边固支的方形薄膜结构,通过在表面粘贴四个压电作动器,达到抑制振动的目的。然而对于超薄的薄膜材料,直接在其表面粘贴作动器,很难实现合理的布线,粘接处很容易出现褶皱,同时粘贴作动器后会改变局部的厚度和刚度,继而改变薄膜应力分布,最终影响薄膜结构动态特性。
另一类是通过在设置面外调节拉索,来抑制膜面振动。Sakamoto提出了索网薄膜结构,并且针对索网结构提出了一种LQR控制方法,通过在边索安装垂直于膜面方向的作动器,同时在拉索和薄膜连接处安装面内作动器,来抑制振动。如果在每一个连接点都安装作动器,会大大增加薄膜结构的整体质量,而且垂直膜面方向的作动器在实际工程应用中不易实现。无论是在薄膜上粘贴作动器,还是设置面外调节索,都是通过施加面外作用力抑制振动,需要安装附加结构,工程实现上难度大。因此,亟需一种合理的控制方法,不需要在薄膜上粘贴作动器,即可达到抑制振动的目的。
综上所述,对于薄膜结构的动态调控方面,主要是通过在薄膜表面或者在边界拉索上安装压电作动器,通过面外控制达到抑制振动的目的。不仅布线和安装困难,直接在薄膜表面粘贴作动器,还会改变薄膜的应力分布,影响结构动态性能。因此不使用面外作动器,仅通过面内控制有效地抑制薄膜结构面外振动成为亟待解决的问题。
发明内容
本发明的目的在于避免上述现有技术分析方法的不足,提供一种基于面内调控的褶皱薄膜主动振动控制方法,达到易于工程实现,同时能够更有效地抑制薄膜面外振动。
本发明所采用的技术方案是:
一种基于面内调控的褶皱薄膜主动振动控制方法,其包括以下步骤:
(1)引入褶皱变形函数,建立褶皱薄膜的振动模型;
(2)用有限元屈曲分析方法确定出褶皱薄膜的褶皱形态,更新有限元分析模型,并进行模态分析,得到步骤(1)振动模型在引入褶皱信息后的振动特性;
(3)根据步骤(2)中得到的振动特性,用模态叠加法得到以固有振型为基础的空间向量,建立不同载荷下褶皱薄膜的动力学模型;
(4)根据各阶模态对系统振动变形的贡献,对步骤(3)不同载荷下褶皱薄膜的动力学模型进行缩聚,建立面向控制的薄膜结构动力学模型;
(5)结合步骤(4)的面向控制的薄膜结构动力学模型以及驱动装置的动力参数,建立驱动模型;
(6)以步骤(5)建立的驱动模型作为被控对象,设计自适应参考控制模型,得到自适应参考控制模型的输出电压,加载到被控对象上,实现性能指标最优条件下的前馈补偿;
(7)根据PID控制模型以及步骤(6)的自适应参考控制模型,得到多模型切换的PID-MRAC控制模型,并对其输出电压进行限幅滤波,加载到被控对象上,实现褶皱薄膜的主动振动控制。
进一步限定,所述步骤(1)的具体方法为:
(1.1)仅考虑薄膜的面外振动,令振动变形函数为wv,褶皱变形函数为wz,得到振动过程中褶皱薄膜弯曲形变势能Ep1和褶皱薄膜的应变势能Ep2以及褶皱薄膜的自由振动动能Ev,分别表达为:
其中,εx是X方向应变,εy是Y方向应变,γxy是切应变,E是薄膜材料的杨氏模量,v是泊松比,hm是薄膜厚度;ρ为薄膜的密度;
(1.2)将振动变形函数wv和褶皱变形函数wz都表示成级数形式:
(1.3)将步骤(1.2)的振动变形函数wv和褶皱变形函数wz带入步骤(1.1)的振动过程中总势能Ep和自由振动动能Ev,则变形为:
其中kijkl,kij,k′ijkl,mij均为中间变量;
(1.4)结合步骤(1.3)的振动过程中总势能Ep和自由振动动能Ev,根据拉格朗日方程,得到非线性振动的微分方程:
(1.5)根据步骤(1.4)的非线性振动的微分方程,得到非线性振动的矩阵表达式为:
其中,M是质量矩阵,K是刚度矩阵。
进一步限定,所述步骤(2)的具体方法为:
2.1)采用有限元屈曲分析法得到褶皱薄膜在一定载荷作用下的褶皱形态;
2.2)将步骤2.1)所得到的褶皱形态更新到有限元分析模型中,得到褶皱变形下的应力分布以及包含褶皱信息的刚度矩阵;
2.3)利用步骤2.2)所得到的褶皱变形下的应力分布以及包含褶皱信息的刚度矩阵,进行模态分析,扩展模态,得到褶皱薄膜结构的各阶固有频率和模态振型。
进一步限定,步骤(2.2)中所述的刚度矩阵K包括初始应力刚度阵Kσ、线性刚度阵KL和褶皱引起的大变形刚度阵KNL。
进一步限定,在步骤(2.3)具体为:结合步骤(2.2)得到的包含褶皱信息的刚度矩阵,对步骤(1.5)非线性振动的矩阵表达式进行求解,得到下述方程:
(K-ω2M)φejωt=0
其中,ω为考虑薄膜表面褶皱后的一组频率值,将其记为ω1,ω2,…,ωn,组成固有频率矩阵Ω为:
进而得到第i阶模态振型:{φ1,φ2,…,φn},每一阶振型中包含nd个节点,各阶振型组成的模态振型矩阵为:
中φij表示第i阶模态第j个节点的面外位移。
进一步限定,所述步骤(3)具体为:
3.1)根据步骤2.3)的褶皱薄膜结构的各阶固有频率和模态振型以及非线性振动的矩阵表达式求解方程,得到第i阶模态质量Mmi为:
Ei为第i阶模态的动能,ωi为第i阶模态的固有频率;
3.2)根据比例阻尼得到的阻尼矩阵,即得到阻尼矩阵与质量矩阵和刚度矩阵的关系式为:
Dm=α1Km+α2Mm
其中α1和α2是瑞丽阻尼系数,且
其中ω1和ω2分别为第一阶和第二阶的固有频率值,ξ1和ξ2为第一阶和第二阶的模态阻尼比系数;
3.3)设定nd为有限元模型的自由度,输入自由度为nr,输出自由度为ns,则褶皱变形下的应力分布在广义坐标下用动力学方程表示为:
其中q是由各节点的面外位移向量组成的nd×1维的矩阵,是nd×1维的节点速度矩阵,是nd×1维的节点加速度矩阵,V=[VX VY]T是输入位移载荷矩阵,即输入自由度nr=2,M是nd×nd维的质量矩阵,D为nd×nd维的阻尼矩阵,K为nd×nd维的刚度矩阵,B为nd×nr维输入矩阵,C为1×nd维面外振动位移输出矩阵,y为输出变量;
3.4)将步骤3.3)广义坐标下用动力学方程转换到以固有振型为基础的空间向量,建立不同载荷下褶皱薄膜的动力学模型,表达式为:
其中:
其中Mm,Km,Dm,Cm,Bm分别为转换后的模态质量矩阵、模态刚度矩阵、模态阻尼矩阵、模态输出矩阵和模态输入矩阵。
进一步限定,步骤(4)具体为:
4.1)根据模态能量法对模型进行缩聚,第i阶模态能量表示为:
则,前G阶中各阶所占的模态能量比ηi为:
4.2)提取前G阶中各阶模态能量,计算出各阶能量占总能量的比率,将各阶能量占总能量的比率按大小顺序排列并相加,当模态能量之和达到总能量的90%时,即可取出相应的阶次;
4.3)根据步骤4.2)确定的阶次,截取步骤3.3)得到的不同载荷下褶皱薄膜的动力学模型中对应的数据,建立面向控制的薄膜结构动力学模型:
其中A′p,B′p,C′p对应为面向控制的薄膜结构动力学模型的状态矩阵、输入矩阵和输出矩阵。
进一步限定,步骤(5)具体为:结合步骤4.3)建立的面向控制的薄膜结构动力学模型以及驱动装置的动力参数,拟合出驱动模型:
进一步限定,步骤(6)具体为:
6.1)令驱动模型的输出跟踪误差ey为:ey=y-yr
yr为参考输入;
其中
6.2)在扩展模型中将输出变量y作为一个状态变量,构建自适应参考控制模型,表达为:
其中参考控制模型的状态矩阵为:
矩阵Px满足条件:
Px(tf)=QT
tf为系统的终端时刻,Q,QT为半正定对称的代价权衡矩阵,R为正定对称的控制输入加权矩阵。
6.3)根据步骤6.2)中构建的自适应参考控制模型,选取自适应控制律为:
进而得到,自适应参考控制模型的输出电压为:
Θ=KxT,Θd是参数Θ=KxT,Θd使用自适应律更新的估计值。
进一步限定,步骤(7)具体为:建立PID控制模型,结合步骤6.2)中建立的自适应参考控制模型,得到多模型切换的PID-MRAC控制模型,其输出电压为:
ua=λ1ua1+λ2ua2
其中,ua2为PID控制模型的输出电压,λ1和λ2为控制输入的权系数。
本发明提供的基于面内调控的褶皱薄膜主动振动控制方法主要是通过分析褶皱形态对薄膜结构动态特性的影响,利用模态叠加法得到以固有振型为基础的空间向量,建立了不同载荷下褶皱薄膜的动力学模型;根据各阶模态对系统振动变形的贡献,对动力学模型进行缩聚,建立了面向控制的薄膜结构动力学模型;在此基础上,采用模型参考自适应控制器,实现性能指标最优条件下的前馈补偿控制,提高系统的抗扰动能力,然后结合PID控制模块,提出了一种基于面内调控多模型切换的PID-MRAC控制方法,与现有技术相比,本发明的有益效果为:1)本发明考虑了褶皱对薄膜结构动态特性的影响,利用模态叠加法建立了不同载荷下褶皱薄膜的动力学模型,然后根据各阶模态对系统振动变形的贡献,对动力学模型进行缩聚,从而建立了面向控制的薄膜结构动力学模型,本发明的缩聚模型可以准确地描述出褶皱薄膜的振动特性,为薄膜结构的主动振动控制奠定了基础。2)本发明的方法无需在薄膜表面粘贴作动器,仅通过调节施加在顶点处的位移载荷,改变薄膜结构本身振动特性,即可有效地抑制薄膜面外振动,无需在薄膜表面粘贴作动器,不会影响薄膜本身的特性,且工程上易于实现。3)本发明基于缩聚后的动力学模型设计了多模型切换自适应控制器,实现性能指标最优条件下的前馈补偿控制,结合PID控制模型,设计了多模型切换的PID-MRAC控制系统,可有效地抑制膜面振动,保证薄膜结构的稳定性。4)本发明的主动振动控制方法可以在薄膜受到单次扰动时将振动时间缩短到无控制时的1/4以下,比传统的PID控制方法下的振动时间缩短了约39.7%,比已有文献中提出的在膜面粘贴作动器的面外控制方法得到的振动时间缩短了约18.9%;当薄膜受到多次扰动时,有效抑制振动的同时,保证最终施加在薄膜结构上的位移载荷达到一定值后不再持续增大,避免因拉力过大而拉断薄膜的情况出现,同时达到平衡时薄膜的抗扰动能力较强,为初始预应力的选取提供依据。
附图说明
图1为扰动施加示意图。
图2为VX=1.0mm,VY=0.25mm时褶皱形态。
图3为P点振动曲线。
图4为PID-MRAC控制框图。
图6为无控制中心点振动曲线。
图7为有控制中心点振动曲线。
图8为膜面粘贴作动器示意图。
图9为振动抑制对比。
图10为多次扰动时位移载荷曲线。
图11为第二次和第三次扰动时振动曲线。
图12为第四次和第五次扰动时振动曲线。
具体实施方式
结合附图和实施例对本发明的技术方案进行详细说明,但是本发明不仅限于下述的实施情形。
本发明的基于面内调控的褶皱薄膜主动振动控制方法可以在静态实验装置的基础上增加了两组音圈电机驱动模组,分别用于调节X方向和Y方向的位移量,通过音圈电机滑块伸缩来调节拉索位移,改变拉索上的拉力,用固连在音圈电机滑块上的拉力传感器实时监控拉力大小,光栅尺实时输出反馈移动位移量。在薄膜背面安装一个敲击装置,通过撞击法模拟面外扰动,利用PSD高精度位置传感器实时反馈薄膜表面特定的点处的面外振动变形量。将测得的所有变量通过数据采集卡传送到计算机,借助labview软件编程实现控制器算法,并采用串口通信向音圈电机发送指令,进而实现褶皱薄膜的主动振动控制。
现以边长L=500mm的正方形薄膜为例进行说明,薄膜材料具体参数在表1中列出,在四个顶点处施加两组大小相等,方向相反的位移载荷VX、VY,边界条件为四角固支。
表1薄膜材料参数
杨氏模量(GPa) | E=2.5 |
厚度(μm) | h<sub>m</sub>=25 |
泊松比 | v=0.34 |
密度(kg/m<sup>3</sup>) | ρ=1400 |
具体的基于面内调控的褶皱薄膜主动振动控制方法是通过下述步骤实现:
S1:引入褶皱变形函数,建立褶皱薄膜的振动模型;具体为:
S1.1)对于具有褶皱的薄膜结构,不考虑外界激励和外界阻尼时,薄膜结构自由振动,假如忽略薄膜的面内振动,仅考虑薄膜的面外振动,令振动变形函数为wv,褶皱变形函数wz,则振动过程中褶皱薄膜弯曲形变势能Ep1可以表示为:
其中
E是薄膜材料的杨氏模量,ν是泊松比,hm是薄膜厚度。
褶皱薄膜的应变势能Ep2可以表示为:
不考虑面内变形时,应变为:
其中εx是X方向应变,εy是Y方向应变,γxy是切应变。
将式(4)代入到式(3)中,得到褶皱薄膜的应变势能为:
则褶皱薄膜的自由振动动能Ev为:
其中ρ为薄膜的密度。
S1.2)将振动变形函数wv和褶皱变形函数wz都表示成级数形式:
S1.3)将式(7)代入式(1)和(5)中,得到振动过程中总势能Ep的表达式为:
将式(7)代入式(6)中,则自由振动动能Ev变形为:
其中kijkl,kij,k′ijkl,mij均为中间变量,其表达式分别如下:
mij=ρhm∫∫wviwvjdxdy (13)
S1.4)根据拉格朗日方程
将式(8)和式(9)代入上式,得到非线性振动的微分方程为:
当不考虑上式中第三项时,得到非线性振动微分方程:
S1.5)根据步骤(S1.4)的非线性振动的微分方程,得到非线性振动的矩阵表达式为:
其中,M是质量矩阵,K是刚度矩阵。
S2:用有限元屈曲分析方法确定出褶皱薄膜的褶皱形态,更新有限元分析模型,并进行模态分析,得到步骤(S1)振动模型在引入褶皱信息后的振动特性;具体为:
S2.1)根据文献(李云良,谭惠丰,王晓华.矩形薄膜和充气管的屈曲及后屈曲行为分析.航空学,2008,29(4):886~892.)中给出的有限元屈曲分析方法,褶皱薄膜在一定载荷作用下的褶皱形态;
S2.2)将步骤(S2.1)所得到的褶皱形态更新到有限元分析模型中,得到褶皱变形下的应力分布以及包含褶皱信息的刚度矩阵;
从步骤(S1.2)可以看到,刚度矩阵K中包括了由褶皱变形引起的大变形刚度,可以将刚度矩阵K分为三部分,分别为初始应力刚度阵Kσ,线性刚度阵KL,和褶皱引起的大变形刚度阵KNL,具体表达式如下:
σ0为初始应力,BL和BNL分别为对应的应变矩阵,De为弹性矩阵。
S2.3)利用步骤(S2.2)所得到的褶皱变形下的应力分布以及包含褶皱信息的刚度矩阵,进行模态分析,扩展模态,得到褶皱薄膜结构的各阶固有频率和模态振型。
步骤(S2.2)相应的矩阵分量分别如式(11)和式(13)所示,带入式(17)中,对步骤(S1.5)非线性振动的矩阵表达式进行求解,即式(17)的解可以写成q=φejωt,得到下述方程:
(K-ω2M)φejωt=0 (18)
ω为考虑薄膜表面褶皱后的一组频率值,将其记为ω1,ω2,…,ωn,ωi表示第i阶固有频率,组成固有频率矩阵Ω为:
将模型转换到模态坐标下,节点位移的绝对坐标值q与模态坐标值qm之间存在关系式:
q=Φqm (22)
Φ代表模态振型矩阵,由各阶振型组成,每一阶振型中包含nd个节点,各阶振型组成的模态振型矩阵为:
其中φij表示第i阶模态第j个节点的面外位移。为了凸显褶皱对薄膜动态特性的影响,令VX=1.0mm,VY=0.25mm时,此时会在薄膜表面产生贯穿对角线的大幅褶皱,引入褶皱信息,更新有限元模型后进行模态分析,去除前几阶中固有频率为零的噪声振型后,提取前六阶模态振型如图1所示,从第三阶和第四阶振型图上也可以体现出褶皱的形态。提取相应的固有频率值,在表2中分别列出不考虑褶皱时与考虑褶皱时的前六阶固有频率,
表2 VX=1.0mm,VY=0.25mm时的固有频率(Hz)
模态 | 第一阶 | 第二阶 | 第三阶 | 第四阶 | 第五阶 | 第六阶 |
考虑褶皱 | 28.44 | 34.40 | 34.70 | 38.75 | 39.36 | 43.82 |
不考虑褶皱 | 18.96 | 27.36 | 28.98 | 32.99 | 33.17 | 33.26 |
通过对比发现,两组数据相差较大,说明了褶皱的存在确实改变了薄膜结构的动态特性,所以在进行动态分析时必须考虑到褶皱的存在。
S3:根据步骤(S2)中得到的振动特性,用模态叠加法得到以固有振型为基础的空间向量,建立不同载荷下褶皱薄膜的动力学模型;具体为:
S3.1)根据步骤(S2.3)的褶皱薄膜结构的各阶固有频率和模态振型以及非线性振动的矩阵表达式求解方程,得到第i阶模态质量Mmi为:
将绝对坐标下的动力学方程转换到模态坐标下的动力学方程:
质量阵、刚度阵和阻尼阵被转换为模态质量阵Mm、模态刚度阵Km和模态阻尼阵Dm,且满足如下关系式:
其中第i阶模态质量Mmi可以表示为:
Ei为第i阶模态的动能,ωi为第i阶模态的固有频率。
S3.2)根据比例阻尼得到的阻尼阵,即得到阻尼矩阵与质量矩阵和刚度矩阵的关系式为:
Dm=α1Km+α2Mm (27)
其中α1和α2是瑞丽阻尼系数,可以根据固有频率和阻尼比系数求出具体数值:
ω1和ω2分别为第一阶和第二阶的固有频率值,ξ1和ξ2为第一阶和第二阶的模态阻尼比系数。
S3.3)设定nd为有限元模型的自由度,输入自由度为nr,输出自由度为ns,则褶皱变形下的应力分布在广义坐标下用动力学方程表示为:
其中q是由各节点的面外位移向量组成的nd×1维的矩阵,是nd×1维的节点速度矩阵,是nd×1维的节点加速度矩阵,V=[VX VY]T是输入位移载荷矩阵,即输入自由度nr=2,M是nd×nd维的质量矩阵,D为nd×nd维的阻尼矩阵,K为nd×nd维的刚度矩阵,B为nd×nr维输入矩阵,C为1×nd维面外振动位移输出矩阵,y为输出变量。
S3.4)将步骤(S3.3)广义坐标下用动力学方程转换到以固有振型为基础的空间向量,建立不同载荷下褶皱薄膜的动力学模型,表达式为:
其中
其中Mm,Km,Dm,Cm,Bm分别为转换后的模态质量矩阵、模态刚度矩阵、模态阻尼矩阵、模态输出矩阵和模态输入矩阵。
S4:根据各阶模态对系统振动变形的贡献,对步骤(S3)不同载荷下褶皱薄膜的动力学模型进行缩聚,建立面向控制的薄膜结构动力学模型;具体为:根据有限元模型得到的式(29)中的振动模型维数等于节点个数,本发明中正方形结构有限元模型的单元个数为10000,如此高的维数使得计算量增大,对动态数值求解困难,更不利于控制器的设计,因此需要对全阶振动模型进行模态降阶,要求降阶后的缩聚模型能够保留全阶模型的固有特性。
S4.1)根据模态能量法对模型进行缩聚是一种常用于柔性结构模型降阶的简单有效的方法,模态能量由两部分组成:模态动能和模态势能,则第i阶模态能量可以表示为:
前G阶中各阶所占的模态能量比ηi为:
如以前15阶为例,前15阶中各阶所占的模态能量比ηi为:
S4.2)提取前15阶中各阶模态能量,计算出各阶能量占总能量的比率,画出相应的条形图,将模态能量比按大小顺序排列并相加,当模态能量之和达到总能量的90%时,即可取出相应的阶次,以截取前六阶模态为例,即步骤(S3)中式(21)变为:
步骤(S3)中式(23)变为:
截取步骤(S3.3)得到的不同载荷下褶皱薄膜的动力学模型中对应的数据,建立面向控制的薄膜结构动力学模型:
其中A′p,B′p,C′p对应为面向控制的薄膜结构动力学模型的状态矩阵、输入矩阵和输出矩阵。
当各拉力大小发生变化后,动态特性发生改变后,由此得到的动力学方程也会随之改变,因此需要对被控对象建立多个模型来逼近系统的动态性能。划分区间越小,模型越逼近实际模型,但是模型太多势必会增大计算复杂度,可根据具体的结构选择合适的子模型来近似逼近实际模型,本发明中划分13个子模型来近似逼近实际模型,得到面向控制的多模型薄膜结构动力学模型。
建立正方形薄膜结构的有限元模型,在薄膜表面一点施加垂直膜面方向的冲击扰动,扰动点位于距右侧顶点250mm的位置,扰动幅值为4mm,作用时间为0.02s,四个角施加相等大小的位移载荷,VX=VY=0.3mm,进行瞬态动力分析,仿真总时间为4s。当位移载荷VX=1.0mm,VY=0.25mm时,此时褶皱将贯穿对角线,褶皱形态如图2所示,在面外扰动不变的情况下,分别采用有限元数值仿真方法和缩聚模型理论分析方法,得到中心线上距离中心点25mm处P点的振动曲线,如图3所示,两种方法得到的振动幅值最大误差为0.025mm,说明缩聚后的降维模型保留了全阶模型的基本属性和固有特性,验证了理论模型的准确性。
S5:结合步骤(S4)的面向控制的薄膜结构动力学模型以及驱动装置的动力参数,建立驱动模型;具体为:
现驱动装置只要能够实现下述功能的电机均可,现选用音圈电机模组为例进行说明,其中ua为电机输入电压,ia为电枢电流,La为电枢电感,Ra为电机内阻,电机运动时产生的反电动势ea为:
ea=Ksva (36)
其中,Ks为电机力常数,va为电枢切割磁力线的速度。
运动状态下,电枢回路的电压平衡方程为:
电磁力克服动子部分的惯性力Fm可以表示为:
式中,ma为动子部分总质量。
令动摩擦系数为ka,电机输出位移为x,则电机的动态力平衡方程式为:
Fe=Fm+kax (39)
根据电机的动态方程以及电的动态方程,可以得到输入电压ua与电机输出位移x之间的传递函数:
结合步骤(S4.3)建立的面向控制的薄膜结构动力学模型以及驱动装置的动力参数,拟合出驱动模型:
S6:以步骤(S5)建立的驱动模型作为被控对象,设计自适应参考控制模型,得到自适应参考控制模型的输出电压,加载到被控对象上,实现性能指标最优条件下的前馈补偿;具体为:
S6.1)控制目标是实现系统的输出y能够一致性跟踪特定的参考输入yr,令系统的输出跟踪误差ey为:
ey=y-yr (42)
其中
S6.2)在扩展模型中将输出变量y作为一个状态变量,构建自适应参考控制模型,定义参考控制模型状态矩阵为:
可以借助线性二次型调节器得到最优状态反馈增益矩阵Kx,选取二次性能指标如下:
上式中t0为系统的初始时刻,tf为系统的终端时刻,Q,QT为半正定对称的代价权衡矩阵,R为正定对称的控制输入加权矩阵。
则参考模型中的最优状态反馈增益矩阵Kx为:
其中,矩阵Px满足条件:
因此,可以得到所要构建的自适应参考控制模型,表达式为:
S6.3)重写扩展模型(43)中的状态方程:
定义新系统的状态跟踪误差e为:
e=x-xm (51)
选取自适应律为:
选取全局二次无界Lyapunov函数:
V(e,ΔΘ)=eTPe+tr(ΔΘTΓΘ -1ΔΘ) (54)
对上式求导,结合式(50)可以得到:
因此,闭环系统的误差模型是一致稳定且有界的。
S7:根据PID控制模型以及步骤(S6)的自适应参考控制模型(MRAC),得到多模型切换的PID-MRAC控制模型,并对其输出电压进行限幅滤波,加载到被控对象上,实现褶皱薄膜的主动振动控制。
多模型切换的PID-MRAC控制模型,如图4所示,控制输入ua由两部分组成:
ua=λ1ua1+λ2ua2 (56)
其中,λ1和λ2为控制输入的权系数,ua1为PID控制模块的输入电压,ua2为PID控制模型的输出电压即MRAC控制模块的输入电压。
当面外振动变形量|y|≥0.5mm时,令λ1=0.8,λ2=0.2,此时PID控制模块起主要作用,能够保证施加在薄膜结构顶点上的位移快速增大,相应地增大薄膜刚度,减小振动幅值;当|y|<0.5mm时,令λ1=0,λ2=1,此时PID控制器不起作用,仅通过MRAC控制模块调节薄膜结构位移载荷大小,是为了保证位移载荷达到一定值时不再持续增大,以免拉力过大破坏薄膜结构。
为了验证本发明的基于面内调控的褶皱薄膜主动振动控制方法的有益效果,下面通过单次扰动实验和多次扰动的实验为例进行说明,具体如下:
1、单次扰动实验
张紧薄膜,调整拉索上的初始拉力为9.8N,约等于在薄膜两个方向上同时施加位移为VX=VY=0.3mm,通过单次撞击,在距右侧顶点250mm的位置施加面外扰动,具体扰动施加示意图见图1所示。提取对角线上各点的振动曲线,如图5所示,可以看到最大振动幅值点位于中心点处。通过PSD高精度位置传感器测得薄膜中心点处振动量如图6所示。
对比实验结果、理论结果和数值仿真结果发现,得到的振动趋势保持一致,从局部放大图也可以看到三种方法得到的结果相差不大,说明本发明的主动振动控制方法具有可行性。
当施加外界扰动后,分别采用本发明的PID-MRAC控制方法和传统的PID控制方法进行振动抑制,图7给出了两种方法得到的中心点处的振动曲线,与不加控制时的自由振动曲线进行对比,可以看到振动得到有效抑制。从图7中的局部放大图可以看到,在2s以后,采用PID-MRAC方法得到的振动幅值明显小于用传统PID方法得到的振动幅值。假设当振动幅值小于0.1mm时认为振动停止,振动持续的时间见表3,无控制时自由振动时间为12.37s,采用本章的控制方法将振动时间缩短了4倍多,仅用了2.96s便将振幅控制在0.1mm以下,而传统的PID控制方法振动持续时间为4.91s,相比两种控制方法,本发明的方法振动时间缩短了约39.7%,说明了PID-MRAC方法对抑制振动的有效性。
表3振动时间(s)
无控制 | PID-MRAC | PID | |
振动时间(s) | 12.37 | 2.96 | 4.91 |
与现有技术中在薄膜上粘贴作动器的方法进行比较,如图8所示,在靠近顶点处粘贴两个面外作动器,通过垂直膜面方法的运动达到抑制薄膜振动的目的,同样以第四章得到的动力学模型作为被控对象,施加相同的面外扰动,设计相应的自适应控制器,调整拉索上的初始拉力为9.8N后保持不变,通过仿真分析得到中心点处振动抑制曲线,如图9所示,与本发明得到的实验结果进行对比,可以看到本文提出的方法得到的结果振动幅值较小,在图中给出Y轴数值为±0.01mm的直线,从局部放大图中可以得到当振幅小于0.01mm时,通过在膜面粘贴作动器的方法得到的振动时间约为3.52s,本发明结果得到的振动时间为2.96s,振动缩短了约15.9%,进一步说明了本发明所提方法对抑制振动的有效性。
2、多次扰动实验
在第一次扰动达到平衡后,施加在薄膜结构上的位移载荷约为0.86mm,在此基础上,再次通过撞击施加外界扰动,图10给出了施加在薄膜结构顶点处的位移载荷曲线,可以看到一共进行了五次扰动,在第三次扰动以后,再次施加扰动时,位移载荷会增大以提高薄膜结构固有频率,但又回到了平衡位置,并不会持续增大拉坏薄膜,最终到达平衡位置时,施加的位移大小约为0.99mm,此时,薄膜的抗扰动能力较强,在受到外界扰动时通过很小的调节即可达到快速抑制振动的目的。
图11和图12显示了每次扰动后中心点处的面外振动曲线,通过施加控制后的振动曲线和无控制的自由振动曲线的对比可以看到,当出现扰动时,控制器开始作用,并且均在3s内将振动幅值控制在0.1mm以下,进一步说明采用本发明的控制方法后,面外振动得到了有效快速地抑制。
从图11中还可以看到,在施加第二次面外扰动时,其最大振动幅值约为1.1mm左右,当施加第三次扰动时,最大振动幅值约为0.85mm,即随着位移载荷增大,膜面刚度随着增大,同样的面外扰动引起的振动幅值随之减小。当位移载荷达到0.99mm时,由图12可以看到,再次施加扰动,最大振动幅值明显减小,约为0.76mm。所以,初始位移越大,即初始预应力越大,膜面刚度越大,当受到外界扰动时,面外振动幅值越小,薄膜结构的抗扰能力越强,但是预应力过大容易破坏薄膜结构。从上面的分析可以得到,最终的到达平衡位置时,施加的位移大小约为0.99mm,对应的薄膜结构固有频率约为36.47Hz,因此设置合理地初始预应力,使薄膜结构固有频率达到36.47Hz左右,可以提高薄膜结构的抗扰能力,在受到外界扰动时,引起的振动幅值较小,且通过很小的调节即可达到快速抑制振动的目的。
以上例举仅仅是对本发明的举例说明,并不构成对本发明的保护范围的限制,凡是与本发明相同或相似的设计均属于本发明的保护范围之内。
Claims (10)
1.一种基于面内调控的褶皱薄膜主动振动控制方法,其特征在于,包括以下步骤:
(1)引入褶皱变形函数,建立褶皱薄膜的振动模型;
(2)用有限元屈曲分析方法确定出褶皱薄膜的褶皱形态,更新有限元分析模型,并进行模态分析,得到步骤(1)振动模型在引入褶皱信息后的振动特性;
(3)根据步骤(2)中得到的振动特性,用模态叠加法得到以固有振型为基础的空间向量,建立不同载荷下褶皱薄膜的动力学模型;
(4)根据各阶模态对系统振动变形的贡献,对步骤(3)不同载荷下褶皱薄膜的动力学模型进行缩聚,建立面向控制的薄膜结构动力学模型;
(5)结合步骤(4)的面向控制的薄膜结构动力学模型以及驱动装置的动力参数,建立驱动模型;
(6)以步骤(5)建立的驱动模型作为被控对象,设计自适应参考控制模型,得到自适应参考控制模型的输出电压,加载到被控对象上,实现性能指标最优条件下的前馈补偿;
(7)根据PID控制模型以及步骤(6)的自适应参考控制模型,得到多模型切换的PID-MRAC控制模型,并对其输出电压进行限幅滤波,加载到被控对象上,实现褶皱薄膜的主动振动控制。
2.根据权利要求1所述的基于面内调控的褶皱薄膜主动振动控制方法,其特征在于,所述步骤(1)的具体方法为:
(1.1)仅考虑薄膜的面外振动,令振动变形函数为wv,褶皱变形函数为wz,得到振动过程中褶皱薄膜弯曲形变势能Ep1和褶皱薄膜的应变势能Ep2以及褶皱薄膜的自由振动动能Ev,分别表达为:
其中,εx是X方向应变,εy是Y方向应变,γxy是切应变,E是薄膜材料的杨氏模量,v是泊松比,hm是薄膜厚度;ρ为薄膜的密度;
(1.2)将振动变形函数wv和褶皱变形函数wz都表示成级数形式:
(1.3)将步骤(1.2)的振动变形函数wv和褶皱变形函数wz带入步骤(1.1)的振动过程中总势能Ep和自由振动动能Ev,则变形为:
其中kijkl,kij,k′ijkl,mij均为中间变量;
(1.4)结合步骤(1.3)的振动过程中总势能Ep和自由振动动能Ev,根据拉格朗日方程,得到非线性振动的微分方程:
(1.5)根据步骤(1.4)的非线性振动的微分方程,得到非线性振动的矩阵表达式为:
其中,M是质量矩阵,K是刚度矩阵。
3.根据权利要求1所述的基于面内调控的褶皱薄膜主动振动控制方法,其特征在于,所述步骤(2)的具体方法为:
2.1)采用有限元屈曲分析法得到褶皱薄膜在一定载荷作用下的褶皱形态;
2.2)将步骤2.1)所得到的褶皱形态更新到有限元分析模型中,得到褶皱变形下的应力分布以及包含褶皱信息的刚度矩阵;
2.3)利用步骤2.2)所得到的褶皱变形下的应力分布以及包含褶皱信息的刚度矩阵,进行模态分析,扩展模态,得到褶皱薄膜结构的各阶固有频率和模态振型。
4.根据权利要求3所述的基于面内调控的褶皱薄膜主动振动控制方法,其特征在于,步骤(2.2)中所述的刚度矩阵K包括初始应力刚度阵Kσ、线性刚度阵KL和褶皱引起的大变形刚度阵KNL。
6.根据权利要求5所述的基于面内调控的褶皱薄膜主动振动控制方法,其特征在于,所述步骤(3)具体为:
3.1)根据步骤2.3)的褶皱薄膜结构的各阶固有频率和模态振型以及非线性振动的矩阵表达式求解方程,得到第i阶模态质量Mmi为:
Ei为第i阶模态的动能,ωi为第i阶模态的固有频率;
3.2)根据比例阻尼得到的阻尼矩阵,即得到阻尼矩阵与质量矩阵和刚度矩阵的关系式为:
Dm=α1Km+α2Mm
其中α1和α2是瑞丽阻尼系数,且
其中ω1和ω2分别为第一阶和第二阶的固有频率值,ξ1和ξ2为第一阶和第二阶的模态阻尼比系数;
3.3)设定nd为有限元模型的自由度,输入自由度为nr,输出自由度为ns,则褶皱变形下的应力分布在广义坐标下用动力学方程表示为:
其中q是由各节点的面外位移向量组成的nd×1维的矩阵,是nd×1维的节点速度矩阵,是nd×1维的节点加速度矩阵,V=[VX VY]T是输入位移载荷矩阵,即输入自由度nr=2,M是nd×nd维的质量矩阵,D为nd×nd维的阻尼矩阵,K为nd×nd维的刚度矩阵,B为nd×nr维输入矩阵,C为1×nd维面外振动位移输出矩阵,y为输出变量;
3.4)将步骤3.3)广义坐标下用动力学方程转换到以固有振型为基础的空间向量,建立不同载荷下褶皱薄膜的动力学模型,表达式为:
其中:
其中Mm,Km,Dm,Cm,Bm分别为转换后的模态质量矩阵、模态刚度矩阵、模态阻尼矩阵、模态输出矩阵和模态输入矩阵。
9.根据权利要求8所述的基于面内调控的褶皱薄膜主动振动控制方法,其特征在于,步骤(6)具体为:
6.1)令驱动模型的输出跟踪误差ey为:ey=y-yr
yr为参考输入;
其中
6.2)在扩展模型中将输出变量y作为一个状态变量,构建自适应参考控制模型,表达为:
其中参考控制模型的状态矩阵为:
矩阵Px满足条件:
Px(tf)=QT
tf为系统的终端时刻,Q,QT为半正定对称的代价权衡矩阵,R为正定对称的控制输入加权矩阵;
6.3)根据步骤6.2)中构建的自适应参考控制模型,选取自适应控制律为:
进而得到,自适应参考控制模型的输出电压为:
10.根据权利要求9所述的基于面内调控的褶皱薄膜主动振动控制方法,其特征在于,步骤(7)具体为:建立PID控制模型,结合步骤6.2)中建立的自适应参考控制模型,得到多模型切换的PID-MRAC控制模型,其输出电压为:
ua=λ1ua1+λ2ua2
其中,ua2为PID控制模型的输出电压,λ1和λ2为控制输入的权系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810523112.3A CN108897213B (zh) | 2018-05-28 | 2018-05-28 | 一种基于面内调控的褶皱薄膜主动振动控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810523112.3A CN108897213B (zh) | 2018-05-28 | 2018-05-28 | 一种基于面内调控的褶皱薄膜主动振动控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108897213A CN108897213A (zh) | 2018-11-27 |
CN108897213B true CN108897213B (zh) | 2021-04-30 |
Family
ID=64343959
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810523112.3A Active CN108897213B (zh) | 2018-05-28 | 2018-05-28 | 一种基于面内调控的褶皱薄膜主动振动控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108897213B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109709796B (zh) * | 2018-12-20 | 2021-11-12 | 武汉科技大学 | 一种惯性轮自平衡方法 |
CN109933882B (zh) * | 2019-03-07 | 2022-12-13 | 中国科学院紫金山天文台 | 环形主镜自由振动模态像差表征方法 |
CN110852015B (zh) * | 2019-11-19 | 2021-11-02 | 南京航空航天大学 | 一种陶瓷基复合材料模态的非线性计算方法 |
CN111475980B (zh) * | 2020-04-09 | 2022-07-26 | 西北工业大学 | 一种综合作动器质量影响的薄壁件动力学参数获取方法 |
CN112069710B (zh) * | 2020-09-07 | 2022-08-26 | 电子科技大学 | 用于气体静压主轴自激振动的预测方法 |
CN112380756B (zh) * | 2020-12-11 | 2022-05-03 | 福耀玻璃工业集团股份有限公司 | 检测调光玻璃的调光膜褶皱的方法及系统 |
CN112799303B (zh) * | 2021-01-06 | 2022-06-10 | 西安电子科技大学 | 一种机械臂的h∞控制方法 |
CN112968650B (zh) * | 2021-03-23 | 2023-02-28 | 武汉理工大学 | 一种基于降阶模型振动反馈控制的电机降噪方法及装置 |
CN115017637B (zh) * | 2022-05-10 | 2023-03-28 | 西北工业大学 | 一种航天用张拉整体模块构件展开过程的动力学特性分析方法 |
CN114996882A (zh) * | 2022-07-18 | 2022-09-02 | 纽威数控装备(苏州)股份有限公司 | 一种数控机床变形协同控制方法 |
CN115357063B (zh) * | 2022-07-27 | 2023-06-16 | 苏州科技大学 | 一种基于参考轨迹强跟踪的伺服系统的振动抑制方法 |
CN115730478B (zh) * | 2022-11-01 | 2023-06-13 | 哈尔滨工业大学 | 一种空间张拉薄膜结构动态响应面型的反馈与估计方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102279566A (zh) * | 2011-06-07 | 2011-12-14 | 广西工学院 | 敷有主动约束层弹性空腔噪声特性一体化自适应控制系统 |
CN102929132A (zh) * | 2012-03-29 | 2013-02-13 | 南京航空航天大学 | 基于新型多片组合驱动器构型的振动非接触主动控制装置及方法 |
CN104536295A (zh) * | 2014-12-17 | 2015-04-22 | 河海大学常州校区 | 一种悬臂梁鲁棒自适应控制方法 |
CN105975649A (zh) * | 2016-03-30 | 2016-09-28 | 西安电子科技大学 | 一种方形太阳帆表面褶皱形态的建模方法 |
CN106200378A (zh) * | 2016-07-08 | 2016-12-07 | 西安工业大学 | 鲁棒非脆弱h∞结构振动控制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070154437A1 (en) * | 2002-09-03 | 2007-07-05 | Pingfan Rao | Peptide-tagged proteins and methods of making and using thereof |
CN107145073A (zh) * | 2017-06-09 | 2017-09-08 | 武汉理工大学 | 一种基于宏纤维复合材料的拉索附加弯矩减振系统 |
-
2018
- 2018-05-28 CN CN201810523112.3A patent/CN108897213B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102279566A (zh) * | 2011-06-07 | 2011-12-14 | 广西工学院 | 敷有主动约束层弹性空腔噪声特性一体化自适应控制系统 |
CN102929132A (zh) * | 2012-03-29 | 2013-02-13 | 南京航空航天大学 | 基于新型多片组合驱动器构型的振动非接触主动控制装置及方法 |
CN104536295A (zh) * | 2014-12-17 | 2015-04-22 | 河海大学常州校区 | 一种悬臂梁鲁棒自适应控制方法 |
CN105975649A (zh) * | 2016-03-30 | 2016-09-28 | 西安电子科技大学 | 一种方形太阳帆表面褶皱形态的建模方法 |
CN106200378A (zh) * | 2016-07-08 | 2016-12-07 | 西安工业大学 | 鲁棒非脆弱h∞结构振动控制方法 |
Non-Patent Citations (3)
Title |
---|
Analysis of wrinkled membrane structures based on a wrinkle-wave model;Mingjun Liu,Jin Huang,Yali Wang;《AIP Advances》;20171231;第1-11页 * |
MODELING AND CONTROL OF WRINKLED MEMBRANES;M.J. Liu, Y.L. Wang, J. Huang;《The Fifth Asia International Symposium on Mechatronics》;20151031;第1-6页 * |
Wrinkling reduction of membrane structure by trimming edges;Mingjun Liu,Jin Huang,Mingyue Liu;《AIP ADVANCES》;20171231;第1-11页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108897213A (zh) | 2018-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108897213B (zh) | 一种基于面内调控的褶皱薄膜主动振动控制方法 | |
Aridogan et al. | A review of active vibration and noise suppression of plate-like structures with piezoelectric transducers | |
Huang et al. | Wind-tunnel tests for active flutter control and closed-loop flutter identification | |
Wang et al. | Model-free control design for multi-input multi-output aeroelastic system subject to external disturbance | |
CN106321719A (zh) | 一种采用正负刚度并联的主被动复合隔振器 | |
Chandiramani | Active control of a piezo-composite rotating beam using coupled plant dynamics | |
Sun et al. | Robust coordinated control for large flexible spacecraft based on consensus theory | |
Jiang et al. | Decentralized robust vibration control of smart structures with parameter uncertainties | |
Zou et al. | Robust active suppression for body-freedom flutter of a flying-wing unmanned aerial vehicle | |
Wang et al. | A dual-loop active vibration control technology with an RBF-RLS adaptive algorithm | |
Zhang et al. | Disturbance rejection control with H∞ optimized observer for vibration suppression of piezoelectric smart structures | |
Jiang et al. | Active vibration control based on modal controller considering structure-actuator interaction | |
Ast et al. | Adaptronic vibration damping for machine tools | |
Li et al. | Multi-modal vibration control for all-clamped plate subjected to periodic disturbances by ESO-based frequency-shaped LQR | |
Qian et al. | New method of modeling uncertainty for robust flutter suppression | |
Li et al. | Suppression of thermal postbuckling and nonlinear panel flutter motions using piezoelectric actuators | |
Kim et al. | Active control of nonlinear panel flutter using aeroelastic modes and piezoelectric actuators | |
Tavakolpour et al. | Active vibration control of a rectangular flexible plate structure using high gain feedback regulator | |
Park et al. | Adaptive control of nonlinear free vibration of shallow shell using piezoelectric actuators | |
Yang et al. | A novel parallel multi-harmonic global multi-channel control algorithm for helicopter active vibration control | |
Qian et al. | Modeling and adaptive vibration control for variable-length strip rolling systems with actuator faults and output constraints | |
Kajiwara et al. | Design of shape and control system for smart structures with piezoelectric films | |
Lee et al. | Composite Adaptive Aeroelastic System Control and Multifold Regressor Integrals-Based Identification | |
Huang et al. | Optimal piezoelectric actuators and sensors configuration for vibration suppression of aircraft framework using particle swarm algorithm | |
Yu et al. | Iterative Learning Control Based on Neural Network and Its Application to Ni-Mn-Ga Alloy Actuator With Local Lipschitz Nonlinearity |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |