留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

二自由度快速控制反射镜系统固有频率优化设计

张玮钒 颜昌翔 高志良 王思宇 申箫 袁静 董有志

张玮钒, 颜昌翔, 高志良, 王思宇, 申箫, 袁静, 董有志. 二自由度快速控制反射镜系统固有频率优化设计[J]. 红外与激光工程, 2021, 50(6): 20200450. doi: 10.3788/IRLA20200450
引用本文: 张玮钒, 颜昌翔, 高志良, 王思宇, 申箫, 袁静, 董有志. 二自由度快速控制反射镜系统固有频率优化设计[J]. 红外与激光工程, 2021, 50(6): 20200450. doi: 10.3788/IRLA20200450
Zhang Weifan, Yan Changxiang, Gao Zhiliang, Wang Siyu, Shen Xiao, Yuan Jing, Dong Youzhi. Optimal design of natural frequency of two-degree-of-freedom fast steering mirror system[J]. Infrared and Laser Engineering, 2021, 50(6): 20200450. doi: 10.3788/IRLA20200450
Citation: Zhang Weifan, Yan Changxiang, Gao Zhiliang, Wang Siyu, Shen Xiao, Yuan Jing, Dong Youzhi. Optimal design of natural frequency of two-degree-of-freedom fast steering mirror system[J]. Infrared and Laser Engineering, 2021, 50(6): 20200450. doi: 10.3788/IRLA20200450

二自由度快速控制反射镜系统固有频率优化设计

doi: 10.3788/IRLA20200450
基金项目: 国家自然科学基金(61627819,61727818,61805235)
详细信息
    作者简介:

    张玮钒,男,硕士生,主要从事机械设计制造及其自动化、光机结构优化设计和结构拓扑优化设计方面的研究

  • 中图分类号: TH122

