共查询到20条相似文献,搜索用时 15 毫秒
1.
Andreas Apostolatos Robert Schmidt Roland Wüchner Kai‐Uwe Bletzinger 《International journal for numerical methods in engineering》2014,97(7):473-504
This paper provides a detailed elaboration and assessment of the most common domain decomposition methods for their application in isogeometric analysis. The methods comprise a penalty approach, Lagrange multiplier methods, and a Nitsche‐type method. For the Nitsche method, a new stabilized formulation is developed in the context of isogeometric analysis to guarantee coercivity. All these methods are investigated on problems of linear elasticity and eigenfrequency analysis in 2D. In particular, focus is put on non‐uniform rational B‐spline patches which join nonconformingly along their common interface. Thus, the application of isogeometric analysis is extended to multi‐patches, which can have an arbitrary parametrization on the adjacent edges. Moreover, it has been shown that the unique properties provided by isogeometric analysis, that is, high‐order functions and smoothness across the element boundaries, carry over for the analysis of multiple domains. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
2.
Ji‐Yeon Kim Sung‐Kie Youn 《International journal for numerical methods in engineering》2012,89(12):1559-1581
In the present work, an isogeometric contact analysis scheme using mortar method is proposed. Because the isogeometric analysis is employed for contact analysis, the geometric exactness of the contact region is maintained without any loss of geometric data because of geometry approximation. Thus, the proposed method can overcome underlying shortcomings that result from the geometric approximation of contact surfaces in the conventional finite element (FE)‐based contact analysis. For an isogeometric contact analysis, the schemes for treating the contact conditions and detecting the real contact surfaces are essentially required. In the proposed method, the mortar method is adopted as a nonconforming contact treatment scheme because it is expected to be in good harmony with the useful characteristics of nonuniform rational B‐spline A new matching algorithm is proposed to combine the mortar method with the isogeometric analysis to guarantee consistent contact surface information with the nonuniform rational B‐spline curve. The present scheme is verified by patch test and the well‐known problems which have theoretical solutions such as interference fit and the Hertzian contact problem. It is shown that the problems with curved contact surfaces which are difficult to treat by conventional approaches can be easily dealt with. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
3.
L. De Lorenzis İ. Temizer P. Wriggers G. Zavarise 《International journal for numerical methods in engineering》2011,87(13):1278-1300
This paper focuses on the application of NURBS‐based isogeometric analysis to Coulomb frictional contact problems between deformable bodies, in the context of large deformations. A mortar‐based approach is presented to treat the contact constraints, whereby the discretization of the continuum is performed with arbitrary order NURBS, as well as C0‐continuous Lagrange polynomial elements for comparison purposes. The numerical examples show that the proposed contact formulation in conjunction with the NURBS discretization delivers accurate and robust predictions. Results of lower quality are obtained from the Lagrange discretization, as well as from a different contact formulation based on the enforcement of the contact constraints at every integration point on the contact surface. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
4.
Margarita Chasapi Wolfgang Dornisch Sven Klinkel 《International journal for numerical methods in engineering》2020,121(14):3206-3226
This contribution is concerned with a coupling approach for nonconforming NURBS patches in the framework of an isogeometric formulation for solids in boundary representation. The boundary representation modeling technique in CAD is the starting point of this approach. We parameterize the solid according to the scaled boundary finite element method and employ NURBS basis functions for the approximation of the solution. Therefore, solid surfaces consist of several sections, which can be regarded as patches and discretized independently. The main objective of this study is to derive an approach for the connection of independent sections in order to allow for local refinement and thus an accurate and efficient discretization of the computational domain. Nonconforming sections are coupled with a mortar approach within a master-slave framework. The coupling of adjacent sections ensures the equality of mutual deformations along the interface in a weak sense and is enforced by constraining the NURBS basis functions on the interface. We apply this approach to nonlinear problems in two dimensions and compare the results with conforming discretizations. 相似文献
5.
S. Li J. Trevelyan W. Zhang D. Wang 《International journal for numerical methods in engineering》2018,114(9):975-998
6.
《International journal for numerical methods in engineering》2018,114(11):1213-1227
In this paper, an approach to blend the Hybrid‐Trefftz Finite Element Method (HTFEM) and the Isogeometric Analysis (IGA) called the Isogeometric Trefftz (IGAT) method is presented. The structure of the isogeometric extension of the Trefftz method is formally the same as for its conventional counterpart, except the approximation of the boundary displacements and geometry that are carried out using the Non‐Uniform Rational B‐Splines (NURBS) instead of polynomials. In other words, only the element boundaries are approximated using NURBS basis while the Trefftz approximation is used in the interior of the elements. For that reason, IGAT can be ranked alongside recently developed Isogeometric Boundary Element Method (IGABEM), the NURBS‐Enhanced Finite Element Method (NEFEM), the Isogeometric Local Maximum Entropy (IGA‐LME) method, and the Isogeometrically enhanced Scaled‐Boundary element method (SBFEM), which all use NURBS approximation at the domain boundary only. Theoretical conjectures made in this paper are accompanied by three examples that show that IGAT leads to excellent results using only a few elements. 相似文献
7.
M. Ruess D. Schillinger Y. Bazilevs V. Varduhn E. Rank 《International journal for numerical methods in engineering》2013,95(10):811-846
Enforcing essential boundary conditions plays a central role in immersed boundary methods. Nitsche's idea has proven to be a reliable concept to satisfy weakly boundary and interface constraints. We formulate an extension of Nitsche's method for elasticity problems in the framework of higher order and higher continuity approximation schemes such as the B‐spline and non‐uniform rational basis spline version of the finite cell method or isogeometric analysis on trimmed geometries. Furthermore, we illustrate a significant improvement of the flexibility and applicability of this extension in the modeling process of complex 3D geometries. With several benchmark problems, we demonstrate the overall good convergence behavior of the proposed method and its good accuracy. We provide extensive studies on the stability of the method, its influence parameters and numerical properties, and a rearrangement of the numerical integration concept that in many cases reduces the numerical effort by a factor two. A newly composed boundary integration concept further enhances the modeling process and allows a flexible, discretization‐independent introduction of boundary conditions. Finally, we present our strategy in the framework of the modeling and isogeometric analysis process of trimmed non‐uniform rational basis spline geometries. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
8.
9.
E. De Luycker D. J. Benson T. Belytschko Y. Bazilevs M. C. Hsu 《International journal for numerical methods in engineering》2011,87(6):541-565
The extended finite element method (X‐FEM) has proven to be an accurate, robust method for solving problems in fracture mechanics. X‐FEM has typically been used with elements using linear basis functions, although some work has been performed using quadratics. In the current work, the X‐FEM formulation is incorporated into isogeometric analysis to obtain solutions with higher order convergence rates for problems in linear fracture mechanics. In comparison with X‐FEM with conventional finite elements of equal degree, the NURBS‐based isogeometric analysis gives equal asymptotic convergence rates and equal accuracy with fewer degrees of freedom (DOF). Results for linear through quartic NURBS basis functions are presented for a multiplicity of one or a multiplicity equal the degree. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
10.
Christian Hesch Peter Betsch 《International journal for numerical methods in engineering》2010,82(3):329-358
The present work deals with transient large‐deformation domain decomposition problems. The tying of dissimilar meshed grids is performed by applying the mortar method. In this connection, a reformulation of the original linear mortar constraints is proposed, which retains frame‐indifference for arbitrary discretizations of the interface. Furthermore, a specific coordinate augmentation technique is proposed to make possible the design of an energy–momentum scheme. Numerical examples demonstrate the robustness and enhanced numerical stability of the newly developed energy–momentum scheme for three‐dimensional problems. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
11.
Reissner–Mindlin shell implementation and energy conserving isogeometric multi‐patch coupling
下载免费PDF全文

