The dynamic response study on thermo-magneto-elastic behavior of shallow conical shell in a time-dependent magnetic field is investigated, and the dynamic responses of displacement of shallow conical shell under mechanical loads, electromagnetic fields and temperature field coupling are analyzed. Based on Maxwell’s equations, heat conduction equation and nonlinear equations of classical plates and shells, the nonlinear dynamic response governing equations are derived. The electromagnetic field and temperature field equations are solved using variable separating technique, the nonlinear elastic field equations are solved by Galerkin method. The variation of temperature, magnetic field intensity and displacement with time under the coupling effect of the applied magnetic field and the surface uniform load were obtained. The influence of frequency of the applied magnetic field on the displacement wave forms is discussed.
Thermo-magneto-elasticity is a new subject to study the strength, stiffness and stability of elastic components under the combined action of electromagnetic, temperature and deformation. Magneto-elasticity is the theory of studying the coupling between electromagnetic and deformation, thermo-elasticity combines elasticity and heat conduction to study the coupling theory between temperature field and elastic field. The thermo-magneto-elasticity includes heat conduction theory, classical elasticity theory and electromagnetic theory. These theories are applied to solve the coupling problems of temperature field, electromagnetic field and elastic field of conductive elastic elements located in magnetic field and considering thermal effect.
With the wide application of electromagnetic structure in the high-tech field, many components work in the environment of temperature change in engineering, the research on the thermos-elastic phenomenon of electromagnetic coupling has a strong engineering background and theoretical value [
The key to the thermo-magneto-elastic analysis is the solution of the thermo-magneto-elastic equation. Xing studied on the dynamic response and quasi-static response of the rectangular plate thermo-magneto-elastic in the temperature-varying magnetic field with difference method. The curves of temperature, magnetic field intensity, induced current, and thermal elastic stress and displacement with time are obtained [
It is very difficult to study the nonlinear dynamic response of a shallow conical shell in an alternating magnetic field and subjected to mechanical loads. The study on the nonlinear dynamic response of the shell is very rare. In this paper, the electromagnetic field equations are derived based on Maxwell’s equations and Ohm’s law. Based on the heat conduction equation and the heat balance equation, temperature field equations are derived. Based on nonlinear equations of classical plates and shells, considering the coupling effect of the Lorentz force and temperature stress, nonlinear magneto-elastic heat equations of shallow conical shell are deduced. Applying Galerkin method, the solution of thermomagneto-elastic equation, the rule of temperature, magnetic field and displacement varying with time under the coupling effect of the applied magnetic field and surface uniform stable mechanical loads are obtained.
Considering shallow conical shell with thickness h, radium a and pyramid dip φ , whose neuter plane is showed in
In the absence of lateral current and the influence of displacement current and volume charge density is not considered, according to Maxwell equation and generalized Ohm’s law, the electrodynamics equations are:
r o t E = − ∂ B ∂ t r o t H = J J = σ ( E − V B ) B = μ H (1)
where E is electric field intensity, B is magnetic induction intensity, H is magnetic field intensity, J is current density, V is velocity, σ is admittance, μ is permeability, rot is rotation.
Ignoring the mechanical electric effect and considering the axial symmetry, the electrodynamics Equation (1) of shallow conical shell can be simplified as
∂ H r ∂ t = ϕ 2 μ σ ∂ 2 H r ∂ z 2 J θ = ϕ ∂ H r ∂ z (2)
As electromagnetic field varying with time induces current in shallow conical shell, which formats Joule heating effect, that is induction current loss. As to shallow conical shell, it can approximately assume that current distribute uniformly in shallow conical shell because of a low frequency of the applied magnetic field and current, then the current loss per unit time per volume can be calculated by the following formula:
Q = J θ 2 σ (3)
Set the initial temperature of shallow conical shell to zero, it is heated by Joule heating effect from t = 0, heat exchange exist among the Inner and outer surface and the bottom of the shell and media whose external temperature is zero. According to Fourier heat transfer law and energy conservation law, the control equation of heat conduction is established, which means that the transient temperature field T ( z , t ) of shallow conical shell should satisfy the followed equation:
∂ T ∂ t = k ∇ 2 T + Q ρ c (4)
where ρ is material mass density, c is specific heat capacity, k is coefficient of thermal conductivity, ∇ 2 is Laplace operator.
In the axial symmetry condition, the control equation of heat conduction of shallow conical shell can be simplified to
∂ T ∂ t = k φ 2 ∂ 2 T ∂ z 2 + Q ρ c (5)
Based on heat exchange law and the condition that the external temperature of media is zero, the heat balance equation of the current-carrying shell’s internal and external surface can be established. Thus, the boundary condition is
± ∂ T ∂ z + Λ T = 0 , at z = ± h 2 (6)
where, Λ is thermal coefficient.
The shallow conical shell in the time dependent electromagnetic field also suffer the temperature stress induced by Joule heat and Lorentz force in addition to the external mechanical loads P, The Lorentz force can be expressed as
f z = − μ H r ∂ H r ∂ z (7)
where f z is z direction component of Lorentz force f.
When discussing the thermo-magneto-elastic nonlinear problem of a shallow conical shell, the Kichhoff-love straight normal hypothesis is adopted, that the normal section perpendicular to the neutral surface before deformation of the shell remains a straight line after deformation and perpendicular to the neutral plane after deformation, and its length remains unchanged. Based on this assumption, the radial, toroidal and normal displacements of any point in a shell with a distance of z from the neutral plane under axial symmetry conditions can be expressed as:
u ( z ) = u − z ∂ w ∂ r v ( z ) = 0 w ( z ) = w (8)
Among them, u and w are the radial displacement and deflection at any point on the neutral surface of the shell respectively.
Considering the Von-Karman type large deflection geometric relation of conical shell, the strain-displacement relation at any point in the shell is obtained:
ε r ( z ) = ∂ u ∂ r + φ ∂ w ∂ r + 1 2 ( ∂ w ∂ r ) 2 − z ∂ 2 w ∂ r 2 ε θ ( z ) = u r − z 1 r ∂ w ∂ r (9)
The physical equation of shallow conical shell considering temperature change can be calculated according to follow formula.
σ r = E 1 − v 2 [ ( ε r + v ε θ ) − ( 1 + v ) α T ] σ θ = E 1 − v 2 [ ( ε θ + v ε r ) − ( 1 + v ) α T ] (10)
The E , v , α in the formula represents the elastic modulus, Poisson’s ratio and linear thermal expansion coefficient respectively, σ r , σ θ are r, θ stress respectively.
In the axial symmetry condition, considering the coupling effect of Lorentz force, temperature stress and mechanical load, according to the classical theory of plates and shells, the control equation of shallow conical shell can be derived as follows:
r 2 ∂ 2 u ∂ r 2 + r ∂ u ∂ r − u + ( 1 − ν ) r φ ∂ w ∂ r + ( 1 − ν ) 2 r ( ∂ w ∂ r ) 2 + r 2 ( φ + ∂ w ∂ r ) ∂ 2 w ∂ r 2 = 0 D ( r ∂ 3 w ∂ r 3 + ∂ 2 w ∂ r 2 − 1 r ∂ w ∂ r ) − r { A [ ∂ u ∂ r + φ ∂ w ∂ r + 1 2 ( ∂ w ∂ r ) 2 + ν u r ] − N T } ( φ + ∂ w ∂ r ) − 1 2 ( P z + F z ) r 2 = ρ h ∫ 0 r r ∂ 2 w ∂ t 2 d r (11)
where A and D are extension rigidity and bending rigidity of shallow conical shell respectively. F z and N T are the resultant forces of Lorentz force in the transverse direction and the film force produced by thermal stress respectively.
F z = ∫ − h / 2 h / 2 f z d z , N T = E α 1 − v ∫ − h / 2 h / 2 T ( z ) d z (12)
To facilitate the calculation, the following dimensionless parameters are introduced
r ¯ = r a , H ¯ r = H r H 0 , τ = 4 φ 2 t μ σ h 2 , u ¯ = 4 a h 2 u , w ¯ = 2 w h , k = 2 a h φ , z ¯ = 2 z h , J ¯ θ = h J θ 2 H 0 , Q ¯ = h 2 σ Q 4 H 0 2 φ 2 , T ¯ = c ρ T μ H 0 2 , h ¯ = h Λ / 2 , ( F ¯ z , P ¯ z ) = a 4 D h 2 ( F z , P z ) , P ¯ z = 8 a 4 A h 3 P z , N ¯ T = a 2 D N T , f ¯ z = h 2 μ H 0 2 f z , ( σ ¯ r , σ ¯ θ ) = 4 a 2 ( 1 − v 2 ) h 2 E ( σ r , σ θ ) , ς 1 = k σ μ , ς 2 = A h 2 4 D , ς 3 = 32 ρ a 4 φ 4 D σ 4 μ 2 h (13)
The non-dimensional forms of Equation (2), (3), (5) and (11) are given as follows
∂ H ¯ r ∂ τ = ∂ 2 H ¯ r ∂ z ¯ 2 (14)
Q ¯ = J ¯ θ 2 (15)
∂ T ¯ ∂ τ = ς 1 ∂ 2 T ¯ ∂ z ¯ 2 + Q ¯ (16)
r ¯ 2 ∂ 2 u ¯ ∂ r ¯ 2 + r ¯ ∂ u ¯ ∂ r ¯ − u ¯ + ( 1 − ν ) k r ¯ ∂ w ¯ ∂ r ¯ + 1 2 ( 1 − ν ) r ¯ ( ∂ w ¯ ∂ r ¯ ) 2 + r ¯ 2 ( k + ∂ w ¯ ∂ r ¯ ) ∂ 2 w ¯ ∂ r ¯ 2 = 0 (17a)
r ¯ ∂ 3 w ¯ ∂ r ¯ 3 + ∂ 2 w ¯ ∂ r ¯ 2 − 1 r ¯ ∂ w ¯ ∂ r ¯ − ς 2 r ¯ { [ ∂ u ¯ ∂ r ¯ + k ∂ w ¯ ∂ r ¯ + 1 2 ( ∂ w ¯ ∂ r ¯ ) 2 + ν u ¯ r ¯ ] − N T } ( k + ∂ w ¯ ∂ r ¯ ) − ( P ¯ z + F ¯ z ) r ¯ 2 = ς 3 ∫ 0 r ¯ r ¯ ∂ 2 w ¯ ∂ τ 2 d r ¯ (17b)
In Equation (17), the expressions of F ¯ z and N T are respectively as follow:
F ¯ z = ∫ − 1 1 f ¯ z d z ¯ , N T = ∫ − 1 1 T ¯ ( z ¯ , τ ) d z ¯ (18)
Consider the following boundary conditions and initial conditions of electromagnetic field, temperature field and elastic field
H ¯ r = φ ( τ ) at z ¯ = 1 H ¯ r = 0 at z ¯ = − 1 H ¯ r = 0 when τ = 0 (19)
± ∂ T ¯ ∂ z ¯ + h ¯ T ¯ = 0 at z ¯ = ± 1 T ¯ ( z ¯ , τ ) = 0 when τ = 0 (20)
w ¯ = 0 , ∂ w ¯ ∂ r ¯ = 0 , u ¯ = 0 at r ¯ = 1 (21)
∂ w ¯ ∂ r ¯ = 0 , u ¯ is finite at r ¯ = 0 (22)
w ¯ = ∂ w ¯ ∂ τ = 0 and u ¯ = ∂ u ¯ ∂ τ = 0 when τ = 0 (23)
To solve the governing Equation (17), the first step is to solve the electromagnetic field Equation (14) to obtain the magnetic field strength H ¯ r and Lorentz force f ¯ z . Solving Equation (16) gets the temperature field distribution T ¯ , and obtains F ¯ z and N T according to formula (18).
Equation (14) is solved by separation of variables. In order to make non-homogeneous boundary Condition (19) homogeneous, new unknown function h r ( z ¯ , τ ) is introduced
h r ( z ¯ , τ ) = H ¯ ( z ¯ , τ ) − z ¯ + 1 2 ϕ ( τ ) (24)
Substituting Formula (24) in to Equation (14), boundary conditions and initial Condition (19) can be written as follows
∂ h r ∂ z ¯ 2 = ∂ h r ∂ τ + 1 + z ¯ 2 d ϕ ( τ ) d τ (25)
h r = 0 at z ¯ = ± 1 (26)
h r ( z ¯ , 0 ) = − ϕ ( 0 ) 1 + z ¯ 2 when τ = 0 (27)
Assume that the solution of Equation (25) satisfy the Conditions (26) and (27) is following series form
h r = ∑ n = 1 ∞ ψ n ( τ ) cos ( k n z ¯ ) (28)
where ψ ( τ ) is the function of τ , k n is the position root of function cos ( k n ) = 0 ,that is
k n = 2 n − 1 2 π ( n = 1 , 2 , 3 , ⋯ ) (29)
Substituting Formula (28) in the Equation (25), both sides of Equation (25) multiply cos ( k m z ¯ ) , using the orthogonality of trigonometric function, it obtains
d ψ ( τ ) d τ + k n 2 ψ ( τ ) + ∫ − 1 1 d ϕ ( τ ) d τ 1 + z ¯ 2 cos ( k n z ¯ ) d z ¯ = 0 (30)
With initial Condition (22), the solution of Equation (30) is
ψ n = ( − 1 ) n k n ψ ¯ n ( τ ) (31)
where ψ ¯ ( τ ) = ∫ 0 τ e − k n 2 ( τ − x ) d φ ( x ) d x d x
Thus, the expressions of non-dimensional magnetic field intensity H ¯ r , induced current and Lorentz force respectively are
H ¯ r = ∑ n = 1 ∞ ( − 1 ) n k n ψ ¯ ( τ ) cos ( k n z ¯ ) + 1 + z ¯ 2 φ ( τ ) J ¯ θ = ∑ n = 1 ∞ ( − 1 ) n + 1 sin ( k n z ¯ ) ψ ¯ ( τ ) + 1 2 φ ( τ ) f ¯ z = [ ∑ m = 1 ∞ ( − 1 ) m ψ ¯ m ( τ ) sin ( k m z ¯ ) − 1 2 φ ( τ ) ] [ ∑ n = 1 ∞ ( − 1 ) n k n ψ ¯ ( τ ) cos ( k n z ¯ ) + 1 + z ¯ 2 φ ( τ ) ] (32)
Based on the separation of variables, assume that the solution of Equation (16) which satisfies the boundary Condition (20) shows the following form
T ¯ = ∑ i = 1 ∞ b i ( τ ) cos ( a i z ¯ ) (33)
where a i is confirmed by tan a i = h ¯ a i . Both sides of Equation (12) multiply
cos ( a j z ¯ ) , using the orthogonality of trigonometric function, it obtains
∂ b i ( τ ) ∂ τ + ς 1 a i 2 b i ( τ ) − h ¯ 2 + a i 2 h ¯ + h ¯ 2 + a i 2 ∫ − 1 1 Q ¯ ( z ¯ , τ ) cos ( a i z ¯ ) d z ¯ = 0 (34)
The solution of Equation (34) is
b i ( τ ) = h ¯ 2 + a i 2 h ¯ + h ¯ 2 + a i 2 [ ∑ m = 1 ∞ ∑ n = 1 ∞ L i m n b ^ i m n ( τ ) + c ( τ ) ] (35)
where
L i m n = h ¯ cos ( a i ) 2 ( k m 2 + k n 2 − a i 2 ) ( 2 k m k n ) 2 − ( k m 2 + k n 2 − a i 2 ) 2 b ^ i m n = ∫ 0 τ e − ς 1 a i 2 ( τ − ξ ) ψ ¯ n ( ξ ) ψ ¯ m ( ξ ) d ξ c ( τ ) = sin ( a i ) 4 a i 4 ω ( 1 − e − ς 1 a i 2 τ ) + a i 4 ς 1 2 ( 1 − cos ( 2 ω τ ) ) − 2 ω a i 2 ς 1 sin ( 2 ω τ ) 4 ω 2 ς 1 a i 2 + a i 6 ς 1 3
Substituting expression (35) into expression (33), it obtains
T ¯ ( z ¯ , τ ) = ∑ i = 1 ∞ h ¯ 2 + a i 2 h ¯ + h ¯ + a i 2 [ ∑ m = 1 ∞ ∑ n = 1 ∞ L i m n b ^ i m n + c ( τ ) ] cos ( a i z ¯ ) (36)
After obtaining Lorentz force f ¯ z and temperature field distribution T ¯ , F ¯ z and N T are obtained according to formula (18), and then substitute it in Equation (17). Applying Galerkin method, the initial boundary value problem (17a), (17b) and (21)-(23). Assume that the deflection w ¯ ( r ¯ , τ ) satisfies with the boundary condition (21) and (22), and has a separable form of time and space as follows
w ¯ = f ( τ ) ( 1 − 2 r ¯ 2 + r ¯ 4 ) (37)
Substituting expression (37) into Equation (17a), the solution of equation (17a) satisfying with boundary Condition (21) and (22) is
u ¯ ( r ¯ , τ ) = k f ( τ ) ( 24 − 16 v 15 r ¯ + 4 v − 2 3 r ¯ 2 + 1 − 4 v 15 r ¯ 4 ) + f 2 ( τ ) [ 6 − 3 v 6 r ¯ + ( 1 − 3 v ) r ¯ 3 − 8 − 2 v 3 r ¯ 5 + 6 − v 6 r ¯ 7 ] (38)
Applying Galerkin integral to Equation (17b), it obtains
∫ 0 1 ∫ 0 2π [ ( r ¯ ∂ 3 w ¯ ∂ r ¯ 3 + ∂ 2 w ¯ ∂ r ¯ 2 − 1 r ¯ ∂ w ¯ ∂ r ¯ ) − ς 2 r ¯ { [ ∂ u ¯ ∂ r ¯ + k ∂ w ¯ ∂ r ¯ + 1 2 ( ∂ w ¯ ∂ r ¯ ) 2 + ν u ¯ r ¯ ] − N T } ( k + ∂ w ¯ ∂ r ¯ ) − ( P ¯ z + F ¯ z ) r ¯ 2 − ς 3 ∫ 0 1 r ¯ ∂ 2 w ¯ ∂ τ 2 d r ¯ ] ∂ w ¯ ∂ f r ¯ d θ d r ¯ = 0 (39)
The dynamic response formula expressed by the deflection of the shell can be obtained
a 1 d 2 f d τ 2 + a 2 f 3 + a 3 f 2 + a 4 f − p 0 − p 1 − p 2 = 0 (40)
In which
a 1 = ς 3 ∫ 0 1 r ¯ ( ∫ 0 1 r ¯ G 1 d r ¯ ) G 1 d r ¯ , a 2 = ς 2 ∫ 0 1 [ G ′ 3 + 1 2 G ′ 1 2 + ν r ¯ ] G ′ 1 G 1 r ¯ 2 d r ¯ , a 3 = ς 2 ∫ 0 1 r ¯ 2 { k [ v r G 3 + G ′ 3 + 1 2 v G ′ 1 2 ] + [ G ′ 2 G ′ 1 + k G ′ 1 2 + v r G 2 G ′ 1 ] } G 1 d r ¯ , a 4 = ∫ 0 1 { − [ r ¯ G ‴ 1 + G ″ 1 − 1 r ¯ G ′ 1 ] + ς 2 r ¯ k [ G ′ 4 + k G ′ 1 + v r ¯ G 2 ] } G 1 r ¯ d r ¯ , p 0 = P ¯ ∫ 0 1 r ¯ 3 G 1 d r ¯ z , p 1 = − ς 2 N ¯ T ∫ 0 1 ( k + G ′ ) G 1 r ¯ 2 d r ¯ , p 2 = F ¯ z ∫ 0 1 r ¯ 3 G 1 d r ¯
where
G 1 ( r ¯ ) = 1 − 2 r ¯ 2 + r ¯ 4 , G 2 ( r ¯ ) = 24 − 16 v 15 r ¯ + 4 v − 2 3 r ¯ 2 + 1 − 4 v 15 r ¯ 4 G 3 ( r ¯ ) = 6 − 3 v 6 r ¯ + ( 1 − 3 v ) r ¯ 3 − 8 − 2 v 3 r ¯ 5 + 6 − v 6 r ¯ 7
Numerical solution of the dynamic response Equation (40) is obtained by Runge-Kutta method, the response relation of deflection and time is obtained.
The shallow conical thin shell is shown in
ϕ ( τ ) = sin ( ω τ ) (41)
where ω is the non-dimensional angular frequency of magnetic field. The physical parameters of the shallow conical thin shell are as follows
μ = 4 π × 10 − 7 H / m , σ = 3.42 × 10 7 s / m , c = 2.7 × 10 3 J / kg ,
ρ = 0.9 × 10 3 kg / m 3 , k = 92.6 × 10 − 6 m 2 / s , ν = 0.3 , E = 70 GPa ,
α = 24 × 10 − 6 k − 1 , h = 2 mm , a = 0.2 m , Λ = 204 1 / m
According to the above analysis, numerical calculation is obtained by Matlab, the results is shown in Figures 2-7.
induced current J ¯ θ in different thickness varying with time in the different thickness. It can be seen from
decreases gradually from surface to central face.
Based on Maxwell’s equations, heat conduction equation and nonlinear equations of classical plates and shells, the dynamic response study on shallow conical shell’s thermo-magneto-elastic behavior in a time-dependent magnetic field is presented. Some conclusions can be obtained through the calculation and analysis of shallow conical shell instances:
1) In the condition that other parameters are constant, mechanical load has an effect on the displacement amplitude of shallow conical shell, but has no effect on the vibration frequency.
2) In the condition that other physical parameters are invariable, the strength of applied magnetic field has influence to the displacement, but has no effect on the vibration frequency.
3) The stress and strain of plate and shell can be controlled when the parameters of magnetic field and mechanical load change appropriately. It has a certain reference value to the practical application of magneto-elastic coupling theory.
The authors declare no conflicts of interest regarding the publication of this paper.
Cen, S.D. and Xu, J.C. (2018) Analysis of Thermo-Magneto-Elastic Nonlinear Dynamic Response of Shallow Conical Shells. Engineering, 10, 837-850. https://doi.org/10.4236/eng.2018.1012058
1. a: radius of the shell
2. Thickness of the shell
3. φ : the cone Angle
4. E θ : θ direction of electric field intensity
5. B r : r direction of magnetic induction intensity
6. H r : r direction of magnetic field intensity
7. J: current density
8. σ: electrical conductivity
9. μ: magnetic permeability
10. ρ: material mass density
11. C: the specific heat capacity of the material
12. K: thermal conductivity
13. Q: current loss per unit time per volume
14. f: Lorentzforce
15. u: radial displacement of the shell in neutral
16. w: deflection of the shell in neutral
17. E: Elastic Modulus
18. v: Poissonratio
19. α : Coefficient of linear thermal expansion
20. σ r : r direction of the stress
21. σ θ : θ direction of the stress
22. A: Tensile rigidity
23. D: Flexural rigidity