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

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

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

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

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

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

各向异性对平板裂纹扩展路径的影响研究

近场动力学理论克服了经典连续介质力学在模拟裂纹扩展方面的不足。本文基于微梁键构建了一般各向异性近场动力学平面模型,并探讨了各向异性对于裂纹扩展路径的影响。首先,建立一般各向异性的微梁键,并基于插值法建立了一般各向异性近场动力学平面模型的应变能密度表达式;然后,利用应变能密度互等原理求解出近场动力学平面本构模型的一般各向异性参数;最后,通过带初始切口的平板断裂实验研究了不同各向异性参数对于平板裂纹扩展路径的影响。     

近场动力学算子方法

提出一种基于非局部思想求解物理学问题的近场动力学算子方法 (peridynamic operator method, PDOM).运用PDOM可将任意阶局部微分及其乘积转化为相应的非局部积分形式,且无需额外地特殊处理间断点与奇异点等问题.近年来研究较多的两种非局部算子:近场动力学微分算子(peridynamic differential operator,PDDO)和非局部算子方法 (nonlocal operator method, NOM),均可视为PDOM的一种特例.以弹性力学问题为例,采用变分原理和拉格朗日方程,建立了适用于分析静/动态弹性力学问题的PDOM模型.理论分析表明,当分别限定相互作用域为与位置无关或位置相关的圆形域时,该PDOM弹性模型即可退化为近年来文献中常见的近场动力学(peridynamics, PD)模型或对偶域近场动力学(dual-horizon peridynamics, DH-PD)模型.通过3个典型实例:杆的拉伸与波动、亥姆霍兹方程和含孔板的拉伸,说明本方法的计算精度、收敛性与数值稳定性.PDOM方法适用于任意均匀或非均匀离散,且能有效避免零能模式以...     

动载作用下复合型裂纹扩展的近场动力学模拟

构建改进的非局部近场动力学模型和粒子系统动力加载算法,开展复合型裂纹动力扩展过程模拟。在常规近场动力学微极模型中引入能反映非局部长程力尺寸效应的核函数修正项,以提高计算精度和收敛稳定性。通过模拟动载作用下含复合型裂纹混凝土三点弯梁的裂纹扩展与破坏过程并与试验结果对比,验证了本文模型和算法的可靠性。在此基础上,通过分组模拟分析了动载作用下,初始裂纹的位置、长度和方向对含复合型裂纹三点弯梁破坏模式、起裂时间与承载能力的影响规律。     

非常规态型近场动力学热黏塑性模型及其应用

在非常规态型近场动力学(non-ordinary state-based peridynamics, NOSB-PD) 理论框架下构建了考虑应变率效应、塑性硬化、热软化效应和材料断裂特征的非局部三维热黏塑性固体本构模型以及相应的非局部空间积分型数值算法, 并应用于金属类材料和构件在冲击载荷作用等工况下的高应变率热黏塑性变形与破坏分析. 通过对经典含初始裂纹Kalthoff-Winkler板冲击试验进行三维近场动力学模拟, 可得到裂纹的起裂角度、扩展路径、扩展速度以及裂纹扩展过程中靶板等效应力和温度分布, 所得结果与已有试验结果和其他数值方法结果吻合较好. 在此基础上, 应用该模型分析了不同冲击速度作用下金属靶板的变形与裂纹扩展过程, 结果表明: 该模型能较好地模拟不同冲击速度(应变率)情况下靶板的变形与破坏全过程. 随着冲击速度变化, 初始裂纹的起裂时间、扩展方向和扩展速度呈一定规律变化. 冲击速度越低, 起裂时间越晚(直至冲击速度低于某值时初始裂纹不扩展), 裂纹扩展速度峰值越低, 冲击过程中靶板温度峰值越低, 完全扩展所需时间越长.      

非常规态型近场动力学热黏塑性模型及其应用

在非常规态型近场动力学(non-ordinary state-based peridynamics,NOSB-PD)理论框架下构建了考虑应变率效应、塑性硬化、热软化效应和材料断裂特征的非局部三维热黏塑性固体本构模型以及相应的非局部空间积分型数值算法,并应用于金属类材料和构件在冲击载荷作用等工况下的高应变率热黏塑性变形与破坏分析.通过对经典含初始裂纹Kalthoff-Winkler板冲击试验进行三维近场动力学模拟,可得到裂纹的起裂角度、扩展路径、扩展速度以及裂纹扩展过程中靶板等效应力和温度分布,所得结果与已有试验结果和其他数值方法结果吻合较好.在此基础上,应用该模型分析了不同冲击速度作用下金属靶板的变形与裂纹扩展过程,结果表明:该模型能较好地模拟不同冲击速度(应变率)情况下靶板的变形与破坏全过程.随着冲击速度变化,初始裂纹的起裂时间、扩展方向和扩展速度呈一定规律变化.冲击速度越低,起裂时间越晚(直至冲击速度低于某值时初始裂纹不扩展),裂纹扩展速度峰值越低,冲击过程中靶板温度峰值越低,完全扩展所需时间越长.     

