# 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.

### 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

### 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 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

### 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

### 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

### 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.

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:

### 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.

### 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

### 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

### 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,

### 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

### 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

### 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

### 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,

### 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,

### 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

### 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

### 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

### 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

### 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

### 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)

### Parametric analysis and torsion design charts for axially restrained RC beams

Structural Engineering and Mechanics, Vol. 55, No. 1 (2015) 1-27 DOI: http://dx.doi.org/10.12989/sem.2015.55.1.001 1 Parametric analysis and torsion design charts for axially restrained RC beams Luís F.A.

### INFLUENCE OF WEB THICKNESS REDUCTION IN THE SHEAR RESISTANCE OF NON-PRISMATIC TAPERED PLATE GIRDERS

INFLUENCE OF WEB THICKNESS REDUCTION IN THE SHEAR RESISTANCE OF NON-PRISMATIC TAPERED PLATE GIRDERS Paulo J. S. Cruz 1, Lúcio Lourenço 1, Hélder Quintela 2 and Manuel F. Santos 2 1 Department of Civil

### 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

### Chapter 12. Static Equilibrium and Elasticity

Chapter 12 Static Equilibrium and Elasticity Static Equilibrium Equilibrium implies that the object moves with both constant velocity and constant angular velocity relative to an observer in an inertial

### 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

### Earth Pressure Theory

Lateral Earth Pressure Page 1 Earth Pressure Theory Examples of Retaining Walls Lateral Earth Pressure Page 2 At-Rest, Active and Passive Earth Pressure Wednesday, August 17, 2011 12:45 PM At-rest condition

### THEME A. Analysis of the elastic behaviour of La Aceña arch-gravity dam

THEME A Analysis of the elastic behaviour of La Aceña arch-gravity dam Gjorgi KOKALANOV, Professor, Faculty of Civil Eng., Skopje, Republic of Macedonia Ljubomir TANČEV, Professor, Faculty of Civil Eng.,

### COMPARISON BETWEEN 2D AND 3D ANALYSES OF SEISMIC STABILITY OF DETACHED BLOCKS IN AN ARCH DAM

COMPARISON BETWEEN 2D AND 3D ANALYSES OF SEISMIC STABILITY OF DETACHED BLOCKS IN AN ARCH DAM Sujan MALLA 1 ABSTRACT The seismic safety of the 147 m high Gigerwald arch dam in Switzerland was assessed for

### Basics of Finite Element Analysis. Strength of Materials, Solid Mechanics

Basics of Finite Element Analysis Strength of Materials, Solid Mechanics 1 Outline of Presentation Basic concepts in mathematics Analogies and applications Approximations to Actual Applications Improvisation

### Some recent developments in computational modelling of concrete fracture

International Journal of Fracture 86: 5 36, 1997. c 1997 Kluwer Academic Publishers. Printed in the Netherlands. Some recent developments in computational modelling of concrete fracture RENÉ DE BORST Delft

### Brittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm

