首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
    
A new technique for treating the mechanical interactions of overlapping finite element meshes is presented. Such methods can be useful for numerous applications, for example, fluid–solid interaction with a superposed meshed solid on an Eulerian background fluid grid. In this work, we consider the interaction of two elastic domains: one mesh is the foreground and defines the surface of interaction, the other is a background mesh and is often a structured grid. Many of the previously proposed methods employ surface defined Lagrange multipliers or penalties to enforce the boundary constraints. It has become apparent that these methods will cause mesh locking under certain conditions. Appropriately applied, the Nitsche method can overcome this locking, but, in its canonical form, is generally not applicable to non‐linear materials such as hyperelastics. The relationship between interior point penalty, discontinuous Galerkin and Nitsche's method is well known. Based on this relationship, a nonlinear theory analogous to the Nitsche method is proposed to treat nonlinear materials in an embedded mesh. Here, a discontinuous Galerkin derivative based on a lifting of the interface surface integrals provides a consistent treatment for non‐linear materials and demonstrates good behavior in example problems. Published 2012. This article is a US Government work and is in the public domain in the USA.  相似文献   

2.
    
Enriched finite element approaches such as the extended finite element method provide a framework for constructing approximations to solutions of non‐smooth problems. Internal features, such as boundaries, are represented in such methods by using discontinuous enrichment of the standard finite element basis. Within such frameworks, however, imposition of interface constraints and/or constitutive relations can cause unexpected difficulties, depending upon how relevant fields are interpolated on un‐gridded interfaces. This work address the stabilized treatment of constraints in an enriched finite element context. Both the Lagrange multiplier and penalty enforcement of tied constraints for an arbitrary boundary represented in an enriched finite element context can lead to instabilities and artificial oscillations in the traction fields. We demonstrate two alternative variational methods that can be used to enforce the constraints in a stable manner. In a ‘bubble‐stabilized approach,’ fine‐scale degrees of freedom are added over elements supporting the interface. The variational form can be shown to have a similar form to a second approach we consider, Nitsche's method, with the exception that the stabilization terms follow directly from the bubble functions. In this work, we examine alternative variational methods for enforcing a tied constraint on an enriched interface in the context of two‐dimensional elasticity. We examine several benchmark problems in elasticity, and show that only Nitsche's method and the bubble‐stabilization approach produce stable traction fields over internal boundaries. We also demonstrate a novel difference between the penalty method and Nitsche's method in that the latter passes the patch test exactly, regardless of the stabilization parameter's magnitude. Results for more complicated geometries and triple interface junctions are also presented. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

3.
    
In this paper we present a strategy for the simulation of a propagating crack under mixed mode linear elastic conditions using a discontinuous finite element method. A key issue is to accurately compute the incremental change in the kink angle of the propagating crack during subsequent steps. We have chosen to work with a domain formulation of the material force vector as a criteria for the propagation direction. We describe the theoretical background together with the numerical implementation in detail and show some results for different loading conditions. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

4.
    
In this paper, we propose a way to weakly prescribe Dirichlet boundary conditions in embedded finite element meshes. The key feature of the method is that the algorithmic parameter of the formulation which allows to ensure stability is independent of the numerical approximation, relatively small, and can be fixed a priori. Moreover, the formulation is symmetric for symmetric problems. An additional element-discontinuous stress field is used to enforce the boundary conditions in the Poisson problem. Additional terms are required in order to guarantee stability in the convection–diffusion equation and the Stokes problem. The proposed method is then easily extended to the transient Navier–Stokes equations. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

5.
    
We analyze several possibilities to prescribe boundary conditions in the context of immersed boundary methods. As basic approximation technique we consider the finite element method with a mesh that does not match the boundary of the computational domain, and therefore Dirichlet boundary conditions need to be prescribed in an approximate way. As starting variational approach we consider Nitsche's methods, and we then move to two options that yield non‐symmetric problems but that turned out to be robust and efficient. The essential idea is to use the degrees of freedom of certain nodes of the finite element mesh to minimize the difference between the exact and the approximated boundary condition. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

6.
    
A key challenge while employing non‐interpolatory basis functions in finite‐element methods is the robust imposition of Dirichlet boundary conditions. The current work studies the weak enforcement of such conditions for B‐spline basis functions, with application to both second‐ and fourth‐order problems. This is achieved using concepts borrowed from Nitsche's method, which is a stabilized method for imposing constraints on surfaces. Conditions for the stability of the system of equations are derived for each class of problem. Stability parameters in the Nitsche weak form are then evaluated by solving a local generalized eigenvalue problem at the Dirichlet boundary. The approach is designed to work equally well when the grid used to build the splines conforms to the physical boundary of interest as well as to the more general case when it does not. Through several numerical examples, the approach is shown to yield optimal rates of convergence. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

