水下爆炸冲击波壁压理论及数值计算方法改进研究

水下爆炸冲击波是舰船抗冲击评估中重要的载荷成分,也是水中结构物毁伤程度快速预报的关键和依据。通过小当量实验发现,由于传统 Taylor 平板理论公式忽略了冲击波波速的非线性变化 ,导致其在预报近距离水下爆炸冲击波壁压脉宽时出现偏差。为此,给出了比例爆距R/W1/3为0.11～5.30 m/kg1/3 (R为爆距,W为炸药质量)下的冲击波速度拟合公式,对传统Taylor理论公式进行修正。修正后,在R/W1/3=0.11 m/kg1/3下,壁压脉宽及冲量偏差大幅减小;在R/W1/3≥0.21 m/kg1/3下,两者偏差均小于12%。此外,在处理水下近场和中远场爆炸问题时,发现数值耗散会导致壁压峰值被明显削弱,于是提出了一种可行的数值策略消除计算中数值耗散导致的削弱效应,结果与修正的Taylor平板理论公式吻合良好,峰值偏差均小于9%。改进后的冲击波壁压理论公式及数值计算方法可为舰船抗爆抗冲击领域提供理论和技术支撑。     

侧向点爆炸作用下地下拱结构的动力响应

首先,利用几何关系与介质中超压、位移衰减公式得出结构表面自由场荷载和位移的分布形式,并对衰减因数、爆距等参数进行讨论;其次,利用MSSI（modified soil-structure interaction）相互作用模型对拱的振动方程进行正交求解,得到任意角度荷载作用下的结构弹性动力响应解析解。基于动力响应解析解,得到了土体的声阻抗对结构位移、速度、加速度时程响应曲线的影响。研究结果表明:侧向爆炸荷载作用下,岩土介质声阻抗越大,埋设其中的地下结构的位移、速度和加速度变化越大。为此,建议地下防护结构应修建在声阻抗小的岩土介质中。     

水下爆炸近场峰值压力的数值模拟

为了研究水下爆炸近场内的压力状态，使用SPH (smoothed particle hydrodynamics)方法，采用C-J爆轰模型，对水下爆炸过程进行了数值模拟。通过与峰值压力规律和中场内经验公式的比较，验证了程序的可靠性。对水下爆炸过程进行了波系分析，与不同维数水下爆炸的数值模拟结果进行了对比研究。结果表明:比距离R/a=6是波形变化的一个分界点，在R/a＜6的近场范围内，峰值压力的拟合分为两段更合适。另外，对数值结果的lnP_m-ln(R/a)曲线进行了分段幂次拟合，得到了与数值模拟结果非常吻合的拟合曲线。     

应用地运动数据反演地下爆炸震源函数

应用近场（震中距为几公里至几十公里）的地表地运动测试数据来研究爆炸震源特性，将近场地震波传播介质近似为均匀半空间模型，通过去除地表响应，反演了硬岩介质中的地下爆炸震源函数折合速度势RVP（reduced velocity potential），得到了描述地下爆炸震源的折合速度势RVP振幅谱及其参数稳态值、特征频率、过冲特征量，通过回归分析，得到稳态值和特征频率与爆炸当量的关系。     

球形爆炸容器内炸药爆炸形成的准静态气体压力

准静态压力是爆炸容器设计的重要参考数据。在球形爆炸容器内开展了爆炸加载实验,压力传感器采用齐平和导孔两种安装方式,均获得了准静态压力数据,且两组数据一致。理论推导了准静态压力的表达式,并通过拟合实测数据得到经验公式。研究结果表明:（1）爆炸冲击波在容器内部往返3次后,容器内气体压力进入准静态;（2）准静态压力与当量容积比近似呈正比例关系,比例系数为1.10 MPa·m3/kg TNT。     

爆炸地震波模拟研究

基于能够描述爆炸地震波的非平稳随机过程模型进行爆炸地震波的模拟研究。基于爆炸地震能量守恒以及经验关系,提出了能够考虑装药量和爆距影响的爆炸地震波功率谱密度和幅值包络线模型,并且利用实测的爆炸试验数据,对模型参数进行标定,给出单点、两点（微差）和多点（微差）爆炸地震波实用模拟方法。结果表明模拟得到的爆炸地震波时程能够很好地再现爆炸地震波的一些特性。     

