Optimum Design of Annular Double Done Shaped Charge Structure
-
摘要: 在环形聚能装药结构中，单锥罩结构形成的射流中间有堆积现象，断裂前射流拉伸长度有限，而双锥罩射流兼顾了上锥小锥角形成高头部速度，下锥大锥角增大射流有效质量的优点，形成的射流更加细长，头部速度高且不易断裂。基于环形切割聚能装药战斗部，综合考虑上锥角大小、上锥罩占药型罩高的比例、药型罩的高度以及药型罩壁厚对射流侵彻能力的影响，并基于灰关联理论对双锥罩环形聚能装药的优化提供依据，通过数值仿真，研究表明：上下锥角对射流成型影响最大，通过比较，当上锥罩为34°、上锥占罩高比例为40%、药型罩高度为70 mm、药型罩壁厚为5 mm时，形成的射流头部速度高，且在空气中能够稳定飞行。相比单锥罩结构，双锥罩射流细长，在空气中飞行时间长，对靶板的侵深大于单锥罩射流。Abstract: There is accumulation in the middle of the jets formed by the single cone liner in the annular shaped charge structure, and the length of the jet is limited before the fracture.The jet of the double cone liner takes both the high head speed formed by upper small angle cone liner and increscent effective mass of the jet formed by the lower cone liner.Therefore, the jet is more slender with high head speed and it is not easy to break.Based on the annular cutting shaped charge warhead, the influence of the upper cone size, the proportion of the upper cone cover, the mask height and the wall thickness of the liner on the penetration capability of jet is considered.It is concluded that the upper cone angle has the greatest influence on the jet molding.By comparison, when the top cone cover is 34°, the proportion of top cone is 40%, the height of the liner is 70 mm and the wall thickness of the liner is 5 mm, the jet head can fly steadily in the air with a high speed.Compared with the single cone liner structure, the jet formed by the double cone liner is long and slender, and has long flight time in the air, its penetration capability is greater than that of the single cone jet.
-
Key words:
- tandem warhead /
- annular double cone liner /
- orthogonal design /
- numerical simulation
-
为应对日益复杂的战场环境，需要不断创新毁伤模式。横向效应增强体(Penetrator with Enhanced Lateral Effect, PELE)是一种新型侵彻弹，最早由法德圣路易斯研究所于1996年提出[1]。PELE弹的典型结构特点是采用密度不同的外壳和内芯，无需引信和装药，其中：外壳一般为高密度材料，如合金钢、钨合金；弹芯为低密度材料，如铝、尼龙、聚乙烯等。弹丸侵彻目标时，由于低密度弹芯的侵彻能力弱，弹芯被压缩产生径向膨胀，弹体内部压力迅速上升；弹丸穿过目标后，内外压力差使壳体断裂，形成破片。由此可见，PELE弹兼具穿甲和破片杀伤效果，具有结构简单、低成本、弹体使用安全等优点，可出色地完成高效毁伤[2-3]。
传统侵彻弹的毁伤机理和作用特点研究已经相当成熟[4]，而PELE弹不同于一般的侵彻弹，需要在靶后形成有效破片，其弹体结构的侵彻规律也有自己的特点，因此学者们开展了大量的研究工作。张谋等[5]通过数值仿真技术分析了PELE弹内芯与横向效应的关联；朱建生等[6]对PELE弹的破碎机理和弹体结构进行了理论分析；杜忠华等[7]通过理论分析与数值仿真，认为PELE弹撞击金属薄板时，壳体的横向速度随着弹丸着靶速度和装填材料声阻抗的增大而增大。
在现代战场上，随着目标防护性能的增强，单一功能的弹丸很难达到理想的破坏效果，这就要求进攻模式多元化。为此，本研究提出分段式PELE弹的概念，期望通过分段式PELE弹体结构增强毁伤能力，尤其是应对多层目标。然而，分段式PELE弹的穿靶过程受多种因素影响，本研究采用数值模拟方法重点探讨弹体的侵彻速度和靶板厚度对其终点效应的影响规律，以实现分段式PELE弹应对不同目标时发挥出最强的优势性能。
1. 数值方法与模型
利用非线性动力学软件LS-DYNA进行数值模拟。弹体采用两段PELE结构，用螺纹连接，其结构如图 1所示，有限元模型见图 2。第1段和第2段弹体的长度分别为32和40 mm，外径10 mm，弹芯尺寸为Ø6 mm×30 mm。靶板为4层，每层靶的材料和厚度均相同，设靠近弹丸一侧为第1层。弹壳均采用钨合金，内芯材料为铝，弹尾采用钨合金。弹体材料参数列于表 1，其中ρ0为密度，E为弹性模量，G为剪切模量，ν为泊松比。
表 1 材料参数Table 1. Material parametersMaterial ρ0/(g·cm-3) E/GPa G/GPa ν Constitutive model Tungsten alloy 17.67 354 138 0.28 Johnson-Cook Aluminum 2.70 69 0.33 Plastic_Kinematic 921 steel 7.85 210 78 0.30 Johnson-Cook 为了提高计算效率，采用1/4模型模拟，约束对称面的位移和转动，设置为轴对称状态。靶板周围边界采用固定约束，限制所有方向的运动。网格尺寸为0.1~0.5 mm，弹体经过区域网格加密，其他区域采用渐变网格。网格单元均采用SOLID164单元，为八节点六面体[8]。靶板和壳体均采用Johnson-Cook模型和Grüneisen状态方程描述，弹芯采用Plastic_Kinematic模型描述。同时添加Add_Erosion，以控制材料的失效。弹体内部采用Contact_Automatic_Surface_To_Surface接触算法，弹体与靶板之间采用Contact_Eroding_Surface_To_Surface接触算法。
2. 模拟结果与分析
2.1 分段PELE弹与一般PELE弹对比
设计一个尺寸与分段PELE弹一致的普通PELE弹进行模拟，入射速度为1.4 km/s，单层靶板厚度为8 mm。两类弹体的轴向速度(va)变化曲线如图 3所示。可见：穿过第1层靶时，两种弹的轴向速度变化基本一致；180 μs左右两类弹体穿过第2层靶，轴向速度产生分离；此后，普通PELE弹的轴向速度降幅明显大于分段PELE弹。两类弹穿过4层靶板后，弹壳未破碎长度Lr如表 2所示。结果显示，分段PELE弹主要在第1层与第2层靶板之间以及第3层与第4层靶板之间形成破片。在侵彻多层靶板过程中，由于分段PELE弹在弹芯之间设置壳体间隔保护，能够在一定程度上限制每次穿靶后弹体的破碎长度，保留其后续侵彻和产生破片的能力，因此分段PELE弹在侵彻多层靶过程中的破片分配更合理。相比之下，普通结构的PELE弹在穿过前两层靶板后，弹体基本全部破碎，难以对后续目标产生有效毁伤。对比表明：分段PELE弹相较普通PELE弹，具有更强的侵彻能力。
表 2 弹壳未破裂长度Table 2. Length of unbroken PELE shell after penetrationProjectile type Lr/mm 1st layer 2nd layer 3rd layer 4th layer Segmented PELE 39.75 22.82 5.12 0 Normal PELE 31.49 8.71 0 0 2.2 靶板厚度对分段PELE弹终点效应的影响
一般来说，PELE弹主要打击较薄的防护壳体及其后方的目标，而靶板厚度对弹丸的终点效应有重要影响[9]。设弹丸的初始侵彻速度为1.4 km/s，改变单层靶板厚度H(4、5、6和8 mm)，进行数值模拟。不同H下，弹丸穿过各层靶板后壳体径向速度峰值vr以及壳体破裂长度L如图 4所示。
图 4 不同靶板厚度条件下壳体的径向速度峰值和破裂长度Figure 4. Peak radial velocity and crushing length of shell under different thicknesses of target plate如图 4(a)所示，随着H的增加，分段PELE弹穿过各层靶板后壳体的径向速度峰值呈现增加趋势；只有当靶板厚度H为8 mm时，弹丸在侵彻第4层靶板过程中完全被侵蚀消耗，无破片产生，速度为零。在H不同的条件下，弹丸穿过多层靶板后壳体径向速度峰值基本逐层递减，但穿过第2层靶时例外，其壳体径向速度峰值最低。由图 4(b)可知，弹丸穿过第1层靶板后第1段弹体几乎完全破碎，表明该阶段产生的破片最多。故而弹体侵彻第2层靶板时，由于破片提前对靶板进行破坏，导致穿靶后弹壳的径向速度峰值相对其他3次较低。
H对弹壳破碎的影响主要体现在后3层靶板的侵彻过程中。增加H会使弹丸在穿过第2层靶板后的破碎长度L增加；穿过第3层靶板后，L多数情况大于穿过第2层靶后，且此时L先随H的增加而增加，当H=8 mm时，因前两次穿靶时L较大，导致第3次穿靶后L反而减少；穿过第4层靶后，L随着H的增加而逐渐变小。由此可见，当H处于5~6 mm区间时，弹丸穿透各层靶后L的分配比较均匀合理。
虽然增加H可使穿靶后壳体的径向速度峰值更大，但也会使弹体在前几次穿靶过程中壳体破碎过多，不利于后续侵彻和毁伤目标。总体而言，H在5~6 mm时，本弹体结构每次穿靶后的壳体径向速度峰值较高，破碎长度均匀，效果最好。
设计一个与分段PELE弹相同尺寸的钨合金杆，以相同的速度侵彻多层靶板。得到两种弹体侵彻不同厚度靶板后的弹孔直径(D)，如图 5所示。从图 5(b)中可以看到，分段PELE弹侵彻不同厚度靶板时，在第2层靶上的弹孔最大。随着H的增加，分段PELE弹在每层靶板上的开孔规律并不一致：对于第1层靶，D随着H的增加略微增大；第2层靶中，H在5~6 mm时D有极小值；第3层靶板中，靶板厚度为6 mm时D达到极大值；而第4层靶板中，D随着H的增加而减小。
图 5 两种弹体侵彻不同厚度靶板的弹孔直径Figure 5. Layer thickness effect on diameter of crater on multi-layer targets图 6中e表示分段PELE弹与钨合金杆侵彻靶板后弹孔直径的相对偏差，e＞0表示分段PELE弹的弹孔直径更大。数据显示：第1层靶板中的e值均较低，H=4 mm时达到最大值，但此时分段PELE弹的D相对钨合金杆仅有8.67%的提升；对于不同厚度的靶板，分段PELE弹侵彻第2层靶板的D值均远大于钨合金杆，即使对于e最低的6 mm厚靶板，D也增加27.68%，而H=8 mm时，D的增幅达到54.1%；第3层靶板的e值在H=6 mm时取极大值；第4层靶板的e值随着H的增加而减小。
图 6 钨合金杆和分段PELE杆在多层靶中造成的开孔大小差异Figure 6. Relative difference in diameters of craters caused by tungsten alloy rod and segmented PELE penetration综合来看，当H为5~6 mm时，分段PELE弹侵彻各层靶板时的弹孔均大于钨合金杆，即就相对开孔能力而言，本研究中的分段PELE弹适应的最佳靶板厚度依然是5~6 mm。
2.3 侵彻速度对分段PELE弹终点效应的影响
根据以上分析结果，选定靶板厚度为5 mm，改变弹丸初始速度v0(0.8、1.1、1.4和1.7 km/s)，进行数值模拟，所得vr和L随v0的变化曲线如图 7所示。
图 7 不同侵彻速度下壳体的径向速度峰值和破裂长度Figure 7. Peak radial velocity (a) and length of broken PELE shell (b) after different velocity penetration图 7(a)显示，弹丸以不同的v0侵彻靶板时，穿过第1层和第4层靶板后，壳体径向速度峰值随着v0的增加而增大；而穿过第2层靶后，壳体径向速度峰值相对较低；穿过第2层和第3层靶板后，壳体径向速度峰值在v0=1.4 km/s时出现极大值，在v0=1.7 km/s时反而降低。从图 7(b)中也可以看到，当v0达到1.7 km/s时，弹丸穿过第2层和第3层靶板后壳体破裂很少，此时弹丸的横向效应不显著。这表明在一定的靶板厚度条件下，分段PELE弹的v0过高反而不利于其在侵彻多层薄靶后产生适量破片。
当H=5 mm、v0在0.8~1.7 km/s范围内时：弹丸穿过第1层靶板后前段弹体完全破裂，L均为30 mm，体现分段PELE弹可在一定程度上控制破碎范围；之后的几层靶板侵彻过程中，随着侵彻速度的变化，L与vr的变化趋势一致，径向速度峰值越高，破碎长度越长。综合考虑后认为弹丸以1.4 km/s的初始速度侵彻最佳。
3. 结论
(1) 相对普通PELE弹，分段PELE弹侵彻多层靶板时，壳体破裂产生的破片分布更合理，侵彻能力更强。
(2) 分段PELE弹在侵彻多层金属薄靶后，各层靶的弹孔直径普遍大于相同金属杆侵彻形成的弹孔直径，特别是第2层靶板，其弹孔直径增大超过26%；分段PELE弹在侵彻不同厚度靶板时，均在第2层靶上的弹孔直径最大。
(3) 随着靶板厚度的增加，弹丸贯穿各层靶板后的壳体径向速度峰值一同增加，弹丸侵彻第2层至第4层靶时壳体随靶板厚度的增加而产生更多破碎。随着弹丸初速度的增加，弹丸在穿过第2层和第3层靶板后壳体径向速度峰值和破碎长度均在初始速度为1.4 km/s时达到极大值。一定范围内改变靶板厚度和弹丸初始速度，弹丸贯穿第2层靶后的壳体径向速度峰值最低，而穿过第1层靶板后的弹体破碎长度基本保持不变。
(4) 此分段PELE弹在单层靶板厚度接近一半弹丸口径、初始速度在1.4 km/s附近时有较好的终点效应。
-
图 2 双锥罩射流成型参数随上锥角的变化曲线
Figure 2. Variation of double cone liner jet parameters with the upper cone angle
图 3 双锥罩射流成型参数随壁厚的变化曲线
Figure 3. Variation of double cone liner jet parameters with the thickness of liner
图 4 双锥罩射流成型参数随药型罩高度的变化曲线
Figure 4. Variation of double cone liner jet parameters with the height of the liner
图 5 双锥罩射流成型参数随上锥罩占罩高比例的变化曲线
Figure 5. Variation of double cone liner jet parameters with the proportion of the upper cone cover
图 8 两种射流侵彻靶板结果图
Figure 8. Result diagrams of two kinds of jet penetrating target plate
表 1 B炸药计算参数
Table 1. Material performance parameters of composition B
ρ/(g·cm-3) D/(m·s-1) pCJ/GPa A/GPa B/GPa R1 R2 ω E0/GPa V0 1.713 7500 28.6 524.2 7.678 4.2 1.1 0.34 8.499 1.0 
下载: 导出CSV
ρ/(g·cm-3) G0/GPa σ0/GPa Tm0/K C/(m·s-1) S1 A β n a 8.96 47.7 0.12 1790 3940 1.49 63.5 36 0.45 1.5
下载: 导出CSV
表 3 30CrMnSiNi2A计算参数
Table 3. Material performance parameters of 30CrMnSiNi2A
A/MPa B/MPa n C m Tm/K T0/K S1 γ0 α 1280 420 0.30 0.03 1.00 1793 294 1.49 2.17 0.46
下载: 导出CSV
ρ/(g·cm-3) E/GPa μ σs/MPa 7.85 207 0.30 600
下载: 导出CSV
表 5 空气的计算参数
Table 5. Material performance parameters of air
ρ/(mg·cm-3) C4 C5 E0/(MJ·m-3) vrel 1.293 0.4 0.4 0.25 1.0
下载: 导出CSV
表 6 正交设计各因子水平值
Table 6. Orthogonal table
Level α/(°) b1/mm H/mm h·H-1/% 1 30 3.8 60 30 2 36 4.6 70 50 3 44 5.4 80 70
下载: 导出CSV
表 7 正交表设计与计算结果
Table 7. Orthogonal design and calculation result
No. α/(°) b1/mm H/mm h·H-1/% V/(m·s-1) L/mm 1 30 3.8 60 30 2807 25 2 30 4.6 70 50 2469 70 3 30 5.4 80 70 2488 70 4 36 3.8 70 70 2297 45 5 36 4.6 80 30 2804 70 6 36 5.4 60 50 2871 90 7 44 3.8 80 30 2834 40 8 44 4.6 60 70 2470 55 9 44 5.4 70 50 2508 70
下载: 导出CSV
表 8 无量纲化的参考序列与比较序列
Table 8. Dimensionless reference sequence and comparison sequence
No. X1′ X2′ X3′ X4′ Y1′ Y2′ 1 1 1 1 1 1 1 2 1 1.211 1.1667 1.6667 0.8796 2.8 3 1 1.421 1.3333 2.3333 0.8864 2.8 4 1.2 1 1.1667 2.3333 0.8183 1.8 5 1.2 1.211 1.3333 1 1.0003 2.8 6 1.2 1.421 1 1.6667 1.0228 3.6 7 1.467 1 1.3333 1 1.0096 1.6 8 1.467 1.211 1 2.3333 0.8799 2.2 9 1.467 1.421 1.1667 1.6667 0.8935 2.8
下载: 导出CSV
表 9 灰关联度矩阵
Table 9. Grey correlation matrix
Reference sequence Comparison sequence α b1 H h·H-1 V1 0.528 0.475 0.554 0.373 L 0.740 0.714 0.755 0.579
下载: 导出CSV
-
[1] 陈闯, 王晓鸣, 李文彬, 等.双锥罩射流侵彻钢靶侵深计算模型[J].兵工学报, 2014, 35(5):604-612. http://d.old.wanfangdata.com.cn/Periodical/bgxb201405005CHEN C, WANG X M, LI W B, et al.Penetration depth calculation model for biconical liner jet penetrating into steel target[J]. Acta Armamentarii, 2014, 35(5):604-612. http://d.old.wanfangdata.com.cn/Periodical/bgxb201405005 [2] 肖强强, 黄正祥, 祖旭东, 等.新结构药型罩射流对等效前舱的抗干扰研究[J].弹道学报, 2012, 24(1):47-51. doi: 10.3969/j.issn.1004-499X.2012.01.010XIAO Q Q, HUANG Z X, ZU X D, et al.Research of anti-interference capability by shaped charge jet into equivalent forward section for new-type line[J]. Journal of Ballistics, 2012, 24(1):47-51. doi: 10.3969/j.issn.1004-499X.2012.01.010 [3] 李磊, 马宏昊, 沈兆武.基于正交设计方法的双锥罩结构优化设计[J].爆炸与冲击, 2013, 33(6):567-573. doi: 10.3969/j.issn.1001-1455.2013.06.002LI L, MA H H, SHEN Z W.Optimal design of biconical liner structure based on orthogonal design method[J]. Explosion and Shock Waves, 2013, 33(6):567-573. doi: 10.3969/j.issn.1001-1455.2013.06.002 [4] 薛建伟, 纪春亮, 李国兰, 等.破甲战斗部药型罩结构设计与试验研究[C]//第15届全国战斗部与毁伤技术学术交流会论文集.北京, 2017: 416-419.XUE J W, JI C L, LI G L, et al.Structure design and experimental study on the liner of the armor piercing warhead[C]//Proceedings of the Fifteenth National Academic Conference on Warhead and Damage Technology.Beijing, 2017: 416-419. [5] 王远飞, 李国兰, 纪春亮, 等.双锥罩聚能战斗部技术研究[C]//第15届全国战斗部与毁伤技术学术交流会论文集.北京, 2017: 443-445.WANG Y F, LI G L, JI C L, et al.Research on the dual cone liner warhead[C]//Proceedings of the Fifteenth National Academic Conference on Warhead and Damage Technology.Beijing, 2017: 443-445. [6] 傅磊.鱼雷新型串联战斗部对舰艇典型舱段结构的毁伤效应研究[D].烟台: 海军航空工程学院, 2015: 50-53.FU L.Damage effect of new tandem warhead-torpedo on typical cabin structure of naval ships[D]. Yantai: Naval Aeronautical Engineering Institute, 2015: 50-53. [7] 尚晓江.ANSYS/LS-DYNA动力分析方法与工程实例[M].北京:中国水利水电出版社, 2006:140-143.SHANG X J.ANSYS/LS-DYNA dynamic analysis method and engineering example[M]. Beijing:China Water Conservancy and Hydropower Press, 2006:140-143. [8] 王成, 付晓磊, 宁建国.起爆方式对聚能射流性能影响的数值分析[J].北京理工大学学报, 2006, 26(5):401-404. doi: 10.3969/j.issn.1001-0645.2006.05.006WANG C, FU X L, NING J G.Numerical simulation of shaped charge jet formation under different ways of initiation[J]. Transactions of Beijing Institute of Technology, 2006, 26(5):401-404. doi: 10.3969/j.issn.1001-0645.2006.05.006 [9] 汪文革, 杨世军, 韩永要, 等.基于ANSYS/LS-DYNA的聚能射流侵彻装甲钢的有限元分析[J].兵工自动化, 2008, 27(3):39-41. doi: 10.3969/j.issn.1006-1576.2008.03.015WANG W G, YANG S J, HAN Y Y, et al.Finite element analysis of shaped charge jet penetrating into target based on ANSYS/LS-DYNA[J]. Ordnance Industry Automation, 2008, 27(3):39-41. doi: 10.3969/j.issn.1006-1576.2008.03.015 [10] TOGAMI T C, BAKER W E, FORRESTAL M J.A split Hopkinson bar technique to evaluate the performance of accelerometers[J]. Journal of Applied Mechanics, 1995, 63(2):353-356. http://adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=1996JAM....63..353T&db_key=PHY&link_type=ABSTRACT [11] FRANZEN R R, SCHNEIDEWIND P N.Observation concerning the penetration mechanics of tubular hypervelocity penetration[J]. International Journal of Impact Engineering, 2001, 11(3):289-303. http://www.sciencedirect.com/science/article/pii/0734743X9190040M [12] WASMUND T L.New model to evaluate weapon effects and platform vulnerability:AJEM[J]. WSTIAC Newsletter, 2001, 2(4):1-11. http://wstiac.alionscience.com/pdf/Vol2Num4.pdf [13] LEE W H, PAINTER J W.Material void-opening computation using particle method[J]. International Journal of Impact Engineering, 1999, 22(1):1-22. doi: 10.1016/S0734-743X(98)00041-4 [14] 蒋志刚, 曾首义, 周建平.刚性尖头弹垂直撞击金属厚靶板极限速度分析[J].固体力学学报, 2004, 25(3):360-364. doi: 10.3969/j.issn.0254-7805.2004.03.021JIANG Z G, ZENG S Y, ZHOU J P.Analysis on ballistic limit of metallic targets struck by rigid sarp-nosed projectiles[J]. Acta Mechanica Solida Sinica, 2004, 25(3):360-364. doi: 10.3969/j.issn.0254-7805.2004.03.021 [15] LEE Y W, WIERZBICKI T.Fracture prediction of thin plates under localized impulsive loading (Part Ⅱ):discing and petalling[J]. International Journal of Impact Engineering, 2005, 31(10):1277-1308. doi: 10.1016/j.ijimpeng.2004.07.011 [16] 李永胜.反舰导弹串级战斗部前级环形切割器的优化设计[D].烟台: 海军航空工程学院, 2011: 42-45.LI Y S.Optimization design of the front annular cutter of the anti ship missile cascade warhead[D]. Yantai: Naval Aeronautical Engineering Institute, 2011: 42-45. [17] 刘思峰.灰色系统理论及其应用[M].第5版.北京:科学出版社, 2010:45-48.LIU S F.Grey system theory and application[M]. 5th Ed.Beijing:Science Press, 2010:45-48. [18] 李永胜, 王伟力, 刘晓红, 等.串联战斗部前级的环形切割器优化设计[J].工程爆破, 2011, 17(4):1-4. doi: 10.3969/j.issn.1006-7051.2011.04.001LI Y S, WANG W L, LIU X H, et al.Design optimization on annular cutter in the front of tandem warhead[J]. Engineering Blasting, 2011, 17(4):1-4. doi: 10.3969/j.issn.1006-7051.2011.04.001 [19] WANG C, YUN S R, HUANG F L, et al.Investigating jet appearance in different background lights with high-speed frame photography[J]. Transactions of Beijing Institute of Technology, 2002, 11(2):113-116. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=bjlgdxxb-e200202001 [20] 王成, 恽寿榕, 黄风雷.W型聚能装药射流形成及侵彻的实验和数值仿真研究[J].兵工学报, 2003, 24(4):451-454. doi: 10.3321/j.issn:1000-1093.2003.04.005WANG C, YUN S R, HUANG F L.An experimental study and numerical simulation on annular jet formation and penetration[J]. Acta Armamentarii, 2003, 24(4):451-454. doi: 10.3321/j.issn:1000-1093.2003.04.005 -
下载:
下载:
点击查看大图
计量
- 文章访问数: 7755
- HTML全文浏览量: 2571
- PDF下载量: 24
