Next Article in Journal
Creation of AlSi12 Alloy Coating by Centrifugal Induction Surfacing with the Addition of Low-Melting Metals
Next Article in Special Issue
Valorization of Aluminum Dross with Copper via High Temperature Melting to Produce Al-Cu Alloys
Previous Article in Journal
Thermodilatometric Study of the Decay of Zeolite-Bearing Building Materials
Previous Article in Special Issue
Research on the Release of Dangerous Compounds from the BTEX and PAHs Groups in Industrial Casting Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Study on X-ray Induced Two-Dimensional Thermal Shock Waves in Carbon/Phenolic

Department of Nuclear Science and Technology, Xi’an Jiaotong University, Xi’an 710000, China
*
Author to whom correspondence should be addressed.
Submission received: 29 April 2021 / Revised: 6 June 2021 / Accepted: 22 June 2021 / Published: 25 June 2021
(This article belongs to the Special Issue Advances in Materials Processing (Second Volume))

Abstract

:
Carbon/Phenolic (C/P), a typical anisotropic material, is an important component of aerospace and often used to protect the thermodynamic effects of strong X-ray radiation. In this paper, we establish the anisotropic elastic-plastic constitutive model, which is embedded in the in-house code “RAMA” to simulate a two-dimensional thermal shock wave induced by X-ray. Then, we compare the numerical simulation results with the thermal shock wave stress generated by the same strong current electron beam via experiment to verify the correctness of the numerical simulation. Subsequently, we discuss and analyze the rules of thermal shock wave propagation in C/P material by further numerical simulation. The results reveal that the thermal shock wave represents different shapes and mechanisms by the radiation of 1 keV and 3 keV X-rays. The vaporization recoil phenomenon appears as a compression wave under 1 keV X-ray irradiation, and X-ray penetration is caused by thermal deformation under 3 keV X-ray irradiation. The thermal shock wave propagation exhibits two-dimensional characteristics, the energy deposition of 1 keV and 3 keV both decays exponentially, the energy deposition of 1 keV-peak soft X-ray is high, and the deposition depth is shallow, while the energy deposition of 3 keV-peak hard X-ray is low, and the deposition depth is deep. RAMA can successfully realize two-dimensional orthotropic elastoplastic constitutive relation, the corresponding program was designed and checked, and the calculation results for inspection are consistent with the theory. This study has great significance in the evaluation of anisotropic material protection under the radiation of intense X-rays.

1. Introduction

Carbon/Phenolic (C/P), as an advanced type of anisotropic composite material, and has been widely used in the field of aerospace [1,2]. This material has the advantages of low density, high strength, oxidation resistance, low coefficient of thermal expansion and mechanical stability at high temperature [3]. Due to the complex external environment of aviation and aerospace, such materials may be facing high-speed collision, radiation, and other dynamic loading environments [4,5]. Under intense pulsed X-ray irradiation, a large amount of energy is rapidly deposited on the surface of the material and decays rapidly inside, which causes a large temperature and pressure gradient in the material. On the other hand, the irradiated material undergoes adiabatic expansion due to a rapid increase in specific internal energy, so with large radiation that in turn has a recoil effect on the material, the material on the light-facing surface would vaporize and then be ejected. The combined effects of these factors will form unsteady stress waves (i.e., X-ray thermal shock waves) within the material [6,7]. When thermal shock waves propagate to the interface or the free surface with low impact impedance, they would unload and form sparse waves that propagate in the opposite direction. They interact with the sparse section of the incident thermal shock waves to produce tensile stress, which would crack the material or delaminate the delamination interface and thus damage the material [8]. With the wide application of C/P composite in spacecraft and missile hulls, the prediction of dynamic response and internal damage of structures caused by explosion, impact or shock wave under the action of intense pulse X-ray deposition, and the evaluation of capability composite materials to resist stress waves, arouse the study of the constitutive relation and the characteristic of stress wave propagating [9]. Reliable numerical simulation results depend on accurate constitutive models as well as effective numerical simulation programs [10]. Among a large number of mechanical numerical simulation programs, the programs used to analyze the material response and structural response related to shock-wave propagation occupy an important position, for example, CTH [11], PUFF-TFT [12], LS-DYNA [13], ABAQUS [14,15], etc. However, since the solid constitutive models of most programs are for isotropic materials, they are not suitable for anisotropic materials. Moreover, the research on the dynamic constitutive model of C/P and its application in numerical simulation is still in the process of update and further perfection. Since the 1990s, many countries have carried out research on the application of anisotropic constitutive models in numerical simulation, and some of the research results have been used in the recent version of the finite element large-scale shock dynamics numerical simulation software, such as LS-DYNA970, that embeds transversely isotropic and orthotropic elastoplastic constitutive models [13]. These models make the software capable of analyzing the dynamic response of fiber-reinforced composites. However, as mentioned above, because of anisotropic mechanical properties, C/P also has a strain rate effect. The interface effect between the fiber and the matrix, damage, temperature, etc., will affect its mechanical properties. Therefore, the establishment of a composite material dynamic constitutive model and improvement of reliability of numerical simulation results are vital to be solved. Lukyanov A et al. [16,17] used the fluid-isotropic ideal elastoplastic constitutive model, and the obtained X-ray plane frontal wave stress peak attenuation curve in the metal is in good agreement with the experiment. However, for C/P, the experimentally measured thermal shock wave stress peak decays faster than numerical simulation. We think this is likely to be the deviation caused by approximating the composite material with anisotropic model in the numerical simulation.
The mechanical properties of C/P composites depend on the properties of their phases and are obviously affected by the manufacturing process. Micromechanics provides a method to calculate the macroscopic mechanical properties of composites based on the properties and geometric shapes of their phases. Therefore, the anisotropy and strain rate correlation is the mechanical properties of most C/P materials, and they will inevitably affect the propagation law of X-ray thermal shock waves. The C/P anisotropic dynamic constitutive model should be introduced in the numerical simulation. In the calculation program. The shell shape of missiles or spacecraft is mostly cylindrical. When exposed to X-ray radiation in the air, the isotropic simplification of composite materials or the one-dimensional simplification of the model will bring errors. Based on the symmetry of the unidirectional X-ray source and the loading shell, it is feasible to adopt a two-dimensional simplified model. However, so far, there is no report about the use of the C/P anisotropic dynamic constitutive model for numerical simulation of two-dimensional X-ray thermal shock wave propagation worldwide. In this study, we used an autonomous program to establish a two-dimensional rate-dependent elastoplastic dynamic constitutive model for orthotropic materials. Taking carbon phenolic (C/P) as an example, the two-dimensional thermal shock waves generated by 1 keV and 3 keV X-ray irradiation were simulated and compared with the measured results. The experimental results show that the anisotropic propagation law of materials is basically consistent with the experimental results.