Lecture 6 Brittle Deformation Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm WW Norton, unless noted otherwise Brittle deformation EarthStructure (2 nd

### Predicting Fatigue Life with ANSYS Workbench

Predicting Fatigue Life with ANSYS Workbench How To Design Products That Meet Their Intended Design Life Requirements Raymond L. Browell, P. E. Product Manager New Technologies ANSYS, Inc. Al Hancq Development

### AN IMPORTANT PITFALL OF PSEUDO-STATIC FINITE ELEMENT ANALYSIS

AN IMPORTANT PITFALL OF PSEUDO-STATIC FINITE ELEMENT ANALYSIS S. Kontoe, L. Pelecanos & D.M. Potts ABSTRACT: Finite Element (FE) pseudo-static analysis can provide a good compromise between simplified

### Modelling Progressive Failure with MPM

Modelling Progressive Failure with MPM A. Yerro, E. Alonso & N. Pinyol Department of Geotechnical Engineering and Geosciences, UPC, Barcelona, Spain ABSTRACT: In this work, the progressive failure phenomenon

### After lecture 16 you should be able to

Lecture 16: Design of paper and board packaging Advanced concepts: FEM, Fracture Mechanics After lecture 16 you should be able to describe the finite element method and its use for paper- based industry

### Numerical Modeling of Direct Shear Tests on Sandy Clay

Numerical Modeling of Direct Shear Tests on Sandy Clay R. Ziaie Moayed, S. Tamassoki, and E. Izadi Abstract Investigation of sandy clay behavior is important since urban development demands mean that sandy

### ENGN 2290: Plasticity Computational plasticity in Abaqus

ENGN 229: Plasticity Computational plasticity in Abaqus The purpose of these exercises is to build a familiarity with using user-material subroutines (UMATs) in Abaqus/Standard. Abaqus/Standard is a finite-element

### MODELING OF NONLINEAR BEHAVIOR OF RC SHEAR WALLS UNDER COMBINED AXIAL, SHEAR AND FLEXURAL LOADING

CD02-003 MODELING OF NONLINEAR BEHAVIOR OF RC SHEAR WALLS UNDER COMBINED AXIAL, SHEAR AND FLEXURAL LOADING B. Ghiassi 1, M. Soltani 2, A. A. Tasnimi 3 1 M.Sc. Student, School of Engineering, Tarbiat Modares

### Flexure: Behavior and Nominal Strength of Beam Sections

4 5000 4000 (increased d ) (increased f (increased A s or f y ) c or b) Flexure: Behavior and Nominal Strength of Beam Sections Moment (kip-in.) 3000 2000 1000 0 0 (basic) (A s 0.5A s ) 0.0005 0.001 0.0015

### EMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading

MA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science Zhe Cheng (2018) 2 Stress & Strain - Axial Loading Statics

### Numerical Characterization of Concrete Heterogeneity

Vol. Materials 5, No. Research, 3, 2002Vol. 5, No. 3, Statistical 309-314, 2002. Characterization of the Concrete Numerical Modeling of Size Effect In Heterogeneity 2002 309 Numerical Characterization

### Modeling the bond of GFRP and concrete based on a damage evolution approach

Modeling the ond of GFRP and concrete ased on a damage evolution approach Mohammadali Rezazadeh 1, Valter Carvelli 2, and Ana Veljkovic 3 1 Dep. Architecture, Built environment and Construction engineering,

### Model tests and FE-modelling of dynamic soil-structure interaction

Shock and Vibration 19 (2012) 1061 1069 1061 DOI 10.3233/SAV-2012-0712 IOS Press Model tests and FE-modelling of dynamic soil-structure interaction N. Kodama a, * and K. Komiya b a Waseda Institute for

### Tuesday, February 11, Chapter 3. Load and Stress Analysis. Dr. Mohammad Suliman Abuhaiba, PE

1 Chapter 3 Load and Stress Analysis 2 Chapter Outline Equilibrium & Free-Body Diagrams Shear Force and Bending Moments in Beams Singularity Functions Stress Cartesian Stress Components Mohr s Circle for

### Mechanics of Solids. Mechanics Of Solids. Suraj kr. Ray Department of Civil Engineering

Mechanics Of Solids Suraj kr. Ray (surajjj2445@gmail.com) Department of Civil Engineering 1 Mechanics of Solids is a branch of applied mechanics that deals with the behaviour of solid bodies subjected

### Análisis Computacional del Comportamiento de Falla de Hormigón Reforzado con Fibras Metálicas

San Miguel de Tucuman, Argentina September 14 th, 2011 Seminary on Análisis Computacional del Comportamiento de Falla de Hormigón Reforzado con Fibras Metálicas Antonio Caggiano 1, Guillermo Etse 2, Enzo

### Numerical simulation of delamination onset and growth in laminated composites

Numerical simulation of delamination onset and growth in laminated composites G. Wimmer, C. Schuecker, H.E. Pettermann Austrian Aeronautics Research (AAR) / Network for Materials and Engineering at the

### The design of a formula student front impact attenuator

The design of a formula student front impact attenuator J.M.J. Schormans MT 10.06 Supervisors: dr. ir. Varvara Kouznetsova Eindhoven University of Technology Department of Mechanical Engineering Computational

### Modeling of Cyclic Load-Deformation Behavior of Shallow Foundations Supporting Rocking Shear Walls. Sivapalan Gajan. Advisor: Bruce Kutter

Modeling of Cyclic Load-Deformation Behavior of Shallow Foundations Supporting Rocking Shear Walls Sivapalan Gajan Advisor: Bruce Kutter Seminar 06. 01. 2005 Overview of Presentation Background Experimental

### Civil Engineering Design (1) Design of Reinforced Concrete Columns 2006/7

Civil Engineering Design (1) Design of Reinforced Concrete Columns 2006/7 Dr. Colin Caprani, Chartered Engineer 1 Contents 1. Introduction... 3 1.1 Background... 3 1.2 Failure Modes... 5 1.3 Design Aspects...

### Introduction to Engineering Materials ENGR2000. Dr. Coates

Introduction to Engineering Materials ENGR2 Chapter 6: Mechanical Properties of Metals Dr. Coates 6.2 Concepts of Stress and Strain tension compression shear torsion Tension Tests The specimen is deformed

### ANSYS Mechanical Basic Structural Nonlinearities

Lecture 4 Rate Independent Plasticity ANSYS Mechanical Basic Structural Nonlinearities 1 Chapter Overview The following will be covered in this Chapter: A. Background Elasticity/Plasticity B. Yield Criteria

### Experimental Study and Numerical Simulation on Steel Plate Girders With Deep Section

6 th International Conference on Advances in Experimental Structural Engineering 11 th International Workshop on Advanced Smart Materials and Smart Structures Technology August 1-2, 2015, University of

### Finite element analyses for limit state design of concrete structures

Finite element analyses for limit state design of concrete structures Posibilities and limitations with SOLVIA Masther s Thesis in the International Master s programme Structural Engineering FERNANDO LOPEZ

### 3-D FINITE ELEMENT CYCLIC ANALYSIS OF RC BEAM/COLUMN JOINT USING SPECIAL BOND MODEL

13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 24 Paper No. 446 3-D FINITE ELEMENT CYCLIC ANALYSIS OF RC BEAM/COLUMN JOINT USING SPECIAL BOND MODEL Kazuki TAJIMA 1,

### Example-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium

Example-3 Title Cylindrical Hole in an Infinite Mohr-Coulomb Medium Description The problem concerns the determination of stresses and displacements for the case of a cylindrical hole in an infinite elasto-plastic

### Transactions on Modelling and Simulation vol 10, 1995 WIT Press, ISSN X

Parameters controlling the numerical simulation validity of damageable composite toughness testing S. Yotte, C. Currit, E. Lacoste, J.M. Quenisset Laboratoire de Genie Meanique - IUT 'A\ Domaine Universitaire,

### Lateral Earth Pressure

1 of 11 6/2/2012 4:28 AM Lateral Earth Pressure The magnitude of lateral earth pressure depends on: 1. Shear strength characteristics of soil 2. Lateral strain condition 3. Pore water pressure 4. State

### Enhancing Prediction Accuracy In Sift Theory

18 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS Enhancing Prediction Accuracy In Sift Theory J. Wang 1 *, W. K. Chiu 1 Defence Science and Technology Organisation, Fishermans Bend, Australia, Department

### A. Belejo, R. Bento & C. Bhatt Instituto Superior Técnico, Lisbon, Portugal 1.INTRODUCTION

Comparison of different computer programs to predict the seismic performance of SPEAR the SPEAR building building by means of by means of Pushover Analysis A. Belejo, R. Bento & C. Bhatt Instituto Superior

### Influence of residual stresses in the structural behavior of. tubular columns and arches. Nuno Rocha Cima Gomes

October 2014 Influence of residual stresses in the structural behavior of Abstract tubular columns and arches Nuno Rocha Cima Gomes Instituto Superior Técnico, Universidade de Lisboa, Portugal Contact:

### Plane Strain Test for Metal Sheet Characterization

Plane Strain Test for Metal Sheet Characterization Paulo Flores 1, Felix Bonnet 2 and Anne-Marie Habraken 3 1 DIM, University of Concepción, Edmundo Larenas 270, Concepción, Chile 2 ENS - Cachan, Avenue

### Anisotropic modeling of short fibers reinforced thermoplastics materials with LS-DYNA

Anisotropic modeling of short fibers reinforced thermoplastics materials with LS-DYNA Alexandre Hatt 1 1 Faurecia Automotive Seating, Simplified Limited Liability Company 1 Abstract / Summary Polymer thermoplastics

### Fracture Mechanics of Non-Shear Reinforced R/C Beams

Irina Kerelezova Thomas Hansen M. P. Nielsen Fracture Mechanics of Non-Shear Reinforced R/C Beams DANMARKS TEKNISKE UNIVERSITET Report BYG DTU R-154 27 ISSN 161-2917 ISBN 97887-7877-226-5 Fracture Mechanics

### Lap splice length and details of column longitudinal reinforcement at plastic hinge region

Lap length and details of column longitudinal reinforcement at plastic hinge region Hong-Gun Park 1) and Chul-Goo Kim 2) 1), 2 Department of Architecture and Architectural Engineering, Seoul National University,

