ICE Paper - Nonlinear Finite Element Analysis of Concrete
ICE Paper - Nonlinear Finite Element Analysis of Concrete
http://dx.doi.org/10.1680/jencm.17.00002
Non-linear analysis of concrete structures: Paper 1700002
Received 17/02/2017 Accepted 13/05/2017
a practical design approach to Eurocode
Keywords: codes of practice & standards/offshore engineering/
Walker, Abhyankar and Gration stress analysis
Reinforced concrete is an inherently non-linear material with non-linear compressive behaviour, concrete cracking and
reinforcement yielding altering the stiffness and causing redistribution of forces within statically indeterminate structures.
A practical design approach is described that has been developed to ensure robust computation of serviceability limit
state and ultimate limit state demand in line with Eurocode BS EN 1992, including explicit modelling of material and
geometric non-linearities and construction imperfections. The approach uses the MAT_CONCRETE_EC2 material model
within the LS-DYNA finite-element analysis package and has been extensively used during the detailed engineering
design of the concrete gravity substructure for the White Rose Extension Project off the coast of Newfoundland, Canada.
The material model and validation of its behaviour in accordance with the Eurocode standard are described, along with
other typical engineering applications where the approach may be of benefit. Specific benefits and implementation
challenges encountered during the design process are also highlighted with a non-linear buckling assessment case study.
1
Engineering and Computational Mechanics Non-linear analysis of concrete structures:
a practical design approach to Eurocode
Walker, Abhyankar and Gration
2
Engineering and Computational Mechanics Non-linear analysis of concrete structures:
a practical design approach to Eurocode
Walker, Abhyankar and Gration
~150 m
~120 m
Figure 1. Indicative render of the White Rose Extension Project CGS and topside
Shared inputs
Linear FE model
• Geometry
• Loading
Sectional forces
and moments
Utilisation calculation
Additional inputs
• Non-linear material
properties
Non-linear FE model
• Reinforcement
layout
• Geometric
imperfections
safety factor’ method for non-linear analysis of concrete The non-linear LS-DYNA model was composed of multilayer
described in the fib Model Code 2010 (Walraven, 2013), in (also referred to as ‘composite’) shell elements. Each layer con-
which further guidance can be found. tained a single integration point defined as either concrete or
3
Engineering and Computational Mechanics Non-linear analysis of concrete structures:
a practical design approach to Eurocode
Walker, Abhyankar and Gration
Concrete
X-direction
reinforcement
Y-direction
reinforcement
Example two-way spanning slab layout Corresponding composite shell element
4
Engineering and Computational Mechanics Non-linear analysis of concrete structures:
a practical design approach to Eurocode
Walker, Abhyankar and Gration
–25 000 The CGS comprises three main sections, as shown in Figure 7.
Axial force: kN/m
The central shaft supports the topsides and contains the main
–20 000 mechanical outfitting. The central shaft and the closed caisson
cells are initially empty and are gradually flooded to lower the
–15 000
CGS onto the seabed at its final location. The open perimeter
cells provide stability during float-out and are filled with solid
ballast once the CGS is in place, to resist overturning moments
–10 000
due to wave loads. The caisson and perimeter cells are divided
into compartments, separated by radial walls. The buckling
–5000
assessment described here was undertaken for the caisson
–2000 –1500 –1000 –500 0 500 1000 1500 2000 radial walls.
0
Moment: kN/m
Linear analysis identified the governing loading condition for
5000 the caisson radial walls to be just prior to touchdown of the
Concrete section solver MAT_CONCRETE_EC2
CGS onto the seabed. At this point in the installation
sequence, the CGS has not yet been fully flooded, and so there
Figure 6. Comparison of force–moment interaction diagram exists a large differential hydrostatic pressure on the external
produced using a MAT_CONCRETE_EC2 shell element and walls of the caisson (Figure 8). The load case also contains
a non-linear concrete section solver
differential water levels between compartments as a result of
potential variation in the rate of flooding (Figure 9). The
differential compartment water levels lead to out-of-plane
deflections, which combined with the in-plane compression
4. Design example – non-linear buckling generate additional moments and necessitate a buckling assess-
assessment ment of the caisson radial walls.
The benefits and challenges associated with this design
approach are illustrated through a non-linear buckling assess- The design of the caisson radial walls was initially based on a
ment undertaken on the White Rose Extension Project. linear buckling analysis. This method relied on underestimating
Caisson
radial wall
Perimeter
radial wall
5
Engineering and Computational Mechanics Non-linear analysis of concrete structures:
a practical design approach to Eurocode
Walker, Abhyankar and Gration
ULS and SLS material parameters and load factors were speci-
fied for separate strength and crack width assessment. The
results quoted here concern the ULS assessment. The load
Shaft
Caisson case previously described, including differential water levels,
water level
water level was applied to both the linear and non-linear models. The
~20 m
~17 m
explicit LS-DYNA analysis was run statically; the loading was
applied as a linear ramp and held at the required load factor.
~17 m During the touchdown load case, the perimeter radial walls
were primarily loaded in tension and thus when non-linear
~20 m material properties were considered, the perimeter radial wall
elements cracked and softened. As the tensile stiffness of the
~20 m perimeter radial walls decreased, the in-plane compression in
the caisson radial wall increased to take the hydrostatic load.
This effect was investigated in the linear analysis through alter-
ing the stiffness of the perimeter radial walls to represent
cracking. The non-linear analysis provided a more accurate
assessment of this effect, shown in Figures 10–12.
6
Engineering and Computational Mechanics Non-linear analysis of concrete structures:
a practical design approach to Eurocode
Walker, Abhyankar and Gration
Figure 10. Elevation view representation of the in-plane behaviour of the caisson and perimeter radial walls during touchdown
MAX_PRINC_STRESS
mid-surface:
MPa
–20·00
–15·00
–10·00
–5·00
5·00
Figure 11. Mid-surface maximum principal stress plot from linear (left) and non-linear (right) analysis
levels between caisson compartments. Both the linear and non- capture material non-linearity explicitly. The non-linear analy-
linear models were subject to the same loading. Force and sis therefore included the effects accounted for in the linear
moment resultants were compared between the models to vali- buckling assessment and could be used to validate its assump-
date the assumptions used in the linear method. tions. It should be noted that as non-linear material properties
were used throughout the CGS in the non-linear analysis, it
The linear buckling model was set up with initial imperfections explicitly captured the stiffness reduction in adjacent structural
and captured the geometric non-linearity (i.e. P-delta effects). elements, which were conservatively assumed to be fully
The stiffnesses of the caisson radial wall and adjacent elements cracked in the linear buckling analysis.
were also modified to account for material non-linearity,
in accordance with the recommendations in Eurocode 2, An illustrative set of sectional resultants are presented below to
Section 5·8. highlight the benefit of the non-linear analysis. The result of
non-linearity (material and geometric) can be seen by compar-
The non-linear model included geometric non-linearity, initial ing Figures 14 and 15, which shows a significant increase in
imperfections and used the MAT_CONCRETE_EC2 to moment resultants, confirming the sensitivity of this structural
7
Engineering and Computational Mechanics Non-linear analysis of concrete structures:
a practical design approach to Eurocode
Walker, Abhyankar and Gration
MIN_PRINC_STRESS
mid-surface:
MPa
–20·00
–15·00
–10·00
–5·00
0
5·00
Figure 12. Mid-surface minimum principal stress plot from linear (left) and non-linear (right) analysis
Uncracked
concrete Horizontal
Cracked moment resultant:
concrete kNm/m
2600
1940
1280
620
–40
–700
–1360
–2200
–2680
–3340
–4000
Figure 13. Concrete cracking pattern
8
Engineering and Computational Mechanics Non-linear analysis of concrete structures:
a practical design approach to Eurocode
Walker, Abhyankar and Gration
9
Engineering and Computational Mechanics Non-linear analysis of concrete structures:
a practical design approach to Eurocode
Walker, Abhyankar and Gration
Figure 17. Thermal load case for liquefied natural gas tank
10
Engineering and Computational Mechanics Non-linear analysis of concrete structures:
a practical design approach to Eurocode
Walker, Abhyankar and Gration
of more standard load cases, in accordance with Eurocode 2. Gibson R (2012) Performance based seismic design of an offshore
The approach uses the explicit LS-DYNA solver and has been platform in the Caspian Sea. The Structural Engineer 90(11):
14–22.
extensively used during the detailed design of the White Rose
Gogus ABE (2016) Seattle civic square tower-engineering a high rise
CGS. The approach enables a less iterative and more realistic for increasing performance goals and seismic hazard in the
consideration of material and geometric non-linearity than Pacific Northwest. Los Angeles Tall Buildings Structural Design
many existing techniques, all within one structural model. On Council Proceedings of the 2016 Annual Meeting, Los Angeles,
the White Rose Extension Project, it has been used to validate CA, USA.
Hoorelbeke P, Izatt C, Bakke JR, Renoult J and Brewerton RW (2006)
linear structural assessments and provided insight into the
Vapor cloud explosion analysis of onshore petrochemical
structural behaviour to inform a better design. The approach facilities. Proceedings of the American Society of Safety
complements simpler techniques based on linear analysis and Engineers Middle East Chapter, 7th Professional Development
offers excellent opportunities to be used more widely as a prac- Conference & Exhibition, Kingdom of Bahrain, paper
tical design tool aligned with Eurocodes. ASSE-MEC-0306-38.
ISO (International Organization for Standardization) (2013) 19903-07:
Petroleum and natural gas industries – fixed concrete offshore
Acknowledgements structures. ISO, Geneva, Switzerland.
The authors acknowledge Nicholas Christie, Maya Shinozaki, Izatt C (2009) Assessing perforation limits of steel section impacts
Mary Ferguson, Babak Pourazarm, Chris Ariyaratana and on reinforced concrete slabs. Proceedings of the Institution of
Civil Engineers – Structures and Buildings 162(1): 37–44,
Alasdair Parkes who have contributed to the development of the
https://doi.org/10.1680/stbu.2009.162.1.37.
techniques and results presented in this paper. They also give Law AK (2015) Engineering geometrically bi-linear columns to deliver
credit to Arnstein Godejord and Richard Sturt who have con- fire resistance: standard heating. Engineering Structures 100:
tributed greatly to the development of this technique. 590–598, https://doi.org/10.1016/j.engstruct.2015.06.046.
The authors are also grateful to Ian Feltham, Gopal Srinivasan LSTC (Livermore Software Technology Corp.) (2013) LS-DYNA Keyword
User’s Manual, Version 971 R7.0. LSTC, CA, USA.
and Mohammed Ezzat, for their assistance with the examples
NEHRP JV (National Earthquake Hazards Reduction Program Consultants
section. All of those previously mentioned are from Arup. Joint Venture) (2014) NIST GCR 14-917-25: Recommendations for
seismic design of reinforced concrete wall buildings based on
REFERENCES studies of the 2010 Maule, Chile Earthquake. NEHRP JV,
Bentz E (2000) Sectional Analysis of Reinforced Concrete Members. Gaithersburg, MD, USA.
PhD thesis, Department of Civil Engineering, University of Nyhus BS (2014) Consistent practical design of concrete structures.
Toronto, Toronto, Canada. Structural Concrete 15(3): 305–316.
BSI (2004) BS EN 1992-1-1:2004, B: Eurocode 2: Design of concrete OL (Oasys Ltd) (2014) Oasys ADSEC User Manual, Version 8.3. OL,
structures part 1-1 general rules and rules for buildings. BSI, Newcastle, UK.
London, UK. Rombach G (2007) Bow-shaped segments in precast segmental bridges.
Cervenka V (2013) Global safety formats in fib Model Code 2010 for Engineering Structures 30(6): 1711–1719.
design of concrete structures. Proceedings of the 11th International Walraven J (2013) fib Model Code for Concrete Structures 2010.
Probabilistic Workshop, Brno, Czech Republic, pp. 31–40. Ernst & Sohn, Berlin, Germany.
11