浅水舰船引起水底附近的扰动速度研究

基于浅水波浪理论,建立了亚临界和超临界航速条件下舰船运动引起水底附近扰动速度的理论计算模型。通过数值计算,揭示了扰动速度分量的变化规律和分布特点。建立了舰船扰动速度的测试系统,对船模运动引起的扰动速度进行了测量。横向扰动速度分量计算结果与实验结果基本吻合,纵向扰动速度分量计算结果与实验结果定性一致。理论与实验研究表明,...     

一种钢管混凝土格构柱极限承载力的快速计算方法

采用半波正弦曲线模拟杆件的变形曲线,在考虑紧箍效应钢管混凝土应力-应变关系和剪切变形影响的基础上,建立了杆件中截面的平衡方程,提出了钢管混凝土格构柱弹塑性极限承载力数值方法,并编制了相应的计算程序。该方法不仅适用于两端偏心相等的偏压构件及两端偏心不等的偏压构件的弹塑性极限承载力的计算,而且,其收敛速度很快,比分段合成数值法更能适用于工程计算。利用提出的计算方法和编制的程序对国内已有两端偏心相同的四肢钢管混凝土格构偏压长柱的试验结果进行了计算,并与现有规程及分段合成法进行了比较,结果表明:规程计算方法结果偏于保守,计算误差大,本文方法的计算结果与试验结果吻合良好。     

预处理方法在低速弹箭流场计算中的应用

基于Weiss-Smith预处理矩阵和全局截断预处理参数,采用有限体积方法对雷诺平均Navier-Stokes方程进行离散。对流项离散采用二阶线性重构和AUSM ＋-up格式,时间推进方法采用多重网格下的LU-SGS方法。结合M PI消息传递方法,建立了一套计算低速流动的并行数值方法。计算了低速椭球体的流场和气动力,压力系数和切应力系数计算结果与文献实验结果对比吻合度较好。生成了末敏弹的流场计算网格,对绕末敏弹流场进行了数值模拟。对多重网格下多进程的加速比和并行效率进行了测试,显示了程序良好的并行效率。计算的气动力结果与实验结果吻合。综合结果表明：本文的数值方法能够用于低速弹箭流场和气动力计算,为新型弹箭的设计和定型提供保证。     

采用NND方法计算三维喷管气流场

本文运用NND显式差分格式，计算了三维喷管气流场。气流场计算的基本方程为一般贴体坐标系下三维守恒型的欧拉方程。采用了时间分裂法和Steger-Warming矢通量分裂技术。在喷管内沿周向的每个由轴线和壁面构成的子午面上根据泊松方程生成贴体网格。本文运用三维程序计算了轴对称JPL喷管，同时与实验结果和前人采用轴对称二维程序所计算的结果做了对比。最后，本文还计算了三维矢量喷管，计算结果与现有的实验结果一致。通过轴对称JPL喷管和三维矢量喷管的计算考核，表明建立的算法和编写的计算程序是正确的。文中提出了采用子午面形式的贴体网格时奇性轴的处理方法。计算结果表明在喷管壁面处，马赫数与压强的计算结果与实验值吻合较好，而在喷管轴线处，只有当网格较密时，才能得出与实验结果接近的计算结果。     

钢管混凝土格构柱偏压承载能力分析的数值方法

采用泰勒级数作为分段插值函数,在考虑多个截面的平衡条件和紧箍效应钢管混凝土应力-应变关系基础上,提出了钢管混凝土格构柱弹塑性极限承载力数值方法,并编制了相应的计算程序。与现有的分段合成法相比,该方法考虑了剪力对柱变形的影响,不仅适用于两端偏心相等的偏压构件,而且适用于两端偏心不等的偏压构件的弹塑性极限承载力的计算。利用提出的计算方法和编制的程序对国内已有两端偏心相同的四肢钢管混凝土格构偏压长柱的试验结果进行了计算,并与现有规程进行了比较,结果表明:现有计算方法结果偏于保守,计算误差大,本文提出的计算结果与试验结果吻合良好。     

无约束气云爆燃压力场的计算

利用球形流体力学方程和燃烧反应模型导出了无约束气云弱点火爆燃过程的压力分布场 ,编制了求解压力场的计算程序 ,可计算不同尺寸气云爆炸时各点的爆炸压力值。进行了可燃气云爆燃实验 ,对计算结果进行了考核。与实验结果相比 ,计算结果的偏差小于2 0 %。     

复合受力下焊接栓钉SRC柱抗扭承载力计算

为研究压、弯、剪、扭复合受力下焊接栓钉型钢混凝土(SRC)柱的抗扭承载力计算方法,基于滞回性能试验结果,分析了试验过程和破坏形态,利用了空间桁架理论与叠加原理,提出了复合受力下焊接栓钉SRC柱抗扭计算模型和强度统一方程。结果表明:所提出的强度统一方程和抗扭承载力计算方法的计算结果与试验结果均值误差均在3%内,吻合较好。     

非规则区域中拉普拉斯方程的有限分析5点格式

建立了非规则区域的有限分析5点格式,增加了有限分析法对不规则边界的适应性。应用所提出的方法对水利工程中常见的有压和无压流动进行了计算,与实验和前人的计算结果相比较,本文的方法都能得到较为满意的结果。本文的计算格式也可以应用到其他非规则区域的计算中。     

DCD格式在破碎发射药床两相流内弹道计算中的应用

利用DCD(dispersion controlled dissipative scheme)格式,提出了一种研究发射装药发射安全性问题的两相流内弹道计算方法。将内弹道气固两相流动力学方程组中与压力有关的项进行变形,实现了用同一种格式对气相和固相统一处理,而无须分别对待,采用DCD格式无须数值粘性和人工滤波,提高了计算精度。实例计算了某榴弹内弹道两相流动力学,计算结果与实验结果吻合较好。把破碎发射药床视为混合装药结构,用DCD格式计算了发射药床不同破碎程度对发射安全性的影响。计算结果表现出了通常计算方法难以反映的破碎发射药床内弹道压力极为剧烈的变化过程和极高的危险膛压。 更多还原     

Euler-Lagrange耦合计算的GEL方法

研究了一种Euler-Lagrange耦合数值方法ghost-fluid Euler-Lagrange(GEL)方法,编写了GEL二维计算程序。其中Euler流场计算采用以SCB格式编制的二阶计算程序,Lagrange域计算采用DEFEL二维动力有限元程序。通过一维黎曼问题的计算结果与高精度PPM方法进行的比较,以及二维移动边界cylinder lift-off problem的计算结果与文献的对比,验证了GEL方法和本文程序的正确性。     

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

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

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


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).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).     

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).