K. Sommerwerk M. Woidt M. C. Haupt P. Horst 《International journal for numerical methods in engineering》2017,109(7):982-1012
A shear‐flexible isogeometric Reissner–Mindlin shell formulation using non‐uniform rational B‐splines basis functions is introduced, which is used for the demonstration of a coupling approach for multiple non‐conforming patches. The six degrees of freedom formulation uses the exact surface normal vectors and curvature. The shell formulation is implemented in an isogeometric analysis framework for computation of structures composed of multiple geometric entities. To enable local model refinement as well as non‐matching domains coupling, a conservative multi‐patch approach using Lagrange multipliers for structured non‐uniform rational B‐splines patches is presented. Here, an additional border frame mesh is used to couple geometries during structural analyses. This frame interface approach avoids the problem of excessive constraints when multiple patches are coupled at one point. First, the shell formulation is verified with several reference cases. Then the influence of the frame interface discretization and frame penalty stiffness on the smoothness of the results is investigated. The effects of the perturbed Lagrangian method in combination with the frame interface approach is shown. In addition, results of models with T‐joint interface connections and perpendicular stiffener patches are presented. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
12.
Isa Kolo René de Borst 《International journal for numerical methods in engineering》2018,114(4):431-453
Implicit gradient plasticity models incorporate higher‐order spatial gradients via an additional Helmholtz type equation for the plastic multiplier. So far, the enrichment has been limited to second‐order spatial gradients, resulting in a formulation that can be discretised using ‐continuous finite elements. Herein, an implicit gradient plasticity model is formulated that includes a fourth‐order gradient term as well. A comparison between the localisation properties of both the implicit gradient plasticity formulations and the explicit second‐order gradient plasticity model is made using a dispersion analysis. The higher‐order continuity requirement for the fourth‐order implicit gradient plasticity model has been met by exploiting the higher‐order continuity property of isogeometric analysis, which uses nonuniform rational B‐splines as shape functions instead of Lagrange polynomials. The discretised variables, displacements, and plastic multiplier may require different orders of interpolation, an issue that is also addressed. Numerical results show that both formulations can be used as a localisation limiter, but that quantitative differences occur, and a different evolution of the localisation band is obtained for 2‐dimensional problems. 相似文献
13.
Tae‐Kyoung Uhm Sung‐Kie Youn 《International journal for numerical methods in engineering》2009,80(4):507-536
A T‐spline surface is a nonuniform rational B‐spline (NURBS) surface with T‐junctions, and is defined by a control grid called T‐mesh. The T‐mesh is similar to a NURBS control mesh except that in a T‐mesh, a row or column of control points is allowed to terminate in the inner parametric space. This property of T‐splines makes local refinement possible. In the present study, shell formulation based on the T‐spline finite element method (FEM) is presented. Shell formulation based on NURBS or T‐splines has fundamental limitations because rotational DOFs, which are necessary in the shell formulation, cannot be defined on control points. In this study, the simple mapping scheme, in which every control point is mapped into one geometric point on the surface, is employed to eliminate the limitations. Using this mapping scheme, T‐spline FEM can be easily extended to the analysis of shells. The proposed shell formulation is verified through various benchmarking problems. This study is a part of the efforts by the authors for the integration of CAD–CAE processes. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
14.
An isogeometric locking‐free NURBS‐based solid‐shell element for geometrically nonlinear analysis
下载免费PDF全文

