含曲线裂纹圆柱扭转问题的新边界元法

研究含曲线裂纹圆柱的Saint-Venant扭转,将问题化归为裂纹上边界积分方程的求解．利用裂纹尖端的奇异元和线性元插值模型,给出了扭转刚度和应力强度因子的边界元计算公式．对圆弧裂纹、曲折裂纹以及直线裂纹的典型问题进行了数值计算,并与用Gauss-Chebyshev求积法计算的直裂纹情形结果进行了比较,证明了方法的有效性和正确性．     

双材料界面裂纹平面问题的半权函数法

应用半权函数法求解双材料界面裂纹的平面问题．由平衡方程、应力应变关系、界面的连续条件以及裂纹面零应力条件推导出裂尖的位移和应力场,其特征值为lambda及其共轭．设置特征值为lambda的虚拟位移和应力场,即界面裂纹的半权函数A·D2由功的互等定理得到应力强度因子KⅠ和KⅡ以半权函数与绕裂尖围道上参考位移和应力积分关系的表达式．数值算例体现了半权函数法精度可靠、计算简便的特点．     

短纤维复合材料的本征应变边界积分方程计算模型

提出了短纤维复合材料的本征应变边界积分方程计算模型,并采用新发展的边界点法进行了求解．模型依据Eshelby等效夹杂物的概念并借助Eshelby张量,采用迭代方法来计算基体中各种性能短纤维的本征应变,其中所需的Eshelby张量不难通过解析或数值方法获得．由于未知量只出现在边界上,与已有的有限元和边界元模型相比,提出的计算模型可极大地减小异质体问题的求解规模,提高计算效率．通过数值算例计算了代表性体积单元上各种短纤维复合材料的整体弹性性能,验证了计算模型和求解方法的正确性和有效性．     

含圆形孔口和水平直裂纹的平面问题的应力分析

利用复变函数方法和积分方程理论研究了既含有圆形孔口又含有水平裂纹的无限大平面的平面弹性问题,将复杂的解析函数的边值问题化成了求解只在裂纹上的奇异积分方程的问题.此外,还给出了裂纹尖端附近的应力场和应力强度因子的公式.     

平面弹性裂纹分析的一种有效边界元方法

提出了一种简单而有效的平面弹性裂纹应力强度因子的边界元计算方法．该方法由Crouch与Starfield建立的常位移不连续单元和闫相桥最近提出的裂尖位移不连续单元构成A·D2在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界．算例(如单向拉伸无限大板中心裂纹、单向拉伸无限大板中圆孔与裂纹的作用)说明平面弹性裂纹应力强度因子的边界元计算方法是非常有效的．此外,还对双轴载荷作用下有限大板中方孔分支裂纹进行了分析．这一数值结果说明平面弹性裂纹应力强度因子的边界元计算方法对有限体中复杂裂纹的有效性,可以揭示双轴载荷及裂纹体几何对应力强度因子的影响．     

压电材料中两平行不相等界面裂纹的动态特性研究

利用Schmidt方法,研究了压电材料中两个平行不相等的可导通界面裂纹对简谐反平面剪切波的散射问题．利用Fourier变换,使问题的求解转换为对两对以裂纹面张开位移为未知变量的对偶积分方程的求解．数值计算结果表明,动态应力强度因子及电位移强度因子受裂纹的几何参数、入射波频率的影响．在特殊情况下,与已有结果进行了比较分析．同时,电位移强度因子远小于不可导通电边界条件下相应问题的结果．     

平面应变Ⅰ型裂纹弹塑性分析的一种方案

本文给出了一种求解平面应变Ⅰ型裂纹弹塑性近似解析解的一种方法.     

法向荷载作用下椭圆孔不等边裂纹的应力强度因子

采用复变函数方法,研究了在法向均布荷载作用下,含两个不等边裂纹椭圆孔的无限大板平面问题,得到了裂纹尖端的应力强度因子的解析解.并通过有限元软件计算了应力强度因子的数值解,与解析解进行对比,吻合较好.另外,研究了随着裂纹和椭圆孔尺寸变化时应力强度因子的变化规律.可以看出应力强度因子随椭圆孔的长短半轴之比和裂纹长度的增大而增大.     

集中力作用下的半平面边缘裂纹问题

讨论了半平面中的倾斜裂纹问题。集中力作用于裂纹表面上,或作用于开裂半平面的边界上。利用有理保角映像函数方法求解这个问题,同时得到了闭合形式的解。最后,给出了二个数值例子和计算结果。     

