Derivation of 3D masonry properties using numerical homogenization technique

Lots of research work has been conducted on homogenization technique, which derives global homogenized properties of masonry from the behaviour of the constitutive materials (brick and mortar). Such a technique mainly focused on two‐dimensional media in the previous studies with the out‐of‐plane properties of masonry material neglected. In this paper, homogenization technique and damage mechanics theory are used to model a three‐dimensional masonry basic cell to numerically derive the equivalent elastic properties, strength envelope, and failure characteristics of masonry material. The basic cell is modelled with distinctive consideration of non‐linear material properties of mortar and brick. Various displacement boundaries are applied on the basic cell surfaces in the numerical simulation. The detailed material properties of mortar and brick are modelled in a finite element program in the numerical analysis. The stress–strain relations of masonry material under various conditions are obtained from the simulation. The homogenized elastic properties and failure characteristics of masonry material are derived from the simulation results. The homogenized 3D model is then utilized to analyse the response of a masonry panel to airblast loads. The same panel is also analysed with distinctive material modelling. The efficiency and accuracy of the homogenized model are demonstrated. The homogenized material properties and failure model can be used to model large‐scale masonry structure response. Copyright © 2006 John Wiley & Sons, Ltd.     

Experimentelle und analytische Untersuchung des out‐of‐plane‐Verhaltens von unbewehrten Mauerwerkswänden unter Erdbebenbelastung

Experimental and analytical investigation of the seismic out‐of‐plane behavior of unreinforced masonry walls In addition to the vertical and horizontal load‐bearing in‐plane, masonry must also withstand out‐of‐plane loads that occur in earthquake scenarios. The out‐of‐plane behavior of unreinforced masonry walls depends on a variety of parameters and is very complex due to the strong non‐linearity. Current design methods in German codes and various international codes have not been explicitly developed for out‐of‐plane behavior and contain considerable conservatism. In the present work, shaking‐table experiments with heat‐insulating masonry walls have been conducted to investigate the out‐of‐plane behavior of vertical spanning unreinforced masonry walls. As shown in previous numerical investigations, important parameters are neglected in existing design and analysis models and the out‐of‐plane capacity is underestimated significantly. In the conducted experiments the results of these numerical investigations are verified. Furthermore, the development of an analytical design model to determine the force‐displacement relationship and the out‐of‐plane load‐bearing capacity considering all significant parameters is presented.     

Masonry infilled reinforced concrete frames under horizontal loading / Stahlbetonrahmen mit Ausfachungen aus Mauerwerk unter horizontalen Belastungen

The behaviour of infilled reinforced concrete frames under horizontal load has been widely investigated, both experimentally and numerically. Since experimental tests represent large investments, numerical simulations offer an efficient approach for a more comprehensive analysis. When RC frames with masonry infill walls are subjected to horizontal loading, their behaviour is highly non‐linear after a certain limit, which makes their analysis quite difficult. The non‐linear behaviour results from the complex inelastic material properties of the concrete, infill wall and conditions at the wall‐frame interface. In order to investigate this non‐linear behaviour in detail, a finite element model using a micro modelling approach is developed, which is able to predict the complex non‐linear behaviour resulting from the different materials and their interaction. Concrete and bricks are represented by a non‐linear material model, while each reinforcement bar is represented as an individual part installed in the concrete part and behaving elasto‐plastically. Each brick is modelled individually and connected taking into account the non‐linearity of a brick mortar interface. The same approach is followed using two finite element software packages and the results are compared with the experimental results. The numerical models show a good agreement with the experiments in predicting the overall behaviour, but also very good matching for strength capacity and drift. The results emphasize the quality and the valuable contribution of the numerical models for use in parametric studies, which are needed for the derivation of design recommendations for infilled frame structures.     

A refined non‐linear non‐conforming triangular plate/shell element

A refined non‐conforming triangular plate/shell element for geometric non‐linear analysis of plates/shells using the total Lagrangian/updated Lagrangian approach is constructed in this paper based on the refined non‐conforming element method for geometric non‐linear analysis. The Allman's triangular plane element with vertex degrees of freedom and the refined triangular plate‐bending element RT9 are used to construct the present element. Numerical examples demonstrate that the accuracy of the new element is quite high in the geometric non‐linear analysis of plates/shells. Copyright © 2003 John Wiley & Sons, Ltd.     

Multi‐time‐step explicit–implicit method for non‐linear structural dynamics

We present a method with domain decomposition to solve time‐dependent non‐linear problems. This method enables arbitrary numeric schemes of the Newmark family to be coupled with different time steps in each subdomain: this coupling is achieved by prescribing continuity of velocities at the interface. We are more specifically interested in the coupling of implicit/explicit numeric schemes taking into account material and geometric non‐linearities. The interfaces are modelled using a dual Schur formulation where the Lagrange multipliers represent the interfacial forces. Unlike the continuous formulation, the discretized formulation of the dynamic problem is unable to verify simultaneously the continuity of displacements, velocities and accelerations at the interfaces. We show that, within the framework of the Newmark family of numeric schemes, continuity of velocities at the interfaces enables the definition of an algorithm which is stable for all cases envisaged. To prove this stability, we use an energy method, i.e. a global method over the whole time interval, in order to verify the algorithms properties. Then, we propose to extend this to non‐linear situations in the following cases: implicit linear/explicit non‐linear, explicit non‐linear/explicit non‐linear and implicit non‐linear/explicit non‐linear. Finally, we present some examples showing the feasibility of the method. Copyright © 2001 John Wiley & Sons, Ltd.     

Coupled mechanics and finite element for non‐linear laminated piezoelectric shallow shells undergoing large displacements and rotations

A theoretical framework is presented for analysing the coupled non‐linear response of shallow doubly curved adaptive laminated piezoelectric shells undergoing large displacements and rotations. The formulated mechanics incorporate coupling between in‐plane and flexural stiffness terms due to geometric curvature, coupling between mechanical and electric fields, and encompass geometric non‐linearity effects due to large displacements and rotations. The governing equations are formulated explicitly in orthogonal curvilinear co‐ordinates and are combined with the kinematic assumptions of a mixed‐field shear‐layerwise shell laminate theory. Based on the above formulation, a finite element methodology together with an incremental‐iterative technique, based on Newton–Raphson method is formulated. An eight‐node coupled non‐linear shell element is also developed. Various evaluation cases on laminated curved beams and cylindrical panels illustrate the capability of the shell finite element to predict the complex non‐linear behaviour of active shell structures including buckling, which is not captured by linear shell models. The numerical results also show the inherent capability of piezoelectric shell structures to actively induce large displacements through piezoelectric actuators, by jumping between multiple equilibrium states. Copyright © 2005 John Wiley & Sons, Ltd.     

Analysis of steel‐reinforced masonry walls regarding maximum shear loads / Traglastberechnung von vertikal bewehrten Mauerwerksscheiben

Loadings on masonry for the earthquake case pose particular challenges for the material. In order to improve the load‐bearing and deformation behaviour, masonry building elements can be strengthened with reinforcement. This article presents an analytical model for the calculation of the load‐bearing capacity of vertically reinforced masonry panels. The masonry is modelled as a homogeneous and anisotropic material and failure conditions are based on the plastic theory. Using uniaxially loaded stress fields and considering the structural constraints, a lower load‐bearing threshold can be given. In order to verify the model, three shear tests on reinforced sand‐lime block masonry were recalculated regarding their load‐bearing capacity. The test panels each contained vertical steel reinforcement in the blocks. The blocks were laid in thin bed mortar.     

Efficient non‐linear reanalysis of skeletal structures using combined approximations

Non‐linear reanalysis of large‐scale structures usually involves much computational effort, because the set of non‐linear equations must be solved repeatedly during the solution process. Various approximations that are often used for linear reanalysis are not sufficiently accurate for non‐linear problems. In this study, solution procedures based on the combined approximations approach are developed and compared in terms of efficiency and accuracy. Various path‐independent non‐linear analysis and reanalysis problems are considered, including material non‐linearity, geometric non‐linearity and buckling analysis. Numerical examples demonstrate the effectiveness of the procedures presented. It is shown that in various cases accurate results can be achieved efficiently. Copyright © 2007 John Wiley & Sons, Ltd.     

