# 6. NON-LINEAR PSEUDO-STATIC ANALYSIS OF ADOBE WALLS

Save this PDF as:

Size: px
Start display at page:

## Transcription

2 114 Sabino Nicola Tarque Ruíz Amongst the different solution procedures used in the implicit finite element solvers, the Newton-Raphson solution procedure is the faster for solving non-linear problems under force control, though no convergence procedure is full proof. When solving quasi static problems, a set of non-linear equations are generally expressed as: T T dv ds G u B u N t (6.2) V S where G is a set of non-linear equations in u which are updated at each iteration and are function of the nodal displacements vector u, is the stresses vector, B is the matrix that relates the strain vector to the displacements, N is the matrix of element shape functions and t is the surface traction vector. V and S represent the volume and surface of the body, respectively. The right-hand side of Equation (6.2) represents the difference between the internal and the external forces. Equation (6.2) is solved for a displacement vector that equilibrates the internal and external forces, as explained in Harewood and McHugh [27]. Besides, Equation (6.2) is solved by incremental methods where load/displacements are applied in time steps t, t t t u is solved from a known state u. In the following equations the subscript denotes iteration number and the superscript denotes an increment step. At the beginning, Equation (6.2) is solved to obtain the displacement correction u t t based on information of u as: u 1 t t G u t t i t t i 1 G ui u i The partial derivative on the right side is the so called Jacobian matrix expressed as the global stiffness matrix K : G u t t Ktan u i = u Thus: t t i tan 1 t t i 1 (6.3) (6.4) t t t t t t ui 1 K tan ui G u i (6.5) The previous equation involves the inversion of the global stiffness matrix (no singular matrix), which ends in a computationally expensive operation, but it ensures that a

3 Numerical modelling of the seismic behaviour of adobe buildings 115 relatively large time increment can be used while maintaining the accuracy of the solution [Harewood and McHugh 27]. t t The displacement correction u i 1 is added to the previous state, thus the improved displacement solution is given with Equation (6.6) and t u t i 1 is used as the current approximation to the solution for the subsequent iteration i 1 in Equation (6.2) until the equilibrium is reached. t t t t t t i 1 i i 1 u u u (6.6) Convergence is measured by ensuring that the difference between external and internal forces G u, displacement increment t t t t ui and displacement correction u i 1 are sufficiently small. The difference between several solution procedures is the way in which t t u i 1 is determined. As an example, Figure 6.1 shows the iteration process followed with a Newton-Rapshon procedure for reaching equilibrium in an implicit procedure. Figure 6.1. Iteration process for an implicit solution. 6.2 DISCRETE MODEL: MODELLING THE PUSHOVER RESPONSE OF AN ADOBE WALL The adobe wall I-1 presented in section 3.5 is now created in Midas FEA using linear elastic solid elements and zero thickness non-linear interface elements. The combined cracking-shearing-crushing model is used for representing the non-linear behaviour of the adobe masonry lumped at the mortar joints. A total of 17 courses were placed to model the adobe masonry. The original dimensions of the unit bricks (.3x.13x.1 and.22x.13x.1 m) were extended to take into account the 13~15 mm mortar thickness. The model includes the reinforced concrete beams (at the top: crown beam, and at the base: foundation), the adobe walls and the timber lintel. The base of the foundation is

4 116 Sabino Nicola Tarque Ruíz fully fixed; the top part of the crown beam is free. Since the test was displacement controlled, in the numerical model a monotonic top-displacement was applied at one vertical edge of the top concrete beam to reach a maximum displacement of 1 mm. The material parameters are changed in the finite element model for a parametric study. a) Complete model view b) View of the interface elements c) In red the applied loads Figure 6.2. Finite element model of the adobe wall subjected to horizontal displacement loads at the top. Discrete model, Midas FEA. The configuration of the numerical model in Midas FEA is shown in Figure All the solid elements are considered elastic and isotropic and the properties, shown in Table 6.1,

6 118 Sabino Nicola Tarque Ruíz are related to the elasticity modulus of the adobe bricks and the mortar joints, as seen in Equation (5.12) and repeated here for convenience: k n Eunit Emortar h E E mortar unit mortar (6.7) Gunit Gmortar ks h G G mortar unit mortar The elasticity modulus for the adobe bricks is assumed as 23 MPa. The elasticity modulus for the mud mortar is assumed to be lower than that for bricks, so values of 79.25, 113.4, and MPa are considered. The mortar thickness, h mortar, is around 15 mm. The shear modulus is taken as.4e for the bricks and mortar. With the previous values, the normal and shear penalty stiffness are those reported in Table 6.2. (6.8) Table 6.2. Variation of elasticity of module of the mortar joints for evaluating the penalty stiffness kn and k s. Mortar thickness (mm) E unit (MPa) E mortar (MPa) k n (N/mm 3 ) K s (N/mm 3 ) In order to make the first round of analysis the other material parameters are defined as shown in Table 6.3, where c is the cohesion, o is the friction angle, is the dilatancy I angle, r is the residual friction angle, f t is the tensile strength, G f is the fracture energy for Mode I (related to tension softening), a and b are factors to evaluate the fracture energy for Mode II (expressed as G II = a. +b and related to the shear behaviour), f c is the compression strength C s is the shear tension contribution factor (which is 9 according to Lourenço [1996], G c f is the compressive fracture energy, and k p is the peak equivalent plastic relative displacement. The parameters which are marked with * are considered for the parametric study and those shown in Table 6.3 are the ones which gave good numerical experimental agreementin terms of global response of the masonry wall (as discussed later). The inelastic parameter values for the adobe masonry were taken lower than those proposed by Lourenço [1996] for clay masonry based on the clear difference in the material strengths.

7 Numerical modelling of the seismic behaviour of adobe buildings 119 k n (N/mm 3 ) Table 6.3. Preliminary material properties for the interface model (mortar joints). k t (N/mm 3 ) Structural c (N/mm 2 ) o (deg) (deg) r (deg) f t (N/mm 2 )* Mode I I G (N/mm)* f Table 6.2 Table Mode II Compression cap a (mm)* b (N/mm)* f c (N/mm 2 ) * C s c G f (N/mm) * k p (mm) * Figure 6.3 shows the effect of the variation of the penalty stiffness k n and k s, on the elastic part of the pushover curve. As expected, it is seen that greater stiffness values lead to a stiffer elastic response. It is shown that k n 8.9 N / mm 3 and ks 3.23 N / mm 3, which implies E unit 23MPa and Emortar 8MPa, match well the elastic part and the yielding initiation, so these values were fixed to analyze other material parameters. 1 8 Force (kn) 6 4 kn= 8.9 N/mm³ 2 kn= N/mm³ kn= N/mm³ Kn= 253 N/mm³ Figure 6.3. Comparison of the pushover curves in models with different penalty stiffness values. Discrete model, Midas FEA. Variation of the tensile strength f t of the mortar was applied to study its influence. The fracture energy G I f was maintained constant and equal to.8 N/mm, which means that the area under each curve is the same. So, due to this the crack displacement values related to low tensile strength are greater than those related to high tensile strength (e.g. see curve for f t.25 MPa in Figure 6.4). The crack displacement is equal to the tensile strain times the thickness of the mortar joint, which is around 15 mm.

8 12 Sabino Nicola Tarque Ruíz.25 Tensile strength (MPa) ft=.2 MPa ft=.1 MPa ft=.5 MPa ft=.25 MPa Crack displacement (mm) Figure 6.4. Variation of tensile strength I f t at the interface, G f =.8 N/mm, h = 15 mm. Figure 6.5 shows the pushover curves obtained varying the tensile strength. It is seen that the global strength of the masonry depends on the crack displacement values in the tensile softening. Greater crack displacement values produce an increment on the lateral resistant at the post-yield peak. When the tensile softening part of the constitutive law descends abruptly (e.g. f t.2 MPa curve in Figure 6.4), the pushover curve stops after the yielding initiation because can not develop greater crack displacements. All the curves, except the one related to f t.2 MPa, have similar global behaviour until 6.2 mm top horizontal displacement. The best match in terms of crack pattern and strength was obtained with f t.1mpa Force (kn) ft=.25 MPa ft=.5 MPa ft=.1 MPa ft=.2 MPa Figure 6.5. Influence of the tensile strength t f on the pushover response. Discrete model, Midas FEA.

9 Numerical modelling of the seismic behaviour of adobe buildings Tensile strength (MPa) Gf=.15 N/mm Gf=.8 N/mm Gf=.5 N/mm Crack displacement (mm) Figure 6.6. Variation of fracture energy I G f at the interface, constant f t =.1 MPa, h = 15 mm. I A variation of the tensile fracture energy G f with constant tensile strength f t was applied to study its influence on the structural response (Figure 6.6). An increment in the lateral strength of the masonry is observed when the tensile softening part does not descent abruptly (Figure 6.7), which partially concludes that greater fracture energy gives greater lateral strength Force (kn) Gf=.5 N/mm Gf=.8 N/mm Gf=.15 N/mm Figure 6.7. Influence of the fracture energy in tension model, Midas FEA. G I f on the pushover response. Discrete A parametric study was carried out for evaluating the compression strength of adobe masonry. Since there is not information about the inelastic material properties, a hardening/softening curve was assumed taking into consideration the experimental results obtained by Lourenço [1996] for clay masonry. It was decided to keep the peak compressive plastic strain around.8 mm/mm (related to k p =.9 mm) and the

10 122 Sabino Nicola Tarque Ruíz compression fracture energy of.2 N/mm, which can be considered a lower bound value of a real compression curve (Figure 6.8)..5 Compression strength (MPa) fc=.4 MPa fc=.25 MPa fc=.15 MPa Strain (mm/mm) C Figure 6.8. Variation of compression strength f c at the interface, constant G f =.2 N/mm, h = 115 mm. The characteristic length for compression, h, is taken as twice the half height of the adobe brick plus the mortar thickness, that is 115 mm. The relative peak compressive displacement, k p, is taken as the peak compressive plastic strain times the characteristic length. Figure 6.9 shows that the failure process does not depend at all on the compression behaviour Force (kn) fc=.15 MPa fc=.25 MPa fc=.4 MPa Figure 6.9. Influence of compression strength FEA. f c on the pushover response. Discrete model, Midas A variation of the compressive fracture energy G was applied to analyze its effect on the global response of the masonry wall. The plastic strain related to f c.25 MPa was kept f c

11 Numerical modelling of the seismic behaviour of adobe buildings 123 constant and close to.8 mm/mm (Figure 6.1). The results showed no variation on the global response of the adobe wall..3 Compression strength (MPa) Gf=.4 N/mm Gf=.2 N/mm Gf=.1 N/mm Figure 6.1. Variation of fracture energy Strain (mm/mm) C G f at the interface, constant f c =.25 MPa, h = 115 mm. Another studied considered a variation of the relative peak compressive displacement k p while maintaining a constant compressive strength (f c.25 MPa ) and a constant f compressive fracture energy ( Gc.2 N / mm), Figure As in the previous cases, no variation in the global response was observed. The pushover curves were more or less the same as the one shown in Figure Compression strength (MPa).2.1 kp=.9 mm kp=.15 mm kp=.2 mm Figure Variation of N/mm, h = 115 mm. Strain (mm/mm) k p for the compression curve, constant f c =.25 MPa and G C f =.2 It should be said that Midas FEA does not include the compression cap model when the analysis refers to a three-dimensional interface model. For this reason no variation on the pushover curves was observed varying the compression strength. Finally, a variation of

12 124 Sabino Nicola Tarque Ruíz the Mode II fracture energy (shear) was applied to study the shear behaviour of the mortar joints. The dilatation angle was assumed zero. In this case it seems that fracture energy of.1 N/mm can be considered for adobe masonry (Figure 6.12) Force (kn) Gf=.2 N/mm Gf=.5 N/mm Gf=.1 N/mm II Figure Influence of the shear fracture energy G f at the interface on the pushover response. Discrete model, Midas FEA Results of the pushover analysis considering a discrete model Figure 6.13 shows the sequence of damage obtained with the selected parameters specified in Table 6.1 and Table 6.3. The crack pattern follows the experimental results: the cracks go from the top left (where the load is applied) to the right bottom of the wall (Figure 6.14). Also, the horizontal cracks at the transversal walls are in agreement with the experimental results. Since the load applied is monotonic, the FE model cannot capture the X-shape failure observed in the experiments. The presence of the two transversal walls prevents the rocking behaviour, representing correctly the tested wall. The maximum displacement reached at the top of the wall is around 6.2 mm, after which the program stopped due to convergence problems. The comparison between the numerical and experimental pushover response is shown in Figure 6.15.

13 Numerical modelling of the seismic behaviour of adobe buildings 125 a) Top displacement= 1 mm b) Top displacement= 2 mm c) Top displacement= 4 mm d) Top displacement= 6.26 mm Figure Damage pattern of the adobe wall subjected to a horizontal top displacement. e) Top displacement= 6.26 mm, isometric view Figure Continuation. Damage pattern of the adobe wall subjected to a horizontal top displacement. Discrete model, Midas FEA.

