杨德庆,马 涛,张梗林
(上海交通大学船舶海洋与建筑工程学院,上海 200030)
舰艇新型宏观负泊松比效应蜂窝舷侧防护结构*
杨德庆,马 涛,张梗林
(上海交通大学船舶海洋与建筑工程学院,上海 200030)
提出一种具有宏观负泊松比效应的新型蜂窝舷侧防护结构,通过对负泊松比效应蜂窝胞元特殊结构构型设计,实现中等弹速下良好抗爆抗冲击性能。利用有限元动力学分析软件,研究鱼雷或导弹水下对舷侧防护结构的撞击侵入和穿透过程,对比研究了不同蜂窝构型、材料、胞元尺寸和胞壁厚度对舷侧结构抗冲击性能的影响。结果表明,蜂窝防护结构具有良好的抗冲击性能,负泊松比蜂窝构型较正泊松比蜂窝构型抗冲击性能更优。
固体力学;舰船防护结构;非线性有限元;蜂窝结构;负泊松比;负泊松比结构;抗冲击
为提高抗爆抗冲击能力,现代舰艇在舷侧设置空舱+液舱+空舱的多层组合防护结构。人们对舰艇舷侧防护结构技术进行了大量研究,D.V.Balandin等[1]、朱锡等[2]、杜志鹏等[3]、姚熊亮等[4]和李青等[5]研究了最佳抗冲击防护结构设计方法与抗爆机理,张延昌等[6]将蜂窝夹层板用于舰艇舷侧防护结构,增加结构塑性吸能。目前,正着力探索防护结构多尺度材料/结构一体化设计。舰艇防护结构设计的原理是,使防护结构具有高效吸收外来能量(穿甲动能、爆炸冲击能)的作用,使结构具有适当刚度来缓冲外载荷造成的变形并具有一定强度,起到主结构作用,强调结构承载性与特殊功能性的综合。工程上兼具上述效能的结构当属多孔固体材料,如泡沫和蜂窝材料等。多孔固体材料同时作为主结构承载材料和吸能材料使用,必须通过特殊的细观与宏观力学设计。特殊多尺度力学设计的多孔材料结构,可表现出宏观负泊松比和零导热率等特性[7-9]。常规蜂窝结构由于蜂窝壁厚过薄,难以承受面内大载荷并保持结构弹性状态,较少作为主结构材料使用[10-11]。将多孔固体材料微观结构构型放大形成大尺度宏观结构构型、解决抗爆抗冲击问题的研究,目前正逐步开展。本文中,提出一种具有宏观蜂窝构型防护层的新型舰艇舷侧防护结构。采用宏观蜂窝主要是解决微观蜂窝无法承受大载荷,蜂窝壁极易屈曲失效的难题。新型蜂窝舷侧防护结构提高抗爆抗冲击性能的机理为:蜂窝构型变化导致的负泊松比效应,使这类结构具有不同于普通材料结构的抗缺口压阻效应、抗断裂性能及大幅提高的回弹韧性等[7]。在穿甲过程中,穿甲破口周围材料由于负泊松比效应,会向破口聚集,将破口填充,封闭或减小弹孔,提高抗爆强度。本文中,将具有正、负泊松比宏观特性的六角形蜂窝结构应用于舰船舷侧防护结构,设计新型防护结构,模拟鱼雷或导弹在水下对舷侧防护结构的侵入和穿透过程,探讨不同蜂窝结构形式的舷侧防护结构抗爆抗冲击性能,研究其抗爆抗冲击设计方法。
新型负泊松比效应蜂窝舷侧防护结构设计如图1~2所示,也是空舱+液舱+空舱组合形式,但在舷侧板背面加装了负泊松比效应蜂窝层。舷侧舱段结构长6 m,高4 m。舱段防护结构由4层钢板构成,里面三层钢板厚均为20 mm,常规防护结构的最外层(第一层)钢板厚48 mm,各层防护板间距为0.3 m。对于新型结构第一层钢板与第二层钢板间填充蜂窝防护结构,第一层钢板厚20 mm,蜂窝胞元初始壁厚5 mm,舷侧防护结构总重量为21 330 kg。正泊松比蜂窝胞元采用等边六角形,立边长度等于斜边长度(H=L),内凹角为15°;负泊松比蜂窝胞元形状为:立边长度两倍于斜边长度(H=2L),内凹角为15°,如图3所示。后文中蜂窝胞元大小均定义为胞元斜边长度。反舰导弹为截锥形圆柱弹体,其中截顶直径为70 mm,弹体直径250 mm,弹体长度1.5 m,半锥角为20°。弹体质量为514.7 kg,弹体初始速度分别为80、200和300 m/s。弹体对舷侧结构作垂直冲击作用,高度方向距舷侧舱段结构底部为2.45 m,水平方向位于舷侧舱段结构中部,撞击部位船体无加强筋。
假设蜂窝结构材料用钢量等于常规防护结构舷侧外板厚度减少的用钢量。常规舷侧防护结构舷侧外板质量,等于负泊松比效应蜂窝舷侧防护结构舷侧外板质量,加上蜂窝层材料质量。
图1 负泊松比蜂窝结构舷侧防护几何模型及有限元模型Fig.1 Geometry and FEM model of defensive structure with re-entrant honeycomb
图2 正、负泊松比效应蜂窝夹芯舷侧防护结构(局部)示意图Fig.2 Defensive structure with honeycomb and re-entrant honeycomb (local)
图3 正、负泊松比蜂窝胞元尺寸示意图Fig.3 Size of honeycomb cell and re-entrant honeycomb cell
采用有限元动力学分析软件LS-DYNA,模拟蜂窝夹芯舷侧防护结构在受亚音速反舰导弹的撞击和侵入作用下的动态响应过程[12]。数值模拟中,考虑水与结构间的流固耦合问题。水采用空材料模型,状态方程采用Grüneisen方程,即:
(1)
式中:密度ρ=1 000 kg/m3,声速c=1 484 m/s,材料常数S1=1.979,S2=0,S3=0,γ0=0.11,α=3.0,单位体积内能E=3.072×105Pa,相对体积V=1.0。船体舷侧外板材料为45钢,蜂窝材料分别选用45钢、921钢和钛合金TC4。计算中采用Johnson-Cook本构模型,该模型是一种与应变率和绝热(忽略热传导)温度相关的塑性模型,适用于很多大应变率的材料,包括绝大多数金属材料。其中流动应力表示如下:
(2)
(3)
T*为相应温度:
(4)
破坏应变定义为:
(5)
σ*为压力与有效压力之比:
(6)
D1~D5为断裂常量,当破坏参数D达到1时即认为发生断裂:
(7)
除上述的失效准则,该材料模型还为壳单元提供了一种基于最大稳定时间步长(Δtmax)的单元删除准则。45钢、921钢和钛合金TC4材料的参数如表1所示。
表1 材料参数Table 1 Material parameters of 45 steel, TC4 and 921 steel
采用弹体剩余速度为指标,对不同蜂窝构型、蜂窝材料、胞元壁厚和蜂窝胞元层数下舷侧防护结构抗冲击性能进行了比较研究(见表2~3)。表中,h为蜂窝胞元厚度,vr为弹体剩余速度,N为胞元层数。这里剩余速度是指穿透第4层防护层的剩余弹速。200 m/s初始弹速下常规、正泊松比蜂窝和负泊松比蜂窝舷侧防护结构破损情况,见图4。
表2 弹体剩余速度Table 2 Residual velocity of missiles
表3 不同胞元层数下弹体剩余速度对比Table 3 Residual velocity of missiles with different cell layers
图4 3种舷侧防护结构破损示意图Fig.4 Crevasse shapes of three kinds of defensive structures
计算结果表明,负泊松比蜂窝构型在200 m/s中等弹速时较正泊松比蜂窝构型抗冲击性能更优,但接近音速时抗冲击性能都不佳。对于负泊松比蜂窝舷侧防护结构,蜂窝材料选用TC4,蜂窝舷侧防护胞元层数分别布置为2、3、5层,调节蜂窝壁厚,使各层数下胞元材料用量(材料总体积)相同,弹体剩余速度计算结果见表3。不同胞元层数(从左至右分别取2、3、5层)负泊松比蜂窝舷侧防护结构破损情况如图5所示。
图5 不同胞元层数下舷侧结构破损图Fig.5 Crevasse shapes of auxetic defensive structure with different layers of honeycomb cell
计算结果表明,等材料用量下蜂窝胞元层数对结构抗冲击性能的影响不是单一的,胞元布置层数并非越密越好,存在抗冲击最佳胞元层数。
针对新型宏观负泊松比效应蜂窝舰艇舷侧防护结构,探讨了蜂窝构型、材料类型、壁厚、胞元大小及蜂窝层数对舷侧结构抗冲击性能的影响。研究表明:
(1) 等材料用量情况下,新型宏观负泊松比效应蜂窝舷侧防护结构抗冲击性能优于常规舷侧防护结构;宏观负泊松比效应蜂窝舷侧防护结构较正泊松比效应蜂窝舷侧防护结构抗冲击性能更优。
(2) 蜂窝胞元材料类型能够影响整体结构抗冲击性能。
(3) 等材料用量条件下,蜂窝胞元大小对结构抗冲击性能的影响不是单一的,胞元布置并非越密越好。对于既定空间的隔舱,存在最优的胞元大小、层数和蜂窝板厚。
本文中主要针对等边六角形蜂窝构型进行了研究,有关蜂窝胞元内角大小、非等间距胞元布置等对舰船舷侧结构抗弹体冲击性能的影响,蜂窝尺寸参数优化设计等有待进行。
[1] Balandin D V, Bolotnik N N, Pilkey W D. Optimal protection from impact, shock and vibration[M]. Amsterdam: Gordon and Breach Science Publishers, 1998.
[2] 朱锡,张振华,刘润泉,等.水面舰艇舷侧防雷舱结构模型抗爆试验研究[J].爆炸与冲击,2004,24(2):134-139. Zhu Xi, Zhang Zhen-hua, Liu Rui-quan, et al. Experimental study on the explosion resistance of cabin near shipboard of surface warship subjected to underwater contact explosion[J]. Explosion and Shock Waves, 2004,24(2):134-139.
[3] 杜志鹏,李晓彬,夏利娟,等.反舰导弹攻击舰船舷侧防护结构过程数值仿真[J].哈尔滨工程大学学报,2006,27(4):484-487. Du Zhi-peng, Li Xiao-bin, Xia Li-Juan, et al. Numerical simulation of anti-ship missile attack warship broadside process[J]. Journal of Harbin Engineering University, 2006,27(6):484-487.
[4] 姚熊亮,侯明亮,李青,等.Y型舷侧结构抗冲击性能数值仿真实验研究[J].哈尔滨工程大学学报,2006,27(6):796-801. Yao Xiong-liang, Hou Ming-liang, Li Qing, et al. Numerical simulation research on counter-impingement capability Y-shape shipboard side structure[J]. Journal of Harbin Engineering University, 2006,27(6):796-801.
[5] 李青,吴广明.水面舰艇舷侧抗冲击防护结构形式初探[J].中国舰船研究,2008,3(3):26-29.
[6] 张延昌,王自力,顾金兰,等.夹层板在舰船舷侧防护结构中的应用[J].中国造船,2009,50(4):36-44. Zhang Yan-chang, Wang Zi-li, Gu Jin-lan, et al. Application of sandwich panel in anti-shock design of warship’s side structure[J]. Shipbuilding of China, 2009,50(4):36-44.
[7] Lorna J G, Michael F A. Cellular solids: Structure and properties[M]. Cambridge University Press, 2005.
[8] Klintworth J W, Stronge W J. Plane punch indentation of a ductile honeycomb[J]. International Journal of Mechanical Sciences, 1989,31(5):359-378.
[9] Kim T, Zhao C Y, Lu T J, et al. Convective heat dissipation with lattice-frame materials[J]. Mechanics of Materials, 2004,36(8):767-780.
[10] Santosa S, Wierzbicki T. Crash behavior of box columns filled with aluminum honeycomb or foam[J]. Computers and Structures, 1998,68(4):343-367.
[11] Alderson A. A triumph of lateral thought[J]. Chemistry & Industry, 1999,17:384-391.
[12] 赵海鸥.LS-DYNA 动力分析指南[M].北京:兵器工业出版社,2003.
(责任编辑 丁 峰)
A novel auxetic broadside defensive structure for naval ships
Yang De-qing, Ma Tao, Zhang Geng-lin
(SchoolofNavalArchitecture,OceanandCivilEngineering,ShanghaiJiaoTongUniversity,Shanghai200030,China)
Broadside defensive structure is important for naval vessels to maintain vitality. A novel broadside defensive structure with macro negative Poisson’s ratio is proposed to achieve higher anti-shock and anti-blast performance by design of auxetic honeycomb configuration. The process of an anti-ship missile impinging on and penetrating broadside structure is simulated by nonlinear finite element software. Effects of different design parameters on broadside structure, like auxetic honeycomb configurations, materials, sizes and thickness of honeycomb cell, are studied and compared. Numerical results indicate that counter-impingement capacity can be improved by adoption of auxetic broadside structure, and honeycomb cell with negative Poisson’s ratio is better than that of common honeycomb cell on anti-blast performance.
solid mechanics; ship defensive structure; nonlinear finite element; honeycomb structure; negative Poisson’s ratio; auxetic structure; anti-shock
10.11883/1001-1455(2015)02-0243-06
2013-07-23;
2014-01-28
国家自然科学基金项目(11072149);高等学校博士学科点专项科研基金项目(20100073110011)
杨德庆(1968— ),男,教授,yangdq@sjtu.edu.cn。
O342;U661.44;TH132.41 国标学科代码: 1301565
A