Nonlinear vibration modeling analysis of viscoelastic vibration isolation platform for a cryocooler of the coronagraph
-
摘要: 中国空间站巡天望远镜(CSST)中的星冕仪模块能够实现系外行星直接成像,然而在轨运行期间其运动模块会产生微振动扰动而导致图像质量下降。文章针对航天载荷常用隔振系统进行非线性建模分析,针对黏弹性材料的Kelvin-Voigt模型进行非线性扩展,通过三阶谐波平衡法对本构方程求解。以微分形式建立了非线性Kelvin-Voigt模型的振动运动方程,对类似的设计假设以及工程应用提供了理论依据;模型分析结果表明,扩展后的非线性实体模型有效补充了非线弹性的微振动隔振系统设计理论,有助于提高工程设计的数据精度。
-
关键词:
- 星冕仪 /
- 巡天望远镜 /
- 非线性振动 /
- Kelvin-Voigt模型
Abstract: Direct imaging of exoplanets can be achieved with the coronagraph module of the China Space Station Telescope (CSST). However, its motion module will generate micro-vibration disturbances during orbital operation, resulting in a decrease in the image quality of the coronagraph. In this paper, the nonlinear modeling analysis on the commonly-used vibration isolation system for spacecraft loads was carried out. The nonlinear expansion of Kelvin-Voigt model for viscoelastic materials was conducted. The constitutive equation was solved using the third-order harmonic balance method. The vibration equation of nonlinear Kelvin-Voigt model was constructed in differential forms, providing a theoretical foundation for analogous design assumptions and engineering applications. The model analysis show that the expanded nonlinear solid model effectively supplements the design theory of nonlinear elastic micro-vibration isolation system, which is conductive to enhance the engineering data accuracy.-
Keywords:
- coronagraph /
- CSST /
- nonlinear vibration /
- Kelvin-Voigt model
-
0. 引言
太阳系外生命及系外行星的探测一直都是国际天文学领域研究的热点之一[1-3]。星冕仪作为系外行星直接成像的关键部件之一[4],其稳定性要求极高,因为如需探测到距离观测仪器10 pc处的类地行星系统,则角分辨率需达0.1″,也就是说,星冕仪模块目标成像对比度将高达10-8及以上。提高星冕仪光学稳定性的方式主要有主动光学矫正、热控和微振动抑制3种。本文主要论及微振动抑制。
在轨航天载荷的微振动抑制技术主要有主动隔振、被动隔振、主被动混合隔振以及半主动隔振技术,其中,被动隔振技术因无需外部能源,在空间任务中能有效压缩载荷质量,且具有可靠性高、稳定性强的优点,故应用较为普遍。目前工程上常用的被动隔振系统物理模型通常为弹簧与阻尼元件的并联组合[5-7],其与黏弹性材料的Kelvin-Voigt本构模型一致,能够在一定程度上反映隔振器的动力学特性,有助于对隔振系统进行理论分析。然而,随着空间载荷的指向精度要求越来越高,对于隔振系统的要求愈发严苛,相关研究从传统的线性隔振系统纷纷转向非线性的振动设计领域[8],并大量使用黏弹性阻尼材料[9]以及记忆合金[10]等新型材料进行微振动抑制。例如:邹元杰等[11]设计了用于大型环形天线的黏弹性阻尼器和金属橡胶阻尼器,试验结果表明被动隔振方案下振动幅值分别下降了78.7%和66%;庞世伟等[12]针对高分辨率遥感卫星设计了黏弹性隔振方案,隔振器在轨工作正常,满足各项功能指标要求;Uchida等[13]研究了超弹性合金在航天器微振动隔离器中的应用,并研制了支柱式隔离器,结果表明其在室温下满足传递率要求。考虑到在诸如黏弹性材料的非线性领域中继续使用线性的振动模型可能会使得研究数据精度下降,因此需要拓展传统的线性物理模型以适应非线性的微振动隔振理论研究。
本文基于中国空间站巡天望远镜(CSST)星冕仪制冷机振动模块,对传统的黏弹性材料弹簧–阻尼并联模型进行非线性扩展,假设模型所含弹簧元件提供三阶多项式的弹性恢复力,通过改变外界激励条件,研究扩展后的模型在受迫振动时的动力学行为,旨在提高黏弹性隔振平台的设计精度,进一步理解空间制冷机隔振器的微振动特性,为愈发复杂的微振动隔振设计提供新思路。
1. 隔振器物理模型
通常,CSST星冕仪制冷机产生的微振动经连接螺栓传递至间隔结构,隔振器即安装于连接螺栓上(如图1所示)。该隔振系统由上平台及基座组成,制冷机刚性安装于上平台,上平台与基座中间由4组隔振单元连接,隔振单元垂直高度为H。假设平台以及制冷机是总质量为m的刚体,则4组隔振单元的质量远小于m,故不考虑隔振单元的质量。P-xyz为原点固定于上平台质心的定坐标系,即该坐标系不随上平台的移动而改变;B-xyz为实验室定参考坐标系。系统内所有坐标系及旋转角都遵循右手法则。
由于所采用的隔振器为硅橡胶材质,故作黏弹性材料假设,物理模型如图2所示,用单自由度运动方程来离散连续系统。
该线性固体黏弹性材料的本构方程为
$$ \sigma = E\varepsilon + \mu \frac{{{\text{d}}\varepsilon }}{{{\text{d}}t}} \text{,} $$ (1) 式中:σ为应力; E为弹簧的弹性模量;ε为应变;μ为阻尼的黏度参数;t为时间。
2. 硅橡胶隔振器准静态试验分析
硅橡胶具有良好的隔振、耐高温及耐低温等特性,因此常被用于空间载荷的隔振系统。使用试验机测试单个硅橡胶隔振垫在不同压缩载荷下的变形(如图3所示),其中,压缩载荷采用位移方法控制,直到位移达到-5 mm。测试得到的载荷–变形曲线如图4所示,可观察到非线性荷载位移特征。
橡胶材料属于一种高分子材料,其分子微观结构是由卷曲的长分子链铰接形成的结构。当受到外界载荷拉伸时,卷曲的长分子链逐渐伸直,从宏观上看,这种长分子链的伸直可以表现为尺寸高达几倍的大变形;当外界拉伸载荷消失后,伸直的长分子链又可以恢复原状[14]。因此,与传统圆柱螺旋弹簧相比,硅橡胶隔振器呈现出明显的非线性特征。
对实验数据进行多项式函数拟合,迭代次数13次,达到10-9的Chi-sqr的容差值,拟合收敛。其负载与形变的关系为y=a0+a1x+a2x2+a3x3+…,其中:a0=0.8,a1=-24.7,a2=177.8,a3=-329。
与线性振动方程相比,非线性的弹性力关系更适宜建立振动控制方程。
3. 非线性振动建模
考虑到材料模型的非线性,对式(1)加以扩展,以导出非线性阻尼的表达式。图5为一个扩展的非线性实体模型,其中弹簧是非线性的,故将其命名为非线性实体模型。非线性弹簧表示系统的刚度是非线性的,提供弹性力k1x+k2x2+k3x3,其中:k1是弹簧的线性刚度;k2是二次刚度系数;k3是立方刚度系数;x是质量块的动态位移。
图5中受迫系统的动态响应 x(t) 可视为周期振动。杨氏模量E可以被视为应变的非线性函数。对于图中的单自由度系统,式(1)可以转化为
$$ F = {k_1}x + {k_2}{x^2} + {k_3}{x^3} + \tau \frac{{\text{d}}}{{{\text{d}}t}}({k_1}x + {k_2}{x^2} + {k_3}{x^3}) \text{,} $$ (2) 式中:τ=μ/E为黏弹性松弛时间常数;x是描述质量位移的坐标。本构方程(2)是一个常微分方程,它给出了施加在质量上的黏弹性动力F(t)。
系统的运动微分方程为
$$ m\ddot x + F(t) = \tilde f\sin (\omega t + \phi ) ,$$ (3) 式中:m表示模块的质量;
$\tilde f$ 为激励力的振幅;ω为激励频率;ϕ为相位角。对于周期激励,通过谐波平衡法可得x的前三阶表达形式如下:
$$x(t) = {a_0} + {a_1}\sin (\omega t) + {a_2}\cos (2\omega t) + {a_3}\sin (3\omega t) + \cdots 。 $$ (4) 由于位移分量并非全部同相,所以不再可能对弹性力和阻尼力进行简单解释,故在此引入所有项同相的简化表达式。同时,式(2)的解有如下形式:
$$\begin{split} F(t) = & {f_0} + {f_{1{\text{s}}}}\sin (\omega t) + {f_{1{\text{c}}}}\cos (\omega t) + {f_{2{\text{s}}}}\sin (2\omega t) +\\ & {f_{2{\text{c}}}}\cos (2\omega t) + {f_{3{\text{s}}}}\sin (3\omega t) + {f_{3{\text{c}}}}\cos (3\omega t) + \cdots 。 \end{split} $$ (5) 式(4)和式(5)中,a0, a1, a2, f0, f1s, f1c, …皆为待确定系数。
忽略式(3)与式(4)中激励频率高于三次的亚谐波与超谐波,将式(4)和(5)代入式(2),则对于式(2)的零阶项,有代数方程
$$\begin{split} {f_0} = & {\text{ }}{k_1}{a_0} + \frac{1}{2}{k_2}\left( {2a_0^2 + a_1^2 + a_2^2 + a_3^2} \right) + \frac{3}{4}{k_3}\times \\ &\left( {\frac{4}{3}a_0^3 + 2{a_0}a_1^2 + 2{a_0}a_2^2 + 2{a_1}{a_2}{a_3} + 2{a_0}a_3^2 - a_1^2{a_2}} \right) 。\end{split} $$ (6) 与sin(ωt)相乘的项展开为
$$\begin{split} {f_{\text {1s}}} = & {\text{ }}{k_1}{a_1} + {k_2}\left( {2{a_0}{a_1} - {a_1}{a_2} + {a_2}{a_3}} \right) + \\ & \frac{3}{4}{k_3}( 4a_0^2{a_1} + a_1^3 - 4{a_0}{a_1}{a_2} + 2{a_1}a_2^2 -\\ & a_1^2{a_3} + 4{a_0}{a_2}{a_3} - a_2^2{a_3} + 2{a_1}a_3^2 ) \text{;} \end{split}$$ (7) 而与cos(ωt)相乘的项为
$$ \begin{split} {f_{\text {1c}}} = & \tau \omega [{k_1}{a_1} + {k_2}(2{a_0}{a_1} - {a_1}{a_2} + {a_2}{a_3}) + \\ & \frac{3}{4}{k_3}(4a_0^2{a_1} + a_1^3 - 4{a_0}{a_1}{a_2} + 2{a_1}a_2^2 - \\ & a_1^2{a_3} + 4{a_0}{a_2}{a_3} - a_2^2{a_3} + 2{a_1}a_3^2)]。 \end{split} $$ (8) 式(7)给出了系统的弹性响应,它是与x同相的机械力。式(8)给出了系统的阻尼力,与频率相关,相对于位移x具有相位差。同理可得出二阶与三阶对应的项:
$$ \begin{split} {f_{\text {2s}}} = & \tau \omega [{\text -} 2{k_1}{a_2} + {k_2}(a_1^2 - 4{a_0}{a_2} - 2{a_1}{a_3}) +\\ & 3{k_3}({a_0}a_1^2 - 2a_0^2{a_2} - a_1^2{a_2} - {{a_2^3} / 2} -\\ & 2{a_0}{a_1}{a_3} + {a_1}{a_2}{a_3} - {a_2}a_3^2)]; \\ {f_{\text {2c}}} = & {k_1}{a_2} + \frac{1}{2}{k_2}( {\text -} a_1^2 + 4{a_0}{a_2} + 2{a_1}{a_3}) +\\ & \frac{3}{4}{k_3}( {\text -} 2{a_0}a_1^2 + 4a_0^2{a_2} + 2a_1^2{a_2} +\\ & a_2^3 + 4{a_0}{a_1}{a_3} - 2{a_1}{a_2}{a_3} + 2{a_2}a_3^2); \\ {f_{\text {3s}}} = & {k_1}{a_3} + {k_2}({a_1}{a_2} + 2{a_0}{a_3}) +\\ & \frac{3}{4}{k_3}( {\text -} a_1^3/3 + 4{a_0}{a_1}{a_2} - {a_1}a_2^2 +\\ & 4a_0^2{a_3} + 2a_1^2{a_3} + 2a_2^2{a_3} + a_3^3);\\ {f_{\text {3c}}} = & \tau \omega [3{k_1}{a_3} + {k_2}(3{a_1}{a_2} + 6{a_0}{a_3}) +\\ & \frac{3}{4}{k_3}( {\text -} a_1^3 + 12{a_0}{a_1}{a_2} - 3{a_1}a_2^2 +\\ & 12a_0^2{a_3} + 6a_1^2{a_3} + 6a_2^2{a_3} + 3a_3^3)] 。 \end{split} $$ (9) 将式(6)~式(9)代入式(3),同时式(3)右端可展开为
$\tilde f\sin (\omega t + \phi ) = \tilde f[\sin (\omega t)\cos \phi + \cos (\omega t)\sin \phi ]$ ,则式(3)的零阶项如下:$$\begin{split} & {k_1}{a_0} + \frac{1}{2}{k_2}[2a_0^2 + a_1^2 + a_2^2 + a_3^2] + \frac{3}{4}{k_3}\times \\ & \left[ {\frac{4}{3}a_0^3 + 2{a_0}a_1^2 + 2{a_0}a_2^2 + 2{a_1}{a_2}{a_3} + 2{a_0}a_3^2 - a_1^2{a_2}} \right] = 0 。\end{split} $$ (10) $\sin (\omega t)$ 与$\cos (\omega t)$ 的控制方程为$$ \begin{split} & {\text{-}}m{a_1}{\omega ^2} + {f_{\text {1s}}} = \tilde f\cos \phi , \\ & {f_{\text {1c}}} = \tilde f\sin \phi 。 \end{split} $$ (11) 同理可得二阶、三阶项的控制方程,通过这6个控制方程,可找到未知参数
$ {a_0},\;{a_1},\;{a_2},\;{a_3} $ 的数值解。阻尼力可写为
$$ \begin{split} {f_{\text {1c}}}\cos (\omega t) & = \tau \omega [{k_1}{a_1} + {k_2}(2{a_0}{a_1} - {a_1}{a_2} + {a_2}{a_3}) +\\ & \frac{3}{4}{k_3}(4a_0^2{a_1} + a_1^3 - 4{a_0}{a_1}{a_2} + 2{a_1}a_2^2 - \\ & a_1^2{a_3} + 4{a_0}{a_2}{a_3} - a_2^2{a_3} + 2{a_1}a_3^2)]\cos (\omega t) = \\ & \tau \omega \left( {{k_1}\frac{{{\text{d}}x}}{{{\text{d}}t}} + {k_2}\frac{{{\text{d}}{x^2}}}{{{\text{d}}t}} + {k_3}\frac{{{\text{d}}{x^3}}}{{{\text{d}}t}}} \right) = \\ & \tau \omega [ {{k_1}\dot x(t) + 2x(t){k_2}\dot x(t) + 3x{{(t)}^2}{k_3}\dot x(t)} ] 。 \end{split} $$ (12) 弹性力可写为
$$ \begin{split} {f_0} + & {f_{1\text{s}}}\sin (\omega t) = {k_1}[ {{a_0} + {a_1}\sin (\omega t)} ] + \frac{1}{2}{k_2}[ 2a_0^2 + a_1^2 + \\ & a_2^2 + a_3^2 + (4{a_0}{a_1} - 2{a_1}{a_2} + 2{a_2}{a_3})\sin (\omega t) ] + \\ & \frac{3}{4}{k_3}[ \frac{4}{3}a_0^3 + 2{a_0}a_1^2 + 2{a_0}a_2^2 + 2{a_1}{a_2}{a_3} + 2{a_0}a_3^2 - a_1^2{a_2} + \\ & (4a_0^2{a_1} + a_1^3 - 4{a_0}{a_1}{a_2} + 2{a_1}a_2^2 - a_1^2{a_3} + 4{a_0}{a_2}{a_3} - \\ & a_2^2{a_3} + 2{a_1}a_3^2)\sin (\omega t) ] = {k_1}x + {k_2}{x^2} + {k_3}{x^3}\text{,} \end{split} $$ (13) 其中,刚度呈现出经典的线性、二次和三次项,与线性黏弹性中获得的动态刚度一致。事实上,目前的方法是将系统的非线性行为考虑在内的线性弹性系数的扩展。
$$ \frac{{{k_1}}}{m} = \omega _n^2 \text{,} \tau = \frac{{2\zeta }}{{{\omega _n}}}\text{,} $$ (14) 式中ζ是阻尼比。将式(12)和式(13)引入式(3),则有
$$\begin{split} & \ddot x + \left(\omega _n^2x + \frac{{{k_2}}}{m}{x^2} + \frac{{{k_3}}}{m}{x^3}\right) + \\ & \quad \frac{{2\zeta \omega }}{{{\omega _n}}}\left[ {{k_1}\dot x(t) + 2x(t){k_2}\dot x(t) + 3x{{(t)}^2}{k_3}\dot x(t)} \right] =\\&\quad \frac{{\tilde f}}{m}\sin (\omega t) 。\end{split} $$ (15) 式(15)是具有二次和三次非线性的单自由度非线性系统对黏弹性行为的推广。涉及二阶时间导数的项是惯性力;第二项表示恒定刚度(具有线性、二次和三次分量);第三项与一阶时间导数相关,表示非线性阻尼;等号右端为外部激励项。式(15)可直接进行数值求解,不必再利用谐波平衡法进行计算。
4. 模型分析
根据硅橡胶隔振器准静态试验分析(如图4所示)可知,线性刚度应用于振动模型时省略了更高阶的弹性分量,显然与实际不符,会导致线性拟合黏弹性材料的载荷–变形特性出现极大偏差(如图6所示),数据分析精度下降。
扩展后的非线性实体模型主要适用于采用黏弹性材料的微振动隔振系统。在黏弹性材料呈现较强的非线性弹性行为时,二次及更高阶次的刚度分量会对该模型的精度产生显著影响;然而更多的高阶项意味着求解时的计算量大幅上升,会对数值求解造成阻碍。因此,选用合适的阶数是平衡精度需求与计算量的关键。
隔振器组件的阻尼特性可通过对硅橡胶隔振垫施加循环应力得出。设计特定工装夹具,首先控制压缩位移到-1 mm,随后释放,再施加拉伸负载到1 mm,通过该位移控制方法循环 0.1 Hz 的准静态载荷,结果如图7所示。
当黏弹性阻尼材料受到交变应力时,其应变滞后于应力德尔塔相位角,从而产生滞后效应。单个循环中拉伸–压缩变形模式下阻尼材料单位体积的能量损失为ΔW=πγ2E1tan δ=πγ2E2,其中:E1为储能模量;E2为损耗模量;γ为应变幅度。
一个循环中阻尼材料每单位体积的总应变能为W=(γ2E1)/2,因此,黏弹性阻尼材料的损耗因子可定义为tan δ=ΔW/(2πW)=2ζ,根据该式即可得出测试材料的阻尼比ζ。
为分析该动力学模型的动力学行为,将式(15)用一阶方程组表示,即
$$\left\{ \begin{aligned} &{{\dot x}_1} = {x_2} \\ & {{\dot x}_2} = f\sin (\omega t) - \omega _n^2(x + {\beta _2}{x^2} + {\beta _3}{x^3}) -\\ & \quad\quad c\omega \dot x(1 + 2{\beta _2}x + 3{\beta _3}{x^2}) \end{aligned}\right.\quad {\text 。} $$ (16) 式中:β2=k2/k1;β3=k3/k1;c=2ζωnm;
$ f = {{\tilde f}/m} $ 。本文取ωn2=35, β2=4, β3=4, m=1, c=0.24, f=10,使用MatCont软件分析该系统的运动特性。
图8为隔振系统随外部激励幅值f变化的相轨迹曲线,可见系统在该条件下表现出复杂的运动行为。
图9为隔振系统随外部激励幅值f变化的幅频曲线,可见:激励幅值的增加导致系统的振幅随之增大。由图9可见特有的跳跃现象,即随着激励频率增大,隔振系统的振幅沿幅频曲线变化:当达到点2时,振幅则从点2跳跃至点3;若激励频率ω逐渐减小,振幅从点4开始沿曲线的下半分支变化至点5处;再减小ω,振幅则从点5跃至点6,然后沿曲线的上半分支向点1。因此受迫振幅在点2~3的振动是不稳定的,这说明外部激励幅值f对整个系统影响很大。在实际隔振系统中,应当尽量控制激励振幅的大小,从而保证系统的平稳运行或工作。
固定激励幅值,初始值取x1=0, x2=0,外部激励频率变化对系统的影响如图10所示,可见:随着外部激励频率的逐渐增大,全部相轨迹渐进趋于一条封闭曲线;而当激励频率继续增大时,轨迹线逐渐退化为以原点为中心的圆。这表明系统有着稳态的周期运动。
5. 结束语
本文以中国空间站巡天望远镜(CSST)星冕仪制冷机的被动隔振系统为研究对象,通过在标准线性固体材料中引入弹性力获得了单自由度的非线性振动模型,从而对Kelvin模型进行了扩展,并从中导出非线性阻尼;然后以微分形式获得了扩展后的几何非线性Kelvin模型振动运动方程。该模型能够补充现如今常用隔振系统的振动理论所涉及的非线性项内容,适用于非线弹性明显的结构假设。
下一步研究将考虑在理论分析的基础上进行更全面的实验以及优化设计,通过设计新的非线性振动试验装置来验证该模型的精确度。
-
-
[1] ZHU Y T, DOU J P, ZHANG X, et al. Portable adaptive optics for exoplanet imaging[J]. Research in Astronomy and Astrophysics, 2021, 21(4): 34-42
[2] DOU J, REN D. Phase quantization study of spatial light modulator for extreme high-contrast imaging[J]. The Astrophysical Journal, 2016, 832(1): 84 DOI: 10.3847/0004-637X/832/1/84
[3] ZHU J, ZHANG H, LIANG E S, et al. Exoplanets in the Antarctic sky: IV Dual-band photometry of variables found by the CSTAR-II commissioning survey at the North sky[J]. The Astronomical Journal, 2020, 159(4): 172 DOI: 10.3847/1538-3881/ab7449
[4] 窦江培, 朱永田, 任德清. 太阳系外行星的研究现状[J]. 自然杂志, 2014, 36(2): 124-128 DOU J P, ZHU Y T, REN D Q. Current research status of exoplanets[J]. Chinese Journal of Nature, 2014, 36(2): 124-128
[5] WEI H, XUN Y, CAI J, et al. Vibration reduction in a pulse tube-cooled high purity germanium detector[J]. Applied Sciences, 2022, 12(4): 1827 DOI: 10.3390/app12041827
[6] CHEN S, XUAN M, XIN J, et al. Design and experiment of dual micro-vibration isolation system for optical satellite flywheel[J]. International Journal of Mechanical Sciences, 2020, 179: 105592 DOI: 10.1016/j.ijmecsci.2020.105592
[7] XU C, XU Z D, HUANG X H, et al. Modeling and analysis of a viscoelastic micro-vibration isolation and mitigation platform for spacecraft[J]. Journal of Vibration and Control, 2018, 24(18): 4337-4352 DOI: 10.1177/1077546317724321
[8] AMABILI M. Derivation of nonlinear damping from viscoelasticity in case of nonlinear vibrations[J]. Nonlinear Dynamics, 2019, 97(3): 1785-1797 DOI: 10.1007/s11071-018-4312-0
[9] KWON S C, JO M S, KO D H, et al. Viscoelastic multilayered blade-type passive vibration isolation system for a spaceborne cryogenic cooler[J]. Cryogenics, 2020, 105: 102982 DOI: 10.1016/j.cryogenics.2019.102982
[10] 渠磊, 闫泽红, 饶智祥, 等. 形状记忆合金在航空航天领域的应用研究综述[J]. 航空动力学报, 2022, 38(10): 1-15 DOI: 10.13224/j.cnki.jasp.20220493 QU L, YAN Z H, RAO Z X, et al. Review on shape memory alloys’ application in field of aerospace[J]. Journal of Aerospace Power, 2022, 38(10): 1-15 DOI: 10.13224/j.cnki.jasp.20220493
[11] 邹元杰, 葛东明, 刘绍奎, 等. 航天器大型天线振动控制方案及其试验验证[J]. 航天器工程, 2018, 27(3): 135-139 ZOU Y J, GE D M, LIU S K, et al. Vibration control scheme and its experimental verification for large spacecraft antenna[J]. Spacecraft Engineering, 2018, 27(3): 135-139
[12] 庞世伟, 潘腾, 范立佳, 等. 一种微振动隔振设计与验证[J]. 强度与环境, 2016, 43(5): 17-23 DOI: 10.19447/j.cnki.11-1773/v.2016.05.003 PANG S W, PAN T, FAN L J, et al. Design and verification of micro-vibration isolator for a satellite[J]. Structure & Environment Engineering, 2016, 43(5): 17-23 DOI: 10.19447/j.cnki.11-1773/v.2016.05.003
[13] UCHIDA H, MIZUTANI T, YASUDA S, et al. Microvibration isolation system with superelastic alloy for space applications[C]∥AIAA SCITECH 2022 Forum, 2022: 1462. doi: 10.2514/6.2022-1462
[14] 张文正. 几类轴对称超弹性橡胶结构的有限变形分析[D]. 大连: 大连理工大学, 2013