14 126 Sabino Nicola Tarque Ruíz Figure damage pattern for wall I-1 due to cyclic displacements applied at the top. Just the adobe wall is shown here, the concrete beam and foundation are hidden, [Blondet et al. 25]. All the models were run in Midas FEA with arc-length method with initial stiffness. The number of load steps was specified as 1, the initial load factor was.1, and the maximum number of iteration per load step was 3. The convergence criteria were given by an energy norm and displacement norm of Force (kn) Numerical Figure Load-displacement diagrams, experimental and numerical. 6.3 TOTAL-STRAIN MODEL: MODELLING THE PUSHOVER RESPONSE In this part the adobe wall tested by Blondet et al. [25] is modelled using a continuum approach. A plane stress finite element model is created in Midas FEA using 4-node rectangular shell elements (Figure 6.16a) and considering drilling DOFs and transverse shear deformation. The size of the mesh is usually kept at 1 x 1 mm, which is related to a characteristic length dimension h= mm, obtained from the square root of the area of the shell element [Bažant and Oh 1983]. The thickness of the shell is 3 mm. The adobe masonry includes the adobe bricks and the mud mortar joints; in this case, a

15 Numerical modelling of the seismic behaviour of adobe buildings 127 homogeneous material is assumed and the cracks are smeared into the continuum. The top and bottom reinforced concrete beams and the timber lintel are considered elastic. The foundation is fully fixed at the base. The crown beam is at the top. The numerical model is subjected to a unidirectional displacement imposed at the two ends of the top crown beam (Figure 6.16b). The elastic material properties for the concrete beam and the timber lintel are given in Table 6.1, while those for the adobe masonry are specified in Table 6.4. The material properties marked with * are calibrated based on the experimental pushover curve as discussed later. Equation (5.29) and Equation (5.37) are considered for computing the inelastic part of the tension and compression constitutive law, respectively. i is taken greater than f c /3 to maintain a parabolic shape of the compression curve. E (N/mm 2 )* Table 6.4. Material properties for the adobe masonry within total-strain model. Elastic Tension Compression m (N/mm 3 ) h (mm) f t (N/mm 2 )* I G f (N/mm)* f c (N/mm 2 )* c G f (N/mm)* p (mm/mm)* 2.2 2e a) Complete view of the model b) Position of horizontal applied loads Figure Finite element model of the adobe wall subjected to horizontal displacement loads at the top. Total-strain model, Midas FEA.

16 128 Sabino Nicola Tarque Ruíz Calibration of material properties The first parameter that was calibrated is the elasticity modulus E of the adobe masonry. According to section and 3.7, the E value can be considered between 2 and 22 MPa. Besides, [Blondet and Vargas 1978] suggests to use E= 17 MPa; however, this value seems to be too conservative. In this work E= 2 MPa has been considered for all the numerical analyses since it yields a good agreement between the numerical and experimental curves. Figure 6.17 shows the variation of the in-plane response of the adobe wall due to the variation of elasticity modulus. The other, elastic and inelastic, material properties used were the ones specified in Table Force (kn) E= 25 MPa E= 22 MPa E= 2 MPa E= 15 MPa Figure Comparison of the pushover curves in models with different E. Total-strain model. The next parameter that was calibrated was the tensile strength f t of the masonry, which can be roughly assumed around 1% of the compression strength f c. The tensile fracture energy G is maintained in all cases as.1 N/mm (Figure 6.18). I f.7.6 Tensile strength (MPa) ft=.2 MPa ft=.4 MPa ft=.6 MPa Crack displacement (mm) I Figure Variation of tensile strength for total-strain model, constant G f =.1 N/mm, h = mm.

