This paper is concerned with the influence of mesh quality on the finite element analysis of geotechnical problems. It studies the relationship between mesh quality and accuracy of results produced. Finite element method (FEM) is commonly used for the analyses of geotechnical problems such as the bearing capacity of foundation(Potts and Zdravković, 2001). In this study, vertical behaviour (such as the load-displacement relationship) of a strip footing will be analysed due to the change in mesh quality by means of a FE model.
For a successful finite element analysis (FEA) of a geotechnical problem, a constitutive soil model is necessary. Finite element method is useful for numerical modelling of appropriate soil bodies. (Md Nujid and Raihan Taha Professor, 2014). The soil in the model is subdivided into elements by meshing it. Mesh sizing is very critical and relates closely to the accuracy. It is eminent that the numerical results depend largely on the complexity of the model and number of elements involved in the analysis. With more elements, computational time also increases. (T. More and Bindu, 2008). The software ABAQUS is chosen to investigate the various effects of the change of parameters on the soil, and to optimise a mesh suitable for Finite Element Analysis.
Results obtained from the simulation of the FE model will then be validated by correlating them to either laboratory experiment data or existing solutions from literature review. Any differences observed will be discussed. The results will then be used to find the perfect mesh size for the finite element analysis.
TABLE OF CONTENTS
The Finite Element Method (FEM) is a numerical method for solving problems with complex geometries, loading and material properties. It is commonly used in analysis of geotechnical problems to obtain approximate solutions of boundary value problems. The FEM is also known as the Finite Element Analysis (FEA), where large and complex problems are discretised into smaller, more manageable parts called the Finite Elements (FE). Different parameters are specified and behaviour of these finite elements are simulated and analysed under a variety of conditions.
The mesh size plays a vital role in the accuracy of results and computing time as it determines the number elements involved.; According to the FEA theories FE models with coarse mesh yields less accurate results and has a reasonably small computing time. A coarse mesh means that there are less number of elements present and that these elements will have a large element size. Coarse mesh is more suitable for problems where quick and rough estimation of design is required. On the other hand, finer mesh is used in problems that require high accuracy. Finer mesh means smaller element size and thus requiring a longer computing time due to the large number of elements involved, hence making it more expensive.
The finite element method is used in this study to analyse a vertically loaded strip footing on soil (which may either be drained or undrained). A strip footing is a component of shallow foundation, it is a continuous strip of concrete that supports the weight of the whole wall across the area of soil in consideration. As per conventional foundation design, the Mohr-Coulomb failure criteria is used to define the plasticity of the soil. Accuracy of the results are investigated through modelling laboratory experiments and comparing with any existing data from literature review.
1.2 OBJECTIVE AND SCOPE
The objective of this paper is to study the influence of mesh size on the vertical loading behaviour of a strip footing and provide guidelines to select the optimal mesh size for different finite element analyses. In this study, the bearing capacity analysed using ABAQUS. Boundary conditions will be defied and loads will be applied to the model and corresponding stress, strain and deformations can be obtained from the analysis. The strip footing on the model will be a rigid body so no deformation will be there due to the vertical loading. The analysis will be run several times using FE models of varying mesh sizes in each of the simulations. The penetration of the strip footing into the soil will be observed and recorded. Bearing capacity of the foundation will then be analysed. This is done to show the effects on stress and displacement as the mesh sizes are varied. This helps achieve the optimal combination of accuracy of results and efficiency of the analysis.
A certain amount of error will be present in the solutions produced by the finite element analysis depending on the size, type and accuracy of the model used relying on the availability of appropriate soil data. To minimise this, validation of the model is compared to the results obtained by other measured tests or any existing results obtained from the literature review.
Foundations are an essential part of a geotechnical structure. It allows the loads from the structure to be transferred into the soil. The type of foundation is determined on the basis of the depth of the soil in which it is constructed. There are two main types of foundation, shallow foundation and deep foundation. Deep foundations include basements, caissons, shafts and piles. They are usually built at depths that are greater than 3m while Shallow foundations are used for depths of 1.5 m. they are usually strong enough to support the structural loads imposed on them. The ratio of the embedment depth of these shallow foundations to their width is usually less than 3 (D/B<3). Shallow foundations include spread footings, strip footings, mat foundation and raft foundation. In this paper the bearing capacity of a strip footing will be looked at. A strip footing is a component of shallow foundation. They are used in case of a load bearing wall to evenly spread out the loads applied on the wall to the ground. It is also used when a row of columns is placed close to each other and their spread footings almost touch one another. Therefore, instead of placing each spread footing individually it appears to be more economical to provide a strip footing. Strip footing also goes by the name continuous footing. The strip footing and soil interaction can be simply illustrated in figure 1. A strip footing having width of B and depth of D is placed on top of a soil model which then slowly embeds shallowly.
Figure 1 Schematic of strip footing and soil interaction as it slowly embeds.
Bearing capacity (q) is defined as the maximum bearing stress of the soil to support the loads that are imposed on the structure after which shear failure occurs in the soil. The ultimate bearing capacity (qult) causes settlement due to the shear failure. To ensure that the foundation is stable and that settlement is tolerable, the bearing pressure needs to be within allowable limits. The allowable bearing capacity can be calculated using a Factor of Safety (Fs) with the ultimate bearing capacity.
The expression proposed by Karl Von Terzaghi is generally used to estimate the bearing capacity of shallow foundation (under both drained and undrained loading (Terzaghi, 1943). The equation is as follows:
= c c + 0 q + 0.5
Where, q = the ultimate bearing capacity (qult). c = the cohesion intercept. 0 = the total stress at the level of the base of the foundation. B = the width of the footing. = the bearing capacity factor that represents the contribution of cohesion of soil. = the bearing capacity factor that represents the contribution of surcharge to the overall bearing capacity. = the bearing capacity factor that represents the contribution of soil self-weight.
Some of the Assumptions Terzaghi used in his analysis were that the strip footing was at a shallow depth with a rough base and that the length, L was greater than five times its width, B (L>5B) and that the depth, D was greater than its width, B (D > B). The soil was assumed to be homogenous and comparatively incompressible. The failure zones were expected not to spread above the horizontal plane over the footing’s base. Lastly, the elastic zone was thought to have straight boundaries with an inclination of џ = φ to the horizontal and the plastic zones being developed fully.
The equation was later attempted to be improved using shape factors Sc and Sy. When restrictions were discovered. For strip footings both Sc =1.0 and Sy =1.0. Therefore, a shallow strip foundation having with B and depth D, will have the following equation for the bearing capacity:
= (c cSc) + ( 0 q) + (0.5 Sy)
= (c c (1.0)) + ( 0 q) + (0.5 (1.0))
= (c c) + ( 0 q) + (0.5 )
The equation used for the ultimate bearing capacity of a strip footing under drained loading is shown below:
= ( 0 ) + (0.5 )
Where, q = ultimate bearing capacity (qult). 0 = the total stress at the level of the base of the foundation. B = the width of the footing. = bearing capacity factor that represents the contribution of surcharge to the overall bearing capacity. = bearing capacity factor that represents the contribution of soil self-weight.
The position of water table is essential. The groundwater level influences the effective stresses that are used in the analysis for drained loading of a shallow foundation. Both shear strength and effective stress of the soil is affected by shallow groundwater. They decrease while pore pressure on the other hand increases. If the water table is located below the strip footing’s base, then water will not have any effect on the ultimate bearing capacity. The values of and can be obtained from the following graph as shown in figure 2. Figure 2 plots the values of and against the soil friction angle.
Figure 2 Bearing capacity factors Nq, Ny (after Terzaghi and Peak, 1967)
The equation used for the ultimate bearing capacity of a strip footing on a saturated undrained clay is shown below in equation. For an undrained analysis, the total stresses are used. The cohesion intercept of the undrained clay is equal to its undrained shear strength, Su and the effective stress friction angle, φ’ is zero.
= c c + 0 q
= × + 0
Where, q = ultimate bearing capacity (qult). = the undrained shear strength, c = bearing capacity factor that represents the contribution of cohesion of soil. It depends on the shape and depth of the foundation. 0 = the total stress at the level of the base of the foundation.
For a strip footing on the surface, Prandtl (1921) used the plasticity theory to determine the bearing capacity factor Nc to be 5.1416.
The bearing capacity factor, Nc in Terzaghi’s equation seems to increase with the increase in depth for a cohesive soil. (‘Module 4 : Design of Shallow Foundations’). This is demonstrated by Skempton in the following figure 3. The figure shows that a strip footing at ground level has depth to width ratio, D/B = 0 and therefore the Nc value is (2+ ) = 5.1416, while =1 and =0 (Skempton, 1951).
Figure 3 also shows the transition of bearing capacity from shallow to deep for soil that is uniform throughout and beyond depth D and has a uniform strength Su.
Figure 3 Skempton’s bearing capacity factor Nc for clay soils
Literature review of the ultimate bearing capacity of a shallow foundation under vertical loading is presented with emphasis on strip footings. A scenario is analysed where a non-uniform mesh is used. A strip footing is vertically loaded on an undrained clay The undrained clay is modelled by means of linear elastic Tresca model. This is done to be consistent with conventional design of foundation. Result obtained in the analyses are compared to the solutions of conventional bearing capacity(Potts and Zdravković, 2001).
In order to obtain the load-displacement curves for rigid strip surface footings, analyses are performed. Elasto-plastic soil with Tresca yield surface and the following properties were considered. The soil had a Young’s Modulus, E of 100 MPa, a Poisson’s Ratio of 0.49 and an undrained strength, Su of 100 KPa.
The kind of mesh used is uniform around the strip footing. As it goes further away the mesh gradually becomes non-uniform as shown in the figure 3:
Figure 4 Finite element mesh for strip footing analysis.
Boundary conditions have been applied to the soil, where the two sides of the mesh have been restricted in the horizontal direction. While the bottom layer has been made fixed. It is restricted in both horizontal and vertical direction. Increments of the vertical displacement are applied to the surface of the soil below the position of the footing as loads are simulated. As the strip footing penetrates the soil, displacement is recorded. The load against displacement curve is shown below:
Figure 5 Load-Displacement curves for (smooth, rough with interface and rough) strip footing.
The bearing capacity of a strip footing loaded vertically on undrained clay is usually expressed as:
Qmax = A NC SU
Where,Qmax = Maximum vertical load applied to the footing. A = Area of the footing. NC = Bearing capacity factor. SU = Undrained strength.
The load is expressed in terms of mobilized bearing capacity factor, NC mob in Figure 4.
NC = =
The displacement is distributed by B.
Settlement = =
δ= displacement and B = half width of the footing
Figure 4 shows results from three analyses modelled. One is a smooth footing while the other two are rough. Interface elements were added between the soil and the footing in one out of the two rough ones. This is illustrated in the figure 6:
Figure 6: detail of mesh with interface elements
A shear strength of 100 KPa and shear and normal stiffness of value, Ks = Kn = 105kN/m3 were given to the interface elements. It can be seen on figure 5, that different load-displacement curves and Nc values were produced for the three analyses.
The theory of conventional bearing capacity indicates that the bearing capacity factor, Nc for both smooth and rough footing should be the same as (2+π) = 5.1416 (Prandtl, 1920). The value is theoretically exact and can be obtained from limit analysis and a closed form plasticity solution. By comparing the results obtained from the finite element analysis and the analytical solution, it is seen that the results differ by 0.94%, 4.6%, 2.8% for the smooth (Nc =5.19), rough (Nc =5.39) and rough with interface element (Nc =5.29) respectively.
Figure 7: Vectors of incremental displacements at failure for three different soil-footing interfaces.
Figure 7 shows the vectors of incremental displacement at failure which is the last increment of each analysis. The direction of movement is indicated by the orientation of these vectors and magnitude of the movement by the length. The failure mechanism is indicated by the orientation and the relative magnitude of the vectors. The failure mechanisms indicate some movement in the horizontal direction at the surface of the soil under the smooth footing. However due to boundary condition of the rough footings, they are not allowed to move so the failure mechanisms are a little deeper and wider in comparison. A bigger ultimate footing load and Nc is obtained due to the larger slip surface area, A which is then combined with Su.
The nature of finite element mesh below and adjacent the footing may contribute to the errors. As per figure 7 it is seen that the direction of displacement vectors next to the edge of the rough footing change abruptly from going vertical downwards at an angle of 45° under the footing corner to horizontal in an upward path next to the footing on the soil surface. This happens over the single element that is located at the corner. Ability of this single element to put up with the rapid change of displacement will determine the accuracy of the results. To obtain accurate solutions, a finer mesh needs to be used. This means that the element sizes should be made even smaller around the edge of the footing. Due to this singularity at the edge of the footing, computation time becomes larger. Nonetheless, more accurate results could be obtained with the help of a refined mesh. To increase accuracy of the results further, the overall mesh size throughout the model can be made smaller instead of just the ones near the corner.
CHAPTER 3: FINITE ELEMENT ANALYSIS
3.1 FINITE ELEMENT METHOD
The Finite Element Analysis (FEA) is the implementation of the Finite Element Method (FEM) that is used to solve certain problems. The Finite Element Method is a numerical method used to solve partial differential equations with boundary conditions by using multiple sub functions in a polynomial interpolation.(WELCOME TO ENGINEERS CADD CENTRE (P) LTD). This method draws out the shape of the original differential equation curve with the help of several polynomial curve shapes in order to simplify the more complex curve. This method is helpful for problems with complicated shapes, loadings, material properties etc. (‘Introduction to Finite Element Analysis (FEA) or Finite Element Method (FEM)’). The intervals are then evaluated with respect to the points that encloses that section of the polynomial curve. The FEA solutions are then evaluated and the results (such as deflection and stress) are only given at these points called nodes.
According to the documentation of Abaqus (2010), it is important to discretize the geometry of a structure into finite element before any finite element simulation is run. Therefore, the structure that will be analysed is subdivided in a mesh of finite elements, breaking down the original shape of the structure or the object into simpler shapes. The finite elements are connected to each other by their nodes. The nodes are a coordinate location which define the degrees of freedom (DOFs) in space. The movement of these points under loading of the structure is represented by the degrees of freedom. They also represent the forces and moments that are being transferred to the next element. The finite elements and the nodes collectively form up the mesh.
It is assumed that the simple polynomial curve shapes and the nodal displacements determine the change in displacement of each of these elements. Stress and strain equations are constructed with respect to the unknown nodal displacements. The equations of equilibrium are then assembled in the form of a matrix so that they can be programmed and solved. After suitable boundary conditions are applied, the stiffness matrix equation can then be solved in order to find nodal displacements. The element stresses and strains can easily be calculated once the nodal displacements are obtained.
Abaqus is a general-purpose simulation software used in different engineering fields and is based on finite element method. It has a graphical user interface to pre-process the input files and post process the results. A lot of literature reviews suggest that Abaqus is frequently used to solve several geotechnical problems. An example of the features which make this software popular to use is the availability of the Mohr-Coulomb yield criteria application. The Mohr – Coulomb yield criterion is frequently used for defining the elasticity of soils in finite element method (Hasan Emre Oktay, 2012). Abaqus is chosen as the software to be used for the finite element analysis due to its several advantages. Numerous ways of modelling and multiple choices of algorithms are offered to find solutions of the modelled problem. A number of powerful mesh generating algorithms can be inserted into each of the modelling schemes. These algorithms (functions) will then determine whether the system will behave in a linear or a non-linear manner. Non-linear systems are more complex when compared to the linear systems. Linear systems usually ignore a few restraints of the model such as loading and behaviour. Non-linear systems on the other hand take more realistic behaviour of the model such as plastic deformation, varying loads etc. into consideration. This allows it to test the model up to its failure. Having said that, for any method chosen, it consists of many predefined formulations for elements and material models which can produce acceptable results.
Abaqus supports scripting functionality and command line can be accessed. It allows the user to create a very detailed model. In addition, it can solve problems that behave in both linear and nonlinear fashion and it can conduct 2D as well as 3D analyses. Furthermore, software license is easily available making it more desirable to use. Three simple steps used in the software are shown in figure 8.
An ASCII text file is required as an input for the Abaqus analysis product, Abaqus/Standard. All aspects of the model such as the coordinates of the nodes, analysis steps and requested output variables are held by this file. This file can be either created manually from scratch, or edited using an existing file. The Abaqus CAE pre-processor software can also be used. The Abaqus CAE uses a graphical user interface for them to create, edit monitor, diagnose and visualize advanced analysis. It is usually very user friendly for beginners. Most of the features of Abaqus can be utilized from the Abaqus CAE. However, input files will be used instead of Abaqus CAE in this study. Abaqus CAE can either store modifications or the input files can be modified manually. Certain variable value of selected entities can be saved by Abaqus as an output. The state of a complete model for a certain variable can also be recorded. Before it was called the history output request but now it has been changed to field output request as per Abaqus documentation of 2010. The requests are written to an output database file. The binary format of the file can be read by Abaqus CAE. Forms of other output formats exist, where history output can be saved in either ASCII file formats or in binary that can be opened by third party software.
Evaluation and simulation
Abaqus/CAE or other products
Abaqus/CAE or other products
Figure 8: Simple tests of Abaqus
3.3.1 Brief overview
To design the Finite Element model of the soil for the strip footing, the geometric domain of the problem needs to be defined along with the element type, material and geometrical properties like the length, area etc. The model can then be meshed using variable mesh sizes in each simulation to determine the element size that is most suitable. Lastly, the boundary conditions will be applied. After the model is created, the simulations will be run and analysed. The results obtained will then be compared to existing solutions from literature review and other laboratory experiment data (if available). By doing this the effects of mesh sizing on the bearing capacity of the strip footing can be found.
The type of soil chosen for the analysis was an undrained soil having a density of 0.1 kg/m3. Young’s modulus of 400×109 Pa and a Poisson’s ratio of 0.495 was assigned to the finite element soil model. The dimension of the soil model was 10m long with a height of 5m. The breadth of the model varied with the size of the mesh for each simulation (for example, 0.2m, 0.1m, 0.05m etc.). A 2m long strip footing model was placed at the centre of the surface of the soil model. The strip footing had a rectangular shape with a height of 0.5m and breadth of 0.2m. The surface of the soil interacting with the base of the strip footing was made rough by setting the cohesion factor, α as 1. This means that the friction was equal to the normal force. The three-dimensional soil model was constructed using an input file and then it was imported into the ABAQUS software. The explicit dynamic finite element module was then used for the analysis.
It was seen that the finite element model could be established using both the ABAQUS CAE and by writing input files based on previous sample work. Based on the conclusion that both files produced same results, input files were chosen to create the FE model to make the process time efficient.
Three finite element models were generated for the soil with mesh size of 0.20m, 0.10m and 0.05m using input files. A sample input file has been attached as appendix. In the input files, the model was divided into two different bodies, a sea-bed and a rigid strip footing. The material for the sea-bed is chosen to be a eulerian material. The coordinates of the sea-bed are defined and the node numbers are assigned. The whole domain is then divided and grouped into node sets (Nset). The nodes on the top and bottom are grouped and named ‘Bottomtop’. The nodes on the left side of the model are grouped and ‘Left’. The nodes on the right side of the model are grouped named ‘Right’. A few more node sets were defined. Similarly, the sea-bed, the void on top of the surface and the pipe (the strip footing) were generated by allocating the node numbers and were grouped into element sets (Elset).
The coordinates of the strip footing were then defined and node numbers were assigned just like the sea-bed. It can be seen on the input files that the node numbers of the strip footing are listed under the heading ‘node’ along with their coordinates. The element numbers for the strip footing is also shown under the heading ‘element’ along with the number of nodes that make up the element. The type of element chosen for the strip footing is the R3D4. This shows that the element is a 3D rigid quadrilateral element. The reference point for the strip footing was set by specifying the node number. After creating the FE model, the material of the sea-bed needs to be assigned. As mentioned earlier, an undrained soil is selected. Apart from the density, young’s modulus, Poisson’s ratio, the cohesion yield stress of the soil is assigned as 1. After assigning all of them in element sets, the boundary conditions are applied. It is usually easier to change parameters or apply boundary conditions to a collection of nodes or elements rather than applying them to each individual node or element. This process can become very difficult and time consuming. In these situations, it is helpful to break down the FE model into smaller sets, making the process much simpler and time saving.
The strip footing did not experience any deformation while it penetrated since it was modelled as a rigid body. The soil was modelled to have a homogenous shear strength all over the model. The edge on the right side was restrained in both vertical and horizontal direction thus making it a fixed end. Similarly, the bottom of the soil domain was also made fixed. Symmetry boundary conditions were assigned to the left edge. This is because It was assumed that the model was symmetrical about the plane normal to the center of the strip footing. Therefore, only the right half of the arrangement was modelled and analysed. The plane of symmetry is shown by the symmetry line in the x-z plane in Figure 9. The strip footing moved in the vertical direction in the downward y plane as it settled.
Figure 9 Half of the arrangement to be analysed
The results of the analysis will be shown and discussed in detail in this chapter. The main objective is to describe the influence of mesh quality on the bearing capacity of the strip footing based on the results obtained from the finite element analysis. Results obtained from the Finite Element Analysis are put into graphical solution so that they can be compared directly. These results are compared among each other and then compared with the analytical value obtained from literature review to see which mesh has the most suitable results.
The outcome of the three input files are shown in this chapter. The figures 9, 10 and 11 show the finite element model with the mesh sizes 0.02m, 0.10m and 0.05m respectively. The strip footing is represented by the red box. The sea-bed is represented by the blue box and the void above the sea-bed is respresented by the purple box.
Figure 10 Finite element model with mesh size 0.20m
Figure 11 Finite element model with mesh size 0.10m
Figure 12 Finite element model with mesh size 0.05m
Analysis part 1
As the strip footing is embedded into the soil bed, the soil-bed undergoes displacement. The displacement increases linearly as the time increases. For all three meshes, the displacement was seen to be 0.1m vertically down for half of the model that was analysed. Figure 12 shows the shape of the displacement against time graph that was obtained for each of the meshes.