Optimal design of natural frequency of two-degree-of-freedom fast steering mirror system

  • 摘要: 在二自由度快速控制反射镜系统设计中,为提高系统的控制带宽,应尽量降低工作方向上的低阶固有频率,提高非工作方向上的高阶固有频率。该课题以某深切口柔性铰链快反镜系统作为研究对象,首先对系统前四阶固有频率的振型运动方向进行了分析,并针对传统刚度计算方法不适用于第三阶振型方向的问题,重新推导了第三阶振型方向上的刚度计算公式;其次,利用能量法和卡式第二定理对深切口柔性铰链上的工作刚度进行了推导,并进行了非线性拟合化简,得出的简化计算公式计算结果与有限元仿真结果误差不超过8.9%,证明了推导的铰链工作刚度理论公式的准确性;然后,将第三阶振型方向刚度计算公式和柔性铰链刚度计算公式代入固有频率计算公式,并进行有限元验证,结果表明理论公式计算结果与有限元仿真结果误差不超过1.7%,证明了新的三阶振型方向上的刚度计算公式的准确性。最后利用遗传算法,对系统前四阶固有频率进行了多目标优化设计,到达设计要求,所求出的优化结构较初始结构有明显优化,工作方向刚度减小19.04%,非工作方向刚度提高297.83%和77.09%。此外还对其进行了有限元仿真验证,结果证明一、二阶固有频率减小8.08%、5.40%,三、四阶固有频率提高了112.59%、16.80%。证明优化结构较初始结构有较大提高,能有效提高系统控制带宽。
  • 图  1  (a)快速控制反射镜总体结构;(b)柔性支撑系统结构

    Figure  1.  (a) Overall structure of fast steering mirror; (b) Flexible support system structure

    图  2  (a)深切口柔性铰链平面图;(b)深切口柔性铰链示意图

    Figure  2.  (a) Plan view of deep-cut flexure hinge; (b) Schematic diagram of deep-cut flexure hinge

    图  3  (a)第一阶振型;(b)第二阶振型;(c)第三阶振型;(d)第四阶振型

    Figure  3.  (a) First-order mode; (b) Second-order mode; (c) Third-order mode; (d) Fourth-order mode

    图  4  (a)一阶振型刚度简化图;(b)二阶振型刚度简化图;(c)三阶振型刚度简化图;(d)四阶振型刚度简化图

    Figure  4.  (a) Simplified diagram of first-order mode shape stiffness; (b) Simplified diagram of second-order mode shape stiffness; (c) Simplified diagram of third-order mode shape stiffness;(d) Simplified diagram of fourth-order mode shape stiffness

    图  5  三阶振型刚度简化模型

    Figure  5.  Simplified model of third-order modal stiffness

    图  6  (a)非线性拟合曲线图;(b)非线性拟合判定系数

    Figure  6.  (a) Non-linear fitting curve diagram; (b) Non-linear fitting determination coefficient

    图  7  深切口柔性铰链工作刚度有限元验证

    Figure  7.  Finite element verification of the working stiffness of the deep-cut flexure hinge

    图  8  第三阶振型仿真图

    Figure  8.  Third-order mode shape simulation diagram

    图  9  遗传算法优化流程图

    Figure  9.  Optimization flow chart of genetic algorithm

    图  10  (a)前视图;(b)右视图;(c)仰视图

    Figure  10.  (a) Front view; (b) Right view; (c) Bottom view

    图  11  (a)一阶模态;(b)二阶模态;(c)三阶模态;(d)四阶模态

    Figure  11.  (a) First-order mode; (b) Second-order mode; (c) Third-order mode; (d) Fourth-order mode

    表  1  深切口柔性铰链工作刚度理论计算与有限元仿真对比

    Table  1.   Comparison of theoretical calculation and finite element simulation of working stiffness of deep-cut flexure hinge

    a
    /mm
    b
    /mm
    t
    /mm
    w
    /mm
    Theoretical
    calculation results
    /N·m·rad1
    Finite element
    simulation results
    /N·m·rad1
    Error
    percentage
    1081825.379926.0417−2.54%
    1061833.839932.65133.64 %
    530.588.46008.25272.45%
    530.555.28755.14662.74%
    520.557.93127.29228.76%
    下载: 导出CSV

    表  2  快反镜系统三阶固有频率理论计算与有限元仿真对比

    Table  2.   Comparison of theoretical calculation and finite element simulation of three-boundary natural frequency of fast mirror system

    a
    /mm
    b
    /mm
    w
    /mm
    t
    /mm
    Mirror
    thickness: h
    /mm
    Theoretical
    calculation
    value of
    third-order
    mode/Hz
    Three-order
    mode finite
    element simulation
    value/Hz
    Error
    percentage
    43.550.823284.40286.250.61%
    1088115218.90220.17−0.58%
    658115434.29427.111.68%
    下载: 导出CSV

    表  3  快反镜系统初始结构参数

    Table  3.   Initial structure parameters of the fast steering mirror system

    a/mmb/mmt/mmw/mmL/mmD/mmh/mm
    108182510015
    下载: 导出CSV

    表  4  遗传算法优化结果

    Table  4.   Optimization results of genetic algorithm

    ParameterValue
    a/mm4.8
    b/mm4
    t/mm0.8
    G(x)0.0181
    下载: 导出CSV

    表  5  多目标优化结构与初始结构系统刚度对比

    Table  5.   Comparison of system stiffness between multi-objective optimized structure and initial structure

    a/mmb/mmt/mmK1/N·m·rad1K2/N·m·rad1K3/N·m·rad1K4/N·m1G(x)
    Initial value 10 8 1 50.04 50.04 844.6296 1.323e+06 0.0279
    Optimized value 4.8 4 0.8 40.51 40.51 3.3602e+03 2.343e+06 0.0181
    Optimization rate −19.04% −19.04% 297.83% 77.09% 41.58%
    下载: 导出CSV

    表  6  优化结构与初始结构四阶固有频率对比

    Table  6.   Comparison of the fourth-order natural frequency between the optimized structure and the initial structure

    First-order mode/HzSecond-order mode/HzThird-order mode/HzFourth-order mode/Hz
    Initial structure60.71560.814236.93549.97
    Optimized structure55.80757.533503.68642.38
    Optimization rate−8.08%−5.40%112.59%16.80%
    下载: 导出CSV
  • [1] 付锦江. 基于椭圆弧柔性铰链支撑的高性能快反镜优化设计[D]. 长春: 中国科学院长春光学精密机械与物理研究所, 2016.

    Fu Jinjiang. Design and optimization on high performance fast steering mirror with flexible platform composed by elliptical hinge[D]. Changchun: Changchun Institute of Optics, Fine Mechanics and Physics Chinese Academy of Sciences, 2016. (in Chinese)
    [2] 艾志伟, 快速反射镜结构分析与控制系统设计研究[D]. 成都: 中国科学院光电技术研究所, 2017.

    Ai Zhiwei. Structure analysis and research on the design of control system for fast steering mirror[D]. Chengdu: Institute of Optics and Electronics, Chinese Academy of Sciences, 2017. (in Chinese)
    [3] 徐宁. 基于柔性机构的快速反射镜研究[D]. 长春: 中国科学院长春光学精密机械与物理研究所, 2018.

    Xu Ning. Research on fast steering mirror based on compliant mechanism[D]. Changchun: Changchun Institute of Optics, Fine Mechanics and Physics Chinese Academy of Sciences, 2018. (in Chinese)
    [4] 吴松航, 董吉洪, 徐抒岩, 等. 快速反射镜椭圆弧柔性铰链多目标优化设计. 红外与激光工程. https://kns.cnki.net/kcms/detail/12.1261.TN.20201010.1621.004.html.

    Wu Songhang, Dong Jihong, Xu Shuyan, et al. Multi-objective optimal design of elliptic flexible hinge in fast steering mirror[J]. Infrared and Laser Engineering, 2021, 50(4): 20200286. (in Chinese)
    [5] Lu Qian, Huang Weiqing, Wang Yin, et al. Optimazation design of deep-notch elliptical flexure hinges [J]. Optics and Precision Engineering, 2015, 23(1): 206-215. (in Chinese) doi:  10.3788/OPE.20152301.0206
    [6] Nicolae Lobontiu, Jeffrey S N Paine, Edward O’Malley, et al. Parabolic and hyperbolic flexure hinges: flexibility, motion precision and stress characterization based on compliance closed-form equations [J]. Precision Engineering, 2002, 26(2): 183-192.
    [7] Ni Yingxue, San Xiaogang, Gao Shijie, et al. Research on flexibility of the novel hybrid flexure hinge [J]. Infrared and Laser Engineering, 2016, 45(10): 1017001. (in Chinese) doi:  10.3788/IRLA201645.1017001
    [8] Fu Jinjiang, Yan Changxiang, Liu Wei, et al. Stiffness calculation and optimal design of elliptical flexure hinges [J]. Optics and Precision Engineering, 2016, 24(7): 1704-1709. (in Chinese)
    [9] Daniele Giannini, Francesco Braghin, Niels Aage. Topology optimization of 2D in-plane single mass MEMS gyroscopes [J]. Structural and Multidisciplinary Optimization, 2020(62): 2069-2089.
    [10] Zhao Lei, Ji Ming, Wang Jia, et al. Design and simulation of fast steering mirrors connected by universal flexure hinges [J]. Infrared and Laser Engineering, 2019, 48(2): 0218002. (in Chinese) doi:  10.3788/IRLA201948.0218002
    [11] 郑树泉. 工业智能技术与应用[M]. 上海: 上海科学技术出版社, 2019: 250-251.

    Zheng Shuquan. Industrial Intelligent Technology and Application[M]. Shanghai: Shanghai Science and Technology Press, 2019: 250-251. (in Chinese)
    [12] 陶经辉, 王声绪, 徐标. 集装箱港口智能化管理优化策略[M]: 中国财富出版社, 2016: 61-63.

    Tao Jinghui, Wang Shengxu, Xu Biao. Optimization Strategy for Intelligent Management of Container Ports[M]. Beijing: China Fortune Press, 2016: 61-63. (in Chinese)
    [13] 王铁方. 计算机基因学 基于家族基因的网格信任模型[M]. 北京: 知识产权出版社, 2016: 93-94.

    Wang Tiefang. Computer Genetics Grid Trust Model Based on Family Genes[M]. Beijing: Intellectual Property Press, 2016: 93-94. (in Chinese)
  • [1] 吴松航, 董吉洪, 徐抒岩, 于夫男, 许博谦.  快速反射镜椭圆弧柔性铰链多目标优化设计 . 红外与激光工程, 2021, 50(4): 20200286-1-20200286-9. doi: 10.3788/IRLA20200286
    [2] 艾志伟, 嵇建波, 王鹏举, 李静, 周皓阳.  两轴柔性支承快速反射镜结构控制一体化设计 . 红外与激光工程, 2020, 49(7): 20190479-1-20190479-8. doi: 10.3788/IRLA20190479
    [3] 费国标, 金勇, 陈新梅.  红外枪瞄减振系统优化设计 . 红外与激光工程, 2020, 49(S1): 20200188-20200188. doi: 10.3788/IRLA20200188
    [4] 汪奎, 辛宏伟, 曹乃亮, 石震.  空间相机快速反射镜的两轴柔性支撑结构设计 . 红外与激光工程, 2019, 48(12): 1214005-1214005(8). doi: 10.3788/IRLA201948.1214005
    [5] 赵磊, 纪明, 王佳, 赵振海, 王虎.  万向柔性铰链连接快速反射镜的设计与仿真 . 红外与激光工程, 2019, 48(2): 218002-0218002(7). doi: 10.3788/IRLA201948.0218002
    [6] 王凯迪, 苏秀琴, 李哲, 吴少博.  像移补偿快速反射镜时频特性优化控制 . 红外与激光工程, 2018, 47(S1): 112-118. doi: 10.3788/IRLA201847.S120003
    [7] 胡亮亮, 米凤文, 金伟其, 盛一成, 何玉婷, 雷琼莹.  基于PI逆模型的快速微摆反射镜的开环控制 . 红外与激光工程, 2017, 46(8): 818001-0818001(7). doi: 10.3788/IRLA201746.0818001
    [8] 吕世良, 刘金国, 周怀得, 梅贵.  星载红外探测器快速反射镜控制系统设计 . 红外与激光工程, 2017, 46(9): 904005-0904005(6). doi: 10.3788/IRLA201746.0904005
    [9] 李林, 王栋, 谭陆洋, 孔林, 杨洪波.  微小卫星天线与飞轮共支撑结构优化设计与试验 . 红外与激光工程, 2016, 45(10): 1018008-1018008(7). doi: 10.3788/IRLA201645.1018008
    [10] 魏传新, 陈洪达, 尹达一.  基于响应面法的交叉簧片铰链微位移机构优化设计 . 红外与激光工程, 2016, 45(10): 1018005-1018005(9). doi: 10.3788/IRLA201645.1018005
    [11] 吴伟彬, 戴一帆, 关朝亮, 范占斌, 钟曜宇.  横向压电效应变形镜优化设计 . 红外与激光工程, 2016, 45(8): 818003-0818003(7). doi: 10.3788/IRLA201645.0818003
    [12] 王夏霄, 王野, 秦祎, 于佳.  光纤电流互感器平均波长漂移变比误差 . 红外与激光工程, 2015, 44(1): 233-238.
    [13] 范磊, 王志, 曹玉岩.  基于力矩主动校正的反射镜支撑分析 . 红外与激光工程, 2015, 44(4): 1273-1277.
    [14] 潘雪涛, 高晓俭, 谷牧, 孟飞, 蔡建文.  接触线几何参数CCD交汇测量分辨率及精度分析 . 红外与激光工程, 2014, 43(11): 3627-3632.
    [15] 陈卫宁, 杨洪涛, 刘伟, 范哲源, 张兆会, 周祚峰, 刘广森, 雷杨杰, 祝青.  变焦凸轮曲线的优化设计方法 . 红外与激光工程, 2014, 43(5): 1535-1539.
    [16] 郭鹏, 张景旭, 杨飞, 张岩, 矫威.  2 m 望远镜K 镜支撑结构优化设计 . 红外与激光工程, 2014, 43(6): 1914-1919.
    [17] 张明月, 杨洪波.  基于大位移柔性结构的动镜支撑机构设计及研究 . 红外与激光工程, 2014, 43(8): 2595-2601.
    [18] 高晓丹.  宽波段7:3 消偏振分光膜的设计及镀制 . 红外与激光工程, 2013, 42(5): 1302-1305.
    [19] 刘齐民, 阮萍, 李福, 潘海俊.  空间光谱仪光栅柔性支撑设计与分析 . 红外与激光工程, 2013, 42(9): 2457-2461.
    [20] 郭万存, 吴清文, 杨近松, 黄勇, 高志良.  2m主镜主动支撑优化设计 . 红外与激光工程, 2013, 42(6): 1480-1484.
  • 加载中
