Influences of material properties of a projectile on hypervelocity penetration depth
-
摘要: 为研究弹体材料参数(主要指屈服强度、韧性等)对超高速侵彻混凝土靶侵彻深度的影响规律,开展了不同材料性能的93W合金柱形弹以
2300 ~3600 m/s的速度侵彻混凝土靶实验,得到了不同材料性能弹体的侵彻深度和残余弹体长度实验数据,并结合已有文献中的实验结果以及数值模拟方法,分析了材料参数对侵彻深度、残余弹体长度的影响规律。得到的结论如下:(1)如果弹体材料的韧性增强而强度不变,残余弹体特征参数并未显著改变,侵彻深度无显著变化,侵彻深度极大值对应的弹速也无显著变化;(2)如果弹体材料的强度提高而韧性不变,则弹体抵抗侵蚀的能力提升,使弹体残余长度增加,侵彻阶段的临界转变速度增加,进而使刚体侵彻深度和总侵深增加,同时使弹体侵彻深度极大值对应的侵彻速度提高。Abstract: In order to study the influence of projectile material parameters (mainly strength, toughness, etc.) on the penetration depth of hypervelocity penetrating concrete targets, experiments of 93W alloy column-shaped projectiles with different material properties penetrating concrete targets at2300 –3600 m/s were carried out on a 57/10 two-stage light gas gun. The projectile velocity was measured by a laser velocimetry system, of which the uncertainty is less than 1%. The experimental data of penetration depth and residual projectile length of different projectiles were obtained by computed tomography (CT) diagnosis technology, which can achieve a measurement accuracy of 0.1 mm. Combined with the experimental results and numerical simulation of Euler type finite element method in the literature, the influences of material parameters on the penetration depth and length of the residual projectile at different impact velocities were analyzed. Numerical simulation was carried out based on the AUTODYN software. In the simulation process, tungsten alloy was described by the Grüneisen equation of state and Steinberg constitutive model, while concrete was described by the pressure-porosity equation of state and RHT dynamic damage constitutive model. The conclusions obtained are as follows. (1) If the toughness of the projectile material increases and the strength does not change, the characteristic parameters of the residual projectile, the penetration depth, and the velocity of the corresponding maximum penetration depth do not change significantly. (2) If the strength of the projectile material increases and the toughness is constant, the ability of the projectile to resist erosion can be enhanced, the residual length of the projectile increases, and the critical transition speed increases, thereby increasing the rigid penetration depth and total penetration depth. At the same time, the velocity corresponding to the maximum value of the projectile penetration depth increases.-
Key words:
- hypervelocity /
- penetration /
- tungsten alloy projectile /
- strength /
- toughness /
- concrete /
- two-stage light gas gun
-
降低集束弹药未爆弹率是目前的关注热点之一, 在《特定常规武器公约》框架下, 对于《集束弹药议定书(草案)》提出的加装自毁装置(不同于发火装置)降低未爆弹率, 目前已经基本达成共识[1-2]。秦栋泽等[1]采用可靠性框图方法探讨了不同时机启动实现高效自毁问题, 认为自毁装置采用一道保险或在抛撒时启动自毁效率高, 而有学者对自毁装置仅一道保险或在抛撒时启动是否会带来弹药引信本身安全性问题尚有疑虑。本文中, 尝试通过理论建模和部分实验结果, 说明经过合理的设计可以保证自毁装置采用一道保险或在抛撒时启动不会降低其安全性, 如可利用起爆信息量大, 能提高其起爆可靠性。
1. 安全性建模
由于自毁装置主要目的是解决未爆弹问题, 需要避免瞎火, 同时不能由此引发安全性问题, 导致可靠起爆和弹药安全性的矛盾非常突出。自毁装置存在起爆信息识别率和起爆信息干扰度之间的矛盾。无论选择何种识别方法设计自毁装置, 总会出现两类错误。第一类是自毁预定条件不存在时, 由于干扰的存在, 自毁启动威胁弹药安全, 这个概率为干扰度Pe0; 第二类是自毁预定条件存在, 而自毁装置判断为不存在, 这个概率为误识率Pe1。第一类错误会导致引信的安全问题, 在自毁装置的错误指令下, 可能出现早炸。第二类错误会导致瞎火。在引信自毁装置设计中, 这两类问题所带来的危害不同, 由于缺乏先验概率, 自毁装置起爆信号检测不宜采用最小错误概率准则和最小平均风险准则, 在设计和评价环境中采用奈曼-皮尔逊(Neyman-Pearson)准则较合理。即, 保持第一类错误概率Pe0为固定的允许值的同时, 使第二类错误概率Pe1最小。这种准则兼顾了安全性和可靠性两大性能, 为各种类型的自毁装置安全性分析建立了统一的评价标准。根据自毁装置的设计准则, 自毁装置的安全性指标规定为必须满足环境干扰度指标α。一般自毁装置由环境识别器、保险器、状态控制器、起爆元件组成, 自毁装置可能有4种安全失效模式。第一种, 环境识别器失效, 保险器、状态控制器、起爆元件均可靠; 第二种, 保险器发生安全性失效、状态控制器工作可靠、起爆元件工作可靠; 第三种, 状态控制器失效、起爆元件工作可靠; 第四种, 起爆元件原发性失效。因此自毁装置安全失效率的计算式为:
PSDe0=PEe0RSRTRD+PSe0RTRD+PTe0RD+PDe (1) 式中:Pe0SD为自毁装置失效率; Pe0E为环境识别器安全失效率; Pe0S为保险器的安全失效率; Pe0T为状态控制器安全失效率; Pe0D为爆炸元件的原发性安全失效率; RS为保险器的可靠度; RT为状态控制器的可靠度; RD为爆炸元件的可靠度。
为了研究不同结构原理的自毁装置, 对自毁装置进行了理论抽象。自毁装置本质输出起爆信息, 起爆信息是从环境信源中提取一定量的信息并转变为信号, 自毁装置要达到规定的可靠性指标, 就有一个必须获得的最小信息量Imin, 若自毁装置获取的信息量大于最小信息量, 则性能可能趋于更优, 自毁可靠性高。
2. 起爆信息输出量模型
2.1 起爆信息输出量模型
自毁装置的实质, 是在引信出厂到战斗部作用于目标的全寿命周期T中, 选择对应的唯一的抛撒主发火失败后, 输出起爆信号起爆爆炸元件。定义自毁所必须处理的最小信息量Imin等于从引信所经历的N次操作中选择不可逆抛撒过程的熵H0, 即:
Imin=H0 (2) H0=−N∑i=1Pilog2Pi (3) 式中:N为引信所经历的操作总次数; Pi为第i次操作为抛撒过程的概率。
由于安全性比可靠性指标苛刻, 所以优先考虑安全性指标[3]。N是随机变量, 要考虑一个自毁装置的安全性总是将问题转化为多个自毁装置的安全失效率, 由此自毁装置的安全性指标α相当于自毁装置在[α-1]次操作中, 其环境识别器只能有一次将环境干扰判断为启动条件,
N=[α−1] (4) 式中:[X]表示不大于X的最大正整数。
考虑对安全性最不利条件, 认为Pi等概率分布, 即Pi=1/N,
Imin=−(log21N)N∑i=11N=log2N=log2[1α] (5) 2.2 顺序识别法模型
现有的自毁装置输出起爆信息识别方法主要有两种, 一种为顺序识别, 即利用M个特定的阈值开关获取信息, 开关按特定的顺序动作所包含的信息量。顺序识别方法包含的开关状态数为M!, 在这M!个状态中, 只有一种状态对应于自毁预定条件的存在, 假设Xi状态发生的概率为P(xi)(i=1, 2, …, M!), 则顺序识别方法所获取的信息量IM为:
IM=−M!∑i=1P(xi)log2P(xi) (6) 考虑最不利条件, 顺序识别方法的所有M!个状态等概率发生:
P(x1)=P(x2)=⋯=P(xM!) (7) M!∑i=1P(xi)=1 (8) P(x1)=1M! (9) 将式(9)代入式(6), 得:
IM=log2M! (10) 2.3 顺序时间窗模型
顺序时间窗识别方法, 即M个开关按预定顺序并在一定的时间区域内闭合, 才判断为自毁预定条件存在。M个开关所具有得状态数为(M-1)!2M-1, 其中(M-1)!为M个开关顺序闭合所拥有的状态数, 2M-1为每一种顺序闭合时, M-1个开关是否处于规定时间区内所处的状态(减1是因为有一个开关为时间基准)。假设Xi状态发生的概率为P(xi)(i=1, 2, …, (M-1)!2M-1), 则顺序时间窗识别方法所获取的信息量IM为:
IM=−(M−1)!22h−1∑i=1P(xi)log2P(xi) (11) 考虑最不利条件, 顺序时间窗识别方法的所有状态等概率发生:
P(x1)=P(x2)=⋯=P(x(M−1)!2M−1) (12) 由于
(M−1)!2M−1∑i=1P(xi)=1 (13) P(x1)=1(M−1)!21−M (14) 将式(14)代入式(11), 得:
IM=M−1+log2(M−1)! (15) 3. 实例分析
选取3个典型引信, 分别为M85子弹药引信(自毁装置二道保险)、XM1161引信(自毁装置一道保险)和M230SD引信(自毁装置抛撒启动), 进行分析。安全性与可靠性的结果见表 1, PT为靶场测试起爆率。
表 1 安全性与可靠性结果比较Table 1. The results of safety and reliability引信 环境识别器 保险器 状态控制器 爆炸元件 Pe0SD/10-6 IM PT/% M85子弹药引信 飘带 飘带 滑块 雷管 3.95 log22 < 94.72 XM1161引信 飘带 飘带 滑块、转子 雷管 3.92 log26 94.72 M230SD引信 电池 电池 处理器、电容 雷管 3.92 2+log22 99.83 M85子弹药引信自毁装置作用原理:子弹抛撒后, 在空气气动力作用下, 飘带打开, 拉力保险解除对滑块的约束, 在离心力作用下, 离心保险解除对滑块的约束, 滑块运动到位后, 自毁锤点燃延期管, 延时后点燃雷管, 子弹药自毁[4]。XM1161引信自毁装置作用原理:子弹抛撒后, 在空气气动力作用下, 飘带打开, 拉力保险解除对滑块的约束, 滑块释放转子, 转子启动后点燃延期管, 延时后点燃雷管, 子弹药自毁[5]。M230SD引信自毁装置作用原理:子弹抛撒后, 电池上电, 处理器定时, 定时时间到后对发火电容充电, 子弹药自毁[6]。
按照相关规定, 要求安全系统失效率不超过10-6, 即假设任一器部件的安全失效率不超过百万分之一, 各元件的可靠度都为0.99。依据公式(1), 可求得XM1161引信、M230SD引信安全失效率低于M85子弹药引信, M85子弹药引信安全失效率为3.95×10-6, XM1161引信、M 230SD引信安全失效率为3.92×10-6。M85子弹药引信在使用过程中安全失效率满足要求, 而限于安全性失效率在10-6这个数量级, 很难通过实验验证, 理论计算说明, 合理设计的自毁装置仅一道保险或在抛撒时启动不一定带来弹药引信本身安全性问题。因此, XM1161引信和M 230SD的引信安全性也应能满足要求。且一般电子元件的可靠度要高于机械元件, 因此实际情况下, M230SD安全失效率应该低于XM1161引信。
M85子弹药引信、XM1161引信起爆信息利用采用顺序识别法(由于延期药管点燃后无法施控, 不属于顺序时间窗方法), M230SD引信采用顺序时间窗识别法具有时间窗口(定时后仍有电容充电过程, 不同于延期管直接起爆雷管), M85子弹药引信开关为飘带和滑块, 采用顺序识别法, XM1161引信开关为飘带、滑块和转子, 采用顺序识别法, M230SD引信采用顺序时间窗识别法, 开关为电池、处理器和电容。简易计算, M230SD引信起爆信息输出量大于XM1161引信, XM1161引信起爆信息输出量大于M85子弹药引信, 与文献[1]采用可靠度框图的起爆效率高低排序结果一致, 说明在器件本身可靠性相同时, 能够采用的起爆信息量大时, 起爆可靠性高。依据真实的靶场测试数据, M230SD引信的可靠性也高于XM1161引信、M85子弹药引信, 一方面由于器件本身的可靠度高, 另一方面也说明, 由于其起爆信息利用量大, 减小了环境干扰, 引信装置起爆度高。
4. 结论
初步完成了自毁装置的抽象, 建立了自毁装置安全性理论计算模型, 在一定假设下计算结果说明, 自毁装置一道保险和自毁装置在抛撒时启动不一定降低弹药引信本身的安全性。在自毁装置的研究中引入了信息论方法, 建立了起爆信息输出量计算模型, 对有关实例进行了分析, 验证了理论模型的合理性, 同时说明, 在器件本身可靠性相同时, 若起爆信息利用量大, 引信装置起爆度高。
-
表 1 3种弹体的材料性能参数
Table 1. Material performance parameters of three kinds of projectiles
实验弹体 材质 σ0.2/MPa σb/MPa δ/% KIC/(MPa·m1/2) Ⅰ型弹 高韧性93W 740 950 26 160 Ⅱ型弹 高强度93W 1222 1252 8 70 Ⅲ型弹[8] 标准93W 731 878 8 130 表 2 Ⅰ型弹体(高韧性93W合金)超高速侵彻混凝土靶成坑数据
Table 2. Crater data of type Ⅰ projectile (high-toughness 93W) penetrating concrete targets at hypervelocities
实验编号 撞击速度/(m∙s−1) 攻角/(°) 侵彻深度/mm 弹坑直径/mm 弹体余长/mm 弹体余长误差/mm 1-1 2390 4 81.0 120.0 4.8 1.2 1-2 2740 6 86.0 112.1 4.6 1.2 1-3 2990 8 75.0 130.0 2.7 1.3 1-4 3310 0 69.9 142.8 0 0 1-5 3580 6 64.1 144.5 0 0 表 3 Ⅱ型弹体(高强度93W合金)超高速侵彻混凝土靶成坑数据
Table 3. Crater data of type Ⅱ projectile (high strength 93W) penetrating concrete targets at hypervelocities
实验编号 撞击速度/(m∙s−1) 攻角/(°) 侵彻深度/mm 弹坑直径/mm 弹体余长/mm 弹体余长误差/mm 2-1 2330 4 79.2 117.0 6.1 1.3 2-2 2680 5 84.6 120.8 5.1 1.2 2-3 2910 0 87.1 125.6 4.1 1.2 2-4 3350 0 82.4 145.3 3.4 1.2 2-5 3500 7 67.6 132.5 0 0 表 4 Ⅲ型弹体(标准93W合金)超高速侵彻混凝土靶成坑数据[8]
Table 4. Crater data of type Ⅲ projectiles (standard 93W) penetrating concrete targets at hypervelocities[8]
实验编号 撞击速度/(m∙s−1) 攻角/(°) 侵彻深度/mm 弹坑直径/mm 弹体余长/mm 弹体余长误差/mm 3-1 1820 7 67.0 3-2 1970 4 69.8 104.5 6.2 1.1 3-3 2020 5 80.6 103.3 6.7 1.2 3-4 2350 0 84.2 101.6 4.9 1.4 3-5 2390 4 82.5 105.5 5.6 0.1 3-6 2610 2 85.9 117.0 4.5 1.1 3-7 2660 0 84.0 115.9 4.2 0.1 3-8 2860 5 84.1 112.0 4.4 1.3 3-9 2900 4 76.7 105.9 3.2 1.4 3-10 3080 8 66.5 127.7 0 0 3-11 3190 0 68.0 128.0 0 0 3-12 3360 0 63.8 131.9 0 0 3-13 3360 4 61.0 144.5 0 0 3-14 3460 5 65.0 136.7 0 0 3-15 3660 7 58.3 141.4 0 0 注:实验3-1因靶体未加钢箍,破碎较严重,无法观测残余弹体 表 5 最大侵深时3种弹体毁伤参数的对比
Table 5. Comparison of the damage parameters for three types of projectiles at the maximum penetration depth
弹体 撞击速度/(m∙s−1) 侵彻深度/mm 弹坑直径/mm 弹体余长/mm Ⅰ型弹 2740 86.0 112.18 4.6 Ⅱ型弹 2910 87.1 125.6 4.1 Ⅲ型弹[8] 2610 85.9 117.0 4.5 表 6 标准93W合金材料模型参数
Table 6. Material model parameters of standard 93W tungsten alloy
ρ/(kg·m−3) G0/GPa σyd/GPa Tm0/K C/(m·s−1) S1 A 17600 160 1.5 2 766 4 040 1.23 183.85 (G′p·G0−1)/GPa−1 (G′T·G0−1)/K−1 β n γ0 a' 0.0094 0.00014 7.7 0.13 1.67 1.3 表 7 混凝土的材料模型参数
Table 7. Material parameters of concrete
G0/GPa fc/MPa ft/fc fs/fc A B ρ/(kg·m−3) M D1 D2 εminf N 16.7 42.7 0.1 0.18 1.4 1.4 2.2 0.5 0.04 1 0.01 0.5 -
[1] GOLD V M, VRADIS G C, PEARSON J C. Concrete penetration by eroding projectiles: experiments and analysis [J]. Journal of Engineering Mechanics, 1996, 122(2): 145–152. DOI: 10.1061/(ASCE)0733-9399(1996)122:2(145). [2] 王明洋, 邱艳宇, 李杰, 等. 超高速长杆弹对岩石侵彻、地冲击效应理论与实验研究 [J]. 岩石力学与工程学报, 2018, 37(3): 564–572. DOI: 10.13722/j.cnki.jrme.2017.1348.WANG M Y, QIU Y Y, LI J, et al. Theoretical and experimental study on penetration in rock and ground impact effects of long rod projectiles of hyper speed [J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(3): 564–572. DOI: 10.13722/j.cnki.jrme.2017.1348. [3] 李干, 宋春明, 邱艳宇, 等. 超高速弹对花岗岩侵彻深度逆减现象的理论与实验研究 [J]. 岩石力学与工程学报, 2018, 37(1): 60–66. DOI: 10.13722/j.cnki.jrme.2017.0584.LI G, SONG C M, QIU Y Y, et al. Theoretical and experimental studies on the phenomenon of reduction in penetration depth of hyper-velocity projectiles into granite [J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(1): 60–66. DOI: 10.13722/j.cnki.jrme.2017.0584. [4] 程怡豪, 邓国强, 李干, 等. 分层地质类材料靶体抗超高速侵彻模型实验 [J]. 爆炸与冲击, 2019, 39(7): 073301. DOI: 10.11883/bzycj-2018-0230.CHENG Y H, DENG G Q, LI G, et al. Model experiments on penetration of layered geological material targets by hypervelocity rob projectiles [J]. Explosion and Shock Waves, 2019, 39(7): 073301. DOI: 10.11883/bzycj-2018-0230. [5] 牛雯霞, 黄洁, 柯发伟, 等. 混凝土房屋结构靶的超高速撞击特性研究 [J]. 实验流体力学, 2014, 28(2): 79–84. DOI: 10.11729/syltlx2014pz38.NIU W X, HUANG J, KE F W, et al. Research on Hypervelocity impact characteristics of concrete building structures target [J]. Journal of Experiments in Fluid Mechanics, 2014, 28(2): 79–84. DOI: 10.11729/syltlx2014pz38. [6] 张浩, 张庆明. 铝弹丸超高速撞击混凝土介质冲击熔化研究 [C]//北京力学会第20届学术年会论文集. 北京: 北京力学会, 2014: 263–264.ZHANG H, ZHANG Q M. Study on hypervelocity impact melting of aluminum projectile into concrete targets [C]//Proceedings of the 20th Annual Meeting of the Beijing Society of Mechanics. Beijing: Beijing Society of Theoretical and Applied Mechanics, 2014: 263–264. [7] 钱秉文, 周刚, 李进, 等. 钨合金弹体超高速撞击混凝土靶成坑特性研究 [J]. 北京理工大学学报, 2018, 38(10): 1012–1017. DOI: 10.15918/j.tbit1001-0645.2018.10.004.QIAN B W, ZHOU G, LI J, et al. Study of the crater produced by hypervelocity tungsten alloy projectile into concrete target [J]. Transactions of Beijing Institute of Technology, 2018, 38(10): 1012–1017. DOI: 10.15918/j.tbit1001-0645.2018.10.004. [8] 钱秉文, 周刚, 李进, 等. 钨合金柱形弹超高速撞击水泥砂浆靶的侵彻深度研究 [J]. 爆炸与冲击, 2019, 39(8): 083301. DOI: 10.11883/bzycj-2019-0141.QIAN B W, ZHOU G, LI J, et al. Penetration depth of hypervelocity tungsten alloy projectile penetrating concrete target [J]. Explosion and Shock Waves, 2019, 39(8): 083301. DOI: 10.11883/bzycj-2019-0141. [9] ANTOUN T H, GLENN L A, WALTON O R, et al. Simulation of hypervelocity penetration in limestone [J]. International Journal of Impact Engineering, 2006, 33(1): 45–52. DOI: 10.1016/j.ijimpeng.2006.09.009. [10] 邓国强, 杨秀敏. 超高速武器对地打击效应数值仿真 [J]. 科技导报, 2015, 33(16): 65–71. DOI: 10.3981/j.issn.1000-7857.2015.16.010.DENG G Q, YANG X M. Numerical simulation of damage effect of hypervelocity weapon on ground target [J]. Science & Technology Review, 2015, 33(16): 65–71. DOI: 10.3981/j.issn.1000-7857.2015.16.010. [11] 章程浩, 沈培辉. 易碎穿甲弹材料性能研究 [J]. 兵器装备工程学报, 2016, 37(7): 144–148. DOI: 10.11809/scbgxb2016.07.031.ZHANG C H, SHEN P H. Study on behavior of materials used in fragile penetrator [J]. Journal of Ordnance Equipment Engineering, 2016, 37(7): 144–148. DOI: 10.11809/scbgxb2016.07.031. [12] 张德志, 唐润棣, 林俊德, 等. 新型气体驱动二级轻气炮研制 [J]. 兵工学报, 2004, 25(1): 14–18. DOI: 10.3321/j.issn:1000-1093.2004.01.004.ZHANG D Z, TANG R D, LIN J D, et al. Development of a new type gas-driven two-stage light gas gun [J]. Acta Armamentarii, 2004, 25(1): 14–18. DOI: 10.3321/j.issn:1000-1093.2004.01.004. [13] RIEDEL W, THOMA K, HIERMAIER S, et al. Penetration of reinforced concrete by BETA2B2500: numerical analysis using a new macroscopic concrete model for hydrocodes [C] // Proceedings of 9th International Symposium on Interaction of the Effects of Munitions with Structures. Berlin-Strausberg: IBMAC, 1999: 315−322. [14] 钱秉文. 钨合金弹体超高速撞击混凝土靶实验研究和机理探索 [D]. 北京: 清华大学, 2016.QIAN B W. Experiment study and mechanism exploration of hypervelocity impact of tungsten alloy projectile into concrete target [D]. Beijing: Tsinghua University, 2016. [15] Livermore Software Technology Corporation. LS-DYNA keywords user’s manual (version 971/Rev5) [M]. California: Livermore Software Technology Corporation, 2010. [16] ORPHAL D L. Phase three penetration [J]. International Journal of Impact Engineering, 1997, 20(6): 601–616. DOI: 10.1016/S0734-743X(97)87448-9. [17] ZUKAS J A. High velocity impact dynamics [M]. New York: Wiley, 1990. [18] TATE A. Long rod penetration models: part Ⅱ: extensions to the hydrodynamic theory of penetration [J]. International Journal of Mechanical Sciences, 1986, 28(9): 599–612. DOI: 10.1016/0020-7403(86)90075-5. 期刊类型引用(6)
1. 娄文忠,李诗怡,吕斯宁,范晨阳,冯恒振,朱豪坤. 基于故障模式影响及危害性分析的引信可靠性分析方法. 探测与控制学报. 2025(01): 49-56+79 . 百度学术
2. 甄国涌,陈晓锦,李祎,王淑琴. 开关脉冲式自毁指令的可靠传输设计. 电子器件. 2021(05): 1047-1052 . 百度学术
3. 程颖,张晓兵. 基于直列式引信的水雷独立自毁装置设计. 舰船电子工程. 2020(07): 178-180 . 百度学术
4. 秦栋泽,何光林. 结合MIL-STD-1316F探讨子弹药引信安全性设计中的若干问题. 弹箭与制导学报. 2019(06): 71-74+82 . 百度学术
5. 于录,宋贵宝,林木. 弹道反舰导弹靶场飞行试验安全分析. 海军航空工程学院学报. 2018(01): 163-167 . 百度学术
6. 吴奇,李晓晨,柳海斌. 基于电荷平衡原理的引信模拟自毁电路. 兵器装备工程学报. 2016(03): 102-105+120 . 百度学术
其他类型引用(2)
-