Robin Bouclier Thomas Elguedj Alain Combescure 《International journal for numerical methods in engineering》2015,101(10):774-808
In this work, we develop an isogeometric non‐uniform rational B‐spline (NURBS)‐based solid‐shell element for the geometrically nonlinear static analysis of elastic shell structures. A single layer of continuous 3D elements through the thickness of the shell is considered, and the order of approximation in that direction is chosen to be equal to two. A complete 3D constitutive relation is assumed. The objective is to develop a highly accurate low‐order element for coarse meshes. We propose an extension of the mixed method of Bouclier et al. [11] to deal with locking in the context of large rotations and large displacements. The main idea is to modify the interpolation of the average through the thickness of the stress components. It is also necessary to stabilize the element in order to avoid the occurrence of spurious zero‐energy modes. This was achieved, for the quadratic version, through the adjunction of artificial elementary stabilization stiffnesses. The result is an element of order 2, which is at least as accurate as standard NURBS shell elements of order 4. Linear and nonlinear test calculations have been carried out along with comparisons with other published NURBS and classical techniques in order to assess the performance of the element. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
15.
K. B. Nakshatrala K. D. Hjelmstad D. A. Tortorelli 《International journal for numerical methods in engineering》2008,75(12):1385-1415
We present a novel partitioned coupling algorithm to solve first‐order time‐dependent non‐linear problems (e.g. transient heat conduction). The spatial domain is partitioned into a set of totally disconnected subdomains. The continuity conditions at the interface are modeled using a dual Schur formulation where the Lagrange multipliers represent the interface fluxes (or the reaction forces) that are required to maintain the continuity conditions. The interface equations along with the subdomain equations lead to a system of differential algebraic equations (DAEs). For the resulting equations a numerical algorithm is developed, which includes choosing appropriate constraint stabilization techniques. The algorithm first solves for the interface Lagrange multipliers, which are subsequently used to advance the solution in the subdomains. The proposed coupling algorithm enables arbitrary numeric schemes to be coupled with different time steps (i.e. it allows subcycling) in each subdomain. This implies that existing software and numerical techniques can be used to solve each subdomain separately. The coupling algorithm can also be applied to multiple subdomains and is suitable for parallel computers. We present examples showing the feasibility of the proposed coupling algorithm. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
16.
Clemens V. Verhoosel Michael A. Scott Thomas J. R. Hughes René de Borst 《International journal for numerical methods in engineering》2011,86(1):115-134
Continuum damage formulations are commonly used for the simulation of diffuse fracture processes. Implicit gradient damage models are employed to avoid the spurious mesh dependencies associated with local continuum damage models. The C0‐continuity of traditional finite elements has hindered the study of higher order gradient damage approximations. In this contribution we use isogeometric finite elements, which allow for the construction of higher order continuous basis functions on complex domains. We study the suitability of isogeometric finite elements for the discretization of higher order gradient damage approximations. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
17.
D. J. Benson Y. Bazilevs E. De Luycker M.‐C. Hsu M. Scott T. J. R. Hughes T. Belytschko 《International journal for numerical methods in engineering》2010,83(6):765-785
Many of the formulations of current research interest, including iosogeometric methods and the extended finite element method, use nontraditional basis functions. Some, such as subdivision surfaces, may not have convenient analytical representations. The concept of an element, if appropriate at all, no longer coincides with the traditional definition. Developing a new software for each new class of basis functions is a large research burden, especially, if the problems involve large deformations, non‐linear materials, and contact. The objective of this paper is to present a method that separates as much as possible the generation and evaluation of the basis functions from the analysis, resulting in a formulation that can be implemented within the traditional structure of a finite element program but that permits the use of arbitrary sets of basis functions that are defined only through the input file. Elements ranging from a traditional linear four‐node tetrahedron through a higher‐order element combining XFEM and isogeometric analysis may be specified entirely through an input file without any additional programming. Examples of this framework to applications with Lagrange elements, isogeometric elements, and XFEM basis functions for fracture are presented. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
18.
B. Flemisch M. A. Puso B. I. Wohlmuth 《International journal for numerical methods in engineering》2005,63(6):813-832
Dual mortar method formulations have shown to be a very effective and efficient way for interfacing (e.g. tying, contacting) dissimilar meshes. On the other hand, we have recently found that they can sometimes perform quite poorly when applied to curved surfaces in some solid mechanics applications. A new modified two‐dimensional dual mortar method for piecewise linear finite elements is developed that overcomes this deficiency and is demonstrated on a model problem. Furthermore, mathematical analysis is provided to demonstrate the optimal convergence and stability of the new method. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
19.
Peter Otto Laura De Lorenzis Jörg F. Unger 《International journal for numerical methods in engineering》2020,121(6):1248-1267
In this paper, the impact problem and the subsequent wave propagation are considered. For the contact discretization an intermediate non-uniform rational B-spline (NURBS) layer is added between the contacting finite element bodies, which allows a smooth contact formulation and efficient element-based integration. The impact event is ill-posed and requires a regularization to avoid propagating stress oscillations. A nonlinear mesh-dependent penalty regularization is used, where the stiffness of the penalty regularization increases upon mesh refinement. Explicit time integration methods are well suited for wave propagation problems, but are efficient only for diagonal mass matrices. Using a spectral element discretization in combination with a NURBS contact layer the bulk part of the mass matrix is diagonal. 相似文献
20.