Mode III fracture behavior of laminated composite with edge crack in torsion

A three-dimensional (3-D) finite element analysis was performed on a [90,(+45/−45)_n,(−45/+45)_n,90]_s class of laminated composites under the edge crack torsion (ECT) test configuration. Finite element delamination models were established and formulas for calculating the Mode III fracture toughness from 3-D finite element models were developed. The relations between the interlaminar fracture behavior and various configuration parameters were investigated and the effects of point loads, ends, geometry, Mode II interference, and friction were evaluated. Results showed that with proper selection of ECT specimen configuration and layup, the delamination could grow in pure Mode III in the middle region of the specimen. Specimen end effect played an important role in the ECT test. A Mode II component occurred in the end regions but it did not interfere significantly with the Mode III delamination state. Specimen dimension ratio, layup, and crack length exhibited significant effect on the interlaminar fracture behavior and the calculated strain energy release rates. However, friction between crackfaces was found to have negligible effect on the interlaminar properties.     

Dual variational formulation for Trefftz finite element method of elastic materials

A dual variational principle is presented for Trefftz finite element analysis. The proof of the stationary conditions of the variational functional and the theorem on the existence of extremum are provided in this paper. They are boundary displacement condition, surface traction condition and interelement continuity condition. Based on the assumed intraelement and frame fields, element stiffness matrix equation is obtained which can easily be implemented into computer programs for numerical analysis with Trefftz finite element method. Two numerical examples are considered to illustrate the effectiveness and applicability of the proposed element model.     

流形Trefftz直接法

基于数值流形方法,结合广义节点的概念,构造了一类节点位移用任意阶多项式展开的流形Trefftz直接法。这种方法融合了数值流形方法和Trefftz直接法的优点。计算结果表明,流形Trefftz直接法与Trefftz直接法以及其他计算方法相比有较高的精度。     

Trefftz method for hydrodynamic pressure on rigid dams with non-vertical upstream face

Based on a two-dimensional potential flow theory, earthquake-induced hydrodynamic pressures on a rigid dam with a non-vertical upstream face are examined by the Trefftz method. The effect of surface waves on the hydrodynamic pressure distribution is discussed in detail. Numerical values are given for different wave effect parameters and different geometries of the dam–water interface.     

Trefftz间接法解自由面渗流问题

渗流问题是岩体水力学研究的重点之一,较常规的有限元法和传统的边界元法而言,Trefftz型边界元法具有程序简单、计算量小及无须奇异积分等特点。本文给出定常和非定常自由面渗流问题的Trefftz解法,计算结果表明,该法收敛性好,结果比较精确。     

随机有限元方法在断裂分析中的应用

在幂律非线性随机有限元基础上,以单边裂纹板为例给出计算含量钢继裂参数,J（J积分）,δ（裂纹张开位移）,Δ（由裂纹引起的裂纹板上下底面相对位移）,θ（由裂纹引起的裂纹板上下底在相对转角）及其对基本随机变量变化率的方法和分析算例。     

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

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

Coarse-graining atomistic dynamics of brittle fracture by finite element method

In this work we present the finite element (FE) implementation of an atomistic formulation of balance equations and its application to coarse-grained (CG) simulation of dynamic fracture. First, we simulate a notched specimen that contains about 1.8 million atoms by the CG-FE method, and we compare the CG-FE results with that by all-atom molecular dynamics (MD) simulations. We find that CG-FE simulations with about 5% degrees of freedom of the MD simulation can capture the essential dynamic features, not in exact correspondence, but qualitatively and quantitatively similar to that obtained by MD simulations. We then proceed to simulate a series of micron-sized specimens by the CF-FE method. We find that it is the interaction of the forward propagating crack with the stress waves being reflected back by the boundaries of the specimen that triggers the dynamic instability and hence the branching of cracks in micron-sized specimens. The potential application of the method and future work for improvements are discussed.     

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.     

Damage analysis by boundary element method

The response of cracked bodies subjected to loading was investigated by the boundary element method in this paper. The two-law elastic-cohesive-softening model was used for crack propagation analysis. The interface conditions for uncracked, craze, open crack, adhesive crack and slid crack parts were discussed and the corresponding incremental iteration algorithm was given. A simplified damage propagation model was presented. The technique has been applied to some specific examples which give the evidence that the method is satisfactory and efficient.     

Electro-mechanical coupling analysis of MEMS structures by boundary element method

In this paper, we present the applications of Boundary Element Method (BEM) to simulate the electro-mechanical coupling responses of Micro-Electro-Mechanical systems (MEMS). The algorithm is programmed in our research group based on BEM modeling for electrostatics and elastostatics. Good agreement is shown while the simulation results of the pull-in voltages are compared with the theoretical/experimental ones for some examples. The project supported by the 973 Program (G199033108) and the national Natural Science Foundation of China (10125211)     

用离散单元法研究混凝土、岩石类材料的拉剪混合型断裂

基于变形体离散单元法,采用固定和旋转的弥散裂缝模型,分析了混凝土、岩石等准脆性材料的拉剪混合型开裂行为,给出了开裂准则以及裂缝扩展机制,初步探讨了从连续介质到非连续介质转化的数值模拟。通过Willam数值试验表明旋转裂缝模型在模拟拉剪混合型开裂问题时更符合客观实际;通过单边切口非对称三点弯梁试验分析,说明在混合形式的荷载条件下,裂缝通常以拉剪Ⅰ/Ⅱ型混合模式起裂,而后以Ⅰ型为控制方式稳定扩展。     

基于单元子分法的结构多尺度边界单元法

