Material Non-linearity (material + non-linearity)

Distribution by Scientific Domains


Selected Abstracts


Cyclic macro-element for soil,structure interaction: material and geometrical non-linearities

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 13 2001
Cécile Cremer
Abstract This paper presents a non-linear soil,structure interaction (SSI) macro-element for shallow foundation on cohesive soil. The element describes the behaviour in the near field of the foundation under cyclic loading, reproducing the material non-linearities of the soil under the foundation (yielding) as well as the geometrical non-linearities (uplift) at the soil,structure interface. The overall behaviour in the soil and at the interface is reduced to its action on the foundation. The macro-element consists of a non-linear joint element, expressed in generalised variables, i.e. in forces applied to the foundation and in the corresponding displacements. Failure is described by the interaction diagram of the ultimate bearing capacity of the foundation under combined loads. Mechanisms of yielding and uplift are modelled through a global, coupled plasticity,uplift model. The cyclic model is dedicated to modelling the dynamic response of structures subjected to seismic action. Thus, it is especially suited to combined loading developed during this kind of motion. Comparisons of cyclic results obtained from the macro-element and from a FE modelization are shown in order to demonstrate the relevance of the proposed model and its predictive ability. Copyright © 2001 John Wiley & Sons, Ltd. [source]


A geometrically and materially non-linear piezoelectric three-dimensional-beam finite element formulation including warping effects

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 5 2008
A. Butz
Abstract This paper is concerned with a three-dimensional piezoelectric beam formulation and its finite element implementation. The developed model considers geometrically and materially non-linear effects. An eccentric beam formulation is derived based on the Timoshenko kinematics. The kinematic assumptions are extended by three additional warping functions of the cross section. These functions follow from torsion and piezoelectrically induced shear deformations. The presented beam formulation incorporates large displacements and finite rotations and allows the investigation of stability problems. The finite element model has two nodes with nine mechanical and five electrical degrees of freedom. It provides an accurate approximation of the electric potential, which is assumed to be linear in the direction of the beam axis and quadratic within the cross section. The mechanical degrees of freedom are three displacements, three rotations and three scaling factors for the warping functions. The latter are computed in a preprocess by solving a two-dimensional in-plane equilibrium condition with the finite element method. The gained warping patterns are considered within the integration through the cross section of the beam formulation. With respect to material non-linearities, which arise in ferroelectric materials, the scalar Preisach model is embedded in the formulation. This model is a mathematical model for the general description of hysteresis phenomena. Its application to piezoelectric materials leads to a phenomenological model for ferroelectric hysteresis effects. Here, the polarization direction is assumed to be constant, which leads to unidirectional constitutive equations. Some examples demonstrate the capability of the proposed model. Copyright © 2008 John Wiley & Sons, Ltd. [source]


Performance-based seismic analysis and design of suspension bridges

EARTHQUAKE ENGINEERING AND STRUCTURAL DYNAMICS, Issue 4-5 2005
Serafim Arzoumanidis
Abstract This paper presents a performance-based seismic analysis and design of a large suspension bridge, the new Tacoma Narrows Parallel Crossing in the State of Washington. The scope of the project included establishment of design criteria, extensive analysis and validation of the design. The analysis was performed using detailed three-dimensional models that included geometric and material non-linearity. The target post-earthquake level of service was verified using stress, deformation and ductility criteria. In the absence of well-established criteria, which relate the structural response of tower shafts to specific levels of performance, capacity analyses were performed to demonstrate that the design fulfills the performance objectives. The seismic analysis and design of this bridge was reviewed throughout the design process. An independent check team also performed separate analysis and validation of the design. Thus, this bridge constitutes an example of a large-scale design project where the performance-based seismic design procedures underwent rigorous assessment. This work demonstrated that the performance-based approach for seismic design is an appropriate way for designing earthquake-resistant structures. Further data that relate the structural response with the performance objectives are necessary. Copyright © 2005 John Wiley & Sons, Ltd. [source]


