清华大学 机械工程系, 北京 100084
收稿日期: 2015-07-02
基金项目: 国家科技重大专项(2015ZX04001002)
作者简介: 杜云松(1983-), 男, 博士研究生。
通讯作者: 李铁民, 副研究员, E-mail: litm@mail.tsinghua.edu.cn
摘要:作为一种高精度的微动机构, 柔性铰链机构已经在多个工业领域得到应用。该文针对平面柔性铰链机构, 对机构设计中柔度计算方法进行了研究。在分析了一般柔性铰链的柔度计算方法的基础上, 利用虚功原理, 推导出单个铰链的变形/受力相较于整个机构的变形/受力的关系; 利用矩阵分析方法, 结合串联机构和并联机构的单元与终端的变形/受力的互逆形式的传递规律, 分别获得了简洁的串并联柔性铰链机构的柔度计算公式。针对3种典型的串并联柔性铰链机构, 进行了柔度计算和不同机构参数与各方向柔度关系的分析, 并利用经典的有限元法进行了对比仿真实验。实验结果表明: 该方法的计算结果与有限元法最大相差7%, 平均相差3%。这表明该方法能够为平面柔性铰链机构的设计和结构优化提供参考。
关键词: 柔性铰链机构 柔度 矩阵分析法
Compliance calculation method for planar flexure-based mechanisms
DU Yunsong, LI Tiemin, JIANG Yao, ZHANG Jinglei
Department of Mechanical Engineering, Tsinghua University, Beijing 100084, China
Abstract:Flexure-based mechanisms are widely used in industry as highly precise micro-motion mechanisms. This paper presents a compliance calculation method for planar flexure-based mechanisms. The compliance of flexure hinges is used to relate the deformations to the loads on flexure member with the end-effector motion obtained using the virtual work principle. Then, a matrix method is used to derive concise compliance equations for serial and parallel flexure mechanisms. The compliance and the relationship between the compliance terms and variable geometric parameters are analyzed for three typical flexure-based mechanisms. The model results compare well with finite element method (FEM) predictions with a maximum difference of 7% and an average difference of 3%. This method provides theoretical and technical support for the design and optimization of flexure-based mechanisms.
Key words: flexure-based mechanismscompliancematrix method
柔性铰链是利用材料的变形使相互连接的刚体产生微小位移及转动的一种柔性单元[1]。若干个柔性铰链可以组成柔性铰链机构,被广泛应用于微电子制造、航空制造及生物工程等领域[2-4]。柔性铰链机构与传统的刚性机构相比,具有许多优点:1) 没有间隙,可提高运动精度; 2) 免于磨损,能够延长寿命; 3) 免于润滑,可避免污染; 4) 采用整体化设计和加工,可简化结构、减小体积和质量等[5]。柔度是柔性铰链机构的重要设计参数,它能够反映载荷与变形量之间的关系。合理的柔度模型能够减小建模误差,进而有利于进行机构的设计优化及精确运动控制。柔性铰链机构包括空间和平面两种结构形式,其中平面柔性铰链机构应用最为广泛。
柔性铰链机构的柔度计算方法有伪刚体模型法、卡式第二定理法、有限元法等[6-7]。伪刚体模型法将弹性杆(或铰链)模型等效为相应的刚性杆模型,用刚体中对机构进行分析与综合的方法,实现模型的简化。Li等[6]采用伪刚体模型法,建立了一种压电陶瓷驱动的XY微定位平台的柔度模型,用于评估其运动学、静力学、刚度、承载能力和动力学等性能,并通过有限元法和实验进行了验证。Pei等[8]提出了具有较大位移特性的杠杆式柔性机构的伪刚体建模方法,在对传统的建模方法进行优化的基础上,建立了相对简洁的伪刚体模型。但是,伪刚体模型法不能进行完整的柔度/刚度建模,只能建立工作方向上的柔度/刚度模型,应用范围受到限制[9]。卡式第二定理法以应变能为基础,将柔性结构的应变能表示为载荷的函数。Lobontiu等[10-11]采用卡式第二定理法,建立了圆弧型和导角型柔性铰链的柔度模型,并采用有限元法进行了对比分析,验证了该柔度模型的准确性和合理性。但是,卡式第二定理涉及到偏微分方程的计算,过程复杂,不利于进行机构的设计优化。有限元法通过变分方法,使得误差函数达到最小值并产生稳定解,是最准确的柔性铰链机构柔度建模方法[12]。Qin等[13]采用有限元法,建立了2自由度对称柔性铰链机构的柔度模型,进而计算出最大耦合比,并通过实验进行了验证。但是,有限元法不能建立机构柔度与机构几何参数之间的解析关系,只适合于形状和尺寸确定以后、 制造和装配之前的性能验证分析[12]。矩阵分析法是根据机构与组成单元之间力和变形量的转换关系计算机构柔度的一种方法。Lobontiu[14]利用矩阵分析法,得到了多种工况下的平面串联柔性铰链机构的柔度模型,并通过对菱形放大机构的分析进行了验证,但是计算过程过于复杂。因此,柔性铰链机构的柔度计算方法,机构几何参数与机构柔度、柔度计算误差之间的关系还有待进一步研究。
本文基于矩阵分析法,得到平面串联、并联柔性铰链机构的简洁柔度计算公式。针对典型柔性铰链机构,采用有限元法对本文得到的公式的准确度进行分析,验证了该柔度计算方法的合理性和有效性; 进而基于该柔度计算方法分析机构几何参数与机构柔度、柔度计算误差之间的关系,为平面柔性铰链机构的设计优化提供参考。
1 柔性铰链的柔度柔性铰链是柔性铰链机构最重要的组成部分。根据自由度的不同,可以分为单轴、两轴和多轴等[1]。其中,单轴柔性铰链主要应用于平面柔性铰链机构。根据截面形状的不同,单轴柔性铰链可以分为圆弧型、直梁型、椭圆型、抛物线型及双曲线型等。
下面以圆弧型柔性铰链为例,分析柔性铰链的柔度。图 1所示为圆弧型柔性铰链结构示意图,主要包括中间的柔性铰链和两端的联接端部(一般只认为中间切口部分为柔性铰链)。结构几何参数包括:铰链厚度w、 铰链半径R、 铰链宽度D、 端部高度H和端部宽度W。
图 1 圆弧型柔性铰链结构示意图 |
图选项 |
如图 2所示,圆弧型柔性铰链单元的一端自由,另一端固定。在柔性铰链自由端a处建立坐标系,x轴沿铰链的轴向指向固定端,y轴沿高度方向向上,z轴垂直纸面向外。载荷作用在柔性铰链自由端a处,一般情况下,包括6个分量:两个弯曲转矩My和Mz,两个剪切力Fy和Fz,一个轴向力Fx,一个扭转力矩Mx。单轴柔性铰链机构主要应用于平面运动,因此只考虑同一个平面内的载荷分量Mz、 Fy和Fx对柔性铰链单元的作用,其他载荷分量Mx、 My和Fz作用在该平面外,在分析平面柔性铰链机构时,可以忽略。
图 2 圆弧型柔性铰链单元受载示意图 |
图选项 |
平面单轴柔性铰链端部变形量可以表示为
$\sigma =CF.$ | (1) |
平面单轴柔性铰链的柔度矩阵可以表示为[1]
$C=\left[ \begin{matrix} {{C}_{x-{{F}_{x}}}} & 0 & 0 \\ 0 & {{C}_{y-{{F}_{y}}}} & {{C}_{y-{{M}_{z}}}} \\ 0 & {{C}_{\alpha -{{F}_{y}}}} & {{C}_{\alpha -{{M}_{z}}}} \\\end{matrix} \right].$ | (2) |
2 平面串联柔性铰链机构的柔度若干个柔性铰链与刚性联接端部相互串联,可以组成平面串联柔性铰链机构,如图 3所示。
图 3 平面串联柔性铰链机构 |
图选项 |
载荷作用点与位移输出点之间的位置关系包括3种情况:载荷作用点靠近机构固定端、位移输出点靠近机构固定端及两作用点位置重合。
如图 3所示,当载荷作用点位于靠近固定端的刚体k上的E点、 位移输出点在刚体m上的P点时,只有固定端和刚体k之间的柔性铰链产生弹性变形,机构的柔度计算只包括前k个柔性铰链。
根据虚功原理,第i个柔性铰链在自身局部坐标系下的变形量δii与在机构位移输出端产生的变形量δi之间的关系可以表示为
${{\delta }_{i}}={{J}_{{{d}_{i}}}}\delta _{i}^{i}\left( i=1,2,3,\cdots k \right).$ | (3) |
作用在机构位移输出端的载荷F与第i个柔性铰链在自身局部坐标系下的载荷Fii之间的关系可以表示为
$F_{i}^{i}={{J}_{{{f}_{i}}}}F\left( i=1,2,3,\cdots k \right).$ | (4) |
将柔性铰链在机构位移输出端产生的变形量进行叠加,得到机构的变形量为
$\delta =\sum\limits_{i=1}^{k}{{{J}_{{{d}_{i}}}}\delta _{i}^{i}=}\sum\limits_{i=1}^{k}{{{J}_{{{d}_{i}}}}C_{i}^{i}F_{i}^{i}=}\left( \sum\limits_{i=1}^{k}{{{J}_{{{d}_{i}}}}C_{i}^{i}{{J}_{_{{{f}_{i}}}}}} \right)F.$ | (5) |
由式(5)得到串联机构的柔度为
$C=\frac{\partial \delta }{\partial F}=\sum\limits_{i=1}^{k}{{{J}_{{{d}_{i}}}}C_{i}^{i}}{{J}_{_{{{f}_{i}}}}}={{J}_{d}}{{C}^{*}}{{J}_{f}}.$ | (6) |
$\begin{align} & {{J}_{d}}=\left[ {{J}_{{{d}_{1}}}}{{J}_{{{d}_{2}}}}{{J}_{{{d}_{3}}}}\cdots {{J}_{{{d}_{k}}}} \right], \\ & {{J}_{f}}={{\left[ {{J}_{{{f}_{1}}}}^{T}{{J}_{{{f}_{2}}}}^{T}{{J}_{{{f}_{3}}}}^{T}\cdots {{J}_{{{f}_{k}}}}^{T} \right]}^{T}}, \\ & {{C}^{*}}=diag\left( C_{1}^{1}C_{2}^{2}C_{3}^{3}\cdots C_{k}^{k} \right). \\ \end{align}$ |
当载荷作用点与位移输出点位置重合,都位于图 3中刚体k上的E点时,满足
$C=\sum\limits_{i=1}^{k}{{{J}_{{{d}_{i}}}}C_{i}^{i}J_{{{d}_{i}}}^{T}.}$ | (7) |
图 4 平面并联柔性铰链机构 |
图选项 |
相较于支链,终端平台具有更大的刚度,可以等效为刚体。因此,终端平台位移输出点P的变形量可以表示为
$\delta ={{J}_{{{p}_{i}}}}{{\delta }_{i}}\left( i=1,2,3,\cdots k \right).$ | (8) |
将作用在串联支链上的载荷进行叠加,得到作用于终端平台的载荷,
$F=\sum\limits_{i=1}^{n}{{{J}_{{{F}_{i}}}}{{F}_{i}}=\sum\limits_{i=1}^{n}{{{J}_{{{F}_{i}}}}{{K}_{i}}{{\delta }_{i}}=}\left( \sum\limits_{i=1}^{n}{{{J}_{{{F}_{i}}}}{{K}_{i}}J_{{{p}_{i}}}^{-1}} \right)\delta .}$ | (9) |
由式(9)得到并联机构的刚度,
$K=\frac{\partial F}{\partial \delta }=\sum\limits_{i=1}^{n}{{{J}_{{{F}_{i}}}}{{K}_{i}}J_{{{p}_{i}}}^{-1}}.$ | (10) |
$C={{K}^{-1}}=\sum\limits_{i=1}^{n}{{{J}_{{{P}_{i}}}}{{C}_{i}}J_{{{F}_{i}}}^{-1}=J_{P}^{*}{{C}^{*}}J_{F}^{*}.}$ | (11) |
4 算例分析不失一般性,采用圆弧型柔性铰链作为机构分析的基本组成单元,结构几何参数如图 1所示,其柔度矩阵通过有限元法计算得到。
4.1 平面串联柔性铰链机构如图 5所示,选取文[15]中的平面串联柔性铰链机构。机构参数L=35 mm,铰链参数R=5 mm,D=10 mm,铰链厚度w分别取1、 3、 5 mm,比较有限元仿真值与本文方法计算值的误差。
图 5 平面串联柔性铰链机构 |
图选项 |
当载荷作用点和位移输出点都位于n点时,结果如表 1所示。
表 1 两作用点位置重合时串联柔性铰链机构误差
w/mm | e/% | ||
Cxn-Fx | Cyn-Fy | Cαn-Mz | |
1 | 2.4 | 0.2 | 0.06 |
3 | 4.4 | 1.02 | 0.47 |
5 | 5.6 | 2.1 | 0.9 |
表选项
当载荷作用点在靠近固定端的m点,位移输出点在n点,或者位移输出点在靠近固定端的m点、 载荷作用点在n点,误差如表 2所示。
表 2 两作用点位置不同时串联柔性铰链机构误差
w/mm | e/% | ||
Cxn-Fx | Cyn-Fy | Cαn-Mz | |
1 | 4.1 | 0.22 | 0.07 |
3 | 5.9 | 1.4 | 0.64 |
5 | 7 | 3.2 | 1.3 |
表选项
从表 1和2可以看出,随着柔性铰链厚度w的增加,轴向柔度、转角柔度和剪切柔度的误差均逐渐增大。其中,轴向柔度误差最大,为7%。总体看来,有限元仿真值与本文方法计算值之间误差较小。
4.2 平面并联柔性铰链机构平行四边形柔性铰链机构是组成并联柔性铰链机构的一种常用结构形式,在实际应用中,经常呈对称结构分布,组成双平行四边形铰链机构。如图 6a所示,平行四边形铰链机构的终端平台在载荷Fx作用下,沿x轴运动时,不会发生转动,但是沿y轴方向会产生不必要的耦合运动误差ey,从而对机构的动力学分析和运动控制带来不利影响。如图 6b所示,两个相同的平行四边形铰链机构沿着驱动轴线呈对称分布组成双平行四边形铰链机构,该对称分布铰链机构可以完全消除机构的耦合运动误差,实现机构的输出解耦,并且可以看作是两个理想的移动副。
图 6 平面并联柔性铰链机构常用结构形式 |
图选项 |
下面分别采用本文方法和有限元法分析典型平面并联柔性铰链机构,包括平行四边形铰链机构和双平行四边形铰链机构。
如图 7所示,平行四边形柔性铰链机构包括两条串联支链,载荷作用点P和位移输出点E位置相同,都位于机构上端部[16]。
图 7 平行四边形柔性铰链机构 |
图选项 |
机构参数B=35 mm,L=54 mm,铰链参数 R=5 mm,D=10 mm,铰链厚度w分别取1、 3、 5 mm 时,比较有限元仿真与本文方法的误差,结果如表 3所示。可以看出,轴向柔度误差最大,转角柔度误差最小,最大误差为6.3%。
表 3 平行四边形柔性铰链机构误差
w/mm | e/% | ||
Cxn-Fx | Cyn-Fy | Cαn-Mz | |
1 | 3 | 2 | 0.6 |
3 | 4.7 | 2.6 | 1.2 |
5 | 6.3 | 3.5 | 2.3 |
表选项
如图 8所示,参考文[15]中的双平行四边形并联柔性铰链机构。该机构包括4条串联支链,左右对称,载荷作用点和位移输出点均位于n点,端部A、 B、 C和D固定约束。
图 8 双平行四边形柔性铰链机构 |
图选项 |
铰链参数R=5 mm,D=10 mm,w=1 mm。当并联机构参数M=30 mm,P=10 mm时,(L-2R)/R取不同值; 当M=30 mm,L=10 mm时,P/R取不同值; 当P=10 mm,L=10 mm时,(M-2R)/R取不同值。分析主要结构尺寸参数 P/R、 (L-2R)/R、 (M-2R)/R与机构柔度之间的关系,比较有限元仿真与本文方法的误差。
由式(7)计算得到并联机构4条串联支链的柔度为
${{C}_{1}}=\sum\limits_{j=1}^{2}{{{J}_{{{d}_{j}}}}C_{j}^{j}J_{{{d}_{j}}}^{T}.}$ | (12) |
$\begin{align} & {{J}_{{{d}_{1}}}}=\left[ \begin{matrix} -1 & 0 & -\left( w/2+R \right) \\ 0 & 1 & -\left( L+B+2R \right) \\ 0 & 0 & -1 \\\end{matrix} \right], \\ & {{J}_{{{d}_{2}}}}=\left[ \begin{matrix} -1 & 0 & -\left( w/2+R \right) \\ 0 & 1 & -B \\ 0 & 0 & -1 \\\end{matrix} \right]. \\ \end{align}$ |
${{C}_{2}}=\sum\limits_{j=3}^{4}{{{J}_{{{d}_{i}}}}C_{i}^{i}J_{{{d}_{j}}}^{T}.}$ | (13) |
$\begin{align} & {{J}_{{{d}_{3}}}}=\left[ \begin{matrix} -1 & 0 & -\left( 3w/2+3R+P \right) \\ 0 & 1 & -\left( L+B+2R \right) \\ 0 & 0 & -1 \\\end{matrix} \right], \\ & {{J}_{{{d}_{4}}}}=\left[ \begin{matrix} -1 & 0 & -\left( 3w/2+3R+P \right) \\ 0 & 1 & -B \\ 0 & 0 & -1 \\\end{matrix} \right]. \\ \end{align}$ |
${{C}_{3}}=\sum\limits_{j=5}^{6}{{{J}_{{{d}_{i}}}}C_{i}^{i}J_{{{d}_{j}}}^{T}.}$ | (14) |
$\begin{align} & {{J}_{{{d}_{5}}}}=\left[ \begin{matrix} -1 & 0 & 3w/2+3R+P \\ 0 & 1 & L+B+2R \\ 0 & 0 & 1 \\\end{matrix} \right], \\ & {{J}_{{{d}_{6}}}}=\left[ \begin{matrix} -1 & 0 & 3w/2+3R+P \\ 0 & 1 & -B \\ 0 & 0 & -1 \\\end{matrix} \right]. \\ \end{align}$ |
${{C}_{4}}=\sum\limits_{j=7}^{8}{{{J}_{{{d}_{i}}}}C_{i}^{i}J_{{{d}_{j}}}^{T}.}$ | (15) |
$\begin{align} & {{J}_{{{d}_{7}}}}=\left[ \begin{matrix} -1 & 0 & w/2+R \\ 0 & 1 & L+B+2R \\ 0 & 0 & 1 \\\end{matrix} \right], \\ & {{J}_{{{d}_{8}}}}=\left[ \begin{matrix} -1 & 0 & w/2+R \\ 0 & 1 & -B \\ 0 & 0 & -1 \\\end{matrix} \right]. \\ \end{align}$ |
$C={{K}^{-1}}={{\left( \sum\limits_{i=1}^{4}{C_{i}^{-1}} \right)}^{-1}}.$ | (16) |
图 9 本文方法计算值与参数P/R的关系 |
图选项 |
图 10 本文方法计算值与参数(L-2R)/R的关系 |
图选项 |
图 11 本文方法计算值与参数(M-2R)/R的关系 |
图选项 |
如图 12a所示,随着P/R比值增大,轴向柔度Cx -Fx和转角柔度Cα -Mz的有限元与本文方法的误差线性增加,剪切柔度Cy -Fy的误差最小,基本稳定在0.1%,轴向柔度x -Fx的误差最大,但小于3%。如图 12b所示,随着(L-2R)/R比值增大,轴向柔度、转角柔度和剪切柔度的误差都缓慢减小。剪切柔度的误差最小,轴向柔度的误差最大,但小于2.5%。如图 12c所示,随着(M-2R)/R比值增大,轴向柔度和转角柔度的误差都减小,而剪切柔度的误差变化不大。总体看来,有限元仿真值与本文方法计算值的误差较小。
图 12 各向柔度误差与参数P/R、 (L-2R)/R、 (M-2R)/R的关系 |
图选项 |
5 结 论本文针对平面柔性铰链机构,基于矩阵分析法计算机构的柔度,在此基础上分析了典型串联柔性铰链机构、平行四边形和双平行四边形并联柔性铰链机构的柔度和柔度计算误差:
1) 将串联柔性铰链的变形量进行叠加,得到串联柔性铰链机构的柔度计算公式; 将并联支链的载荷进行叠加,得到并联柔性铰链机构的柔度计算公式。
2) 采用有限元法进行对比分析,结果表明典型柔性铰链机构的柔度计算误差较小,验证了本文方法的合理性和有效性。
3) 本文方法能够用于分析铰链机构几何参数与建模误差、机构柔度之间的关系,为铰链机构的设计优化提供了参考。
参考文献
[1] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Lobontiu N. Compliant Mechanisms: Design of Flexure Hinges[M].Boca Raton, FL: CRC Press, 2010. |
[2] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Yao Q, Dong J, Ferreira P M. Design, analysis, fabrication and testing of a parallel-kinematic micropositioning XY stage[J].International Journal of Machine Tools and Manufacture,2007, 47(6): 946–961. |
[3] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Zubir M N M, Shirinzadeh B, Tian Y. Development of a novel flexure-based microgripper for high precision micro-object manipulation[J].Sensors and Actuators A: Physical,2009, 150(2): 257–266. |
[4] | Journal of Central South University(Science and Technology), 41(2):649-654.-->贠远, 徐青松, 李杨民. 并联微操作机器人技术及应用进展[J].机械工程学报,2009, 44(12): 12–23.YUN Yuan, XU Qingsong, LI Yangmin, et al. Survey on parallel manipulators with micro/nano manipulation technology and applications[J].Journal of Mechanical Engineering,2009, 44(12): 12–23.(in Chinese) |
[5] | Journal of Central South University(Science and Technology), 41(2):649-654.-->于靖军, 裴旭, 毕树生, 等. 柔性铰链机构设计方法的研究进展[J].机械工程学报,2010, 46(13): 2–13.YU Jingjun, PEI Xu, BI Shusheng, et al. State-of-arts of design method for flexure mechanisms[J].Journal of Mechanical Engineering,2010, 46(13): 2–13.(in Chinese) |
[6] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Li Y, Xu Q. A novel piezoactuated XY stage with parallel, decoupled, and stacked flexure structure for micro-/ nanopositioning[J].IEEE Transactions on Industrial Electronics,2011, 58(8): 3601–3615. |
[7] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Xiao S, Li Y. Development of a large working range flexure-based 3-DOF micro-parallel manipulator driven by electromagnetic actuators[J].International Journal of Precision Engineering and Manufacturing,2014, 15(4): 735–744. |
[8] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Pei X, Yu J, Zong G, et al. An effective pseudo-rigid-body method for beam-based compliant mechanisms[J].Precision Engineering,2010, 34(3): 634–639. |
[9] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Li Y, Xu Q. Design and analysis of a totally decoupled flexure-based XY parallel micromanipulator[J].IEEE Transactions on Robotics,2009, 25(3): 645–657. |
[10] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Lobontiu N, Paine J S N, Garcia E, et al. Design of symmetric conic-section flexure hinges based on closed-form compliance equations[J].Mechanism and Machine Theory,2002, 37(5): 477–498. |
[11] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Lobontiu N, Paine J S N, Garcia E, et al. Corner-filleted flexure hinges[J].Journal of Mechanical Design,2001, 123(3): 346–352. |
[12] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Pham H H, Chen I M. Stiffness modeling of flexure parallel mechanism[J].Precision Engineering,2005, 29(4): 467–478. |
[13] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Qin Y, Tian Y, Zhang D. Design and dynamic modeling of a 2-DOF decoupled flexure-based mechanism[J].Chinese Journal of Mechanical Engineering,2012, 25(4): 688–696. |
[14] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Lobontiu N. Compliance-based matrix method for modeling the quasi-static response of planar serial flexure-hinge mechanisms[J].Precision Engineering,2014, 38(3): 639–650. |
[15] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Qin Y, Shirinzadeh B, Zhang D, et al. Compliance modeling and analysis of statically indeterminate symmetric flexure structures[J].Precision Engineering,2013, 37(2): 415–424. |
[16] | Journal of Central South University(Science and Technology), 41(2):649-654.-->Jia X, Liu J, Tian Y, et al. Stiffness analysis of a compliant precision positioning stage[J].Robotica,2012, 30(6): 925–939. |