Chapter 5 CENTRIC TENSION OR COMPRESSION ( AXIAL LOADING ) 5.1 DEFINITION A construction member is subjected to centric (axial) tension or compression if in any cross section the single distinct stress

### UNIT 1 STRESS STRAIN AND DEFORMATION OF SOLIDS, STATES OF STRESS 1. Define stress. When an external force acts on a body, it undergoes deformation.

UNIT 1 STRESS STRAIN AND DEFORMATION OF SOLIDS, STATES OF STRESS 1. Define stress. When an external force acts on a body, it undergoes deformation. At the same time the body resists deformation. The magnitude

### MECE 3321 MECHANICS OF SOLIDS CHAPTER 3

MECE 3321 MECHANICS OF SOLIDS CHAPTER 3 Samantha Ramirez TENSION AND COMPRESSION TESTS Tension and compression tests are used primarily to determine the relationship between σ avg and ε avg in any material.

### Chapter 3. Load and Stress Analysis

Chapter 3 Load and Stress Analysis 2 Shear Force and Bending Moments in Beams Internal shear force V & bending moment M must ensure equilibrium Fig. 3 2 Sign Conventions for Bending and Shear Fig. 3 3

### Crash and Impact Simulation of Composite Structures by Using CAE Process Chain

Crash and Impact Simulation of Composite Structures by Using CAE Process Chain Madhukar Chatiri 1, Thorsten Schütz 2, Anton Matzenmiller 3, Ulrich Stelzmann 1 1 CADFEM GmbH, Grafing/Munich, Germany, mchatiri@cadfem.de