爆炸冲击波作用下围岩与被覆结构的动力相互作用

爆炸作用下围岩与被覆结构的动力相互作用对于合理确定防护结构荷载、科学设计被覆结构具有重要意义。运用ANSYS/LS-DYNA非线性显式动力有限元程序和流-固耦合算法,对垂直爆炸作用下不同爆距、不同跨度的地下结构与围岩的动力相互作用进行了数值模拟,应用波动理论进行了动力相互作用力分析,讨论了相互作用动载计算公式在岩石结构中的适用性,得到了围岩与被覆结构的最大相互作用力变化规律。研究结果表明:在距拱顶1~25 m垂直爆炸作用下,14~25 m跨地下结构都发生了拱顶局部破坏,整个拱的混凝土均会产生震动裂缝;当爆距为4 m时,围岩与结构的动力耦合作用最大,可以作为确定最大荷载的依据。     

黄土和砂土岩中填实爆炸辐射弹性波的对比研究

研究地下爆炸弹性区的震动特性,关键是获得场地介质与爆炸能量耦合作用下辐射弹性波的实验参数。对于不易加工成大尺寸模型的砂土岩,为研究其填实爆炸下辐射弹性波的特征,采用0.125 g TNT微型炸药球作为爆炸源,以塑性区可置换的?1 370 mm×1 200 mm黄土样品作为提供应力波传播路径的载体,用波阻抗近似相等的重塑黄土和砂土岩样品分别作为源区介质,对比分析了两种介质中微药量填实爆炸辐射的弹性波传播特征。实验结果表明:在测试范围内,两种介质中填实爆炸激发的弹性应力波粒子速度(位移)峰值的衰减规律、波形的主频变化规律一致;砂土岩中爆炸辐射的弹性波粒子速度(位移)峰值整体高于黄土、粒子速度波形的半高宽和主频低于黄土;砂土岩中爆炸耦合的向外传播的弹性波能量比黄土大。实测结果反映,黄土和砂土岩中填实爆炸弹性波能量耦合强度的差别。     

黏弹性固体中地下爆炸辐射地震波能量的演化

地下爆炸与介质的能量耦合和介质中的波传播机制是理解地下爆炸源物理的重要基础。为研究地下爆炸辐射地震波能量的传播衰减规律,分析了黏弹性介质中地下爆炸地震波能量的组成。基于无限介质中黏弹性球面波理论,给出了速度、位移、应力、应变等物理量Laplace域的理论解。利用Laplace数值逆求解方法,建立了黏弹性介质中地下爆炸辐射地震波场的计算方法。以干黄土作为典型黏弹性材料,计算给出了地震波能量的传播特征,分析了地下爆炸辐射能量的传播衰减规律。结果表明:（1）在黏弹性介质中,某球面处流入的能量随半径增加而逐渐降低。在理想弹性介质中,某球面处流入的能量在几倍弹性半径外即可稳定到某一定值;（2）在某一固定的有限观测区域内,当观测时间足够长时,势能和耗散能均趋于某一定值,辐射动能趋于零;（3）当有限的观测区域能容纳一个完整波长的地震波时,地震波辐射动能的稳态值随波传播距离的增大而减小,总体上可以用指数函数和幂函数进行分段拟合。     

由邻近爆炸引起的岩洞振动规律

本文对一个有空腔耦合的爆炸多次试验结果进行了总结。重点叙述了地下岩洞运动参数(加速度、速度、位移)的衰减规律及沿空腔断面分布特征。给出了可供工程应用参考的衰减指数。     

Preface: Current research progress on mechanics of high speed rail

The construction and operation of high-speed rail（HSR）grid within the past two decades in China,in terms of the scale,may be possibly comparable to any national-wide construction in the history of China,even the Great Wall.By counting railways with commercial train service at the speed of200 km/h,China has the world’s longest HSR network with over 19 369.8 km of track in service today and this number     

ACTA MECHANICA SINICA 2014 ANNUAL SUBJECT INDEX

正~~     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

人造地震动合成中的位移误差分析

三角函数级数法是合成人工地震动常用算法之一,但是通过对加速度积分求取位移时,却存在与零线漂移相类似的位移漂移现象。