2. Theoretical Basis and Modelling Derivation

2.1. Introduction of Numerical Calculation Programme

The thermal shock wave propagation problems are caused by pulsed X-ray irradiation of anisotropic materials. The existing finite element analysis software cannot deal with these problems well. We use the inhouse code “RAMA” to program the two-dimensional finite element dynamics program [8,18] for numerical simulation. The “RAMA” program can deal with material fracture and vaporization problems and embeds a variety of constitutive models and equations of state, including orthotropic dynamic elastoplastic constitutive models. At present, it can simulate the stress wave propagation in the two-dimensional flat plate collision problem of anisotropic and isotropic materials and the two-dimensional X-ray thermal shock wave propagation in various shapes. In the following text, we will discuss in detail the important issues in the program design, such as the spindle rotation of orthotropic materials, the correction of objective stress rate in two-dimensional problems, and the simple algorithm of X-ray energy deposition and so on.

2.2. Two-Dimensional Orthotropic Elastoplastic Constitutive Model

Considering the plane strain (taking 1–2 plane strain as an example), there are three main axis directions (expressed by 1, 2, and 3) for orthotropic materials: the fiber warp (warp), the fiber weft (fill), and the thickness direction (thickness) ε 33 = ε 13 = ε 23 = 0 ,   σ 13 = σ 23 = 0 . The stress–strain relationship is described by Hooke’s law in the elastic denaturation.
{ σ 11 = c 11 ε 11 + c 12 ε 22 σ 22 = c 12 ε 11 + c 22 ε 22 σ 33 = c 13 ε 11 + c 23 ε 22 σ 12 = c 44 ε 12
c i j is the stiffness matrix coefficient related to the elastic modulus, Poisson ratio, and shear modulus of the material.
There is no corresponding relationship between both strain and stress whose state is related to the deforming path or process in plastic deformation, but the Hooke law is satisfied between stress increment and elastic strain increment. Thus, the constitutive relation is expressed in an incremental form. The strain increment can be decomposed into an elastic strain increment and plastic strain increment, d ε i j = d ε i j e + d ε i j p , as shown below:
[ d σ 11 d σ 22 d σ 33 d σ 12 ] = [ c 11 , c 12 , c 13 , 0 c 12 , c 22 , c 23 , 0 c 13 , c 23 , c 33 , 0 0 , 0 , 0 , c 44 ] [ d ε 11 e d ε 22 e d ε 33 e d ε 12 e ] = [ c 11 , c 12 , c 13 , 0 c 12 , c 22 , c 23 , 0 c 13 , c 23 , c 33 , 0 0 , 0 , 0 , c 44 ] [ d ε 11 d ε 11 p d ε 22 d ε 22 p d ε 33 p d ε 12 d ε 12 p ]
For isotropic materials, the decoupling of the tolerance and distortion laws by most impingement dynamics programs could be carried out by algorithm, i.e., calculating the static pressure and the partial stress by the state equation and the constitutive relation, respectively. Because of the possible changes in shape and volume by hydrostatic pressure and stress bias, respectively, anisotropic materials cannot be simply decoupled. To reflect the difference from isotropic conditions, the hydrostatic pressure in isotropic materials would be replaced by the mean normal stress p ( p = ( σ 11 + σ 22 + σ 33 ) / 3 ) will be used to replace the concept of hydrostatic pressure in isotropic materials. In this paper, stress and strain are defined as positive pressure and negative tension; that is, the average normal stress and deviator stress of orthotropic materials. It can be decoupled in this form to facilitate its application in calculation programs. The stress ( σ i j ) is decomposed into mean positive stress ( p ) and partial stress ( S i j ), as represented by   σ i j = p δ i j + S i j . The strain ( ε i j ) is decomposed into body strain ( θ ( θ = ε 11 + ε 22 ) ) and partial strain ( e i j ), as represented by   ε i j = 1 3 θ δ i j + e i j . The basic form of the average positive stress at the stage of elastic deformation is shown below:
p = 1 3 ( σ 11 + σ 22 + σ 33 ) p = ( c 11 + 3 c 12 + c 13 + 2 c 22 + 2 c 23 ) θ 9 + ( c 11 + c 13 c 22 c 23 ) e 11 3
The basic form of average positive stress increment during plastic deformation is as follows:
d p = ( c 11 + 3 c 12 + c 13 + 2 c 22 + 2 c 23 ) d θ 9 + ( c 11 + c 13 c 22 c 23 ) d e 11 3 ( c 11 + c 12 + c 13 ) d ε 11 p 3 ( c 12 + c 22 + c 23 ) d ε 22 p 3 ( c 13 + c 23 + c 33 ) d ε 33 p 3
In order to make the expression of average normal stress reflect the nonlinear effect of volume change and the anisotropic characteristics of the material, a modified equation of state can be introduced for calculation [19]. Since the X-ray irradiation process not only describes the shock compression state with relatively low temperature but also deals with the state formed by the coupling of X-ray energy deposition and hydrodynamic motion, this work uses the PUFF Equations (5) and (6) of state to describe the specific form [20,21].
Compression zone:
{ p = p H + ρ Γ ( E E H ) p H = ρ 0 c 0 2 ( 1 v / v 0 ) [ 1 s ( 1 v / v 0 ) ] 2 , E H = 1 2 p H ( v 0 v )
Expansion zone:
p = ρ [ γ 1 + ( Γ 0 γ + 1 ) ρ ρ 0 ] [ E E S ( 1 exp ( N ρ 0 ρ ( 1 ρ 0 ρ ) ) ) ]
v is the specific volume, ρ is the density, ρ = 1 / v , Γ is the Grüneisen coefficient, and c 0 ,   s are the fitting experimental parameters in the rational expression of D = c 0 + s u of the material shock wave velocity D and the post-wave particle velocity u. γ is the specific heat ratio, and E s is the sublimation energy, N = c 0 2 / Γ E s .
Expand Equations (5) and (6) into a polynomial series about θ ( θ = ε 11 + ε 22 ( v 0 v ) / v ) and combine Equation (3) to obtain a correction that can not only reflect the nonlinear characteristics of the material volume change but also reflect the anisotropic strength effect of the material’s equation of state [4]. From this, the corrected average normal stress of the compression zone and expansion zone in the elastic deformation stage is as follows:
p = A 1 θ + [ A 2 Γ 2 A 1 ] θ 2 + [ A 3 Γ 2 A 2 ] θ 3 + ( ρ 0 Γ + ρ 0 Γ θ ) E + 1 3 ( c 11 + c 13 c 22 c 23 ) e 11
p = B 1 θ + B 2 θ 2 + B 3 θ 3 + ( ρ 0 Γ + 3 2 ρ 0 Γ θ ) E γ 1 2 ρ 0 θ E + 1 3 ( c 11 + c 13 c 22 c 23 ) e 11
Therefore, the modified mean normal stress increment of the compression zone and the expansion zone in the plastic deformation stage is obtained by Equation (4).
d p = A 1 d θ + 2 [ A 2 Γ 2 A 1 ] θ d θ + 3 [ A 3 Γ 2 A 2 ] θ 2 d θ + ρ 0 Γ E d θ + ( ρ 0 Γ + ρ 0 Γ θ ) d E + 1 3 ( c 11 + c 13 c 22 c 23 ) d e 11 1 3 ( c 11 + c 12 + c 13 ) d e 11 p 1 3 ( c 12 + c 22 + c 23 ) d e 22 p 1 3 ( c 13 + c 23 + c 33 ) d e 33 p
d p = B 1 d θ + 2 B 2 θ d θ + 3 B 3 θ 2 d θ + ( ρ 0 Γ + 3 2 ρ 0 Γ θ γ 1 2 ρ 0 θ ) d E ( 3 2 ρ 0 Γ E γ 1 2 ρ 0 E ) d θ + 1 3 ( c 11 + c 13 c 22 c 23 ) d e 11 1 3 ( c 11 + c 12 + c 13 ) d e 11 p 1 3 ( c 12 + c 22 + c 23 ) d e 22 p 1 3 ( c 13 + c 23 + c 33 ) d e 33 p
A 1 = ρ 0 c 0 2 , A 2 = ρ 0 c 0 2 ( 2 s 1 ) , A 3 = ρ 0 c 0 2 ( 3 s 2 4 s + 1 ) , B 1 = ρ 0 c 0 2 , B 2 = B 1 2 γ 1 2 Γ B 1 + B 1 2 N , B 2 = 5 B 1 24 + 5 8 γ 1 Γ [ 5 4 + 1 4 γ 1 Γ ] B 1 N + 1 6 B 1 N 2 A 1 = B 1 = 1 9 ( c 11 + 3 c 12 + c 13 + 2 c 22 + 2 c 33 )
A 1 ( B 1 ) reflects the material anisotropy or equivalent volume contact and is reduced to volume modulus ( A 1 ( B 1 ) ) under the isotropic limit condition.
The partial stress ( S i j ) of an elastic stage can be obtained by S i j = σ i j [ ( σ 11 + σ 22 + σ 33 ) / 3 ] δ i j and (1), for example,
S 11 = ( 2 c 11 + 3 c 12 c 13 2 c 22 2 c 23 ) θ 9 + ( 2 c 11 3 c 12 c 13 + c 22 + c 23 ) e 11 9
The increment of partial stress ( d S i j ) in the plastic stage can be obtained according to Equation (2) and d S i j = d σ i j [ ( d σ 11 + d σ 22 + d σ 33 ) / 3 ] δ i j , as shown as follows:
d S 11 = ( 2 c 11 + 3 c 12 c 13 2 c 22 2 c 23 ) d θ 9 + ( 2 c 11 3 c 12 c 13 + c 22 + c 23 ) d e 11 9 ( 2 c 11 c 12 c 13 ) d e 11 p 9 ( 2 c 12 2 c 22 2 c 23 ) d e 22 p 9 ( 2 c 13 c 23 c 33 ) d e 33 p 9
Other expressions of partial stress components can be obtained in the same way. The plastic strain increment of the above formula needs to be solved by yield criterion, and the Tsai–Hill yield criterion can be adopted to determine whether the material has entered the plastic stage for orthogonal anisotropic materials.
F = σ 11 2 Y 11 2 + σ 22 2 Y 22 2 + σ 33 2 Y 33 2 + σ 12 2 Y 12 2 + Y ¯ 33 σ 11 σ 22 + Y ¯ 22 σ 11 σ 33 + Y ¯ 11 σ 22 σ 33 1 = 0
Y ¯ 11 = 1 Y 11 2 1 Y 22 2 1 Y 33 2 , Y ¯ 22 = 1 Y 22 2 1 Y 11 2 1 Y 33 2 , Y ¯ 33 = 1 Y 33 2 1 Y 22 2 1 Y 11 2 , Y 11 , Y 22 , Y 33 , Y 12 are the yield strength in 3 main material directions and the shear yield strength in 1–2 planes.
By the strain rate correlation [17,18], [ ( F + 1 ) / R 2 ] 1 = 0 is obtained. R ( ε ) = 1 + β ln ( ε / ε 0 ) is the strain rate factor, β is the experiment parameters, and ε 0 is the reference strain rate. The strain rate ( ε ) is considered the equivalent plastic strain rate ( ε p = 2 3 ε i j p ). Thus, the yield criterion of strain rate separation is f = ε 0 exp [ 1 β F + 1 1 ) ] ε p = 0 . According to the orthogonal law, it is shown as follows:
ε ˙ i j p = λ ˙ f σ i j = λ ˙ ε 0 exp [ 1 β F + 1 1 ) ] ( 1 2 β F + 1 ) D i j
where D i j = F σ i j , λ is the plastic flow factor. Based on ε p = 2 3 ε i j p = 2 3 λ 2 ε i j p , we can obtain ( λ = ε p / 2 3 f σ i j ). According to the consistency law, the stress state is always on the yield surface when plastic deformation occurs, and we can obtain ε p = ε 0 exp [ 1 β F + 1 1 ) ] by f = 0 , as shown below:
λ = ε ˙ p 2 3 f σ i j = 2 β F + 1 2 3 D i j   and   D i j 2 = D 11 2 + D 22 2 + D 33 2 + 2 D 12 2
The plastic strain rate can be solved based on Equation (14), which is equivalent to the plastic strain increment. Based on this, the average positive stress and partial stress in the elastoplastic stage can be obtained according to corresponding formulas. Such an elastic-plastic constitutive model takes into account not only the rate-related effect but also the nonlinear characteristics of volume change and the anisotropic strength effect of materials in the process of compression and expansion. Thus, it can objectively reflect real material changes.

