Processing math: 100%
FU Qiang, YAN Lei, TAN Shuang-long, LIU Yang, WANG Ling-jie, ZHANG Xin. Light-and-small optical systems by metal-based additive manufacturing[J]. Chinese Optics, 2022, 15(5): 1019-1028. doi: 10.37188/CO.2022-0128
Citation: ZHAO Chen-hang, LU Qi-peng, SONG Yuan, GONG Xue-peng, WANG Yi, XU Bin-hao. Design and analysis of stress-free clamping of mirrors used in free-electron laser beamlines[J]. Chinese Optics, 2020, 13(4): 787-794. doi: 10.37188/CO.2019-0131

Design and analysis of stress-free clamping of mirrors used in free-electron laser beamlines

Funds:  Supported by the National Natural Science Foundation of China (No. 11079035); National Science and Technology Major Project (No. 2012ZX02702001); National Key R&D Program of China (No. 2018YFD0401003-04)
More Information
  • Corresponding author: luqp@ciomp.ac.cn
  • Received Date: 27 Jun 2019
  • Rev Recd Date: 07 Aug 2019
  • Publish Date: 01 Aug 2020
  • The reflector is an important optical element in free-electron laser beamlines. Deformation error caused by gravity can seriously affect the image quality of a beamline. To reduce deformation error, a gravity compensation scheme based on the Bessel point theory is proposed and a stress-free clamping device is designed. Taking a 440 mm × 50 mm × 50 mm mirror as an example, the analysis results indicate that the deformation error in the bottom surface of a mirror clamped with the traditional support method is 1.647 μrad. Adopting the newly designed device proposed in this paper, the results of a finite element analysis showed that the deformation error reduced to 0.085 7 μrad, which is better than the engineering index of 0.1 μrad. To prevent the mirror from moving when switching modes, a small clamping force of no more than 2 N can be added to the mirror, at which point the surface error of the mirror becomes 0.093 9 μrad. Additionally, a dynamic analysis of the device is also carried out, which indicates that the device mutes the low natural frequency, which means that resonance will not occur during operation. Therefore, this scheme satisfies our requirements for the beamline.

     

  • 深空低冷背景环境下目标探测是一个相对广阔的领域,如深空探测、天基预警、空间攻防等任务都属于该领域范畴,是各大国必争之地[1]。折反射式光学系统被广泛用于低冷目标探测领域[2]。随着科学技术的不断发展,折反射光学系统向着灵巧型、轻质、低辐射、大视场和高分辨率等方向发展。

    对于折反射式光学系统前端的反射光学元件,随着单点金刚石加工技术的发展,金属反射镜因其具有加工工艺性好、材料价格较低等优势逐渐获得认可和应用,典型的金属基材料包括铝合金、镁铝合金、铍铝合金等[3-5]。随着轻质高刚度指标要求的进一步提高,传统机械加工方式加工的铝合金反射镜越来越不能满足需求,这在一定程度上限制了铝合金反射镜的应用。为解决上述问题,基于增材制造的铝合金反射镜技术应运而生,成为近年来的研究热点[6-7]。增材制造金属基反射镜的加工工艺与传统金属反射镜相同,同时可以将拓扑优化思想应用于金属反射镜的设计,设计封闭式金属反射镜,从而获得传统金属反射镜难以实现的轻质和高刚度。

    2015年,美国康宁公司通过增材制造技术制备了一个蜂窝状轻质高性能铝镜,与传统金属反射镜基体制备相比,提高了加工成形效率[8-9]

    德国弗琅禾费研究所提出一种夹层式封闭蜂窝状结构,通过内部加强筋上的孔使所有蜂窝状空间连通,反射镜刚度高于传统金属反射镜[10]

    2017年,Hilpert E等人通过对比5种结构形式的金属反射镜,分析了增材制造技术制备的金属反射镜的优势[11]。2019年,Hilpert E等人对反射镜的轻量化形式进行了进一步的优化设计,使其在保证金属反射镜刚度的情况下,达到了60.5%的轻量化效果[12]

    本文针对深空低冷目标的探测需求,提出一种光学系统局部制冷方法,将传统常温折反射式光学系统中的透镜组部分,放置到探测器杜瓦内部,与探测器芯片一起封装到杜瓦中制冷,实现光学自辐射抑制。在此基础上,设计完成了局部制冷型折反射式光学系统,光学系统口径为55 mm,焦距为110 mm,视场达到4°×4°。利用拓扑优化方法对前组反射式系统进行了优化设计,利用增材制造方法打印完成主镜组件、次镜组件和连接筒。利用单点金刚石车削的方法进行光学加工,针对打印表面存在的缺陷,采用镍磷改性工艺加以解决,最终完成的主镜和次镜面形精度高。最后,对光机装调后的系统性能进行了测试。

    拓扑优化方法本质上是对于特定的设计区间,寻求最佳的分布,以实现最优设计构型。目前常用的拓扑表达形式包括:变厚度法、均匀化方法和相对密度法等[13]

    采取相对密度法进行拓扑优化设计,相对密度法是在均匀化方法的基础上改进的一种方法,以单元相对密度作为设计变量,相对密度设置为0~1,相对密度越接近1,说明该单元在设计中越发重要,表明该单元应该保留;而当相对密度接近零时,表明该单元材料可以去除。目前相对密度法中最为常用的是各向正交惩罚材料密度(SIMP)法,该方法在许多通用有限元分析软件中应用广泛,其单元弹性模量和相对密度关系式为:

    Ei=(α0+(α1α0)αiP)E0i=1,2,,N, (1)

    其中,$ {E_i} $为弹性模量,$ {\alpha _i} $为单元相对密度,$ {\alpha _0} $为单元相对密度下限,$ {\alpha _1} $为单元相对密度上限,$ {E_0} $为材料弹性模量,$ N $为设计区间的单元数量,$ P $为惩罚因子。

    通过设定拓扑优化单元阈值$ {\alpha _c} $,并对比$ {\alpha _i} $$ {\alpha _c} $大小对设计区间内单元进行取舍,可以将设计变量表述为:

    α=(α1,α2,,αN)T. (2)

    基于此,对反射镜支撑结构展开拓扑优化设计,通过分析传力路径,寻求设计空间内最优材料分布,以实现支撑区域高刚度、高轻量化设计。以体积为约束条件,最小应变能为设计目标,建立数学优化模型,具体如式(3)所示:

    min:J(α)=UTKUis.t.{KU=FV(α)=fV00αi1, (3)

    其中,$ J $为反射镜支撑区域应变能,$\boldsymbol\alpha$为设计变量向量,$ {\alpha _i} $为第$ i $个设计变量,$\boldsymbol{K}$为全局刚度矩阵,$\boldsymbol U$为全局位移向量,$\boldsymbol F$为全局载荷向量,$V$为设计域体积约束,$ f $为体积分数,$ {V_0} $设计域体积。

    光学系统具体设计指标要求见表1。光学系统的谱段为长波红外,选择的长波红外探测器指标见表2。由于所选探测器为制冷型探测器,故有冷光阑效率的要求,需要达到100%。主次镜组件重量要求较苛刻,要求不大于100 g。

    表  1  光学系统设计指标
    Table  1.  Design index of the optical system
    参数指标要求
    波段/μm8~10
    相对孔径1∶2
    焦距/mm110
    视场/(°)4×4
    冷光阑效率100%
    主次镜组件重量不大于100 g
    下载: 导出CSV 
    | 显示表格
    表  2  长波红外探测器指标
    Table  2.  Index of the long-wave infrared detector
    参数指标数据
    阵列尺寸256×256
    像元间距30 μm×30 μm
    F数2
    下载: 导出CSV 
    | 显示表格

    为实现光学系统的轻小型化,在光学系统结构形式上,采用同轴折反射式结构形式。指标中光学系统的视场较大,达到4°×4°。考虑将由透镜组成的中继镜组放置在红外探测器的杜瓦中,光学系统元件数尽量少,数量少于等于3片为宜。考虑中继镜组的第一片透镜同时作为探测器杜瓦的窗口片,中继透镜组的口径和间隔尺寸应尽量紧凑,上述的一系列要求给光学系统设计带来了一定挑战。

    设计后的光学系统二维图如图1所示。光学系统由主镜、次镜以及中继镜组组成,其中主镜和次镜为非球面设计,均为高次非球面;中继镜组由三片透镜组成,正光焦度材料为锗,负光焦度材料为硒化锌,冷阑放置在透镜2的后表面上。系统为了冷阑匹配采用二次成像结构形式,中间像面位于主镜和次镜之间。为尽量减小镜体自身辐射对探测性能的影响,光阑匹配是设计的重点,冷光阑效率达到100%。中心遮拦影响能量集中度和光学传递函数,设计后系统面遮拦比为18.3%。

    图  1  光学系统二维图
    Figure  1.  2-D diagram of the optical system

    图2给出了光学调制传递函数曲线,由图2可以看出,全视场范围内各视场点像质均达到衍射极限。

    图  2  光学系统调制传递函数曲线
    Figure  2.  Modulation transfer function curves of the optical system

    对于点目标探测系统来说,能量集中度是一个重要的衡量指标,图3给出了光学系统的能量集中度曲线。可见,在一个像元范围内,能量集中度大于50%,受中心遮拦的影响,能量集中度低于传统透射式系统。

    图  3  光学系统能量集中度曲线
    Figure  3.  Enclosed energy curves of the optical system

    轻小型金属基光学系统结构包含前组和后组两部分。其中,前组包括:主镜组件、遮光筒、次镜组件;后组包括:透镜组和透镜壳等。为了实现结构超轻设计,同时考虑到光学系统与总体的接口关系,在设计上,采用拓扑优化方法对前组进行设计,工艺上采用3D打印方式实现。设计完成的系统结构如图4所示。

    图  4  轻小型金属基光机系统图
    Figure  4.  Diagram of light-and-small metal-based optical system
    3.3.1   主镜组件设计

    为了使零部件最少,同时尽量避免安装过程中引入的外力等,将主镜及其背板进行一体式设计,在镜体和支撑之间添加柔性设计,有效起到隔绝外力及温度变化产生的应力的作用。运行Hypermesh软件中的Optistruct优化模块,对初始结构进行拓扑优化设计,为保证主镜的面形精度,结合点阵结构对镜体内部进行适度填充,在满足轻量化设计思想下,有效提高镜面刚度。拓扑优化完成的主镜组件如图5所示,在自重条件下的变形情况如图6(彩图见期刊电子版)所示,对主镜的面形进行拟合分析,主镜面形RMS为3 nm(约λ/210,λ=632.8 nm),对成像质量影响很小,可以忽略。

    图  5  主镜组件的(a)正视图和(b)侧视图
    Figure  5.  (a) Front view and (b) side view of primary mirror assembly
    图  6  主镜在自重下的变形情况
    Figure  6.  Deformation of the primary mirror under its own gravity
    3.3.2   次镜组件设计

    与主镜组件形式类似,次镜及其支撑进行一体化设计,对次镜组件支撑筋形式进行多轮优化,组件内部多处采用中空设计,采用放射形三根筋支撑方案,三根筋呈环切式分布,可有效减小组件弯曲应力,提高结构承受强冲击和振动能力。经优化后,每根筋径向宽度为1.8 mm,在满足光学通光要求下,具备足够刚度。

    同时,次镜背部结构采用共面设计,可为光加和装调阶段提供高定位基准。次镜组件结构形式如图7所示。

    图  7  次镜组件
    Figure  7.  Secondary mirror assembly
    3.3.3   遮光筒设计

    遮光筒作为连接主镜组件和次镜组件的转接结构,需承载次镜组件的重量,同时起到遮挡杂光的作用。

    在设计上,遮光筒与主次组件类似,经拓扑优化,确定了遮光筒外形尺寸,结合光学系统要求,同时考虑加工工艺性,对遮光筒筒壁进行封闭处理,优化后遮光筒最薄壁厚0.8 mm,在拓扑优化结果的基础上,进一步进行轻量化设计,采用中空结构结合局部支撑的方式,进一步实现超轻量化设计。遮光筒结构如图8所示。

    图  8  遮光筒结构
    Figure  8.  Shading baffle structure
    3.3.4   系统模态及静力分析

    对光学系统进行整机模态分析,以验证整机及各组件刚度分布。整机模态分析结果如图9(彩图见期刊电子版)所示,系统三阶和四阶模态达到1213.7 Hz,整机刚度高。

    图  9  光学系统整体模态分析结果
    Figure  9.  Results of mode analysis of the whole optial system

    为了考察系统在自身重力下,因重力导致的系统变形情况,开展系统自重下静力分析,模拟实际安装状态,对主镜背板3个安装孔进行6个自由度全约束处理,模拟系统在1 g重力下的变形情况,分析结果如图10(彩图见期刊电子版)所示。由分析可知,系统在自身重力下,最大合成位移为5.86×10−4 mm,属于亚微米量级,结合模态分析结果,结构具备足够刚度。

    图  10  系统自重下的静力分析结果
    Figure  10.  Results of system static analysis under its own grarity
    3.3.5   重量估算

    对主镜组件、次镜组件、遮光筒和螺钉进行估算,总质量为91.1 g,小于100 g,满足指标要求,详见表3

    表  3  主、次镜组件、遮光筒和螺钉的重量估算
    Table  3.  Weight estimation of primary and secondary mirror assemblies, shadding baffle and screws
    名称重量g
    主镜组件40.6
    次镜组件16.8
    遮光筒32.2
    螺钉1.50
    合计91.1
    下载: 导出CSV 
    | 显示表格

    主镜组件、次镜组件和连接筒均由增材制造3D打印完成,反射镜的材料为AlSi10Mg。打印完成后对组件进行了高低温时效处理,图11~图13分别为打印完成的主镜组件、次镜组件和连接筒。

    图  11  打印完成的主镜组件
    Figure  11.  Primary mirror assembly by additive manufacturing
    图  12  打印完成的次镜组件
    Figure  12.  Secondary mirror assembly by additive manufacturing
    图  13  打印完成的连接筒
    Figure  13.  Shading baffle by additive manufacturing

    用单点金刚石车床对主镜和次镜进行光学加工,光学加工后的主镜和次镜分别如图14(a)和图15(a)所示,主镜面形RMS值达到0.044 μm,见图14(b)(彩图见期刊电子版),次镜面形RMS值达到0.028 μm,见图15(b)(彩图见期刊电子版)。虽然面形满足使用要求,但从图中可以看出,镜子表面有明显的环带和点状斑点。

    图  14  单点车削后的主镜及其面形数据
    Figure  14.  Primary mirror after SPDT and it′s surface shape data
    图  15  单点车削后的次镜及其面形数据
    Figure  15.  Secondary mirror after SPDT and it′s surface shape data

    由于光学加工后的主镜和次镜表面有明显的环带和点状斑点,与通常非3D打印的铝合金表面有明显差异,故用高倍相机对表面进行了测试,测试图见图16。从图16可以看出,环带非常明显,且点状斑点均为小坑,对图中最大的圆形坑尺寸进行测量,直径达到了0.3 mm。预计上述表面对光学系统的能量集中度影响会很大,不能直接使用。为解决该问题,对光学元件表面进行镀镍磷合金改性。

    图  16  高倍相机下的主镜表面
    Figure  16.  Primary mirror surface observed by a high magnification camera

    对光学加工后的主镜和次镜光学表面进行改性处理,在表面镀制一层镍磷。为防止镍磷进到3D打印金属镜内部,对金属镜的排粉孔表面进行了保护处理。为了让改性层致密,对镜面进行喷砂处理。图17给出了镀镍磷改性完成的主镜。

    图  17  镍磷改性完成的主镜
    Figure  17.  Primary mirror modified by Ni-P coating

    镀完镍磷后,对主镜和次镜重新进行光学加工,加工后的主镜和次镜分别如图18(a)和图19(a)所示。从图中可以看出,镜子表面质量得到明显改善。主镜面形RMS值为0.044 μm,见图18(b)(彩图见期刊电子版)。次镜面形RMS值达到了0.018 μm,见图19(b)(彩图见期刊电子版)。

    图  18  改性后光学加工完成的主镜及面形数据
    Figure  18.  Primary mirror and surface shape data after optical processing
    图  19  改性后光学加工完成的次镜及其面形数据
    Figure  19.  Second mirror and surface shape data after optical processing

    光学加工完成后对主镜和次镜进行镀膜处理,镀膜为金膜。镀膜后的主镜和次镜见图20。接着,对镀膜后的主镜和次镜进行面形测试,主镜的面形测试结果见图21(彩图见期刊电子版),主镜面形RMS值为0.076λ,为48 nm。与镀膜前相比,RMS值改变量为4 nm,考虑到测量误差,这个变化可忽略。次镜的面形测试结果见图22(彩图见期刊电子版),次镜面形RMS值为0.03λ,为19 nm,与镀膜前18 nm相比,RMS值改变量为1 nm,考虑到测量误差,面形镀膜前后可视为无变化。

    图  20  镀金后的主镜和次镜
    Figure  20.  Primary and secondary mirrors after gold plating
    图  21  镀金后的主镜面形数据
    Figure  21.  Surface quality of the primary mirror after gold plating
    图  22  镀金后的次镜面形数据
    Figure  22.  Surface quality of the second mirror after gold plating

    利用定心仪对主镜组件和次镜组件进行定心装调,见图23

    图  23  主次镜组件定心装调实验图
    Figure  23.  Centering alignment experiment of the primary and secondary mirror assembly

    整机装配完成后,利用传递函数测试仪对核心指标进行测试,测试现场如图24所示。光学系统焦距测试结果为108.26 mm,与设计值110 mm偏差为1.74 mm,在±3%的公差范围内,见图25

    图  24  整机性能测试实验图片
    Figure  24.  The performance test of the optical system
    图  25  焦距测量结果
    Figure  25.  The measurement result of the focal length

    对光学系统不同视场的MTF进行了测试,图26分别给出了中心视场、0.7视场和1视场下的MTF曲线,在Nyquist频率16.7 lp/mm处的测试值均大于0.35,接近衍射极限。

    图  26  不同视场调制传递函数测试结果
    Figure  26.  MTF test results for different FOVs

    对主次镜组件进行称重,见图27,重量为96.04 g。各组件详细重量数据见表4,比表3预估的重量(91.1 g)增加了4.94 g。其中,主镜组件增加了4.04 g,重量增加主要是由3D打印偏差、镍磷改性和镀金导致的。遮光筒重量增加了0.87 g,3D打印的实际尺寸与理论设计略有差别,而次镜组件和螺钉质量基本一致。满足小于100 g的技术要求。

    图  27  重量测试
    Figure  27.  Weight test
    表  4  各部件的实测重量
    Table  4.  Weight test results of each assembly
    名称实测质量(g)
    主镜组件44.64
    次镜组件16.8
    遮光筒33.07
    螺钉1.53
    合计96.04
    下载: 导出CSV 
    | 显示表格

    本文针对深空低冷目标探测需求,设计完成局部制冷型折反射式光学系统,将透镜组放置在红外探测器的杜瓦中,光学系统口径为55 mm,焦距为110 mm,视场达到4°×4°。利用拓扑优化方法对主次镜组件进行了优化设计,利用增材制造方法打印完成主镜组件、次镜组件和连接筒。利用单点金刚车削方法进行光学加工,针对打印表面存在缺陷的问题,采用镍磷改性工艺,再次单点加工完成镜面加工,测试结果表明其面形精度较高。定心装调后,对整机性能进行了测试。测试结果如下:光学系统焦距为108.26 mm,满足指标要求;全视场范围内调制传递函数均达到衍射极限;重量仅为96.04 g,小于100 g,满足要求。通过整机研制,充分说明金属基增材制造方法可以作为提升光学系统性能的有效手段。

  • [1]
    赵振堂, 冯超. X射线自由电子激光[J]. 物理,2018,47(8):481-490. doi: 10.7693/wl20180801

    ZHAO ZH T, FENG CH. X-ray free electron lasers[J]. Physics, 2018, 47(8): 481-490. (in Chinese) doi: 10.7693/wl20180801
    [2]
    ASSOUFID L, HIGNETTE O, HOWELLS M, et al. Future metrology needs for synchrotron radiation grazing-incidence optics[J]. Nuclear Instruments and Methods in Physics Research Section A:Accelerators,Spectrometers,Detectors and Associated Equipment, 2001, 467-468: 267-270.
    [3]
    卢启鹏, 高飒飒, 彭忠琦. 同步辐射水平偏转压弯镜面形误差分析与补偿[J]. 光学精密工程,2011,19(11):2644-2650. doi: 10.3788/OPE.20111911.2644

    LU Q P, GAO S S, PENG ZH Q. Analysis and compensation of slope error for synchrotron radiation horizontal deflected mirror[J]. Optics and Precision Engineering, 2011, 19(11): 2644-2650. (in Chinese) doi: 10.3788/OPE.20111911.2644
    [4]
    孙福权, 傅远, 祝万钱, 等. 压弯镜系统自重平衡多点调节方法的研究[J]. 核技术,2011,34(4):246-250.

    SUN F Q, FU Y, ZHU W Q, et al. A study on multi-point gravity compensation of mirror bending system[J]. Nuclear Techniques, 2011, 34(4): 246-250. (in Chinese)
    [5]
    MORI Y, YAMAMURA K, SANO Y. The study of fabrication of the x-ray mirror by numerically controlled plasma chemical vaporization machining: development of the machine for the x-ray mirror fabrication[J]. Review of Scientific Instruments, 2000, 71(12): 4620. doi: 10.1063/1.1322580
    [6]
    SENBA Y, KISHIMOTO H, MIURA T, et al. Upgrade of surface profiler for x-ray mirror at SPring-8[J]. Proceedings of SPIE,SPIE, 2014, 9206: 920604.
    [7]
    黄智超, 程建高, 李飞, 等. 同步辐射压弯镜重力补偿方法[J]. 强激光与粒子束,2018,30(8):085102. doi: 10.11884/HPLPB201830.180066

    HUANG ZH CH, CHENG J G, LI F, et al. Gravity compensation for bent mirror of synchrotron radiation[J]. High Power Laser and Particle Beams, 2018, 30(8): 085102. (in Chinese) doi: 10.11884/HPLPB201830.180066
    [8]
    李庆祥, 王东生, 李玉和. 现代精密仪器设计[M]. 北京: 清华大学出版社, 2004: 64-65.

    LI Q X, WANG D SH, LI Y H. Design of Modern Precision Instruments[M]. Beijing: Tsinghua University Press, 2004: 64-65. (in Chinese)
    [9]
    REED R R. A glass beam reference surface for quality control measurements[J]. International Journal of Mechanical Sciences, 1966, 8(11): 703-715. doi: 10.1016/0020-7403(66)90049-X
    [10]
    邵明振, 邵春雷, 卢启鹏, 等. 高功率TEA CO2激光器主机结构优化设计[J]. 发光学报,2013,34(3):388-393. doi: 10.3788/fgxb20133403.0388

    SHAO M ZH, SHAO CH L, LU Q P, et al. Design on mainframe of high power TEA CO2 laser and optimization[J]. Chinese Journal of Luminescence, 2013, 34(3): 388-393. (in Chinese) doi: 10.3788/fgxb20133403.0388
    [11]
    程显超. 同步辐射束线工程中相关技术的研究[D]. 合肥: 中国科学技术大学, 2016.

    CHENG X CH. Research on the beamline engineering technologies for the synchrotron radiation[D]. Hefei: University of Science and Technology of China, 2016. (in Chinese)
    [12]
    翟岩, 梅贵, 江帆, 等. Ф2020 mm口径空间红外相机主反射镜设计[J]. 发光学报,2018,39(8):1170-1176. doi: 10.3788/fgxb20183908.1170

    ZHAI Y, MEI G, JIANG F, et al. Ф2020 mm aperture space infrared camera main reflector design[J]. Chinese Journal of Luminescence, 2018, 39(8): 1170-1176. (in Chinese) doi: 10.3788/fgxb20183908.1170
    [13]
    NIJSSE G J P. Linear motion systems: a modular approach for improved straightness performance[D]. Delft: Delft University of Technology, 2001: 39-40.
  • Cited by

    Periodical cited type(3)

    1. 刘承麟,战俊彤,张肃,王超,付强,李英超,段锦,姜会林. 基于改进型Blinn遮蔽函数的目标材料表面可见光偏振反射研究. 中国光学(中英文). 2024(05): 1162-1174 .
    2. 陈宝华,吴泉英,唐运海,范君柳,陈晓翌,余浩墨,孙毅. 产生环形激光的光学系统设计. 中国光学(中英文). 2023(06): 1365-1375 .
    3. 卜和阳,虞林瑶,田浩南,王健. 中波红外成像系统冷反射抑制. 中国光学(中英文). 2023(06): 1414-1423 .

    Other cited types(3)

  • 加载中

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Figures(17)  / Tables(2)

    Article views(2347) PDF downloads(108) Cited by(6)
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return