### Elasto-plastic concrete beam analysis by 1-dimensional Finite Element Method

AALBORG UNIVERSITY MASTER S THESIS Elasto-plastic concrete beam analysis by 1-dimensional Finite Element Method Authors: Niels F. Overgaard Martin B. Andreasen Supervisors: Johan Clausen Lars V. Andersen

### Tensile behaviour of anti-symmetric CFRP composite

Available online at www.sciencedirect.com Procedia Engineering 1 (211) 1865 187 ICM11 Tensile behaviour of anti-symmetric CFRP composite K. J. Wong a,b, *, X. J. Gong a, S. Aivazzadeh a, M. N. Tamin b

### **********************************************************************

Department of Civil and Environmental Engineering School of Mining and Petroleum Engineering 3-33 Markin/CNRL Natural Resources Engineering Facility www.engineering.ualberta.ca/civil Tel: 780.492.4235

### STRAIN ASSESSMENT USFOS

1 STRAIN ASSESSMENT IN USFOS 2 CONTENTS: 1 Introduction...3 2 Revised strain calculation model...3 3 Strain predictions for various characteristic cases...4 3.1 Beam with concentrated load at mid span...

### PLANES OF WEAKNESS IN ROCKS, ROCK FRCTURES AND FRACTURED ROCK. Contents

PLANES OF WEAKNESS IN ROCKS, ROCK FRCTURES AND FRACTURED ROCK Contents 7.1 Introduction 7.2 Studies On Jointed Rock Mass 7.2.1 Joint Intensity 7.2.2 Orientation Of Joints 7.2.3 Joint Roughness/Joint Strength

### Nonlinear analysis in ADINA Structures

Nonlinear analysis in ADINA Structures Theodore Sussman, Ph.D. ADINA R&D, Inc, 2016 1 Topics presented Types of nonlinearities Materially nonlinear only Geometrically nonlinear analysis Deformation-dependent

### Game Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost

Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit

### DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS

DEVELOPMENT OF A CONTINUUM PLASTICITY MODEL FOR THE COMMERCIAL FINITE ELEMENT CODE ABAQUS Mohsen Safaei, Wim De Waele Ghent University, Laboratory Soete, Belgium Abstract The present work relates to the

### COLUMN BASE WEAK AXIS ALIGNED ASYMMETRIC FRICTION CONNECTION CYCLIC PERFORMANCE

8 th International Conference on Behavior of Steel Structures in Seismic Areas Shanghai, China, July 1-3, 2015 COLUMN BASE WEAK AXIS ALIGNED ASYMMETRIC FRICTION CONNECTION CYCLIC PERFORMANCE J. Borzouie*,

### Aim of the study Experimental determination of mechanical parameters Local buckling (wrinkling) Failure maps Optimization of sandwich panels

METNET Workshop October 11-12, 2009, Poznań, Poland Experimental and numerical analysis of sandwich metal panels Zbigniew Pozorski, Monika Chuda-Kowalska, Robert Studziński, Andrzej Garstecki Poznan University

### Geometric nonlinear sensitivity analysis for nonparametric shape optimization with non-zero prescribed displacements