17 Numerical modelling of the seismic behaviour of adobe buildings 129 Since the area under the tensile softening curve is fixed for all values of the tensile strength, the crack displacement values are greater for lower tensile strengths (Figure 6.18). Figure 6.19 shows that the most accurate pushover curve is obtained when f t.4 MPa. Lower values of f t reduce the in-plane strength of the masonry. However, they also give more stable pushover curves due to the large values of the crack displacement (see Figure 6.18). On the other hand, larger values of f t increase the seismic in-plane capacity, but the pushover curve stops due to convergence problems Force (kn) ft=.2 MPa ft=.4 MPa ft=.6 MPa Figure Influence of the tensile strength FEA. ft on the pushover response. Total-strain model, Midas Tensile strength (MPa) Gf=.5 N/mm Gf=.1 N/mm Gf=.15 N/mm Crack displacement (mm) Figure 6.2. Variation of fracture energy G I f for total-strain model, constant f t =.4 MPa, h = mm.

18 13 Sabino Nicola Tarque Ruíz I A variation of the tensile fracture energy G f was considered.5,.1 and.15 N/mm (Figure 6.2). The tensile strain values are equal to the crack displacements divided by the characteristic element length h (141.4 mm). It is seen that greater the fracture energy, the larger the crack displacement values. The best fit of the experimental pushover curve is obtained with G.1MPa (Figure 6.21), in terms of both wall strength and crack pattern. This preliminary study concludes that even though adobe is a very brittle material, it still retains some tension fracture energy, which controls the crack formation process. I f Force (kn) Gf=.5 N/mm Gf=.1 N/mm Gf=.15 N/mm I Figure Influence of the fracture energy in tension G f on the pushover response. Total-strain model, Midas FEA. The smeared crack model takes into account the effect of shear through a reduction factor that multiplies the shear stiffness. This is only possible when a fixed crack model is used, as is the case for this work (see section 4.2.2). A variation of is analyzed to see how much this can influence the global numerical response. As shown in Figure 6.22, the best numerical result was obtained considering.5. The compression strength f c of the adobe masonry was varied from.3 to.8 MPa. The hardening/softening curve is similar to the one used by Lourenço [1996] for clay masonry C but proportionally scaled for adobe masonry. For this reason the ratio Gf / fc is kept at about.344 mm. The peak plastic compression strain p is kept at.2 mm/mm, and the plastic strain at 5% of the compression strength m is kept at.5 mm/mm in all cases. The different compression curves are shown in Figure According to the experimental data, it is seen that a reasonable f c value can be greater than.5 MPa; however, this value is calibrated for within a smeared crack approach.

19 Numerical modelling of the seismic behaviour of adobe buildings Force (kn) bheta=.2 N/mm bheta=.5 N/mm bheta=.1 N/mm Figure Influence of the shear retention factor on the pushover response. Total-strain model, Midas FEA. The total peak compression strain values shown in Figure 6.23, which are the sum of the elastic plus the peak plastic strain, are.345,.41,.5275 and.56 mm/mm for f =.3,.45,.7 and.8 MPa, respectively. c.9 Compression strength (MPa) fc=.3 MPa fc=.45 MPa fc=.7 MPa fc=.8 MPa Strain (mm/mm) Figure Variation of compression strength f c for total-strain model. Relation mm in all cases, h = mm. G C f / f =.344 c Figure 6.24 shows the numerical pushover curves obtained with different compression strength values. It can be observed that the best results are obtained with f c =.3 MPa. The pushover curves obtained with other compression strengths are superimposed and give higher strength, but stop earlier due to convergence problems due to the

20 132 Sabino Nicola Tarque Ruíz concentration of compression stress at the right top window corner. If convergence is reached, so the pushover curve will down and continues closes to the experimental curve Force (kn) fc=.45 MPa fc=.3 MPa fc=.7 MPa fc=.8 MPa Figure Influence of the compression strength Midas FEA. fc on the pushover response. Total-strain model, Results of the pushover analysis considering a total-strain model Figure 6.25 shows the sequence of damage obtained with the parameters specified in Table 6.4. Only the adobe walls are shown here; the ring concrete beams and the lintel are hidden. The maximum top displacement reached was around 9.34 mm. Similar to the experimental response (Figure 6.13), the numerical results show a diagonal crack forming from the corners of the opening. Horizontal cracks are also detected in the perpendicular walls. a) Top displacement= 1 mm b) Top displacement= 4 mm Figure Damage evolution of the wall subjected to a horizontal top load (4 displacement levels). The top and bottom concrete beam and the lintel are hidden. Total-strain model, Midas FEA.

21 Numerical modelling of the seismic behaviour of adobe buildings 133 c) Top displacement= 6 mm d) Top displacement= 9.34 mm Figure Continuation. Damage evolution of the wall subjected to a horizontal top load (4 displacement levels). The top and bottom concrete beam and the lintel are hidden. Total-strain model, Midas FEA. Figure 6.26 shows the evolution of the maximum principal stresses on the in-plane adobe wall (the concrete beams, lintel and perpendicular walls are not shown). In the legend the maximum tensile stress value is limited to.5 MPa. A good agreement is seen between the maximum tensile stresses and the experimental damage pattern (Figure 6.14). Large tensile stresses correspond to large crack openings. The principal stresses reach their maximum values at the opening corners and then travel to the wall corners. The white zones inside the adobe wall show the parts where the tensile strength has been exceeded; this is possible when dealing with a fixed crack [de Borst and Nauta 1985; Feenstra and Rots 21; Noghabai 1999]. The good global agreement between the numerical results and the experimental tests leads to the conclusion that the total-strain model can be successfully applied to the analysis of adobe masonry, and the assumption of a homogeneous material is reasonable. Figure 6.27 shows the deformation pattern at the last computation step. The models were run in Midas FEA with arc-length method with initial stiffness. The number of load steps was specified as 1, the initial load factor was.1, and the maximum number of iterations per load step was 8. The convergence criteria were given by an energy norm equal to.1 and a displacement norm equal to.5.

22 134 Sabino Nicola Tarque Ruíz a) Top displacement= 1 mm b) Top displacement= 4 mm c) Top displacement= 6 mm d) Top displacement= 9.34 mm Figure Evolution of maximum principal stresses in the adobe wall subjected to a horizontal top displacement. Total-strain model, Midas FEA. Figure Deformation of the adobe wall due to a maximum horizontal top displacement of 9.34 mm. Total-strain model, Midas FEA.

23 Numerical modelling of the seismic behaviour of adobe buildings CONCRETE DAMAGED PLASTICITY: MODELLING THE PUSHOVER RESPONSE. A finite element model was created in Abaqus/Standard using 4-node rectangular shell elements without integration reduction (Figure 6.28a). For element controls, a finite membrane strain and a default drilling hourglass scaling factors were selected. The numerical model is similar to the one created in Section 6.3 with Midas FEA, so the shell elements are 1 x 1 mm with 3 mm thick, and the characteristic length is equal to the diagonal of the shell element. The reinforced concrete beams (top and bottom) and the wooden lintel are modelled using linear material properties. The adobe masonry is represented by the concrete damaged plasticity model, which takes into account the tension and compression constitutive laws for adobe. The displacement history is applied at one edge of the top concrete beam as seen in Figure 6.28b. The base of the foundation is fully fixed, while the top part of the crown concrete beam is free of movement. a) Complete view of the model b) Position of horizontal applied loads at the top beam Figure Finite element model of the adobe wall subjected to horizontal displacement loads at the top. Concrete Damaged Plasticity model, Abaqus/Standard. The material parameters used in Abaqus/Standard are essentially the ones used for Midas FEA, though the compression strength is increased to.45 MPa as shown in Table 6.5. Table 6.5. Material properties for the adobe masonry within concrete damaged plasticity model. E (N/mm 2 ) Elastic Tension Compression m (N/mm 3 ) h (mm) f t (N/mm 2 ) I G f (N/mm) f c (N/mm 2 ) c G f (N/mm) p (mm/mm) 2.2 2e

24 136 Sabino Nicola Tarque Ruíz The following default additional parameters are required for the concrete damaged plasticity model: dilatation angle= 1, eccentricity=.1, ratio of initial equibiaxial compressive yield stress to initial uniaxial compressive yield stress= 1.16, k parameter related to yield surface= 2/3, and null viscosity parameter. Again, a parametric study for selection of the tensile and compression strength is carried out. The tensile strength is between.2 to 6 MPa and the compression strength is varied between.3 to.8 MPa (see Figure 6.18 and Figure 6.23). As in the previous analysis, it is seen from the numerical pushover curves that the tensile strength is the parameter that controls the global behaviour of the adobe masonry; low values of tensile strength allows to a fast inelastic excursion and will end with convergence problems, high values of tensile strength make more brittle the adobe masonry but increase its lateral strength. According to Figure 6.29 a value of f t equal to.4 MPa should be selected Force (kn) ft=.2 MPa ft=.4 MPa ft=.6 MPa Figure Influence of tensile strength model, Abaqus/Standard. f t on the pushover response. Concrete damaged plasticity Figure 6.3 shows the numerical pushover curves analyzed with different compression strength values. It is seen that the compression strength increment influences on the maximum lateral strength of the adobe wall, but it maintains similar post peak behaviour and failure pattern. The main difference in lateral strength is seen from 2 to 4 mm of top displacement and it is due to the biaxial interaction between tensile and compressive strength. Less difference it is seen for the pushover curves computed with f c =.45 MPa to.8 MPa. A lower bound of f c =.3 MPa can be considered without loss of accuracy on the global response for the plasticity damage model implemented in Abaqus. In this case, there was a not convergence problem as seen in Midas FEA.