2.3. Energy Deposition of Unit Mass

The interaction between X-rays and matter is essentially the interaction between X-ray photons and matter atoms. The main interaction process is photoelectric effect and Compton scattering, accompanied by fluorescence effect and secondary scattering effect. A parallel X-ray incident vertically along the x-direction onto a 2-D plate should be considered. If φ 0 is the initial energy flux, the spectral energy after the distance x is as follows:
φ ( x ) = φ 0 0 f ( λ , T ) exp [ μ ( λ ) ρ x ] d λ
where f ( λ , T ) = c 1 λ 5 1 exp ( c 2 / λ T ) 1 is the Plank function. c 1 = 3.7416 × 10 16   W m 2 and c 2 = 0.014388   m K are the first and second radiation constants. λ is the X-ray wavelength, T is the blackbody spectrum temperature, ρ is the material density, and μ ( λ ) as the mass absorption coefficient for the material wavelength ( λ ) depending on the photoelectric effect and Compton scattering.
The continuous energy spectrum is decomposed into several monochromatic lights, and the upper and lower wavelengths should be cut off reasonably for the convenience of numerical calculation. Therefore, Equation (16) can be translated into φ ( x ) = φ 0 w j exp [ μ ( λ ) ρ x ] , where w j is the percentage of the total energy spectrum by the group j monochromatic light. So, the energy deposition of unit mass ( x ( x + Δ x ) , y ( y + Δ y ) ) can be calculated by:
Δ E = [ φ ( x ) φ ( x + Δ x ) ] Δ y ρ Δ x Δ y = φ ( x ) φ ( x + Δ x ) ρ Δ x
The X-ray energy spectrum and time spectrum are shown in Figure 1. B(E) represents the normalized energy spectrum, and η(t) is the time spectrum, which provides input conditions for subsequent numerical simulation calculations of strong pulse X-rays.

