Finite Element Analysis of a New Type of Self-Insulating Concrete Masonry Wall System ()
1. Introduction
Design and selection of building materials and their components are an efficient way to reduce energy consumption. Thus, enveloping buildings with thermal insulation materials is one of the most effective ways to reduce energy consumption [1] [2]. Several studies have been carried out on this subject, and it is thought that the building type, shape, construction materials, insulation materials, and costs affect the thickness of an insulator [3]-[7].
In general, external insulation methods are preferred around the world. However, the use of insulation blocks composed of Expanded Polystyrene Foam (EPS), which is an alternative to building external insulation, has become more widespread with the intention of thermal insulation. EPS used for heat and sound insulation or packaging is composed of small, white and interconnected beads and exhibits superior engineering properties due to its structure such as lightweight, versatile, energy efficient, and cost effective. Therefore, it is used as insulation material in buildings and can be molded into many shapes for different purposes [8] [9].
Most previous researches focus on light weight concrete masonry units as good thermal insulation materials in buildings because they have a lower thermal conductivity than normal weight concrete. Unfortunately, masonry units made from light weight concrete have low mechanical properties compared with normal weight concrete [10]-[14]. By inserting insulation material such as EPS into normal weight concrete with a special configuration of concrete masonry units (CMUs) leads to increase in their thermal resistance without an effect of their mechanical performance [15].
2. Modelling Verification
Firstly, the FE using Abaqus software was verified with experimental results which published by Abu-Bakre et al. (2016, 2017) [16] [17] and the details of the tested specimens are shown in Table 1. The comparison of FE model results with the experimental results was aimed to ensure that the material and elements properties, and convergence criteria are suitable to model the behavior of the wall and that the simulation process is correct. Therefore, the walls which tested experimentally were simulated for verification study.
Table 1. Self-insulated concrete masonry shear wall specimens’ details [16] [17].
Wall |
Height (mm) |
Length (mm) |
Thickness (mm) |
Aspect ratio |
Reinforcement |
Axial stress (N/mm2) |
vertical |
horizontal |
ρv% |
ρh% |
SW1 |
1590 |
1590 |
240 |
1.0 |
Ø16@200 mm |
Ø8@200 mm |
0.46 |
0.115 |
0.00 |
SW2 |
3190 |
1590 |
240 |
2.0 |
Ø16@200 mm |
Ø8@200 mm |
0.46 |
0.115 |
0.50 |
SW3 |
3190 |
1590 |
240 |
2.0 |
Ø16@200 mm |
Ø8@200 mm |
0.46 |
0.115 |
1.0 |
SW4 |
5780 |
1590 |
240 |
3.6 |
Ø16@200 mm |
Ø8@200 mm |
0.46 |
0.115 |
1.0 |
SW5 |
5780 |
1590 |
240 |
3.6 |
Ø20@200 mm |
Ø8@200 mm |
0.718 |
0.115 |
1.0 |
SW6 |
5780 |
2190 |
240 |
2.6 |
Ø16@200 mm |
Ø8@200 mm |
0.667 |
0.115 |
1.0 |
As shown in Table 2, there was good relationship between the FE and experimental load-deflection curves. This confirmed the validity of the settled FE models and reliability of the FE analysis.
3. Finite Element Analysis
The finite element method (FEM) is one of the powerful tools for modeling a structure with a very complicated geometry and materials. There are many strategies as shown in Figure 1 to model masonry structure with FEM, which includes macro and micro modeling. The macro model is based on the assumption of homogenous materials, and the mortar joints and units can be smeared into one isotropic or anisotropic material. This procedure may be preferred for the analysis of large masonry structures due to the reduced time and memory requirements as well as a user-friendly mesh generation. In addition, this type of modeling is most valuable when a compromise between accuracy and efficiency is needed [18] [19].
Finite element analysis was performed using Abaqus software to validate the experimental results of SCMSW specimens.
Figure 1. Modeling strategies for block masonry. (a) Typical masonry sample; (b) Detailed micro modeling; (c) Simplified micro modeling; and (d) Macro modeling.
3.1. Material Constitutive Models
The damage plasticity model from the Abaqus software was applied to simulate the concrete constitutive law. Using the same methods described in chapter three. The masonry compressive strength obtained was extracted from compression test data for masonry prism, and Poisson’s ratio was taken as 0.2. The concrete material used in the foundation and loading beams was obtained from concrete compression test, and Poisson’s ratio was taken as 0.2.
The “plasticity” model from Abaqus software was adopted to simulate the reinforcement law. The stress-strain curve is plotted in Figure 2. Each stress-strain curve, made up of two linear portions, represents the character of a bare mild steel bar; where the modulus Es is the elastic modulus of the reinforcement and the modulus E' is the deformation modulus at the strain hardening stage, E' = 0.01Es. The test yield strength and elastic modulus of the reinforcement were adopted in chapter five, and Poisson’s ratio was taken as 0.3.
Figure 2. Stress-strain curve of reinforcement.
3.2. Finite Element Analysis Model
Solid element C3D8R was used for concrete and masonry wall, and truss element T3D2 was used for reinforcement. The interface model between masonry wall and loading beam was composed by contact in the lateral axis and stick-slip along the tangential axis. The “hard contact” model was adopted for the contact on the lateral axis and the “Coulomb friction” model for the stick-slip along the tangential axis. The friction coefficient was taken as 0.7 according to the code GB5003-2011 [20]. Reinforcements were all embedded in the concrete. The bottom of the base was fixed rigidly. A vertically distributed load was applied to the top of the loading beam first and kept constant. Then, horizontal displacement was applied to a coupling node in the middle plane of the loading beam. A typical finite element analysis model of SW2 is shown in Figure 3.
Figure 3. Finite element analysis model of SW2: (a) Loading and boundary conditions (BCs); and (b) Element mesh.
4. Finite Element Simulation Results and Analysis
Figure 4 shows the Von Mises stresses distribution in SW3 specimen at maximum load. While Table 2 shows the maximum load values measured by experiment and predicted by finite element software Abaqus. A good agreement is shown between the numerical and experimental results.
Figure 4. SW3 F.E. results: (a) Von Mises stresses distribution at wall reinforcement; and (b) Von Mises stresses at masonry wall.
Table 2. Comparison of predicted load by F.E and experimental capacities.
Wall |
Aspect ratio |
Axial stress (N/mm2) |
Vertical
reinforcement |
F.E. maximum load (kN) |
Average observed capacity (kN) |
Ratio of observed
capacity to F.E. results |
SW1 |
1.0 |
0.0 |
Ø16@200 mm |
290 |
327 |
1.13 |
SW2 |
2.0 |
0.5 |
Ø16@200 mm |
194.4 |
159 |
0.82 |
SW3 |
2.0 |
1.0 |
Ø16@200 mm |
208.6 |
209.7 |
0.95 |
SW4 |
3.6 |
1.0 |
Ø16@200 mm |
125.4 |
112.6 |
1.0 |
SW5 |
3.6 |
1.0 |
Ø20@200 mm |
154.5 |
146.3 |
0.95 |
SW6 |
2.6 |
1.0 |
Ø16@200 mm |
320 |
359.6 |
1.12 |
5. Design Example
The 18-storey residential building as shown in Figure 5 is composed of self-insulated concrete masonry shear walls was analyzed by using ETABS (v9.5.0) software. The building is designed for the seismic intensity of 8 according to Chinese standard.
(a)
(b)
Figure 5. Model example geometry (a) Typical architecture floor plan (b) reinforcement details.
The structural model was generated using ETABS software to give a practical application for the self-insulating concrete masonry shear walls in general engineering practice. Self-weight of structural elements, finishing load of 3 kN/m2 and live load of 2 kN/m2, act at each floor and at roof level were applied to the structural model. The maximum horizontal displacements for X and Y direction under earthquake excitations were found 15 and 30 mm, respectively. Moreover, the first mode shape time period is 0.97 second, while the second and third mode shapes time period are 0.14 and 0.07 second, respectively. According to this example the new type of self-insulating concrete masonry shear walls can be used in areas have seismic intensity of 8 or less and it is application valid until 18-storey height building. The FE model and shape of the studied Building are shown in Figure 6.
(a) (b)
Figure 6. (a) Rendered view of the model (b) Building mode shapes.
6. Conclusion
This paper evaluated the effects that the height-to-length aspect ratio, axial compressive stress, and reinforcement ratio have on the behavior of self-insulating concrete masonry shear walls under in-plane cyclic loading. SCMSW specimens’ performance was evaluated based on a comparison of predicted vs. actual load capacity, drift capacity, displacement ductility, height of plasticity, equivalent plastic hinge length, amount of energy dissipated, and value of equivalent hysteretic damping. Moreover, finite element results obtained from the Abaqus program were compared to those obtained experimentally. Finally, a design example of the 18-storey building under an earthquake excitation was provided to give a practical application of the self-insulating concrete masonry shear walls.