25 Numerical modelling of the seismic behaviour of adobe buildings Force (kn) fc=.3 MPa fc=.45 MPa fc=.7 MPa fc=.8 MPa Figure 6.3. Influence of compression strength f c on the pushover response. Concrete damaged plasticity model, Abaqus/Standard. For the cyclic analysis done in Section 6.5 a compression strength of.3 MPa is assumed for the adobe masonry, while for the dynamic analysis performed in Section 7 the compression strength was increased to.45 MPa. In both cases the relation G C f / fc is maintained as.344 mm Results of the pushover response considering the concrete damaged plasticity model The displacement pattern at the last stage is shown in Figure The global behaviour of the numerical analysis on the adobe wall represents well the experimental test in terms of crack pattern and lateral capacity (see Figure 6.14). It is preliminarily concluded that the calibrated material properties can be used for further numerical analyses. Figure Deformation of the adobe wall due to a maximum horizontal top displacement of 1 mm. Concrete damaged plasticity model, Abaqus/Standard.

26 138 Sabino Nicola Tarque Ruíz Furthermore, the damage pattern is analyzed based on the formation of plastic strain in tension at different levels of top displacement (Figure 6.32). It is observed that the formation of cracks starts at the opening corners and at the contact zone of the lintel with the adobe masonry. Horizontal cracks are also observed at the perpendicular walls, similarly to the ones observed in the experimental test. a) Top displacement= 1 mm b) Top displacement= 4 mm c) Top displacement= 6 mm d) Top displacement= 1 mm Figure Evolution of maximum in-plane plastic strain in the adobe wall subjected to a horizontal top displacement. Concrete damaged plasticity model, Abaqus/Standard. Figure 6.33 shows the evolution of the maximum principal stresses in the in-plane adobe wall, without the two concrete beams, the timber lintel and the perpendicular walls. It is observed that the maximum tensile zones (shown in red) are reached first at the opening corners and evolve diagonally to the wall corners. After any integration point reaches f t, so the tensile stress value descends but increasing the crack displacement (softening part of the tensile constitutive law). The models were run in Abaqus/Standard specifying a direct method -for equation solver- with full Newton solution technique. The total displacement load is applied in 1s, having a minimum increment size of.1 s with a maximum of.5 s. The maximum number of increments is 2. Non linear geometric effects are considered for the

27 Numerical modelling of the seismic behaviour of adobe buildings 139 analysis of equilibrium even though they are not expected to affect the results of such a stiff wall. Two control parameters are also specified, the automatic stabilization for a dissipated energy fraction of.1, and the adaptive stabilization with maximum ratio of stabilization to strain energy of.1. a) Top displacement= 1 mm b) Top displacement= 4 mm c) Top displacement= 6 mm d) Top displacement= 1 mm Figure Evolution of maximum principal stresses in the adobe wall subjected to a horizontal top displacement. Concrete damaged plasticity model, Abaqus/Standard. 6.5 CONCRETE DAMAGED PLASTICITY: MODELLING THE CYCLIC BEHAVIOUR The finite element model used for the pushover analysis in Abaqus/Standard was then used for calibration of parameters for cyclic behaviour. The idea is to calibrate the damage factors d t and d c (which control the closing of cracking and reduces the elastic stiffness during unloading), and the stiffness recovery w t and w c, for tension and compression respectively, to reproduce the behaviour of adobe masonry under reversal loads. The tensile and compression strength are kept as.4 and.3 MPa, respectively. According to Abaqus 6.9 SIMULIA [29], it was seen from experimental tests on concrete that the compressive stiffness can be recovered upon crack closure as the load changes from tension to compression. On the other hand, the tensile stiffness is not

28 14 Sabino Nicola Tarque Ruíz recovered as the load changes from compression to tension once crushing micro-cracks have developed. This assumption is kept for the adobe masonry. The sequence of applied displacements consists two cycles for top displacement of 1 mm and one cycle for top displacements of 2, 5 and 1 mm, as shown in Figure In the experimental test, each displacement limit was repeated twice Time step (s) Figure History of static horizontal displacement load applied to the numerical model Calibration of stiffness recovery and damage factors Since no data tests exist for adobe masonry under reversal loads, stiffness recovery and damage factors were assumed. The idea was to numerically represent the closing of cracking during the change of tension to compression by selecting appropiate damage factors. The numerical force-displacement curves are compared with the results of the cyclic experimental test. By default, Abaqus assigns w c = 1 and w t =, which indicate full stiffness recovery when the integration point is under compression stress, and no stiffness recovery when it is subjected again to tensile stress (see Figure 5.18). Much attention was given to the tensile behaviour rather than to the compression one because it seems that the tension controls the global response of the in-plane adobe walls. As Table 6.6 to Table 6.16 show it is not possible to reach zero crack displacement when the load changes from tension to compression -also with large values of the tensile damage factors- so a residual crack deformation remains when the load is revearsed. Table 6.6. Proposed compression damage factor: Dc-1. Damage factor dc Plastic strain (mm/mm)

32 144 Sabino Nicola Tarque Ruíz problems. Models 18 to 21 maintain the tensile damage factor specified in Table 6.7, which are close to the ones specified in Table 6.9, but with a variation of the compression stiffness recovery from.5 to.8. The best match is obtained with Model 2, which considered wc.6, despite the large residual deformations Numerical model Numerical model 2 a) w c = 1, w t =, no damage factor b) w c = 1, 5 w t =.9, no damage factor Numerical model Numerical model c) w c =.75, w t =, no damage factor d) w c =.5, w t =, no damage factor Numerical model Numerical model e) w c = 1, w t =, d c =Table 6.6, d t =Table 6.7 f) w c =.5, w t =.5, d c =Table 6.6, d t =Table 6.8

33 Numerical modelling of the seismic behaviour of adobe buildings Numerical model 7 g) w c =.5, -4 w t =.5, d c = Table 6.6, d t = Table Numerical model Numerical model 9 h) w c = 1, w t =, d t = Table 6.7 i) w c = 1, w t =, d t = Table Numerical model Numerical model j) w c = 1, w t =, d t = Table 6.11 k) w c = 1, w t =, d t = Table 6.12 Figure Comparison of the experimental and numerical cyclic behaviour of the adobe wall taking into account variability in the recovery stiffness and damage factors, in tension and compression. Horizontal load applied on the left part of the top concrete beam.

34 146 Sabino Nicola Tarque Ruíz Numerical model 12 l) w c = 1, -5 w t =, d t = Table m) w c =.9, Numerical model 13 w t =, d t = Table 6.7 n) w c =.9, Numerical model 14 w t =, d t = Table Numerical model Numerical model o) w c =.9, w t =, d t = Table 6.14 p) w c =.8, w t =, d t = Table 6.4 Figure Continuation. Comparison of the experimental and numerical cyclic behaviour of the adobe wall taking into account variability in the recovery stiffness and damage factors, in tension and compression. Horizontal load applied on the left part of the top concrete beam.

35 Numerical modelling of the seismic behaviour of adobe buildings Numerical model 17 q) w c =.7, -5 w t =, d t = Table Numerical model Numerical model r) w c =.8, w t =, d t = Table 6.7 s) w c =.7, w t =, d t = Table Numerical model Numerical model 21 t) w c =.6, w t =, d t = Table 6.7 u) w c =.5, w t =, d t = Table 6.7 Figure Continuation. Comparison of the experimental and numerical cyclic behaviour of the adobe wall taking into account variability in the recovery stiffness and damage factors, in tension and compression. Horizontal load applied on the left part of the top concrete beam.