3. Experimental Comparison of C/P Materials Irradiated by Low Energy X-rays

In order to verify the two-dimensional orthotropic elastoplastic constitutive model and the energy deposition derivation, a low-energy X-ray irradiation experiment of C/P materials was carried out. The X-ray experimental device is a combined multi-purpose low-energy high-current pulsed electron beam accelerator, which is mainly used to generate pulsed gamma rays and X-rays. According to the energy spectrum and energy fluence output, it can be divided into three X-ray working states. The experiment uses an X-ray mode with an average energy of 0.1–1.5 keV, a pulse width of 40 ns, and a pulse beam energy of 65 KJ. Use Polyvinylidene Fluoride (PVDF) to measure the thermal shock wave of X-rays on carbon phenolic materials. PVDF is a semi-crystalline polymer. The PVDF is directly placed in the target component to directly measure the evolution of the thermal shock wave. As shown in Figure 2, the thickness of the self-made PVDF is about 70 μm, and the influence on the propagation behavior of the thermal shock wave is negligible. PVDF adopts the current mode and connects a load resistance directly to the output end of the membrane, and the signal is transmitted to the oscilloscope to record by the load resistance cable. The oscilloscope adopts high resistance gear, and the shunt of the signal can be ignored. The load R = 50 Ω, the characteristic impedance of the signal cable is 50 Ω, and the two are matched. The relationship between the voltage V(t) recorded by the oscilloscope and the thermal shock wave stress output of PVDF is as follows:
σ ( t ) = 1 AKR 0 t V ( t ) d t
where σ ( t ) is the PVDF thermal shock wave stress output, V(t) is the experimental recording voltage, A is the effective area of the PVDF, K is the dynamic piezoelectric coefficient given by the calibration experiment, and R is the load resistance. At this point, the thermal shock wave waveform can be obtained by integrating the numerical value of the voltage signal recorded by the oscilloscope.
Figure 3 shows the ablation of the material surface before and after irradiation. After the light-facing surface of the material is irradiated by the electron beam, a large amount of energy is deposited on the surface and converted into heat energy of the substance. The carbon phenolic matrix does not melt after being heated but transforms into a layer of porous carbon. The carbon particles will fly away from the target surface under the action of the sparse wave on the facing surface. Carbon fiber does not melt and only sublimates when its ratio can reach sublimation energy. Therefore, some sublimated carbon fibers will remain on the edge of the target under the ’scouring’ of the jet of carbon particles (as shown in Figure 3). After the experiment, it was found that, from the recovered target, there was spallation on the back side, and some sensors were degummed at the bonding point within the energy fluence range of the experiment. Figure 3 verifies the propagation and effect analysis of thermal shock waves in C/P materials and provides data for the theoretical derivation of numerical calculations.
Table 1 shows the results of five experiments. E is the X-ray average energy, F is the energy fluence, σ is the thermal shock wave stress, and the PVDF installation position is 3 mm, 6 mm, and 9 mm away from the illuminated surface.
Table 1 and Figure 4 show the measured peak stress of the thermal shock wave and the fitted peak attenuation curve and reveal the relationship between the peak stress of the thermal shock wave and the energy fluence at 3, 6, and 9 mm from the target surface in the C/P material. The results show the peak value of the thermal shock wave increases with the increase in energy fluence, but its growth rate not consistent. When the energy fluence is low, the peak stress of the thermal shock wave increases slowly with the increase in the energy fluence, and when the energy fluence is high, its growth accelerates significantly.
Figure 5 shows the comparison of thermal shock wave stress experiments and numerical simulation data at 3, 6, and 9 mm in the X direction of C/P materials. The numerical simulation data in Figure 4 uses the calculation results of 1 keV X-ray with similar energy fluence. The experimental energy fluence in the Figure is 383 J/cm2, and the numerical simulation energy fluence is 380 J/cm2. There is not much difference in between, and the experimental thermal shock wave stress peak value of the same C/P material position is greater than the simulation result. The arrival time of the thermal shock wave coincides with the peak time. When the high current electron accelerator X-ray simulation source is working, the average X-ray energy production is low, which is only 0.4–0.7 keV, while the average energy of the strong pulse 1 keV blackbody spectrum X-ray in the numerical simulation is 2.7 keV. That is because the surface layer of the material has vaporization, the combined effect of vaporization recoil and thermal deformation, the thermal shock wave is mainly manifested as a compression wave. Therefore, the peak pressure of the surface layer is very large, and the stress decays quickly, causing the peak value of the experimental result to be larger than the calculation result. Although the experimental data is limited, it can be determined that the 1 keV blackbody spectrum acts on the surface thermal shock wave stress of the C/P material to be larger (Figure 5). In the experiment, the thermal shock wave stress decay rate is relatively low, so that the experimental thermal shock wave stress is larger than the calculated result when the depth is 3, 6, and 9 mm.
Theoretical analysis shows that under the same energy fluence, the simulated coupling coefficient of X-ray and intense pulsed X-ray in the experiment is greater than 1 keV blackbody spectrum X-ray and intense pulsed X-ray. This is related to the shallower X-ray energy deposition and less vaporized recoil material in the experiment. Therefore, the stress peak value of the thermal shock wave at the deeper position of the X-ray generated by the strong current electron accelerator is greater than the simulation result, which is in line with the law of thermal shock wave generation and propagation.
Through low-energy X-ray experimental research, a qualitatively consistent result with the 1 keV simulation result was obtained, indicating that the calculation method is credible.

