Published: 27 September 2024

Investigation of the load’s variation rate effect on the efforts and moment intensity factors values in a cylindrical shell with a longitudinal crack

Mykola Makhorkin1
Tetiana Makhorkina2
Ihor Demkiv3
Andrii Kunynets4
1Pidstryhach Institute for Applied Problems of Mechanics and Mathematics National Academy of Sciences of Ukraine, 3-b, Naukova Str. 79060, Lviv, Ukraine
2, 3, 4, 1Department of Computational Mathematics and Programming, Lviv Polytechnic National University, Lviv, Ukraine
Corresponding Author:
Mykola Makhorkin
Views 18
Reads 9
Downloads 32

Abstract

The limit equilibrium of a cylindrical shell with a longitudinal crack under the action of a time-varying load is considered. Dynamic coefficients of intensity of forces and moments near the tips of cracks in shells made of steel, copper, and cast iron are calculated. The influence of the mechanical characteristics of the shell material and the rate of load variation on their behavior is analyzed.

1. Introduction

Cylindrical shells are important structural elements that are widely used in various industries. They are used in civil engineering, mechanical engineering, aircraft and rocketry and the chemical industry, particularly in the design of tanks, reservoirs, pipelines, machinery covers, etc., that are subject to various external loads. Such structural elements in modern engineering usually operate under the influence of a time-varying load, which requires taking into account the inertia of the shell material when studying its limit equilibrium.

At present, methods of solving problems of the stress-strain state of different kinds of shells with cracks under static loading are described in detail in numerous works [1-3]. Many fewer scientific papers deal with the limiting equilibrium of shell structures subjected to dynamic loading [4-8]. This situation is caused by the need to consider the inertia of the material in the equilibrium equations under dynamic loading, which leads to difficulty in constructing their solutions. That is why the study of the limiting equilibrium of cylindrical shells under the action of a time-varying load is an actual problem of shell theory.

This article is based on the author’s previous research. In the works of M. Nykolyshyn and M. Makhorkin [9, 10] the problem of the limit equilibrium of a long cylindrical shell with a longitudinal crack under the action of a load varying in time with regard to the inertia of the material was considered. In [9], the methodology for constructing a system of singular integral equations (SIE) to determine the coefficients of intensity of forces and moments at the tips of a crack is described, and in [10], such a system is constructed and ways to solve it are proposed. The following is a study of the effect of load's variation rate on the values of the intensity coefficients of forces and moments in cylindrical shells with a longitudinal crack based on the relations obtained in [10] and the influence of shell parameters (physical and geometrical) on its limit equilibrium is analyzed.

2. Problem formulation

An infinite cylindrical shell of radius R with a through longitudinal crack of length 2l was considered (Fig. 1). The shell is under the action of a surface loading that varies according to an exponential law eγ~τ and the location of the crack does not violate the Saint-Vernand principle

According to [2, 3], as it was realized in [9, 10], we write:

– The corresponding efforts that will take place in the shell:

1
N2iα,β,τ=N2i*α,βeγ~τ, Siα,β,τ=Si*α,βeγ~τ,
M2i(α,β,τ)=M2i*(α,β)eγ~τ, Q2i*(α,β,τ)=Q2i**(α,β)eγ~τ, i=0,1.

– Stress-free deformation field:

2
εββ0=εββ*eγ~τ=R-1vα,τδβ, εαβ0=εαβ*eγ~τ=R-1uα,τδα,
æββ0=æββ*eγ~τ=-R-1θβα,τδβ-R-2wα,τβδβ,
æαβ0=æαβ*eγ~τ=-R-2α[w(α,τ)] δ(β), εαα0=æαα0=0.

– Conditions on the crack edges [10]:

3
N2+α,0,τ=N2-α,0,τ=f1αeγ~τ, S+α,0,τ=S-α,0,τ=f2αeγ~τ,
M2+α,0,τ=M2-α,0,τ=f3αeγ~τ, Q2+*α,0,τ=Q2-*α,0,τ=f4αeγ~τ,

where:

4
f1α=N21*α,+ 0-N20*α,+ 0, f2α=S1*α,+ 0-S0*α,+ 0,
f3α=M21*α,+ 0-M20*α,+ 0, f4α=Q21**α,+ 0-Q20**α,+ 0,

