Research on time domain parameters measurement of single pulse femtosecond laser
-
摘要: 钛宝石飞秒激光放大系统具有重复频率低、脉冲波形复杂等特点，为准确测量其峰值功率，提出对单次脉冲波形和脉冲宽度测量的需求。介绍了单脉冲飞秒激光时域波形和脉冲宽度的测量原理和测量方法，设计了基于频率分辨光学开关法的单脉冲飞秒激光时域参数测量装置。讨论了单脉冲飞秒激光时域参数测量校准面临的问题及解决方法。Abstract: The Ti sapphire femtosecond laser amplification system has the characteristics of low repetition frequency and complex pulse waveform. In order to accurately measure the peak power of the femtosecond laser amplification system, the measurement demand of single pulse waveform and pulse width was proposed. The measurement principle and method of single pulse femtosecond laser time domain waveform and pulse width were introduced, the time domain parameters measuring device of single pulse femtosecond laser based on frequency-resolved optical gating method was designed, and the technical problems and solutions in the calibration of time domain parameters measurement of single pulse femtosecond laser were discussed.
-
Keywords:
- femtosecond laser /
- pulse width /
- pulse waveform /
- single pulse measurement
-
引言
光学系统由多个光学元件组成，其中以各类透明元件的应用最为广泛。由于光学元件在生产过程中受到制造工艺的制约，成品会有一定的残余应力存在，致使其在使用过程中产生收缩、变形、翘曲或裂纹。不仅影响制品的表面形状和光学性能，而且直接关系到产品的使用寿命和产品安全，因此应力测量对于保证光学元件的强度和光学性能具有重要意义[1]。
应力的测量方法主要分为机械法与光学检测法。机械法如钻孔法、取条法，在测量残余应力时一般将部件切割出来，然后测量其应力的变化，从而求出残余应力，是一种间接测量手段，会对元件造成一定的损坏；光学检测法如光弹法、拉曼法及X射线衍射法等，这些方法均属于无损检测，其中光弹法是对应力导致的双折射相位差的测量，具有光学关联性更直接、测量精度更高的特点[2]。用光弹效应测量相位的方法可分为偏振光调制法[3-4]和相移法[5]。偏振光调制法一般采用补偿器、四分之一波片等偏振元件作补偿得到相位延迟量，往往导致测量时间长、测量次数多，很难按要求的精度来复现各个等差线及等倾线的位置；相移法又称光强法，利用待测样品的相位延迟量和出射光的光强之间的关系求解相位延迟。基于这一理论，近年来研究者们不断提出新的测量系统，主要包括基于五步相移法获得光弹特征的方法[6]，基于四步相移法并利用自准直光路系统测量玻璃应力的方法 [7]，基于六步相移法的红外光弹系统以计算硅片应力分布的方法[8-9]。这些方法都需要旋转检偏器和四分之一波片，难以实现对应力的实时监控。为了同时获取不同方向上的偏振信息，有研究者采用分振幅法[10]。但由于使用了多个光电传感器，使得捕获的偏振图像存在采集不同步和图像不匹配的问题。为解决上述问题，有学者提出可利用偏振相机测量应力[11]。偏振相机应用范围广泛，对于划痕缺陷检测、成像效果增强及应力检测等方面具有显著优势，因此这一技术得到了越来越多的关注，但是在应力测量领域的研究仍有不足。
针对上述研究情况，本文将像素偏振相机应用在应力测量系统中。首先利用Stokes矩阵和Mueller矩阵求解透明光学元件的相位延迟量分布及应力方向的计算公式，其次分析了影响系统测量精度的主要因素，并根据后续分析提出减小误差的方法，随后搭建出一套测量装置。实验过程简便，重复性好，最后通过四分之一波片的标定和车灯透镜应力分布的测量验证了方法的有效性。
1 测量原理与方案
1.1 光路模型
光路示意图如图1所示，光源沿z轴出发，经过一块透光轴为竖直y方向的起偏器成为线偏振光，然后垂直入射到待测样品表面，经透射后其偏振状态发生改变，最后光强信号由偏振相机接收。
由于偏振相机内置检偏方向为0°、45°、90°、135°，为方便与偏振相机组合使用，本文将采用四步相移法求解样品的应力分布。选用线偏光而非圆偏光的原因在于考虑到圆偏光的四步相移法在不涉及到旋转四分之一波片的前提下容易产生波片与光源波长不匹配的情况，从而影响测量结果[12]。
1.2 理论推导
光路原理图如图2所示，对整个光路运用Stokes矢量和Mueller矩阵表示法推导如下：
光强为I0的单色光光源S0为
$$ {{\boldsymbol{S}}_0} = \left[ {\begin{array}{*{20}{c}} {{{\boldsymbol{I}}_0}} \\ 0 \\ 0 \\ 0 \end{array}} \right] $$ 主轴方向与x轴成$\text{γ} $的起偏器Ρ0的Mueller矩阵为
$$ {{\boldsymbol{P}}_0} = \frac{1}{2}\left[ {\begin{array}{*{20}{c}} 1&{\cos 2\gamma }&{\sin 2\gamma }&0 \\ {\cos 2\gamma }&{{{\cos }^2}2\gamma }&{\sin 2\gamma \cos 2\gamma }&0 \\ {\sin 2\gamma }&{\sin 2\gamma \cos 2\gamma }&{{{\sin }^2}2\gamma }&0 \\ 0&0&0&0 \end{array}} \right] $$ 与x轴的夹角为θ、相位延迟量为φ的待测样品的Mueller矩阵为
$$ {\boldsymbol{M}} = \left[ {\begin{array}{*{20}{c}} 1&0&0&0 \\ 0&{{{\cos }^2}2\theta + {{\sin }^2}2\theta \cos \varphi }&{(1 - \cos \varphi )\sin 2\theta \cos 2\theta }&{ - \sin 2\theta \sin \varphi } \\ 0&{(1 - \cos \varphi )\sin 2\theta \cos 2\theta }&{{{\sin }^2}2\theta + {{\cos }^2}2\theta \cos \varphi }&{\cos 2\theta \sin \varphi } \\ 0&{\sin 2\theta \sin \varphi }&{ - \cos 2\theta \sin \varphi }&{\cos \varphi } \end{array}} \right] $$ 与x轴成β的检偏器分为0°、45°、90°、135° 等4个检偏方向，可表示为
$$ \begin{split} {{\boldsymbol{P}}_i} &= \frac{1}{2}\left[ {\begin{array}{*{20}{c}} 1&{\cos 2(\beta + {k_i})}&{\sin 2(\beta + {k_i})}&0 \\ {\cos 2(\beta + {k_i})}&{{{\cos }^2}2(\beta + {k_i})}&{\sin 2(\beta + {k_i})\cos 2(\beta + {k_i})}&0 \\ {\sin 2(\beta + {k_i})}&{\sin 2(\beta + i)\cos 2(\beta + {k_i})}&{{{\sin }^2}2(\beta + {k_i})}&0 \\ 0&0&0&0 \end{array}} \right], &(i = 1,2,3,4;{k_i} = 0^\circ ,45^\circ ,90^\circ ,135^\circ ) \end{split}$$ 则经过检偏器后的出射光的Stokes矢量可以由（2）式推出：
$$ {{\boldsymbol{I}}_i} = {{\boldsymbol{P}}_i}{\boldsymbol{M}}{{\boldsymbol{P}}_0}{{\boldsymbol{S}}_0},\;\;\;\;\;(i = 1,2,3,4) $$ (1) 根据Stokes矢量的定义可知，出射光矢量第1行的式子即为探测器CCD检测到的光强，则：
$$ {{\boldsymbol{I}}_1} = \frac{{{{\boldsymbol{I}}_0}}}{4}\left[ {1 + \cos 2(\beta - \gamma ) - 2{{\sin }^2}\frac{\varphi }{2}\sin 2(\theta - \gamma )\sin 2(\theta - \beta )} \right] $$ (2) 取$\gamma = 0$，$\beta = 0$,可得：
$$ {{\boldsymbol{I}}_1} = \frac{{{{\boldsymbol{I}}_0}}}{2} - \frac{{{{\boldsymbol{I}}_0}}}{2}{\sin ^2}\frac{\varphi }{2}{\sin ^2}(2\theta ) $$ (3) 同理可得${{\boldsymbol{I}}}_{2}、{{\boldsymbol{I}}}_{3}、{{\boldsymbol{I}}}_{4}$为
$$ \left\{ {\begin{array}{*{20}{l}} {{{\boldsymbol{I}}_1} = \dfrac{{{{\boldsymbol{I}}_0}}}{2} - \dfrac{{{{\boldsymbol{I}}_0}}}{2}{{\sin }^2}\dfrac{\varphi }{2}{{\sin }^2}(2\theta )} \\ {{{\boldsymbol{I}}_2} = \dfrac{{{{\boldsymbol{I}}_0}}}{4} + \dfrac{{{{\boldsymbol{I}}_0}}}{2}{{\sin }^2}\dfrac{\varphi }{2}\sin 2\theta \cos 2\theta } \\ {{{\boldsymbol{I}}_3} = \dfrac{{{{\boldsymbol{I}}_0}}}{2}{{\sin }^2}\dfrac{\varphi }{2}{{\sin }^2}(2\theta )} \\ {{{\boldsymbol{I}}_4} = \dfrac{{{{\boldsymbol{I}}_0}}}{4} - \dfrac{{{{\boldsymbol{I}}_0}}}{2}{{\sin }^2}\dfrac{\varphi }{2}\sin 2\theta \cos 2\theta } \end{array}} \right. $$ (4) 由（4）式可知：
$$ \varphi = 2\arcsin \sqrt {\frac{{{{\left( {{{\boldsymbol{I}}_3} + {{\boldsymbol{I}}_3}} \right)}^2} + {{({{\boldsymbol{I}}_2} - {{\boldsymbol{I}}_4})}^2}}}{{4{{\boldsymbol{I}}_3}({{\boldsymbol{I}}_1} + {{\boldsymbol{I}}_3})}}} $$ (5) $$ \theta = \frac{1}{2}\arctan \left(\frac{{{{\boldsymbol{I}}_3} - {{\boldsymbol{I}}_1} + {{\boldsymbol{I}}_2} + {{\boldsymbol{I}}_4}}}{{{{\boldsymbol{I}}_2} - {{\boldsymbol{I}}_4}}}\right) $$ (6) 从（5）式、（6）式可以看出，该方法在求解相位延迟量φ或应力方向角θ时，只与光强大小相关，二者互不影响，能够得到较好的结果。
1.3 误差模拟
影响精度的主要因素是器件的方位角误差，其中偏振相机的4个检偏方向可看作严格满足条件，因此误差主要在于起偏器与光轴夹角和偏振相机与光轴夹角不满足（3）式。由于γ=0°、β=0°，即起偏器、检偏器与x轴没有偏差角，是等式成立的先决条件。当此条件不成立时， I1产生的误差可表示为(2)式，其中cos2(β－γ)为常数项。空测时将I3调至光强最小处使得β－γ足够小，可令cos2(β－γ)=1－δ，此时sin2(β－γ)=(2δ－δ2)1/2，对I2、I3、I4做同样处理，可得到以下两种误差。
1） 方位角偏差给相位测量带来的误差
将无偏差角时测量的相位延迟量记作φ，有偏差角时测量的相位延迟量记作φ1，两者之差为Δφ, 经运算可得φ的最大误差Δφ/φ为
$$ \begin{gathered} \frac{{\Delta \varphi }}{\varphi } = 1 - {\left[ {\frac{{\sin 2(\theta - \gamma )}}{{\sin 2(\theta - \beta )}}} \right]^{1/2}} \\ \\ \end{gathered} $$ (7) 2） 方位角偏差给应力方向测量带来的误差
将无偏差角时测量的应力方向角记作θ，有偏差角时测量的应力方向角记作θ1，两者之差为Δθ, 经运算可得θ的最大误差Δθ/θ为
$$ \frac{{\Delta \theta }}{\theta } = 1 - \frac{{\tan 2\left( {\theta - \beta } \right)}}{{\tan 2\theta }} $$ (8) 2 实验流程
偏振相机是利用像素偏振片阵列对光束的偏振状态进行像素化的空间调制，使得4个方向上的偏振信息可以被同时捕获。使用最为广泛的偏振相机是由4个线偏振片单元、2×2周期排布形成的像素偏振片阵列组成[13]。从偏振片阵列中可分别提取出4个方向上的偏振信息，组成4幅偏振图像。图3为偏振相机原理示意图，以0°偏振方向为例，显示出偏振信息的获取过程。本文采用的偏振相机型号为MER-502-79U3M POL，分辨率为2 448 (H)×2 048 (V)，帧率为79 帧/s，传感器为2/3" Sony IMX250MZR，像素尺寸为3.45 μm×3.45 μm。
参考上述光路及分析计算，设计了一款实验系统，如图4（a）所示。该系统选择波长为620 nm的红光背光源，谱线宽度为5 nm。起偏器为线偏振片，波段为400 nm～700 nm，消光比>1000∶1。载物台用以承载样品，样品以车灯透镜为例。调节光学镜头可使样品的成像质量达到偏振相机的分辨率要求。
使用时先调节光源强度及曝光时间，使得偏振相机接收到的光强落在其线性相应区域；其次调节焦距及光圈大小，使样品能清晰成像；最后旋转起偏器在计算机中显示偏振相机捕获的4幅图像。当代表90°方向的I3光最暗、0°方向的I1光最亮，代表45°、135°方向的I2、I4光强相等时，可认为偏差角β－γ=0°。系统标定如图4（b）所示。
3 实验结果与分析
3.1 四分之一波片的测量及分析
选取一块波长为633 nm的四分之一波片并绘制误差曲线，当光源波长为620 nm时，该波片的相位延迟量为0.255 λ，即158.1 nm，精度为λ/100。图5显示了四分之一波片处于水平中位线上的各点相位延迟量的分布情况，滤波去噪后基本是一条水平线，显示出良好的均一性。
由上述误差的理论仿真及测量结果可知，应力值和应力方向的精度受到系统偏差角的影响且与样品本身的应力方向有关。将波片快轴方向在0°～90°内每5°旋转一次，测量误差如表1所示。多次测量的结果显示，精度与快轴方向有关。图6（a）和图6（b）模拟了相位延迟量与应力方向在不同偏差角下的误差曲线。当系统偏差角为2°时，模拟得到的理论误差曲线与实验误差曲线对比如图6（c）和图6（d）所示，二者相关系数大于0.85。
误差曲线显示出随着系统偏差角增大，误差会升高，偏差角存在的原因是空测时I3未调至最小。由于本例采用的是目测法，故I3仍有一定剩余量，使用光功率计调节可降低此误差。此外，注意到图6（a）和图6（b）中当快轴方向处于40°～60°范围内时，波片误差较小，在测量样品时可利用这一特性进一步减小误差。
表 1 四分之一波片的测量分析结果Table 1. Measurement analysis results of quarter wave plate快轴方向/（°） 延迟量/nm 误差/nm 30 151.53 6.57 35 160.41 2.31 45 158.96 0.86 55 159.64 1.54 60 152.65 5.45 3.2 透镜应力分布的测量及分析
透明元件的应力在各个区域不同，但整体分布是连续的，其应力分布特征影响着光学系统的性能。本文实际测量了一块LED车灯透镜，其尺寸规格为20 mm×15 mm×5 mm。图7显示了由（5）式计算出的透镜应力相位延迟量分布。从图7中可以看出，该透镜的相位延迟在中心菱形区域分布较均匀，选取菱形中心位置处的一点Ρ，该点φ=45°。计算其应力双折射值以此反映该透镜的应力大小。
该透镜由硼硅玻璃制成，其应力光学常数为3.8×10−6 mm2/N。取中心一小块区域测量多点，延迟量取平均值得46.06 nm，厚度为5 mm。由应力-双折射公式可知[14]，该中心区域应力双折射值为9.21 nm/mm、应力值为2.45 MPa。使用德国ILIS的Strain Scope应力测量仪检测的延迟量为44.24 nm，结果相差在2 nm范围以内。
图8（a）显示的是应力方向的全场测量结果，可以看出有明显的跳变，不满足应力连续性原则。造成这种现象的原因是应力方向角θ是由（6）式反正切函数解出，反正切函数的值域限制了应力方向角的范围。若要θ∈[0,π]，则需：
$$ \arctan \left(\frac{{{{\boldsymbol{I}}_3} - {{\boldsymbol{I}}_1} + {{\boldsymbol{I}}_2} + {{\boldsymbol{I}}_4}}}{{{{\boldsymbol{I}}_2} - {{\boldsymbol{I}}_4}}}\right) \subset \left[ {0,2\pi } \right] $$ (9) 根据参考文献 [15]中使用的符号规则（1），处理后得到真实的应力方向角全场分布图，满足应力连续性原则，如图8（b）所示。
4 结论
使用相移法确定透明元件内部应力的分布，相比于传统方法，不需要冻结应力和切片，符合实际应用的需求。本文基于应力双折射原理，提出了一种利用偏振相机实现四步相移可快速测量透明元件应力分布的测量方案。提出的理论误差曲线图符合实验误差曲线图，二者相似系数大于0.85。四分之一波片的测量结果显示，测量值为158.96 nm，误差小于1%。对一块车灯透镜的残余应力进行实测，得出透镜的全场应力延迟量分布，对测量得到的应力方向信息使用符号规则后得到了连续的应力方向分布图。与德国ILIS公司的应力测量仪进行了比较，结果相差在2 nm以内，验证了该方法的有效性。本文所述方法的主要优势在于无需旋转偏振元件即可实现全场应力测量，可适用于应力的实时测量。
-
![]()
图 2 单脉冲SHG-FROG信号光与光开关延迟示意图
Figure 2. Schematic diagram of delay time for single pulse SHG-FROG signal light and optical gating
![]()
图 4 单脉冲飞秒激光时域参数测量装置组成图
Figure 4. Composition block diagram of time domain parameters measuring device for single pulse femtosecond laser
![]()
图 5 热释电激光能量计探头硬件组成图
Figure 5. Composition block diagram of pyroelectric laser energy meter probe hardware
![]()
图 6 FROG脉冲波形测量单元组成图
Figure 6. Composition block diagram of FROG pulse waveform measurement unit
-
[1] MAINE P, STRICKLAND D, BADO P, et al. Generation of ultrahigh peak power pulses by chirped pulse amplification[J]. IEEE Journal of Quantum Electronics,1988,24(2):398-403. doi: 10.1109/3.137
[2] PERRY M D, PENNINGTON D, STUART B C, et al. Petawatt laser pulses[J]. Optics Letters,1999,24(3):160-162. doi: 10.1364/OL.24.000160
[3] WANG Z H, LIU C, SHEN Z W, et al. High-contrast 116 PW Ti: sapphire laser system combined with a doubled chirped-pulse amplification scheme and a femtosecond optical-parametric amplifier[J]. Optics Letters,2011,36(16):3194-3196. doi: 10.1364/OL.36.003194
[4] YU T J, LEE S K, SUNG J H, et al. Generation of high-contrast, 30 fs, 15 PW laser pulses from chirped-pulse amplification Ti: sapphire laser[J]. Optics Express,2012,20(10):10807-10815. doi: 10.1364/OE.20.010807
[5] CHU Y X, GAN Z B, LIANG X Y, et al. High-energy large-aperture Ti: sapphire amplifier for 5 PW laser pulses[J]. Optics Letters,2015,40(21):5011-5014. doi: 10.1364/OL.40.005011
[6] GAN Z B, YU L H, LI S, et al. 200 J high efficiency Ti: sapphire chirped pulse amplifier pumped by temporal dual-pulse[J]. Optics Express,2017,25(5):5169-5178. doi: 10.1364/OE.25.005169
[7] LI M, WANG Z H, TENG H, et al. Prospects for femtosecond ultrahigh intensity laser system towards Exawatt level[J]. Scientia Sinica (Physica, Mechanica & Astronomica),2018,48(2):024201.
[8] WANG Z H, FANG S B, TENG H, et al. Femtosecond laser user facility for application research on ultrafast science[J]. Chinese Physics B,2018,27(7):074204. doi: 10.1088/1674-1056/27/7/074204
[9] 张志刚. 飞秒激光技术[M]. 北京: 科学出版社, 2011. ZHANG Zhigang. Femtosecond laser technology[M]. Beijing: Science Press, 2011.
[10] DENG Yuqiang, SUN Qing, CAO Shiying, et al. Measurement of ultrashort optical pulses, modern metrology concerns[M]. Rijeka: InTech, 2012.
[11] 吴磊, 阴万宏, 俞兵, 等. 飞秒激光脉冲宽度和脉冲波形测试技术[J]. 应用光学,2019,40(2):291-299. WU Lei, YIN Wanhong, YU Bin, et al. Research on femto-second laser pulse width and pulse waveform measurement technology[J]. Journal of Applied Optics,2019,40(2):291-299.
[12] MONMAYRANT A, WEBER S, CHATEL B. A newcomer’s guide to ultrashort pulse shaping and characterization[J]. Journal of Physics B: Atomic, Molecular and Optical Physics,2010,43(10):103001. doi: 10.1088/0953-4075/43/10/103001
[13] GALLMANN L, SUTTER D H, MATUSCHEK N, et al. Techniques for the characterization of sub-10-fs optical pulses: a comparison[J]. Applied Physics B,2000,70(S1):S67-S75. doi: 10.1007/s003400000307
[14] TREBINO R. Frequency-resolved optical gating: the measurement of ultrashort laser pulses[M]. Boston, MA: Springer US, 2000.
[15] 王兆华, 魏志义, 滕浩, 等. 飞秒激光脉冲的谐波频率分辨光学开关法测量研究[J]. 物理学报,2003,52(2):362-366. doi: 10.3321/j.issn:1000-3290.2003.02.022 WANG Zhaohua, WEI Zhiyi, TENG Hao, et al. Measurement of femtosecond laser pulses using SHG frequency-resolved optical gating technique[J]. Acta Physica Sinica,2003,52(2):362-366. doi: 10.3321/j.issn:1000-3290.2003.02.022
[16] 姚尧. 超短脉冲单发频率分辨光学开关法测量技术研究[D]. 哈尔滨: 哈尔滨理工大学, 2012. YAO Yao. Research on ultrashort pulse measurement with GRENOUILLE[D]. Harbin: Harbin University of Science and Technology, 2012.
[17] 鲁琳琳, 贾豫东, 张晓青. 频率分辨光学开关法飞秒脉宽测量技术研究进展[J]. 激光与光电子学进展,2015,52(7):40-46. LU Linlin, JIA Yudong, ZHANG Xiaoqing. Research progress of frequency resolved optical gating femtosecond pulse width measurement technology[J]. Laser & Optoelectronics Progress,2015,52(7):40-46.
[18] IACONIS C, WALMSLEY I A. Self-referencing spectral interferometry for measuring ultrashort optical pulses[J]. IEEE Journal of Quantum Electronics,1999,35(4):501-509. doi: 10.1109/3.753654
[19] 王鹏, 王兆华, 魏志义, 等. 用SPIDER法测量飞秒激光脉冲的光谱相位[J]. 物理学报,2004,53(9):3004-3009. doi: 10.3321/j.issn:1000-3290.2004.09.036 WANG Peng, WANG Zhaohua, WEI Zhiyi, et al. Measurement of spectral phase of femotosecond laser pulse using SPIDER technique[J]. Acta Physica Sinica,2004,53(9):3004-3009. doi: 10.3321/j.issn:1000-3290.2004.09.036
[20] 何铁英, 柴路, 王清月, 等. SPIDER光谱相位干涉议的实现及相关理论分析[J]. 中国激光,2004,31(s1):152-154. HE Tieying, CHAI Lu, WANG Qingyue, et al. Realization and analysis for a SPIDER apparatus[J]. Chinese Journal of Lasers,2004,31(s1):152-154.
[21] 文汝红, 刘天夫. 用SPIDER法还原飞秒脉冲振幅和相位误差[J]. 中国计量学院学报,2006,17(2):114-118. WEN Ruhong, LIU Tianfu. Error of spectral intensity and phase of femtosecond laser pulse using SPIDER technique[J]. Journal of China University of Metrology,2006,17(2):114-118.
[22] 文汝红, 刘天夫. 用SPIDER法重构飞秒脉冲位相中参数的选择[J]. 应用光学,2006,27(4):340-343. WEN Ruhong, LIU Tianfu. Optimization of parameters on reconstruction of phase of femtosecond laser pulse with SPIDER[J]. Journal of Applied Optics,2006,27(4):340-343.
[23] 张素侠, 李方家, 刘军. 自参考光谱干涉法测量飞秒脉冲的研究进展[J]. 激光与光电子学进展,2014,51(1):010001. ZHANG Suxia, LI Fangjia, LIU Jun. Research progress of self-referenced spectral interferometry on the characterization of femtosecond pulses[J]. Laser & Optoelectronics Progress,2014,51(1):010001.
-
期刊类型引用(0)
其他类型引用(1)
下载:
计量
- 文章访问数: 3583
- HTML全文浏览量: 669
- PDF下载量: 56
- 被引次数: 1
陕公网安备 61011302001501号 