Analysis of reinforced concrete with FEM

Background
The FEM work bench has the capability of estimating the level of reinforcement required in a concrete structure to prevent brittle failure under tension or shear.



This is done with the method described here:. In essence it is a post processing routine for Calculix, which calculates the principal tensile stresses in the concrete from an elastic analysis and uses those to determine the minimum reinforcement in the three coordinate directions required to prevent failure. In the analysis it is assumed that the concrete material cannot carry tensile stresses, whereas the steel is utilized to its maximum capacity (i.e. reaches yield).

The required reinforcement is expressed in terms of a reinforcement ratio. This is the ratio of steel to concrete area. For example a reinforcement ratio of 0.01 in x-direction (rx=0.01) means that the total cross sectional area of reinforcement bars running in x direction should be 1% of the concrete cross sectional area they are passing through. A hypothetical cross section of 1mx1m should in that case contain 0.01 m2 steel, which could be achieved by using 90 reinforcement bars of 12mm diameter each (steel area = 90*PI*(0.012)^2/4=0.0102 m^2). If the required reinforcement ratio over this concrete cross section is uniform then the bars could be placed at an equidistance 9x10 grid with a center-to-center distance of approximately 10cm. This is still a practical number where sufficient space between bars is left for the concrete to pass through and ensure a high quality cover. Much higher values would lead to a very dense reinforcement grid with potential quality issues, whereas much lower values could lead to large tension cracks in the cross section between the bars. A typical range in practice is from 0.002 to 0.02 (= 0.2% to 2%). Further guidance can be found in design codes.

If the required reinforcement ratio is not uniform over the full cross section then the cross section can be divided into pragmatic sub-sections with more or less uniform ratio and reinforcement applied to those cross-sections. An example will be given later on.

As a word of caution, it takes much more to design a safe and durable concrete structure than what the FEM work bench can currently provide. For example, the method does not calculate crack width (important for durability and functionality), accurate deformations (FEM results for concrete are simply linear-elastic) or take account of reinforcement anchoring requirements (causing an increase of required reinforcement ratios in anchoring zones). It also doesn’t predict concrete crushing (although an indication of that can be obtained by plotting the Mohr Coulomb stress - see further), which could mean that the concrete fails before the reinforcement yields, causing brittle failure of the overall structure. These and other limitations mean that the FEM concrete functionality can only be used to assess conceptual designs, whereas detailed design decisions critical to safety and performance should be left to qualified professionals.

Model geometry, loads and supports
Although the FEM concrete routine does not have any additional requirements for geometry, loads and supports, it should be borne in mind that sharp corners and a support on an edge or vertex can introduce stress concentrations that will lead to extremely high and unrealistic reinforcement ratios at or near those locations.

Material Parameters
FEM workbench has a special material object for reinforced materials, which combines a matrix material (e.g. concrete) and a reinforcement material (e.g.steel). For the analysis of reinforced concrete with FEM, the following parameters need to be specified, as a minimum:

for concrete:

- Young’s modulus (used in the Calculix analysis to calculate elastic deformations and stresses) - Poisson ratio (same) - uniaxial compressive strength (used during post-processing in FEM to calculate the Mohr Coulomb stress as an indicator for crushing or shear failure in concrete) - friction angle (same)

for steel:

- yield strength (used during post-processing in FEM to calculate reinforcement ratios)

Please note that three types of analysis are performed: 1) An elastic analysis using Calculix (only utilising the elastic parameters for concrete); 2) A post processing step to analyse the required reinforcement (only utilising the yield strength of steel) and 3) Calculation of the Mohr Coulomb stress (only using the strength parameters of concrete, i.e. uniaxial compressive strength and friction angle). The Mohr Coulomb stress can be reviewed in the VTK pipeline.

Application
In the remainder of this article a few practical cases will be analysed to discuss application of the method.

Simply supported beam with uniform load
A 4.0x0.1x0.3m concrete beam is loaded by self weight and a 100kN (25kN/m) distributed load.

The material parameters are as follows:

for concrete:

- Young’s modulus = 32 GPa (as per Calculix default for concrete) - Poisson ratio = 0.17 (as per Calculix default for concrete) - uniaxial compressive strength = 30 MPa (concrete type C30/37) - friction angle = 30 degrees

for steel:

- yield strength = 500 MPa

The specific weight of the concrete is taken as 24kN/m^3

The required reinforcement in x-direction is very high (5.4%) and exceeds typical maximum percentages allowed by code to prevent brittle failure. The high shear stresses at the supports also lead to a requirement of high reinforcement:



The Mohr Coulomb plot shows that beam is indeed prone to crushing on the compression side (Mohr Coulomb stress > 0.0), as would be expected with a very high reinforcement percentage:



Both the reinforcement ratio and Mohr Coulomb stress indicate that we have an issue and that we need to rethink our conceptual design. Potential solutions are to increase the beam dimensions or use pre-stressed concrete. Further details can be found in the following post:

https://forum.freecadweb.org/viewtopic.php?f=18&t=28821&start=10#p235003

Beam with mid-span support
A 8.0x0.2x0.4m concrete beam is loaded by self weight and a 160kN (20kN/m) distributed load.

The material parameters are as follows:

for concrete:

- Young’s modulus = 32 GPa (as per Calculix default for concrete) - Poisson ratio = 0.17 (as per Calculix default for concrete) - uniaxial compressive strength = 25 MPa (concrete type B25) - friction angle = 30 degrees

for steel:

- yield strength = 286 MPa (reduced from 500 MPa to account for a safety factor of 1.75)

The specific weight of the concrete is taken as 24 kN/m^3

The ParaView plot of the exported VTK file shows that the reinforcement requirement is largest at the top of the beam near the central support. Here the highest bending moment occurs. The maximum reinforcement ratio is with 0.02 at the high end of the practical range quoted earlier:



The required area of steel at the central support can be obtained with a ParaView integration filter applied to the mid section of the beam:



The panel at the bottom of this picture shows that the total required steel area at this cross section is 389.6 mm^2. As one reinforcement bar of diameter 12mm has a cross sectional area of 113mm^2, it means that 4 bars would be required, giving a cross sectional area of 452 mm^2. These would need to be placed near the top of the beam, while maintaining sufficient concrete cover. In general, the center of gravity of the reinforcement ratio plot gives a good indication of where the reinforcement should be centered:



The reinforcement requirement found above agrees well with that obtained using traditional methods:

https://forum.freecadweb.org/viewtopic.php?f=18&t=28821&start=20#p235063

Finally, a Mohr Coulomb stress check should be performed to check potential crushing of the concrete. For this check, the characteristic compressive strength of concrete (25MPa) should be divided by an appropriate material factor (>1.0).

Shear wall with uniform load
Summary, picture and link (under construction).

https://forum.freecadweb.org/viewtopic.php?f=18&t=28821&start=10#p234673



Deep beam with opening
Summary, picture and link (under construction).

https://forum.freecadweb.org/viewtopic.php?f=18&t=33049