where εij*, æij*, i,j=α,β are the coordinate functions; [u(α,τ)]=[u*(α)] eγ~τ, [v(α,τ)]=[v*(α)] eγ~τ, [w(α,τ)]=[w*(α)] eγ~τ, [θβ(α,τ)]=[θβ*(α)] eγ~τ are functions of displacement jump and rotation angles, respectively; τ=t/T is the dimensionless time; T=πR/vτ is the time for an elastic wave to traverse half of the shell circumference; vτ=E/ρ is the velocity of elastic wave propagation in the shell material; E is the Young’s modulus of the shell material; ρ is the density of the shell material; γ~ is constant coefficient that characterizes the rate of growtes of the load; δ( ) is the Dirac delta function; jn=n/jn, n=1,2,3,, j=α,β,τ; N20, S 0, Q20*, M20 – are the normal, shear, and generalized shear forces and the bending moment in a shell without a crack respectively; N21, S 1, Q21*, M21 are the similar efforts applied to the real crack edges; corresponding forces and moments on the edges of a notch β=± 0 are marked by the superscripts “+” and “-”.

In the case of a complex value of γ~=Reγ~+i Imγ~, we will come to the formulation of the problem for a load that varies harmonically according to the law of cosine or sine. In this formulation of the problem, it is necessary to consider the equality of their real (cosine law) or imaginary (sinusoidal law) parts in Eqs. (1-3).

The system of equilibrium equations in displacements that takes into account the inertia of the material in accordance with [2, 4, 9] is as follows:

5
Lk1u+Lk2v+Lk3w-R2cτ-2g¨k=qk0*, k=1,2,3,

where cτ2=vτ(1-ν2); ν is the Poisson’s ratio; we do not give the expressions for Lk1, Lk2, Lk3 because of their cumbersomeness (they can be found in [10]).

3. Construction of a problem solutions

In [9] is given detailed presentation of the methodology for solving the formulated problem and the procedure of constructing a system of singular integral equations, the solution of which makes it possible to calculate the dynamic coefficients of intensity of forces and moments near the tips of the crack. Since the above expressions are extremely cumbersome, we do not present them in this paper.

Fig. 1Сylindrical shell with a longitudinal crack

Сylindrical shell with a longitudinal crack

The intensity factors (IF) of forces and moments are given as follows [10]:

6
KN=12l KN0eγ~t, KM=12Rcl KM0eγ~t,

where:

7
KN0=limu11-u2Ω1u, KM0=limu11-u2Ω3u,

where the functions Ω1(u), Ω3(u) are solutions of the singular integral equations system [10]:

8
-11Ω1(u)u-sdu+α0-11Ω1(u)K011(z) du+α0-11Ω3(u)K013(z) du=f1*α0s,
-11Ω3(u)u-sdu+α0-11Ω1(u)K031(z) du+α0-11Ω3(u)K033(z) du=f3*α0s.

Also, [10] describes in detail all the steps of constructing these singular integral equations (SIE) system and provides expressions for all the necessary functions.

In Eq. (6) the factor eγ~t characterizes the variation of the value of the intensity factors of forces and moments depending on the time of load action, the parameter γ~ describes the rate of the load’s variation, and 12lKN0, 12RclKM0 take into account the geometric and physical parameters of the shell as well as the load’s variation rate.

Similarly to [2, 3] to estimate the intensity of forces and moments, we will use relative intensity factors for dynamic loading KN0 = 2KNle-γ~t and KM0=2KMRcle-γ~t. These factors describe, respectively, the intensity of efforts and moments in a period of time t<T.

4. Numerical results and discussion

Based on the constructed in [10] solutions of the SIE system Eq. (8) and on the relationships Eq. (7), the dependence of the relative IF of forces and moments on the rate of load variation and crack length was studied.

The study was carried out for a cylindrical shell with a longitudinal crack under the following assumptions:

1) The radius of the shell and its thickness are respectively equal to R = 15×10-2 m, h = 15×10-4 m;

2) To the sides of the crack are applied the following forces and moments with constant intensity N2(α,0,τ)=- N0eγ~τ=eγ~τ and M2(α,0,τ)=-M0eγ~τ=0.

3) The shell can be made of one of the following materials: steel (E= 2×1011 Pa, ν=0.25, ρ= 8×103 kg/m3), copper (E= 1.1×10^11 Pa, ν=0.34, ρ= 8.9×103 kg/m3), cast iron (E= 1×1011 Pa, ν=0.22, ρ= 7×103 kg/m3).

4) The approximate value of the time for which an elastic wave traverses the perimeter of such shell, depending on its material, is equal to Tst=0.9×10-4 for steel, Tcop=1.3×10-4 for copper and Tc-iron=1.2×10-4 for cast iron.