单裂纹问题的虚边界无网格伽辽金法分析

使用含裂纹复变基本解,虚边界无网格伽辽金法被进一步推广应用于弹性材料的单裂纹问题求解.为了清晰地说明单裂纹问题的虚边界元法实现过程,单裂纹问题的虚边界元法示意图、复变坐标平面下含裂纹问题的复变位移和复变面力基本解示意图被展示.含裂纹复变基本解,因自动满足裂纹处边界条件,故使用虚边界无网格伽辽金法计算单裂纹问题,无需在裂纹处布置节点或单元.给出含裂纹复变基本解中的Φ'(x)的详细表达式、裂纹左右裂尖应力强度因子的虚边界无网格离散公式,方便了其他学者使用本方法计算裂纹问题.数值计算两端受拉长方形钢板中心含有裂纹的应力强度因子的算例,计算结果证明了本方法的精确性与稳定性.     

A new model for the analysis of reinforced concrete members with a coupled HdBNM/FEM

As a truly boundary-type meshless method, the hybrid boundary node method (HdBNM) does not require ‘boundary element mesh’, either for the purpose of interpolation of the solution variables or for the integration of ‘energy’. In this paper, the HdBNM is coupled with the finite element method (FEM) for predicting the mechanical behaviors of reinforced concrete. The steel bars are considered as body forces in the concrete. A bond model is presented to simulate the bond-slip between the concrete and steels using fictitious spring elements. The computational scale and cost for meshing can be further reduced. Numerical examples, in 2D and 3D cases, demonstrate the efficiency of the proposed approach.     

Research on the Fictitious Source Points of the Hybrid Fundamental Solution⁃Based Finite Element Method for Heat Conduction Problems北大核心CSCD

针对热传导问题,提出了杂交基本解有限元法.首先,假设两个独立场:一个为利用基本解线性组合近似的单元域内温度场,另一个为使用与传统有限元法相同形式的辅助网线温度场.然后,利用修正变分泛函将上述两个独立场关联起来,并导出有限元列式.然而,该方法的准确性很大程度上取决于源点的分布和数量,通常将源点布置在单元外部两种虚拟边界上:与单元相似的边界和圆形边界.此外,还提出了双重虚拟边界,并与上述两种源点布局方式进行对比.通过两个典型数值算例,验证了该文方法在不同源点布局下的有效性和对网格畸变的不敏感性.     

Effect of shear stiffness and concrete wedge cone failure on fracture of a fiber-reinforced plastic rod at an intersection with a concrete crack

Experimental results of model speciments in which FRP rods fractured due to local deformation at a crack intersection in a concrete member were analyzed by a 3D nonlinear finite element method in which orthogonal anisotropy of the FRP rod was considered. The analytical results indicated that accurate prediction of shear modulus of the FRP rod and size of concrete wedge cone failure around the FRP rod was significant to predict deformation and fracture of the FRP rod. FRP rods as reinforcement in concrete members, the small shear modulus, because of the orthogonal anisotropy and the wedge cone failure, may prevent the FRP rod from fracturing at a very low tensile stress due to the local deformation at the crack intersection.Presented at the Ninth International Conference on the Mechanics of Composite Materials, Riga, October, 1995.Published in Mekhanika Kompozitnykh Materialov, Vol. 21, No. 2, pp. 158–166, March–April, 1996.     

Hexahedral finite element with an open crack

A method of creating the stiffness matrix of a hexahedral eight-node finite element with a single, nonpropagating, transverse, one-edge crack at half of its length is presented in this paper. The crack was modelled by adding an additional flexibility matrix to that of the noncracked element. The terms of the additional matrix have been calculated by use of the laws of fracture mechanics. Employing the elaborated element a numerical test has been worked out, the results of which are compared with the data of analytical solutions accessible in the literature, and a high conformity with them has been obtained. The element presented in the paper may be applied to the static and dynamic analysis of different types of structural elements with material defects in the form of cracks. The described method of creating the stiffness matrix of the element allows to create different kinds of finite elements with cracks provided that the stress intensity factors for a given type of crack are known.     

线性蠕变平面问题中的摄动有限元方法

