求解三维混凝土骨料模型分层高次元方程的多水平方法

三维混凝土随机骨料模型是由基体、骨料及界面层所组成的一种复杂弱不连续问题。高阶有限元方法是确保各骨料界面附近数值解保持较高精度的计算方法之一,但与低阶单元(如线性元)相比,高阶单元具有更高的计算复杂性,需要更多的计算机存储单元。本文针对混凝土骨料模型四面体分层高次有限元方程,提出了相应的两水平方法及多水平方法,这些方法本质性地将高次有限元方程化归为相应线性有限元方程的求解。在程序实现中,由于采用了分层基,因而不需建立判定未知数变量指标和所属几何节点类型对应关系的代数判据,构造相应网格层之间的转换算子也变得简单有效,从而显著提高了运算效率。数值例子验证了方法的有效性,为求解实际三维混凝土骨料模型提供了相应的快速计算方法。     

算子自定义小波弹性板单元构造及自适应分析

提出一类用于分析弹性板问题的算子自定义小波弹性板单元构造方法。该方法的优点在于根据工程问题的求解需要灵活构造具有解耦特性的算子自定义小波基，使得系统多尺度刚阵具有沿对角线的强稀疏性，从而实现了该算法在每个尺度上独立、快速求解，系统方程的求解效率得到较大提高。建立多分辨Lagrange有限元空间和多尺度计算理论，提出基于稳定完备法的算子自定义小波弹性板单元构造方法及解耦条件。依据两尺度相对误差估计，提出自适应算子自定义小波有限元算法。数值算例证明，算子自定义小波弹性板单元具有求解精度与计算效率高等特点。     

三维非均匀介质中弹性波传播的数值模拟

提出了一种三维非均匀介质中弹性波传播数值模拟的方法,文中称为三维格子法。该算法是二维格子法（一种二维非均匀介质中P－SV波传播的数值模拟算法）向三维非均匀介质情况的推广。在空间离散上该文方法与有限元方法类似,容许根据连续体的形状和介质分界面任意剖面网格,且自然满足自由表面边界条件。不同于常规有限差分法在各个节点上满足动力学微分方程,该算法通过满足各节点周围格子的整体平衡（积分平衡方程）来对问题进行求解,三维格子法所需的计算机内存及计算耗时与同阶精度的规则网格有限差分法相当。算例表明,该文提出的三维格子法具有较高的精度且可很好地模拟三维复杂形状地表对弹性波的反射和绕射。     

弱不连续问题高阶有限元离散系统的GAMG法

弱不连续问题(如含夹杂问题)是固体力学计算中的一类重要问题。高阶有限元方法由于其具有更好的逼近效果,是确保数值解在界面保持较高精度的计算方法之一。但与线性元相比,高阶单元需要更多的计算机存储单元,具有更高的计算复杂性。本文利用两水平算法的思想,将高阶有限元离散系统化归于线性元离散系统的求解,为弱不连续问题高阶有限元离散系统设计了一种新的基于几何与分析信息的代数多重网格(GAMG)法,并应用于圆形求解域含单夹杂问题的高阶有限元离散系统的求解。数值试验结果表明,相比于常用GAMG法,新方法的迭代次数基本不依赖于问题规模、单元阶次以及杨氏模量的间断性,CPU计算时间得到明显改善,具有更好的计算效率和鲁棒性,可大大提高弱不连续问题有限元分析的整体效率。     

计算无限弹性空间中共面Ⅲ型裂纹动应力强度因子的一种新方法

本文应用一种含有黎曼意义下不可积核的新的边界积分方程——位移积分方程,研究无限弹性空间内共面裂纹在瞬态SH波入射下动应力强度因子的计算。本文采用线性样条插值全离散伽辽金法数值求解此积分方程。通过对简单问题的计算,我们发现:应用位移积分方程求解本问题具有提高解的稳定性和节省机时的优点。本方法可看作冯康、余德浩提出的正则积分方程法从静断裂域到动断裂领域的推广。     

求解三维裂纹前缘SIF分布时间历程的通用权函数法和有限变分法

将三维热权函数法扩展为适用于表面力、体积力和温度载荷的通用权函数法(UWF).推导出以变分型积分方程表达的UWF法基本方程,从变分的角度,将求解三维热权函数法基本方程的多虚拟裂纹扩展法(MVCE)改造为可以适用于一般的变分型积分方程的一类新型数值方法--有限变分法(FVM).在FVM中可以引入无穷多种线性无关的局部变分模式,可以根据计算要求在求解域中插入任意多个计算节点,单一型裂纹问题FVM所得到的最终方程组的系数矩阵总是一个对称的窄带矩阵,而且对角元总是大数,具有良好的数值计算性能.FVM对于SIF沿裂纹前缘急剧变化的复杂情况具有较好的数值模拟能力和较高的计算精度,利用自身一致性,可以求得三维裂纹前缘SIF的高精度解.     

二维弹性结构入水冲击过程中的流固耦合效应

描述了一个研究弹性结构入水冲击过程中水弹性效应的数值方法,在弹性结构入水冲击过程中,流体域作用在结构上的水动力载荷由边界元法获得,而结构的弹性动力响应则由有限元方法求解,通过线性给离散Ｂｅｒｎｏｕｌｌｉ方程将有限元方程和边界元方程耦合到一起,从而获得了求解流场和结构动力响应的相互耦合的运动方程。在数值考虑了自由表面的非线性边界条件,通过引入射流单元以及最大射流厚度,较好地处理了冲击引起的射流问题。     

三维势流场的比例边界有限元求解方法

比例边界有限元法(SBFEM)是线性偏微分方程的一种新的数值求解方法。该方法只对计算域边界利用Galerkin方法进行数值离散,相对于有限元方法(FEM)减少了一个空间坐标的维数,而在减少的空间坐标方向利用解析方法进行求解;相对于边界元法(BEM),比例边界有限元方法不需要基本解,避免了奇异积分的计算,所以它结合了有限元和边界元方法的优点。本文建立了利用比例边界有限元法求解三维Laplace方程的数值模型并用于计算三维物体周围的水流场,将计算结果与解析解和边界元方法进行了对比,结果表明此方法可以很好地模拟水流场,且具有较高的计算精度。     

无网格MLPG法求解轴对称弹性体扭转问题

应用无网格局部彼得洛夫-伽辽金法(MLPG)研究轴对称弹性体扭转问题,给出了矩阵形式的控制方程,发展了MLPG求解轴对称体弹性扭转问题的数值计算方法。算例分析表明:此方法对求解此类问题具有良好的适应性,数值解能达到理想的计算精度。     

功能梯度材料矩形板的三维静动态响应分析

基于线弹性理论的基本方程,选用3个位移分量和3个应力分量作为状态变量,利用状态空间法建立了功能梯度矩形板的三维状态方程.考虑四边简支的边界条件,采用打靶法数值求解了材料常数沿板厚按幂率变化的弯曲问题和自由振动问题,为求解功能梯度材料三维弹性响应提供了一种方法.并且给出了功能梯度材料三维矩形板的静动态响应受组分材料分布以及板厚长比变化的影响规律.     

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

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

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

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