首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
    
A model which allows the introduction of displacements jumps to conventional finite elements is developed. The path of the discontinuity is completely independent of the mesh structure. Unlike so‐called ‘embedded discontinuity’ models, which are based on incompatible strain modes, there is no restriction on the type of underlying solid finite element that can be used and displacement jumps are continuous across element boundaries. Using finite element shape functions as partitions of unity, the displacement jump across a crack is represented by extra degrees of freedom at existing nodes. To model fracture in quasi‐brittle heterogeneous materials, a cohesive crack model is used. Numerical simulations illustrate the ability of the method to objectively simulate fracture with unstructured meshes. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

2.
    
Discontinuous failure is simulated via the introduction of a geometrical discontinuity. The cohesive zone is modelled via the use of the partition‐of‐unity property of the finite element interpolation. By this approach, a crack can pass through the elements without any restriction to the underlying mesh. Despite such a feature, it has been confirmed that a sufficiently fine mesh discretization still needs to be ensured in order to obtain a correct crack path and mechanical response. The p‐adaptive scheme, driven by error in an energy norm measure or a goal‐oriented measure, has been examined due to its implementational simplicity. The results have shown that, if considering only increasing the polynomial degree, the p‐approach can greatly improve the results. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

3.
    
A recent approach to fracture modeling has combined the extended finite element method (XFEM) with cohesive zone models. Most studies have used simplified enrichment functions to represent the strong discontinuity but have lacked an analytical basis to represent the displacement gradients in the vicinity of the cohesive crack. In this study enrichment functions based upon an existing analytical investigation of the cohesive crack problem are proposed. These functions have the potential of representing displacement gradients in the vicinity of the cohesive crack and allow the crack to incrementally advance across each element. Key aspects of the corresponding numerical formulation and enrichment functions are discussed. A parameter study for a simple mode I model problem is presented to evaluate if quasi‐static crack propagation can be accurately followed with the proposed formulation. The effects of mesh refinement and mesh orientation are considered. Propagation of the cohesive zone tip and crack tip, time variation of the cohesive zone length, and crack profiles are examined. The analysis results indicate that the analytically based enrichment functions can accurately track the cohesive crack propagation of a mode I crack independent of mesh orientation. A mixed mode example further demonstrates the potential of the formulation. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

4.
    
A new two‐level multiscale enrichment methodology for analysis of heterogeneous plates is presented. The enrichments are applied in the displacement and strain levels: the displacement field of a Reissner–Mindlin plate is enriched using the multiscale enrichment functions based on the partition of unity principle; the strain field is enriched using the mathematical homogenization theory. The proposed methodology is implemented for linear and failure analysis of brittle heterogeneous plates. The eigendeformation‐based model reduction approach is employed to efficiently evaluate the non‐linear processes in case of failure. The capabilities of the proposed methodology are verified against direct three‐dimensional finite element models with full resolution of the microstructure. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

5.
    
In the first part of this contribution, a brief theoretical revision of the mechanical and variational foundations of a Failure‐Oriented Multiscale Formulation devised for modeling failure in heterogeneous materials is described. The proposed model considers two well separated physical length scales, namely: (i) the macroscale where nucleation and evolution of a cohesive surface is considered as a medium to characterize the degradation phenomenon occurring at the lower length scale, and (ii) the microscale where some mechanical processes that lead to the material failure are taking place, such as strain localization, damage, shear band formation, and so on. These processes are modeled using the concept of Representative Volume Element (RVE). On the macroscale, the traction separation response, characterizing the mechanical behavior of the cohesive interface, is a result of the failure processes simulated in the microscale. The traction separation response is obtained by a particular homogenization technique applied on specific RVE sub‐domains. Standard, as well as, Non‐Standard boundary conditions are consistently derived in order to preserve objectivity of the homogenized response with respect to the micro‐cell size. In the second part of the paper, and as an original contribution, the detailed numerical implementation of the two‐scale model based on the finite element method is presented. Special attention is devoted to the topics, which are distinctive of the Failure‐Oriented Multiscale Formulation, such as: (i) the finite element technologies adopted in each scale along with their corresponding algorithmic expressions, (ii) the generalized treatment given to the kinematical boundary conditions in the RVE, and (iii) how these kinematical restrictions affect the capturing of macroscopic material instability modes and the posterior evolution of failure at the RVE level. Finally, a set of numerical simulations is performed in order to show the potentialities of the proposed methodology, as well as, to compare and validate the numerical solutions furnished by the two‐scale model with respect to a direct numerical simulation approach. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

6.
    
The performance of partition‐of‐unity based methods such as the generalized finite element method or the extended finite element method is studied for the simulation of cohesive cracking. The focus of investigation is on the performance of bilinear quadrilateral finite elements using these methods. In particular, the approximation of the displacement jump field, representing cohesive cracks, by extended finite element method/generalized finite element method and its effect on the overall behavior at element and structural level is investigated. A single element test is performed with two different integration schemes, namely the Newton‐Cotes/Lobatto and the Gauss integration schemes, for the cracked interface contribution. It was found that cohesive crack segments subjected to a nonuniform opening in unstructured meshes (or an inclined crack in a structured finite element mesh) result in an unrealistic crack opening. The reasons for such behavior and its effect on the response at element level are discussed. Furthermore, a mesh refinement study is performed to analyze the overall response of a cohesively cracked body in a finite element analysis. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

