三维Laplace方程边界元中线性单元的精确积分法

边界元中的边界积分计算影响计算精度和计算速度。非奇异积分一般采用数值积分，当配置点接近积分单元时，计算精度降低。未知函数线性插值得到的解是连续解，但计算难度增大。本文采用积分区域变换，将三维Laplace问题的二维积分化为一维积分，这样奇异积分和非奇异积分能采用精确积分的方法计算，使求解精度，计算速度都得到提高。     

三维非规划非均匀边界元网格的简便的高精度算法

对三维直接边界元中一阶奇异积分，一阶近奇异积分以及非奇异积分进行统一处理，给出了一种提高积分计算精度的简便有效的方法，对非规则非均匀边界元网格可获得比一般方法高得多的计算精度，非常适合边界形状比较复杂的三维实际问题的边界元分析。     

三维非规则非均匀边界元网格的简便的高精度算法

对三维直接边界元中一阶奇异积分、一阶近奇异积分以及非奇异积分进行统一处理，给出了一种提高积分计算精度的简便有效的方法，对非规则非均匀边界元网格可获得比一般方法高得多的计算精度，非常适合边界形状比较复杂的三维实际问题的边界元分析．     

基于奇异分解和自适应BEM积分算法的水下航行体机械噪声预报

提出一种改进的声学边界元法（M-BEM）用于准确计算水下航行体发动机振动引起的近场辐射噪声。分别采用奇异分解技术和自适应边界元积分算法解决了Helmholtz积分方程在求解近场声压时出现的超奇异积分和奇异积分问题。采用一脉动球源的声辐射算例对方法进行验证,数值解与精确解误差小于1．5dB。结合有限元方法并考虑流固耦合作...     

三维常数势边界元中的精确积分

对三维问题边界元方法中应用最广泛的常数边界元的积分提出一种精确积分方法。借助于一个假想的闭合曲面,将特定的势场应用于边界积分方程,发现对于三维问题,常数势项的积分可以化作球面三角型的面积计算,而导数项的积分则可在平面域用极坐标进行。本文方法结果精确,公式简单,同一计算公式可以用来计算非奇异、几乎奇异和奇异积分,统一了积分算法。     

弹性地基板的无奇异边界元解法

本文用无奇异边界单元法,分析弹性地基上薄板的弯曲问题,考虑Winkler和双参数地基模型,选择第三类复变量的Bessel函数作为该问题的基本解,由此导出了一组权函数,在数值解法中,对面荷载积分项均统一化为边界积分,避免了在域内划分网格,文中还给出了一种考虑桩支承的边界元分析方法及域内点弯矩和接触应力的计算公式。     

非连续边界元积分的精确表达式及相关问题

以二维位势问题边界元分析为例,给出了利用线性非连续边界元离散边界积分方程时系数矩阵积分计算的精确表达式,通过和利用Gauss积分方法计算系数矩阵所得数值结果的比较表明：配位点选择不同对数值计算结果精度影响的主要原因是积分计算的精度,尤其当配位因子选择较大时,存在的准奇异积分（Nearly Singular Integrals)很难利用常规Gauss积分方法准确求得。     

梯度材料中矩形裂纹的对偶边界元方法分析

采用对偶边界元方法分析了梯度材料中的矩形裂纹. 该方法基于层状材料基本解,以非裂纹边界的位移和面力以及裂纹面的间断位移作为未知量. 位移边界积分方程的源点配置在非裂纹边界上,面力边界积分方程的源点配置在裂纹面上. 发展了边界积分方程中不同类型奇异积分的数值方法. 借助层状材料基本解,采用分层方法逼近梯度材料夹层沿厚度方向力学参数的变化. 与均匀介质中矩形裂纹的数值解对比,建议方法可以获得高精度的计算结果. 最后,分析了梯度材料中均匀张应力作用下矩形裂纹的应力强度因子,讨论了梯度材料非均匀参数、夹层厚度和裂纹与夹层之间相对位置对应力强度因子的影响.      

pFFT快速边界元方法模拟三维声散射

研究了用pFFT快速边界元方法模拟声散射问题的关键技术。采用Burton—Miller方程消除了声学边界元方法中外问题解的不唯一现象。为此,文中研究了采用常量元时该方程中超奇异积分的计算方法。最后,通过对平面声波的刚性圆球声散射的数值模拟,验证了建立的声学pFFT快速边界元方法。     

各向异性位势平面问题的规则化边界元法

基于转化域方程为边界积分方程的极限定理及一个新颖的基本解分解技术, 建立间接变量规则化边界积分方程, 它有效地避免了奇异积分的直接计算. 与已有方法比,该方法不将问题变换为各向同性的问题去处理, 因而无需反演运算, 也有别于Galerkin方法, 无需计算重积分. 可计算任意边界位势梯度, 而不仅限于法向通量. 针对椭圆边界的边值问题, 提交一种精确单元来描述边界几何. 数值算例表明, 所提算法稳定且效率高, 所得数值结果与精确解吻合较好.      

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

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

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

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