基于Lasserre算法的自然单元法形函数计算

基于Lasserre体积算法推导了两种插值方案下自然单元法形函数及其导数 的具体计算方法,特别是对计算点处于某些特殊位置时可能造成计算失败的原因和处理方法 进行了较为深入的研究. 算例结果验证了Sibson与non-Sibson插值形函数在三角形外接圆 的圆周上具有不同的连续性,自然单元法形函数在凸区域的边界结点间是线性变化的,因而 可以方便地施加本质边界条件.     

弹性力学的杂交自然单元法

为了解决自然单元法无法直接求解节点应力以及应力解精度不高的问题, 将应力杂交的思想引入自然单元法中, 与弹性问题的Hellinger-Reissner变分原理结合, 提出了弹性问题的杂交自然单元法. 算例表明: 杂交自然单元法的计算结果与解析解吻合, 证明该方法是可行的; 在求解应力方面, 杂交自然单元法比自然单元法有着更高的计算精度, 而且可以直接求解出节点的应力.     

自然单元法在Winkler地基薄板计算中的应用

自然单元法是一种基于Voronoi图及Delaunay三角形剖分图,以自然邻接点插值为试函数的一种无网格数值方法.本文以目前该方法中自然邻接点的Laplace插值形函数为基础,求出了其一阶及二阶导函数,建立了Winkler地基上薄板弯曲挠度的自然单元法求解控制方程,并编制了相应的计算程序,通过算例分析表明了本文方法的可行性和有效性.     

虚拟裂纹闭合法在结构断裂分析中的应用

基于三维虚拟裂纹技术(3DVCCT),利用ABAQUS用户单元子程序(UEL)编写裂纹界面单元,使3DVCCT集成于ABAQUS软件中,直接计算出裂纹的断裂参数.采用此方法对连杆杆身表面裂纹进行研究,得到了连杆裂纹的应力强度因子的分布规律.     

Winkler地基上厚板分析的自然单元法

以自然单元法中自然邻接点的Laplace插值形函数为基础,基于Mindlin厚板理论,建立了Winkler地基上厚板弯曲挠度的自然单元法求解控制方程,并进行了相应的程序实现,最后通过算例分析,表明了该文方法的可行性和有效性.     

自然单元法研究进展

自然单元法是一种基于Voronoi图和Delaunay三角化几何结构,以自然邻点插值为试函数的一种新型数值方法.其既具有无网格方法和经典有限元方法的优点,又克服了两者的一些缺陷,是一种发展前景广阔的求解微分方程的数值方法.自然单元法的形函数满足插值性质,可以像有限元法一样直接施加本质边界条件,不存在基于移动最小二乘拟合的无网格方法不能直接施加本质边界条件的难题.由于自然单元法是无网格方法,可以方便处理有限元方法较难处理的一些问题,例如移动边界和大变形等问题.自然单元法与其他数值方法的最根本区别于其插值格式的不同.将自然邻点插值用于Galerkin过程,就得到基于Voronoi结构的自然单元Galerkin法.自然邻点插值有自然邻点Sibson插值和Laplace插值(非Sibson插值)两种.Laplace插值比Sibson插值在计算上要简单的多,并且不论对凸的或非凸的区域都能精确施加本质边界条件.以Laplace插值为试函数的自然单元法在数值实施上比以Sibson插值为试函数的自然单元法简单.本文对基于Voronoi结构的自然邻点插值和自然单元法的基本思想作了介绍,综述了国内外关于自然单元法的研究成果,总结了自然单元法的优点和尚需解决的问题.      

基于自然单元法的极限上限分析

自然单元法是一种基于离散点集的Voronoi图和Delaunay三角化几何信息,以自然邻近插值为试函数的新型数值方法.相对于一般无网格法中常采用的移动最小二乘近似而言,自然邻近插值不涉及到复杂的矩阵求逆运算,更不需要任何人为的参数,可以提高计算效率.采用该方法构造的形函数满足Delta函数的性质,可以像有限元一样准确地施加边界条件,可以方便处理场函数及其导数的不连续性的问题.论文将自然单元法应用到极限上限分析中,编制了相应的计算程序,通过极限分析的几个经典算例进行了验证,同时采用类似于分片应力磨平的方式,编制相应的磨平程序,由计算点上的塑性耗散功外推得到了节点上的塑性耗散功的值,从而画出了极限状态下结构的塑性耗散功的分布云图.计算结果表明采用自然单元法求解极限上限分析具有稳定性好,精度高,收敛快等优点.     

断裂力学判据的评述

从Inglis 和Griffith 的著名论文到Irwin 和Rice 等的奠基性贡献,对断裂力学中的线弹性断裂力学的K判据,界面断裂力学的G判据,和弹塑性断裂力学的J 判据作了扼要的综述. 介绍了在界面断裂力学G判据的基础上提出的界面断裂力学的K判据,以说明断裂力学的判据存在改进的可能性. 在综述中归纳出断裂力学判据中目前还没有较好解决的几个问题. 在总结以往断裂力学研究经验的基础上,指出裂纹端应力奇异性的源是对断裂力学判据存在的问题作进一步研究的切入点. 探讨了裂纹端应变间断的奇点是裂纹端应力奇异性的源的问题,从而对裂纹端应力强度因子的物理意义进行了讨论. 最后,阐述了进行可靠的裂纹端应力场的弹塑性分析是改进弹塑性断裂力学判据的关键,而进行可靠的裂纹端应力场的弹塑性分析的前提是要通过裂纹端应力奇异性的源的研究来获得作用在裂纹端的造成裂纹端应变间断的有限值应力.      

