Applied element method

Last updated January 17, 2023

The applied element method (AEM) is a numerical analysis used in predicting the continuum and discrete behavior of structures. The modeling method in AEM adopts the concept of discrete cracking allowing it to automatically track structural collapse behavior passing through all stages of loading: elastic, crack initiation and propagation in tension-weak materials, reinforcement yield, element separation, element contact and collision, as well as collision with the ground and adjacent structures.

History

Exploration of the approach employed in the applied element method began in 1995 at the University of Tokyo as part of Dr. Hatem Tagel-Din's research studies. The term "applied element method" itself, however, was first coined in 2000 in a paper called "Applied element method for structural analysis: Theory and application for linear materials".^[1] Since then AEM has been the subject of research by a number of academic institutions and the driving factor in real-world applications. Research has verified its accuracy for: elastic analysis;^[1] crack initiation and propagation; estimation of failure loads at reinforced concrete structures;^[2] reinforced concrete structures under cyclic loading;^[3] buckling and post-buckling behavior;^[4] nonlinear dynamic analysis of structures subjected to severe earthquakes;^[5] fault-rupture propagation;^[6] nonlinear behavior of brick structures;^[7] and the analysis of glass reinforced polymers (GFRP) walls under blast loads.^[8]

Technical discussion

In AEM, the structure is divided virtually and modeled as an assemblage of relatively small elements. The elements are then connected through a set of normal and shear springs located at contact points distributed along with the element faces. Normal and shear springs are responsible for the transfer of normal and shear stresses from one element to the next.

Element generation and formulation

The modeling of objects in AEM is very similar to modeling objects in FEM. Each object is divided into a series of elements connected and forming a mesh. The main difference between AEM and FEM, however, is how the elements are joined together. In AEM the elements are connected by a series of non-linear springs representing the material behavior.

There are three types of springs used in AEM:

Matrix Springs: Matrix springs connect two elements together representing the main material properties of the object.
Reinforcing Bar Springs: Reinforcement springs are used to implicitly represent additional reinforcement bars running through the object without adding additional elements to the analysis.
Contact Springs: Contact Springs are generated when two elements collide with each other or the ground. When this occurs three springs are generated (Shear Y, Shear X and Normal).

Automatic element separation

When the average strain value at the element face reaches the separation strain, all springs at this face are removed and elements are no longer connected until a collision occurs, at which point they collide together as rigid bodies.

Separation strain represents the strain at which adjacent elements are totally separated at the connecting face. This parameter is not available in the elastic material model. For concrete, all springs between the adjacent faces including reinforcement bar springs are cut. If the elements meet again, they will behave as two different rigid bodies that have now contacted each other. For steel, the bars are cut if the stress point reaches ultimate stress or if the concrete reaches the separation strain.

Automatic element contact/collision

Contact or collision is detected without any user intervention. Elements are able to separate, contract and/or make contact with other elements. In AEM three contact methods include Corner-to-Face, Edge-to-Edge, and Corner-to-Ground.

Stiffness matrix

The spring stiffness in a 2D model can be calculated from the following equations:

K_{n}={\frac {E\cdot T\cdot d}{a}}

K_{s}={\frac {G\cdot T\cdot d}{a}}

Where d is the distance between springs, T is the thickness of the element, a is the length of the representative area, E is the Young's modulus, and G is the shear modulus of the material. The above equation's indicate that each spring represents the stiffness of an area (T·d) within the length of the studied material.

To model reinforcement bars embedded in concrete, a spring is placed inside the element at the location of the bar; the area (T·d) is replaced by the actual cross section area of the reinforcement bar. Similar to modeling embedded steel sections, the area (T·d) may be replaced by the area of the steel section represented by the spring.

