-
目前常用光纤埋入槽型有四种,分别为矩型槽、U型槽、梯型槽和三角型槽,考虑到实际制造的难易程度,选用矩形槽结构建立了有/无涂覆层光纤埋入基板结构的数值模型,这两种光纤埋入结构如图1所示。
图 1 (a)无涂覆层光纤埋入结构示意图;(b) 有涂覆层光纤埋入结构示意图
Figure 1. (a) Structure diagram of uncoated optical fibers embedding; (b) Structure diagram of coated optical fibers embedding
对于无涂覆层多模裸光纤埋入结构,其光纤纤芯半径为62.5 μm,半槽宽度值为65 μm,槽间距为250 μm,上下聚酰亚胺(Polyimide,PI)保护层厚度为25 μm,中间PI基板厚度为100 μm,上下填充胶层厚度为50 μm,上下导电层铜箔厚度为18 μm,挠性基板宽度为2 mm。对于有涂覆层多模裸光纤埋入结构,其光纤纤芯半径为62.5 μm,涂覆层为聚酰亚胺材料,厚度为15 μm,能够长期承受300 ℃高温,兼容层压高温工艺;半槽宽度值为80 μm,槽间距为250 μm,上下PI保护层厚度为25 μm,中间PI基板厚度为100 μm,上下填充胶层厚度为50 μm,上下导电层铜箔厚度为18 μm,挠性基板宽度为2 mm。各层材料参数如表1所示[3,15,17]。
网格的剖分对于有限元的求解非常关键,根据挠性基板不同的结构采用不同的网格形状,可在提高仿真精度的同时加快计算结果收敛。对PI保护层与铜箔等形状规则的几何结构采用映射法来划分网格;对于光纤、中间PI基板层及填充胶结构单元采用自由三角形网格来划分,对最大网格单元和最小网格单元尺寸进行控制,最小单元为0.06 μm,最大单元尺寸为10 μm。统计三角形网格数为13141,四边形网格数为10970,平均网格质量为0.8423,求解自由度数为190388。图2为有涂覆层光纤埋入基板结构网格划分后的全局视图。
表 1 材料属性
Table 1. Properties of materials
Material Elastic
modulus/GPaPoisson's
ratioThermal expansion
coefficient/℃Density/
kg∙m−3Heat capacity/
J∙(Kg∙℃)−1Thermal conductivity/
W∙(m∙℃)−1Copper layer 110 0.33 18×10−6 8940 384 398 Protective layer 3.20 0.36 22×10−6 1420 1090 0.12 Silica fiber 71.9 0.16 0.55×10−6 2200 745 1.50 Filled resin 7.84 0.30 27×10−6 970 1600 0.21 Coating layer 3.20 0.36 22×10−6 1420 1090 0.12 -
在实际生产制造过程中,需先将挠性基板固定,再对挠性基板施加时间、温度、压力载荷,故在仿真中对基板下表面施加固定约束,对基板上下表面施加温度载荷,对基板上表面施加压力载荷,如图3所示。图4为层压工艺曲线[17],首先对挠性光电基板从室温(25 ℃)上升到130 ℃,此时压力为0 MPa;然后在15 min内将压力增加到1 MPa (温度为130 ℃);再将温度升到180 ℃并加压到2 MPa保压0.5 h;接着缓慢升温到220 ℃,在2 MPa的环境下保温固化2 h;最后卸载压力和温度。升温速率和降温速率均为5 ℃/min。
-
图5所示为有/无涂覆层光纤埋入双面覆铜矩型槽内受到的应力。从图5(a)可知,无涂覆层光纤的最大应力约为125.04 MPa;从图5(b)可知,有涂覆层光纤的最大应力约为81.30 MPa;图5(c)为层压工艺下涂覆层应力分布图,最大应力约为93.09 MPa。
图 5 (a) 无涂覆层光纤应力分布;(b) 有涂覆层光纤应力分布;(c) 涂覆层应力分布
Figure 5. (a) Stress distribution of uncoated optical fibers; (b) Stress distribution of coated optical fibers; (c) Stress distribution of coatings
图6为中间PI层厚度对埋入光纤应力的影响。由图6可知,随着PI层厚度值从90 μm增大至160 μm,埋入光纤所受应力不断减小,无涂覆层光纤所受应力从129.72 MPa减小至116.80 MPa,有涂覆层光纤所受应力从89.47 MPa减小至52.02 MPa。当PI层厚度值超过140 μm后,有/无涂覆层光纤所受应力趋于稳定。
图 6 中间PI层厚度对光纤应力的影响
Figure 6. The influence of the thickness of the middle PI layer on optical fibers stress
在实际制造过程中,很难保证光纤与铜层槽底直接接触,通常会有填充胶在层压工艺下流动至槽底将光纤与槽底隔开。图7所示为光纤底部填充胶厚度对光纤应力的影响。由图7可知,随着填充胶厚度由0 μm增加到10 μm,有涂覆层光纤应力呈现不断增大的趋势,从81.30 MPa增至84.52 MPa,但趋势较缓;对于无涂覆层光纤,填充胶厚度由0 μm 变为1 μm时,光纤所受应力从125.04 MPa骤减至86.82 MPa,由1 μm增加到10 μm时光纤所受应力不断增大至93.53 MPa。
-
FEOPCB制造主要分为两步:第一步是在双面覆铜PI基材层上制作高精度光纤矩形槽,以实现埋入光纤的定位;第二步是高温层压,将布置好有涂覆层光纤的PI基材层、半固化片和PI保护层进行层叠定位,然后放到热压设备的层压窗口,加热加压。光纤矩形定位槽与半固化片的使用可降低埋入光纤在层压过程中所受热应力和偏移量,提升FEOPCB的使用寿命和耦合效率。
-
由于FEOPCB采用挠性PI基材,无法使用机械方式制作的矩形槽。利用激光的光热烧蚀机理,在双面覆铜PI基材上制作矩形槽结构。制造厂商提供的PI基材尺寸参数如下:基材层厚度约为138 μm (铜层厚度约为18 μm,中间PI层厚度约为102 μm)。首先通过传统的电路板图形刻蚀工艺,在顶层铜上制作出精确矩形槽尺寸宽度Fpw,间距为P,呈现出激光烧蚀PI层的窗口。由于铜箔烧蚀能量大于PI材料,通过设置CO2激光参数值,在未刻蚀掉的顶层铜箔层保护下对PI材料层进行光热烧蚀,直至露出基材层的下表面铜层,从而制造出光纤高精度矩型定位槽,槽深为H。
在挠性PI基材上制作的四条光纤矩形槽在光学显微镜下的横截面如图8所示,光纤矩形槽表面光洁,槽底也无残留基材。光纤埋入精度取决于制造的光纤矩形槽尺寸参数精度。选取某一矩形槽横截面,利用Olympus STM6光学显微镜测量了不同位置的宽度值、深度值及间隔值,如图8中所示,测量结果如表2所示。
表 2 光纤定位槽参数测量结果
Table 2. Measurement results of structural parameters of optical fiber positioning groove
Measurement item Measurement value/μm Design value/μm Error value/μm Fpw (copper layer) 155.95 155 0.95 Fpw (position 1) 155.78 155 0.78 Fpw (position 2) 140.88 155 −14.12 Pitch (P) 244.87 250 −5.13 Depth (H) 118.29 120 −1.71 从测量结果可知,在位置2处,矩形槽宽度误差值最大为−14.12 μm,由于光纤为圆形,该误差对光纤的埋入精度影响不大。位置1宽度误差值为0.78 μm,顶层铜采用传统化学刻蚀工艺,误差值仅为0.95 μm。矩形槽间距P为244.87 μm,误差值为−5.13 μm;矩形槽深度H为118.29 μm,误差值为−1.71 μm。
-
将有涂覆层的光纤布置到光纤定位槽中,然后在埋入光纤的双面覆铜PI基材层的上下表面放置半固化片和PI保护层进行叠层与对位,并按照图4所示的层压工艺曲线进行层压。半固化片在高温高压条件下将变为液态,并填充光纤与定位槽中间的空隙,包裹光纤,降低光纤在层压过程所受热应力,保证其不被损坏,半固化片厚度为50 μm。制作的FEOPCB如图9(a)所示,埋入光纤长度为10 cm和15 cm。对FEOPCB进行截面切片金相观察分析,如图9(b)所示,内埋有涂覆层光纤端面光洁,无高温降解、无裂纹及断裂等缺陷。
-
虽然半固化片在层压过程中变为液态,包裹光纤,降低其所受热应力,同时也会导致埋入光纤偏移。图10所示为埋入光纤偏移量测量示意图,假设光纤布置在矩形定位槽中心线位置,光纤与槽底层铜相切。光纤中心位置在O点,通过高温层压后的位置变为O',x和y分别为光纤横向和纵向偏移量。横向偏移量为光纤中心到矩形定位槽中心线距离值,纵向偏移量为光纤中心到底层铜距离值减去77.5 μm (光纤半径值),通过光学显微镜测量并计算FEOPCBA端口和D端口光纤偏移量,如图9(b)所示,结果如表3所示。
表 3 埋入光纤偏移量测量结果
Table 3. Measurement values of buried optical fiber offset
Path
numberSection
numberLateral
offset/μmVertical
offset/μm1 A 2.98 7.15 D 3.87 6.94 2 A −0.56 6.26 D −1.24 5.74 3 A 0.89 5.70 D 0.66 6.07 4 A −2.58 6.84 D 2.27 6.71 横向偏移量正负号分别表示埋入光纤相对于矩形槽中心线向左和向右偏移。从测量计算结果可知,埋入光纤横向偏移量最大值为3.87 μm,最小值为0.56 μm,说明制作的光纤矩形定位槽具有较高精度;纵向偏移量最大值为7.15 μm,最小值为5.70 μm,差值为1.45 μm,说明埋入光纤底部填充胶厚度较均匀。根据参考文献[18-19],埋入光纤横向和纵向偏移量小于±10 μm,均能够保证VCESL激光器和PIN探测器阵列与埋入光纤之间具有较高的耦合效率。在实际应用过程中,可对光纤端面依次经过粒度5、3、1 μm研磨砂纸和抛光垫进行研磨抛光处理,并配合粒度为25~35 nm抛光液,获得平整光洁的端面,可进一步提升耦合效率。
-
为了验证环境适应性试验和10万次弯曲疲劳试验后FEOPCB样件的性能,对其弯曲损耗、通道间串扰及高速信号传输性能展开了测试。试验测试设备包括:Agilent示波器,型号为86100D;Anritsu误码分析仪,型号为MP2100A;中电科思仪光功率计,型号为6334B;Finisar公司光模块,波长为850 nm,速率为10 Gbps,功耗约为1.5 W。
-
对FEOPCB在不同弯曲半径下的损耗开展了测试,其测试示意图如图13所示。首先采用850 nm的光模块作为光源,通过标准LC-MT转接线与FEOPCB一端相连,再通过标准MT-FC转接线与光功率计相连,形成光通路,测量FEOPCB输出功率P。将FEOPCB在不同弯曲半径下分别弯曲90°和180°,分别测量输出功率P0,通过公式(1)计算FEOPCB中埋入光纤所有模式总功率的损耗[20],每组数据测量10次,并取平均值。
$$ B L=10 \lg \left(P / P_0\right) $$ (1) 由于4个通道测量计算结果值接近,仅展示了其中某一通道的弯曲损耗值,如图14所示。从图中可知,FEOPCB最小弯曲半径可达2 mm,在90°和180°弯曲损耗约为0.57 dB和1.12 dB。弯曲损耗随着弯曲半径增加而降低,在弯曲半径为5 mm时,90°和180°弯曲损耗约为0.16 dB和0.36 dB。
-
相邻光纤之间的串扰必须足够低,才能保证FEOPCB埋入光纤中信息的稳定传输。目前,串扰的计算和测试方法主要基于模式耦合理论和功率耦合理论,其中功率耦合理论更适用于FEOPCB中相邻光纤串扰的计算和测试[21]。则相邻光纤之间的串扰计算表达式为:
$$ {XT}_{i j}=10 \lg \left(P_j / P_i\right) $$ (2) 式中:Pi为第i通道的光纤输出功率;Pj为第i通道光纤耦合至相邻第j通道光纤输出功率,其测试示意图如图15所示。
为了与VCSEL激光器和PIN探测器阵列进行对准耦合,FEOPCB埋入裸光纤纤芯间距为250 μm,通过测量,相邻裸光纤之间无芯间串扰。参考文献[22]中多芯光纤的芯间距为38 μm,其串扰值已符合传输要求。目前,国际上多芯光纤的芯间距一般为42 μm[22]。
-
图16为FEOPCB高速信号传输性能测试原理示意图,图17为搭建的测试平台实物图。将误码仪输出码型设置为PRBS31,传输速率为10 Gbps。FEOPCB 4个通道在30 min测试时间的条件下,误码数均为0,误码率均小于10−16,其中某一通道测试结果如图17所示,眼图良好。
Simulation and test of optical fiber characteristics in flexible electro-optical printed circuit board process
-
摘要: 挠性光电印制电路板(Flexible Electro-Optical Printed Circuit Board, FEOPCB)在高温层压制作过程中,埋入光纤会产生热应力,造成光纤损坏等缺陷,影响其可靠性和高速信号传输性能。为了降低FEOPCB弯曲半径并提升其可靠性,将在双面覆铜聚酰亚胺(PI)基板上设计制作高精度矩形光纤定位槽。首先建立有/无涂覆层光纤埋入挠性基板有限元仿真模型,对FEOPCB制造工艺进行模拟仿真,并对埋入光纤应力及影响因素进行分析。结果表明,有涂覆层光纤所受应力远小于无涂覆层光纤。针对有涂覆层光纤,采用激光刻蚀技术在双面覆铜PI基板上制作了高精度矩形定位槽,通过高温层压工艺完成了FEOPCB制作。FEOPCB完成了温度冲击、低温、高温、湿热和10万次弯曲疲劳可靠性试验,利用光学显微镜观察分析,埋入光纤无高温降解和破裂等缺陷。FEOPCB最小弯曲半径小至2 mm,弯曲损耗分别为0.57 dB (90°)和1.12 dB (180°),且相邻光纤之间无串扰,在850 nm波长条件下通信速率可达10 Gbps,误码率小于10−16。Abstract:
Objective The traditional electric interconnection method has become the bottleneck to limit the rapid development of high-speed communication electronic products for its inherent physical characteristics in the case of high frequency. Optical interconnection technology can be used instead of electric interconnection technology to realize high-speed, large-capacity, high-density and flexible information transmission in electronic products and eliminate the technical bottleneck. As a new development direction of board-level optoelectronic interconnection technology, FEOPCB can realize flexible interconnection among different subsystems and meet the development trend of lightweight, miniaturization and high performance of high-speed electronic systems. However, during the high-temperature lamination manufacturing process of FEOPCB, the embedded fibers will generate thermal stress, which will cause damage to the embedded fibers, affecting high-speed signal transmission performance and reliability of FEOPCB. Therefore, it is necessary to establish the finite element model with bare optical fibers embedded for simulating and analyzing the thermal force to guide the design and manufacturing of FEOPCB. For this purpose, the research work on simulation and test of fiber characteristics in FEOPCB process was carried out in this paper. Methods In order to reduce the bending radius and improve its reliability, high-precision rectangular positioning grooves for the fibers were designed and fabricated on polyimide substrate with double-sided copper-clad. Firstly, finite element simulation models of fibers with and without coating layer embedded in PI substrate were established and the manufacturing process of FEOPCB was simulated and analyzed with the influence factors of stress for embedded fibers (Fig.1). The results indicate that the stress of the coated fibers is much smaller than that of the uncoated fiber (Fig.5-7). Then, the laser-etching technology was used to fabricate the high-precision rectangular positioning grooves on the double-sided copper-clad PI substrate for the coated fibers (Fig.8). FEOPCB was fabricated by the high-temperature lamination process (Fig.9). Results and Discussions FEOPCB has completed the reliability tests of temperature shock, low temperature, high temperature, wet heat and bending fatigue for 100 000 times (Fig.11). Through the observation and analysis with optical microscopy, the result shows that the embedded fibers have no high temperature degradation and cracking defects under high temperature lamination process (Fig.12). The minimum bending radius of FEOPCB is as small as 2 mm, and the bending loss is 0.57 dB and 1.12 dB respectively under 90° and 180° bending conditions (Fig.14). The measurement results show that there is no crosstalk between adjacent fibers. Finally, the high-speed signal transmission performance was measured which indicated that a 10 Gbps communication rate with bit error rate of 10-16 could be reached at the wavelength of 850 nm (Fig.17). Conclusions In this study, the finite element analysis method is used to establish the model with coated or uncoated optical fiber embedded in rectangular groove of PI substrate, and the stress and influence factors of embedded optical fiber are analyzed. The results show that the stress of uncoated fiber decreases from 129.72 MPa to 116.80 MPa, and the stress of coated fiber decreases from 89.47 MPa to 52.02 MPa with the increase of intermediate PI layer thickness. The stress value tends to be stable, when the thickness value is greater than 140 μm. With the increase of the thickness of the filling adhesive at the bottom of the optical fiber, the stress on the uncoated optical fiber decreases from 125.04 MPa to 86.82 MPa, and then increases to 93.53 MPa, and the stress on the coated optical fiber increases from 81.30 MPa to 84.52 MPa. The transverse and longitudinal offsets of the embedded optical fibers at both ends of FEOPCB were measured, and the maximum values were 3.87 μm and 7.15 μm, respectively. It can ensure high coupling efficiency between bare optical fiber and photoelectric chip. FEOPCB has completed the reliability experiments and performance tests perfectly. The research results show that the coated optical fiber can be compatible with the printed circuit board lamination process. FEOPCB has high reliability and can meet the requirements of high-speed signal transmission. -
表 1 材料属性
Table 1. Properties of materials
Material Elastic
modulus/GPaPoisson's
ratioThermal expansion
coefficient/℃Density/
kg∙m−3Heat capacity/
J∙(Kg∙℃)−1Thermal conductivity/
W∙(m∙℃)−1Copper layer 110 0.33 18×10−6 8940 384 398 Protective layer 3.20 0.36 22×10−6 1420 1090 0.12 Silica fiber 71.9 0.16 0.55×10−6 2200 745 1.50 Filled resin 7.84 0.30 27×10−6 970 1600 0.21 Coating layer 3.20 0.36 22×10−6 1420 1090 0.12 表 2 光纤定位槽参数测量结果
Table 2. Measurement results of structural parameters of optical fiber positioning groove
Measurement item Measurement value/μm Design value/μm Error value/μm Fpw (copper layer) 155.95 155 0.95 Fpw (position 1) 155.78 155 0.78 Fpw (position 2) 140.88 155 −14.12 Pitch (P) 244.87 250 −5.13 Depth (H) 118.29 120 −1.71 表 3 埋入光纤偏移量测量结果
Table 3. Measurement values of buried optical fiber offset
Path
numberSection
numberLateral
offset/μmVertical
offset/μm1 A 2.98 7.15 D 3.87 6.94 2 A −0.56 6.26 D −1.24 5.74 3 A 0.89 5.70 D 0.66 6.07 4 A −2.58 6.84 D 2.27 6.71 -
[1] Zhou Dejian, Cheng Lei. The domestic and international research situation of photoelectric interconnection technology [J]. Advanced Materials Research, 2013, 760-762: 383-387. doi: 10.4028/www.scientific.net/AMR.760-762.383 [2] 周德俭, 成磊, 汪智奇. 光电互联技术及其光电集成器件/组件发展概述[C]//第十一届中国高端SMT学术会议论文集, 2017: 3-16. Zhou Dejian, Cheng Lei, Wang Zhiqi. The Development overview of photoelectric interconnection technology and photoelectric integrated devices/components [C]//Proceeding of 11th China Upscale SMT Academic Conference, 2017: 3-16. (in Chinese) [3] 成磊. 埋入光纤的光电互联PCBA关键技术研究[D]. 西安: 西安电子科技大学, 2017. Cheng Lei. Research on the key technologies of photoelectric interconnection PCBA embedded silica fibers [D]. Xi’an: Xidian University, 2017. (in Chinese) [4] Yang Wei, Mao Jiubing, Feng Xiaojuan. Research status and developing tread for waveguide-based board-level optical interconnects technology [J]. Laser & Optoelectronics Progress, 2016, 53(6): 060004. (in Chinese) [5] Mao Jiubing, Yang Wei, Feng Xiaojuan, et al. Research status of flexible electro-optical circuit for interconnection [J]. Laser & Optoelectronics Progress, 2016, 53(8): 080004. (in Chinese) [6] 王一鸣. 基于光波导的光电互联技术的设计与实现[D]. 西安: 西安电子科技大学, 2019. Wang Yiming. The design and implementation of board-level optoelectronic interconnects technology based on optical waveguides[D]. Xi’ an: Xidian University, 2019. (in Chinese) [7] 方云山. 基于光波导的板级光电互联技术的研究[D]. 西安: 西安电子科技大学, 2018. Fang Yunshan. Research on board-level optoelectronic interconnection technology based on optical waveguide [D]. Xi’an: Xidian University, 2018. (in Chinese) [8] Nasu H, Nagashima K, Uemura T, et al. >1.3-Tb/s VCSEL-based on-board parallel-optical transceiver module for high-density optical interconnects [J]. Journal of Lightwave Technology, 2018, 36(2): 159-167. doi: 10.1109/JLT.2017.2766081 [9] Nasu H, Nagashima K, Uemura T, et al. >1-Tb/s on-board optical engine for high-density optical interconnects [C]//IEEE Optical Fiber Communication Conference, 2017: 1-3. [10] Bamiedakis N, Shi F, Chu D, et al. Mode-mixing in multimode polymer waveguides for on-board optical interconnects [C]//Optical Interconnects XIX of SPIE, 2019, 10924: 1092404. [11] Brusberg L, Weber D, Perbtgaker D , et al. Silicon photonics packaging on board-level [C]//IEEE Photonics Conference, 2015. [12] She Yulai, Zhou Dejian, Chen Xiaoyong. Optimization design of bending performance for ultralow loss hole assisted fiber [J]. Infrared and Laser Engineering, 2019, 48(9): 0918006. (in Chinese) doi: 10.3788/IRLA201948.0918006 [13] She Yulai, Zhou Dejian, Chen Xiaoyong, et al. Study on ultralow bending loss of bend-insensitive single mode optical fiber [J]. Optical Fiber Technology, 2019, 50(1): 225-232. (in Chinese) doi: 10.1016/j.yofte.2019.03.025 [14] Cheng Lei, Zhou Dejian, Wu Zhaohua. Design of the optical fiber embedded structure in PCB-level circuit based on the fracture mechanism [J]. Journal of Mechanical Engineering, 2016, 52(9): 158-164. (in Chinese) doi: 10.3901/JME.2016.09.158 [15] She Yulai, Zhou Dejian, Chen Xiaoyong, et al. Effect analysis of embedded fiber transmission performance by lamination process [J]. Acta Photonica Sinica, 2020, 49(3): 0306001. (in Chinese) [16] Norio Chujo, Rika Nomura, Toru Yazaki, et al. 7.2-Tb/s Compact optical backplane using ribbbon fiber sheet and high-density connector [C]//IEEE Optical Interconnects Conference, 2014: 93-94. [17] 郑强. 聚酰亚胺树脂基和环氧树脂基两类覆铜板的研究[D]. 厦门: 厦门大学, 2014. Zheng Qiang. Studies on two kinds of copper clad laminates based on polyimide and epoxy resin [D]. Xiamen: Xiamen University, 2014. (in Chinese) [18] Mao Jiubing, Yang Wei, Feng Xiaojuan, et al. Simulation and analysis of indirect coupling for laser and optical fiber in electro-optical interconnection circuit [J]. Electro-Mechanical Engineering, 2019, 35(1): 50-54. (in Chinese) doi: 10.19659/j.issn.1008-5300.2019.01.012 [19] Marc S, Thomas K, Jurgen M, et al. Fibers in printed circuit boards with passively aligned coupling [J]. Journal of Lightwave Technology, 2010, 28(15): 2121-2128. doi: 10.1109/JLT.2010.2053019 [20] Dong Fanlong, Zhao Fangzhou, Ge Tingwu, et al. Fiber bending impacts on beam quality of Yb-doped fiber laser [J]. Infrared and Laser Engineering, 2014, 43(11): 3565-3569. (in Chinese) doi: 10.3969/j.issn.1007-2276.2014.11.011 [21] Li Zenghui, Li Shuguang, Li Jianshe, et al. Double-trench assisted thirteen-core five-mode fibers with low crosstalk and low non-linearity [J]. Acta Physica Sinica, 2021, 70(10): 104208. (in Chinese) doi: 10.7498/aps.70.20201825 [22] Ma Liling, Li Shuguang, Li Jianshe, et al. A kind of single trench 19-core single-mode heterogeneous fiber with low crosstalk and anti-bending performance [J]. Acta Physica Sinica, 2022, 71(10): 104206. (in Chinese)