4. Calculation and Analysis of Thermal Shock Wave for Carbon Phenolic (C/P) Materials

Take carbon phenolic (hereinafter referred to as C/P) plate as an example to simulate the propagation of thermal shock waves under the action of pulsed X-rays. The C/P material is a fiber-reinforced composite material. Its three main axis directions are the warp, fill, and thickness of the fiber cloth. 1, 2, and 3 are used to indicate the material main axis and x, y, and z to represent the system coordinate system. Initially, the two coordinate systems coincide, and the x-direction is taken as the material thickness direction (Thickness), the y-direction is taken as the material warp direction (Warp), and the z-direction is taken as material weft (Fill). As shown in Figure 6a. The size of the C/P board in the x-direction is 1 cm, the size in the y-direction is 2 cm, and the size in the z-direction is much larger than the size in the x and y-directions. X-rays are incident parallel to the x-direction, the irradiation problem can be simplified to a plane strain problem, and the model is shown in Figure 6b.
The parameters of the state equation are shown in Table 2. The main constitutive parameters are shown in Table 3. The mesh size is 0.0025 cm × 0.01 cm. According to the X-ray spectrum, the blackbody spectra of kT = 1 keV and 3 keV are obtained. The initial energy flux is 418 J/cm2, and the irradiation time spectrum is a rectangular spectrum with a width of 0.1 μs.
Figure 7 shows the energy deposition profile in carbon phenolic calculated by Equation (17). It can be seen that the energy deposition of 1 keV and 3 keV of X-ray blackbody spectrum decays exponentially, and the peak of energy deposition of soft X-ray of 1 keV is high, the deposition depth is shallow. The 3 keV hard X-ray energy deposition peak is low, but the deposition depth is large.
Figure 8 shows the time and space distribution diagrams of σxx along the x-direction on the symmetry axis y = 1 cm in the C/P material, and Figure 9 and Figure 10 show the σxx contour diagrams in the two-dimensional plane.
As shown in Figure 8, the waveform of thermal shock waves is almost triangular, and the full width at half maximum (FWHM) of thermal shock wave expands slightly with wave propagation, with a higher pressure near the material surface and the peak stress value decreasing with wave propagation under different hardness spectrum X-ray irradiation. However, due to the shallow ray penetration under the soft X-ray irradiation of 1 keV, the energy is mainly concentrated near the surface layer. The result shows that the surface layer of the material is vaporized, and the thermal shock wave is mainly expressed as a compression wave due to the combined action of vaporization recoil and thermal deformation. Thus, the peak pressure on the surface is so high, and the attenuation is extremely fast. When the X-ray penetration is deeper, the material would not vaporize under the hard X-ray irradiation of 3 keV. The thermal shock wave is mainly caused by thermal deformation, and the tensile wave is produced after the compression wave. If the X-ray spectrum is hard enough and the flux is large enough, the tensile wave would be much stronger. Therefore, its propagation can cause tensile fracture at any position in the target.
It can be seen from Figure 9 and Figure 10 that the penetration of 3 keV hard X-ray in the material is much deeper than that of 1 keV soft X-ray. The material vaporization takes place under the X-ray irradiation of 1 keV instead of 3 keV, and its tensile strength is much greater. Because it is a 2-D model with an obvious sparse free boundary on its upper and lower sides, the stress wave propagation shows obvious 2-D characteristics.