7.
    
A stabilized variational formulation, based on Nitsche's method for enforcing boundary constraints, leads to an efficient procedure for embedding kinematic boundary conditions in thin plate bending. The absence of kinematic admissibility constraints allows the use of non‐conforming meshes with non‐interpolatory approximations, thereby providing added flexibility in addressing the C1‐continuity requirements typical of these problems. Work‐conjugate pairs weakly enforce kinematic boundary conditions. The pointwise enforcement of corner deflections is key to good performance in the presence of corners. Stabilization parameters are determined from local generalized eigenvalue problems, guaranteeing coercivity of the discrete bilinear form. The accuracy of the approach is verified by representative computations with bicubic C2 B‐splines, exhibiting optimal rates of convergence and robust performance with respect to values of the stabilization parameters. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

8.
A method for removing the numerical instability of the conventional Green's formula in the vicinity of the boundary surface is proposed. The approach uses a pair of Green's formulae written both inside and outside the volume of interest, even when used for solving a single-phase problem. It is shown that cancellation between the pair yields an alternative, singularity-free integral representation that is amenable to numerical calculation. The derivation of the formula is given explicitly and is accompanied by several numerical test results for validation.  相似文献   

9.
    
Simulation approaches for fluid-structure-contact interaction, especially if requested to be consistent even down to the real contact scenarios, belong to the most challenging and still unsolved problems in computational mechanics. The main challenges are 2-fold—one is to have a correct physical model for this scenario, and the other is to have a numerical method that is capable of working and being consistent down to a zero gap. Moreover, when analyzing such challenging setups of fluid-structure interaction, which include contact of submersed solid components, it gets obvious that the influence of surface roughness effects is essential for a physical consistent modeling of such configurations. To capture this system behavior, we present a continuum mechanical model that is able to include the effects of the surface microstructure in a fluid-structure-contact interaction framework. An averaged representation for the mixture of fluid and solid on the rough surfaces, which is of major interest for the macroscopic response of such a system, is introduced therein. The inherent coupling of the macroscopic fluid flow and the flow inside the rough surfaces, the stress exchange of all contacting solid bodies involved, and the interaction between fluid and solid are included in the construction of the model. Although the physical model is not restricted to finite element–based methods, a numerical approach with its core based on the cut finite element method, enabling topological changes of the fluid domain to solve the presented model numerically, is introduced. Such a cut finite element method–based approach is able to deal with the numerical challenges mentioned above. Different test cases give a perspective toward the potential capabilities of the presented physical model and numerical approach.  相似文献   

10.
Cut finite element method–based approaches toward challenging fluid-structure interaction (FSI) are proposed. The different considered methods combine the advantages of competing novel Eulerian (fixed grid) and established arbitrary Lagrangian-Eulerian (moving mesh) finite element formulations for the fluid. The objective is to highlight the benefit of using cut finite element techniques for moving-domain problems and to demonstrate their high potential with regard to simplified mesh generation, treatment of large structural motions in surrounding flows, capturing boundary layers, their ability to deal with topological changes in the fluid phase, and their general straightforward extensibility to other coupled multiphysics problems. In addition to a pure fixed-grid FSI method, advanced fluid-domain decomposition techniques are also considered, leading to highly flexible discretization methods for the FSI problem. All stabilized formulations include Nitsche-based weak coupling of the phases supported by the ghost penalty technique for the flow field. For the resulting systems, monolithic solution strategies are presented. Various two- and three-dimensional FSI cases of different complexity levels validate the methods and demonstrate their capabilities and limitations in different situations.  相似文献   

11.
12.
    
We propose a new embedded finite element method to simulate partial differential equations over domains with internal interfaces. Our approach belongs to the family of surrogate/approximate interface methods and relies on the idea of shifting the location and value of jump interface conditions. This choice has the goal of preserving optimal convergence rates while avoiding small cut cells and related problematic issues, typical of traditional embedded methods. The proposed approach is accurate, robust, efficient, and simple to implement. We apply this concept to internal interface computations in the context of the mixed Poisson problem, also known as the Darcy flow problem. An extensive set of numerical tests is provided to demonstrate the performance of the proposed approach.  相似文献   

