用于应力双折射分布测量的多构型扫描系统设计

岐临凯, 李克武, 李坤钰, 王志斌, 李晋华

岐临凯, 李克武, 李坤钰, 王志斌, 李晋华. 用于应力双折射分布测量的多构型扫描系统设计[J]. 应用光学, 2024, 45(3): 608-615. DOI: 10.5768/JAO202445.0301003
引用本文: 岐临凯, 李克武, 李坤钰, 王志斌, 李晋华. 用于应力双折射分布测量的多构型扫描系统设计[J]. 应用光学, 2024, 45(3): 608-615. DOI: 10.5768/JAO202445.0301003
QI Linkai, LI Kewu, LI Kunyu, WANG Zhibin, LI Jinhua. Design of multi-configuration scanning system for stress birefringence distribution measurement[J]. Journal of Applied Optics, 2024, 45(3): 608-615. DOI: 10.5768/JAO202445.0301003
Citation: QI Linkai, LI Kewu, LI Kunyu, WANG Zhibin, LI Jinhua. Design of multi-configuration scanning system for stress birefringence distribution measurement[J]. Journal of Applied Optics, 2024, 45(3): 608-615. DOI: 10.5768/JAO202445.0301003

用于应力双折射分布测量的多构型扫描系统设计

基金项目: 国家自然科学基金(62205310);山西省自然科学基金(201901D211234)
详细信息
    作者简介:

    岐临凯(1997—),男,硕士,主要从事光电信息检测研究。E-mail:858515048@qq.com

    通讯作者:

    王志斌(1966—),男,硕士,教授,主要从事光电信息物理、测试计量技术及仪器、光电子技术与仪器、光学工程研究。E-mail:wangzhibin@nuc.edu.cn

  • 中图分类号: TN206

