共查询到20条相似文献，搜索用时 0 毫秒
1.
Ran Ma;WaiChing Sun;Tong Guo; 《International journal for numerical methods in engineering》2024,125(13):e7462
We present a phase-field fracture model for a stress resultant geometrically exact shell in finite deformation regime where the configuration manifold evolves according to deformation and fracture. The Reissner–Mindlin shell problem is first solved via the finite element method, where the independent unknown fields are the displacement and director. The phase-field ductile fracture model is then coupled with the verified geometrically exact shell by enriching the three-field elasto-plastic free energy with a regularized crack surface energy. To capture the geometrical nonlinearity due to the large plastic deformation exhibited in the processing zone, the corresponding finite-strain stress resultant elasto-plasticity model is coupled with the phase field model. This coupling between the stress resultant elasto-plasticity model and the phase-field fracture model enables us to predict the different amounts of energy dissipation attributed to plastic deformation and fracture and hence simulate the crack propagation in the ductile regime properly. We introduce a mixed finite element model with displacement, director, and phase-field order parameters as the nodal degree of freedoms formulated on the mid-surface. An energy-based arc-length method is generalized to track the equilibrium path and mitigate instabilities arising from material nonlinearity. Four numerical examples are presented to validate the implementation of the model and demonstrate its capability to simulate ductile fracture in shell structures. These examples include a plane-stress tension/shear fracture model, the Muscat-Fenech and Atkins plate, axial tension of a notched cylindrical shell, and ductile fracture of a simply supported plate under uniform pressure. 相似文献
2.
初始缺陷是影响网壳稳定性的主要因素之一。技术规程建议的一致缺陷模态法综合考虑了结点的位置偏差与杆件的初弯曲,但现有计算方法难于便捷地单独全面体现初弯曲的影响。该文基于几何精确梁理论构造了一种弱形式求积元模型。通过引入纤维模型模拟材料非线性,结合柱面弧长法实现了对空间曲梁结构的弹塑性大位移分析,通过算例验证了提出模型的有效性。该模型构建不需具体指定位移形函数,避免了使用有限元软件中的复杂建模过程。运用该模型计算了几类典型网壳在杆件初弯曲方向随机分布,不同弯曲挠度下的极限承载力。计算结果表明,在现有钢结构生产工艺下杆件初弯曲的缺陷对于网壳稳定性的影响较小,不起控制作用。通过与一致缺陷模态法的计算结果和空间网格结构技术规程计算的承载力进行对比,对规范的承载力计算公式提出了考虑特殊网壳形式及材料屈服强度的改进建议。 相似文献
3.
Jia Lu Xianlian Zhou Madhavan L. Raghavan 《International journal for numerical methods in engineering》2007,69(6):1239-1261
The paper presents a computational method for predicting the initial geometry of a finitely deforming anisotropic elastic body from a given deformed state. The method is imperative for a class of problem in stress analysis, particularly in biomechanical applications. While the basic idea has been established elsewhere Comput. Methods Appl. Mech. Eng. 1996; 136 :47–57; Int. J. Numer. Meth. Engng 1998; 43 : 821–838), the implementation in general anisotropic solids is not a trivial exercise, but comes after a systematic development of Eulerian representations of constitutive equations. In this paper, we discuss the general representation in the context of fibrous hyperelastic solids, and provide explicit stress functions for some commonly used soft tissue models including the Fung model and the Holzapfel model. A three‐field mixed formulation is introduced to enforce quasi‐incompressibility constraints. The practical utility of this method is demonstrated using an example of aneurysm stress analysis. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
4.
Murillo V. B. Santana Carlo Sansour Mohammed Hjiaj Hugues Somja 《International journal for numerical methods in engineering》2022,123(2):444-464
This article describes a novel equilibrium-based geometrically exact beam finite element formulation. First, the spatial position and rotation fields are interpolated by nonlinear configuration-dependent functions that enforce constant strains along the element axis, completely eliminating locking phenomena. Then, the resulting kinematic fields are used to interpolate the spatial sections force and moment fields in order to fulfill equilibrium exactly in the deformed configuration. The internal variables are explicitly solved at the element level and closed-form expressions for the internal force vector and tangent stiffness matrix are obtained, allowing for explicit computation, without numerical integration. The objectivity and absence of locking are verified and some important numerical and theoretical aspects leading to a computationally efficient strategy are highlighted and discussed. The proposed formulation is successfully tested in several numerical application examples. 相似文献
5.
6.
R. Hauptmann K. Schweizerhof 《International journal for numerical methods in engineering》1998,42(1):49-69
In the present contribution we propose a so-called solid-shell concept which incorporates only displacement degrees of freedom. Thus, some major disadvantages of the usually used degenerated shell concept are overcome. These disadvantages are related to boundary conditions—the handling of soft and hard support, the need for special co-ordinate systems at boundaries, the connection with continuum elements—and, in geometrically non-linear analyses, to a complicated update of the rotation vector. First, the kinematics of the so-called solid-shell concept in analogy to the degenerated shell concept are introduced. Then several modifications of the solid-shell concept are proposed to obtain locking-free solid-shell elements, leading also to formulations which allow the use of general three-dimensional material laws and which are also able to represent the normal stresses and strains in thickness direction. Numerical analyses of geometrically linear and non-linear problems are finally performed using solely assumed natural shear strain elements with a linear approximation in in-plane direction. Although some considerations are needed to get comparable boundary conditions in the examples analysed, the solid-shell elements prove to work as good as the degenerated shell elements. The numerical examples show that neither thickness nor shear locking are present even for distorted element shapes. © 1998 John Wiley & Sons, Ltd. 相似文献
7.
B. Skallerud B. Haugen 《International journal for numerical methods in engineering》1999,46(12):1961-1986
Due to the very non‐linear behaviour of thin shells under collapse, numerical simulations are subject to challenges. Shell finite elements are attractive in these simulations. Rotational degrees of freedom do, however, complicate the solution. In the present study a co‐rotated formulation is employed. The deformation of the shell is decomposed in to a contribution from large rigid body rotation and a strain producing term. A triangular assumed strain shell finite element is used. Hence, a high performance elastic element is combined with the co‐rotated formulation. In the co‐rotated co‐ordinate system the plasticity is accounted for by a simplifyed Ilyushin stress resultant yield surface. The stress update is determined from the backward Euler difference, and a consistent geometrical and material tangent stiffness is derived. Comparison with other published analysis results show that the present formulation gives acceptable accuracy. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
8.
G. M. Kulikov S. V. Plotnikova A. O. Glebov 《International journal for numerical methods in engineering》2020,121(17):3795-3823
In this work, the finite rotation exact geometry four-node solid-shell element using the sampling surfaces (SaS) method is developed for the analysis of the second Piola-Kirchhoff stresses in laminated piezoelectric shells. The SaS method is based on choosing inside the layers the arbitrary number of SaS parallel to the middle surface and located at Chebyshev polynomial nodes in order to introduce the displacements and electric potentials of these surfaces as fundamental shell unknowns. The outer surfaces and interfaces are also included into a set of SaS. To circumvent shear and membrane locking, the hybrid-mixed solid-shell element on the basis of the Hu-Washizu variational principle is proposed. The tangent stiffness matrix is evaluated by 3D analytical integration throughout the finite element. This novelty provides a superior performance in the case of coarse meshes. A comparison with the SOLID226 element showed that the developed exact geometry SaS solid-shell element allows the use of load increments, which are much larger than possible with existing displacement-based finite elements. Thus, it can be recommended for the 3D stress analysis of doubly-curved laminated piezoelectric shells because the SaS formulation gives the opportunity to obtain the 3D solutions of electroelasticity with a prescribed accuracy. 相似文献
9.
This paper presents a robust non-linear piezoelectric exact geometry (EG) four-node solid-shell element based on the higher-order 9-parameter equivalent single-layer (ESL) theory, which permits one to utilize 3D constitutive equations. The term EG reflects the fact that coefficients of the first and second fundamental forms of the reference surface are taken exactly at each element node. The finite element formulation developed is based on a new concept of interpolation surfaces (I-surfaces) inside the shell body. We introduce three I-surfaces and choose nine displacements of these surfaces as fundamental shell unknowns. Such choice allows us to represent the finite rotation piezoelectric higher-order EG solid-shell element formulation in a very compact form and to utilize in curvilinear reference surface coordinates the strain-displacement relationships, which are objective, that is, invariant under arbitrarily large rigid-body shell motions. To avoid shear and membrane locking and have no spurious zero energy modes, the assumed displacement-independent strain and stress resultant fields are introduced. In this connection, the Hu-Washizu variational equation is invoked. To implement the analytical integration throughout the element, the modified ANS method is applied. As a result, the present finite rotation piezoelectric EG solid-shell element formulation permits the use of coarse meshes and very large load increments. 相似文献
10.
G. M. Kulikov S. V. Plotnikova 《International journal for numerical methods in engineering》2011,88(13):1363-1389
This paper presents the finite rotation exact geometry (EG) 12‐node solid‐shell element with 36 displacement degrees of freedom. The term ‘EG’ reflects the fact that coefficients of the first and second fundamental forms of the reference surface and Christoffel symbols are taken exactly at each element node. The finite element formulation developed is based on the 9‐parameter shell model by employing a new concept of sampling surfaces (S‐surfaces) inside the shell body. We introduce three S‐surfaces, namely, bottom, middle and top, and choose nine displacements of these surfaces as fundamental shell unknowns. Such choice allows one to represent the finite rotation higher order EG solid‐shell element formulation in a very compact form and to derive the strain–displacement relationships, which are objective, that is, invariant under arbitrarily large rigid‐body shell motions in convected curvilinear coordinates. The tangent stiffness matrix is evaluated by using 3D analytical integration and the explicit presentation of this matrix is given. The latter is unusual for the non‐linear EG shell element formulation. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
11.
Chen Wanji Zheng Shijie 《International journal for numerical methods in engineering》1998,41(7):1195-1213
Based on a variational principle with relaxed inter-element continuity requirements, a refined hybrid quadrilateral degenerated shell element GNRH6, which is a non-conforming model with six internal displacements, is proposed for the geometrically non-linear analysis. The orthogonal approach and non-conforming modes are incorporated into the geometrically non-linear formulation. Numerical results show that the orthogonal approach can improve computational efficiency while the non-conforming modes can eliminate the shear/membrane locking phenomenon and improve the accuracy. © 1998 John Wiley & Sons, Ltd. 相似文献
12.
13.
几何非线性高性能复合材料筋混凝土梁Heterosis组合壳单元 总被引：2，自引：0，他引：2
对于高性能碳纤维增强聚合物复合材料(CFRP)筋混凝土梁,研究几何非线性组合壳单元模型,对预应力CFRP筋混凝土梁进行了全过程分析。引入Von Karman理论,推导了局部坐标系下Piola2Kirchhoff 应力矩阵和几何刚度矩阵;分别采用组合壳单元和分层壳单元模拟预应力CFRP 筋和玻璃纤维增强聚合物复合材料(GFRP)筋,并推导了CFRP筋对组合壳单元刚度矩阵的贡献,同时采用Heterosis选择积分技术以避免剪切锁定和零能量模式,研制了相应的非线性计算程序。计算结果与试验数据对比可知,挠度发展规律和预应力CFRP筋应变发展规律均吻合良好,说明了研究单元的有效性及研制程序的正确性;CFRP筋具有高强度性能,梁试件破坏时CFRP筋均未失效;利用预应力CFRP筋应变重分布系数研究了梁的刚度退化规律,表明采用GFRP筋代替普通钢筋在加载后期会使梁的刚度退化减小。 相似文献
14.
I. Romero F. Armero 《International journal for numerical methods in engineering》2002,54(7):1043-1086
This paper presents a new family of time‐stepping algorithms for the integration of the dynamics of non‐linear shells. We consider the geometrically exact shell theory involving an inextensible director field (the so‐called five‐parameter shell model). The main characteristic of this model is the presence of the group of finite rotations in the configuration manifold describing the deformation of the solid. In this context, we develop time‐stepping algorithms whose discrete solutions exhibit the same conservation laws of linear and angular momenta as the underlying physical system, and allow the introduction of a controllable non‐negative energy dissipation to handle the high numerical stiffness characteristic of these problems. A series of algorithmic parameters for the different components of the deformation of the shell (i.e. membrane, bending and transverse shear) fully control this numerical dissipation, recovering existing energy‐momentum schemes as a particular choice of these algorithmic parameters. We present rigorous proofs of the numerical properties of the resulting algorithms in the full non‐linear range. Furthermore, it is argued that the numerical dissipation is introduced in the high‐frequency range by considering the proposed algorithm in the context of a linear problem. The finite element implementation of the resulting methods is described in detail as well as considered in the final arguments proving the aforementioned conservation/dissipation properties. We present several representative numerical simulations illustrating the performance of the newly proposed methods. The robustness gained over existing methods in these stiff problems is confirmed in particular. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
15.
Ezgi Günay 《Applied Composite Materials》1999,6(6):381-395
A finite-element solution for a geometrically nonlinear fiber composite shallow shell with Stringer-type stiffeners is presented. A laminated, anisotropic thin/thick shallow stiffened shell finite element is developed and applied for the solution of several static problems. Geometrically, a nonlinear finite-element model is based on the nonlocking shear deformable theory. Stiffened composite cylinders are subjected to mechanical loading and the shallow shell theory is used for the geometric representation and formulation of the axially stiffened cylinders. A new, two-sided meshing system is generated to represent a cylindrical shell with stiffeners in a three-dimensional coordinate system. 相似文献
16.
17.
S. A. Hosseini Kordkheili R. Naghdabadi 《International journal for numerical methods in engineering》2007,72(8):964-986
A finite element formulation governing the geometrically non‐linear thermoelastic behaviour of plates and shells made of functionally graded materials is derived in this paper using the updated Lagrangian approach. Derivation of the formulation is based on rewriting the Green–Lagrange strain as well as the 2nd Piola–Kirchhoff stress as two second‐order functions in terms of a through‐the‐thickness parameter. Material properties are assumed to vary through the thickness according to the commonly used power law distribution of the volume fraction of the constituents. Within a non‐linear finite element analysis framework, the main focus of the paper is the proposal of a formulation to account for non‐linear stress distribution in FG plates and shells, particularly, near the inner and outer surfaces for small and large values of the grading index parameter. The non‐linear heat transfer equation is also solved for thermal distribution through the thickness by the Rayleigh–Ritz method. Advantages of the proposed approach are assessed and comparisons with available solutions are presented. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
18.
Frdric Boyer Dominique Primault 《International journal for numerical methods in engineering》2004,59(5):669-702
This article is devoted to the modelling of thin beams undergoing finite deformations essentially due to bending and torsion and to their numerical resolution by the finite element method. The solution proposed here differs from the approaches usually implemented to treat thin beams, as it can be qualified as ‘geometrically exact’. Two numerical models are proposed. The first one is a non‐linear Euler–Bernoulli model while the second one is a non‐linear Rayleigh model. The finite element method is tested on several numerical examples in statics and dynamics, and validated through comparison with analytical solutions, experimental observations and the geometrically exact approach of the Reissner beam theory initiated by Simo. The numerical result shows that this approach is a good alternative to the modelling of non‐linear beams, especially in statics. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
19.
An isogeometric locking‐free NURBS‐based solid‐shell element for geometrically nonlinear analysis 
下载免费PDF全文
下载免费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. 相似文献
20.
Adnan Ibrahimbegovi&#x; Mazen Al Mikdad 《International journal for numerical methods in engineering》1998,41(5):781-814
We examine theoretical and computational aspects of three-dimensional finite rotations pertinent to the dynamics of beams. The model problem chosen for consideration is the Reissner beam theory capable of modelling finite strains and finite rotations in geometrically exact manner. Special emphasis is placed on clarifying the geometry aspects, finite rotation updates and the associated linearization procedure pertaining to different choices of rotation parameters. The latter is shown to play an important role in constructing the optimal implementation of a time-stepping scheme. © 1998 John Wiley & Sons, Ltd. 相似文献