图(11) / 表(6)
计量
  • 文章访问数:  34
  • HTML全文浏览量:  13
  • PDF下载量:  2
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-11-16
  • 修回日期:  2021-03-03
  • 刊出日期:  2021-06-30

二自由度快速控制反射镜系统固有频率优化设计

doi: 10.3788/IRLA20200450
    作者简介:

    张玮钒,男,硕士生,主要从事机械设计制造及其自动化、光机结构优化设计和结构拓扑优化设计方面的研究

基金项目:  国家自然科学基金(61627819,61727818,61805235)
  • 中图分类号: TH122

摘要: 在二自由度快速控制反射镜系统设计中,为提高系统的控制带宽,应尽量降低工作方向上的低阶固有频率,提高非工作方向上的高阶固有频率。该课题以某深切口柔性铰链快反镜系统作为研究对象,首先对系统前四阶固有频率的振型运动方向进行了分析,并针对传统刚度计算方法不适用于第三阶振型方向的问题,重新推导了第三阶振型方向上的刚度计算公式;其次,利用能量法和卡式第二定理对深切口柔性铰链上的工作刚度进行了推导,并进行了非线性拟合化简,得出的简化计算公式计算结果与有限元仿真结果误差不超过8.9%,证明了推导的铰链工作刚度理论公式的准确性;然后,将第三阶振型方向刚度计算公式和柔性铰链刚度计算公式代入固有频率计算公式,并进行有限元验证,结果表明理论公式计算结果与有限元仿真结果误差不超过1.7%,证明了新的三阶振型方向上的刚度计算公式的准确性。最后利用遗传算法,对系统前四阶固有频率进行了多目标优化设计,到达设计要求,所求出的优化结构较初始结构有明显优化,工作方向刚度减小19.04%,非工作方向刚度提高297.83%和77.09%。此外还对其进行了有限元仿真验证,结果证明一、二阶固有频率减小8.08%、5.40%,三、四阶固有频率提高了112.59%、16.80%。证明优化结构较初始结构有较大提高,能有效提高系统控制带宽。