直接增强自然单元法计算应力强度因子

自然单元法是一种新兴的无网格数值计算方法,但应用于裂纹问题计算时,其近似函数并不能准确反映裂纹尖端渐进应力场的奇异性,为获得足够的计算精度,需要在缝尖附近增大结点的布置密度。针对裂纹问题提出一种增强的自然单元法,将缝尖渐近位移场函数嵌入到自然单元法近似函数中,给出了增强试函数的构造方法,推导了总体刚度矩阵和荷载列阵的相关列式。应力强度因子可以作为附加未知量直接算得,也可用J积分或相互作用能量积分方法进行计算,对增强区域的选择和影响进行了分析。算例结果表明,基于增强自然单元法采用围线积分方法计算应力强度因子具有很高的精度,但直接以附加结点自由度形式计算则精度有所降低。     

人工神经网络在实验断裂力学中的应用

提出一种基于人工神经网络的多裂纹柱体扭转问题的数据处理新方法，运用MATLAB中的神经网络工具箱，对柱体的抗扭刚度实验值进行快速拟合，从而得到了裂纹尖端的应力强度因子，实例应用证明该方法是简便、有效的。     

广义扩展有限元法及其在裂纹扩展分析中的应用

结合广义有限元法(GFEM)和扩展有限元法(XFEM)的特点,提出了一种新的数值方法——广义扩展有限元法(GXFEM)。阐述了广义扩展有限元法的基本原理,对相关公式进行推导,探讨数值实施中需注意的重要问题,给出利用广义扩展有限元法进行断裂分析时应力强度因子的计算方法,编写了广义扩展有限元法程序。通过算例进行了应力强度因子的计算,模拟了结构裂纹的扩展过程。算例结果表明,利用广义扩展有限元法计算裂纹扩展问题,不需要进行过密的网格划分,且网格在裂纹扩展后无需重新剖分,具有相当高的计算精度。     

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method (SBFEM) is a recently developed numerical method combining advantages of both finite element methods (FEM) and boundary element methods (BEM) and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors (SIFs) directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.The project supported by the National Natural Science Foundation of China (50579081) and the Australian Research Council (DP0452681)The English text was polished by Keren Wang.     

Semi-analytical finite element method for fictitious crack model in fracture mechanics of concrete

Based on the Hamiltonian governing equations of plane elasticity for sectorial domain, the variable separation and eigenfunction expansion techniques were employed to develop a novel analytical finite element for the fictitious crack model in fracture mechanics of concrete. The new analytical element can be implemented into FEM program systems to solve fictitious crack propagation problems for concrete cracked plates with arbitrary shapes and loads. Numerical results indicate that the method is more efficient and accurate than ordinary finite element method.     

Photoelastic experimental hybrid method for fracture mechanics of anisotropic materials

In this paper, the isotropic and anisotropic photoelastic experimental hybrid methods for fracture mechanics are developed. Using the photoelastic experimental hybrid method, it is demonstrated that one can precisely obtain stress intensity factors and separate the stress components of isotropic and anisotropic plate problems from the only isochromatics.     

自然单元法的自适应研究

自然单元法(NEM)是一种新兴的无网格数值计算方法,具有前处理简单和易于准确施加本质边界条件等优点.本文基于Z-Z后验误差估计方法,给出了一种自然单元法的误差估计因子和自适应分析细化方案.采用结点处的光滑应变计算相应的恢复应力,并用于构造全域上的恢复应力场.通过结点Voronoi单胞内的能量范数相对误差指示需要进行结点...     

基于有限断裂法和比例边界有限元法的裂纹扩展模拟

基于有限断裂法和比例边界有限元法提出了一种裂缝开裂过程模拟的数值模型。采用基于有限断裂法的混合断裂准则作为起裂及扩展的判断标准,当最大环向应力和能量释放率同时达到其临界值时,裂缝扩展。结合多边形比例边界有限元法,可以半解析地求解裂尖区域附近的应力场和位移场,在裂尖附近无需富集即可获得高精度的解。计算能量释放率时,只需将裂尖多边形内的裂尖位置局部调整,无需改变整体网格的分布,网格重剖分的工作量降至最少。裂缝扩展步长通过混合断裂准则确定,避免了人为假设的随意性,并可以实现裂缝变步长扩展的模拟,更符合实际情况。通过对四点剪切梁的复合型裂缝扩展过程的模拟,对本文模型进行了验证,并应用于重力坝模型的裂缝扩展模拟,计算结果表明,本文提出的模型简单易行且精度较高。     

扩展有限元裂尖场精度研究

论述了扩展有限元方法和基本原理,研究了单元类型(四边形单元和三角形单元、线性单元和二次单元)、网格密度、J积分区域半径等因素对裂尖局部应力场(应力强度因子)计算精度的影响。研究发现,上述因素对裂尖应力强度因子计算的收敛速度与稳定性影响不大,证实了XFEM可以用较少的节点获得较高的裂尖场精度,并提出了通过固定裂尖附加区半径可以进一步改善XFEM的收敛速度。