5. Conclusions

In this paper, we use an in-house program containing an orthotropic two-dimensional dynamic elastoplastic constitutive model, take C/P materials as an example, and simulate pulsed X-ray induced thermal shock waves in anisotropic materials. The following conclusions are obtained.
(1)
The adopted elastoplastic constitutive model not only considers rate-dependent effects but also considers the nonlinear characteristics of volume changes during compression, expansion, and the anisotropic strength effects of materials, which can objectively reflect the true changes of materials. The theoretical predictions are in good agreement with the experimental results.
(2)
Using the X-rays generated by the high current electron accelerator to irradiate the surface of the C/P material, the stress peaks at the thermal shock wave measurement points are larger than the numerical simulation results, which are related to the X-ray energy deposition and stress attenuation rate.
(3)
X-ray-induced thermal shock waves in C/P materials have two mechanisms, namely thermal deformation and vaporization recoil. The vaporization recoil phenomenon appears under 1 keV X-ray irradiation, which is mainly manifested as a compression wave. Under 3 keV X-ray irradiation, X-ray penetration is deeper, and the induced thermal shock wave is caused by thermal deformation, resulting in a strong stretching phenomenon, which becomes the main cause of damage to the material.
(4)
The thermal shock waveforms in C/P under 1 keV and 3 keV X-ray irradiation, such as the peak value of thermal shock wave, X-ray penetration depth, vaporization phenomenon, tensile strength, etc., are very different. Thermal shock wave propagation exhibits two-dimensional characteristics. The energy deposition of 1 keV and 3 keV both decays exponentially. The energy deposition peak of 1 keV soft X-ray is high, and the deposition depth is shallow; the hard X-ray energy deposition peak of 3 keV is low, and the deposition depth is large.

Author Contributions