The research results are presented graphically in Figs. 2-4, where the dashed curves correspond to the values calculated for the shell made of cast iron, dotted curves - of copper, solid curves – of steel.

The relative dynamic IFs of forces (Fig. 2(a), Fig. 3(a)) and moments (Fig. 2(b), Fig. 3(b)) at low load growth rates γ~<0.01 are different from the corresponding relative IFs at static load [2, 3] no more than 1 %, which indirectly confirms the reliability of the performed studies. With an increase in the loading rate, the relative force intensity factor decreases approaching a certain limit value, which depends on the crack length. The critical value of the loading rate, after which the values of the dynamic IF begin to differ by more than 5 % from the values of the static IF, depends on the length of the crack (decreases with its prolongation).

Fig. 2Dependence of relative intensity factors of forces KN0 and moments KM0, on the relative crack length α0=l/R for values of the loading rate γ~

Dependence of relative intensity factors of forces KN0 and moments KM0,  on the relative crack length α0=l/R for values of the loading rate γ~

a) Relative IF of the efforts

Dependence of relative intensity factors of forces KN0 and moments KM0,  on the relative crack length α0=l/R for values of the loading rate γ~

b) Relative IF of the moments

Fig. 3Dependence of the relative intensity factors of forces KN0 and moments KM0 on the rate of load variation γ~ for different values of the relative crack length α0=l/R

Dependence of the relative intensity factors of forces KN0 and moments KM0  on the rate of load variation γ~ for different values of the relative crack length α0=l/R

а) Relative IF of the efforts

Dependence of the relative intensity factors of forces KN0 and moments KM0  on the rate of load variation γ~ for different values of the relative crack length α0=l/R

b) Relative IF of the moments

Based on the analysis of SIE kernels in [9], was made an assumption, that it is possible to neglect the influence IF of moments when estimating the strength of the shell and predicting its fracture. The results of calculations of the values of relative IF of forces (Fig. 2(a), Fig. 3(a)) and values of IF of moments (Fig. 2(b), Fig. 3(b)) confirm this assumption. The diagrams show that the moments IF are more than three orders of magnitude smaller than the corresponding forces IF, and therefore have little effect on the overall picture of the limit equilibrium

The value of the relative dynamic IF of moments (Fig. 2(b), Fig. 3(b)) increases with increasing loading rate. The graphs clearly show local maxima at the critical load rate, after which the IF values decrease towards a certain asymptote. Note that the critical rate is greater if the longer the crack length is, and the asymptotic value of the relative dynamic IF is higher than the value of the corresponding relative IF of moments under static loading (at low rates of loading γ~<0.01).

The respective diagrams show that the relative force and moment intensity factors for shells made of materials such as copper and cast iron are similar in their values. At the same time, they are twice as much as the corresponding values for the steel shell. A comparison of the mechanical characteristics of the respective materials allows us to conclude that the value of the Young’s modulus has a predominant influence on the magnitude of the relative dynamic IF of forces and moments.

An illustration of the influence of the relative dynamic factor on the overall pattern of changes in the values of the dynamic force factor KN=0.5lKN0eγ~τ is shown in Fig. 4. It shows the dependence on the relative time (0τ1) of the actual dynamic intensity factors of forces (Fig. 4(a)) and their exponential time dependence eγ~τ (Fig. 4(b)) in a steel shell with a longitudinal crack (crack length α0=3.85) at different loading rates γ~=0.5, γ~=1, γ~=2. The analysis of the presented graphs demonstrates that for an adequate qualitative and quantitative assessment of the value of dynamic intensity factor of forces and moments, it is necessary to take into account the value of the relative IF at the initial stage of loading.

Fig. 4Time dependence of dynamic intensity factors of forces and their exponential time dependence in a steel shell with a longitudinal crack

Time dependence of dynamic intensity factors of forces and their exponential  time dependence in a steel shell with a longitudinal crack

a) Dynamic intensity factors of effort

Time dependence of dynamic intensity factors of forces and their exponential  time dependence in a steel shell with a longitudinal crack

b) Exponential dependence

5. Conclusions

It is confirmed that the written in the explicit form system of singular integral equations with the Cauchy kernel [10] can be successfully used to determine the stress-strain state of cylindrical shells with a longitudinal crack under the action of a time-varying load. The reliability of the research is confirmed by the coincidence of the results of calculations at low rates of load variation with the known results for the case of static load [2, 3].

