Because of the safety concerns, we have to prove by means of safety assessment and calculations and build confidence in the safety of the CC and show that the consequences of an aircraft impact on the environment and humans are negligible. These consequences are traditionally categorized into local and global effects. Accordingly, we must distinguish between global and local damage (or local and global response). The global effects include overall axial, bending, and shear effects in the structural elements between the impact area and support locations.
The global damage identifies the collapse of large portions of the building walls, load carrying members, etc., apart from vibration, that can potentially shake the whole building structure. The extent of such a “global failure” will depend primarily on the dynamic characteristics of the CC and missile, and the relative damage is, in general, associated with the excessive deformation of the entire structural system, assuming that local perforation does not occur. The behavior of concrete in compression and its compressive strength become very important in assessing the local behavior, whereas the tensile strength properties and the stress–strain response of reinforcement steel, in many cases, dominate the deflection and vibration behavior of the reinforced concrete structural members.
As for the local effects concerned, according to [
19] these are due to hard (or semi-hard) missile impacts, resulting in punching failure mode. The sequence of localized loading effects consists of (
Figure 6): (a) penetration into the target; (b) cone cracking and plugging due to the inertial stress wave; (c) spalling with the ejection of target material from the proximal face of the target; (d) radial cracking behind initial wave front; (e) scabbing with the ejection of fragments from the distal face of the target; (f) potentially, perforation completely through the target; and failure. Although these would not, in general, determine a structural collapse, they have to be considered because of the secondary effects on the safety-related systems or components.
For the purposes of this study, a deterministic method which uses dynamic simulations to verify “the acceptance criteria” are respected, as required by the applicants, has been adopted.
The numerical simulation of a reinforced concrete structure is a challenging task due to the inhomogeneity of the material (non-linear behavior) which is also dependent on the strain rate. In order to correctly consider the dynamic characteristics of an impacted structure, a missile–target interaction method with explicit integration techniques has been adopted.
3.1. CC Model
The reinforced concrete containment of this study is that of a typical pressurized water reactor (60 m tall, about 45 m diameter and 1 m thickness) resting on a rock soil (“rigid foundation”), represented by means of clamped restraints. To comply with the ASCE rules and ACI standards, the CC reinforcement is made of #8 reinforcing steel bars in the dome, distributed about 30 mm from the center of each direction and each face, #18 vertical bars in the cylindrical part at about 30 mm spacing for each face, and #18 horizontal bars on both sides of the vertical reinforcing grid. The basement has instead steel plates of 13 mm thickness [
21]. The material properties of reinforcement are provided in
Table 1.
Figure 7 shows the FE model of CC that was made of more than 68,000 solid elements. The concrete structure and the main body of F4 plane are modelled with SOLID-3D Lagrange elements and the internal structures using 3D thick shell elements. The TRUSS-3D elements steel is used for the reinforcement (rebar).
To achieve sufficient accuracy of the dynamic analysis results, sensitivity studies were performed to identify the allowable maximum element size. In determining this, the frequency of the highest mode and the longitudinal propagation velocity of the concrete material were considered. The mesh size was chosen to be uniform in each 15° angular sector of CC discretization, with 20 elements through the thickness (element vertical dimension of 1 m) and about 72 elements across the airplane cross-section.
The connection of the reinforcement and concrete may be a crucial issue. All the models presented in this paper utilize a perfect bond with fully kinematic coupling, therefore the effects of bond slipping cannot be taken into account.
Adequate initial and boundary conditions (continuum mechanics) have been assumed so as to preserve energy and momentum.
The mathematical model is able to describe the characteristic of contact deformation during crushing. The contact force is automatically calculated by a code based on the conservation of the impulse. Moreover, assuming a continuous material, contact theory and structure mechanics allowed us to describe the CC mechanical behavior through its mechanical properties such as Young’s modulus, yielding stress value, and contact type, etc. A proper contact table describing the interaction between the two bodies—i.e., missile and target—was also defined. The contact detection method was described by a stress-based contact algorithm (which allows us to decide whether or not bodies in contact should separate).
In simulations dealing with impact loading, pronounced material nonlinearities occur. The constitutive laws, especially for concrete, must be able to reproduce the relevant phenomena.
The behavior of concrete is assumed to be linear elastic within an elastic limit surface and up to the limit stress. With the onset of plastic strain and its further increase, material damage will rise until failure occurs. This is considered by damage parameters representing the progressive cracking and crushing of concrete caused by the abrupt stiffness changes (due to the propagation of the aircraft missile’s kinetic energy), which scale the failure surface down to a residual strength surface.
The main challenge in non-linear analyses of concrete structures is due to the relatively low tensile strength. The tensile stress causing cracking in concrete is about one tenth of its crushing strength; in reinforced concrete structures, this “weakness” is compensated by the use of reinforcement. In addition, it should be considered that the concrete material strength may increase at elevated strain rates. To account for this effect, which for concrete is particularly due to the inertial effects, a dynamic increase factor (DIF) equal to 1 is considered for steel and concrete, as the dynamic load factor associated with the impactive or impulsive loading is less than 1.2 [
22].
The implemented concrete ‘plasticity’ model accounts for material degradation in compression and in tension as impact progresses. Damage effect types are associated with cracking properties, such as critical stress, softening modulus, crushing strain, etc. It is necessary to remark that in the analyses, the compressive stiffness was recovered upon crack closure as the load changed from tension to compression, but the tensile stiffness was not recovered when the load changed from compression to tension. Equations (1) and (2) provide the stress–strain relationships for uniaxial tension and compression:
where
dt and
dtc are the tensile and compression damage parameters, respectively.
E0 is the undamaged modulus of elasticity.
εt and
εc are the tensile and compression strains, whereas
and
are the equivalent plastic strains in tension and compression.
The airplane and reinforcement steel non-linear material behavior can be modelled using an elastic-plastic material model with the von Mises yield criteria and the Cowper–Symonds strain rate method. To simulate extreme loading conditions, such as impact phenomena, material models should be strain rate-dependent. Moreover, the dynamic yield strength is considered in the Cowper–Symonds formula:
σy and
σys are the static and the dynamic yield stress, respectively;
is the strain rate; D and q are the steel material parameters, equal to 40 1/s and 5, respectively, and are taken from the literature [
23]. Further restrictions and simplification have been also imposed on the airplane velocity (assumed constant), impact direction, and internal components.
In preliminarily analyzing the ultimate structural capacity, varying material properties with age (realistic values) are taken into account. Such a variation/degradation is considered in the assessment assuming that the values of the mechanical and physical properties of both the concrete and steel—e.g., yielding strength and Young’s modulus—are equal to 80% of the nominal value.
For the purposes of the simulation, a horizontal impact onto the CC dome is assumed in order to consider the worst accident scenario (
Figure 8).
The aircraft crashing was simulated by performing a nonlinear transient analysis and assuming the updated Lagrangian procedure with follower force/stiffness command. For the current analysis, the large strain plasticity procedure is activated. The “large displacement” model option, which causes distributed loads to be based on the current deformed geometry, contributes to a stiffness effect on the tangential stiffness matrix. Erosion is not thus implemented because the model procedure allows us to account for the failure process and handle the large distortion problem (large deformation), such that the actual strength of the structure is calculated by the code by immediately deactivating the damaged elements.
The time step size was set to be equal to 10−5 s to properly capture the damage of the material as the impact progresses.
3.2. Results Discussion
Because of the strike against the target, a part of the aircraft is crushed while the remaining portion undergoes elastic deformation. The kinetic energy transferred during the impact is, thus, partially dissipated by the crushing and buckling of the plane and partially through the deformation of the walls. The impact force at the time of the crushing of the fuselage (30 ms) is about 89 MN; this value is in good agreement with the Riera impact load (represented with a dashed line in the trends of
Figure 9) and with the Reisemann experimental data [
18,
19,
20,
21,
22,
23,
24].
Upon impact, a compressive wave propagates away from the impact point; the radial stress built up is tensile in nature. When the CC material stress is greater than the ultimate dynamic strength, radial and/or circumferential cracks occur. As the compressive wave propagates outwards, the hoop or circumferential stresses generate radial cracks due to the Poisson effect.
The containment structure deforms beyond the point of permanent yielding; pronounced damage on the front and rear face occurs as a consequence of the impact.
The extensive penetration of the wall (critical condition) is reached in less time for an aged containment than for an unaged one due to the reduced strength of the materials.
Figure 10 provides the contour bands of the stress; with the same impact force transmitted to the structure, the resistant capacity of the containment reduces by about 30%. Furthermore as the results of the aircraft impact, a narrow cone forms at the entry location, confirming the wall penetration (
Figure 11). The penetration depth ranges between about 5.8 and 7.7 cm for the unaged and aged structures, respectively. In both the two conditions, the penetration is accompanied by cone cracking and some spalling of the concrete wall, while acceleration propagates faster through the structure. Even if perforation is precluded, the concentrated loading can produce significant scabbing, weakening a portion of the structure up to ten times the missile diameter in each direction and affecting the global behavior.
Figure 12 shows the acceleration propagation inside the CC structure. It can be stated that the observed dynamic response due to the aircraft impact is characterized by much higher acceleration values, especially in the high frequency range (approximately from 25 to 80 Hz), in comparison to further external event loading cases. The high accelerations at the outer CC wall in the vicinity of the impact location are transferred directly to the inner structures.
If the inner and outer structures are separated, the travelling shock waves are attenuated significantly and filtered during propagation all the way down to the foundation and, from there, to the supporting point of the equipment in the inner structure, as shown by the A3int and A2int trends.
As for vibration, the calculated acceleration of about 40 g is not sufficient to damage the internal structures, and particularly the containment vessel.
Finally, it is possible to say that, despite the 20% reduction in the materials’ properties, the strength capacity/integrity of the CC seemed not to be compromised. In fact, the phenomena of instability induced by overturning moments do not appear, nor does the strain exceed the limit value of 0.005, which is the acceptance criteria suggested by the IAEA for the compressive strain of concrete.