13.
The paper introduces a methodology for numerical simulation of landslides experiencing considerable deformations and topological changes. Within an interface capturing approach, all interfaces are implicitly described by specifically defined level‐set functions allowing arbitrarily evolving complex topologies. The transient interface evolution is obtained by solving the level‐set equation driven by the physical velocity field for all three level‐set functions in a block Jacobi approach. The three boundary‐coupled fluid‐like continua involved are modeled as incompressible, governed by a generalized non‐Newtonian material law taking into account capillary pressure at moving fluid–fluid interfaces. The weighted residual formulation of the level‐set equations and the fluid equations is discretized with finite elements in space and time using velocity and pressure as unknown variables. Non‐smooth solution characteristics are represented by enriched approximations to fluid velocity (weak discontinuity) and fluid pressure (strong discontinuity). Special attention is given to the construction of enriched approximations for elements containing evolving triple junctions. Numerical examples of three‐fluid tank sloshing and air–water‐liquefied soil landslides demonstrate the potential and applicability of the method in geotechnical investigations. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

14.
15.
    
In this paper, several boundary element regularization methods, such as iterative, conjugate gradient, Tikhonov regularization and singular value decomposition methods, for solving the Cauchy problem associated to the Helmholtz equation are developed and compared. Regularizing stopping criteria are developed and the convergence, as well as the stability, of the numerical methods proposed are analysed. The Cauchy problem for the Helmholtz equation can be regularized by various methods, such as the general regularization methods presented in this paper, but more accurate results are obtained by classical methods, such as the singular value decomposition and the Tikhonov regularization methods. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

16.
    
A new higher‐order accurate method is proposed that combines the advantages of the classical p‐version of the FEM on body‐fitted meshes with embedded domain methods. A background mesh composed by higher‐order Lagrange elements is used. Boundaries and interfaces are described implicitly by the level set method and are within elements. In the elements cut by the boundaries or interfaces, an automatic decomposition into higher‐order accurate sub‐elements is realized. Therefore, the zero level sets are detected and meshed in a first step, which is called reconstruction. Then, based on the topological situation in the cut element, higher‐order sub‐elements are mapped to the two sides of the boundary or interface. The quality of the reconstruction and the mapping largely determines the properties of the resulting, automatically generated conforming mesh. It is found that optimal convergence rates are possible although the resulting sub‐elements are not always well‐shaped. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

17.
    
We discuss explicit coupling schemes for fluid‐structure interaction problems where the added mass effect is important. In this paper, we show the close relation between coupling schemes by using Nitsche's method and a Robin–Robin type coupling. In the latter case, the method may be implemented either using boundary integrals of the stresses or the more conventional discrete lifting operators. Recalling the explicit method proposed in Comput. Methods Appl. Mech. Engrg. 198(5‐8):766–784, 2009, we make the observation that this scheme is stable under a hyperbolic type CFL condition, but that optimal accuracy imposes a parabolic type CFL conditions because of the splitting error. Two strategies to enhance the accuracy of the coupling scheme under the hyperbolic CFL‐condition are suggested, one using extrapolation and defect‐correction and one using a penalty‐free non‐symmetric Nitsche method. Finally, we illustrate the performance of the proposed schemes on some numerical examples in two and three space dimensions. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

18.
    
In this paper various formulations for the eddy current problem are presented. The formulations are based on solving directly for the magnetic field h, and they differ from each other mainly by how the field on the boundary is treated. The electromagnetic problem is studied in connection with the fivefold decomposition of the space of square integrable vector fields within a bounded region. This provides us with numerical approaches with clear signposts about how to solve the eddy current problem in multiply connected domains. Besides the fivefold decomposition, another essential tool in our approach is Whitney elements, as they provide the structure needed to retain consistency between the continuous and discrete problems. The paper demonstrates the usefulness of these mathematical tools in solving electromagnetic field problems. © 1998 John Wiley & Sons, Ltd.  相似文献   

19.
    
The boundary knot method is an inherently meshless, integration‐free, boundary‐type, radial basis function collocation technique for the solution of partial differential equations. In this paper, the method is applied to the solution of some inverse problems for the Helmholtz equation, including the highly ill‐posed Cauchy problem. Since the resulting matrix equation is badly ill‐conditioned, a regularized solution is obtained by employing truncated singular value decomposition, while the regularization parameter for the regularization method is provided by the L‐curve method. Numerical results are presented for both smooth and piecewise smooth geometry. The stability of the method with respect to the noise in the data is investigated by using simulated noisy data. The results show that the method is highly accurate, computationally efficient and stable, and can be a competitive alternative to existing methods for the numerical solution of the problems. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

20.
    
This paper presents a numerical method to solve the forward position problem in spatial mechanisms. The method may be incorporated in a software for the kinematic analysis of mechanisms, where the procedure is systematic and can be easily implemented, achieving a high degree of automation in simulation. The procedure presents high computational efficiency, enabling its incorporation in the control loop to solve the forward position problem in the case of a velocity control scheme. Also, in this paper preliminary results on the convergence of the proposed procedure are shown, and efficiency results of the method applied to representative spatial mechanisms are presented. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

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

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