An operator-split ALE model for large deformation analysis of geomaterials

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 12 2007
Y. Di
Abstract Analysis of large deformation of geomaterials subjected to time-varying load poses a very difficult problem for the geotechnical profession. Conventional finite element schemes using the updated Lagrangian formulation may suffer from serious numerical difficulties when the deformation of geomaterials is significantly large such that the discretized elements are severely distorted. In this paper, an operator-split arbitrary Lagrangian,Eulerian (ALE) finite element model is proposed for large deformation analysis of a soil mass subjected to either static or dynamic loading, where the soil is modelled as a saturated porous material with solid,fluid coupling and strong material non-linearity. Each time step of the operator-split ALE algorithm consists of a Lagrangian step and an Eulerian step. In the Lagrangian step, the equilibrium equation and continuity equation of the saturated soil are solved by the updated Lagrangian method. In the Eulerian step, mesh smoothing is performed for the deformed body and the state variables obtained in the updated Lagrangian step are then transferred to the new mesh system. The accuracy and efficiency of the proposed ALE method are verified by comparison of its results with the results produced by an analytical solution for one-dimensional finite elastic consolidation of a soil column and with the results from the small strain finite element analysis and the updated Lagrangian analysis. Its performance is further illustrated by simulation of a complex problem involving the transient response of an embankment subjected to earthquake loading. Copyright © 2007 John Wiley & Sons, Ltd. [source]


Displacement-controlled method and its applications to material non-linearity

INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS, Issue 3 2005
H. Zheng
Abstract For the analysis of non-linear problems, the displacement-controlled method (DCM) has a more extensive application scope and more powerful abilities than the load-controlled method (LCM). However, difficulties of the DCM's procedure not amenable to most finite element implementations of the conventional LCM have restricted its applications in geomechanics. By means of Sherman,Morrison's theorem, the solution of DCM is improved. The improved procedure is characterized by high efficiency, good numerical stability and a programme structure similar to LCM. Two aspects of applications of DCM are illustrated. The first application is to compute the response of a structure under a given load level like the conventional finite element analysis. The second application is to trace the equilibrium path of a structure under a given load distribution type. A simple but effective algorithm is presented for automatically adjusting the step length in tracing the equilibrium path. Examples illustrate that the proposed procedures are suited for modelling complicated non-linear problems in geomechanics. Copyright © 2005 John Wiley & Sons, Ltd. [source]


A new triangular layered plate element for the non-linear analysis of reinforced concrete slabs

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING, Issue 7 2006
Y. X. Zhang
Abstract A new 3-node, 18-DOF triangular layered plate element is developed in this paper for the geometric and material non-linear analysis of isotropic plates and reinforced concrete slabs under service loads. The proposed model is a combination of Allman's 3-node, 9-DOF triangular membrane element with drilling degrees of freedom and the refined non-conforming 3-node, 9-DOF triangular plate-bending element RT9 in order to account for the coupling effects between membrane and bending actions. The element is modelled as a layered system of concrete and equivalent smeared steel reinforcement layers, and perfect bond is assumed between the concrete layers and the smeared steel layers. The maximum normal stress criterion is employed to detect cracking of the concrete, and a smeared fixed crack model is assumed. Both geometric non-linearity with large displacements but moderate rotations and material non-linearity, which incorporates tension, compression, concrete cracking and tension stiffening, are included in the model. An updated Lagrangian approach is employed as a solution strategy for the non-linear finite element analysis and a numerical example of reinforced concrete slab is given to demonstrate the efficacy of this robust element. Copyright © 2005 John Wiley & Sons, Ltd. [source]


Adaptive superposition of finite element meshes in non-linear transient solid mechanics problems

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 9 2007
Z. Yue
Abstract An s-adaptive finite element procedure is developed for the transient analysis of 2-D solid mechanics problems with material non-linearity due to progressive damage. The resulting adaptive method simultaneously estimates and controls both the spatial error and temporal error within user-specified tolerances. The spatial error is quantified by the Zienkiewicz,Zhu error estimator and computed via superconvergent patch recovery, while the estimation of temporal error is based on the assumption of a linearly varying third-order time derivatives of the displacement field in conjunction with direct numerical time integration. The distinguishing characteristic of the s-adaptive procedure is the use of finite element mesh superposition (s-refinement) to provide spatial adaptivity. Mesh superposition proves to be particularly advantageous in computationally demanding non-linear transient problems since it is faster, simpler and more efficient than traditional h-refinement schemes. Numerical examples are provided to demonstrate the performance characteristics of the s-adaptive method for quasi-static and transient problems with material non-linearity. Copyright © 2007 John Wiley & Sons, Ltd. [source]