建立在基于单元子分法的一种有效自适应格式以及多区域边界元三步求解技术基础上提出了一种计算结构多尺度问题的多区域边界元法。首先,通过高斯积分误差分析公式确定边界单元在满足精度要求下所需要的高斯点数,当所需高斯点数超过规定数目时该单元就被自动划分成一定数量的子单元,从而消除结构多尺度所引起的近奇异性。在单元子分技术的基础上采用多区域边界元三步求解技术来处理材料非均质问题:第一步消除各子域的内部未知量,第二步消除各子域独自拥有的边界未知量,第三步根据位移相容性条件和面力平衡条件建立系统方程组并求解公共界面节点位移以及每个子域的其他未知量。数值算例结果表明本方法可以用较少的计算时间得到满意的结果,是处理结构多尺度问题的一种有效方法。     

基于离散元的杆件断裂全过程数值模拟

杆件的断裂会涉及到大变形、非线性以及不连续等问题,通常的数值计算方法模拟这种复杂力学行为具有局限性。本文基于颗粒离散元法DEM,将接触粘结处的分布式弹簧用梁纤维进行等效,提出了一种适于结构弹塑性分析的DEM纤维梁模型,然后在此基础上构建了构件断裂模拟算法以及纤维破环准则。将该模型应用于悬臂梁结构,模拟了悬臂梁从弹性到弹塑性阶段,再到断裂破坏的全过程,数值模拟得到的结构响应和截面开裂破坏形态均较合理。最后将该方法应用于单层网壳倒塌破坏模拟,并与网壳振动台倒塌试验进行对比,结果表明,数值模拟得到的杆件断裂过程及结构倒塌模式与试验现象一致,验证了该模型的正确性和适用性。     

板壳问题的三维无网格伽辽金直接分析法

对于板壳问题,共有三种数值模拟方案:线性或非线性的板壳理论、退化连续体方案和直接三维连续体方案。无网格法近似函数可具有C1甚至更高的连续性,便于在K irchhoff-Love理论中应用。但当各种无网格法用于M ind lin-R e issner板理论时,会遇到数值锁死的困扰。对比之下,三维连续体方案是最简单,最精确但并不常用的一种方案。无网格法近似函数具有高度光滑性,在板壳的厚度方向仅布置2~5层点就可以很好地捕捉此方向场的梯度,同时还可以在一定参数范围内避免剪切和体积锁死,在处理复杂本构关系、非线性板壳等问题中更是具有很大优势。本文采用无网格伽辽金法(EFG)和三维连续体方案分析了线性板壳问题,与有限单元法做了对比,并讨论了数值锁死等问题。     

基于等几何分析的比例边界有限元方法

提出了一种具有比例边界有限元的半解析特性和等几何分析的几何特性的新方法。该新方法是在比例边界有限元框架中用NURBS曲线或曲面精确描述域边界几何形状,同时域边界位移场采用描述几何形状的NURBS形函数等参构造。这种新方法具有比例边界有限元固有的径向解析特性和NURBS的高阶连续性的优点。数值算例显示,与传统的比例边界有限元相比,基于等几何分析的比例边界有限元方法提高了域边界单元和域内应力场的连续性,减少了计算自由度。应用此方法可以用较少的计算自由度获得更高连续阶和更高精度的位移、应力和应变场。     

Some applications of boundary element methods, boundary element-finite element coupling techniques in elasto-plastic stress analysis

In this paper some research works of author's group on boundary element method for linear elastic and elasto-plastic problems are briefly summarized. The BE subrogions method and BE-FE coupling techniques are introduced, and some idea of the error estimation of BE solutions is presented.This work is supported by the National Natural Science Foundation of China.     

APPLICATION OF TREFFTZ BEM TO ANTI-PLANE PIEZOELECTRIC PROBLEM

Anti-plane electroelastic problems are studied by the Trefftz boundary elementmethod (BEM) in this paper. The Trefftz BEM is based on a weighted residual formulation andindirect boundary approach. In particular the point-collocation and Galerkin techniques, in whichthe basic unknowns are the retained expansion coefficients in the system of complete equations,are considered, Furthermore, special Trefftz functions and auxiliary functions which satisfy ex-actly the specified boundary conditions along the slit boundaries are also used to derive a specialpurpose element with local defects. The path-independent integral is evaluated at the tip of acrack to determine the energy release rate for a mode Ⅲ fracture problem. In final, the accuracyand efficiency of the Trefftz boundary element method are illustrated by an example and thecomparison is made with other methods.     

基于正交基无单元Galerkin法和非线性规划的安定分析方法

基于安定分析的下限定理,用正交基无单元Galerkin法建立了交交载荷作用下理想弹塑性结构安定分析的下限计算格式.在给定载荷域的载荷角点所对应的载荷作用下,采用正交基无单元Galerkin法计算相应的虚拟弹性应力场.并且利用结构在正交基无单元Galerkin法弹塑性增量分析中平衡迭代结果计算得到自平衡应力场基矢量,然后由这些基矢量的线性组合模拟自平街应力场.安定分析问题最终被归结为一系列未知变量较少的非线性数学规划子问题,通过复合形法求解.算例表明该方法的计算结果是令人满意的,并且对初始复合形顶点和用于构造自平衡应力场基矢量的载荷增量是非常不敏感的.     

Discontinuous Galerkin spectral element lattice Boltzmann method on triangular element

Discontinuous Galerkin spectral element method is used to solve the lattice Boltzmann equation (LBE) in the discrete velocity space. The triangular elements are adopted because of their flexibility to deal with complex geometries. The flow past a circular cylinder is simulated by the proposed scheme. The results are consistent with those obtained from the previous numerical methods and experiments. Copyright © 2003 John Wiley & Sons, Ltd.