7.
8.
    
Present extended finite element method (XFEM) elements for cohesive crack growth may often not be able to model equal stresses on both sides of the discontinuity when acting as a crack‐tip element. The authors have developed a new partly cracked XFEM element for cohesive crack growth with extra enrichments to the cracked elements. The extra enrichments are element side local and were developed by superposition of the standard nodal shape functions for the element and standard nodal shape functions for a sub‐triangle of the cracked element. With the extra enrichments, the crack‐tip element becomes capable of modelling variations in the discontinuous displacement field on both sides of the crack and hence also capable of modelling the case where equal stresses are present on each side of the crack. The enrichment was implemented for the 3‐node constant strain triangle (CST) and a standard algorithm was used to solve the non‐linear equations. The performance of the element is illustrated by modelling fracture mechanical benchmark tests. Investigations were carried out on the performance of the element for different crack lengths within one element. The results are compared with previously obtained XFEM results applying fully cracked XFEM elements, with computational results achieved using standard cohesive interface elements in a commercial code, and with experimental results. The suggested element performed well in the tests. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

9.
Mesh-independent discrete numerical representations of cohesive-zone models   总被引:2,自引:0,他引:2  
The importance of the cohesive-zone approach to analyse localisation and fracture in engineering materials is emphasised and ways to incorporate the cohesive-zone methodology in computational methods are discussed. Until recently, numerical implementations of cohesive-zone models have suffered from a certain mesh bias. For discrete representations this is caused by the initial mesh design, while for smeared representations it is rooted in the ill-posedness of the rate boundary value problem that arises upon the introduction of decohesion. A proper representation of the discrete character of cohesive-zone formulations which avoids any mesh bias is obtained by exploiting the partition-of-unity property of finite element shape functions. The recently developed cohesive segments method, which is well-suited for simulating the entire process of crack nucleation, growth and coalescence is reviewed. The effectiveness of this approach is illustrated by some examples for crack nucleation and growth in heterogeneous materials and for fast crack growth.  相似文献   

10.
    
A novel numerical method is proposed for modelling time‐harmonic acoustic propagation of short wavelength disturbances on non‐uniform potential flows. The method is based on the partition of unity finite element method in which a local basis of discrete plane waves is used to enrich the conventional finite element approximation space. The basis functions are local solutions of the governing equations. They are able to represent accurately the highly oscillatory behaviour of the solution within each element while taking into account the convective effect of the flow and the spatial variation in local sound speed when the flow is non‐uniform. Many wavelengths can be included within a single element leading to ultra‐sparse meshes. Results presented in this article will demonstrate that accurate solutions can be obtained in this way for a greatly reduced number of degrees of freedom when compared to conventional element or grid‐based schemes. Numerical results for lined uniform two‐dimensional ducts and for non‐uniform axisymmetric ducts are presented to indicate the accuracy and performance which can be achieved. Numerical studies indicate that the ‘pollution’ effect associated with cumulative dispersion error in conventional finite element schemes is largely eliminated. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

11.
    
In this paper, we propose a novel image interpolation method by usingGaussian-Sinc automatic interpolators with partition of unity property. A comprehensivecomparison is made with classical image interpolation methods, such as the bicubicinterpolation, Lanczos interpolation, cubic Schaum interpolation, cubic B-splineinterpolation and cubic Moms interpolation. The experimental results show theeffectiveness of the improved image interpolation method via some image quality metricssuch as PSNR and SSIM.  相似文献   

12.
This paper presents a new computational approach dedicated to the fracture of nonlinear heterogeneous materials. This approach extends the standard periodic homogenization problem to a two field cohesive-volumetric finite element scheme. This two field finite element formulation is written as a generalization Non-Smooth Contact Dynamics framework involving Frictional Cohesive Zone Models. The associated numerical platform allows to simulate, at finite strain, the fracture of nonlinear composites from crack initiation to post-fracture behavior. The ability of this computational approach is illustrated by the fracture of the hydrided Zircaloy under transient loading.  相似文献   

13.
The paper focuses on the Continuum Strong Discontinuity Approach (CSDA) to fracture mechanics, and the traction-separation cohesive laws induced from continuum dissipative models as their projections onto the failure interface. They are compared with the cohesive laws commonly used for the fracture simulation in quasi-brittle materials, typically concrete. Emphasis is placed in the analysis of the mechanical stress-strain states induced by the CSDA into the fracture process zone: first when the damage mechanism is initiated and, after, when the cohesive model determines the crack response. The influence of the material parameters, particularly the fracture energy and the initial continuum softening modulus, in the obtained phenomenological responses is also analyzed. Representative numerical solutions of fracture problems are finally presented.  相似文献   

14.
    
The numerical analysis of large numbers of arbitrarily distributed discrete thin fibres embedded in a continuum is a computationally demanding process. In this contribution, we propose an approach based on the partition of unity property of finite element shape functions that can handle discrete thin fibres in a continuum matrix without meshing them. This is made possible by a special enrichment function that represents the action of each individual fibre on the matrix. Our approach allows to model fibre‐reinforced materials considering matrix, fibres and interfaces between matrix and fibres individually, each with its own elastic constitutive law. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

15.
    
A new approach for the computation of the forced vibrations up to the medium frequency range is formulated for thin plates. It is based on the partition of unity method (PUM), first proposed by Babu?ka, and used here to solve the elastodynamic problem. The paper focuses on the introduction of local information in the basis of the PUM coming from previous approximations, in order to enhance the accuracy of the solution. The method may be iterative and generates a PUM approximation leading to smaller models compared with the finite element ones required for a same accuracy level. It shows very promising results, in terms of frequency range, accuracy and computational time. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

16.
    
Enhanced transformation field analysis (E-TFA), recently proposed for reduced-order modeling, is here formulated for and applied to multiscale analysis. The approach is able to reproduce a highly complex nonlinear macroscale behavior, resulting from a microstructure with cohesive interfaces embedded in an elasto-plastic bulk. E-TFA features a consistent tangent matrix in its solution procedure, which enables a straightforward definition of the upscaled tangent stiffness tensor. Numerical tests show that, compared to FE 2$$ {}^2 $$, the proposed approach yields accurate solutions at a lower computational cost.  相似文献   

17.
沈珉  郝培 《复合材料学报》2016,33(1):189-197
采用内聚力模型(CZM)描述颗粒增强复合材料(PRCs)颗粒/基体非理想界面的力学行为,采用细观力学的Mori-Tanaka(M-T)方法和稀释解方法研究非理想界面刚度对该复合材料有效模量的影响。结果表明:复合材料某一体积含量下的有效模量和其颗粒/基体非理想界面刚度之间存在单调递增的关系曲线。同一复合材料不同体积含量的有效模量和非理想界面刚度关系曲线存在一个临界交点(CP),该点对应的临界界面刚度控制着颗粒体积分数对有效模量的影响。研究了基体和增强相的力学性能以及增强颗粒的尺寸对CP点临界界面刚度的影响。采用有效模量-界面刚度关系曲线,结合实验测得的有效模量,提出了估算PRCs的非理想界面刚度的方法,进而估计了复合材料的宏观有效模量。   相似文献   

18.
    
In this paper, the Polytope Finite Element Method is employed to model an embedded interface through the body, independent of the background FEM mesh. The elements that are crossed by the embedded interface are decomposed into new polytope elements which have some nodes on the interface line. The interface introduces discontinuity into the primary variable (strong) or into its derivatives (weak). Both strong and weak discontinuities are studied by the proposed method through different numerical examples including fracture problems with traction‐free and cohesive cracks, and heat conduction problems with Dirichlet and Dirichlet–Neumann types of boundary conditions on the embedded interface. For traction‐free cracks which have tip singularity, the nodes near the crack tip are enriched with the singular functions through the eXtended Finite Element Method. The concept of Natural Element Coordinates (NECs) is invoked to drive shape functions for the produced polytopes. A simple treatment is proposed for concave polytopes produced by a kinked interface and also for locating crack tip inside an element prior to using the singularity enrichment. The proposed method pursues some implementational details of eXtended/Generalized Finite Element Methods for interfaces. But here the additional DOFs are constructed on the interface lines in contrast to X/G‐FEM, which attach enriched DOFs to the previously existed nodes. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

19.
    
We present a function space‐oriented coupling approach for the multiscale simulation of non‐linear processes in mechanics using finite elements and molecular dynamics concurrently. The key idea is to construct a transfer operator between the different scales on the basis of weighted local averaging instead of using point wise taken values. The local weight functions are constructed by assigning a partition of unity to the molecular degrees of freedom (Shepard's approach). This allows for decomposing the micro scale displacements into a low‐frequency and a high‐frequency part by means of a weighted L2‐projection. Numerical experiments illustrating the stabilizing effect of our coupling approach are given. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

20.
    
An enriched partition of unity FEM is developed to solve time‐dependent diffusion problems. In the present formulation, multiple exponential functions describing the spatial and temporal diffusion decay are embedded in the finite element approximation space. The resulting enrichment is in the form of a local asymptotic expansion. Unlike previous works in this area where the enrichment must be updated at each time step, here, the temporal decay in the solution is embedded in the asymptotic expansion. Thus, the system matrix that is evaluated at the first time step may be decomposed and retained for every next time step by just updating the right‐hand side of the linear system of equations. The advantage is a significant saving in the computational effort where, previously, the linear system must be reevaluated and resolved at every time step. In comparison with the traditional finite element analysis with p‐version refinements, the present approach is much simpler, more efficient, and yields more accurate solutions for a prescribed number of DoFs. Numerical results are presented for a transient diffusion equation with known analytical solution. The performance of the method is analyzed on two applications: the transient heat equation with a single source and the transient heat equation with multiple sources. The aim of such a method compared with the classical FEM is to solve time‐dependent diffusion applications efficiently and with an appropriate level of accuracy. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号