36 148 Sabino Nicola Tarque Ruíz A new cyclic analysis was performed in Abaqus considering a variation of the zones where the load displacements are applied. In the previous cases the load was applied at one vertical edge of the top concrete beam (Figure 6.28b), which can be good for a monotonic test but probably not good representative for a cyclic one. The same load pattern was later applied at both vertical edges of the top concrete beam and part of the adobe masonry, as shown in Figure 6.36, in order to simulate better the experimental test (see Figure 3.9). The results of the parametric study are shown in Figure 6.37 and demonstrate some improvements for the numerical results. Figure 6.36 Finite element model of the adobe wall considering both ends of the top concrete beam for application of the cyclic horizontal displacement. Concrete Damaged Plasticity model, Abaqus/Standard. Models 22 to 25 evaluate the variation of the tensile damage factors. Again, the needs to reduce the compression stiffness when the stress goes from tension to compression are observed to match the experimental curve, especially for the loading branch at 1 mm displacement. The new tensile factors specified in Table 6.15 and Table 6.16 do not show improvement in the reduction of residual deformations for unloading. Models 26 to 28 consider the tensile damage factors given in Table 6.7 and consider compression stiffness factors from.5 to.8. The best results are obtained with models 27 and 28, concluding that the compression stiffness factors w c should be specified between.5 and.6. The tensile damage factor can not be significantly different from those given in Table 6.7 or Table 6.9; otherwise, convergence problems may stop the analysis before the last stage, always in Abaqus.

37 Numerical modelling of the seismic behaviour of adobe buildings Numerical model Numerical model 23 a) w c = 1, w t =, d t = Table 6.7 b) w c = 1, w t =, d t = Table Numerical model Numerical model c) w c = 1, w t =, d t = Table 6.15 d) w c = 1, w t =, d t = Table Numerical model Numerical model 27 e) w c =.8, w t =, d t = Table 6.7 f) w c =.6, w t =, d t = Table 6.7 Figure Comparison of the experimental and numerical cyclic behaviour of the adobe wall taking into account variability in the recovery stiffness and damage factors, in tension and compression. Horizontal load applied at both ends of the top concrete beam.

38 15 Sabino Nicola Tarque Ruíz Numerical model 28 g) w c =.5, -5 w t =, d t = Table 6.7 Figure Continuation. Comparison of the experimental and numerical cyclic behaviour of the adobe wall taking into account variability in the recovery stiffness and damage factors. Horizontal load applied at both ends of the top concrete beam. The Model 28 (Figure 6.37g) is used for showing the cracking process. From the analysis of the plastic strain it is seen that after the first 2 cycles of 1 mm some regions of the adobe masonry already exceed the maximum elastic strain. This effect is seen at the opening corners, where a concentration of tensile stresses is expected to occur (Figure 6.38a). Figure 6.38 shows the formation process of tensile plastic strains at different values of the top displacement load in Model 28, where the most important aspect is the formation of the X-diagonal cracks, typical of the in-plane behaviour of masonry. The numerical results match the failure pattern seen in the experimental test (Figure 6.14). Horizontal cracks at the perpendicular walls are also formed due to bending. a) Plastic strain values at the end of the 2 cycles of 1 mm Figure Formation process of the tensile plastic strain on the adobe wall under cyclic loads. A non unique legend in placed each to each figure to visualize better the plastic strain. Concrete damaged plasticity model, Abaqus/Standard.

39 Numerical modelling of the seismic behaviour of adobe buildings 151 c) Plastic strain values at the end of the cycle of 5 mm d) Plastic strain values at the end of the cycle of 1 mm Figure Continuation. Formation process of the tensile plastic strain on the adobe wall under cyclic loads. A non unique legend in placed each to each figure to visualize better the plastic strain. Concrete damaged plasticity model, Abaqus/Standard. Another way for interpretation of the tensile damage occurred in the adobe masonry is to show the tensile damage factor (Figure 6.39). Figure Tensile damage factor for Model 28 at the end of the history of cyclic horizontal displacement load.

40 152 Sabino Nicola Tarque Ruíz The tensile damage factor is a non-decreasing quantity associated with the tensile failure of the material. In Figure 6.39, the zones which are not in blue ( d t = ) indicate the zones which already are in the softening part of the tensile constitutive law and can be interpreted as damage zones. 6.6 VIBRATION MODES An eigenvalue analysis is done to compute the vibration modes of the model, especially in the direction of the applied load (X-X direction). The reinforced concrete beam placed as foundation of the wall is removed, so the total weight of the model is 1.63 kn. The base of the wall is fully fixed. The analysis is performed with Abaqus/Standard through the linear perturbation option and considering Lanczos method for extraction of the frequency values. A 5% of the elasticity of modulus has been used according to Tarque [28] to take into account early cracking into the material. Figure 6.4 shows the effective mass related to the first 11 vibration modes, represented here by the frequency values. In theory the sum of all the effective masses should be equal to the total mass of the model. It is seen that 11 modes of vibration are required to reach the 9% of the total mass (Table 6.17), being the fundamental one the first mode. 8 7 Effective mass (%) Frequency (Hz) Figure 6.4. Contribution of the modes of vibration in the X-X direction until reaches the 9% of the total mass of the model. The deflected shapes given for each mode of vibration are shown in Figure The first vibration mode, which involves 74.24% of the total mass, is a translational mode; while the others are basically out-of-plane deformations of the flange walls. This analysis considers the use of the elastic material properties. However, the adobe material is brittle and goes into the inelastic range very early; therefore the frequencies are expected to shorten.

41 Numerical modelling of the seismic behaviour of adobe buildings 153 a) Mode 1. T 1 =.1157s b) Mode 2. T 2 =.533s c) Mode 3. T 3 =.391s d) Mode 4. T 4 =.344s e) Mode 5. T 5 =.275s f) Mode 6. T 6 =.242s Figure Modes of vibration in the X-X direction.

### Structural behaviour of traditional mortise-and-tenon timber joints

Structural behaviour of traditional mortise-and-tenon timber joints Artur O. Feio 1, Paulo B. Lourenço 2 and José S. Machado 3 1 CCR Construtora S.A., Portugal University Lusíada, Portugal 2 University

9th International Conference on Fracture Mechanics of Concrete and Concrete Structures FraMCoS-9 Chenjie Yu, P.C.J. Hoogenboom and J.G. Rots DOI 10.21012/FC9.288 ALGORITHM FOR NON-PROPORTIONAL LOADING

### Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling

Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling F. Portioli, L. Cascini, R. Landolfo University of Naples Federico II, Italy P. Foraboschi IUAV University,

### Abstract. 1 Introduction

Contact analysis for the modelling of anchors in concrete structures H. Walter*, L. Baillet** & M. Brunet* *Laboratoire de Mecanique des Solides **Laboratoire de Mecanique des Contacts-CNRS UMR 5514 Institut

### MASONRY MICRO-MODELLING ADOPTING A DISCONTINUOUS FRAMEWORK

MASONRY MICRO-MODELLING ADOPTING A DISCONTINUOUS FRAMEWORK J. Pina-Henriques and Paulo B. Lourenço School of Engineering, University of Minho, Guimarães, Portugal Abstract Several continuous and discontinuous

### Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models

Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models I. Rhee, K.J. Willam, B.P. Shing, University of Colorado at Boulder ABSTRACT: This paper examines the global

### NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS

NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS ABSTRACT : P Mata1, AH Barbat1, S Oller1, R Boroschek2 1 Technical University of Catalonia, Civil Engineering

### Fluid driven cohesive crack propagation in quasi-brittle materials

Fluid driven cohesive crack propagation in quasi-brittle materials F. Barpi 1, S. Valente 2 Department of Structural and Geotechnical Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129

### University of Sheffield The development of finite elements for 3D structural analysis in fire

The development of finite elements for 3D structural analysis in fire Chaoming Yu, I. W. Burgess, Z. Huang, R. J. Plank Department of Civil and Structural Engineering StiFF 05/09/2006 3D composite structures

### EDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d

Institute of Industrial Science, University of Tokyo Bulletin of ERS, No. 48 (5) A TWO-PHASE SIMPLIFIED COLLAPSE ANALYSIS OF RC BUILDINGS PHASE : SPRING NETWORK PHASE Shanthanu RAJASEKHARAN, Muneyoshi

### Numerical Modelling of Blockwork Prisms Tested in Compression Using Finite Element Method with Interface Behaviour

13 th International Brick and Block Masonry Conference Amsterdam, July 4-7, 2004 Numerical Modelling of Blockwork Prisms Tested in Compression Using Finite Element Method with Interface Behaviour H. R.

Essentially, all models are wrong, but some are useful George E.P. Box, (Box and Draper, 1987) 5 ADVANCED FRACTURE MODELS In the previous chapter it was shown that the MOR parameter cannot be relied upon

### Chapter Two: Mechanical Properties of materials

Chapter Two: Mechanical Properties of materials Time : 16 Hours An important consideration in the choice of a material is the way it behave when subjected to force. The mechanical properties of a material