基于近场动力学的功能梯度材料动态断裂分析

近场动力学(Peridynamics)是一种基于非局部理论的数值计算方法,通过域内积分建立物质基本运动方程。本文根据功能梯度材料的特点提出了一种基于近场动力学理论的动态断裂分析方法,模拟了在冲击载荷作用下功能梯度材料的动态响应,并分析了微模量函数对裂纹扩展行为的影响。通过与已有实验结果的对比发现,近场动力学模型所得裂纹扩展路径及破坏时间均与实验结果吻合较好,从而证明近场动力学方法能够用于分析功能梯度材料的动态断裂行为。     

基于态型近场动力学的淬火陶瓷破坏过程分析

基于态型近场动力学理论进行复合陶瓷板淬火破坏过程力学建模与分析。引入热膨胀项反映材料的热致变形,基于物质点对的断裂定义损伤,构建三维非局部常规态型近场动力学热弹脆性模型,结合多速率显式积分法计算热-力耦合,实现热冲击荷载作用下陶瓷板起裂和裂纹扩展全过程的模拟。通过模拟典型淬火Al₂O₃陶瓷板的破坏过程,并与实验和其他数值结果对比,验证了该模型和方法。使用该模型研究了含不同α-Al₂O₃纤维掺杂比率的Al₂O₃复合陶瓷板淬火破坏过程,结果表明,断裂韧性的提高显著减少了裂纹数量并延后起裂时间;含5%α-Al₂O₃纤维掺杂比率的复合陶瓷材料综合力学性能最优,能够有效提升构件抵抗热冲击荷载能力。     

On the asymptotic crack-tip stress fields in nonlocal orthotropic elasticity

The crack-tip stress fields in orthotropic bodies are derived within the framework of Eringen’s nonlocal elasticity via the Green’s function method. The modified Bessel function of second kind and order zero is considered as the nonlocal kernel. We demonstrate that if the localisation residuals are neglected, as originally proposed by Eringen, the asymptotic stress tensor and its normal derivative are continuous across the crack. We prove that the stresses attained at the crack tip are finite in nonlocal orthotropic continua for all the three fracture modes (I, II and III). The relative magnitudes of the stress components depend on the material orthotropy. Moreover, non-zero self-balanced tractions exist on the crack edges for both isotropic and orthotropic continua. The special case of a mode I Griffith crack in a nonlocal and orthotropic material is studied, with the inclusion of the T-stress term.     

考虑相变的近场动力学热?力耦合模型及多孔介质冻结破坏模拟

多孔介质的传热传质现象广泛存在于自然界和工业领域中. 低温条件可能导致多孔介质中的组分发生相变, 并由此诱发材料损伤, 甚至导致结构失效破坏. 对这类破坏现象的预测需要精细化建模, 以能够反映物质的相变过程和材料的破坏特征. 本文采用热焓法改写经典的热传导方程, 在近场动力学框架下, 建立了一种考虑物质相变的热?力耦合模型, 发展了交错显式求解的数值计算方法, 进行了方板角冻结、热致变形和多孔介质冻结破坏等问题的模拟, 得到了方板的冻结特征、温度场和变形场的分布规律以及多孔介质的冻结破坏过程, 与试验和其他数值方法的结果具有较好的一致性. 研究表明, 本文所建立的考虑物质相变的近场动力学热?力耦合模型能够反映材料的非局部效应和物质相变潜热的影响, 准确捕捉相变过程中液固界面的演化特征, 再现多孔介质中材料相变、基质热致变形和冻结破坏过程, 突破了传统连续性模型求解这类破坏问题时面临的瓶颈, 为深入研究多孔介质冻融破坏过程和破坏机理提供了有效途径.      

混凝土复合型裂纹扩展的非局部近场动力学建模分析

基于非局部近场动力学理论,构建了修正的能反映混凝土宏观拉压异性和断裂特征的近场动力学本构模型,开发了相应的离散、加载和时间积分算法,实现典型混凝土构件中复合型裂纹扩展过程模拟。在物质点对尺度上定义局部损伤并考虑物质点对的相对转动,通过求解时空微-积分方程实现裂纹的自然萌生与扩展,避免裂尖不连续带来的求解奇异性、网格依赖性和网格重构以及常规近场动力学本构模型的泊松比限制。通过含单边和双边初始裂纹四点剪切混凝土梁裂纹扩展破坏全过程模拟,得到破坏形态、破坏荷载以及完整的荷载-裂纹开口滑移曲线,并与试验和其他数值模拟结果对比,验证了模型的精确性和算法的稳定性。     