Design of multi-configuration scanning system for stress birefringence distribution measurement

  • 摘要:

    为了实现对不同构型光学样品进行应力双折射分布测量,在使用双弹光调制方法的基础上,设计了一种应用于双折射分布测量的多构型扫描系统。该系统在保证测量的高分辨率的同时,通过保持激光器静止,同时使样品进行快速移动,提高了测量精度与广度。在样品测量方面,采用633 nm(1/4)玻片测试,测试相对误差的范围为0.79%~0.95%,波动范围为0.12 nm,标准差为0.035 2;采用BK7玻璃样品测试,波动范围为0.25 nm,标准差为0.038 9。在扫描精度方面,连续扫描精度误差不超过0.05 mm,连续寸动扫描精度误差不超过0.009 mm。对比实验结果可得出,该多构型扫描系统可有效解决对样品任一区域实现高精度应力双折射测量的问题。

    Abstract:

    In order to measure the stress birefringence distribution of optical samples with different configurations, a multi-configuration scanning system for birefringence distribution measurement was designed based on the double elastic light modulation method. The system not only ensured the high resolution of the measurement, but improved the measurement accuracy and breadth by keeping the laser stationary and making the sample move quickly. In terms of sample measurement, the 633 nm (1/4) glass was adopted to test, the relative error of the detection results was 0.79%~0.95%, with a variation range of 0.12 nm, and the standard deviation was 0.035 2. The variation range of the BK7 glass sample was 0.25 nm, and the standard deviation of the experimental results was 0.038 9. For scanning accuracy, the continuous scanning accuracy error does not exceed 0.05 mm, and the continuous inching scanning accuracy error is less than 0.009 mm. According to the above two experimental results, it can be concluded that the multi-configuration scanning system can effectively solve the problem of high precision stress birefringence measurement in any area of the sample.

  • 随着光学镀膜技术的不断发展进步,广泛应用于高能激光系统、引力波探测、激光陀螺和腔衰荡光谱测试等的低损耗、高反射率光学元件性能已大大提高,其反射率已达到99.9%以上,有的甚至达到99.99%以上。基于分光光度技术的常规检测方法有光透射法、多次反射法等方法,但由于测量精度只能达到0.1%左右[1],无法准确测量高于99.9%光学元件的反射率。目前能精确测量反射率高于99.9%的唯一方法是光腔衰荡技术[2]。光腔衰荡技术是一种基于高精细度谐振腔的高灵敏度探测技术,利用光在谐振腔内来回振荡,其出射光的强度随时间呈指数形式衰减,通过测量衰荡时间来确定待测镜的反射率,是一种绝对测量方法,不需要标定。

    光腔衰荡技术的测试原理为用合适波长的激光器发出激光,通过耦合光路进入2个凹面镜组成的谐振腔,激光在衰荡腔中经过反射会形成衰荡,用探测器探测衰荡腔中光强信号,即可监视腔内光强的损耗情况。通过理论分析,光强信号的衰减过程符合一次指数函数衰减规律,衰减速率的快慢直接反映光腔损耗。得到初始腔的损耗后,将待测光学元件放入衰荡腔内,就构成了测试腔。光腔衰荡法是通过测量激光在谐振腔中的衰荡时间来确定腔的损耗或者待测腔镜的反射率。一般采用直腔和折叠腔相结合的方法,根据原理公式计算得到腔镜和待测光学元件的反射率。

    理论上来说,对于光腔衰荡技术,初始腔的腔内激光损耗越小,测量精度就越高。减小腔内损耗的方法包括升级系统硬件和优化腔结构。其中,硬件升级主要依赖于超高反射率腔镜(例如R=99.999 8%)的使用。但是在实际调试过程中,机械结构的组装加工都会带来不可避免的误差,会导致实际的旋转轴与理想位置不重合,从而引起扫描系统在扫描过程中腔长发生变化,因此就必须对扫描系统精度进行分析计算,根据系统的具体指标提出对扫描系统的要求。本文通过理论计算分析,结合实际测得的光腔衰荡数据得到对应的反射率,通过对比检验误差结果的影响判断是否准确。为提升光腔衰荡法测光学元件反射率提供了分析支撑。

    光腔衰荡技术依据所采用的光源不同,分脉冲光衰荡技术[3]、窄谱连续波光腔衰荡技术[4](NBCW-CRD)、相移光腔衰荡技术[5](PS-CRD)和自混合光腔衰荡技术[6](OF-CRD)等,在这些方法中,脉冲光衰荡技术的优点是装置简单[7-9],物理模型最为直观。

    脉冲光衰荡技术的测量原理如图1所示,脉冲激光光束透过耦合镜组将脉冲激光入射到衰荡腔,高反射率反射镜M1和M2构成一个光学谐振腔[10],当一束单色脉冲激光沿着光轴入射到光腔内,激光脉冲在2个腔镜之间往返振荡,每经过一次循环透射出部分光,M2后面的探测器接收脉冲信号。由于脉冲光入射到高反腔内时,只有很少的光进入光腔内衰荡,每次振荡透过镜片的光就更少,透射光的强度与输出镜的反射率有关,是时间的指数函数[11-13]

    图  1  脉冲光腔衰荡法测量原理图
    Figure  1.  Schematic diagram of measurement for pulse cavity ring-down method

    当激光脉冲进入衰荡腔后,光腔输出信号可表示为

    $$ I(t) = {I_0}\exp \Bigg(\dfrac{{ - t}}{\tau }\Bigg) $$ (1)

    直形衰荡光腔的衰荡时间τ与反射率R的关系为

    $$ \tau = \dfrac{{nL}}{{c(\alpha L - \ln R)}} $$ (2)

    式中:L为腔长;n为腔内介质的折射率;α为腔内吸收系数。当腔内介质为空气时,折射率n近似为1;c为光速;$\alpha $为损耗系数,当腔内无吸收介质时$\alpha $=0。

    将待测反射镜置于折叠腔的折点处,以一定角度将直腔改变为折叠腔,如图2所示。

    图  2  折叠衰荡腔示意图
    Figure  2.  Schematic diagram of folding ring-down cavity

    折叠衰荡光腔的衰荡时间$t$与反射率Rx的关系为

    $$ \dfrac{1}{t} = \dfrac{1}{{{t_0}}} - \dfrac{{c\ln {R_x}}}{L} $$ (3)

    脉冲光衰荡方法的测量结果不受光源的光强波动的影响,测量精度高,可进行不同入射角度的反射率测量,测量精度可达到${10^{ - 5}}$量级以上。这种测量法适用于高反射率元件测量,若元件的反射率低于97%,激光脉冲信号衰荡时间很短,从而影响测量精度。因此,该测量法更适合99.5%以上高反射率镜片的测量,且镜片的反射率越高,测量精度越高[14]

    R1R2分别为2个腔镜的反射率。

    若衰荡腔内介质为空气,忽略腔内空气吸收,当谐振腔镜反射率(R1=R2)相等时,则直腔结构谐振腔镜反射率为

    $$ R = \sqrt {{R_1}{R_2}} = 1 - \dfrac{L}{{c\tau }} $$ (4)

    对折叠腔,则待测反射镜的反射率为

    $$ {R_x} = \exp \Bigg[ - \Bigg(\dfrac{1}{t} - \dfrac{1}{{{t_0}}}\Bigg)\dfrac{L}{c}\Bigg] $$ (5)

    因此,本方法通过测量光束在衰荡腔内的衰荡时间,进而得到待测腔镜的反射率。

    系统采用脉冲光腔衰荡技术,其系统装置如图3所示,光源为1 064 nm固体激光器,TEM00模输出,光斑大小约2 mm,能量约2 mJ。激光光束经2个45°反射镜对光束进行折转后,经光束整形单元对光束变焦后,与衰荡腔本征模进行匹配,最后入射到由平凹反射腔镜 M1和M2(曲率半径为−1 m)构成的初始光学谐振腔,激光在腔中来回反射而形成振荡,在腔镜 M2后用光电探测器探测透射光,采集卡对测试数据进行记录并保存。实验装置放置在洁净度万级的超净实验室环境,空调净化系统始终保持开启状态。通过对初始腔的光腔衰荡信号进行单指数拟合,得到初始腔衰荡时间t0。将待测反射光学元件Rx按使用角度插入衰荡腔空腔中,形成测试腔,拟合得到测试腔衰荡时间t,并由(5)式计算得到待测反射光学元件的反射率。

    图  3  光腔衰荡法测试系统
    Figure  3.  Test system of cavity ring-down method

    在本系统中,腔镜M2与旋转扫描平台在结构上采用共轴的设计,目的是减小扫描结构的系统误差[15]

    由(5)式可知,待测镜的反射率主要影响因素有腔长L和衰荡时间t,而光学衰荡腔的影响因素包括腔长失调、模式匹配、腔镜倾斜。衰荡时间的影响因素包括探测器孔径大小、探测器响应时间、数据拟合算法等因素[16]。上述因素的影响规律相关文献已经进行描述分析,若采用大口径旋转扫描平台对大尺寸光学元件的反射率进行测量,扫描结构如图4所示。扫描系统涉及2个旋转轴、2个平移轴,扫描系统的系统精度对反射率测量精度的影响规律目前尚未见相关报道[17],因此,本文主要从该方面进行分析。

    图  4  大口径光学元件精密扫描系统结构图
    Figure  4.  Structure diagram of precision scanning system for large-aperture optical element

    图5所示,从理论上讲,当入射激光的初始传播方向保持不变时,在扫描机构扫描过程中,入射激光与待测光学元件交汇点的位置应保持不变,即激光在待测镜上的入射点的空间位置恒定;然而实际上,由于机械结构的加工、装调引入的各类误差的累积,不可避免地会导致实际的旋转轴与理想位置不重合,从而导致扫描系统在扫描过程中腔长会发生变化,从而增大测量的误差。因此,必须对扫描系统的精度进行分析计算,根据系统的具体指标提出对扫描系统精度的设计/调试要求。

    图  5  光学偏差示意图
    Figure  5.  Schematic diagram of optical deviation

    由于腔镜M2与旋转扫描平台在结构上采用共轴的设计,通过系统设计及装调,能够保证2个旋转轴的旋转中心偏差在10 μm以内,若衰荡腔的腔长为1 m,待测反射镜放置在衰荡腔中心位置,则腔镜(M2)的旋转半径为r=0.5 m,所对应的角度偏移在20 μrad以内,约为人眼的最小可辨别失调角580 μrad的3.5%,因此可以忽略不计。

    当旋转轴与理想旋转轴不重合时,如图5所示,实际旋转轴为点(ab),理想旋转轴为o点,反射镜旋转$ \alpha $角度,经反射镜的实际反射光线与理想反射光线是平行的,当腔镜M2的旋转半径为r时,引起的衰荡腔长变化量为$ a + \sqrt {{r^2} - {b^2}} - r $,入射角度变化量为$\arcsin \Bigg(\dfrac{b}{r}\Bigg)$。由上述2个公式可知,腔长变化量和入射角度变化量均与扫描系统的旋转角度无关,只与实际旋转中心的位置有关,由于a、b的误差量级在mm级,当r为0.5 m时,$ a + \sqrt {{r^2} - {b^2}} - r \approx a $,即腔长变化量只与a相关。由文献[18]衰荡光腔中腔镜倾斜分析可知,若b的值小于0.29 mm时,由入射角度变化引入的测量相对误差为10−6;腔长失调量约为a值;若其量值在mm量级,可通过调整反射镜M2的前后位置来改变衰荡腔长来修正,通过示波器监视衰荡信号,对衰荡腔腔长调节实现精确调节,减小腔长的变化对测量误差的影响。

    通过对旋转扫描平台进行精密调试,获得了在大口径光学元件在旋转时不同角度下的衰荡曲线,通过最小二乘法对衰荡曲线进行拟合,得到衰荡时间,从而计算得到反射率,对其单点进行N次(N不小于6)测量得到$ {R_i} $,计算出测量的平均值$ \mathop R\limits^ - $=0.999 815 3,即可求出标准偏差$ s(R) $=1.032 8×$10_{}^{ - 6}$ 。如表1所示。

    表  1  光腔衰荡时间和对应的反射率
    Table  1.  Cavity ring-down time and corresponding reflectivities
    实验次数123456
    τ18.01818.01218.08317.99218.12618.051
    R0.999 8150.999 8140.999 8150.999 8170.999 8160.999 815
    下载: 导出CSV 
    | 显示表格

    把光电探测器采集得到的时间与光照强度的数据通过一次指数函数进行拟合,得到拟合数据曲线如图6所示。通过原始数据点与拟合图像的关系,计算得到对应函数关系式的统计值R2的均值为0.994 3,说明一次指数函数拟合得到的值拟合度较为准确。测得的待测镜反射率为0.999 815,将待测镜旋转一定的角度,得到的反射率为0.999 817,可知所设计的旋转扫描平台对测试结果影响很小,可忽略不计。

    图  6  一次指数函数拟合曲线图
    Figure  6.  Fitting curve of first-order exponential function

    本文从腔衰荡大口径光学元件反射率测试系统的工程化角度,分析了旋转扫描单元的设计及调试误差对系统测试误差的影响。在垂直于光的传播方向上,实际旋转轴与理想旋转轴的水平偏差在0.29 mm时,测量相对误差为10−6;在光的传播方向上,实际旋转轴与理想旋转轴的水平偏差在1 mm以内,可通过调节衰荡腔腔镜的前后距离,实现实际腔长与理论腔长相近,从而降低系统的测试误差。通过试验测试,所设计的测试系统的测量误差在10−6以内,能够满足较高测试系统的测试精度需求。

  • 图  1   双折射测量系统结构图

    Figure  1.   Structure diagram of birefringence measurement system

    图  2   系统整体设计框图

    Figure  2.   Block diagram of overall system design

    图  3   扫描模块流程图

    Figure  3.   Flow chart of scanning module

    图  4   数据处理模块流程图

    Figure  4.   Flow chart of data processing module

    图  5   人机交互界面

    Figure  5.   Man-machine interaction interface

    图  6   x轴连续扫描误差

    Figure  6.   X-axis continuous scanning error

    图  7   y轴连续扫描误差

    Figure  7.   Y-axis continuous scanning error

    图  8   x轴连续寸动扫描误差

    Figure  8.   X-axis continuous inching scanning error

    图  9   y轴连续寸动扫描误差

    Figure  9.   Y-axis continuous inching scanning error

    图  10   1/4玻片延迟量重复性测试结果

    Figure  10.   Test results of 1/4 slide retardation repeatability

    图  11   BK7玻璃延迟量重复性测试结果

    Figure  11.   Test results of BK7 glass retardation repeatability

    图  12   1/4玻片快轴方位角重复性测试结果

    Figure  12.   Test results of 1/4 slide fast axis azimuth repeatability

    图  13   BK7玻璃快轴方位角重复性测试结果

    Figure  13.   Test results of BK7 glass fast axis azimuth repeatability

    图  14   1/4玻片扫描结果

    Figure  14.   Scanning results of 1/4 slide

    图  15   BK7玻璃扫描结果

    Figure  15.   Scanning results of BK7 glass

    表  1   不同样品的重复性测量结果

    Table  1   Repeatability measurement results of different samples

    样品测量数据类型波动范围标准差
    633 nm(1/4)玻片延迟量0.12 nm0.035 2 nm
    快轴方位角0.13°0.022 5°
    BK7玻璃延迟量0.25 nm0.038 9 nm
    快轴方位角0.13°0.037 0°
    下载: 导出CSV
  • [1] 肖昊苏, 张运强, 范志刚, 等. 偏振干涉法测量晶体应力双折射精度分析[J]. 红外与激光工程,2011,40(2):271-276.

    XIAO Haosu, ZHANG Yunqiang, FAN Zhigang, et al. Precision analysis of polarization interference method for measuring stress birefringence of crystal[J]. Infrared and Laser Engineering,2011,40(2):271-276.

    [2] 左芬, 陈磊, 徐晨. 1/4波片相位延迟分布的动态测量[J]. 光子学报,2008,37(11):2296-2299.

    ZUO Fen, CHEN Lei, XU Chen. Dynamic phase-shifting interferometry for full field retardation distributing of quarter wave plate[J]. Acta Photonica Sinica,2008,37(11):2296-2299.

    [3] 李金鹏, 陈磊, 乌兰图雅, 等. 基于移相算法的玻璃应力延迟量快速测量方法[J]. 光学学报,2013,33(6):125-129.

    LI Jinpeng, CHEN Lei, WU Lantuya, et al. A method based on phase shifting algorithm for fast measurement of stress retardation of optical glass[J]. Acta Optica Sinica,2013,33(6):125-129.

    [4] 王金伟, 李克武, 景宁, 等. 波片相位延迟量测量和快轴标定系统[J]. 红外与激光工程,2019,48(2):211-216.

    WANG Jinwei, LI Kewu, JING Ning, et al. Phase retardation measurement and fast axis calibration system for wave plate[J]. Infrared and Laser Engineering,2019,48(2):211-216.

    [5] 王雷, 许荣国, 阴万宏, 等. 光学元件残余偏振测量方法与装置[J]. 应用光学,2021,42(6):1080-1085. doi: 10.5768/JAO202142.0603004

    WANG Lei, XU Rongguo, YIN Wanhong, et al. Testing method and equipment for residual polarization of optical lens[J]. Journal of Applied Optics,2021,42(6):1080-1085. doi: 10.5768/JAO202142.0603004

    [6] 张瑞, 陈媛媛, 景宁, 等. 基于双弹光差频调制的中红外波片相位延迟高精度测量[J]. 光学学报,2019,39(3):184-190.

    ZHANG Rui, CHEN Yuanyuan, JING Ning, et al. High-precision measurement of mid-infrared waveplate phase retardation based on dual photoelastic difference frequency modulation[J]. Acta Optica Sinica,2019,39(3):184-190.

    [7] 古兆兵, 郝淑杰, 王雷, 等. 一种基于旋转1/4波片法的激光偏振度测量仪[J]. 应用光学,2018,39(6):936-941.

    GU Zhaobing, HAO Shujie, WANG Lei, et al. Laser polarimeter based on rotary 1/4 waveplate[J]. Journal of Applied Optics,2018,39(6):936-941.

    [8] 宋连科, 牛明生, 韩培高, 等. 基于电光补偿测量的光相位延迟拓展测量法[J]. 光学技术,2017,43(3):199-202.

    SONG Lianke, NIU Mingsheng, HAN Peigao, et al. Optical phase delay expansion measuring method based on electro-optical compensation measurement[J]. Optical Technique,2017,43(3):199-202.

    [9] 李继扬, 谈宜东, 吴季, 等. 基于激光回馈效应的液晶双折射特性测量[J]. 红外与激光工程,2017,46(3):114-119.

    LI Jiyang, TAN Yidong, WU Ji, et al. Birefringence measurement of liquid crystals based on laser feedback effect[J]. Infrared and Laser Engineering,2017,46(3):114-119.

    [10] 肖石磊, 李斌成. 光学元件残余应力无损检测技术概述[J]. 光电工程,2020,47(8):51-61.

    XIAO Shilei, LI Bincheng. Residual stress measurement methods of optics[J]. Opto-Electronic Engineering,2020,47(8):51-61.

    [11] 马响, 邓勇, 张书练. 激光回馈半钢化玻璃应力双折射测量技术[J]. 激光技术,2020,44(3):371-376. doi: 10.7510/jgjs.issn.1001-3806.2020.03.018

    MA Xiang, DENG Yong, ZHANG Shulian. Stress birefringence measurement technology of heat strengthened glass based on laser feedback[J]. Laser Technology,2020,44(3):371-376. doi: 10.7510/jgjs.issn.1001-3806.2020.03.018

    [12] 李克武, 王黎明, 王志斌, 等. 弹光调制和电光调制联合测剩余双折射[J]. 中国激光,2016,43(5):208-215.

    LI Kewu, WANG Liming, WANG Zhibin, et al. Measurement of residual birefringence combined photo-elastic modulation with electro-optic modulation[J]. Chinese Journal of Lasers,2016,43(5):208-215.

    [13] 王爽, 韩燮, 李晓, 等. 弹光调制测椭偏参量的数字锁相数据处理[J]. 光学 精密工程,2018,26(6):1314-1321. doi: 10.3788/OPE.20182606.1314

    WANG Shuang, HAN Xie, LI Xiao, et al. Digital phase-locked data processing for ellipsometric parameter measurements based on photoelastic modulation[J]. Optics and Precision Engineering,2018,26(6):1314-1321. doi: 10.3788/OPE.20182606.1314

    [14] 刘梓良, 王志斌, 李克武, 等. 一种弹光调制器精确定标微系统设计[J]. 光学技术,2022,48(1):8-13.

    LIU Ziliang, WANG Zhibin, LI Kewu, et al. A precise calibration microsystem for photoelastic modulator[J]. Optical Technique,2022,48(1):8-13.

    [15] 杜洽伟, 李克武, 王志斌. 基于弹光调制的双折射测试系统数据处理研究[J]. 光学技术,2023,49(1):64-69.

    DU Qiawei, LI Kewu, WANG Zhibin. Study on data processing of birefringence measurement system based on photoelastic modulation[J]. Optical Technique,2023,49(1):64-69.

  • 期刊类型引用(12)

    1. 袁利,姜甜甜,魏春岭,杨孟飞. 空间控制技术发展与展望. 自动化学报. 2023(03): 476-493 . 百度学术
    2. 覃振鹏,徐金敏,吕德深. 基于三维光学测量的光纤微缺陷无损检测研究. 激光杂志. 2023(09): 233-237 . 百度学术
    3. 张志鹏,邵学君,庞庆. 基于三维激光扫描的铁路罐车点云优化处理方法. 应用光学. 2022(03): 460-465 . 本站查看
    4. 万青. 基于数据挖掘的光纤通信网络异常数据检测. 中国新通信. 2022(18): 4-6 . 百度学术
    5. 蒋心学,唐飞笼,杨迪. 机载三维激光点云数据分类数学模型. 激光杂志. 2021(07): 142-146 . 百度学术
    6. 王宪伦,丁文壮,孙旭祥. 场景点云中小孔洞边界提取算法. 机械制造与自动化. 2021(05): 39-41+52 . 百度学术
    7. 王洪. 基于关联规则算法的不动产边界三维测量方法. 经纬天地. 2021(06): 14-17 . 百度学术
    8. 林婷婷,江晟,李荣华,葛研军,周颖. 非合作目标视觉位姿测量与地面验证方法. 大连交通大学学报. 2020(03): 34-40 . 百度学术
    9. 李敏,金晶. 基于VR技术的舰船航行环境快速重构研究. 舰船科学技术. 2020(14): 79-81 . 百度学术
    10. 田九玲,杨永菊. 基于Rodrigues参数交互模型航天器相对位姿测量技术. 计算机测量与控制. 2020(09): 19-22+33 . 百度学术
    11. 陈立胜,于丽娜. 基于虚拟光学的室内三维空间视觉合理性设计. 激光杂志. 2020(11): 183-187 . 百度学术
    12. 马莉莉,刘江平. 基于数据挖掘的光纤通信网络异常数据检测研究. 应用光学. 2020(06): 1305-1310 . 本站查看

    其他类型引用(5)

图(15)  /  表(1)
计量
  • 文章访问数:  78
  • HTML全文浏览量:  13
  • PDF下载量:  29
  • 被引次数: 17
出版历程
  • 收稿日期:  2023-07-26
  • 修回日期:  2023-09-18
  • 网络出版日期:  2024-04-02
  • 刊出日期:  2024-05-14

目录

/

返回文章
返回