Although the element motion moves as a rigid body, its internal deformations are represented by the spring deformation around each element. This means the element shape does not change during analysis, but the behavior of assembly of elements is deformable. The two elements are assumed to be connected by only one pair of normal and shear springs. To have a general stiffness matrix, the locations of element and contact springs are assumed in a general position. The stiffness matrix components corresponding to each degree of freedom are determined by assuming a unit displacement in the studied direction and by determining forces at the centroid of each element. The 2D element stiffness matrix size is 6 × 6; the components of the upper left quarter of the stiffness matrix are shown below:

{\begin{bmatrix}\sin ^{2}(\theta +\alpha )K_{n}&-K_{n}\sin(\theta +\alpha )\cos(\theta +\alpha )&\cos(\theta +\alpha )K_{s}L\sin(\alpha )\\+\cos ^{2}(\theta +\alpha )K_{s}&+K_{s}\sin(\theta +\alpha )\cos(\theta +\alpha )&-\sin(\theta +\alpha )K_{n}L\cos(\alpha )\\\\-K_{n}\sin(\theta +\alpha )\cos(\theta +\alpha )&\sin ^{2}(\theta +\alpha )K_{s}&\cos(\theta +\alpha )K_{n}L\cos(\alpha )\\+K_{s}\sin(\theta +\alpha )\cos(\theta +\alpha )&+\cos ^{2}(\theta +\alpha )K_{n}&+\sin(\theta +\alpha )K_{s}L\sin(\alpha )\\\\\cos(\theta +\alpha )K_{s}L\sin(\alpha )&\cos(\theta +\alpha )K_{n}L\cos(\alpha )&L^{2}\cos ^{2}(\alpha )K_{n}\\-\sin(\theta +\alpha )K_{n}L\cos(\alpha )&+\sin(\theta +\alpha )K_{s}L\sin(\alpha )&+L^{2}\sin ^{2}(\alpha )K_{s}\end{bmatrix}}

The stiffness matrix depends on the contact spring stiffness and the spring location. The stiffness matrix is for only one pair of contact springs. However, the global stiffness matrix is determined by summing up the stiffness matrices of individual pairs of springs around each element. Consequently, the developed stiffness matrix has total effects from all pairs of springs, according to the stress situation around the element. This technique can be used in both load and displacement control cases. The 3D stiffness matrix may be deduced similarly.

Applications

The applied element method is currently being used in the following applications:

Structural vulnerability assessment
- Progressive collapse
- Blast analysis
- Impact analysis
- Seismic analysis
Forensic engineering
Performance based design
Demolition analysis
Glass performance analysis
Visual effects

Related Research Articles

In optics, polarized light can be described using the Jones calculus, discovered by R. C. Jones in 1941. Polarized light is represented by a Jones vector, and linear optical elements are represented by Jones matrices. When light crosses an optical element the resulting polarization of the emerging light is found by taking the product of the Jones matrix of the optical element and the Jones vector of the incident light. Note that Jones calculus is only applicable to light that is already fully polarized. Light which is randomly polarized, partially polarized, or incoherent must be treated using Mueller calculus.

In mathematical physics and mathematics, the Pauli matrices are a set of three $2 \times 2$ complex matrices which are Hermitian, involutory and unitary. Usually indicated by the Greek letter sigma, they are occasionally denoted by tau when used in connection with isospin symmetries.

Kinematics is a subfield of physics, developed in classical mechanics, that describes the motion of points, bodies (objects), and systems of bodies without considering the forces that cause them to move. Kinematics, as a field of study, is often referred to as the "geometry of motion" and is occasionally seen as a branch of mathematics. A kinematics problem begins by describing the geometry of the system and declaring the initial conditions of any known values of position, velocity and/or acceleration of points within the system. Then, using arguments from geometry, the position, velocity and acceleration of any unknown parts of the system can be determined. The study of how forces act on bodies falls within kinetics, not kinematics. For further details, see analytical dynamics.

<span class="mw-page-title-main">Angular displacement</span>

