近场动力学与有限元方法耦合求解热传导问题

求解含裂纹等不连续问题一直是计算力学的重点研究课题之一,以偏微分方程为基础的连续介质力学方法处理不连续问题时面临很大的困难. 近场动力学方法是一种基于积分方程的非局部理论,在处理不连续问题时有很大的优越性. 本文提出了求解含裂纹热传导问题的一种新的近场动力学与有限元法的耦合方法. 结合近场动力学方法处理不连续问题的优势以及有限元方法计算效率高的优势,将求解区域划分为两个区域,近场动力学区域和有限元区域. 包含裂纹的区域采用近场动力学方法建模,其他区域采用有限元方法建模. 本文提出的耦合方案实施简单方便,近场动力学区域与有限元区域之间不需要设置重叠区域. 耦合方法通过近场动力学粒子与其域内所有粒子(包括近场动力学粒子和有限元节点)以非局部方式连接,有限元节点与其周围的所有粒子以有限元方式相互作用. 将有限元热传导矩阵和近场动力学粒子相互作用矩阵写入同一整体热传导矩阵中,并采用Guyan缩聚法进一步减小计算量. 分别采用连续介质力学方法和近场动力学方法对一维以及二维温度场算例进行模拟,结果表明,本文的耦合方法具有较高的计算精度和计算效率. 该耦合方案可以进一步拓展到热力耦合条件下含裂纹材料和结构的裂纹扩展问题.      

近场动力学与有限元方法耦合求解热传导问题

求解含裂纹等不连续问题一直是计算力学的重点研究课题之一,以偏微分方程为基础的连续介质力学方法处理不连续问题时面临很大的困难.近场动力学方法是一种基于积分方程的非局部理论,在处理不连续问题时有很大的优越性.本文提出了求解含裂纹热传导问题的一种新的近场动力学与有限元法的耦合方法.结合近场动力学方法处理不连续问题的优势以及有限元方法计算效率高的优势,将求解区域划分为两个区域,近场动力学区域和有限元区域.包含裂纹的区域采用近场动力学方法建模,其他区域采用有限元方法建模.本文提出的耦合方案实施简单方便,近场动力学区域与有限元区域之间不需要设置重叠区域.耦合方法通过近场动力学粒子与其域内所有粒子(包括近场动力学粒子和有限元节点)以非局部方式连接,有限元节点与其周围的所有粒子以有限元方式相互作用.将有限元热传导矩阵和近场动力学粒子相互作用矩阵写入同一整体热传导矩阵中,并采用Guyan缩聚法进一步减小计算量.分别采用连续介质力学方法和近场动力学方法对一维以及二维温度场算例进行模拟,结果表明,本文的耦合方法具有较高的计算精度和计算效率.该耦合方案可以进一步拓展到热力耦合条件下含裂纹材料和结构的裂纹扩展问题.     

基于常规态型近场动力学的热力耦合变形破坏分析

基于常规态型近场动力学(Peridynamics, PD)理论,考虑物质点之间的非局部热传导,采用拉格朗日坐标系下的欧拉-拉格朗日方程,推导出常规态型近场动力学热传导方程;结合热力学第一定律和自由能函数建立了常规态型近场动力学热力耦合方程。模拟了含单孔单裂纹和多孔多裂纹板的热传导问题以及典型混凝土构件的热致变形和热力耦合破坏问题,并根据模拟出的温度分布及裂纹扩展情况,验证了本文模型和方法的有效性。基于提出的模型进一步分析了含初始裂缝混凝土三点弯梁的破坏过程及温度梯度和初始裂缝位置对其的影响。将提出的PD热力耦合模型运用于混凝土材料,使PD理论与工程应用相结合。     

近场动力学方法及其应用

近场动力学(peridynamics,PD)是一种新兴的基于非局部作用思想建立模型并通过求解空间积分方程描述物质力学行为的方法.它兼有分子动力学方法和无网格方法的优点,避免了基于连续性假设建模和求解空间微分方程的传统宏观方法在面临不连续问题时的奇异性,又突破了经典分子动力学方法在计算尺度上的局限,在宏/微观不连续力学问题分析中均表现出很高的求解精度和效率.首先概述了PD方法的理论基础、建模思路和计算体系;进而介绍了PD方法在不同尺度不连续力学问题中的应用,包括均匀与非均匀材料和结构的大变形、损伤、断裂、冲击、穿透和失稳问题,结晶相变动力学问题以及纳米材料和结构的破坏问题;最后讨论了PD方法在理论、计算和应用等方面值得进一步研究的问题.     

近场动力学研究进展