Conceptualization, D.W.; methodology, J.W.; validation, J.W. and H.L.; formal analysis, D.W. and Y.G.; investigation, D.W. and Y.G.; data curation, D.W.; writing—original draft preparation, D.W.; writing—review and editing, D.W., Y.G., S.W., J.W. and H.L.; supervision, S.W.; project administration, S.W.; funding acquisition, S.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the key project of the Intergovernmental International Scientific and Technological Innovation Cooperation in China under Grant No.2016YFE0128900 and the National Natural Science Foundation of China under Grant No.11775166.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data included in this study are available upon request by contact with the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. White, S. Radiation Testing of PICA at the Solar Power Tower. In Proceedings of the 10th AIAA/ASME Joint Thermophysics and Heat Transfer Conference, Chicago, IL, USA, 28 June–1 July 2010; p. 28. [Google Scholar]
  2. Du, Z.; Zhu, M.; Wang, Z.; Yang, J. Design and application of composite platform with extreme low thermal deformation for satellite. Compos. Struct. 2016, 152, 693–703. [Google Scholar]
  3. Xu, Y.Y. Thermal-Chemical Coupling Analysis of Low-Density Carbon/Phenolic Composites under Ablation Condition. Master’s Thesis, Harbin Institute of Technology, Harbin, China, 2016. [Google Scholar]
  4. Ryan, S.; Schaefer, F.; Riedel, W. Numerical simulation of hypervelocity impact on CFRP/Al HC SP spacecraft structures causing penetration and fragment ejection. Int. J. Impact Eng. 2006, 33, 703–712. [Google Scholar] [CrossRef]
  5. Min, S.; Chen, X.; Chai, Y.; Lowe, T. Effect of reinforcement continuity on the ballistic performance of composites reinforced with multiply plain weave fabric. Compos. Part B Eng. 2016, 90, 30–36. [Google Scholar] [CrossRef]
  6. Wang, D.; Zhang, K.; Tang, W. Numerical Simulation of Thermal Shock Waves Induced by Pulsed X-ray in C/PF Materials. J. Phys. Conf. Ser. 2021, 1865, 022066. [Google Scholar] [CrossRef]
  7. Tang, W.H.; Zhang, R.Q.; Zhao, G.M. Thermal shock wave induced by impulsive X-ray. Chin. J. High Press. Phys. 1995, 9, 107–110. [Google Scholar]
  8. Zhou, N.; Qiao, D. Materials Dynamics under Pulse Beam Radiation; National Defence Industry Press: Beijing, China, 2002; pp. 401–404. [Google Scholar]
  9. Huang, X.; Tang, W.H.; Jiang, B.H. Influences of constitutive models on numerically simulated results of X-ray thermal shock waves in composite materials. Explos. Shock Waves 2011, 31, 600–605. [Google Scholar]
  10. Zhang, R.Q.; Tang, W.H.; Zhao, G.M. Several influential factors ON numerical simulated results for the X-ray thermal shock wave. Chin. J. High Press. Phys. 1998, 12, 161–167. [Google Scholar]
  11. Hertel, E.S.; Bell, R.L.; Elrick, M.G.; Farnsworth, A.V.; Kerley, G.I.; McGlaun, J.M.; Petney, S.V.; Silling, S.A.; Taylor, P.A.; Yarrington, L. CTH: A Software Family for Multi-Dimensional Shock Physics Analysis; Springer: Berlin/Heidelberg, Germany, 1995; pp. 377–382. [Google Scholar]
  12. Flock, R.A.; Liu, D.T. Numerical analysis of “thin film” Gruneisen test measurements. In Proceedings of the American Physical Society Topical Conference, Williamsburg, VA, USA, 17–20 June 1991; pp. 403–406. [Google Scholar]
  13. Grujicic, M.; Pandurangan, B.; Koudela, K.; Cheeseman, B. A computational analysis of the ballistic performance of light-weight hybrid composite armors. Appl. Surf. Sci. 2006, 253, 730–745. [Google Scholar] [CrossRef]
  14. Clegg, R.; White, D.; Riedel, W.; Harwick, W. Hypervelocity impact damage prediction in composites: Part I—Material model and characterisation. Int. J. Impact Eng. 2006, 33, 190–200. [Google Scholar] [CrossRef]
  15. Riedel, W.; Nahme, H.; White, D.; Clegg, R. Hypervelocity impact damage prediction in composites: Part II—Experimental investigations and simulations. Int. J. Impact Eng. 2006, 33, 670–680. [Google Scholar] [CrossRef]
  16. Lukyanov, A.A. Constitutive behaviour of anisotropic material under shock loading. Int. J. Plast. 2008, 24, 140–167. [Google Scholar] [CrossRef]
  17. Lukyanov, A.A. An equation of state for anisotropic solids under shock loading. Eur. Phys. J. B 2008, 64, 159–164. [Google Scholar] [CrossRef]
  18. Qiao, D. Thermodynamic Effect and Reinforcing Technology under Pulse X-ray Radiation; National Defense Industry Press: Beijing, China, 2012; pp. 118–196. [Google Scholar]
  19. Nagai, M.; Ikeda, T.; Miyazaki, N. Stress intensity factor analysis of an interface crack between dissimilar anisotropic materials under thermal stress using the finite element analysis. Int. J. Fract. 2007, 146, 233–248. [Google Scholar] [CrossRef]
  20. Khan, A.S.; Liu, H. Strain rate and temperature dependent fracture criteria for isotropic and anisotropic metals. Int. J. Plast. 2012, 37, 1–15. [Google Scholar] [CrossRef]
  21. Castres, M.; Berthe, J.; Brieu, M.; Deletombe, E. A strain rate and temperature dependent criterion to describe the linear nonlinear-behaviour’s transition of organic matrix composite materials in shear: Application to T700GC/M21. Mech. Mater. 2018, 124, 100–105. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The normalized energy spectrum and time spectrum of intense pulse X-ray. (a) energy spectrum; (b) time spectrum.