The effectiveness of the constructed algorithm for numerical solution of the system of equations obtained in [10] is confirmed.

The intensity factors of forces and moments near the tips of a longitudinal crack in the cylindrical shell under the action of a varying in time according to an exponential law load (this type of load occurs, for example, at a water hammer) are calculated. The corresponding intensity factors were calculated for shells made of different materials (copper, steel, cast iron).

The analysis of the obtained results is carried out, the main results of which are as follows:

– The Young’s modulus has a predominant influence on the magnitude of the dynamic intensity factor.

– At low rates of load variation, the effect of taking into account the inertia of the material when calculating the relative coefficients of force intensity can be neglected.

– The intensity factors of moments are more than three orders of magnitude smaller than the corresponding forces intensity factors, and therefore have a small effect on the overall picture of the limit equilibrium.

– With an increase in the rate of load variation, the IF of forces and moments approach to a certain asymptotic value.

The results of the studies confirm the effectiveness of the methods proposed in [9, 10] for determining the stress-strain state of cylindrical shells with longitudinal cuts under the action of a time-varying load. In the future, they can be used to determine the opening of the crack and estimate the strength of the shell.

References

  • H. S. Kit, R. M. Kushnir, V. V. Mykhas’Kiv, and M. M. Nykolyshyn, “Methods for the determination of static and dynamic stresses in bodies with subsurface cracks,” Materials Science, Vol. 47, No. 2, pp. 177–187, Nov. 2011, https://doi.org/10.1007/s11003-011-9382-9
  • R. M. Kushnir and M. M. Nykolyshyn, “Stress state and limit equilibrium of piecewise homogeneous cylindrical shells withcracks,” (in Ukrainian), Mathematical Methods and Physicomechanical Fields, Vol. 46, No. 1, pp. 60–74, 2003.
  • V. A. Osadchuk, Stress-Strain State and Limit Equilibrium of Shells with Notches. (in Russian), Kiev: Naukova Dumka, 1985.
  • Y. S. Podstrigach and R. N. Shvets, Thermoelasticity of Thin Shells. (in Russian), Kiev: Naukova Dumka, 1978.
  • I. D. Breslavsky, M. Amabili, and M. Legrand, “Static and dynamic behavior of circular cylindrical shell made of hyperelastic arterial material,” Journal of Applied Mechanics, Vol. 83, No. 5, p. 05100, May 2016, https://doi.org/10.1115/1.4032549
  • A. A. Hamzah, H. K. Jobair, O. I. Abdullah, E. T. Hashim, and L. A. Sabri, “An investigation of dynamic behavior of the cylindrical shells under thermal effect,” Case Studies in Thermal Engineering, Vol. 12, No. 5, pp. 537–545, Sep. 2018, https://doi.org/10.1016/j.csite.2018.07.007
  • S. M. Nabavi and A. R. Shahani, “Dynamic stress intensity factors for a longitudinal semi-elliptical crack in a thick-walled cylinder,” International Journal of Engineering, Science and Technology, Vol. 6, No. 5, pp. 57–77, Jan. 1970, https://doi.org/10.4314/ijest.v6i5.6
  • A. N. Guz, I. A. Guz, A. V. Men’Shikov, and V. A. Men’Shikov, “Stress-intensity factors for materials with interface cracks under harmonic loading,” International Applied Mechanics, Vol. 46, No. 10, pp. 1093–1100, Mar. 2011, https://doi.org/10.1007/s10778-011-0401-1
  • M. Makhorkin and M. Nykolyshyn, “Construction of integral equations describing limit equilibrium of cylindrical shell with a longitudinal crack under time-varying load,” Econtechmod, Vol. 5, No. 3, pp. 141–146, 2016.
  • M. I. Makhorkin and M. M. Nykolyshyn, “Limit equilibrium of a cylindrical shell with longitudinal crack with regard for the inertia of the material,” Journal of Mathematical Sciences, Vol. 249, No. 3, pp. 462–477, Aug. 2020, https://doi.org/10.1007/s10958-020-04953-4

About this article

Received
20 June 2024
Accepted
02 August 2024
Published
27 September 2024
SUBJECTS
Mathematical models in engineering
Keywords
cylindrical shell
longitudinal crack
time-varying load
exponential dependence on time
singular integral equations
dynamic intensity factors
Acknowledgements

The authors have not disclosed any funding.

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflict of interest

The authors declare that they have no conflict of interest.