近场动力学(Peridynamics或PD)理论基于非局部作用思想,采用空间积分描述物质内部作用,对于从连续到非连续、微观到宏观的力学行为具有统一的表述,数值上天然具有无网格属性和不连续求解功能,在分析不连续,多尺度等问题时展现出了具有优势的适用性和可靠性.本文介绍了近场动力学的发展背景;概述了其理论基础、数值实现过程和计算体系,并在此基础上探讨了近场动力学理论和数值模型的适定性,以及与传统连续介质模型和分子动力学模型进行耦合的可行性;系统分析了近场动力学方法在各个领域上的应用发展现状和趋势,包括静态、动态破坏问题,基于近场动力学的材料模型,以及新兴的疲劳问题研究和多尺度、多物理场的耦合问题;最后对近场动力学方法目前存在的局限性和将来的研究进行了探讨.     

近场动力学在断裂力学领域的研究进展

近场动力学采用非局部积分计算节点内力, 利用统一数学框架描述空间连续与非连续, 避免了非连续区局部空间导数引起的应力奇异, 数值上具有无网格属性, 可自然模拟材料结构的断裂问题. 本文概述了近场动力学的弹性本构力模型, 系统介绍了近场动力学临界伸长率、临界能量密度以及材料强度相关的键失效准则. 详细介绍了近场动力学在断裂力学领域的研究进展, 包括断裂参数能量释放率与应力强度因子的求解、J积分、混合型裂纹、弹塑性断裂、黏聚力模型、动态断裂、材料界面断裂以及疲劳裂纹扩展等. 最后讨论了断裂问题近场动力学研究的发展方向.      

钢筋混凝土结构破坏过程的近场动力学模拟

采用近场动力学(Peridynamics,PD)方法对钢筋混凝土结构破坏过程进行模拟,在"键"型近场动力学模型的基础上,考虑物质点对间的转动以突破泊松比的限制,采用能够描述混凝土材料的拉压异性和断裂特征的损伤模型,引入动态松弛、分级加载、平衡收敛准则和冲击接触等算法,建立了能准确描述钢筋混凝土结构破坏的近场动力学模型。模拟了钢筋混凝土梁在不同荷载作用下破坏的全过程,裂纹扩展路径以及最终破坏形式均与试验结果吻合,为解决工程结构破坏问题提供了一种有效的方法。     

多物理场耦合作用分析的近场动力学理论与方法

广义来说, 近场动力学(peri-dynamics,PD)是假设每个物质点在承受一定范围内的非接触相互作用下,研究整个物理系统演化过程的理论,为涉及非连续和非局部相互作用的问题提供了一个统一的数学框架,具有广泛的适用性.在简要介绍诸多工程对于多物理场模型和数值计算软件的迫切需求后,针对现有商用软件在处理结构非连续演化问题时遇到的瓶颈,引入近场动力学理论和方法. 概述近场动力学固体力学模型,系统阐述近场动力学扩散模型和近场动力学多物理场耦合建模的研究现状和进展,主要涉及电子元器件、电子封装和岩土工程领域的多物理场耦合建模,包括热--力、湿--热--力、热--氧、热--力--氧、力--电、热--电、力--热--电、多孔介质的水--力流固相互作用等非耦合、半耦合与完全耦合模型,强调发展耦合方程数值解法的重要性.最后对扩散问题和多物理场耦合问题的近场动力学理论模型、数值算法和工程应用做进一步展望.      

Review of peridynamics for multi-physics coupling modeling

广义来说, 近场动力学(peri-dynamics,PD)是假设每个物质点在承受一定范围内的非接触相互作用下,研究整个物理系统演化过程的理论,为涉及非连续和非局部相互作用的问题提供了一个统一的数学框架,具有广泛的适用性.在简要介绍诸多工程对于多物理场模型和数值计算软件的迫切需求后,针对现有商用软件在处理结构非连续演化问题时遇到的瓶颈,引入近场动力学理论和方法. 概述近场动力学固体力学模型,系统阐述近场动力学扩散模型和近场动力学多物理场耦合建模的研究现状和进展,主要涉及电子元器件、电子封装和岩土工程领域的多物理场耦合建模,包括热--力、湿--热--力、热--氧、热--力--氧、力--电、热--电、力--热--电、多孔介质的水--力流固相互作用等非耦合、半耦合与完全耦合模型,强调发展耦合方程数值解法的重要性.最后对扩散问题和多物理场耦合问题的近场动力学理论模型、数值算法和工程应用做进一步展望.     

功能梯度材料热传导问题的近场动力学模型

当材料中存在不连续性缺陷的时候,传统的基于连续介质理论的方法在研究热传导问题时存在诸多不便．我们基于近场动力学方法(Peridynamics)建立了功能梯度材料的热传导模型,并且研究了在温度荷载作用下功能梯度材料的温度场的变化表现．该文概述了PD方法应用于热传导问题时的详细理论基础,描述了其建模思路以及计算体系,给出了使用PD方法模拟结构承受温度荷载时的计算格式,讨论了不同梯度形式对功能梯度材料内热传导的影响．算例结果表明：通过PD模型所得到的温度场与解析解吻合较好,证明了PD方法在分析功能梯度材料热传导行为等问题时的可行性．     

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

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

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

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

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

正~~     

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