Figure 1. The normalized energy spectrum and time spectrum of intense pulse X-ray. (a) energy spectrum; (b) time spectrum.
Materials 14 03553 g001
Figure 2. A schematic diagram of PVDF thermal shock wave stress measurement.
Figure 2. A schematic diagram of PVDF thermal shock wave stress measurement.
Materials 14 03553 g002
Figure 3. Experimental photos before and after X-ray irradiation.
Figure 3. Experimental photos before and after X-ray irradiation.
Materials 14 03553 g003
Figure 4. The measured peak stress of thermal shock wave and fitted peak attenuation curve.
Figure 4. The measured peak stress of thermal shock wave and fitted peak attenuation curve.
Materials 14 03553 g004
Figure 5. A comparison of C/P material experiment and calculation data.
Figure 5. A comparison of C/P material experiment and calculation data.
Materials 14 03553 g005
Figure 6. Schematic diagrams of the C/P material model: (a) principal directions and lamina of the composite and (b) a simplified model for the numerical simulation.
Figure 6. Schematic diagrams of the C/P material model: (a) principal directions and lamina of the composite and (b) a simplified model for the numerical simulation.
Materials 14 03553 g006
Figure 7. The deposited energy distribution along the symmetry axis under irradiation of X-ray with different blackbody spectra. (a) kT = 1 keV; (b) kT = 3 keV.
Figure 7. The deposited energy distribution along the symmetry axis under irradiation of X-ray with different blackbody spectra. (a) kT = 1 keV; (b) kT = 3 keV.
Materials 14 03553 g007
Figure 8. The stress wave along symmetry axis of C/P under radiation of X-ray with different blackbody spectra. (a) time distribution, kT = 1 keV; (b) time distribution, kT = 3 keV; (c) space distribution, kT = 1 keV; (d) space distribution, kT = 3 keV.
Figure 8. The stress wave along symmetry axis of C/P under radiation of X-ray with different blackbody spectra. (a) time distribution, kT = 1 keV; (b) time distribution, kT = 3 keV; (c) space distribution, kT = 1 keV; (d) space distribution, kT = 3 keV.
Materials 14 03553 g008
Figure 9. σxx (GPa) contours of C/P under radiation of X-ray with 1 keV blackbody spectrum. (a) t = 1 us; (b) t = 2 us.
Figure 9. σxx (GPa) contours of C/P under radiation of X-ray with 1 keV blackbody spectrum. (a) t = 1 us; (b) t = 2 us.
Materials 14 03553 g009
Figure 10. σxx (GPa) contours of C/P under radiation of X-ray with 3 keV blackbody spectrum. (a) t = 1 us; (b) t = 2 us.
Figure 10. σxx (GPa) contours of C/P under radiation of X-ray with 3 keV blackbody spectrum. (a) t = 1 us; (b) t = 2 us.
Materials 14 03553 g010
Table 1. The X-ray thermal shock wave experiment results.
Table 1. The X-ray thermal shock wave experiment results.
NO.E /keVF /(J/cm2)σ/GPa
3 mm6 mm9 mm
10.4973831.5491.2091.087
20.4223120.8070.7030.598
30.4613580.8760.6670.512
40.6975312.2231.5331.270
50.5874291.8211.3521.139
Table 2. State parameters of C/P composites.
Table 2. State parameters of C/P composites.
ρ 0   ( g / cm 3 ) c 0   ( km / s ) s Γ 0 γ e s   ( kJ / g )
1.382.351.662.321.45.15
Table 3. Constitutive model parameters of C/P composites.
Table 3. Constitutive model parameters of C/P composites.
Elastic Module
(GPa)
Shear Module (GPa)Poisson RatioYield Strength
(GPa)
Rate Parameter
E x E y E z σ x y v x y v x z v y z Y x x 0 Y y y 0 Y z z 0 Y x y 0 β
6.965.454.873.50.300.400.3130.120.0630.170.070.0218
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, D.; Gao, Y.; Wang, S.; Wang, J.; Li, H. Study on X-ray Induced Two-Dimensional Thermal Shock Waves in Carbon/Phenolic. Materials 2021, 14, 3553. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14133553

AMA Style

Wang D, Gao Y, Wang S, Wang J, Li H. Study on X-ray Induced Two-Dimensional Thermal Shock Waves in Carbon/Phenolic. Materials. 2021; 14(13):3553. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14133553

Chicago/Turabian Style

Wang, Dengwang, Yong Gao, Sheng Wang, Jie Wang, and Haipeng Li. 2021. "Study on X-ray Induced Two-Dimensional Thermal Shock Waves in Carbon/Phenolic" Materials 14, no. 13: 3553. https://0-doi-org.brum.beds.ac.uk/10.3390/ma14133553

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop