Stability of incompressible formulations enriched with X-FEM
|
|
- Deborah Miller
- 5 years ago
- Views:
Transcription
1 Legrain, N. Moes and A. Huerta,Stability of incompressible formulations enriched with X-FEM, G., Computer Methods in Applied Mechanics and Engineering, Vol. 197, Issues 21-24, pp , 2008 Stability of incompressible formulations enriched with X-FEM G. Legrain a N. Moës a, A. Huerta b a GeM Institute - École Centrale de Nantes / Université de Nantes / CNRS 1 Rue de la Noë, Nantes, France. b Laboratori de Càlcul Numèric (LaCàN) - Edifici C2, Campus Nord, Universitat Politècnica de Catalunya - E Barcelona, Spain. Abstract The treatment of (near-)incompressibility is a major concern for applications involving rubber-like materials, or when important plastic flows occurs as in forming processes. The use of mixed finite element methods is known to prevent the locking of the finite element approximation in the incompressible limit. However, it also introduces a critical condition for the stability of the formulation, called the infsup or LBB condition. Recently, the finite element method has evolved with the introduction of the partition of unity. The extended Finite Element Method (X- FEM) uses the partition of unity to remove the need to mesh physical surfaces or to remesh them as they evolve. The enrichment of the displacement field makes it possible to treat surfaces of discontinuity inside finite elements. In this paper, some strategies are proposed for the enrichment of mixed finite element approximations in the incompressible setting. The case of holes, material interfaces and cracks are considered. Numerical examples show that for well chosen enrichment strategies, the finite element convergence rate is preserved and the inf-sup condition is passed. Key words: Mixed formulation, X-FEM, Partition of Unity, inf-sup condition, incompressibility, holes, inclusions, fracture mechanics Corresponding author addresses: gregory.legrain@ec-nantes.fr (G. Legrain), nicolas.moes@ec-nantes.fr (N. Moës), antonio.huerta@upc.edu (A. Huerta). Preprint submitted to Elsevier Science 13 September 2005
2 1 Introduction Displacement based finite element methods are nowadays abundantly used in engineer analysis. Indeed, they can solve a wide variety of problems, and have now been deeply mathematically investigated. However there still exists two main drawbacks for these methods. First, the treatment of incompressible or nearly incompressible problems necessitates the use of adapted formulations. If not, incompressibility constraint locks the approximation, leading for instance to non-physical displacement fields. Second, the generation and especially the update of the mesh in complex 3D settings for evolving boundaries such as cracks, material interfaces and voids still lacks robustness, and involves important human effort. Several techniques have been developed to respond to the locking issue. For instance, the selective-reduced-integration procedures [1 3] or the Bbar approach of Hugues [4] in which the volumetric part of the strain tensor is evaluated at the center of the element. Another way to avoid locking is to enhance the strain tensor in order to enlarge the space on which the minimization is performed, and meet the divergence-free condition (enhanced assumed strain methods, see [5 9]). Here, we will focus on two-field mixed finite element methods. The incompressibility constraint is weakened by the introduction of the pressure field. This alleviates locking at the price of additional pressure unknowns. However, mixed finite element methods are not stable in all cases, some of them showing spurious pressure oscillations if displacement and pressure spaces are not chosen carefully. To be stable, a mixed formulation must verify consistency, ellipticity and the so called inf-sup (or LBB) condition. The later is a severe condition which depends on the connection between the displacement and pressure approximation spaces. Stable mixed formulations can be obtained by stabilizing non-stable formulations with the use of parameters whose values may depend on the problem at hand. Otherwise, one has to work with approximation spaces which passe the inf-sup condition. To prove that a displacement-pressure pair satisfies the inf-sup condition is not a trivial task. However, a numerical test has been proposed by Brezzi and Fortin [10] then by Chapelle et al. [11], in order to draw a prediction on the fulfilment of the inf-sup. This test proved to be useful for the study of the stability of various mixed elements [12]. The second drawback of classical finite element methods (evolving boundaries) has been overcome by the development of alternative methods such as meshless methods in which the connectivity between the nodes is no longer obtained by the mesh, but by domains of influence which can be split by the boundaries. Moreover, the approximation basis can be enriched with functions coming from the physical knowledge of the problem. Note that in the incompressible limit, meshless methods are now known to lock [13,14] as classical finite elements. Thus, some strategies have been developed to circumvent this 2
3 issue [13,15]. Apart from mixing meshless methods and finite elements [16] another alternative to overcome the re-meshing issue in finite elements is to use the extended Finite Element Method (X-FEM) based on the partition of unity framework introduced by Babuška and Melenk [17]. Proper enrichment of the finite element basis makes it possible to model crack, material inclusions and holes with non-conforming meshes. The X-FEM method has been used for the simulation of a wide variety of problems such as fracture mechanics problems (2D[18 20], 3D[21 23], plates [24,25], cohesive zone modeling [26,27], dynamic fracture [28], nonlinear fracture mechanics[29 31]), holes [32,33], but also material inclusions [33,34] or multiple phase flows [35]. Here, we focus on the application of this method to mixed formulations for the treatment of holes, material inclusions and cracks in the incompressible limit. Bbar or selective-reduced formulations are not considered, because they do not seem to be generalized easily to enriched displacement fields. The main contribution of this paper is the design of enrichment strategies for the pressure and displacement fields, so that it leads to a stable formulation. The enrichment of mixed finite element approximations has already been used by Dolbow et al. [25] and Areias et al.[31] for fracture mechanics in plates and shells, and by Wagner et al. [36] for rigid particles in Stokes flow. However, the stability and the convergence of these approaches was not studied. The latest work concerning volumetric incompressibility was proposed by Dolbow and Devan [29]. In this papers, the authors focus on the application of the enhanced assumed strain method to X-FEM in large strain. This approach seems to lead to a stable low order formulation in the case of nearly incompressible nonlinear fracture mechanics. However, the stability of the method was not shown, and the influence of the near-tip enrichment was not studied. More precisely, it is not clear whether the near-tip enrichment could make this approach unstable, as the construction of an orthogonal enhanced strain field becomes difficult with non polynomial functions. The paper is organized as follows: first, the governing equations of incompressible linear elasticity are recalled. The conditions for the stability of mixed formulations are also reviewed. Next, some strategies are proposed to keep the stability of enriched finite elements. The case of holes, material interfaces and 2D cracks are presented. Finally, in a last section the stability of these strategies is investigated. 2 Governing Equations In this section, we focus on the design of stable mixed formulations for the treatment of incompressible elasticity. First, the equations governing incompressible linear elasticity are recalled. Then the inf-sup condition is presented together with a numerical test. 3
4
5 and ε D is the deviatoric strain operator: ε D = ε ε V 3 I (5) When the material tends to incompressibility, the bulk modulus tends to infinity. This means that ε V must tend to zero in order to meet condition (2) (the displacement field must be divergence free in the incompressible limit). ε V = div(u) 0 as ν 0.5 (6) The strong form (1) is equivalent to the stationarity of a displacement potential Π: Π(u) = 1 ε(u) : C : ε(u) dω u b dω u T d dγ (7) 2 Ω Ω Ω t In order to model incompressible or almost incompressible problems, a two field principle is considered by introducing a second variable (the hydrostatic pressure p) in the potential (7): p = κε V (u) = 1 T r(σ) (8) 3 When κ increases, the volumetric strain ε V decreases and becomes very small. For total incompressibility, the bulk modulus is infinite, the volumetric strain is zero, and the pressure remains finite (of the order of the applied boundary tractions). The stress tensor is then expressed as: σ = p I + 2µ ε D in Ω (9) The solution of the governing differential equations (1) now involves two variables: the displacement field and the pressure field. Writing the two fields variational principle, the total potential for the u p formulation is expressed as: χ(u, p) = 1 2 Ω ε D (u) : C : ε D (u) dω 1 2 u b dω u T d dγ (10) Ω t p 2 Ω κ Ω dω p ε V (u) dω Invoking the stationarity of χ(u, p) with respect to the two independent variables u and p, we obtain: Ω δε D : C : ε D dω p δε V dω = R(δv) (11) Ω ( ) p κ + ε V δp dω = 0 (12) Ω where R(δv) represents the virtual work of the external loads. In the case of Ω 5
6
7 expresses as: The existence of a stable finite element approximate solution (u h, p h ) depends on choosing a pair of spaces V h and Q h such that the following condition holds: inf sup q h Q h v h V h Ω qh div v h dω v h 1 q h 0 β > 0 (17) where 1 and 0 indicates H 1 and L 2 norms respectively and β is independent of the mesh size h. If the inf-sup compatibility condition is satisfied, then there exists a unique u h V h and a p h Q h (determined up to an arbitrary constant in the case of purely Dirichlet boundary conditions) Numerical assessment of the inf-sup condition As seen before, the prediction of the stability of a mixed formulation involves the fulfillment of the inf-sup criterion. This criterion is however impossible to prove for practical situations. This is why the numerical evaluation of the inf-sup condition has received considerable attention [11,10]. This numerical evaluation, although not equivalent to the analytical inf-sup, gives indications on whether (17) is fulfilled or not for a given set of finite element discretizations. The numerical inf-sup test is based on the following theorem. Proposition 1 Let M uu and M pp be the mass matrices associated to the scalar products of V h and Q h respectively and let µ min be the smallest non zero eigenvalue defined by the following eigenproblem: then the value of β is simply µ min. K up T M uu 1 K up q = µ 2 M pp q (18) The proof can be found in [37] or [10]. The numerical test proposed in [11] consists in testing a particular formulation by calculating β using meshes of increasing refinement. On the basis of three or four results it can be predicted whether the inf-sup value is probably bounded from underneath or, on the contrary, goes down to zero when the mesh is refined. The reliability of this test is demonstrated on several examples of elements for incompressible elasticity problems in [11]. In the following section this test is used to check the behavior of proposed enrichment strategies. However, we follow [11] and use only S uu = Ω u : u dω instead of M uu in (18). In order to perform the numerical infsup test, a sequence of successive refined meshes is considered. The objective is to monitor the inf-sup values, β, when h decreases. If a steady decrease in log(β) is observed when h goes to zero, the element is predicted to violate the inf-sup condition and said to fail the numerical test. But, if the log(β) value is stable as the number of elements increases, the test is numerically passed. 7
8 3 X-FEM Discretization 3.1 Displacement field With classical finite elements, the approximation of a vector field u on an element Ω e is written as: u(x) Ωe = n u α=1 u α N α u(x) (19) where n u is the number of coefficients describing the approximation of the displacement over the element, u α is the α th coefficient of this approximation and N α u is the vectorial shape function associated to the coefficient u α. Within the partition of unity, the approximation is enriched as: u(x) Ωe = n u N α u α=1 u α + n enr β=1 a α β φ u β(x) (20) where n enr is the number of enrichment modes, a α β is the additional dof associated to dof α and φ u β stands for the β th scalar enrichment function. The number and the expression of the enrichment functions vary with the problem to model. The expression of this enrichment function will be recalled for holes, inclusions and fracture mechanics in the next sections. 3.2 Pressure field Using the same scheme, the pressure approximation is written as: p(x) Ωe = n p α=1 N α p (pα + a α φ p (x)) (21) where n p is the number of coefficients describing the approximation over the element, p α is the α th coefficient of this approximation and Np α is the scalar shape function associated to the coefficient p α, a α is the additional dof associated to dof α and φ p stands for the scalar pressure enrichment function. The key issue is the combined choice of enrichment functions φ u and φ p such that the whole enriched approximation (displacement and pressure) passes the inf-sup condition. 8
9
10
11
12
13
14
15
16
17
18
19
20 function are the most effective since they pass the inf-sup test. On the contrary, Heaviside based strategies do not pass this test and should not be considered. The formulations where only the linear part of the displacement is enriched (N 4 and 8) with the ridge are interesting, since less degrees of freedom are involved in the approximation. This should be important in the context of an extension to three dimensional studies. 4.3 Incompressible fracture mechanics The resolution of compressible fracture mechanics problems has been extensively studied in the context of the X-FEM for both 2D [39,18,40,32,24] and 3D fracture mechanics [21,22]. The most common enrichment strategy consists in using the asymptotic displacement field as an enrichment for the displacement finite element approximation. In the context of incompressible media, the analytical asymptotic displacement field (Westergaard solution) is shown to be identical to the limit of the compressible one. The asymptotic evolution of the pressure field can be obtained also using the Westergaard solution. p(r, θ) = 2K I 3 2 π r cos ( ) θ + 2K II π r sin ( ) θ 2 (25) φ u = { ( ) r θ sin, r cos 2 ( ) θ, r sin 2 r cos ( θ 2 ( ) θ sin(θ), 2 ) } sin(θ) (26) We use these expressions as an enrichment for the pressure field in the near-tip region. Thus, the enrichment basis for the pressure is expressed as: φ p I = 1 r cos ( ) θ 2 ) (27) φ p II = 1 r sin ( θ 2 Note that a classical Heaviside enrichment is considered for both pressure and displacement for nodes whose support is fully cut by the crack, and that only the Mini element is considered hereunder Convergence study Consider a domain Ω = [ 1, 1] [ 1, 1] under tension (see Figure 23). The tensions applied on the boundary of the domain are related to the exact ten- 20
21
22
23
24 lation. Moreover, the use of the geometrical enrichment leads to an improved convergence case, similar to the compressible rate. Finally, degrees of freedom can be saved by only enriching the linear part of the approximation. 5 Conclusion Some strategies for enriching existing mixed finite element methods have been presented. These strategies are natural extensions of the displacement-based X-FEM, and are shown to preserve the classical finite element convergence rate. The stability of these strategies has been shown through the numerical inf-sup test. However, quadratic-based elements could not be tested completely, as in some cases the linear interpolation of the level-set leads to a degraded rate of convergence. The construction and the validation of isoparametric quadratic elements will be the subject of a forthcoming paper. The method should also be applied to finite strain mechanics, as the fulfilment of the inf-sup condition seems to be a prerequisite to build efficient large strain formulations [42]. 24
25 References [1] D.J. Naylor. Stress in nearly incompressible materials for finite elements with application to the calculation of excess pore pressure. international journal for numerical methods in engineering, 8: , [2] T.J.R. Hughes, R.L. Taylor, and J.F. Levy. High Reynolds number, steady, incompressible flows by a finite element method. John Wiley & Sons, [3] S.F. Pawsey and R.W. Clough. Improved numerical integration of thick slab finite elements. International journal for numerical methods in engineering, 3: , [4] TJR Hughes. Generalization of selective reduced integration procedures to anisotropic and nonlinear media. International Journal for Numerical Methods in Engineering, 15: , [5] RL. Taylor, PJ. Beresford, and EL. Wilson. A nonconforming element for stress analysis. International Journal for Numerical Methods in Engineering, 10(6): , [6] EL. Wilson, RL. Taylor, WP. Doherty, and J. Ghaboussi. Incompatible displacement models In Numerical and Computer Models in Structural Mechanics [7] JC. Simo and MS. Rifai. A class of mixed assumed strain methods and the method of incompatible modes. International Journal for Numerical Methods in Engineering, 29: , [8] JC. Simo and F. Armero. Geometrically non-linear enhanced strain mixed methods and the method of incompatible modes. International Journal for Numerical Methods in Engineering, 33: , [9] B.D. Reddy and J.C. Simo. stability and convergence of a class of enhanced strain methods. SIAM J. Numer. Anal., 32(6): , [10] F. Brezzi and M. Fortin. Mixed and Hybrid Finite Element Methods. Springer, New York, [11] D. Chapelle and K.J. Bathe. The inf-sup test. Computers & structures, 47(4-5): , [12] K.J. Bathe. The inf-sup condition and its evaluation for mixed finite element methods. Computers & Structures, 79: , [13] A. Huerta and S. Fernández-Méndez. Locking in the incompressible limit for the element free galerkin method. international journal for numerical methods in engineering, 50, [14] J. Dolbow and T. Belytschko. Volumetric locking in the element free galerkin method. Int. J. Numer. Methods Eng., 46(6): ,
26 [15] A. Huerta, Y. Vidal, and P. Villon. Pseudo-divergence-free element free galerkin method for incompressible fluid flow. Computer Methods in Applied Mechanics and Engineering, 193(12-14): , [16] A. Huerta and S. Fernández-Méndez. Enrichment and coupling of the finite element and meshless methods. International Journal for Numerical Methods in Engineering, 48(11): , [17] J.M. Melenk and I. Babuška. The partition of unity finite element method : Basic theory and applications. Comput. Methods Appl. Mech. Engrg., 139: , [18] N. Moës, J. Dolbow, and T. Belytschko. A finite element method for crack growth without remeshing. International Journal for Numerical Methods in Engineering, 46: , [19] N. Sukumar, Z.Y. Huang, J.H. Prévost, and Z. Suo. Partition of unity enrichment for bimaterial interfacial cracks. International Journal for Numerical Methods in Engineering, 59: , [20] E. Budyn, G. Zi, N. Moës, and T. Belytschko. A model for multiple crack growth in brittle materials without remeshing. International journal for numerical methods in engineering, [21] N. Moës, A. Gravouil, and T. Belytschko. Non-planar 3D crack growth by the extended finite element and level sets. part I: Mechanical model. International Journal for Numerical Methods in Engineering, 53: , [22] A. Gravouil, N. Moës, and T. Belytschko. Non-planar 3d crack growth by the extended finite element and level sets. part II: level set update. International Journal for Numerical Methods in Engineering, 53: , [23] N. Sukumar, N. Moës, T. Belytschko, and B. Moran. Extended Finite Element Method for three-dimensional crack modelling. International Journal for Numerical Methods in Engineering, 48(11): , [24] J. Dolbow, N. Moës, and T. Belytschko. An extended finite element method for modeling crack growth with frictional contact. Comp. Meth. in Applied Mech. and Engrg., 190: , [25] J. Dolbow, N. Moës, and T. Belytschko. Modeling fracture in Mindlin-Reissner plates with the extended finite element method. Int. J. Solids Structures, 37: , [26] N. Moës and T. Belytschko. Extended finite element method for cohesive crack growth. Engineering Fracture Mechanics, 69: , [27] G.N. Wells and L.J. Sluys. A new method for modelling cohesive cracks using finite elements. International Journal for Numerical Methods in Engineering, 50: , [28] J. Réthoré, A. Gravouil, and A. Combescure. An energy-conserving scheme for dynamic crack growth using the extended finite element method. International Journal for Numerical Methods in Engineering,
27 [29] J.E. Dolbow and A. Devan. Enrichment of enhanced assumed strain approximations for representing strong discontinuities : Addressing volumetric incompressibility and the discontinuous patch test. International journal for numerical methods in engineering, 59:47 67, [30] G. Legrain, N. Moës, and E. Verron. Stress analysis around crack tips in finite strain problems using the extended finite element method. International Journal for Numerical Methods in Engineering, 63(2): , [31] Pedro M. A. Areias and Ted Belytschko. Non-linear analysis of shells with arbitrary evolving cracks using xfem. International Journal for Numerical Methods in Engineering, 62(3): , [32] C. Daux, N. Moës, J. Dolbow, N. Sukumar, and T. Belytschko. Arbitrary branched and intersecting cracks with the extended Finite Element Method. International Journal for Numerical Methods in Engineering, 48: , [33] N. Sukumar, D. L. Chopp, N. Moës, and T. Belytschko. Modeling holes and inclusions by level sets in the extended finite element method. Comp. Meth. in Applied Mech. and Engrg., 190: , [34] N. Moës, M. Cloirec, P. Cartraud, and J.-F. Remacle. A computational approach to handle complex microstructure geometries. Comp. Meth. in Applied Mech. and Engrg., 192: , [35] J. Chessa and T. Belytschko. An extended finite element method for two-phase fluids. Journal of Applied Mechanics (ASME), 70:10 17, [36] G.J. Wagner, N. Moës, W.K. Liu, and T. Belytschko. The extended finite element method for Stokes flow past rigid cylinders. International Journal for Numerical Methods in Engineering, 51: , [37] D. S. Malkus. Eigenproblems associated with the discrete LBB condition for incompressible finite elements. Int. J. Eng. Sci., 19(10): , [38] D.N. Arnold, F. Brezzi, and M. Fortin. Stable finite element for stokes equations. Calcolo, 21: , [39] T. Belytschko and T. Black. Elastic crack growth in finite elements with minimal remeshing. International Journal for Numerical Methods in Engineering, 45(5): , [40] D.L. Chopp and N. Sukumar. Fatigue crack propagation of multiple coplanar cracks with the coupled extended finite element/fast marching method. International journal of engineering science, 41: , [41] E. Bechet, H. Minnebo, N. Moës, and B. Burgardt. Improved implementation and robustness study of the x-fem method for stress analysis around cracks. International Journal for Numerical Methods in Engineering, Accepted, [42] D. Pantuso and K.J. Bathe. On the stability of mixed finite elements in large strain analysis of incompressible solids. Finite elements in analysis and design, 28:83 104,
28 A Analytical solution for curved interface We did construct a specific analytical solution to investigate the convergence of the X-FEM. It represents two perfectly bounded rings made from different materials (see figure 7). The loading is chosen such that: b 1r (r) = 3 W 3 c W 2 c + W 1 b 2r (r) = 3 V 3 c V 2 c + V 1 T r = W 3 c 3 + W 2 c 2 + W 1 c + W 0 T θ = 2 µ 1µ 2 a 2 b 2 ( ) a 2 µ 1 c 2 + a 2 µ 2 c 2 µ 2 b 2 c 2 + a 2 b 2 µ 1 c Under these boundary conditions, the stress field components are: σ 1 rr = ( W 3 r 3 + W 2 r 2 + W 1 r + W 0 ) σ 1 θθ = σ1 rr σ 1 rθ = 2µ 1 µ 2 c a 2 b 2 ( a 2 µ 1 c 2 + a 2 µ 2 c 2 µ 2 b 2 c 2 + a 2 b 2 µ 1 ) r 2 (A.1) (A.2) (A.3) (A.4) (A.5) (A.6) (A.7) σ 2 rr = ( V 3 r 3 + V 2 r 2 + V 1 r + V 0 ) σθθ 2 = σ2 rr σrθ 2 2µ 1 µ 2 c a 2 b 2 = ( a 2 µ 1 c 2 + a 2 µ 2 c 2 µ 2 b 2 c 2 + a 2 b 2 µ 1 ) r 2 (A.8) (A.9) (A.10) The displacement field components are: u 1 θ (r) = c b 2 µ 2 (r 2 a 2 ) ( a 2 µ 1 c 2 + a 2 µ 2 c 2 µ 2 b 2 c 2 + a 2 b 2 µ 1 )r u 1 r (r) = 0 u 2 θ(r) = c r2 a 2 µ 2 + r 2 a 2 µ 1 + r 2 µ 2 b 2 a 2 b 2 µ 1 ( a 2 µ 1 c 2 + a 2 µ 2 c 2 µ 2 b 2 c 2 + a 2 b 2 µ 1 )r u 2 r (r) = 0 (A.11) (A.12) (A.13) (A.14) And the pressure evolution is: p 1 (r) = W 3 r 3 + W 2 r 2 + W 1 r + W 0 (A.15) p 2 (r) = V 3 r 3 + V 2 r 2 + V 1 r + V 0 (A.16) 28
29 Moreover, in the example section 4.2.1, we have considered: µ 1 = 1/3 µ 2 = 10/3 a = 0.4 b = 1.0 c = 2.0 W 3 = 10.0 W 2 = 20.0 W 1 = 5.0 W 0 = 10.0 V 3 = 0.0 V 2 = 10.0 V 1 = 5.0 V 0 = 0.0 B Analytical solution for a straight interface The problem (see figure 15) is built so that the interface between the two materials is straight. The loading is chosen such that: b 1x = 6µ 1 x b 1y = 6 µ 1 y 2 µ y 2 b 2x = 24 µ 2 x b 2y = 24 µ 2 y 2 µ y 2 (B.1) (B.2) (B.3) (B.4) T 1x = 6 µ 1 x T 1y = 6 µ 1 y 2 µ y 2 T 2x = 24 µ 2 x T 2y = 24 µ 2 y 2 µ y 2 (B.5) (B.6) (B.7) (B.8) Under these boundary conditions, the stress field components are: σ 1 = 2 µ 1 ( 3 y 2 2 µ 2y µ 1 ) y 3 µ 1 ( ( 6 y + 2 µ 2 ) ) µ 1 x + 2 µ 1 ( (6 y + 2 µ 2 /µ 1 ) x + 2) 2 µ 1 ( 3 y µ 2y µ 1 ) y 3 (B.9) σ 2 = 2 µ ( 2 (12 y 2 2 y) 2 y 3 µ 2 ( 24 y + 2) x + 2 µ 1 ( ) µ 2 ( 24 y + 2) x + 2 µ 1 µ 2 2 µ2 ( 12 y y) 2 y 3 ) µ 2 (B.10) The displacement field is: 29
30
31 x y Fig. B.3. Bimaterial problem with straight interface: pressure field. 31
Some improvements of Xfem for cracked domains
Some improvements of Xfem for cracked domains E. Chahine 1, P. Laborde 2, J. Pommier 1, Y. Renard 3 and M. Salaün 4 (1) INSA Toulouse, laboratoire MIP, CNRS UMR 5640, Complexe scientifique de Rangueil,
More informationPseudo-divergence-free element free Galerkin method for incompressible fluid flow
International Workshop on MeshFree Methods 23 Pseudo-divergence-free element free Galerkin method for incompressible fluid flow Y. Vidal () and A. Huerta (2) Abstract: Incompressible modelling in finite
More informationModélisation d'interfaces linéaires et non linéaires dans le cadre X-FEM
Modélisation d'interfaces linéaires et non linéaires dans le cadre X-FEM Nicolas MOËS Ecole Centrale de Nantes, FRANCE GeM Institute - CNRS Journées Interface Numérique, GDR-IFS, CNAM, PARIS, 14-15 Mai
More informationDevelopment of X-FEM methodology and study on mixed-mode crack propagation
Acta Mech. Sin. (2011) 27(3):406 415 DOI 10.1007/s10409-011-0436-x RESEARCH PAPER Development of X-FEM methodology and study on mixed-mode crack propagation Zhuo Zhuang Bin-Bin Cheng Received: 2 February
More informationLocking phenomena in Computational Mechanics: nearly incompressible materials and plate problems
Locking phenomena in Computational Mechanics: nearly incompressible materials and plate problems C. Lovadina Dipartimento di Matematica Univ. di Pavia IMATI-CNR, Pavia Bologna, September, the 18th 2006
More informationStress analysis around crack tips in finite strain problems using the extended finite element method
Stress analysis around crack tips in finite strain problems using the etended finite element method Grégory Legrain, Nicolas Moes, Erwan Verron To cite this version: Grégory Legrain, Nicolas Moes, Erwan
More informationMulti-scale digital image correlation of strain localization
Multi-scale digital image correlation of strain localization J. Marty a, J. Réthoré a, A. Combescure a a. Laboratoire de Mécanique des Contacts et des Strcutures, INSA Lyon / UMR CNRS 5259 2 Avenue des
More informationENHANCED STRAIN METHODS FOR ELASTICITY PROBLEMS
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, S. Korotov, E. Oñate, J. Périaux, and D. Knörzer (eds.) Jyväskylä, 24 28 July 2004
More informationLocking in the incompressible limit: pseudo-divergence-free element free Galerkin
Locking in the incompressible limit: pseudodivergencefree element free Galerkin Yolanda Vidal Pierre Villon and Antonio Huerta Departament de Matemàtica Aplicada III, ES de Ingenieros de Caminos, Canales
More informationDevelopment of discontinuous Galerkin method for linear strain gradient elasticity
Development of discontinuous Galerkin method for linear strain gradient elasticity R Bala Chandran Computation for Design and Optimizaton Massachusetts Institute of Technology Cambridge, MA L. Noels* Aerospace
More informationINNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD
INNOVATIVE FINITE ELEMENT METHODS FOR PLATES* DOUGLAS N. ARNOLD Abstract. Finite element methods for the Reissner Mindlin plate theory are discussed. Methods in which both the tranverse displacement and
More informationAn exact reanalysis algorithm using incremental Cholesky factorization and its application to crack growth modeling
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 01; 91:158 14 Published online 5 June 01 in Wiley Online Library (wileyonlinelibrary.com). DOI: 10.100/nme.4 SHORT
More information3D XFEM Modelling of Imperfect Interfaces
3D XFEM Modelling of Imperfect Interfaces Elena Benvenuti 1, Giulio Ventura 2, Antonio Tralli 1, Nicola Ponara 1 1 Dipartimento di Ingegneria, Universit di Ferrara, Italy E-mail: elena.benvenuti@unife.it,
More informationSome New Elements for the Reissner Mindlin Plate Model
Boundary Value Problems for Partial Differential Equations and Applications, J.-L. Lions and C. Baiocchi, eds., Masson, 1993, pp. 287 292. Some New Elements for the Reissner Mindlin Plate Model Douglas
More informationStudies of Bimaterial Interface Fracture with Peridynamics Fang Wang 1, Lisheng Liu 2, *, Qiwen Liu 1, Zhenyu Zhang 1, Lin Su 1 & Dan Xue 1
International Power, Electronics and Materials Engineering Conference (IPEMEC 2015) Studies of Bimaterial Interface Fracture with Peridynamics Fang Wang 1, Lisheng Liu 2, *, Qiwen Liu 1, Zhenyu Zhang 1,
More informationA truly meshless Galerkin method based on a moving least squares quadrature
A truly meshless Galerkin method based on a moving least squares quadrature Marc Duflot, Hung Nguyen-Dang Abstract A new body integration technique is presented and applied to the evaluation of the stiffness
More informationMIXED FINITE ELEMENTS FOR PLATES. Ricardo G. Durán Universidad de Buenos Aires
MIXED FINITE ELEMENTS FOR PLATES Ricardo G. Durán Universidad de Buenos Aires - Necessity of 2D models. - Reissner-Mindlin Equations. - Finite Element Approximations. - Locking. - Mixed interpolation or
More informationCohesive band model: a triaxiality-dependent cohesive model for damage to crack transition in a non-local implicit discontinuous Galerkin framework
University of Liège Aerospace & Mechanical Engineering Cohesive band model: a triaxiality-dependent cohesive model for damage to crack transition in a non-local implicit discontinuous Galerkin framework
More informationA Stabilized Formulation for Incompressible Elasticity Using Linear Displacement and Pressure Interpolations
A Stabilized Formulation for Incompressible Elasticity Using Linear Displacement and Pressure Interpolations M. Chiumenti y, Q. Valverde Λ, C. Agelet de Saracibar y and M. Cervera y ydept. de Resistència
More informationCheckerboard instabilities in topological shape optimization algorithms
Checkerboard instabilities in topological shape optimization algorithms Eric Bonnetier, François Jouve Abstract Checkerboards instabilities for an algorithm of topological design are studied on a simple
More informationCRACK ANALYSIS IN MAGNETOELECTROELASTIC MEDIA USING THE EXTENDED FINITE ELEMENT METHOD
International Conference on Extended Finite Element Methods Recent Developments and Applications XFEM 2009 T.P. Fries and A. Zilian (Eds) c RWTH Aachen, Germany, 2009 CRACK ANALYSIS IN MAGNETOELECTROELASTIC
More informationComparative study of FEM and XFEM
et International Journal on Emerging Technologies 4(2): 47-55(2013) ISSN No. (Print) : 0975-8364 ISSN No. (Online) : 2249-3255 Comparative study of FEM and XFEM Ajay Singh*, Jagdeep Kumar**, Vinod Dhull***
More informationModel-independent approaches for the XFEM in fracture mechanics
Model-independent approaches for the XFEM in fracture mechanics Safdar Abbas 1 Alaskar Alizada 2 and Thomas-Peter Fries 2 1 Aachen Institute for Computational Engineering Science (AICES), RWTH Aachen University,
More informationChapter 2. General concepts. 2.1 The Navier-Stokes equations
Chapter 2 General concepts 2.1 The Navier-Stokes equations The Navier-Stokes equations model the fluid mechanics. This set of differential equations describes the motion of a fluid. In the present work
More informationA super convergent hybridisable discontinuous Galerkin method for linear elasticity
A super convergent hybridisable discontinuous Galerkin method for linear elasticity arxiv:82.223v [math.na] 6 Feb 28 Ruben Sevilla Zienkiewicz Centre for Computational Engineering, College of Engineering,
More informationA weak penalty formulation remedying traction oscillations in interface elements
c 2016. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/. A weak penalty formulation remedying traction oscillations in interface
More informationGeneralized Finite Element Methods for Three Dimensional Structural Mechanics Problems. C. A. Duarte. I. Babuška and J. T. Oden
Generalized Finite Element Methods for Three Dimensional Structural Mechanics Problems C. A. Duarte COMCO, Inc., 7800 Shoal Creek Blvd. Suite 290E Austin, Texas, 78757, USA I. Babuška and J. T. Oden TICAM,
More informationAn Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation
An Energy Dissipative Constitutive Model for Multi-Surface Interfaces at Weld Defect Sites in Ultrasonic Consolidation Nachiket Patil, Deepankar Pal and Brent E. Stucker Industrial Engineering, University
More informationPowerful Modelling Techniques in Abaqus to Simulate
Powerful Modelling Techniques in Abaqus to Simulate Necking and Delamination of Laminated Composites D. F. Zhang, K.M. Mao, Md. S. Islam, E. Andreasson, Nasir Mehmood, S. Kao-Walter Email: sharon.kao-walter@bth.se
More informationCAST3M IMPLEMENTATION OF THE EXTENDED FINITE ELEMENT METHOD FOR COHESIVE CRACK
Vietnam Journal of Mechanics, VAST, Vol. 33, No. 1 (2011), pp. 55 64 CAST3M IMPLEMENTATION OF THE EXTENDED FINITE ELEMENT METHOD FOR COHESIVE CRACK Nguyen Truong Giang, Ngo Huong Nhu Institute of Mechanics,
More informationJ. Liou Tulsa Research Center Amoco Production Company Tulsa, OK 74102, USA. Received 23 August 1990 Revised manuscript received 24 October 1990
Computer Methods in Applied Mechanics and Engineering, 94 (1992) 339 351 1 A NEW STRATEGY FOR FINITE ELEMENT COMPUTATIONS INVOLVING MOVING BOUNDARIES AND INTERFACES THE DEFORMING-SPATIAL-DOMAIN/SPACE-TIME
More information3D crack propagation with XFEM cohesive elements
Transactions, SMiRT-23, Paper ID 472 3D crack propagation with XFEM cohesive elements Patrick Massin 1, Guilhem Ferté 1,2, and Nicolas Moës 2 1 Institute of Mechanical Sciences and Industrial Applications,
More informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationANALYSIS OF MIXED FINITE ELEMENTS FOR LAMINATED COMPOSITE PLATES
ANALYSIS OF MIXED FINITE ELEMENTS FOR LAMINATED COMPOSITE PLATES F. Auricchio Dipartimento di Meccanica Strutturale Università di Pavia, Italy C. Lovadina Dipartimento di Ingegneria Meccanica e Strutturale
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More informationCohesive Band Model: a triaxiality-dependent cohesive model inside an implicit non-local damage to crack transition framework
University of Liège Aerospace & Mechanical Engineering MS3: Abstract 131573 - CFRAC2017 Cohesive Band Model: a triaxiality-dependent cohesive model inside an implicit non-local damage to crack transition
More informationA NEW CLASS OF MIXED FINITE ELEMENT METHODS FOR REISSNER MINDLIN PLATES
A NEW CLASS OF IXED FINITE ELEENT ETHODS FOR REISSNER INDLIN PLATES C. LOVADINA Abstract. A new class of finite elements for Reissner indlin plate problem is presented. The family is based on a modified
More informationFinite Element Computations of Complex Stress Intensity Factor Magnitude of Interfacial Crack in Bi-Material Composites
16 th International Conference on AEROSPACE SCIENCES & AVIATION TECHNOLOGY, ASAT - 16 May 26-28, 215, E-Mail: asat@mtc.edu.eg Military Technical College, Kobry Elkobbah, Cairo, Egypt Tel : +(22) 2425292
More informationGeneralized Finite Element Method with global-local enrichments for nonlinear fracture analysis
Generalized Finite Element Method with global-local enrichments for nonlinear fracture analysis D.-J. Kim, C.A. Duarte Department of Civil and Environmental Engr., University of Illinois at Urbana-Champaign,
More information1 Exercise: Linear, incompressible Stokes flow with FE
Figure 1: Pressure and velocity solution for a sinking, fluid slab impinging on viscosity contrast problem. 1 Exercise: Linear, incompressible Stokes flow with FE Reading Hughes (2000), sec. 4.2-4.4 Dabrowski
More informationThis document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore.
This document is downloaded from DR-NTU, Nanyang Technological University Library, Singapore. Title An enriched 6-node MITC plate element for yield line analysis Author(s) Xu, Jin; Lee, Chi King; Tan,
More informationA new formulation for imposing Dirichlet boundary conditions on non-matching meshes
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Meth. Engng 2015; 103:430 444 Published online 19 March 2015 in Wiley Online Library (wileyonlinelibrary.com)..4898 A new formulation
More informationSome applications of discontinuous Galerkin methods in solid mechanics
Some applications of discontinuous Galerkin methods in solid mechanics Adrian Lew, Alex Ten Eyck, and Ramsharan Rangarajan Mechanical Engineering, Stanford University, USA. Email: {lewa,alabute,rram}@stanford.edu
More informationComputers and Structures
Computers and Structures 87 (2009) 1451 1460 Contents lists available at ScienceDirect Computers and Structures journal homepage: www.elsevier.com/locate/compstruc A triangular six-node shell element Do-Nyun
More informationDynamic crack propagation analysis of orthotropic media by the extended finite element method
Int J Fract (2010 161:21 39 DOI 10.1007/s10704-009-9423-7 ORIGINAL PAPER Dynamic crack propagation analysis of orthotropic media by the extended finite element method D. Motamedi S. Mohammadi Received:
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationFLOATING NODE METHOD AND VIRTUAL CRACK CLOSURE TECHNIQUE FOR MODELING MATRIX CRACKING- DELAMINATION MIGRATION
THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS FLOATING NODE METHOD AND VIRTUAL CRACK CLOSURE TECHNIQUE FOR MODELING MATRIX CRACKING- DELAMINATION MIGRATION N. V. De Carvalho 1*, B. Y. Chen
More informationA Locking-Free MHM Method for Elasticity
Trabalho apresentado no CNMAC, Gramado - RS, 2016. Proceeding Series of the Brazilian Society of Computational and Applied Mathematics A Locking-Free MHM Method for Elasticity Weslley S. Pereira 1 Frédéric
More informationThe Finite Element Method for Solid and Structural Mechanics
The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in
More informationOn Mixed Methods for Signorini Problems
Annals of University of Craiova, Math. Comp. Sci. Ser. Volume 30, 2003, Pages 45 52 ISS: 1223-6934 On Mixed Methods for Signorini Problems Faker Ben Belgacem, Yves Renard, and Leila Slimane Abstract. We
More informationMixed-mode crack propagation using a Hybrid Analytical and extended finite element method
Mixed-mode crack propagation using a Hybrid Analytical and extended finite element method Julien Réthoré a, Stéphane Roux b, François Hild b a LaMCoS, Université de Lyon/INSA Lyon/CNRS, 20 Avenue des Sciences,
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationIntroduction to fracture mechanics
Introduction to fracture mechanics Prof. Dr. Eleni Chatzi Dr. Giuseppe Abbiati, Dr. Konstantinos Agathos Lecture 6-9 November, 2017 Institute of Structural Engineering, ETH Zu rich November 9, 2017 Institute
More informationIMPLEMENTATION OF FUNCTIONAL-TYPE A POSTERIORI ERROR ESTIMATES FOR LINEAR PROBLEMS IN SOLID MECHANICS. Maksim Frolov
IMPLEMENTATION OF FUNCTIONAL-TYPE A POSTERIORI ERROR ESTIMATES FOR LINEAR PROBLEMS IN SOLID MECHANICS Maksim Frolov Peter the Great St.Petersburg Polytechnic University AANMPDE11, FINLAND 06 10 August
More informationComputational homogenization of microfractured continua using weakly periodic boundary conditions
c 2015. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/. Computational homogenization of microfractured continua using weakly
More informationError Analyses of Stabilized Pseudo-Derivative MLS Methods
Department of Mathematical Sciences Revision: April 11, 2015 Error Analyses of Stabilized Pseudo-Derivative MLS Methods Don French, M. Osorio and J. Clack Moving least square (MLS) methods are a popular
More informationJournal of Theoretical and Applied Mechanics, Sofia, 2012, vol. 42, No. 1, pp
Journal of Theoretical and Applied Mechanics, Sofia, 2012, vol. 42, No. 1, pp. 23 40 SOLID MECHANICS CALCULATION OF STRESS INTENSITY FACTORS BASED ON FORCE-DISPLACEMENT CURVE USING ELEMENT FREE GALERKIN
More informationAdaptive C1 Macroelements for Fourth Order and Divergence-Free Problems
Adaptive C1 Macroelements for Fourth Order and Divergence-Free Problems Roy Stogner Computational Fluid Dynamics Lab Institute for Computational Engineering and Sciences University of Texas at Austin March
More informationDISPLACEMENT/PRESSURE BASED MIXED FINITE ELEMENT FORMULATIONS FOR ACOUSTIC FLUID STRUCTURE INTERACTION PROBLEMS
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, VOL. 40, 2001 2017 (1997) DISPLACEMENT/PRESSURE BASED MIXED FINITE ELEMENT FORMULATIONS FOR ACOUSTIC FLUID STRUCTURE INTERACTION PROBLEMS XIAODONG
More informationOn the implementation of the extended Finite Element Method (XFEM) for interface problems
On the implementation of the extended Finite Element Method (XFEM) for interface problems Thomas Carraro and Sven Wetterauer Institute for Applied Mathematics, Heidelberg University July 16, 2015 arxiv:1507.04238v1
More informationDiscontinuous Galerkin methods for the Stokes equations using divergence-free approximations
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 27; : 9 [Version: 22/9/8 v.] Discontinuous Galerkin methods for the Stokes equations using divergence-free approximations
More informationThe Virtual Element Method: an introduction with focus on fluid flows
The Virtual Element Method: an introduction with focus on fluid flows L. Beirão da Veiga Dipartimento di Matematica e Applicazioni Università di Milano-Bicocca VEM people (or group representatives): P.
More informationFig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double cantilever beam (DCB) specimens.
a). Cohesive Failure b). Interfacial Failure c). Oscillatory Failure d). Alternating Failure Fig. 1. Different locus of failure and crack trajectories observed in mode I testing of adhesively bonded double
More informationAdvances in Engineering Software
Advances in Engineering Software 41 (21) 712 728 Contents lists available at ScienceDirect Advances in Engineering Software journal homepage: www.elsevier.com/locate/advengsoft The quadratic MITC plate
More informationA posteriori error estimates applied to flow in a channel with corners
Mathematics and Computers in Simulation 61 (2003) 375 383 A posteriori error estimates applied to flow in a channel with corners Pavel Burda a,, Jaroslav Novotný b, Bedřich Sousedík a a Department of Mathematics,
More informationLEAST-SQUARES FINITE ELEMENT MODELS
LEAST-SQUARES FINITE ELEMENT MODELS General idea of the least-squares formulation applied to an abstract boundary-value problem Works of our group Application to Poisson s equation Application to flows
More informationPSEUDO ELASTIC ANALYSIS OF MATERIAL NON-LINEAR PROBLEMS USING ELEMENT FREE GALERKIN METHOD
Journal of the Chinese Institute of Engineers, Vol. 27, No. 4, pp. 505-516 (2004) 505 PSEUDO ELASTIC ANALYSIS OF MATERIAL NON-LINEAR PROBLEMS USING ELEMENT FREE GALERKIN METHOD Raju Sethuraman* and Cherku
More informationA NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION
A NOTE ON THE LADYŽENSKAJA-BABUŠKA-BREZZI CONDITION JOHNNY GUZMÁN, ABNER J. SALGADO, AND FRANCISCO-JAVIER SAYAS Abstract. The analysis of finite-element-like Galerkin discretization techniques for the
More informationFast multipole boundary element method for the analysis of plates with many holes
Arch. Mech., 59, 4 5, pp. 385 401, Warszawa 2007 Fast multipole boundary element method for the analysis of plates with many holes J. PTASZNY, P. FEDELIŃSKI Department of Strength of Materials and Computational
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationQUASI-STATIC AND DYNAMIC SIMULATION OF SHEET METAL FORMING PROCESSES USING LINEAR AND QUADRATIC SOLID- SHELL ELEMENTS
Proceedings of the 6th International Conference on Mechanics and Materials in Design, Editors: J.F. Silva Gomes & S.A. Meguid, P.Delgada/Azores, 26-30 July 2015 PAPER REF: 5456 QUASI-STATIC AND DYNAMIC
More informationarxiv: v1 [math.na] 27 Jan 2016
Virtual Element Method for fourth order problems: L 2 estimates Claudia Chinosi a, L. Donatella Marini b arxiv:1601.07484v1 [math.na] 27 Jan 2016 a Dipartimento di Scienze e Innovazione Tecnologica, Università
More informationINCOMPATIBLE BUBBLES: A GENERALIZED VARIATIONAL FRAMEWORK INCLUDING INCOMPATIBLE MODES AND ENHANCED ASSUMED STRAINS
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2004 P. Neittaanmäki, T. Rossi, S. Korotov, E. Oñate, J. Périaux, and D. Knörzer (eds.) Jyväskylä, 24 28 July 2004
More informationOverview of Extended Finite Element Method
Chapter 1 Overview of Extended Finite Element Method Chapter Outline 1.1 Significance of Studying Computational Fracture Mechanics 1 1.2 Introduction to X-FEM 2 1.3 Research Status and Development of X-FEM
More informationA RATE-DEPENDENT MULTI-SCALE CRACK MODEL FOR CONCRETE
VIII International Conference on Fracture echanics of Concrete and Concrete Structures FraCoS-8 J.G.. Van ier, G. Ruiz, C. Andrade, R.C. Yu and X.X. Zhang (Eds) A RATE-DEPENDENT ULTI-SCALE CRACK ODEL FOR
More informationAn Extended Finite Element Method for a Two-Phase Stokes problem
XFEM project An Extended Finite Element Method for a Two-Phase Stokes problem P. Lederer, C. Pfeiler, C. Wintersteiger Advisor: Dr. C. Lehrenfeld August 5, 2015 Contents 1 Problem description 2 1.1 Physics.........................................
More informationTHE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER QUADRILATERAL MESHES. Zhong-Ci Shi and Xuejun Xu.
DISCRETE AND CONTINUOUS Website: http://aimsciences.org DYNAMICAL SYSTEMS SERIES B Volume 9, Number, January 2008 pp. 63 82 THE PATCH RECOVERY FOR FINITE ELEMENT APPROXIMATION OF ELASTICITY PROBLEMS UNDER
More informationAn Eulerian-Lagrangian approach for fluid-structure coupled systems using X-FEM Antoine Legay / CNAM-Paris
An Eulerian-Lagrangian approach for fluid-structure coupled systems using X-FEM Antoine Legay / CNAM-Paris Collaborations: A. Tralli, P. Gaudenzi, Università di Roma La Sapienza, Italia 8th USNCCM, 25-27th
More informationStatic & Dynamic. Analysis of Structures. Edward L.Wilson. University of California, Berkeley. Fourth Edition. Professor Emeritus of Civil Engineering
Static & Dynamic Analysis of Structures A Physical Approach With Emphasis on Earthquake Engineering Edward LWilson Professor Emeritus of Civil Engineering University of California, Berkeley Fourth Edition
More informationGlasgow, Rankine building, G12 8LT, Glasgow, UK; Tel: +44(0) ;
Three-dimensional crack initiation, propagation, branching and junction in non-linear materials by an extended meshfree method without asymptotic enrichment Stéphane Bordas a Timon Rabczuk b Goangseup
More informationA Moving Least Squares weighting function for the Element-free Galerkin Method which almost fulfills essential boundary conditions
Structural Engineering and Mechanics, Vol. 21, No. 3 (2005) 315-332 315 A Moving Least Squares weighting function for the Element-free Galerkin Method which almost fulfills essential boundary conditions
More informationCommon pitfalls while using FEM
Common pitfalls while using FEM J. Pamin Instytut Technologii Informatycznych w Inżynierii Lądowej Wydział Inżynierii Lądowej, Politechnika Krakowska e-mail: JPamin@L5.pk.edu.pl With thanks to: R. de Borst
More informationTREBALL FINAL DE MÀSTER
A continuous-discontinuous model to simulate crack branching in quasi-brittle failure Treball realitzat per: Jordi Feliu Fabà Dirigit per: Antonio Rodríguez Ferran Màster en: Enginyeria de Camins, Canals
More informationSMOOTHED PARTICLE HYDRODYNAMICS METHOD IN MODELING OF STRUCTURAL ELEMENTS UNDER HIGH DYNAMIC LOADS
The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China SMOOTHE PARTICLE HYROYAMICS METHO I MOELIG OF STRUCTURAL ELEMETS UER HIGH YAMIC LOAS. Asprone *, F. Auricchio, A. Reali,
More informationHomework Problems. ( σ 11 + σ 22 ) 2. cos (θ /2), ( σ θθ σ rr ) 2. ( σ 22 σ 11 ) 2
Engineering Sciences 47: Fracture Mechanics J. R. Rice, 1991 Homework Problems 1) Assuming that the stress field near a crack tip in a linear elastic solid is singular in the form σ ij = rλ Σ ij (θ), it
More informationA TIME-DEPENDENT DAMAGE LAW IN DEFORMABLE SOLID: A HOMOGENIZATION APPROACH
9th HSTAM International Congress on Mechanics Limassol, Cyprus, - July, A TIME-DEPENDENT DAMAGE LAW IN DEFORMABLE SOLID: A HOMOGENIZATION APPROACH Cristian Dascalu, Bertrand François, Laboratoire Sols
More informationLarge sliding contact along branched discontinuities with X-FEM
DOI 10.1007/s00466-012-0807-6 ORIGINAL PAPER Large sliding contact along branched discontinuities with X-FEM Maximilien Siavelis Martin L. E. Guiton Patrick Massin Nicolas Moës Received: 19 April 2012
More informationAlternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering
University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering
More informationArchetype-Blending Multiscale Continuum Method
Archetype-Blending Multiscale Continuum Method John A. Moore Professor Wing Kam Liu Northwestern University Mechanical Engineering 3/27/2014 1 1 Outline Background and Motivation Archetype-Blending Continuum
More informationPROPAGATION OF CURVED CRACKS IN HOMOGENEOUS AND GRADED MATERIALS
PROPAGATION OF CURVED CRACKS IN HOMOGENEOUS AND GRADED MATERIALS Abstract Matthew T. Tilbrook, Robert J. Moon and Mark Hoffman School of Materials Science and Engineering University of New South Wales,
More informationAbstract. 1 Introduction
Contact analysis for the modelling of anchors in concrete structures H. Walter*, L. Baillet** & M. Brunet* *Laboratoire de Mecanique des Solides **Laboratoire de Mecanique des Contacts-CNRS UMR 5514 Institut
More informationStabilized and Coupled FEM/EFG Approximations for Fluid Problems
COMPUTATIONAL MECHANICS WCCM VI in conjunction with APCOM4, Sept. 5-, 24, Beijing, China c 24 Tsinghua University Press & Springer-Verlag Stabilized and Coupled FEM/EFG Approximations for Fluid Problems
More informationComputational non-linear structural dynamics and energy-momentum integration schemes
icccbe 2010 Nottingham University Press Proceedings of the International Conference on Computing in Civil and Building Engineering W Tizani (Editor) Computational non-linear structural dynamics and energy-momentum
More information2008 by authors and 2008 Springer Science+Business Media
Antti H. Niemi, Harri Hakula, and Juhani Pitkäranta. 28. Point load on a shell. In: Karl Kunisch, Günther Of, and Olaf Steinbach (editors). Numerical Mathematics and Advanced Applications. Proceedings
More informationSTABILIZED FINITE ELEMENT METHODS TO DEAL WITH INCOMPRESSIBILITY IN SOLID MECHANICS IN FINITE STRAINS
European Congress on Computational Methods in Applied Sciences and Engineering ECCOMAS 2012) J. Eberhardsteiner et.al. eds.) Vienna, Austria, September 10-14, 2012 STABILIZED FINITE ELEMENT METHODS TO
More informationPSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS
PSEUDO-COMPRESSIBILITY METHODS FOR THE UNSTEADY INCOMPRESSIBLE NAVIER-STOKES EQUATIONS Jie Shen Department of Mathematics, Penn State University University Par, PA 1680, USA Abstract. We present in this
More informationA local projection stabilization of fictitious domain method for elliptic boundary value problems
A local projection stabilization of fictitious domain method for elliptic boundary value problems Saber Amdouni, Maher Moakher, Yves Renard. Oct. 25, 2013 Abstract In this paper, a new consistent method
More information