Angular displacement of a body is the angle through which a point revolves around a centre or a specified axis in a specified sense. When a body rotates about its axis, the motion cannot simply be analyzed as a particle, as in circular motion it undergoes a changing velocity and acceleration at any time (t). When dealing with the rotation of a body, it becomes simpler to consider the body itself rigid. A body is generally considered rigid when the separations between all the particles remains constant throughout the body's motion, so for example parts of its mass are not flying off. In a realistic sense, all things can be deformable, however this impact is minimal and negligible. Thus the rotation of a rigid body over a fixed axis is referred to as rotational motion.

In linear algebra, an invertible complex square matrix $U$ is unitary if its conjugate transpose $U *$ is also its inverse, that is, if

A composite material is a material which is produced from two or more constituent materials. These constituent materials have notably dissimilar chemical or physical properties and are merged to create a material with properties unlike the individual elements. Within the finished structure, the individual elements remain separate and distinct, distinguishing composites from mixtures and solid solutions.

In mechanics and geometry, the 3D rotation group, often denoted SO(3), is the group of all rotations about the origin of three-dimensional Euclidean space $under the operation of composition.$

In linear algebra, linear transformations can be represented by matrices. If $is a linear transformation mapping to and is a column vector with entries, then$

In linear algebra, a rotation matrix is a transformation matrix that is used to perform a rotation in Euclidean space. For example, using the convention below, the matrix

In numerical analysis, the Clenshaw algorithm, also called Clenshaw summation, is a recursive method to evaluate a linear combination of Chebyshev polynomials. The method was published by Charles William Clenshaw in 1955. It is a generalization of Horner's method for evaluating a linear combination of monomials.

The contact angle is the angle, conventionally measured through the liquid, where a liquid–vapor interface meets a solid surface. It quantifies the wettability of a solid surface by a liquid via the Young equation. A given system of solid, liquid, and vapor at a given temperature and pressure has a unique equilibrium contact angle. However, in practice a dynamic phenomenon of contact angle hysteresis is often observed, ranging from the advancing (maximal) contact angle to the receding (minimal) contact angle. The equilibrium contact is within those values, and can be calculated from them. The equilibrium contact angle reflects the relative strength of the liquid, solid, and vapour molecular interaction.

In geometry, various formalisms exist to express a rotation in three dimensions as a mathematical transformation. In physics, this concept is applied to classical mechanics where rotational kinematics is the science of quantitative description of a purely rotational motion. The orientation of an object at a given instant is described with the same tools, as it is defined as an imaginary rotation from a reference placement in space, rather than an actually observed rotation from a previous placement in space.

The Wigner D-matrix is a unitary matrix in an irreducible representation of the groups SU(2) and SO(3). It was introduced in 1927 by Eugene Wigner, and plays a fundamental role in the quantum mechanical theory of angular momentum. The complex conjugate of the D-matrix is an eigenfunction of the Hamiltonian of spherical and symmetric rigid rotors. The letter $D$ stands for Darstellung, which means "representation" in German.

In the differential geometry of surfaces, a Darboux frame is a natural moving frame constructed on a surface. It is the analog of the Frenet–Serret frame as applied to surface geometry. A Darboux frame exists at any non-umbilic point of a surface embedded in Euclidean space. It is named after French mathematician Jean Gaston Darboux.

In mechanical engineering, the Denavit–Hartenberg parameters are the four parameters associated with a particular convention for attaching reference frames to the links of a spatial kinematic chain, or robot manipulator.

Slope stability analysis is a static or dynamic, analytical or empirical method to evaluate the stability of earth and rock-fill dams, embankments, excavated slopes, and natural slopes in soil and rock. Slope stability refers to the condition of inclined soil or rock slopes to withstand or undergo movement. The stability condition of slopes is a subject of study and research in soil mechanics, geotechnical engineering and engineering geology. Analyses are generally aimed at understanding the causes of an occurred slope failure, or the factors that can potentially trigger a slope movement, resulting in a landslide, as well as at preventing the initiation of such movement, slowing it down or arresting it through mitigation countermeasures.