Model validation and parametric study on the blast response of unreinforced brick masonry walls

Numerical simulations are carried out to estimate the response and damage of unreinforced brick masonry walls subjected to explosive blast loading based on the transient dynamic finite element program LS-DYNA. A previously developed dynamic plastic damage model was used for brick and mortar. A new model for strain rate effects of bricks and mortar is included in the numerical analysis. The results obtained from the numerical models are compared with field test data and good agreement can be found. Parametric studies are conducted to evaluate the effect of material strength, boundary conditions, and thickness of the wall on the blast response of unreinforced brick masonry walls. It was found that boundary conditions and wall thickness significantly affect the blast response, while the effect of material strength is relatively small.     

Classification of SIM infill panels / Klassifikation von SIM‐Ausfachungswänden

SIM is an innovative building system for mortar‐less walls. It utilises a special method of interlocking SIM bricks that allows relative sliding of brick courses in‐plane of a wall and prevents out‐of‐plane relative movement of bricks. One of its structural applications is in multistorey frame buildings as earthquake resistant masonry infill panels. It improves energy dissipation of frame structures during earthquakes. The energy dissipation occurs through friction between bricks as they engage in relative sliding by the frame vibrating during earthquake. This paper explains the novelty of SIM and offers classification of SIM panels based on the gap width between the frame and the top of the panel.     

Investigation of the seismic out‐of‐plane behaviour of unreinforced masonry walls

The structural stability of unreinforced masonry (URM) walls has to be guaranteed not only under static (permanent and live) loads but also under earthquake loads. Loads transverse to the plane (out‐of‐plane) often have a decisive influence on the load‐bearing capacity. In practical applications, simplified methods from codes, guidelines and literature are often used to analyse and evaluate the out‐of‐plane capacity of load‐bearing and non‐load‐bearing URM walls. The results of these simplified methods can be significantly conservative and inaccurate since essential influencing effects are neglected. For many existing buildings, the simplified methods underestimate the capacity, which leads to cost‐intensive retrofitting and strengthening measures or complete replacement by other wall systems. In order to realistically estimate the out‐of‐plane capacity, parameters such as wall geometry, boundary conditions, vertical loads and especially dynamic effects (e.g. inertia forces) have to be taken into account. In this paper, non‐linear time history simulations are presented to investigate the influence of these effects. The numerically determined maximum acceptable earthquake acceleration is compared with results from simplified analysis models. The comparison shows that the out‐of‐plane capacity is significantly higher than the values predicted by simplified models. Finally, several initial experimental seismic tests conducted on the shaking table of the TU Kaiserslautern are presented, together with the planned extensive experimental test program on the out‐of‐plane capacity of masonry walls.     

c‐Type method of unified CAMG and FEA. Part 1: Beam and arch mega‐elements—3D linear and 2D non‐linear

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.     

On the design of energy–momentum integration schemes for arbitrary continuum formulations. Applications to classical and chaotic motion of shells

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.     

A non‐linear programming method approach for upper bound limit analysis

This paper presents a finite element model based on mathematical non‐linear programming in order to determine upper bounds of colapse loads of a mechanical structure. The proposed formulation is derived within a kinematical approach framework, employing two simultaneous and independent field approximations for the velocity and strain rate fields. The augmented Lagrangian is used to establish the compatibility between these two fields. In this model, only continuous velocity fields are used. Uzawa's minimization algorithm is applied to determine the optimal kinematical field that minimizes the difference between external and dissipated work rate. The use of this technique allows to bypass the complexity of the non‐linear aspects of the problem, since non‐linearity is addressed as a set of small local subproblems of optimization for each finite element. The obtained model is quite versatile and suitable for solving a wide range of collapse problems. This paper studies 3D strut‐and‐tie structures, 2D plane strain/stress and 3D solid problems. Copyright © 2007 John Wiley & Sons, Ltd.     

Buckling of thin‐walled cylindrical shells under axial compression