### Introduction and Background

Introduction and Background Itasca Consulting Group, Inc. (Itasca) has been participating in the geomechanical design of the underground 118-Zone at the Capstone Minto Mine (Minto) in the Yukon, in northwestern

### Modelling the nonlinear shear stress-strain response of glass fibrereinforced composites. Part II: Model development and finite element simulations

Modelling the nonlinear shear stress-strain response of glass fibrereinforced composites. Part II: Model development and finite element simulations W. Van Paepegem *, I. De Baere and J. Degrieck Ghent

### Seismic Response Analysis of Structure Supported by Piles Subjected to Very Large Earthquake Based on 3D-FEM

Seismic Response Analysis of Structure Supported by Piles Subjected to Very Large Earthquake Based on 3D-FEM *Hisatoshi Kashiwa 1) and Yuji Miyamoto 2) 1), 2) Dept. of Architectural Engineering Division

### THE BEHAVIOUR OF REINFORCED CONCRETE AS DEPICTED IN FINITE ELEMENT ANALYSIS.

THE BEHAVIOUR OF REINFORCED CONCRETE AS DEPICTED IN FINITE ELEMENT ANALYSIS. THE CASE OF A TERRACE UNIT. John N Karadelis 1. INTRODUCTION. Aim to replicate the behaviour of reinforced concrete in a multi-scale

### Heterogeneous structures studied by interphase elasto-damaging model.

Heterogeneous structures studied by interphase elasto-damaging model. Giuseppe Fileccia Scimemi 1, Giuseppe Giambanco 1, Antonino Spada 1 1 Department of Civil, Environmental and Aerospace Engineering,

### NON-LINEAR STATIC ANALYSIS OF MASONRY STRUCTURES

Non-linear static analysis of masonry structures XIII International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XIII E. Oñate, D.R.J. Owen, D. Peric & M. Chiumenti (Eds)

### PREDICTION OF THE CYCLIC BEHAVIOR OF MOMENT RESISTANT BEAM-TO-COLUMN JOINTS OF COMPOSITE STRUCTURAL ELEMENTS

SDSS Rio 21 STABILITY AND DUCTILITY OF STEEL STRUCTURES E. Batista, P. Vellasco, L. de Lima (Eds.) Rio de Janeiro, Brazil, September 8-1, 21 PREDICTION OF THE CYCLIC BEHAVIOR OF MOMENT RESISTANT BEAM-TO-COLUMN

### A Nonlinear Interface Element for 3D Mesoscale Analysis. of Brick-Masonry Structures

A Nonlinear Interface Element for 3D Mesoscale Analysis of Brick-Masonry Structures L. Macorini 1, B.A. Izzuddin 2 Abstract This paper presents a novel interface element for the geometric and material

### Nonlinear static analysis PUSHOVER

Nonlinear static analysis PUSHOVER Adrian DOGARIU European Erasmus Mundus Master Course Sustainable Constructions under Natural Hazards and Catastrophic Events 520121-1-2011-1-CZ-ERA MUNDUS-EMMC Structural

### FINITE ELEMENT ANALYSIS OF TAPERED COMPOSITE PLATE GIRDER WITH A NON-LINEAR VARYING WEB DEPTH

Journal of Engineering Science and Technology Vol. 12, No. 11 (2017) 2839-2854 School of Engineering, Taylor s University FINITE ELEMENT ANALYSIS OF TAPERED COMPOSITE PLATE GIRDER WITH A NON-LINEAR VARYING

### Discrete Element Modelling of a Reinforced Concrete Structure

Discrete Element Modelling of a Reinforced Concrete Structure S. Hentz, L. Daudeville, F.-V. Donzé Laboratoire Sols, Solides, Structures, Domaine Universitaire, BP 38041 Grenoble Cedex 9 France sebastian.hentz@inpg.fr

### Finite element analysis of diagonal tension failure in RC beams

Finite element analysis of diagonal tension failure in RC beams T. Hasegawa Institute of Technology, Shimizu Corporation, Tokyo, Japan ABSTRACT: Finite element analysis of diagonal tension failure in a

### A Performance Modeling Strategy based on Multifiber Beams to Estimate Crack Openings ESTIMATE in Concrete Structures CRACK

A Performance Modeling Strategy based on Multifiber Beams to Estimate Crack Openings ESTIMATE in Concrete Structures CRACK A. Medjahed, M. Matallah, S. Ghezali, M. Djafour RiSAM, RisK Assessment and Management,

### Non-linear and time-dependent material models in Mentat & MARC. Tutorial with Background and Exercises

Non-linear and time-dependent material models in Mentat & MARC Tutorial with Background and Exercises Eindhoven University of Technology Department of Mechanical Engineering Piet Schreurs July 7, 2009

### 3D Finite Element analysis of stud anchors with large head and embedment depth

3D Finite Element analysis of stud anchors with large head and embedment depth G. Periškić, J. Ožbolt & R. Eligehausen Institute for Construction Materials, University of Stuttgart, Stuttgart, Germany

### EQUIVALENT FRACTURE ENERGY CONCEPT FOR DYNAMIC RESPONSE ANALYSIS OF PROTOTYPE RC GIRDERS

EQUIVALENT FRACTURE ENERGY CONCEPT FOR DYNAMIC RESPONSE ANALYSIS OF PROTOTYPE RC GIRDERS Abdul Qadir Bhatti 1, Norimitsu Kishi 2 and Khaliq U Rehman Shad 3 1 Assistant Professor, Dept. of Structural Engineering,

### Nonlinear Analysis of Reinforced Concrete Shells Subjected to Impact Loads

Transactions of the 7 th International Conference on Structural Mechanics in Reactor Technology (SMiRT 7) Prague, Czech Republic, August 7, 00 Paper # J0- Nonlinear Analysis of Reinforced Concrete Shells

### Centrifuge Shaking Table Tests and FEM Analyses of RC Pile Foundation and Underground Structure

Centrifuge Shaking Table s and FEM Analyses of RC Pile Foundation and Underground Structure Kenji Yonezawa Obayashi Corporation, Tokyo, Japan. Takuya Anabuki Obayashi Corporation, Tokyo, Japan. Shunichi

### Cracked concrete structures under cyclic load

Cracked concrete structures under cyclic load Fabrizio Barpi & Silvio Valente Department of Structural and Geotechnical Engineering, Politecnico di Torino, Torino, Italy ABSTRACT: The safety of cracked

### Modelling Localisation and Spatial Scaling of Constitutive Behaviour: a Kinematically Enriched Continuum Approach

Modelling Localisation and Spatial Scaling of Constitutive Behaviour: a Kinematically Enriched Continuum Approach Giang Dinh Nguyen, Chi Thanh Nguyen, Vinh Phu Nguyen School of Civil, Environmental and

### Using the Abaqus CDP Model in Impact Simulations