0 th World Congress on Structural and Multidisciplinary Optimization May 9-24, 203, Orlando, Florida, USA Geometric nonlinear sensitivity analysis for nonparametric shape optimization with non-zero prescribed

### DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION

October 1-17,, Beijing, China DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION Mohammad M. Ahmadi 1 and Mahdi Ehsani 1 Assistant Professor, Dept. of Civil Engineering, Geotechnical Group,

### CDPM2: A damage-plasticity approach to modelling the failure of concrete

CDPM2: A damage-plasticity approach to modelling the failure of concrete Peter Grassl 1, Dimitrios Xenos 1, Ulrika Nyström 2, Rasmus Rempling 2, Kent Gylltoft 2 arxiv:1307.6998v1 [cond-mat.mtrl-sci] 26

### PACIFIC EARTHQUAKE ENGINEERING RESEARCH CENTER

PACIFIC EARTHQUAKE ENGINEERING RESEARCH CENTER Shear-Flexure Interaction Modeling for Reinforced Concrete Structural Walls and Columns under Reversed Cyclic Loading Kristijan Kolozvari Department of Civil

### INFLUENCE OF FLANGE STIFFNESS ON DUCTILITY BEHAVIOUR OF PLATE GIRDER

International Journal of Civil Structural 6 Environmental And Infrastructure Engineering Research Vol.1, Issue.1 (2011) 1-15 TJPRC Pvt. Ltd.,. INFLUENCE OF FLANGE STIFFNESS ON DUCTILITY BEHAVIOUR OF PLATE

Bilinear Quadrilateral (Q4): CQUAD4 in GENESIS The Q4 element has four nodes and eight nodal dof. The shape can be any quadrilateral; we ll concentrate on a rectangle now. The displacement field in terms

### Comparison of 2D vs. 3D Modeling Approaches for the Analyses of Concrete Faced Rockfill Dams

Comparison of 2D vs. 3D Modeling Approaches for the Analyses of Concrete Faced Rockfill Dams H. F. Özel & Y. Arici Dept. of Civil Engineering, Structural Mechanics Laboratory,K2-303, Middle East Technical

### The Dynamic Response Analysis of Concrete Gravity Dam under the Earthquake

Copyright 2013 Tech Science Press SL, vol.9, no.1, pp.23-36, 2013 The Dynamic Response Analysis of Concrete Gravity Dam under the Earthquake Yang Lu 1, Li Shi-Min 1, Cao Peng 2 and Shen Xin-Pu 1 Abstract:

### A NEW SIMPLIFIED AND EFFICIENT TECHNIQUE FOR FRACTURE BEHAVIOR ANALYSIS OF CONCRETE STRUCTURES

Fracture Mechanics of Concrete Structures Proceedings FRAMCOS-3 AEDFCATO Publishers, D-79104 Freiburg, Germany A NEW SMPLFED AND EFFCENT TECHNQUE FOR FRACTURE BEHAVOR ANALYSS OF CONCRETE STRUCTURES K.

### Size Effects In the Crushing of Honeycomb Structures

45th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics & Materials Conference 19-22 April 2004, Palm Springs, California AIAA 2004-1640 Size Effects In the Crushing of Honeycomb Structures Erik C.

### DEVELOPMENT OF AUTOMATIC CONTROL OF MULTI-STAGE TRIAXIAL TESTS AT THE UNIVERSITY OF MISKOLC

Geosciences and Engineering, Vol. 2, No. 3 (2013), pp. 37 43. DEVELOPMENT OF AUTOMATIC CONTROL OF MULTI-STAGE TRIAXIAL TESTS AT THE UNIVERSITY OF MISKOLC BALÁZS CSUHANICS ÁKOS DEBRECZENI Institute of Mining

### A FINITE ELEMENT MODEL FOR THE ANALYSIS OF DELAMINATIONS IN FRP SHELLS

TRENDS IN COMPUTATIONAL STRUCTURAL MECHANICS W.A. Wall, K.-U. Bletzinger and K. Schweizerhof (Eds.) c CIMNE, Barcelona, Spain 2001 A FINITE ELEMENT MODEL FOR THE ANALYSIS OF DELAMINATIONS IN FRP SHELLS

### The Frictional Regime

The Frictional Regime Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 1/25/2016 10:08 AM We Discuss The Frictional Regime Processes of Brittle Deformation

### Harmonized European standards for construction in Egypt

Harmonized European standards for construction in Egypt EN 1998 - Design of structures for earthquake resistance Jean-Armand Calgaro Chairman of CEN/TC250 Organised with the support of the Egyptian Organization