Constant strain-rate loading of liquid-driving expanding ring
-
摘要: 膨胀环实验技术主要包括爆炸膨胀环实验技术和电磁膨胀环实验技术,实验过程中膨胀环的加载应变率在达到峰值后会随着圆环的膨胀而迅速降低,给研究应变率敏感材料的拉伸碎裂带来极大的不便。在前期提出的液压膨胀环实验技术的基础上,发展了一种恒应变率加载技术。首先,从理论上获得了实现金属圆环恒应变率膨胀所需的液压加载曲线的近似表达式;然后,采用有限元流固耦合数值模拟了液压膨胀环装置中1060-O铝环的膨胀碎裂过程,在给定液压加载曲线下,膨胀环的环向应变率在应变率稳定阶段上下波动范围最大不超过20%;并进一步研究了加载曲线对碎裂过程中应变率的影响规律。在液压膨胀环实验装置上对1060-O 铝环开展了膨胀环实验,验证了恒应变率加载技术的可行性。Abstract: The expanding ring experimental technology mainly refers to the explosion expanding ring and the electromagnetic expanding ring experimental technology. During the experiment, the loading strain rate of the expansion ring decreases rapidly with the expansion of the ring after reached the peak value, which creates great inconvenience to the study of tension fragmentation of strain-rate sensitive solids. In this paper, a constant strain-rate loading technology is developed on the basis of the liquid-driving expanding ring experimental technology. Since it is not possible to apply sudden loading to the expansion ring during the experiment, it is assumed that the strain rate of the expansion ring during the expansion process is divided into linear growth stage and stable stage of the strain rate. By reasonably controlling the loading velocity and loading time of the liquid, an approximate expression of the liquid-driving loading curve required to realize the constant strain-rate expansion of the metal ring is deduced theoretically. The tension fragmentation process of the 1060-O aluminum ring under liquid-driving loading is simulated by the fluid-solid coupling numerical simulation. Under the liquid-driving loading curve, the hoop strain rate of the expanding ring fluctuates within a maximum of 20% in the stable stage of the strain rate. Before occurring of the significant necking of the expansion ring, the circumferential velocity of the expansion ring is basically zero, indicating that the expansion ring is under uniform tensile loading and there is no stress wave propagation in the circumferential direction. When the expansion ring is significantly necked, an obvious sudden change in the circumferential velocity will take place, indicateing that a Mott wave from the fracture site propagates to the corresponding position. The influence of the loading curve on the strain rate during the fracture process is further studied. Then an expanding ring experiment was carried out on the 1060-O aluminum ring on the liquid-driving expanding ring experimental device to verify the feasibility of the constant strain rate loading technology.
-
在火箭、导弹的发射过程中,固体火箭发动机会产生高温、高速的燃气射流,存在着声、光、焰、烟等明显的发射特征,不仅对周围环境及工作人员产生影响,而且容易暴露发射位置,威胁发射平台的安全性。为了克服固体火箭发动机工作过程中的固有缺陷,研究一种有效的方式来减弱发射特征具有重要意义。而通过注水方式来降低发射特征已经在欧洲、美国的大型发射场得到应用,M.M.Molnar[1]就喷水对固体火箭发动机排气羽流流场参数的影响进行了数值计算。P.Giordan等[2]采用FLUENT5对火箭发射时燃气与水的相互作用进行了数值模拟。周帆等[3]采用FLUENT软件中耦合了水的汽化方程的Mixture模型对超声速高温燃气射流的注水降温机理进行了数值模拟和实验研究。马艳丽等[4]采用FLUENT中的Mixture多相流模型求解气液两相流流场,使用离散坐标法求解辐射传输方程,就喷水对羽流红外辐射的降低效果进行了研究。张磊等[5]以单兵火箭为背景,提出了液体水柱放置在尾管中的设计,研究了注水对火箭燃气射流噪声的抑制作用。为了在此基础上更精确地计算单兵固体火箭发动机燃气射流与液体水柱之间的相互作用,本文中采用FLUENT软件中耦合了水的汽化方程的VOF模型对超声速燃气射流驱动液柱的流场进行了数值模拟。
为了对该模型的有效性进行验证,本文中将计算结果与实验现象进行对比,同时将燃气射流驱动液柱的流场计算结果与无液柱平衡体的流场计算结果进行对比分析,对液柱平衡体的降温效果进行初步分析和研究。
1. 数学模型
固体火箭发动机燃气射流驱动液柱平衡体涉及复杂的多相流问题,为了实现有效地计算,采用简化的物理模型,并提出以下假设:(1)将实际三维流动问题简化为二维湍流流动,视为非稳态膨胀过程;(2)假设火药燃烧产物为单一组分的CO2气体,采用冻结等熵流动模型,即气相燃烧产物的成分不随温度和压强的变化而变化,始终保持与喷管入口条件相一致的气体成分不变;(3)由于燃气自喷管射出在尾管内高速流动,故不考虑燃气与尾管壁面之间的热量交换,即尾管壁面视为绝热壁面;(4)视尾管壁面为无滑移壁面;(5)采用带旋流修正的k-ε模型(Realizable模型)来描述燃气射流与液体工质相互作用的过程。
1.1 控制方程
根据非稳态气液两相流场的物理过程描述,本文气液两相流动的控制方程采用分相流动模型(Separated Flow Model),在此基础上结合改进的带旋流修正k-ε模型(Realizable模型)两方程湍流控制方程,建立非稳态气液两相流场的数学模型。
(1)连续性方程:
∂∂t=(ρm)+∇⋅(ρmvm)=0 (1) 式中:vm为平均流速,vm=∑kαkρkvk/ρm,k=1, 2,分别代表气相和液相,αk为第k相的体积分数;ρm为混合密度,ρm=∑kαkρk,ρk为第k相的密度。
(2)动量守恒方程:
∂∂t(ρmvm)+∇⋅(ρmvm×vm)=−∇p+∇×(μ∇×vm)+ρmg+F+∇(∑kαmρmv2dr,k) (2) 式中:F为体积力;μm为混合粘度;μm=∑kαkμk;vdr, k为气、液相的漂移速度,vdr, k=vk-vm;g为重力加速度。
(3)能量守恒方程:
∂∂t∑k(αkρkEk)+∇⋅∑k[αkvk(ρkEk+ρ)]=∇⋅(keff∇T)+Sk (3) 式中:T为热力学温度;keff有效热传导率,keff=k+kt,kt为湍流热传导率;对可压缩相,Ek=hk-p/ρk+vk2/2,对不可压缩相,Ek=hk,hk为第k相的显焓;Sk为粘性耗散项,包含流体内热源及粘性耗散热。
(4) 气体状态方程:
p=ρRT (4) 式中:p为压强;ρ为燃气密度,ρ=1.977×10-3 g/m3;R为理想气体常数,R=8.134 J/(mol·K)。
(5) 辅助项体积含率方程:
液体工质为辅助项,体积含率方程为
∂∂t(α2ρ2)+∇⋅(α2ρ2vm)=−∇⋅(α2ρ2vdr,2) (5) 气体体积含率为α1=1-α2。
2. 计算模型
2.1 计算方法
本文中采用有限体积法来离散控制方程,湍流模型选用Realizable k-ε模型,壁面附近采用标准壁面函数。气液两相流模型采用耦合了液态水汽化方程的VOF多相流计算模型。利用PRESTO方法对压力项离散,扩散项和对流项的离散采用二阶迎风格式,体积分数离散采用QUICK格式,压力与速度耦合采用PISO算法。
2.2 计算几何模型及边界条件
单兵固体火箭发动机燃气射流驱动液柱平衡体的实验中在喷管的喉部处会放置一压力膜片,当燃烧室内的压力达到破膜压力后,燃气从喷管的扩张段喷出,所以计算中仅选取喷管的扩张段结构。为了消除重力对射流扩展形态的影响,将实验装置竖直放置,喷管向上喷射。本文中采用与实验结构1:1比例的计算模型,包含喷管扩张段、尾管以及外部大气环境空间,计算区域具体尺寸如图 1所示,喷管扩张段结构尺寸如图 2中所示,图中所示尺寸单位均为毫米。
本文中模拟计算了2种情况下的燃气射流流场:(1)无液柱平衡体——入口采用压力入口边界条件,根据实验过程中所测燃烧室内最大压力10 MPa,设置入口初始压力为10 MPa,总温为3 000 K,出口为外部大气环境,所以将出口选取为压力出口,压力设置为环境压力101 325 Pa,温度设置为300 K。壁面采用绝热壁面条件;(2)有液柱平衡体——入口、出口及壁面初始条件与无液柱平衡体时的初始条件相同,在尾管内选取长度为150 mm,宽为93 mm的长方形区域,将其设置为液体区域,液柱平衡体距喷管扩张段出口50 mm,距尾管末端20 mm,如图 1中所示。
3. 实验与计算结果分析
3.1 实验结果对比分析
图 3所示为实验所采用的发动机喷管结构,数值计算中采用的喷管扩张段尺寸与实验结构相同。为了分析燃气射流驱动液柱的降温、消烟、消焰效果,选取了实验中部分比较有代表性的照片,如图 4~5所示。图 4和5分别代表单兵固体火箭发动机在有无液柱作为平衡体2种状态下的射流流场分布。可以看出,在无液柱的情况下,燃气射流在轴向的扩展速度明显比径向的扩展速度更快,由于燃气的温度较高,所以射流流场的内部呈现出高亮特性。在t=35.00 ms时,在尾管出口出现大量白烟,这是由推进剂燃烧产物中的凝相微粒形成的。当尾管中放置液体水柱平衡体时,燃气射流会对液体水柱产生冲击,燃气射流与液体柱在强烈掺混的过程中共同从尾管中喷出。由于液体水柱的阻碍作用,气液混合射流在t=34.00 ms时才开始从尾管中喷出,射流流场的轴向发展速度也变得缓慢。如图 5所示,在t=36.00 ms时,射流流场外轮廓逐渐开始呈现锯齿状,在气液交界面处出现强烈的Richtmyer-Meshkov不稳定性,这种流动不稳定使液体发生破碎,增加了液体汽化表面,有助于液体水汽化吸收更多的燃气能量,降低了射流流场的温度,与无液体水柱的射流流场相比高亮区域不再可见。由于液体水柱的雾化,整个流场区域被大量水雾所覆盖,使得微粒对光的散射和吸收作用不足以被察觉,在无液体水柱状态下尾管出口出现的白烟也不再出现。
通过对有无液体水柱状态下的射流流场对比可以发现,在尾管后端放置液体水柱平衡体,在一定程度上改变了射流流场的机构,对射流流场有明显的降温、消焰、消烟效果,验证了该方案的可行性。
3.2 计算结果与分析
3.2.1 气液两相射流流场结果分析
为了更好的反映燃气射流与液体水柱掺混共同喷射到大气环境中的过程,对此进行了数值模拟计算。图 6所示为非稳态气液两相流场在不同时刻的液相体积分数分布云图。从图 6中可以看出,在射流流场发展的初始阶段,燃气射流驱动液柱向前发展,在尾管内形成空腔,液体在燃气推动下沿着两侧壁面运动,汽化作用即在空腔的气液交界面上进行。在t=3 ms时,燃气射流已推动液柱从尾管中喷出,尾管内有部分液体出现回流现象,粘附在两侧壁面上。随着射流流场的发展,燃气驱动大部分液体从尾管中喷出,在这个过程中气液交界面呈现出了不规则的轮廓形状,表现出了Richtmyer-Meshkov不稳定。这是由于射流流场在激波的作用下,气液界面两侧的气体与液体之间存在较大的压力梯度和密度梯度,导致了Richtmyer-Meshkov不稳定,使得两相流体间界面不光滑,具有随机脉动性,存在较强的湍流掺混,有助于液相吸收燃气中的能量蒸发雾化。在混合介质内部,燃气是被卷吸的连续相介质,而水滴是被燃气动力“撕破”的离散相介质。由于湍流脉动、燃气动力与水的表面张力复合作用,被燃气驱动的高速水射流内出现波纹式扰动,并由此破碎成水滴。这些水滴是不稳定的,在燃气动力作用下会变形、二次破碎,进一步雾化成很小的雾滴。当进行到9 ms时,射流流场中已无液相,全为气相和介于液态与气态之间的高温水蒸气。根据尾喷管的结构,认为在尾喷管内液体首次破碎,这个过程是液体整个雾化过程的开始。随着气液混合射流逐渐喷射到大气环境中,液滴发生二次破碎,与高温燃气相互作用汽化吸收热量。
图 7所示为实验与数值模拟2种条件下得到的射流流场沿轴向发展速度曲线。从图 7可以看出数值模拟所得速度发展趋势与实验所得结果吻合较好。在射流流场发展初始阶段沿轴向的发展速度较慢,随着尾管内燃气不断的推动液体与其共同喷出尾管,射流流场轴向发展速度逐渐增加,在t=6 ms时速度达到最大值,实验所得最大速度为265 m/s,数值计算所得最大速度为220 m/s,这时液体已全部从尾管中喷出。由于在数值计算时将燃气假设为单一组分的CO2气体,但实际燃气中包含多种组分,燃气在流动的过程中各组分间会进行一定的化学反应,释放出一定的能量,加快了射流流场的发展速度,这可能是导致数值计算所得射流流场沿轴向发展速度略小于实验所得发展速度的原因。随着射流流场在大气环境中发展,燃气对液体的推动作用逐渐减弱,同时液体与燃气相互作用汽化吸收了大部分能量,降低了气液混合物的动能,导致射流流场沿轴向的发展速度在达到最大速度之后逐渐减小。
3.2.2 气液射流流场与燃气射流流场对比分析
燃气驱动液柱掺混喷射过程是一个非稳态过程,射流流场中各项参数随时间和空间的变化而改变,图 8所示为有、无液柱平衡体2种状态下以喷管喉部为起点的完全发展流场沿轴线上压力(p)、温度(T)、速度(v)分布曲线。从图 8可以看出,激波的存在造成波前、波后流动参数突跃,并随着位置外移,射流流场的压强和温度在下游逐渐趋于周围大气环境的压强和温度。有液体水柱为平衡体的流场参数波动幅度小于无液柱平衡体的流场参数波动幅度,液柱减弱了流场的湍流脉动强度。从图中可以看出2种状态下在起始段射流流场均有明显的加速过程,燃气热能大部分转化为燃气的动能,所以射流流场的压力和温度在起始段均有所降低。经过喷管扩张段后,射流流场速度有一个跌落过程,而压力和温度有一个回升过程,有液柱平衡体的波动位置要比无液柱平衡体的波动位置靠前。两种状态下射流流场完全发展时,有液柱为平衡体的射流流场轴向发展位移小于无液柱平衡体的射流流场轴向发展位移,这是由于液柱对燃气的喷射过程有阻碍作用,液柱在汽化吸收燃气能量的同时降低了流场的特征参数。通过对比2种状态下射流流场各参数所达到的峰值可知,在有液柱情况下压力峰值减小了0.9 MPa,温度峰值减小了503 K,速度峰值减小了291 m/s。
图 9所示为在有、无液体水柱平衡体两种条件下射流流场发展相同位移时温度场的等值线图,从图 9可以看出,2种射流流场的温度分布明显不同,无液柱平衡体条件下温度场与周围环境发生的扰动更剧烈,射流流场头部的温度更高。在尾管出口之后的射流流场轴向核心区域内,有液柱为平衡体的射流流场温度有所降低,相对于无液柱平衡体的射流温度场降低了8.8%~23.5%,验证了实验中液体水柱对单兵固体火箭发动机工作过程中射流流场的降温效果。
4. 结论
本文中对单兵固体火箭发动机燃气射流驱动液柱的过程进行了实验研究并利用FLUENT软件对该过程进行了数值模拟,得到如下结论:
(1) 利用Realizable k-ε模型以及耦合了液态水汽化方程的VOF多相流计算模型可以对燃气驱动液柱过程的气液两相射流进行较好的描述,数值模拟结果与实验结果吻合较好;
(2) 通过实验研究,有液柱平衡体的射流流场沿轴向发展速度减慢,气液掺混喷射过程中在气液交界面上产生了Richtmyer-Meshkov不稳定,液体水柱的汽化作用减弱了流场的特征参数,起到了消烟、消焰、降温的作用;
(3) 通过数值模拟发现,有液柱平衡体的完全发展射流流场中的压力、温度、速度波动幅度要小于无液柱平衡体的完全发展射流流场波动幅度,减弱了射流流场的湍流脉动强度,射流流场中的压力峰值降低了0.9 MPa,温度峰值降低了503 K,速度峰值降低了291 m/s;
(4) 当有液柱为平衡体的射流流场与无液柱平衡体的射流流场发展到相同位置时,射流核心区域内的温度降低了8.8%~23.5%,验证了实验中液体水柱对射流流场的降温效果。
-
ρ/(kg·m−3) c/(J·kg−1·K−1) β θt/K θm/K 弹性参数 E/GPa µ 2 770 900 0.9 298 1 048 70 0.34 塑性参数 损伤演化参数 A/MPa B/MPa C n m ˙ε0/s−1 d1 d2 d3 d4 d5 Gc/(J·m−2) 27 43 0.025 0.34 1 1 0.13 0.13 −1.5 0.01 1 5 700 表 2 加载曲线中的基本物理参数
Table 2. Physical parameters in the loading curve
r0/mm h/mm R/mm ˙ε1/s−1 t1/μs 17.5 1.5 15 4 000~10 000 40 -
[1] JOHNSON P C, STEIN B A, DAVH R S. Measurement of dynamic plastic flow properties under uniform stress [C] // Symposium on the Dynamic Behavior of Materials. Albuquerque, USA: American Society for Testing and Materials, 1963: 195–198. DOI: 10.1520/STP42030S. [2] NIORDSON F I. A unit for testing materials at high strain rates [J]. Experimental Mechanics, 1965, 5(1): 29–32. DOI: 10.1007/BF02320901. [3] WARNES R H, DUFFEY T A, KARPP R R, et al. Improved technique for determining dynamic material properties using the expanding ring [C] // International Conference on the Meteallurgical Effects of High Strain Rate Deformation and Fabrication. Albuquerque, NM, USA, 1980. [4] GRADY D E, BENSON D A. Fragmentation of metal rings by electromagnetic loading [J]. Experimental Mechanics, 1983, 23(4): 393–400. DOI: 10.1007/BF02330054. [5] GOURDIN W H. Analysis and assessment of electromagnetic ring expansion as a high-strain-rate test [J]. Journal Applied Physics, 1989, 65: 411–422. DOI: 10.1063/1.343121. [6] GOURDIN W H, WEINLAND S L, BOLING R M. Development of the electromagnetically launched expanding ring as a high-strain-rate test technique [J]. Review of Scientific Instruments, 1989, 60(3): 427–432. DOI: 10.1063/1.1140395. [7] 桂毓林, 孙承纬, 李强, 等. 实现金属环动态拉伸的电磁加载技术研究 [J]. 爆炸与冲击, 2006, 26(6): 481–485. DOI: 10.11883/1001-1455(2006)06-0481-05.GUI Y L, SUN C W, LI Q, et al. Experimental studies on dynamic tension of metal ring by electromagnetic loading [J]. Explosion and Shock Waves, 2006, 26(6): 481–485. DOI: 10.11883/1001-1455(2006)06-0481-05. [8] 桂毓林. 电磁加载下金属膨胀环的动态断裂与碎裂研究 [D]. 四川绵阳: 中国工程物理研究院, 2007.GUI Y L. The studies on the dynamic fracture and fragmentation of metal freely expanding ring driven by electromagnetically loading [D]. Mianyang, Sichuan, China: China Academy of Engineering Physics, 2007. [9] 汤铁钢, 李庆忠, 陈永涛, 等. 实现材料高应变率拉伸加载的爆炸膨胀环技术 [J]. 爆炸与冲击, 2009, 29(5): 546–549. DOI: 10.11883/1001-1455(2009)05-0546-04.TANG T G, LI Q Z, CHEN Y T, et al. An improved technique for dynamic tension of metal ring by explosive loading [J]. Explosion and Shock Waves, 2009, 29(5): 546–549. DOI: 10.11883/1001-1455(2009)05-0546-04. [10] 汤铁钢, 桂毓林, 李庆忠, 等. 爆炸膨胀环实验数据处理方法讨论 [J]. 爆炸与冲击, 2010, 30(5): 505–510. DOI: 10.11883/1001-1455(2010)05-0505-06.TANG T G, GUI Y L, LI Q Z, et al. A discussion of data processing techniques for expanding ring tests [J]. Explosion and Shock Waves, 2010, 30(5): 505–510. DOI: 10.11883/1001-1455(2010)05-0505-06. [11] MOTT N F. Fragmentation of shell cases [J]. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 1947, 189(1018): 300–308. DOI: 10.1098/rspa.1947.0042. [12] GRADY D E. Fragmentation of rings and shells: the legacy of N. F. Mott [M]. Berlin, Germany: Springer, 2006. [13] KIPP M E, GRADY D E. Dynamic fracture growth and interaction in one dimension [J]. Journal of the Mechanics and Physics of Solids, 1985, 33(4): 399–415. DOI: 10.1016/0022-5096(85)90036-5. [14] 郑宇轩, 周风华, 胡时胜. 一种基于SHPB的冲击膨胀环实验技术 [J]. 爆炸与冲击, 2014, 34(4): 483–488. DOI: 10.11883/1001-1455(2014)04-0483-06.ZHENG Y X, ZHOU F H, HU S S. An SHPB-based experimental technique for dynamic fragmentations of expanding rings [J]. Explosion and Shock Waves, 2014, 34(4): 483–488. DOI: 10.11883/1001-1455(2014)04-0483-06. [15] 张佳, 郑宇轩, 周风华. 立式液压膨胀环实验技术研究 [J]. 宁波大学学报(理工版), 2017, 30(2): 35–38. DOI: 10.3969/j.issn.1001-5132.2017.02.007.ZHANG J, ZHENG Y X, ZHOU F H. Experimental technique for fragmentation of liquid-driven expanding ring [J]. Journal of Ningbo University (Natural Science and Engineering), 2017, 30(2): 35–38. DOI: 10.3969/j.issn.1001-5132.2017.02.007. [16] 李天密, 张佳, 方继松, 等. PMMA膨胀环动态拉伸碎裂实验研究 [J]. 力学学报, 2018, 50(4): 820–827. DOI: 10.6052/0459-1879-18-016.LI T M, ZHANG J, FANG J S, et al. Experimental study of the high velocity expansion and fragmentation of PMMA rings [J]. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(4): 820–827. DOI: 10.6052/0459-1879-18-016. [17] 汤佳妮, 徐便, 郑宇轩, 等. 脆性膨胀环动态拉伸碎裂实验研究 [J]. 爆炸与冲击, 2021, 41(1): 014101. DOI: 10.11883/bzycj-2020-0049.TANG J N, XU B, ZHENG Y X, et al. Experimental study for dynamic fragmentation of brittle expansion rings [J]. Explosion and Shock Waves, 2021, 41(1): 014101. DOI: 10.11883/bzycj-2020-0049. [18] 卢思凡, 张佳, 王珠, 等. 液压膨胀环动态拉伸碎裂的有限元模拟 [J]. 固体力学学报, 2019(4): 372–380. DOI: 10.19636/j.cnki.cjsm42-1250/o3.2019.012.LU S F, ZHANG J, WANG Z, et al. FEM Simulation of dynamic fragmentation of liquid-driving expanding ring [J]. Chinese Journal of Solid Mechanics, 2019(4): 372–380. DOI: 10.19636/j.cnki.cjsm42-1250/o3.2019.012. [19] 张佳. 基于SHPB的液压膨胀环实验研究 [D]. 浙江宁波: 宁波大学, 2017.ZHANG J. Research on experiment of hydraulic expanding ring technology based on SHPB [D]. Ningbo, Zhejiang, China: Ningbo University, 2017. 期刊类型引用(4)
1. 朱家萱,丁宁,蔡振涛,郭保全,潘玉田,李魁武. 燃气射流驱动液柱和堵片相互作用特性分析. 火炮发射与控制学报. 2022(03): 1-7 . 百度学术
2. 丁宁,郭保全,黄通,栾成龙,李鑫波. 筒式武器有限空间内发射流场特性研究. 弹箭与制导学报. 2020(06): 30-33+40 . 百度学术
3. 庞春桥,陶钢,闻鹏,李智宇,任保祥,王小峰. 某单兵无后坐力炮后喷区域危险性综合评估. 兵工学报. 2019(07): 1546-1552 . 百度学术
4. 郭保全,黄通,丁宁,栾成龙,张彤. 燃气射流在变截面充液管道内作用特性研究. 火炮发射与控制学报. 2019(03): 24-28 . 百度学术
其他类型引用(4)
-