SAFIR project (http://safir2018.vtt.fi/) The Finnish Research Programme on Nuclear Power Plant Safety 2015-2018 Using the Abaqus CDP Model in Impact Simulations Alexis Fedoroff Technical reseach centre

### Mechanics of Materials

Mechanics of Materials Notation: a = acceleration = area (net = with holes, bearing = in contact, etc...) SD = allowable stress design d = diameter of a hole = calculus symbol for differentiation e = change

### Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics

Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics Tobias Gasch *1 and Anders Ansell 1 1 KTH Royal Institute of Technology, Department of Civil and Architectural Engineering *Corresponding

### Numerical Analysis of Composite Panels in the Post-Buckling Field taking into account Progressive Failure

Copyright c 007 ICCES ICCES, vol.1, no.3, pp.93-98, 007 Numerical Analysis of Composite Panels in the Post-Buckling Field taking into account Progressive Failure C. Bisagni 1 Summary The research here

### NUMERICAL SIMULATIONS OF CORNERS IN RC FRAMES USING STRUT-AND-TIE METHOD AND CDP MODEL

Numerical simulations of corners in RC frames using Strut-and-Tie Method and CDP model XIII International Conference on Computational Plasticity. Fundamentals and Applications COMPLAS XIII E. Oñate, D.R.J.

### MESOSCOPIC MODELLING OF MASONRY USING GFEM: A COMPARISON OF STRONG AND WEAK DISCONTINUITY MODELS B. Vandoren 1,2, K. De Proft 2

Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm MESOSCOPIC MODELLING OF MASONRY USING GFEM: A COMPARISON OF STRONG AND WEAK DISCONTINUITY

### A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials

Dublin, October 2010 A Constitutive Framework for the Numerical Analysis of Organic Soils and Directionally Dependent Materials FracMan Technology Group Dr Mark Cottrell Presentation Outline Some Physical

### VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS

The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.

### Impact and Crash Modeling of Composite Structures: A Challenge for Damage Mechanics

Impact and Crash Modeling of Composite Structures: A Challenge for Damage Mechanics Dr. A. Johnson DLR Dr. A. K. Pickett ESI GmbH EURO-PAM 99 Impact and Crash Modelling of Composite Structures: A Challenge

### Bruno Silva 1, Athanasios Pappas 2, Maria R. Valluzzi 3, Francesca da Porto 4,Claudio Modena 5 ABSTRACT

Structural Analysis of Historical Constructions Jerzy Jasieńko (ed) 2012 DWE, Wrocław, Poland, ISSN 08602395, ISBN 9788371252167 CALIBRATION OF A NUMERICAL MATERIAL BEHAVIOUR MODEL FOR THE SIMULATION OF

### Course in. Geometric nonlinearity. Nonlinear FEM. Computational Mechanics, AAU, Esbjerg

Course in Nonlinear FEM Geometric nonlinearity Nonlinear FEM Outline Lecture 1 Introduction Lecture 2 Geometric nonlinearity Lecture 3 Material nonlinearity Lecture 4 Material nonlinearity it continued

Constitutive Modeling of Reinforced Concrete Panel Behavior under Cyclic Loading K. Orakcal & D. Ulugtekin Bogazici University, Istanbul, Turkey L. M. Massone University of Chile, Santiago, Chile SUMMARY:

### Sensitivity and Reliability Analysis of Nonlinear Frame Structures

Sensitivity and Reliability Analysis of Nonlinear Frame Structures Michael H. Scott Associate Professor School of Civil and Construction Engineering Applied Mathematics and Computation Seminar April 8,

### Module 5: Failure Criteria of Rock and Rock masses. Contents Hydrostatic compression Deviatoric compression

FAILURE CRITERIA OF ROCK AND ROCK MASSES Contents 5.1 Failure in rocks 5.1.1 Hydrostatic compression 5.1.2 Deviatoric compression 5.1.3 Effect of confining pressure 5.2 Failure modes in rocks 5.3 Complete

### An orthotropic damage model for crash simulation of composites

High Performance Structures and Materials III 511 An orthotropic damage model for crash simulation of composites W. Wang 1, F. H. M. Swartjes 1 & M. D. Gan 1 BU Automotive Centre of Lightweight Structures

### Theory at a Glance (for IES, GATE, PSU)

1. Stress and Strain Theory at a Glance (for IES, GATE, PSU) 1.1 Stress () When a material is subjected to an external force, a resisting force is set up within the component. The internal resistance force

### Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield Surface

Transaction A: Civil Engineering Vol. 16, No. 6, pp. 512{519 c Sharif University of Technology, December 2009 Research Note Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield

### EFFECT OF SHEAR REINFORCEMENT ON FAILURE MODE OF RC BRIDGE PIERS SUBJECTED TO STRONG EARTHQUAKE MOTIONS

EFFECT OF SHEAR REINFORCEMENT ON FAILURE MODE OF RC BRIDGE PIERS SUBJECTED TO STRONG EARTHQUAKE MOTIONS Atsuhiko MACHIDA And Khairy H ABDELKAREEM SUMMARY Nonlinear D FEM was utilized to carry out inelastic

### Card Variable MID RO E PR ECC QH0 FT FC. Type A8 F F F F F F F. Default none none none 0.2 AUTO 0.3 none none

Note: This is an extended description of MAT_273 input provided by Peter Grassl It contains additional guidance on the choice of input parameters beyond the description in the official LS-DYNA manual Last

Girão-Coelho AM. Finite element guidelines for simulation of delamination dominated failures in composite materials validated by case studies. Archives of Computational Methods in Engineering 2016, 1-26.

### Numerical simulations of a full scale load test on a stone masonry arch bridge

Historical Constructions, P.B. Lourenço, P. Roca (Eds.), Guimarães, 21 739 Numerical simulations of a full scale load test on a stone masonry arch bridge J. Alfaiate Prof. Auxiliar, Instituto Superior

### Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics

Cracking in Quasi-Brittle Materials Using Isotropic Damage Mechanics Tobias Gasch, PhD Student Co-author: Prof. Anders Ansell Comsol Conference 2016 Munich 2016-10-12 Contents Introduction Isotropic damage

### Finite Element Modelling with Plastic Hinges

01/02/2016 Marco Donà Finite Element Modelling with Plastic Hinges 1 Plastic hinge approach A plastic hinge represents a concentrated post-yield behaviour in one or more degrees of freedom. Hinges only

### A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE

A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE Daxu Zhang and D. R. Hayhurst School of Mechanical, Aerospace and Civil Engineering,

### Pillar strength estimates for foliated and inclined pillars in schistose material

Pillar strength estimates for foliated and inclined pillars in schistose material L.J. Lorig Itasca Consulting Group, Inc., Minneapolis, MN, USA A. Cabrera Itasca S.A., Santiago, Chile ABSTRACT: Pillar

### A 3D Discrete Damage Modeling Methodology for Abaqus for Fatigue Damage Evaluation in Bolted Composite Joints

A 3D Discrete Damage Modeling Methodology for Abaqus for Fatigue Damage Evaluation in Bolted Composite Joints Eugene Fang 1, Ling Liu 1, Michael Stuebner 1, Jim Lua 1 1 Global Engineering and Materials,

### ENERGY DIAGRAM w/ HYSTERETIC

ENERGY DIAGRAM ENERGY DIAGRAM w/ HYSTERETIC IMPLIED NONLINEAR BEHAVIOR STEEL STRESS STRAIN RELATIONSHIPS INELASTIC WORK DONE HYSTERETIC BEHAVIOR MOMENT ROTATION RELATIONSHIP IDEALIZED MOMENT ROTATION DUCTILITY

### Use of negative stiffness in failure analysis of concrete beams

Use of negative stiffness in failure analysis of concrete beams A. Salam Al-Sabah Research Scientist, Debra F. Laefer Associate Professor* Urban Modelling Group, School of Civil Engineering, University

### Sabah Shawkat Cabinet of Structural Engineering Walls carrying vertical loads should be designed as columns. Basically walls are designed in

Sabah Shawkat Cabinet of Structural Engineering 17 3.6 Shear walls Walls carrying vertical loads should be designed as columns. Basically walls are designed in the same manner as columns, but there are

### Numerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation

Numerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation P. Sushma Ph D Scholar, Earthquake Engineering Research Center, IIIT Hyderabad, Gachbowli,

### Mechanics of Materials Primer

Mechanics of Materials rimer Notation: A = area (net = with holes, bearing = in contact, etc...) b = total width of material at a horizontal section d = diameter of a hole D = symbol for diameter E = modulus

### Nonlinear seismic analysis of masonry buildings

Erlenbach, September 12 th, 2013 Nonlinear seismic analysis of masonry buildings Andrea Penna andrea.penna@unipv.it Department of Civil Engineering and Architecture University of Pavia, Italy EUCENRE Foundation

### Seismic Pushover Analysis Using AASHTO Guide Specifications for LRFD Seismic Bridge Design

Seismic Pushover Analysis Using AASHTO Guide Specifications for LRFD Seismic Bridge Design Elmer E. Marx, Alaska Department of Transportation and Public Facilities Michael Keever, California Department

### A Modified Response Spectrum Analysis Procedure (MRSA) to Determine the Nonlinear Seismic Demands of Tall Buildings

Fawad A. Najam Pennung Warnitchai Asian Institute of Technology (AIT), Thailand Email: fawad.ahmed.najam@ait.ac.th A Modified Response Spectrum Analysis Procedure (MRSA) to Determine the Nonlinear Seismic

### Exercise: concepts from chapter 8

Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic

### Static and Time Dependent Failure of Fibre Reinforced Elastomeric Components. Salim Mirza Element Materials Technology Hitchin, UK

Static and Time Dependent Failure of Fibre Reinforced Elastomeric Components Salim Mirza Element Materials Technology Hitchin, UK Introduction Fibre reinforced elastomers are used in many applications,

### MODELING SLAB-COLUMN CONNECTIONS REINFORCED WITH GFRP UNDER LOCALIZED IMPACT

MODELING SLAB-COLUMN CONNECTIONS REINFORCED WITH GFRP UNDER LOCALIZED IMPACT QI ZHANG and AMGAD HUSSEIN Faculty of Engineering, Memorial University of Newfoundland St. John s, Newfoundland, Canada, A1B

### Mesoscopic Simulation of Failure of Mortar and Concrete by 3D RBSM

Journal of Advanced Concrete Technology Vol., No., 85-4, October 5 / Copyright 5 Japan Concrete Institute 85 Scientific paper Mesoscopic Simulation of Failure of Mortar and Concrete by D RBSM Kohei Nagai,

### FEA A Guide to Good Practice. What to expect when you re expecting FEA A guide to good practice

FEA A Guide to Good Practice What to expect when you re expecting FEA A guide to good practice 1. Background Finite Element Analysis (FEA) has transformed design procedures for engineers. Allowing more

### Nonlinear Analysis Of An EPDM Hydraulic Accumulator Bladder. Richard Kennison, Race-Tec

Nonlinear Analysis Of An EPDM Hydraulic Accumulator Bladder Richard Kennison, Race-Tec Agenda Race-Tec Overview Accumulator Experimental Testing Material Testing Numerical Analysis: 1. Linear Buckling

### A FINITE ELEMENT MODEL FOR SIZE EFFECT AND HETEROGENEITY IN CONCRETE STRUCTURES

A FINITE ELEMENT MODEL FOR SIZE EFFECT AND HETEROGENEITY IN CONCRETE STRUCTURES Roque Luiz Pitangueira 1 and Raul Rosas e Silva 2 1 Department of Structural Engineering -Federal University of Minas Gerais

### Role of Force Resultant Interaction on Ultra-High Performance Concrete

First International Interactive Symposium on UHPC 216 Role of Force Resultant Interaction on Ultra-High Performance Concrete Author(s) & Affiliation: Titchenda Chan (1), Kevin R. Mackie (1), and Jun Xia

### Practical methodology for inclusion of uplift and pore pressures in analysis of concrete dams

Practical methodology for inclusion of uplift and pore pressures in analysis of concrete dams Michael McKay 1 and Francisco Lopez 2 1 Dams Engineer, GHD Pty 2 Principal Dams/Structural Engineer, GHD Pty

### Settlement and Bearing Capacity of a Strip Footing. Nonlinear Analyses

Settlement and Bearing Capacity of a Strip Footing Nonlinear Analyses Outline 1 Description 2 Nonlinear Drained Analysis 2.1 Overview 2.2 Properties 2.3 Loads 2.4 Analysis Commands 2.5 Results 3 Nonlinear

### NORMAL STRESS. The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts.

NORMAL STRESS The simplest form of stress is normal stress/direct stress, which is the stress perpendicular to the surface on which it acts. σ = force/area = P/A where σ = the normal stress P = the centric

### 14. LS-DYNA Forum 2016

14. LS-DYNA Forum 2016 A Novel Approach to Model Laminated Glass R. Böhm, A. Haufe, A. Erhart DYNAmore GmbH Stuttgart 1 Content Introduction and Motivation Common approach to model laminated glass New

### On Nonlinear Buckling and Collapse Analysis using Riks Method

Visit the SIMULIA Resource Center for more customer examples. On Nonlinear Buckling and Collapse Analysis using Riks Method Mingxin Zhao, Ph.D. UOP, A Honeywell Company, 50 East Algonquin Road, Des Plaines,

### BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS

BIAXIAL STRENGTH INVESTIGATION OF CFRP COMPOSITE LAMINATES BY USING CRUCIFORM SPECIMENS H. Kumazawa and T. Takatoya Airframes and Structures Group, Japan Aerospace Exploration Agency 6-13-1, Ohsawa, Mitaka,

### Chapter 3. Load and Stress Analysis. Lecture Slides

Lecture Slides Chapter 3 Load and Stress Analysis 2015 by McGraw Hill Education. This is proprietary material solely for authorized instructor use. Not authorized for sale or distribution in any manner.

### strain appears only after the stress has reached a certain critical level, usually specied by a Rankine-type criterion in terms of the maximum princip

Nonlocal damage models: Practical aspects and open issues Milan Jirasek LSC-DGC, Swiss Federal Institute of Technology at Lausanne (EPFL), Switzerland Milan.Jirasek@ep.ch Abstract: The purpose of this

### DETERMINING THE STRESS PATTERN IN THE HH RAILROAD TIES DUE TO DYNAMIC LOADS 1

PERIODICA POLYTECHNICA SER. CIV. ENG. VOL. 46, NO. 1, PP. 125 148 (2002) DETERMINING THE STRESS PATTERN IN THE HH RAILROAD TIES DUE TO DYNAMIC LOADS 1 Nándor LIEGNER Department of Highway and Railway Engineering

### Effect of buttress on reduction of rock slope sliding along geological boundary

Paper No. 20 ISMS 2016 Effect of buttress on reduction of rock slope sliding along geological boundary Ryota MORIYA *, Daisuke FUKUDA, Jun-ichi KODAMA, Yoshiaki FUJII Faculty of Engineering, Hokkaido University,

### Displacement ductility demand and strength reduction factors for rocking structures

Earthquake Resistant Engineering Structures VI 9 Displacement ductility demand and strength reduction factors for rocking structures M. Trueb, Y. Belmouden & P. Lestuzzi ETHZ-Swiss Federal Institute of

### Supplementary Figures

Fracture Strength (GPa) Supplementary Figures a b 10 R=0.88 mm 1 0.1 Gordon et al Zhu et al Tang et al im et al 5 7 6 4 This work 5 50 500 Si Nanowire Diameter (nm) Supplementary Figure 1: (a) TEM image

### five Mechanics of Materials 1 ARCHITECTURAL STRUCTURES: FORM, BEHAVIOR, AND DESIGN DR. ANNE NICHOLS SUMMER 2017 lecture

ARCHITECTURAL STRUCTURES: FORM, BEHAVIOR, AND DESIGN DR. ANNE NICHOLS SUMMER 2017 lecture five mechanics www.carttalk.com of materials Mechanics of Materials 1 Mechanics of Materials MECHANICS MATERIALS

### five Mechanics of Materials 1 ARCHITECTURAL STRUCTURES: FORM, BEHAVIOR, AND DESIGN DR. ANNE NICHOLS SUMMER 2014 lecture

ARCHITECTURAL STRUCTURES: FORM, BEHAVIOR, AND DESIGN DR. ANNE NICHOLS SUMMER 2014 lecture five mechanics www.carttalk.com of materials Mechanics of Materials 1 Mechanics of Materials MECHANICS MATERIALS

### Recent Research on EPS Geofoam Seismic Buffers. Richard J. Bathurst and Saman Zarnani GeoEngineering Centre at Queen s-rmc Canada

Recent Research on EPS Geofoam Seismic Buffers Richard J. Bathurst and Saman Zarnani GeoEngineering Centre at Queen s-rmc Canada What is a wall (SEISMIC) buffer? A compressible inclusion placed between

### Computational Analysis for Composites

Computational Analysis for Composites Professor Johann Sienz and Dr. Tony Murmu Swansea University July, 011 The topics covered include: OUTLINE Overview of composites and their applications Micromechanics

### D : SOLID MECHANICS. Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown.

D : SOLID MECHANICS Q. 1 Q. 9 carry one mark each. Q.1 Find the force (in kn) in the member BH of the truss shown. Q.2 Consider the forces of magnitude F acting on the sides of the regular hexagon having

### Members Subjected to Torsional Loads

Members Subjected to Torsional Loads Torsion of circular shafts Definition of Torsion: Consider a shaft rigidly clamped at one end and twisted at the other end by a torque T = F.d applied in a plane perpendicular

### REPRESENTING MATRIX CRACKS THROUGH DECOMPOSITION OF THE DEFORMATION GRADIENT TENSOR IN CONTINUUM DAMAGE MECHANICS METHODS

20 th International Conference on Composite Materials Copenhagen, 19-24 th July 2015 REPRESENTING MATRIX CRACKS THROUGH DECOMPOSITION OF THE DEFORMATION GRADIENT TENSOR IN CONTINUUM DAMAGE MECHANICS METHODS

### Numerical investigation of EDZ development around a deep polymetallic ore mine

Paper No. 198 ISMS 2016 Numerical investigation of EDZ development around a deep polymetallic ore mine Mountaka Souley a *, Marwan Al Heib a, Vincent Renaud a a INERIS, c/o Ecole des Mines de Nancy, Campus

### Experimental and theoretical characterization of Li 2 TiO 3 and Li 4 SiO 4 pebbles

Experimental and theoretical characterization of Li 2 TiO 3 and Li 4 SiO 4 s D. Aquaro 1 N. Zaccari ABSTRACT Dipartimento di Ingegneria Meccanica Nucleare e della Produzione University of Pisa (Italy)