Extreme Loading for Structures (ELS) is commercial structural-analysis software based on the applied element method (AEM) for the automatic tracking and propagation of cracks, separation of elements, element collision, and collapse of structures under extreme loads. AEM combines features of Finite element method and Discrete element method simulation with its own solver capabilities for the generation of PC-based structural analysis.

The concept of angles between lines in the plane and between pairs of two lines, two planes or a line and a plane in space can be generalized to arbitrary dimension. This generalization was first discussed by Jordan. For any pair of flats in a Euclidean space of arbitrary dimension one can define a set of mutual angles which are invariant under isometric transformation of the Euclidean space. If the flats do not intersect, their shortest distance is one more invariant. These angles are called canonical or principal. The concept of angles can be generalized to pairs of flats in a finite-dimensional inner product space over the complex numbers.

The fracture of soft materials involves large deformations and crack blunting before propagation of the crack can occur. Consequently, the stress field close to the crack tip is significantly different from the traditional formulation encountered in the Linear elastic fracture mechanics. Therefore, fracture analysis for these applications requires a special attention. The Linear Elastic Fracture Mechanics (LEFM) and K-field are based on the assumption of infinitesimal deformation, and as a result are not suitable to describe the fracture of soft materials. However, LEFM general approach can be applied to understand the basics of fracture on soft materials. The solution for the deformation and crack stress field in soft materials considers large deformation and is derived from the finite strain elastostatics framework and hyperelastic material models.

References

1 2 Meguro, K.; Tagel-Din, H. (2000). "Applied element method for structural analysis: Theory and application for linear materials". Structural Engineering/Earthquake Engineering. Japan: Japan Society of Civil Engineers. 17 (1): 21–35. F0028A. Archived from the original on 2012-02-29. Retrieved 2009-08-10.
↑ Tagel-Din, H.; Meguro, K (2000). "Applied Element Method for Simulation of Nonlinear Materials: Theory and Application for RC Structures". Structural Engineering/Earthquake Engineering. Japan: Japan Society of Civil Engineers. 17 (2): 137–148. Retrieved 2009-08-10.
↑ Tagel-Din, H.; Meguro, Kimiro (November 2001). "Applied Element Simulation of RC Structures under Cyclic Loading". Journal of Structural Engineering. Japan: ASCE. 127 (11): 137–148. doi:10.1061/(ASCE)0733-9445(2001)127:11(1295). ISSN 0733-9445 . Retrieved 2009-08-10.
↑ Tagel-Din, H.; Meguro, K (2002). "AEM Used for Large Displacement Structure Analysis" (PDF). Journal of Natural Disaster Science. Japan. 24 (1): 25–34. Retrieved 2009-08-10.^{[ permanent dead link ]}
↑ Tagel-Din, Hatem; Kimiro Meguro, K (January 30 – February 4, 2000). Analysis of a Small Scale RC Building Subjected to Shaking Table Tests using Applied Element Method. New Zealand: Proceedings of the 12th World Conference on Earthquake Engineering. pp. 25–34.
↑ HATEM, Tagel-Din; Kimiro MEGURO, K (August 1–6, 2004). Dynamic Modeling of Dip-Slip Faults for Studying Ground Surface Deformation Using Applied Element Method. Vancouver, Canada: Proceedings of the 13th World Conference on Earthquake Engineering.
↑ Mayorka, Paola; Kimiro Meguro, K (October 2003). "Modeling Masonry Structures using the Applied Element Method". Seisan Kenkyu. Japan: Institute of Industrial Science, The University of Tokyo. 55 (6): 123–126. ISSN 1881-2058 . Retrieved 2009-08-10.
↑ Mayorka, Paola; Kimiro Meguro, K (2005). Blast Testing and Research Bridge at the Tenza Viaduct. Japan: University of Missouri-Rolla, TSWG Contract Number N4175-05-R-4828, Final Report of Task 1.