English Abstract

    • 快速控制反射镜(Fast Steering Mirror, FSM)是一种可以在光源和接收器之间控制光束方向的反射镜装置,具有体积小、功耗低、响应速度快,精度高等特点,最早应用于自适应光学系统中,现广泛应用在精密捕获、瞄准和跟踪、光束稳定、目标指向、空间激光通信等领域[1]。其中,FSM系统的控制带宽表征了系统对输入信号的跟踪能力,系统控制带宽越高,系统响应速度越快,反应时间越短,FSM系统的跟踪能力就越精确、及时[1]。而国内现有的大口径快反镜系统控制带宽普遍有待提高,所以大口径、高控制带宽的FSM系统是未来FSM系统研究的重点领域。提高系统控制带宽需要提高系统的高阶固有频率,减小系统低阶固有频率[1],而系统的固有频率与系统在该振型方向上的刚度有关,故有必要对系统的工作刚度和固有频率计算公式展开研究。

      国内外已有多名学者对快反镜系统工作刚度和固有频率计算公式开展过分析研究。付锦江等[1]对二自由度椭圆弧柔性支撑系统的前四阶振型、运动刚度和固有频率进行了分析,并通过反复计算得到了椭圆弧柔性支撑系统的优化设计结果;艾志伟等[2]对FSM系统进行建模分析,并讨论了系统机械结构与控制系统的关系;徐宁等[3]对快反镜系统固有频率计算公式进行了推导,并进行了有限元仿真测试。上述工作主要集中于对快反镜系统刚度组成以及固有频率的分析,并以固有频率计算公式得出各方向上的刚度设计范围,缺乏以提高系统控制带宽为目标的定量优化计算。故有必要利用优化算法对快反镜系统结构固有频率进行定量的多目标优化计算。

      对于单个柔性铰链的工作刚度以及其多目标优化计算,也有多位学者对其进行了研究。吴松航等[4]利用NSGA-II算法对单个椭圆的柔性铰链进行了减小z轴扭转刚度,提高xy轴扭转刚度的多目标优化计算;卢倩等[5]对单个深切口柔性铰链的刚度矩阵进行了推导,并进行了减小两个工作方向刚度、提高七个非工作方向刚度的多目标优化计算,优化率在10%~68%;以上工作对于单个柔性铰链的分析优化有很大的指导意义,但缺乏将柔性铰链代入整个系统进行优化分析。

      针对以上问题,文中对二自由度快反镜系统的固有频率计算进行了理论分析,并基于减小系统低阶固有频率,提高系统高阶固有频率提出了一种多目标优化函数。首先,文中对深切口柔性铰链组成的二自由度快反镜系统进行了有限元仿真分析,得出了前四阶固有频率的振型,并将其简化为弹簧的串并联模型,针对传统刚度分析不适用于第三阶振型运动方向的问题,对其振型方向刚度与单个柔性铰链刚度之间的关系进行了分析,得出了新的快反镜第三阶振型方向的刚度理论公式;然后,文中利用能量法以及卡式第二定理对深切口柔性铰链工作刚度进行了推导,并利用非线性拟合对其进行了简化,解决了传统深切口柔性铰链工作刚度计算公式过于复杂的问题;之后,文中结合以上的刚度计算公式,对快反镜固有频率计算公式进行了有限元仿真验证,证明了三阶固有频率计算公式的准确性。最后,文中提出了一种减小系统低阶固有频率,提高系统高阶固有频率的多目标优化函数,并对四个不同的目标赋予不同的权值,利用遗传算法对其进行优化计算,并进行了有限元仿真验证,结果证明,优化计算有效地减小了系统低阶固有频率,提高了系统高阶固有频率,设计结果良好,理论计算与仿真结果吻合良好。文中提出的基于系统固有频率的多目标优化函数以及推导简化的深切口柔性铰链工作刚度计算公式对于进一步提高快反镜系统控制带宽、优化快反镜结构设计具有一定意义。

    • 文中讨论的快速控制反射镜以及柔性支撑系统三维结构模型如图1所示。

      图  1  (a)快速控制反射镜总体结构;(b)柔性支撑系统结构

      Figure 1.  (a) Overall structure of fast steering mirror; (b) Flexible support system structure

      图1(a)中柔性支撑结构与隔离板之间采用螺栓连接;隔离板与反射镜之间采用胶粘连接。图1(b)中的坐标系为O-xyzO为柔性支撑系统转动中心,工作方向为绕xy轴的转动。该结构由四个结构完全一样的柔性铰链组成,其中铰链1、3组合,可以实现结构在x轴上的偏转;铰链2、4组合,可以实现结构在y轴上的偏转[1]。单个深切口柔性铰链模型如图2所示。

      图  2  (a)深切口柔性铰链平面图;(b)深切口柔性铰链示意图

      Figure 2.  (a) Plan view of deep-cut flexure hinge; (b) Schematic diagram of deep-cut flexure hinge

      图2 中的坐标系为${{{O}}'} - {x'}{y'}{{\textit{z}}'}$,该坐标系适用于单个深切口柔性铰链。图中,h为柔性铰链的高度,w为柔性铰链的宽度,t为柔性铰链最小切割厚度。同时椭圆切口的长半轴和短半轴分别为ab

      对FSM系统进行模态分析,结果如图3所示。

      图  3  (a)第一阶振型;(b)第二阶振型;(c)第三阶振型;(d)第四阶振型

      Figure 3.  (a) First-order mode; (b) Second-order mode; (c) Third-order mode; (d) Fourth-order mode

      图3可知,FSM系统前四阶模态振型为绕xy轴偏转,绕z轴扭转以及沿z轴的拉压。所以可以将二自由度柔性支撑系统简化为弹簧串并联模型,并进行振型方向上的刚度计算[1]

      (1)第一阶振型:绕x轴转动

      x轴的转动刚度是由铰链1、3的刚度${K_{{\theta _{{z'}}},{M_{{z'}}}}}$并联,铰链 2,4的刚度${K_{{\theta _{{y'}}},{M_{{y'}}}}}$并联,再由二者串联而成[1],简化模型如图4(a)所示,计算公式如下:

      图  4  (a)一阶振型刚度简化图;(b)二阶振型刚度简化图;(c)三阶振型刚度简化图;(d)四阶振型刚度简化图

      Figure 4.  (a) Simplified diagram of first-order mode shape stiffness; (b) Simplified diagram of second-order mode shape stiffness; (c) Simplified diagram of third-order mode shape stiffness;(d) Simplified diagram of fourth-order mode shape stiffness

      $$\frac{1}{{{K_x}}} = \frac{1}{{{K_{1,{\theta _{{{\textit{z}}'}}},{M_{{{\textit{z}}'}}}}} + {K_{3,{\theta _{{{\textit{z}}'}}},{M_{{{\textit{z}}'}}}}}}} + \frac{1}{{{K_{2,{\theta _{{y'}}},{M_{{y'}}}}} + {K_{4,{\theta _{{y'}}},{M_{{y'}}}}}}}$$ (1)

      计算结果为:

      $${K_x} = \frac{{({K_{1,{\theta _{{z'}}},{M_{{z'}}}}} + {K_{3,{\theta _{{z'}}},{M_{{z'}}}}})({K_{2,{\theta _{{y'}}},{M_{{y'}}}}} + {K_{4,{\theta _{{y'}}},{M_{{y'}}}}})}}{{({K_{1,{\theta _{{z'}}},{M_{{z'}}}}} + {K_{3,{\theta _{{z'}}},{M_{{z'}}}}} + {K_{2,{\theta _{{y'}}},{M_{{y'}}}}} + {K_{4,{\theta _{{y'}}},{M_{{y'}}}}})}}$$ (1)

      (2)第二阶振型:绕y轴转动

      y轴的转动刚度是由铰链2、4的刚度${K_{{\theta _{{{\textit{z}}'}}},{M_{{{\textit{z}}'}}}}}$并联,铰链 1、3的刚度${K_{{\theta _{{y'}}},{M_{{y'}}}}}$并联,再由二者串联而成[1],简化模型如图4(b)所示,计算公式如下:

      $$\frac{1}{{{K_y}}} = \frac{1}{{{K_{2,{\theta _{{{\textit{z}}'}}},{M_{{{\textit{z}}'}}}}} + {K_{4,{\theta _{{{\textit{z}}'}}},{M_{{{\textit{z}}'}}}}}}} + \frac{1}{{{K_{1,{\theta _{{y'}}},{M_{{y'}}}}} + {K_{3,{\theta _{{y'}}},{M_{{y'}}}}}}}$$ (2)

      计算结果为:

      $${K_y} = \frac{{({K_{2,{\theta _{{{\textit{z}}'}}},{M_{{{\textit{z}}'}}}}} + {K_{4,{\theta _{{{\textit{z}}'}}},{M_{{{\textit{z}}'}}}}})({K_{1,{\theta _{{y'}}},{M_{{y'}}}}} + {K_{3,{\theta _{{y'}}},{M_{{y'}}}}})}}{{({K_{2,{\theta _{{{\textit{z}}'}}},{M_{{{\textit{z}}'}}}}} + {K_{4,{\theta _{{{\textit{z}}'}}},{M_{{{\textit{z}}'}}}}} + {K_{1,{\theta _{{y'}}},{M_{{y'}}}}} + {K_{3,{\theta _{{y'}}},{M_{{y'}}}}})}}$$ (3)

      (3)第三阶振型:绕z轴扭转

      常用的二自由度快反镜系统三阶振型计算公式为:

      $${K_{{\theta _{\rm{{\textit{z}}}}}}} = 4L{K_{{\theta _{{{\textit{z}}'}}},{F_{{y'}}}}}$$ (4)

      在柔性支撑系统面对绕z轴的转矩${M_{\textit{z}}}$时,${M_{\textit{z}}}$分布在四个铰链上的剪切力为$F = \dfrac{{{M_{\textit{z}}}}}{L}$,其中,$L$为柔性铰链偏转中心到柔性支撑系统偏转中心的距离,即系统偏转半径。但是分析快反镜系统第三阶运动形式可知,${K_{{\theta _{\textit{z}}},{F_y}}}$只能表述柔性铰链在受切向力的情况下绕y轴偏转的角度,而与快反镜系统绕z轴的偏转不在一个平面上。其简化模型如图5所示。

      图  5  三阶振型刚度简化模型

      Figure 5.  Simplified model of third-order modal stiffness

      图5中,O点为柔性铰链偏转中心,b为柔性铰链短半轴长度,L为柔性铰链至支撑系统转动中心的长度,P为支撑系统转动中心,F为柔性铰链所受切向力。传统公式中的${K_{{\theta _z},{F_y}}}$只能表示在受到切向力时产生的偏转角$\theta $,而系统三阶振型产生的扭转角为$\delta $,故传统公式的$4L{K_{{\theta _z},{F_y}}}$无法表示系统扭转刚度。

      根据图5中的几何关系可以得出:

      $$b\tan \theta = L\tan \delta $$ (6)

      因为是小变形情况,有:

      $$\begin{array}{l} \tan \theta \approx \theta \\ \tan \delta \approx \delta \\ \end{array} $$ (5)

      故:

      $${K_{{\theta _{{{\textit{z}}'}}},{F_{{y'}}}}} = \frac{F}{\theta } = \frac{F}{{\left(\dfrac{L}{b} \cdot \delta \right)}} = \frac{b}{L} \cdot \frac{F}{\delta } = \frac{b}{L} \cdot {K_{{\delta _{\textit{z}}},{F_{{y'}}}}}$$ (8)

      式中:${K_{{\delta _z},{F_{{y'}}}}}$表示在单位切向力F的 作用下,单个柔性铰链围绕支撑系统偏转中心产生的偏转角。所以,绕z轴扭转刚度是由四个铰链的刚度${K_{{\delta _z},{F_{{y'}}}}}$并联组合而成,简化模型如图4(c)所示,故其计算公式为:

      $$ \begin{split} {{{K}}_{{{{\theta }}_{{z}}}}}{\rm{ = }}\dfrac{{{{{M}}_{{z}}}}}{{{{{\delta }}_{{z}}}}} {\rm{ = }}\frac{{{{FL}}}}{{{{{\delta }}_{{z}}}}} {{ = L(}}{{{K}}_{{\rm{1,}}{{{\delta }}_{{z}}},{{{F}}_{{{{y}}'}}}}}+ {{{K}}_{{\rm{2,}}{{{\delta }}_{{z}}},{{{F}}_{{{{y}}'}}}}}{\rm{ + }}{{{K}}_{{\rm{3,}}{{{\delta }}_{{z}}},{{{F}}_{{{{y}}'}}}}}{\rm{ + }}{{{K}}_{{\rm{4,}}{{{\delta }}_{{z}}},{{{F}}_{{{{y}}'}}}}}{\rm{)}} \\ \end{split} $$ (9)

      计算结果为:

      $${K_{{\theta _{\rm{{\textit{z}}}}}}} = 4L{K_{{\delta _{\textit{z}}},{F_y}'}} = 4\frac{{{L^2}}}{b}{K_{{\theta _{{{\textit{z}}'}}},{F_{{y'}}}}}$$ (10)

      (4)第四阶:沿z轴拉压刚度

      沿z轴的转动刚度是由铰链2、4的刚度${K_{{\delta _{{x'}}},{F_{{x'}}}}}$并联,铰链 1、3的${K_{{\delta _{{x'}}},{F_{{x'}}}}}$并联,再由二者串联而成[1],简化模型如图4(d)所示,计算公式如下:

      $$\frac{1}{{{K_{{\Delta _z}}}}} = \frac{1}{{{K_{1,{\delta _{{x'}}},{F_{{x'}}}}} + {K_{3,{\delta _{{x'}}},{F_{{x'}}}}}}} + \frac{1}{{{K_{2,{\delta _{{x'}}},{F_{{x'}}}}} + {K_{4,{\delta _{{x'}}},{F_{{x'}}}}}}}$$ (11)

      计算结果为:

      $${K_{{\Delta _{\textit{z}}}}} = {K_{{\delta _{{x'}}},{F_{{x'}}}}}$$ (12)
    • 定义柔性铰链上的力和弯矩为$F = [ {{F_{{x'}}}} {F_{{y'}}}{F_{{{\textit{z}}'}}} $$ {M_{{x'}}}{M_{{y'}}}{ {{M_{{{\textit{z}}'}}}} ]^{\rm T}}$,对应的铰链变形为$D = {\left[ {{\delta _{{x'}}}{\rm{ }}{\delta _{{y'}}}{\rm{ }}{\delta _{{{\textit{z}}'}}}{\theta _{{x'}}}{\theta _{{y'}}}{\theta _{{{\textit{z}}'}}}} \right]^{\rm T}}$,根据材料力学理论可得它们的关系为[1]

      $$D = CF$$ (13)

      其中,C为柔性铰链的柔度矩阵,公式为:

      $$C = \left[ {\begin{array}{*{20}{c}} {\dfrac{{{\delta _{{x'}}}}}{{{F_{{x'}}}}}}&0&0&0&0&0\\ 0&{\dfrac{{{\delta _{{y'}}}}}{{{F_{{y'}}}}}}&0&0&0&{\dfrac{{{\delta _{{y'}}}}}{{{M_{{z'}}}}}}\\ 0&0&{\dfrac{{{\delta _{{z'}}}}}{{{F_{{z'}}}}}}&0&{\dfrac{{{\delta _{{z'}}}}}{{{M_{{y'}}}}}}&0\\ 0&0&0&{\dfrac{{{\theta _{{x'}}}}}{{{M_{{x'}}}}}}&0&0\\ 0&0&{\dfrac{{{\theta _{{y'}}}}}{{{F_{{z'}}}}}}&0&{\dfrac{{{\theta _{{y'}}}}}{{{M_{{y'}}}}}}&0\\ 0&{\dfrac{{{\theta _{{z'}}}}}{{{F_{{y'}}}}}}&0&0&0&{\dfrac{{{\theta _{{z'}}}}}{{{M_{{z'}}}}}} \end{array}} \right]$$ (14)

      根据参考文献[6],柔性铰链工作方向上的弯曲应变能公式为:

      $${{U}} = \frac{1}{2}\left(\int_t {\frac{{F_{{x'}}^2}}{{EA({x'})}} {\rm{d}}{x'} + \int_t {\frac{{M_{{{\textit{z}}'}1}^2}}{{EI({x'})}} {\rm{d}}{x'}} } \right)$$ (15)
      $$\left\{ \begin{array}{l} {M_{{{\textit{z}}'}1}} = {M_{{{\textit{z}}'}}} + {F_{{y'}}}{x'} \\ A({x'}) = wt({x'}) \\ \end{array} \right.$$ (16)

      式中:E为柔性铰链材料的杨氏模量;$I({x'})$为柔性铰链截面对中心轴的惯性矩;${F_{{x'}}}$为铰链自由端上沿${x'}$轴向的拉力;${M_{{{\textit{z}}'}}}$为铰链自由端上沿${y'}$轴切向力引起的转矩以及绕${{\textit{z}}'}$轴的转矩之和。先对旋转力矩${M_{{{\textit{z}}'}1}}$引起的绕${{\textit{z}}'}$ 轴的转矩进行分析,根据卡式第二定理可得,柔性铰链的转角计算公式为[7]

      $${\theta _{{{\textit{z}}'}1}} = \frac{{\partial U}}{{\partial {M_{{{\textit{z}}'}1}}}}$$ (17)

      根据图2(b)可得:

      $$I({{{x}}'}) = \frac{{wt{{({{{x}}'})}^3}}}{{12}}$$ (18)
      $$\begin{array}{l} t({x'}) = 2a + t - 2\sqrt {{a^2} - \dfrac{{{a^2}}}{{{b^2}}}{x'}^2} = 2a + t - \cos \varphi \\ \end{array} $$ (19)
      $$ {\rm{d}}{x}^{\rm{'}}={\rm{d}}({{b}}\sin\varphi )={{b}} \cdot \cos\varphi {\rm{d}}\varphi $$ (20)

      将公式(19)代入公式(15)可得柔性铰链应变能为:

      $$\begin{array}{l} U = \dfrac{{{M_{{{\textit{z}}'}}}^2}}{{2E}}\displaystyle\int_t^{} {\dfrac{{12}}{{wt{{({x'})}^3}}}{\rm{d}}{x'}} = \dfrac{{6{M_{{{\textit{z}}'}}}^2}}{{wE}}\displaystyle\int_t^{} {\dfrac{1}{{t{{({x'})}^3}}}{\rm{d}}{x'}} \\ \end{array} $$ (21)

      再将公式(21)代入公式(17)可得转角${\theta _{{z'}1}}$为:

      $${\theta _{{{\textit{z}}'}1}} = \dfrac{{\partial U}}{{\partial {M_{{{\textit{z}}'}1}}}} = \dfrac{{12{M_{{{\textit{z}}'}}}}}{{wE}}\int_t^{} {\dfrac{1}{{t{{({x'})}^3}}}{\rm{d}}{x'}} $$ (22)

      将公式(19)、(20)和公式(22)代入公式(14)中,可计算出柔性铰链工作方向上的柔度为:

      $$ \begin{split}{C}_{{\theta }_{{{\textit{z}}}^{\rm{'}}},{M}_{{{\textit{z}}}^{\rm{'}}}}=&\dfrac{{\theta }_{{{\textit{z}}}^{\rm{'}}}}{{M}_{{{\textit{z}}}^{\rm{'}}}}=\frac{12}{wE}{\displaystyle {\int }_{t}^{}\frac{1}{t{({x}^{\rm{'}})}^{3}}{\rm{d}}{x}^{\rm{'}}} =\\ &\frac{12b}{wE}{\displaystyle {\int }_{-\frac{\pi }{2}}^{\frac{\pi }{2}}\dfrac{\mathrm{cos}\varphi }{{(2a+t-2a·{\rm{cos}}\varphi )}^{3}}{\rm{d}}\varphi } \end{split}$$ (23)

      $p = \dfrac{a}{t}$,则

      $$ {C}_{{\theta }_{{{\textit{z}}}^{\rm{'}}},{M}_{{{\textit{z}}}^{\rm{'}}}}=\frac{12b}{Ew{t}^{3}}{\displaystyle {\int }_{-\frac{\pi }{2}}^{\frac{\pi }{2}}\frac{\mathrm{cos}\varphi }{{(2p+1-2p \cdot \cos\varphi )}^{3}}{\rm{d}}\varphi }$$ (24)

      $ Y={\displaystyle {\int }_{-\frac{\pi }{2}}^{\frac{\pi }{2}}\frac{\mathrm{cos}\varphi }{{(2p+1-2p\cdot \cos\varphi )}^{3}}{\rm{d}}\varphi }$,因为Y求积分结果过于复杂,不利于工程上的优化设计应用,所以采用非线性曲线拟合的方法,对Y进行化简。曲线拟合的优度判断采用残差平方和(SSE)、判定系数(R-Square)、调整的判定系数(Adjusted R-Square)、均方根误差(RMSE)和残差分布图进行判定,残差平方和、均方根误差越接近0,判定系数、调整的判定系数越接近1,则曲线拟合程度越好[8]

      Y函数采用了幂函数拟合,如图6所示,由残差分布图以及各项判定系数可知,幂函数拟合结果较好,满足误差要求。

      图  6  (a)非线性拟合曲线图;(b)非线性拟合判定系数

      Figure 6.  (a) Non-linear fitting curve diagram; (b) Non-linear fitting determination coefficient

      拟合结果为:

      $$Y = 1.079{p^{ - 0.475\;3}}$$ (25)

      则深切口柔性铰链工作方向上的柔度为:

      $${C_{{\theta _{{z'}}},{M_{{z'}}}}} = \frac{{12b \times 1.079{p^{ - 0.475\;3}}}}{{Ew{t^3}}}$$ (26)

      根据参考文献[5],计算系统前四阶固有频率所需的深切口柔性铰链在其他次要方向上的柔度理论计算公式如下:

      $${C_{{\theta _{{{\textit{z}}'}}},{F_{{y'}}}}} = \frac{{5.027{b^2}}}{{Ew}}\left[\frac{{3.771}}{{{{(0.58a + t)}^3}}} + \frac{1}{{{t^3}}}\right]$$ (27)
      $${C_{{\theta _{{y'}}},{M_{{y'}}}}} = \frac{{5.027b}}{{E{w^3}}}\left[\frac{{3.771}}{{0.58a + t}} + \frac{1}{t}\right]$$ (28)
      $${C_{{\delta _{{x'}}},{F_{{x'}}}}} = \frac{{0.419b}}{{Ew}}\left[\frac{{3.771}}{{0.58a + t}} + \frac{1}{t}\right]$$ (29)
    • 根据模态分析结果可知,FSM系统前两阶固有频率振型分别为绕x,y两轴的扭振运动,故可将其简化为单自由度扭振模型,根据参考文献[19-10]可知,该模型在扭振方向上的运动刚度和该方向上固有频率的关系可表示为:

      $${f_{n1,2}} = \frac{1}{{2\pi }}\sqrt {\frac{{{K_{\theta n}}}}{{{J_n}}}} $$ (30)

      式中:${K_{\theta n}}$为运动方向的扭转刚度;${J_n}$为运动方向的转动惯量。

      快反镜系统第三阶振型绕z轴扭转,因为是小变形运动,故可将其简化为单自由度的扭转振动,将新推导出的三阶振型方向上的刚度计算公式代入计算,则三阶固有频率计算公式为:

      $${f_{n3}} = \frac{1}{{2\pi }}\sqrt {\frac{{{K_{{\theta _\text{z}}}}}}{{{J_\text{z}}}}} $$ (31)

      式中:${K_{{\theta _z}}}$为系统绕z轴的扭转刚度;${J_z}$为运动部分绕z轴扭转刚度。

      快反镜系统第四阶为沿z轴拉压,根据参考文献[8]可知,第四阶固有频率计算公式为:

      $${f_{n4}} = \frac{1}{{2\pi }}\sqrt {\frac{{{K_{{\Delta _\text{z}}}}}}{M}} $$ (32)

      式中:${K_{{\delta _z}}}$为系统沿z轴的拉压刚度;$M$为运动部分质量。

    • 利用ug软件对深切口柔性铰链进行建模,并导入Ansys软件进行了有限元仿真验证,如图7所示。

      仿真结果对比如表1所示。文中推导出的深切口理论计算公式计算结果与有限元仿真结果误差不超过9%,证明了理论推导结果的准确性。

      图  7  深切口柔性铰链工作刚度有限元验证

      Figure 7.  Finite element verification of the working stiffness of the deep-cut flexure hinge

      表 1  深切口柔性铰链工作刚度理论计算与有限元仿真对比

      Table 1.  Comparison of theoretical calculation and finite element simulation of working stiffness of deep-cut flexure hinge

      a
      /mm
      b
      /mm
      t
      /mm
      w
      /mm
      Theoretical
      calculation results
      /N·m·rad1
      Finite element
      simulation results
      /N·m·rad1
      Error
      percentage
      1081825.379926.0417−2.54%
      1061833.839932.65133.64 %
      530.588.46008.25272.45%
      530.555.28755.14662.74%
      520.557.93127.29228.76%
    • 利用ug软件对快反镜系统进行建模,并导入Ansys软件进行了模态分析,如图8所示。

      图  8  第三阶振型仿真图

      Figure 8.  Third-order mode shape simulation diagram

      有限元仿真验证结果如表2所示。根据有限元仿真结果可知,该理论公式与有限元仿真结果误差不超过1.7%,证明了三阶振型计算公式以及新推导系统绕z轴扭转刚度计算公式的准确性。

      表 2  快反镜系统三阶固有频率理论计算与有限元仿真对比

      Table 2.  Comparison of theoretical calculation and finite element simulation of three-boundary natural frequency of fast mirror system

      a
      /mm
      b
      /mm
      w
      /mm
      t
      /mm
      Mirror
      thickness: h
      /mm
      Theoretical
      calculation
      value of
      third-order
      mode/Hz
      Three-order
      mode finite
      element simulation
      value/Hz
      Error
      percentage
      43.550.823284.40286.250.61%
      1088115218.90220.17−0.58%
      658115434.29427.111.68%
    • 文中要求设计的快反镜系统镜面直径D为100 mm,反射镜厚度h为15 mm,隔离板厚度为3 mm。反射镜系统偏转范围为$ \pm {\rm{3\;mrad}}$,重复定位精度${\delta _p} \leqslant 8\;{\rm{urad}}$。柔性支持系统材料为TC4,隔离板材料为铟钢,反射镜材料为微晶玻璃。初始结构参数如表3所示。

      表 3  快反镜系统初始结构参数

      Table 3.  Initial structure parameters of the fast steering mirror system

      a/mmb/mmt/mmw/mmL/mmD/mmh/mm
      108182510015

      图3中,参数abt为自变量,参与函数优化设计;wLDh为定值,不参与优化设计。

    • 在快速控制反射镜的系统设计中,为提高系统控制带宽,要求系统结构的固有频率高于系统的控制带宽的2~4倍。但是这样的设计会导致固有频率过大,在FSM系统转动惯量不变的情况下,系统工作方向上的刚度就会过大,这会极大地增大系统驱动器的负担[1]。故所采用的设计方法是尽可能地减小系统工作方向上的低阶固有频率,并通过控制系统对其加以抑制,而另一方面提高系统非工作方向上的高阶固有频率,并且越大越好,以提高系统的控制带宽[2]

      根据公式(30)~(32)可知,系统的四阶固有频率大小${f_n}$与柔性支撑系统三轴转动刚度的平方根$\sqrt {{K_{\theta n}}} $z轴拉伸刚度的平方根$\sqrt {{K_{{\Delta _n}}}} $成正比关系,因为其比值随系统负载的质量和转动惯量不同而变化。为适应不同负载的快反镜系统,故将转动惯量和质量值简化为各阶固有频率前的加权系数,构建如下多目标优化函数:

      $$\min G(x) = \frac{{{\beta _1}\sqrt {{K_{\theta_x}}} + {\beta _2}\sqrt {{K_{\theta_y}}} }}{{{\beta _3}\sqrt {{K_{\theta_\text{z}}}} + {\beta _4}\sqrt {{K_{\Delta_\text{z}}}} }}$$ (33)

      式中:${\;\beta _1} \sim {\;\beta _4}$为柔性支撑系统各运动方向刚度平方根的加权系数,反映了各运动方向刚度平方根在固有频率优化中的重要性。

      在多目标优化问题中,采用不同的加权值会产生不同的优化结果,故一般会产生最优解集,需根据各子目标的重要程度不同以及多次的参数调整来最终确认加权值,从而最终计算出一组最优解[4]

      文中采用固定权重值法来解决该多目标优化问题。在固有频率的优化计算中,振型为绕y轴扭转的二阶固有频率决定了低阶工作方向固有频率的上限,振型为沿z轴扭转的三阶固有频率决定了高阶非工作方向固有频率的下限。根据FSM系统设计提高非工作方向上的固有频率,同时减小工作方向上的固有频率这一设计思想,着重在降低二阶固有频率的同时提高三阶固有频率是该优化计算中比较重要的部分,故$\sqrt {{K_{\theta_y}}} $$\sqrt {{K_{\theta_\text{z}}}} $两个子目标在优化函数中占据比较重要的地位,在分配权重时需要予以更大的权重值。根据以上需求,并通过调参计算,选取的加权系数为:β1= 0.2,β2= 0.3,β3=0.4,β4=0.1。

    • (1)柔性铰链强度不等式约束

      柔性铰链弯曲时,铰链的最大应力出现在铰链的最小切割厚度t处,同时有应力集中的影响,故根据材料力学中的纯弯曲理论,柔性铰链弯曲时最大应力为[8]

      $${\sigma _{\max }} = {K_t}\frac{{6{M_{\max }}}}{{{t^2}w}}$$ (34)

      式中:${K_t}$为应力集中系数[8]${M_{\max }}$为铰链在最大偏转角时外部施加的转矩。其计算公式为:

      $$\left\{ \begin{array}{l} {K_t} = {\left(1 + \dfrac{{at}}{{2{b^2}}}\right)^{\frac{9}{{20}}}} \\ {M_{\max }} = \dfrac{{Ew{t^3} \cdot {{10}^{ - 3}}}}{{4.316\;b{p^{ - 0.475\;3}}}} \\ \end{array} \right.$$ (35)

      将公式(35)代入公式(34)可得:

      $${\sigma _{\max }} = {\left(1 + \frac{{at}}{{2{b^2}}}\right)^{\frac{9}{{20}}}} \times \frac{{6Et \cdot {{10}^{ - 3}}}}{{4.316b{{\left(\dfrac{a}{t}\right)}^{ - 0.475\;3}}}}$$ (36)

      TC4屈服强度为$4.3 \times {10^8}$ Pa,安全系数为1.5,将屈服强度和安全系数代入公式(36)可得:

      $${c_{\rm{1}}}(x) = {(1 + \frac{{at}}{{2{b^2}}})^{\frac{9}{{20}}}} \times \frac{{6Et \cdot {{10}^{ - 3}}}}{{4.316b{{\left(\dfrac{a}{t}\right)}^{ - 0.475\;3}}}} - \frac{{4.3 \times {{10}^2}}}{{1.5}} \leqslant 0$$ (37)

      (2)柔性铰链运动精度不等式约束

      在对柔性铰链运动精度的分析中,一般将柔性铰链几何中心点的位移,即柔性铰链偏转中心的位移作为衡量柔性铰链运动精度的定量指标[8],根据材料力学的公式可得:

      $$yc = \iint {\frac{{M({x'})}}{{EI({x'})}} {\rm{d}}{x'}{\rm{d}}{x'}}$$ (38)

      其中,$\left\{ \begin{array}{l} I({{{x}}'}) = \dfrac{{{{w}}{{{t}}^{\rm{3}}}{\rm{(}}{{{x}}'}{\rm{)}}}}{{{\rm{12}}}} \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;{\text{。}}\\ M({{{x}}'}) = M \\ \end{array} \right.$

      将公式(38)化为极坐标形式可得:

      $$ \begin{split} &\dfrac{{{{12Mb}}}}{{{{Ew}}}} \displaystyle\iint {\frac{{{\rm{cos\varphi }}}}{{{{{{(2a + t - 2a \cos\varphi )}}}^{\rm{3}}}}}{{{\rm{d}}\varphi b\cos\varphi {\rm{d}}\varphi }}} =\\ & \dfrac{{{{12Mb}}}}{{Ew{t^3}}} \displaystyle\int_{ - \frac{\pi }{2}}^0 {\int_{ - \frac{\pi }{2}}^\varphi {\dfrac{{\cos \beta }}{{{{(2p + 1 - 2p\cos \beta )}^3}}}} } {\rm{d}}\beta b\cos \varphi {\rm{d}}\varphi = \\ & \dfrac{{6{{Mb}}{^2}}}{{Ew{t^3}}} \cdot \dfrac{1}{{1 + 2p}} \\ \end{split} $$ (39)

      将公式(19)代入公式(39),并代入柔性铰链重复定位精度${\delta _p} \leqslant 8 $ μrad,可得出柔性铰链运动精度约束表达式为:

      $${c_2}(x) = \frac{{6 \times {{10}^{ - 3}}b}}{{4.316 \cdot \left(1 + \dfrac{{2a}}{t}\right) \cdot {{\left(\dfrac{a}{t}\right)}^{ - 0.475\;3}}}} - 8 \times {10^{ - 3}} \leqslant 0$$ (40)
    • 优化模型为:

      $$\begin{array}{l} {f_{obj}} = \min G(x) = \dfrac{{{\beta _1}\sqrt {{K_{\theta_x}}} + {\beta _2}\sqrt {{K_{\theta_y}}} }}{{{\beta _3}\sqrt {{K_{\theta_{\textit{z}}}}} + {\beta _4}\sqrt {{K_{\Delta_{\textit{z}}}}} }} \\ {{suppose:x1 = a,x2 = b,x3 = t,x = (x1,x2,x3}}{{\rm{)}}^{\rm{T}}} \\ s.t.\left\{ \begin{array}{l} c(x) \leqslant 0 \\ A \cdot x \leqslant b \\ lb < x < ub \\ \end{array} \right. \\ \end{array} $$ (41)

      将柔性铰链的长轴a、短轴b和最小切割厚度t作为变量。同时根据几何约束条件可得出线性不等式约束矩阵A为:[0 −1 5 0; −1 0 5 0; 1 0 −50 0; −1 1.2 0 0],b=[0; 0]。根据现有的加工精度,得出边界约束lb=[0,0,0.8],ub=[20; 20; 5]。非线性不等式约束为柔性铰链强度不等式约束c1,柔性铰链运动精度不等式约束c2。

    • 遗传算法(Genetic Algorithm,GA)最早是由美国的 John holland于20世纪70年代提出的,该算法的设计灵感来自于大自然中生物体进化的原理,将达尔文生物进化论的自然选择和遗传学机理的生物进化过程的计算模型作为搜索最优解的方法[11]。该算法通过数学的方式,将问题的求解过程转换成类似生物进化中的染色体基因的交叉、变异等过程[11]。与一般的优化算法相比,遗传算法从问题解的串集开始搜索,而不是从单个解开始。传统优化算法是从单个初始值迭代求最优解的,容易误入局部最优解。而遗传算法从串集开始搜索,覆盖面大,利于全局择优,不易陷入局部最优解[12]。在求解较为复杂的组合优化问题时,相对一些常规的优化算法,通常能够较快地获得较好的优化结果,遗传算法的基本运算过程如图9所示[13]

      图  9  遗传算法优化流程图

      Figure 9.  Optimization flow chart of genetic algorithm

    • 文中采用matlab遗传算法工具箱进行优化计算,优化初始值为a=10,b=8,t=1。设定的优化迭代次数上限为100代,终止计算条件为两次迭代的适应度函数度值的平均变化小于1×10−6并且不违反约束。优化计算结果为a=4.798,b=3.999,t=0.8。为便于加工,对以上优化结果进行近似处理,所得的最终优化结果参数如表4所示。

      表 4  遗传算法优化结果

      Table 4.  Optimization results of genetic algorithm

      ParameterValue
      a/mm4.8
      b/mm4
      t/mm0.8
      G(x)0.0181

      利用以上优化后的结构参数,并结合之前所展示的系统结构其他关键部件参数,对该结构进行了ug 建模,优化结构三视图如图10所示。

      图  10  (a)前视图;(b)右视图;(c)仰视图

      Figure 10.  (a) Front view; (b) Right view; (c) Bottom view

      利用之前推导的刚度计算公式计算该优化结构四个振型方向上的刚度,并与初始结构的四个振型方向刚度做对比,以初步分析优化计算的效果,对比结果如表5所示。

      表 5  多目标优化结构与初始结构系统刚度对比

      Table 5.  Comparison of system stiffness between multi-objective optimized structure and initial structure

      a/mmb/mmt/mmK1/N·m·rad1K2/N·m·rad1K3/N·m·rad1K4/N·m1G(x)
      Initial value 10 8 1 50.04 50.04 844.6296 1.323e+06 0.0279
      Optimized value 4.8 4 0.8 40.51 40.51 3.3602e+03 2.343e+06 0.0181
      Optimization rate −19.04% −19.04% 297.83% 77.09% 41.58%

      从以上计算结果可以看出,优化后的快反镜结构工作方向刚度减19.04%,非工作方向刚度提高了297.83%和77.09%,优化函数值下降41.58%。因为系统固有频率与系统刚度成正相关,与系统负载成负相关,故在反射镜负载不变的情况下,前二阶固有频率会下降,而三、四阶固有频率会上升,低级固有频率与高阶固有频率的差距会进一步提高,故该优化结构初步符合设计预期。

      为进一步分析优化结构的固有频率,并得到更加直观的结果,文中将优化后的结构进行Ansys模态仿真,结果如图11所示,仿真结果以及与初始结构对比如表6所示,从以上仿真结果可以看出,经过优化计算,快反镜一、二阶固有频率较初始结构分别下降8.08%和5.40%,分别为55.8 Hz和57.5 Hz。三、四阶固有频率较初始结构上升112.59%和16.80%,达到503.7 Hz和642.4 Hz,可见优化计算达到其设计目的,满足快反镜系统的设计要求。其中最为关键的二、三阶固有频率得到了较大的优化,证明了加权系数的有效性。该优化计算有效降低了结构工作方向上的低阶固有频率,为下一步加入控制系统打下了良好的基础,使系统的低阶固有频率可以得到有效抑制。同时,该优化计算大大提高了结构非工作方向上的三、四阶固有频率,使系统系统控制带宽可以得到显著提高,预计在加入控制系统抑制低阶频率后系统控制带宽可达300 Hz左右。

      图  11  (a)一阶模态;(b)二阶模态;(c)三阶模态;(d)四阶模态

      Figure 11.  (a) First-order mode; (b) Second-order mode; (c) Third-order mode; (d) Fourth-order mode

      表 6  优化结构与初始结构四阶固有频率对比

      Table 6.  Comparison of the fourth-order natural frequency between the optimized structure and the initial structure

      First-order mode/HzSecond-order mode/HzThird-order mode/HzFourth-order mode/Hz
      Initial structure60.71560.814236.93549.97
      Optimized structure55.80757.533503.68642.38
      Optimization rate−8.08%−5.40%112.59%16.80%
    • 为对二自由度深切口柔性铰链快速控制反射镜系统进行控制带宽优化设计,文中首先对快速控制反射镜系统四阶固有频率进行了有限元仿真分析,并对四个振型方向上工作刚度进行了理论分析推导,解决了传统计算公式不适用于三阶振型方向的问题;然后利用能量法和卡式第二定理,对深切口柔性铰链的工作刚度计算公式进行推导,并利用非线性拟合方法对计算公式进行了拟合简化,得出的简化计算公式与有限元仿真结果误差不超过9%,解决了传统公式过于复杂的问题;接着利用以上推导出的刚度公式进行了快反镜系统第三阶固有频率理论计算和有限元仿真验证,计算结果与仿真结果误差不超过1.7%,证明了新推导的系统三阶振型方向刚度和三阶固有频率计算公式的可靠性;最后以系统一、二阶固有频率的平方根最小,三、四阶固有频率的平方根最大为目标进行了多目标优化设计,得出的最优化结构有效降低了一、二阶固有频率,提高了三、四阶固有频率,并进行了有限元验证,证明了优化结果的可靠性和有效性。文中所推导的公式和提出的优化方法对二自由度柔性铰链快反镜系统的固有频率理论分析和提高控制带宽有一定的工程意义。

参考文献 (13)

目录

    /

    返回文章
    返回