On the design of energy,momentum integration schemes for arbitrary continuum formulations.

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 15 2004
Applications to classical, chaotic motion of shells
Abstract The construction of energy,momentum methods depends heavily on three kinds of non-linearities: (1) the geometric (non-linearity of the strain,displacement relation), (2) the material (non-linearity of the elastic constitutive law), and (3) the one exhibited in displacement-dependent loading. In previous works, the authors have developed a general method which is valid for any kind of geometric non-linearity. In this paper, we extend the method and combine it with a treatment of material non-linearity as well as that exhibited in force terms. In addition, the dynamical formulation is presented in a general finite element framework where enhanced strains are incorporated as well. The non-linearity of the constitutive law necessitates a new treatment of the enhanced strains in order to retain the energy conservation property. Use is made of the logarithmic strain tensor which allows for a highly non-linear material law, while preserving the advantage of considering non-linear vibrations of classical metallic structures. Various examples and applications to classical and non-classical vibrations and non-linear motion of shells are presented, including (1) chaotic motion of arches, cylinders and caps using a linear constitutive law and (2) large overall motion and non-linear vibration of shells using non-linear constitutive law. Copyright © 2004 John Wiley & Sons, Ltd. [source]


c-Type method of unified CAMG and FEA.

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 9 2003
2D non-linear, 3D linear, Part 1: Beam, arch mega-elements
Abstract Computer-aided mesh generation (CAMG) dictated solely by the minimal key set of requirements of geometry, material, loading and support condition can produce ,mega-sized', arbitrary-shaped distorted elements. However, this may result in substantial cost saving and reduced bookkeeping for the subsequent finite element analysis (FEA) and reduced engineering manpower requirement for final quality assurance. A method, denoted as c-type, has been proposed by constructively defining a finite element space whereby the above hurdles may be overcome with a minimal number of hyper-sized elements. Bezier (and de Boor) control vectors are used as the generalized displacements and the Bernstein polynomials (and B-splines) as the elemental basis functions. A concomitant idea of coerced parametry and inter-element continuity on demand unifies modelling and finite element method. The c-type method may introduce additional control, namely, an inter-element continuity condition to the existing h-type and p-type methods. Adaptation of the c-type method to existing commercial and general-purpose computer programs based on a conventional displacement-based finite element method is straightforward. The c-type method with associated subdivision technique can be easily made into a hierarchic adaptive computer method with a suitable a posteriori error analysis. In this context, a summary of a geometrically exact non-linear formulation for the two-dimensional curved beams/arches is presented. Several beam problems ranging from truly three-dimensional tortuous linear curved beams to geometrically extremely non-linear two-dimensional arches are solved to establish numerical efficiency of the method. Incremental Lagrangian curvilinear formulation may be extended to overcome rotational singularity in 3D geometric non-linearity and to treat general material non-linearity. Copyright © 2003 John Wiley & Sons, Ltd. [source]


A vertex-based finite volume method applied to non-linear material problems in computational solid mechanics

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 4 2003
G. A. Taylor
Abstract A vertex-based finite volume (FV) method is presented for the computational solution of quasi-static solid mechanics problems involving material non-linearity and infinitesimal strains. The problems are analysed numerically with fully unstructured meshes that consist of a variety of two- and three-dimensional element types. A detailed comparison between the vertex-based FV and the standard Galerkin FE methods is provided with regard to discretization, solution accuracy and computational efficiency. For some problem classes a direct equivalence of the two methods is demonstrated, both theoretically and numerically. However, for other problems some interesting advantages and disadvantages of the FV formulation over the Galerkin FE method are highlighted. Copyright © 2002 John Wiley & Sons, Ltd. [source]