Heat conduction and thermal stress induced by an electric current in an infinite thin plate containing an elliptical hole with an edge crack

A uniform electric current at infinity was applied to a thin infinite conductor containing an elliptical hole with an edge crack. The electric current gives rise to two states, i.e., uniform and uneven Joule heat. These two states must be considered to analyze the heat conduction problem. The uneven Joule heat gives rise to uneven temperature and thus to heat flux, and to thermal stress.Using a rational mapping function, problems of the electric current, the Joule heat, the temperature, the heat flux, the thermal stress are analyzed, and each of their solutions is obtained as a closed form. The distributions of the electric current, the Joule heat, the temperature, the heat flux and the stress are shown in figures.The heat conduction problem is solved as a temperature boundary value problem. Solving the thermal stress problem, dislocation and rotation terms appear, which complicates this problem. The solutions of the Joule heat, the temperature, the heat flux and the thermal stress are nonlinear in the direction of the electric current. The crack problems are also analyzed, and the singular intensities at the crack tip of each problem are obtained. Mode II (sliding mode) stress intensity factor (SIF) is produced as well as Mode I (opening mode) SIF, for any direction of the electric current. The relations between the electric current density and the melting temperature and between the electric current density and SIF are investigated for some crack lengths in an aluminum plate.     

Nonlocal Constrained Value Problems for a Linear Peridynamic Navier Equation

In this paper, we carry out further mathematical studies of nonlocal constrained value problems for a peridynamic Navier equation derived from linear state-based peridynamic models. Given the nonlocal interactions effected in the model, constraints on the solution over a volume of nonzero measure are natural conditions to impose. We generalize previous well-posedness results that were formulated for very special kernels of nonlocal interactions. We also give a more rigorous treatment to the convergence of solutions to nonlocal peridynamic models to the solution of the conventional Navier equation of linear elasticity as the horizon parameter goes to zero. The results are valid for arbitrary Poisson ratio, which is a characteristic of the state-based peridynamic model.     

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

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

Static and dynamic behaviour of nonlocal elastic bar using integral strain-based and peridynamic models

The static and dynamic behaviour of a nonlocal bar of finite length is studied in this paper. The nonlocal integral models considered in this paper are strain-based and relative displacement-based nonlocal models; the latter one is also labelled as a peridynamic model. For infinite media, and for sufficiently smooth displacement fields, both integral nonlocal models can be equivalent, assuming some kernel correspondence rules. For infinite media (or finite media with extended reflection rules), it is also shown that Eringen's differential model can be reformulated into a consistent strain-based integral nonlocal model with exponential kernel, or into a relative displacement-based integral nonlocal model with a modified exponential kernel. A finite bar in uniform tension is considered as a paradigmatic static case. The strain-based nonlocal behaviour of this bar in tension is analyzed for different kernels available in the literature. It is shown that the kernel has to fulfil some normalization and end compatibility conditions in order to preserve the uniform strain field associated with this homogeneous stress state. Such a kernel can be built by combining a local and a nonlocal strain measure with compatible boundary conditions, or by extending the domain outside its finite size while preserving some kinematic compatibility conditions. The same results are shown for the nonlocal peridynamic bar where a homogeneous strain field is also analytically obtained in the elastic bar for consistent compatible kinematic boundary conditions at the vicinity of the end conditions. The results are extended to the vibration of a fixed–fixed finite bar where the natural frequencies are calculated for both the strain-based and the peridynamic models.     

Analytic solutions to crack tip plastic zone under various loading conditions

Based on stress field equations and Hill yield criterion, the crack tip plastic zone is determined for orthotropic materials and isotropic materials under small-scale yielding condition. An analytical solution to calculating the crack tip plastic zone in plane stress states is presented. The shape and size of the plastic zone are analyzed under different loading conditions. The obtained results show that the crack tip plastic zones present “butterfly-like” shapes, and the elastic–plastic boundary is smooth. The size of the plastic zone for orthotropic composites is less at the crack tip for various loading conditions, compared with the case of isotropic materials. Crack inclination angle and loading conditions affect greatly the size and shape of crack tip plastic zone. The mode I crack has a crucial effect on the plastic zone for mixed mode case in plane stress state. The plastic zone for pure mode I crack and pure mode II crack have a symmetrical distribution to the initial crack plane.     

Buckling analysis of functionally graded nanobeams under non-uniform temperature using stress-driven nonlocal elasticity

In this work, the size-dependent buckling of functionally graded(FG)Bernoulli-Euler beams under non-uniform temperature is analyzed based on the stressdriven nonlocal elasticity and nonlocal heat conduction. By utilizing the variational principle of virtual work, the governing equations and the associated standard boundary conditions are systematically extracted, and the thermal effect, equivalent to the induced thermal load, is explicitly assessed by using the nonlocal heat conduction law. The ...