Lightweight thin‐walled cylindrical shells subjected to external loads are prone to buckling rather than strength failure. The buckling of an axially compressed shell is studied using analytical, numerical and semi‐empirical models. An analytical model is developed using the classical shell small deflection theory. A semi‐empirical model is obtained by employing experimental correction factors based on the available test data in the theoretical model. Numerical model is built using ANSYS finite element analysis code for the same shell. The comparison reveals that the analytical and numerical linear model results match closely with each other but are higher than the empirical values. To investigate this discrepancy, non‐linear buckling analyses with large deflection effect and geometric imperfections are carried out. These analyses show that the effects of non‐linearity and geometric imperfections are responsible for the mismatch between theoretical and experimental results. The effect of shell thickness, radius and length variation on buckling load and buckling mode has also been studied. Copyright © 2009 John Wiley & Sons, Ltd.     

A two‐scale approximation of the Schur complement and its use for non‐intrusive coupling

This paper presents a two‐scale approximation of the Schur complement of a subdomain's stiffness matrix, obtained by combining local (i.e. element strips) and global (i.e. homogenized) contributions. This approximation is used in the context of a coupling strategy that is designed to embed local plasticity and geometric details into a small region of a large linear elastic structure; the strategy consists in creating a local model that contains the desired features of the concerned region and then substituting it into the global problem by the means of a non‐intrusive solver coupling technique adapted from domain decomposition methods. Using the two‐scale approximation of the Schur complement as a Robin condition on the local model enables to reach high efficiency. Examples include a large 3D problem provided by our industrial partner Snecma. Copyright © 2011 John Wiley & Sons, Ltd.     

Shear tests on full scale storey‐height specimens constructed with thermally insulating clay unit masonry with thin‐layer mortar

The paper presents results of a series of 6 in‐plane shear tests on storey‐height clay unit masonry panels [1] with thin‐layer mortar, carried out in addition to previous test campaigns [2], [3], and [4]. The walls were constructed with unfilled thermally insulating clay units with a thermal conductivity of λ = 0.09 W/(m · K). The current design rules for clay unit masonry according to DIN EN 1996‐1‐1/NA [5] are conservative compared to the presented test results for thermally insulating clay unit masonry.     

Mixed‐mode delamination in 2D layered beam finite elements

In this work, a 2D finite element (FE) formulation for a multi‐layer beam with arbitrary number of layers with interconnection that allows for mixed‐mode delamination is presented. The layers are modelled as linear beams, while interface elements with embedded cohesive‐zone model are used for the interconnection. Because the interface elements are sandwiched between beam FEs and attached to their nodes, the only basic unknown functions of the system are two components of the displacement vector and a cross‐sectional rotation per layer. Damage in the interface is modelled via a bi‐linear constitutive law for a single delamination mode and a mixed‐mode damage evolution law. Because in a numerical integration procedure, the damage occurs only in discrete integration points (i.e. not continuously), the solution procedure experiences sharp snap backs in the force‐displacements diagram. A modified arc‐length method is used to solve this problem. The present model is verified against commonly used models, which use 2D plane‐strain FEs for the bulk material. Various numerical examples show that the multi‐layer beam model presented gives accurate results using significantly less degrees of freedom in comparison with standard models from the literature. Copyright © 2015 John Wiley & Sons, Ltd.     

A novel non‐linearly explicit second‐order accurate L‐stable methodology for finite deformation: hypoelastic/hypoelasto‐plastic structural dynamics problems

A novel non‐linearly explicit second‐order accurate L‐stable computational methodology for integrating the non‐linear equations of motion without non‐linear iterations during each time step, and the underlying implementation procedure is described. Emphasis is placed on illustrative non‐linear structural dynamics problems employing both total/updated Lagrangian formulations to handle finite deformation hypoelasticity/hypoelasto‐plasticity models in conjunction with a new explicit exact integration procedure for a particular rate form constitutive equation. Illustrative numerical examples are shown to demonstrate the robustness of the overall developments for non‐linear structural dynamics applications. Copyright © 2004 John Wiley & Sons, Ltd.     

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

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.