本文提出一种利用摄动有限元求解线性蠕变问题的方法(PFMC).它可用于平面蠕变问题,诸如钢筋混凝土梁,预应力混凝土梁,钢筋混凝土圆筒或位于弹性或粘弹性介质中的钢筋混凝土隧洞及地下建筑等.本方法不采用一般增量法中关于在一个时段内各物理量保持不变的假设,提高了分析精度,加大了计算步长,减少了机器存储,提高了运算效率.本文构造了包含钢筋在内的四节点四边形等参单元的有限元摄动格式.并给出五个算例,与解析解相比,有令人满意的精度.     

Improved Convergence Results for the Finite Element Method with Holomorphic Functions

The purpose of this paper is to present extended results on convergence of a method for coupling of an analytical and a finite element solution for a boundary value problem with a singularity. As an example of a singular problem we consider a solution to the Lamé-Navier equation with a singularity caused by a crack. The main idea of such an approach is to construct a continuous coupling between an analytical solution near a singularity and a finite element solution through the whole interaction interface. In this paper we present a basic theory of convergence for a method of coupling, particularly we discuss some points which show a difference to the standard theory of the finite element method.     

Crack identification of the sprocket wheel in sintering machine based on WFEM

In order to identify the crack of the sprocket wheel correctly, the wavelet finite element method is studied in depth. Firstly, the progress of study on the wavelet finite element method is summarized, and then the basic property of wavelet analysis is analyzed, and then the wavelet finite element theory model of sprocket wheel in sintering machine is studied, and the Daubechies wavelet plate element and isoparametric plate element of crack tip are established, and then the theory of constructing crack identification database of the crack for sprocket wheel is studied, and finally the effective of this method is verified by identify the sprocket wheels with five kinds of cracks based on the vibration test, and results showed that this method can identify the crack of the sprocket wheel correctly.     

考虑流固耦合效应的重力坝水力劈裂模拟

裂缝的高压水力劈裂是混凝土高坝安全评估的重要部分,研究其过程中的流固耦合作用是准确预测在各种情况下裂纹扩展路径和危险程度的关键.该文利用扩展有限元法在模拟裂纹扩展方面的优势,对大坝的裂纹进行水力劈裂模拟研究.裂纹中的水压分布模型采用Brühwiler和Saouma水力劈裂试验的成果,体现了水压和裂纹宽度的耦合关系,给出了扩展有限元在裂纹面上施加水压力荷载的实施方法,对一典型重力坝裂纹的水力劈裂进行了数值模拟分析.研究结果表明:采用扩展有限元法模拟水力劈裂,克服了常规有限元法存在的缺点,裂纹扩展时不用重新划分网格,裂纹的实时宽度可以由加强节点的附加自由度得到,裂纹面上水压的施加也变得简单易行.当考虑裂纹内的流固耦合效应时,裂纹的扩展路径相比不考虑耦合效应时的扩展路径(均布全水头水压),扩展角变大,扩展距离变短.     

Analysis of optimal preconditioners for CutFEM

In this article, we consider a class of unfitted finite element methods for scalar elliptic problems. These so-called CutFEM methods use standard finite element spaces on a fixed unfitted triangulation combined with the Nitsche technique and a ghost penalty stabilization. As a model problem we consider the application of such a method to the Poisson interface problem. We introduce and analyze a new class of preconditioners that is based on a subspace decomposition approach. The unfitted finite element space is split into two subspaces, where one subspace is the standard finite element space associated to the background mesh and the second subspace is spanned by all cut basis functions corresponding to nodes on the cut elements. We will show that this splitting is stable, uniformly in the discretization parameter and in the location of the interface in the triangulation. Based on this we introduce an efficient preconditioner that is uniformly spectrally equivalent to the stiffness matrix. Using a similar splitting, it is shown that the same preconditioning approach can also be applied to a fictitious domain CutFEM discretization of the Poisson equation. Results of numerical experiments are included that illustrate optimality of such preconditioners for the Poisson interface problem and the Poisson fictitious domain problem.     

Formulation and Application of the Initially Rigid Cohesive Finite Element Method

The presented procedure for cohesive crack propagation is based on an adaptive finite element (FE) implementation, which enables the introduction of cohesive surfaces in dependence on the current crack state. In contrast to already existing formulations, the focus of the present model lies on failure processes that can be described at quasi-static conditions within an implicit framework. Furthermore, an extension for mesh independent crack propagation in terms of an additional mesh adaptive formulation is presented. By the evaluation of the failure criterion considering the preferred crack direction, a new crack tip coordinate is computed and the discretization is accordingly adjusted. The remaining mesh is modified for the new boundary representation. The application of the proposed method is shown by the numerical investigation of a concrete fracture specimen from an experimental research project. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)