**Engineering**

Vol.10 No.12(2018), Article ID:89137,14 pages

10.4236/eng.2018.1012058

Analysis of Thermo-Magneto-Elastic Nonlinear Dynamic Response of Shallow Conical Shells

Shende Cen, Jiachu Xu^{ }

School of Mechanics and Construction Engineering, Jinan University, Guangzhou, China

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: September 13, 2018; Accepted: December 11, 2018; Published: December 14, 2018

ABSTRACT

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.

**Keywords:**

Thermo-Magneto-Elastic, Shallow Conical Shell, Nonlinear Dynamic Response, Galerkin Method

1. Introduction

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 [1] [2] . The elastic elements and structures in high energy-varying magnetic field under mechanical loads can produce various stresses. In addition to mechanical stress, there are the thermal stress generated by the induced eddy current losses, and the magnetic stress generated by the Lorentz force. These stresses affect each other, and to be high nonlinear. Previous studies on the thermo-magneto-elastic problem of plate and shell were mainly based on the simplified theory and the linear theories. However, in actual situation, most plate and shell structure are in the temperature-varying high-energy electromagnetic fields, which is a highlycoupled nonlinearity.

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 [3] . Higuchi et al. studied the cylinder affected by the change magnetic field and the thermo-magneto-elastic stress of the hollow cylindrical shell with separation of variables [4] [5] . Wang Ping studied the chaotic motion of the large deflection simple support plate under the coupling effect of mechanical loads, electromagnetic and temperature field [6] [7] . Alberto Milazzo presented a new one-dimensional model for the dynamic problem of magneto-electro-elastic generally laminated beams [8] . Based on the successful validation of the model, new results for free vibrations of functionally graded magneto-electro-elastic beams are presented. Zhang Lang analyzed the buckling and vibration of functionally graded magneto-electro-thermo-elastic circular cylindrical shells. Kattimani studied the Control of geometrically nonlinear vibrations of functionally graded magneto-electro-elastic plates [9] [10] .

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.

2. Basic Equations

Considering shallow conical shell with thickness h, radium a and pyramid dip $\phi $ , whose neuter plane is showed in Figure 1 Assume that shallow conical shell in alternating magnetic field works under axisymmetric state, whose outer surface is subjected to normal stable mechanical load. The orthogonal curvilinear coordinate $\left(r,\theta ,z\right)$ is established in Figure 1, where r, θ, z are radium, annular and normal coordinate of shallow conical shell respectively. The applied magnetic field intensity is $H\left({H}_{r},0,0\right)$ and Mechanical load is $P\left(0,0,{P}_{z}\right)$ . ${H}_{r}$ is a function of the coordinate z and time t, and ${P}_{z}$ is constant.

2.1. Electro Dynamics Equations

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:

$\begin{array}{l}rotE=-\frac{\partial B}{\partial t}\\ rotH=J\\ J=\sigma \left(E-VB\right)\\ B=\mu H\end{array}$ (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

$\begin{array}{l}\frac{\partial {H}_{r}}{\partial t}=\frac{{\varphi}^{2}}{\mu \sigma}\frac{{\partial}^{2}{H}_{r}}{\partial {z}^{2}}\\ {J}_{\theta}=\varphi \frac{\partial {H}_{r}}{\partial z}\end{array}$ (2)

Figure 1. The diagram of shallow conical shell.

2.2. Temperature Field

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=\frac{{J}_{\theta}^{2}}{\sigma}$ (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\left(z,t\right)$ of shallow conical shell should satisfy the followed equation:

$\frac{\partial T}{\partial t}=k{\nabla}^{2}T+\frac{Q}{\rho c}$ (4)

where ρ is material mass density, c is specific heat capacity, k is coefficient of thermal conductivity, ${\nabla}^{2}$ is Laplace operator.

In the axial symmetry condition, the control equation of heat conduction of shallow conical shell can be simplified to

$\frac{\partial T}{\partial t}=k{\phi}^{2}\frac{{\partial}^{2}T}{\partial {z}^{2}}+\frac{Q}{\rho 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

$\pm \frac{\partial T}{\partial z}+\Lambda T=0$ , at $z=\pm \frac{h}{2}$ (6)

where, $\Lambda $ is thermal coefficient.

2.3. Elastic Field

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}=-\mu {H}_{r}\frac{\partial {H}_{r}}{\partial 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:

$\begin{array}{l}u\left(z\right)=u-z\frac{\partial w}{\partial r}\\ v\left(z\right)=0\\ w\left(z\right)=w\end{array}$ (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:

$\begin{array}{l}{\epsilon}_{r}\left(z\right)=\frac{\partial u}{\partial r}+\phi \frac{\partial w}{\partial r}+\frac{1}{2}{\left(\frac{\partial w}{\partial r}\right)}^{2}-z\frac{{\partial}^{2}w}{\partial {r}^{2}}\\ {\epsilon}_{\theta}\left(z\right)=\frac{u}{r}-z\frac{1}{r}\frac{\partial w}{\partial r}\end{array}$ (9)

The physical equation of shallow conical shell considering temperature change can be calculated according to follow formula.

$\begin{array}{l}{\sigma}_{r}=\frac{E}{1-{v}^{2}}\left[\left({\epsilon}_{r}+v{\epsilon}_{\theta}\right)-\left(1+v\right)\alpha T\right]\\ {\sigma}_{\theta}=\frac{E}{1-{v}^{2}}\left[\left({\epsilon}_{\theta}+v{\epsilon}_{r}\right)-\left(1+v\right)\alpha T\right]\end{array}$ (10)

The $E,v,\alpha $ in the formula represents the elastic modulus, Poisson’s ratio and linear thermal expansion coefficient respectively, ${\sigma}_{r},{\sigma}_{\theta}$ 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:

$\begin{array}{l}{r}^{2}\frac{{\partial}^{2}u}{\partial {r}^{2}}+r\frac{\partial u}{\partial r}-u+\left(1-\nu \right)r\phi \frac{\partial w}{\partial r}+\frac{\left(1-\nu \right)}{2}r{\left(\frac{\partial w}{\partial r}\right)}^{2}+{r}^{2}\left(\phi +\frac{\partial w}{\partial r}\right)\frac{{\partial}^{2}w}{\partial {r}^{2}}=0\\ D\left(r\frac{{\partial}^{3}w}{\partial {r}^{3}}+\frac{{\partial}^{2}w}{\partial {r}^{2}}-\frac{1}{r}\frac{\partial w}{\partial r}\right)-r\left\{A\left[\frac{\partial u}{\partial r}+\phi \frac{\partial w}{\partial r}+\frac{1}{2}{\left(\frac{\partial w}{\partial r}\right)}^{2}+\nu \frac{u}{r}\right]-{N}^{T}\right\}\left(\phi +\frac{\partial w}{\partial r}\right)\\ -\frac{1}{2}\left({P}_{z}+{F}_{z}\right){r}^{2}=\rho h{\displaystyle {\int}_{\text{}0}^{\text{}r}r\frac{{\partial}^{2}w}{\partial {t}^{2}}\text{d}r}\end{array}$ (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}={\displaystyle {\int}_{-h/2}^{h/2}{f}_{z}\text{d}z},\text{}{N}^{T}=\frac{E\alpha}{1-v}{\displaystyle {\int}_{-h/2}^{h/2}T\left(z\right)\text{d}z}$ (12)

2.4. Non-Dimensionalization

To facilitate the calculation, the following dimensionless parameters are introduced

$\begin{array}{l}\stackrel{\xaf}{r}=\frac{r}{a},\text{\hspace{0.17em}}{\stackrel{\xaf}{H}}_{r}=\frac{{H}_{r}}{{H}_{0}},\text{\hspace{0.17em}}\tau =\frac{4{\phi}^{2}t}{\mu \sigma {h}^{2}},\text{\hspace{0.17em}}\text{}\stackrel{\xaf}{u}=\frac{4a}{{h}^{2}}u,\text{\hspace{0.17em}}\stackrel{\xaf}{w}=\frac{2w}{h},\text{\hspace{0.17em}}\text{\hspace{0.17em}}k=\frac{2a}{h}\phi ,\text{}\stackrel{\xaf}{z}=\frac{2z}{h},\\ {\stackrel{\xaf}{J}}_{\theta}=\frac{h{J}_{\theta}}{2{H}_{0}},\text{\hspace{0.17em}}\stackrel{\xaf}{Q}=\frac{{h}^{2}\sigma Q}{4{H}_{0}{}^{2}{\phi}^{2}},\text{\hspace{0.17em}}\stackrel{\xaf}{T}=\frac{c\rho T}{\mu {H}_{0}{}^{2}},\text{\hspace{0.17em}}\stackrel{\xaf}{h}=h\Lambda /2,\text{\hspace{0.17em}}\left({\stackrel{\xaf}{F}}_{z},{\stackrel{\xaf}{P}}_{z}\right)=\frac{{a}^{4}}{D{h}^{2}}\left({F}_{z},{P}_{z}\right),\\ {\stackrel{\xaf}{P}}_{z}=\frac{8{a}^{4}}{A{h}^{3}}{P}_{z},\text{\hspace{0.17em}}{\stackrel{\xaf}{N}}^{T}=\frac{{a}^{2}}{D}{N}^{T},\text{\hspace{0.17em}}{\stackrel{\xaf}{f}}_{z}=\frac{h}{2\mu {H}_{0}^{2}}{f}_{z},\text{\hspace{0.17em}}\left({\stackrel{\xaf}{\sigma}}_{r},{\stackrel{\xaf}{\sigma}}_{\theta}\right)=\frac{4{a}^{2}\left(1-{v}^{2}\right)}{{h}^{2}E}\left({\sigma}_{r},{\sigma}_{\theta}\right),\\ {\varsigma}_{1}=k\sigma \mu ,\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\varsigma}_{2}=\frac{A{h}^{2}}{4D},\text{\hspace{0.17em}}{\varsigma}_{3}=\frac{32\rho {a}^{4}{\phi}^{4}}{D{\sigma}^{4}{\mu}^{2}h}\end{array}$ (13)

The non-dimensional forms of Equation (2), (3), (5) and (11) are given as follows

$\frac{\partial {\stackrel{\xaf}{H}}_{r}}{\partial \tau}=\frac{{\partial}^{2}{\stackrel{\xaf}{H}}_{r}}{\partial {\stackrel{\xaf}{z}}^{2}}$ (14)

$\stackrel{\xaf}{Q}={\stackrel{\xaf}{J}}_{\theta}^{2}$ (15)

$\frac{\partial \stackrel{\xaf}{T}}{\partial \tau}={\varsigma}_{1}\frac{{\partial}^{2}\stackrel{\xaf}{T}}{\partial {\stackrel{\xaf}{z}}^{2}}+\stackrel{\xaf}{Q}$ (16)

${\overline{r}}^{2}\frac{{\partial}^{2}\overline{u}}{\partial {\overline{r}}^{2}}+\overline{r}\frac{\partial \overline{u}}{\partial \overline{r}}-\overline{u}+\left(1-\nu \right)k\overline{r}\frac{\partial \overline{w}}{\partial \overline{r}}+\frac{1}{2}\left(1-\nu \right)\overline{r}{\left(\frac{\partial \overline{w}}{\partial \overline{r}}\right)}^{2}+{\overline{r}}^{2}\left(k+\frac{\partial \overline{w}}{\partial \overline{r}}\right)\frac{{\partial}^{2}\overline{w}}{\partial {\overline{r}}^{2}}=0$ (17a)

$\begin{array}{l}\overline{r}\frac{{\partial}^{3}\overline{w}}{\partial {\overline{r}}^{3}}+\frac{{\partial}^{2}\overline{w}}{\partial {\overline{r}}^{2}}-\frac{1}{\overline{r}}\frac{\partial \overline{w}}{\partial \overline{r}}-{\varsigma}_{2}\overline{r}\left\{\left[\frac{\partial \overline{u}}{\partial \overline{r}}+k\frac{\partial \overline{w}}{\partial \overline{r}}+\frac{1}{2}{\left(\frac{\partial \overline{w}}{\partial \overline{r}}\right)}^{2}+\nu \frac{\overline{u}}{\overline{r}}\right]-{N}^{T}\right\}\left(k+\frac{\partial \overline{w}}{\partial \overline{r}}\right)\\ -\left({\overline{P}}_{z}+{\overline{F}}_{z}\right){\overline{r}}^{2}={\varsigma}_{3}{\displaystyle {\int}_{0}^{\overline{r}}\overline{r}\frac{{\partial}^{2}\overline{w}}{\partial {\tau}^{2}}d\overline{r}}\end{array}$ (17b)

In Equation (17), the expressions of ${\stackrel{\xaf}{F}}_{z}$ and ${N}^{T}$ are respectively as follow:

${\stackrel{\xaf}{F}}_{z}={\displaystyle {\int}_{-1}^{1}{\stackrel{\xaf}{f}}_{z}\text{d}\stackrel{\xaf}{z}},\text{}{N}^{T}={\displaystyle {\int}_{-1}^{1}\stackrel{\xaf}{T}\left(\stackrel{\xaf}{z},\tau \right)\text{d}\stackrel{\xaf}{z}}$ (18)

Consider the following boundary conditions and initial conditions of electromagnetic field, temperature field and elastic field

$\begin{array}{l}{\stackrel{\xaf}{H}}_{r}=\phi \left(\tau \right)\text{at}\text{\hspace{0.17em}}\stackrel{\xaf}{z}=1\\ {\stackrel{\xaf}{H}}_{r}=0\text{at}\text{\hspace{0.17em}}\stackrel{\xaf}{z}=-1\\ {\stackrel{\xaf}{H}}_{r}=0\text{when}\text{\hspace{0.17em}}\tau =0\end{array}$ (19)

$\begin{array}{l}\pm \frac{\partial \stackrel{\xaf}{T}}{\partial \stackrel{\xaf}{z}}+\stackrel{\xaf}{h}\stackrel{\xaf}{T}=0\text{at}\stackrel{\xaf}{z}=\pm 1\\ \stackrel{\xaf}{T}\left(\stackrel{\xaf}{z},\tau \right)=0\text{when}\tau =0\end{array}$ (20)

$\stackrel{\xaf}{w}=0,\text{}\frac{\partial \stackrel{\xaf}{w}}{\partial \stackrel{\xaf}{r}}=0,\text{}\stackrel{\xaf}{u}=0\text{at}\text{\hspace{0.17em}}\stackrel{\xaf}{r}=1$ (21)

$\frac{\partial \stackrel{\xaf}{w}}{\partial \stackrel{\xaf}{r}}=0$ , $\stackrel{\xaf}{u}$ is finite at $\stackrel{\xaf}{r}=0$ (22)

$\stackrel{\xaf}{w}=\frac{\partial \stackrel{\xaf}{w}}{\partial \tau}=0$ and $\stackrel{\xaf}{u}=\frac{\partial \stackrel{\xaf}{u}}{\partial \tau}=0$ when $\tau =0$ (23)

3. Solution

To solve the governing Equation (17), the first step is to solve the electromagnetic field Equation (14) to obtain the magnetic field strength ${\stackrel{\xaf}{H}}_{r}$ and Lorentz force ${\stackrel{\xaf}{f}}_{z}$ . Solving Equation (16) gets the temperature field distribution $\stackrel{\xaf}{T}$ , and obtains ${\stackrel{\xaf}{F}}_{z}$ and ${N}^{T}$ according to formula (18).

3.1. Electromagnetic Field

Equation (14) is solved by separation of variables. In order to make non-homogeneous boundary Condition (19) homogeneous, new unknown function ${h}_{r}\left(\stackrel{\xaf}{z},\tau \right)$ is introduced

${h}_{r}\left(\stackrel{\xaf}{z},\tau \right)=\stackrel{\xaf}{H}\left(\stackrel{\xaf}{z},\tau \right)-\frac{\stackrel{\xaf}{z}+1}{2}\varphi \left(\tau \right)$ (24)

Substituting Formula (24) in to Equation (14), boundary conditions and initial Condition (19) can be written as follows

$\frac{\partial {h}_{r}}{\partial {\stackrel{\xaf}{z}}^{2}}=\frac{\partial {h}_{r}}{\partial \tau}+\frac{1+\stackrel{\xaf}{z}}{2}\frac{\text{d}\varphi \left(\tau \right)}{\text{d}\tau}$ (25)

${h}_{r}=0$ at $\stackrel{\xaf}{z}=\pm 1$ (26)

${h}_{r}\left(\stackrel{\xaf}{z},0\right)=-\varphi \left(0\right)\frac{1+\stackrel{\xaf}{z}}{2}$ when $\tau =0$ (27)

Assume that the solution of Equation (25) satisfy the Conditions (26) and (27) is following series form

${h}_{r}={\displaystyle \underset{n=1}{\overset{\infty}{\sum}}{\psi}_{n}\left(\tau \right)\mathrm{cos}\left({k}_{n}\stackrel{\xaf}{z}\right)}$ (28)

where $\psi \left(\tau \right)$ is the function of $\tau $ , ${k}_{n}$ is the position root of function $\mathrm{cos}\left({k}_{n}\right)=0$ ,that is

${k}_{n}=\frac{2n-1}{2}\text{\pi}$ ( $n=1,2,3,\cdots $ ) (29)

Substituting Formula (28) in the Equation (25), both sides of Equation (25) multiply $\mathrm{cos}\left({k}_{m}\stackrel{\xaf}{z}\right)$ , using the orthogonality of trigonometric function, it obtains

$\frac{\text{d}\psi \left(\tau \right)}{\text{d}\tau}+{k}_{n}^{2}\psi \left(\tau \right)+{\displaystyle {\int}_{-1}^{1}\frac{\text{d}\varphi \left(\tau \right)}{\text{d}\tau}}\frac{1+\stackrel{\xaf}{z}}{2}\mathrm{cos}\left({k}_{n}\stackrel{\xaf}{z}\right)\text{d}\stackrel{\xaf}{z}=0$ (30)

With initial Condition (22), the solution of Equation (30) is

${\psi}_{n}=\frac{{\left(-1\right)}^{n}}{{k}_{n}}{\stackrel{\xaf}{\psi}}_{n}\left(\tau \right)$ (31)

where $\stackrel{\xaf}{\psi}\left(\tau \right)={\displaystyle {\int}_{0}^{\tau}{\text{e}}^{-{k}_{n}^{2}\left(\tau -x\right)}\frac{\text{d}\phi \left(x\right)}{\text{d}x}\text{d}x}$

Thus, the expressions of non-dimensional magnetic field intensity ${\stackrel{\xaf}{H}}_{r}$ , induced current and Lorentz force respectively are

$\begin{array}{l}{\stackrel{\xaf}{H}}_{r}={\displaystyle \underset{n=1}{\overset{\infty}{\sum}}\frac{{\left(-1\right)}^{n}}{{k}_{n}}}\stackrel{\xaf}{\psi}\left(\tau \right)\mathrm{cos}\left({k}_{n}\stackrel{\xaf}{z}\right)+\frac{1+\stackrel{\xaf}{z}}{2}\phi \left(\tau \right)\\ {\stackrel{\xaf}{J}}_{\theta}={\displaystyle \underset{n=1}{\overset{\infty}{\sum}}{\left(-1\right)}^{n+1}\mathrm{sin}\left({k}_{n}\stackrel{\xaf}{z}\right)\stackrel{\xaf}{\psi}\left(\tau \right)}+\frac{1}{2}\phi \left(\tau \right)\\ {\stackrel{\xaf}{f}}_{z}=\left[{\displaystyle \underset{m=1}{\overset{\infty}{\sum}}{\left(-1\right)}^{m}{\stackrel{\xaf}{\psi}}_{m}\left(\tau \right)\mathrm{sin}\left({k}_{m}\stackrel{\xaf}{z}\right)}-\frac{1}{2}\phi \left(\tau \right)\right]\left[{\displaystyle \underset{n=1}{\overset{\infty}{\sum}}\frac{{\left(-1\right)}^{n}}{{k}_{n}}\stackrel{\xaf}{\psi}\left(\tau \right)\mathrm{cos}\left({k}_{n}\stackrel{\xaf}{z}\right)}+\frac{1+\stackrel{\xaf}{z}}{2}\phi \left(\tau \right)\right]\end{array}$ (32)

3.2. Temperature Field

Based on the separation of variables, assume that the solution of Equation (16) which satisfies the boundary Condition (20) shows the following form

$\stackrel{\xaf}{T}={\displaystyle \underset{i=1}{\overset{\infty}{\sum}}{b}_{i}}\left(\tau \right)\mathrm{cos}\left({a}_{i}\stackrel{\xaf}{z}\right)$ (33)

where ${a}_{i}$ is confirmed by $\mathrm{tan}{a}_{i}=\frac{\stackrel{\xaf}{h}}{{a}_{i}}$ . Both sides of Equation (12) multiply

$\mathrm{cos}\left({a}_{j}\stackrel{\xaf}{z}\right)$ , using the orthogonality of trigonometric function, it obtains

$\frac{\partial {b}_{i}\left(\tau \right)}{\partial \tau}+{\varsigma}_{1}{a}_{i}^{2}{b}_{i}\left(\tau \right)-\frac{{\stackrel{\xaf}{h}}^{2}+{a}_{i}^{2}}{\stackrel{\xaf}{h}+{\stackrel{\xaf}{h}}^{2}+{a}_{i}^{2}}{\displaystyle {\int}_{-1}^{1}\stackrel{\xaf}{Q}\left(\stackrel{\xaf}{z},\tau \right)\mathrm{cos}\left({a}_{i}\stackrel{\xaf}{z}\right)\text{d}\stackrel{\xaf}{z}}=0$ (34)

The solution of Equation (34) is

${b}_{i}\left(\tau \right)=\frac{{\stackrel{\xaf}{h}}^{2}+{a}_{i}^{2}}{\stackrel{\xaf}{h}+{\stackrel{\xaf}{h}}^{2}+{a}_{i}^{2}}\left[{\displaystyle \underset{m=1}{\overset{\infty}{\sum}}{\displaystyle \underset{n=1}{\overset{\infty}{\sum}}{L}_{imn}{\stackrel{^}{b}}_{imn}\left(\tau \right)+c\left(\tau \right)}}\right]$ (35)

where

$\begin{array}{l}{L}_{imn}=\stackrel{\xaf}{h}\mathrm{cos}\left({a}_{i}\right)\frac{2\left({k}_{m}^{2}+{k}_{n}^{2}-{a}_{i}^{2}\right)}{{\left(2{k}_{m}{k}_{n}\right)}^{2}-{\left({k}_{m}^{2}+{k}_{n}^{2}-{a}_{i}^{2}\right)}^{2}}\\ {\stackrel{^}{b}}_{imn}={\displaystyle {\int}_{0}^{\tau}{\text{e}}^{-{\varsigma}_{1}{a}_{i}^{2}\left(\tau -\xi \right)}{\stackrel{\xaf}{\psi}}_{n}\left(\xi \right){\stackrel{\xaf}{\psi}}_{m}\left(\xi \right)\text{d}\xi}\\ c\left(\tau \right)=\frac{\mathrm{sin}\left({a}_{i}\right)}{4{a}_{i}}\frac{4\omega \left(1-{\text{e}}^{-{\varsigma}_{1}{a}_{i}^{2}\tau}\right)+{a}_{i}^{4}{\varsigma}_{1}^{2}\left(1-\mathrm{cos}\left(2\omega \tau \right)\right)-2\omega {a}_{i}^{2}{\varsigma}_{1}\mathrm{sin}\left(2\omega \tau \right)}{4{\omega}^{2}{\varsigma}_{1}{a}_{i}^{2}+{a}_{i}^{6}{\varsigma}_{1}^{3}}\end{array}$

Substituting expression (35) into expression (33), it obtains

$\stackrel{\xaf}{T}\left(\stackrel{\xaf}{z},\tau \right)={\displaystyle \underset{i=1}{\overset{\infty}{\sum}}\frac{{\stackrel{\xaf}{h}}^{2}+{a}_{i}^{2}}{\stackrel{\xaf}{h}+\stackrel{\xaf}{h}+{a}_{i}^{2}}\left[{\displaystyle \underset{m=1}{\overset{\infty}{\sum}}{\displaystyle \underset{n=1}{\overset{\infty}{\sum}}{L}_{imn}{\stackrel{^}{b}}_{imn}+c\left(\tau \right)}}\right]}\mathrm{cos}\left({a}_{i}\stackrel{\xaf}{z}\right)$ (36)

3.3. Elastic Field

After obtaining Lorentz force ${\stackrel{\xaf}{f}}_{z}$ and temperature field distribution $\stackrel{\xaf}{T}$ , ${\stackrel{\xaf}{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 $\stackrel{\xaf}{w}\left(\stackrel{\xaf}{r},\tau \right)$ satisfies with the boundary condition (21) and (22), and has a separable form of time and space as follows

$\stackrel{\xaf}{w}=f\left(\tau \right)\left(1-2{\stackrel{\xaf}{r}}^{2}+{\stackrel{\xaf}{r}}^{4}\right)$ (37)

Substituting expression (37) into Equation (17a), the solution of equation (17a) satisfying with boundary Condition (21) and (22) is

$\begin{array}{c}\stackrel{\xaf}{u}\left(\stackrel{\xaf}{r},\tau \right)=kf\left(\tau \right)\left(\frac{24-16v}{15}\stackrel{\xaf}{r}+\frac{4v-2}{3}{\stackrel{\xaf}{r}}^{2}+\frac{1-4v}{15}{\stackrel{\xaf}{r}}^{4}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{f}^{2}\left(\tau \right)\left[\frac{6-3v}{6}\stackrel{\xaf}{r}+\left(1-3v\right){\stackrel{\xaf}{r}}^{3}-\frac{8-2v}{3}{\stackrel{\xaf}{r}}^{5}+\frac{6-v}{6}{\stackrel{\xaf}{r}}^{7}\right]\end{array}$ (38)

Applying Galerkin integral to Equation (17b), it obtains

$\begin{array}{l}{\displaystyle {\int}_{0}^{1}{\displaystyle {\int}_{0}^{\text{2\pi}}[\left(\stackrel{\xaf}{r}\frac{{\partial}^{3}\stackrel{\xaf}{w}}{\partial {\stackrel{\xaf}{r}}^{3}}+\frac{{\partial}^{2}\stackrel{\xaf}{w}}{\partial {\stackrel{\xaf}{r}}^{2}}-\frac{1}{\stackrel{\xaf}{r}}\frac{\partial \stackrel{\xaf}{w}}{\partial \stackrel{\xaf}{r}}\right)}}\\ \text{\hspace{0.05em}}-{\varsigma}_{2}\stackrel{\xaf}{r}\left\{\left[\frac{\partial \stackrel{\xaf}{u}}{\partial \stackrel{\xaf}{r}}+k\frac{\partial \stackrel{\xaf}{w}}{\partial \stackrel{\xaf}{r}}+\frac{1}{2}{\left(\frac{\partial \stackrel{\xaf}{w}}{\partial \stackrel{\xaf}{r}}\right)}^{2}+\nu \frac{\stackrel{\xaf}{u}}{\stackrel{\xaf}{r}}\right]-{N}^{T}\right\}\left(k+\frac{\partial \stackrel{\xaf}{w}}{\partial \stackrel{\xaf}{r}}\right)\\ \text{\hspace{0.05em}}-\left({\stackrel{\xaf}{P}}_{z}+{\stackrel{\xaf}{F}}_{z}\right){\stackrel{\xaf}{r}}^{2}-{\varsigma}_{3}{\displaystyle {\int}_{0}^{1}\stackrel{\xaf}{r}}\frac{{\partial}^{2}\stackrel{\xaf}{w}}{\partial {\tau}^{2}}\text{d}\stackrel{\xaf}{r}]\frac{\partial \stackrel{\xaf}{w}}{\partial f}\stackrel{\xaf}{r}\text{d}\theta \text{d}\stackrel{\xaf}{r}=0\end{array}$ (39)

The dynamic response formula expressed by the deflection of the shell can be obtained

${a}_{1}\frac{{\text{d}}^{2}f}{\text{d}{\tau}^{2}}+{a}_{2}{f}^{3}+{a}_{3}{f}^{2}+{a}_{4}f-{p}_{0}-{p}_{1}-{p}_{2}=0$ (40)

In which

$\begin{array}{l}{a}_{1}={\varsigma}_{3}{\displaystyle {\int}_{0}^{1}\stackrel{\xaf}{r}\left({\displaystyle {\int}_{0}^{1}\stackrel{\xaf}{r}}{G}_{1}\text{d}\stackrel{\xaf}{r}\right){G}_{1}\text{d}}\stackrel{\xaf}{r},\text{}{a}_{2}={\varsigma}_{2}{\displaystyle {\int}_{0}^{1}\left[{{G}^{\prime}}_{3}+\frac{1}{2}{{G}^{\prime}}_{1}^{2}+\frac{\nu}{\stackrel{\xaf}{r}}\right]{{G}^{\prime}}_{1}{G}_{1}{\stackrel{\xaf}{r}}^{2}\text{d}}\stackrel{\xaf}{r},\\ {a}_{3}={\varsigma}_{2}{\displaystyle {\int}_{0}^{1}{\stackrel{\xaf}{r}}^{2}\left\{k\left[\frac{v}{r}{G}_{3}+{{G}^{\prime}}_{3}+\frac{1}{2}v{{G}^{\prime}}_{1}^{2}\right]+\left[{{G}^{\prime}}_{2}{{G}^{\prime}}_{1}+k{{G}^{\prime}}_{1}^{2}+\frac{v}{r}{G}_{2}{{G}^{\prime}}_{1}\right]\right\}{G}_{1}\text{d}}\stackrel{\xaf}{r},\\ {a}_{4}={\displaystyle {\int}_{0}^{1}\left\{-\left[\stackrel{\xaf}{r}{{G}^{\u2034}}_{1}+{{G}^{\u2033}}_{1}-\frac{1}{\stackrel{\xaf}{r}}{{G}^{\prime}}_{1}\right]+{\varsigma}_{2}\stackrel{\xaf}{r}k\left[{{G}^{\prime}}_{4}+k{{G}^{\prime}}_{1}+\frac{v}{\stackrel{\xaf}{r}}{G}_{2}\right]\right\}{G}_{1}\stackrel{\xaf}{r}\text{d}}\stackrel{\xaf}{r},\\ {p}_{0}=\stackrel{\xaf}{P}{\displaystyle {\int}_{0}^{1}{\stackrel{\xaf}{r}}^{3}{G}_{1}\text{d}}{\stackrel{\xaf}{r}}_{z},\text{}{p}_{1}=-{\varsigma}_{2}{\stackrel{\xaf}{N}}^{T}{\displaystyle {\int}_{0}^{1}\left(k+{G}^{\prime}\right){G}_{1}{\stackrel{\xaf}{r}}^{2}\text{d}}\stackrel{\xaf}{r},\text{}{p}_{2}={\stackrel{\xaf}{F}}_{z}{\displaystyle {\int}_{0}^{1}{\stackrel{\xaf}{r}}^{3}{G}_{1}\text{d}}\stackrel{\xaf}{r}\end{array}$

where

$\begin{array}{l}{G}_{1}\left(\stackrel{\xaf}{r}\right)=1-2{\stackrel{\xaf}{r}}^{2}+{\stackrel{\xaf}{r}}^{4},\text{}{G}_{2}\left(\stackrel{\xaf}{r}\right)=\frac{24-16v}{15}\stackrel{\xaf}{r}+\frac{4v-2}{3}{\stackrel{\xaf}{r}}^{2}+\frac{1-4v}{15}{\stackrel{\xaf}{r}}^{4}\\ {G}_{3}\left(\stackrel{\xaf}{r}\right)=\frac{6-3v}{6}\stackrel{\xaf}{r}+\left(1-3v\right){\stackrel{\xaf}{r}}^{3}-\frac{8-2v}{3}{\stackrel{\xaf}{r}}^{5}+\frac{6-v}{6}{\stackrel{\xaf}{r}}^{7}\end{array}$

Numerical solution of the dynamic response Equation (40) is obtained by Runge-Kutta method, the response relation of deflection and time is obtained.

4. Numerical Example

The shallow conical thin shell is shown in Figure 1, which is made of aluminum and subjected to mechanical load and time-varying magnetic field. Assume that the expression of the function $\varphi \left(\tau \right)$ of time-varying applied magnetic field is shown as follows.

$\varphi \left(\tau \right)=\mathrm{sin}\left(\omega \tau \right)$ (41)

where $\omega $ is the non-dimensional angular frequency of magnetic field. The physical parameters of the shallow conical thin shell are as follows

$\mu =4\text{\pi}\times {10}^{-7}\text{\hspace{0.17em}}\text{H}/\text{m}$ , $\sigma =3.42\times {10}^{7}\text{\hspace{0.17em}}\text{s}/\text{m}$ , $c=2.7\times {10}^{3}\text{\hspace{0.17em}}\text{J}/\text{kg}$ ,

$\rho =0.9\times {10}^{3}\text{\hspace{0.17em}}\text{kg}/{\text{m}}^{\text{3}}$ , $k=92.6\times {10}^{-6}\text{\hspace{0.17em}}{\text{m}}^{\text{2}}/\text{s}$ , $\nu =0.3$ , $E=70\text{\hspace{0.17em}}\text{GPa}$ ,

$\alpha =24\times {10}^{-6}{k}^{-1}$ , $h=2\text{\hspace{0.17em}}\text{mm}$ , $a=0.2\text{\hspace{0.17em}}\text{m}$ , $\Lambda =204\text{\hspace{0.17em}}\text{1}/\text{m}$

According to the above analysis, numerical calculation is obtained by Matlab, the results is shown in Figures 2-7.

Figure 2 and Figure 3 demonstrate the curve of magnetic field ${\stackrel{\xaf}{H}}_{r}$ and

Figure 2. Curve of Magnetic field intensity ${\stackrel{\xaf}{H}}_{r}$ vary with time.

Figure 3. Curve of induced current ${\stackrel{\xaf}{J}}_{\theta}$ vary with time.

Figure 4. Curve of Lorentz force ${\stackrel{\xaf}{f}}_{z}$ vary with time.

Figure 5. Curve of temperature $\stackrel{\xaf}{T}$ vary with time.

induced current ${\stackrel{\xaf}{J}}_{\theta}$ in different thickness varying with time in the different thickness. It can be seen from Figure 2 that magnetic field ${\stackrel{\xaf}{H}}_{r}$ and induced current ${\stackrel{\xaf}{J}}_{\theta}$ show sinusoidal variation. The ${\stackrel{\xaf}{H}}_{r}$ amplitude gradually increases to 1 when z is in the interval [0, 1] in Figure 2. However, the induced current ${\stackrel{\xaf}{J}}_{\theta}$ amplitude deceases first then increases when z is in the interval [0, 1] in Figure 3. the outermost surface of .shallow conical shell has the maximum induced current amplitude. Figure 4 is the curve of Lorentz force ${\stackrel{\xaf}{f}}_{z}$ varying with time in different thickness when $\stackrel{\xaf}{r}=0.5$ , it can be seen that Lorentz force shows sinusoidal variation, and vibration amplitude of shallow conical shell

Figure 6. Curve of deflection $\omega $ vary with time in different mechanical load ( $\stackrel{\xaf}{r}=0$ , ${H}_{0}=0.01/\mu $ ).

Figure 7. Curve of deflection $\omega $ vary with time in different magnetic field intensity ( $\stackrel{\xaf}{r}=0$ , ${P}_{z}={10}^{4}\text{\hspace{0.17em}}\text{N}$ ).

decreases gradually from surface to central face.

Figure 5 is the curve of temperature $\stackrel{\xaf}{T}$ vary with time in different thickness when $\stackrel{\xaf}{s}=0.5$ , it can be seen that temperature reaches a steady state when the time is long enough,

Figure 6 shows the curve of deflection $\stackrel{\xaf}{w}$ vary with time in different mechanical load which is expressed by non-dimension mechanical load ${\stackrel{\xaf}{P}}_{z}={10}^{-3}{P}_{z}$ . It can be seen that mechanical load has an effect on the vibration amplitude of deflection $\stackrel{\xaf}{w}$ , but has no effect on vibration frequency. The vibration amplitude of deflection $\omega $ increases as mechanical load increases.

Figure 7 shows the curve of deflection $\stackrel{\xaf}{w}$ vary with time in different applied magnetic field intensity which is expressed by non-dimension magnetic field intensity ${\stackrel{\xaf}{H}}_{0}=100{H}_{0}$ . It can be seen that vibration amplitude increases with the increase of the magnetic field strength.

5. Conclusions

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.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite 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

References

- 1. Zhou, Y.H. and Zheng, X.J. (1999) Electromagnetic Solid Structural Mechanics. Science Press, Peking.
- 2. Bai, X.Z. (2006) Magnetic Base Plate and Shell Elasticity. Science Press, Peking.
- 3. Xing, Y.F. and Liu, B. (2010) A Differential Quadrature Analysis of Dynamic and Quasi-Static Magneto-Thermo-Elastic Stresses in a Conducting Rectangular Plate Subjected to an Arbitrary Variation of Magnetic Field. International Journal of Engineering Science, 48, 1944-1960. https://doi.org/10.1016/j.ijengsci.2010.06.010
- 4. Higuchi, M., Tanigawa, Y. and Kawamura, R. (2007) Magneto-Thermo-Elastic Stresses Induced by a Transient Magnetic Field in a Conducting Solid Circular Cylinder. Journal of Thermal Stresses, 44, 5316-5335.
- 5. Tanigawa, Y., Higuchi, M. and Kawamura, R. (2008) Dynamic and Quasi-Static Behaviors of Magneto-Thermo-Elastic Stresses in a Conducting Hollow Circular Cylinder Subjected to an Arbitrary Variation of Magnetic Field. International Journal of Mechanical Sciences, 50, 365-379. https://doi.org/10.1016/j.ijmecsci.2007.11.001
- 6. Tu, J.X., Wang, Z.R., Li, Y.Z. and Wang, P. (2014) Nonlinear Random Vibration of a Current Carrying Circular Plate in a Magnetic Field. Journal of Mechanical Strength, 36, 504-509.
- 7. Wang, P., Chen, S.M. and Wang, Z.R. (2013) Bifurcation and Chaos of a Thin Rectangular Plate Simply Supported with Large Deflection in a Coupled Environment of Heating, Force and Magnetic Field. Journal of Vibration and Shock, 32, 129-134.
- 8. Milazzo, A. (2013) A One-Dimensional Model for Dynamic Analysis of Generally Layered Magneto-Electro-Elastic Beams. Journal of Sound and Vibration, 332, 465-483. https://doi.org/10.1016/j.jsv.2012.09.004
- 9. Lang, Z. and Li, X. (2013) Buckling and Vibration Analysis of Functionally Graded Magneto-Electro-Thermo-Elastic Circular Cylindrical Shells. Applied Mathematical Modelling, 37, 2279-2292. https://doi.org/10.1016/j.apm.2012.05.023
- 10. Kattimani, S.C. and Ray, M.C. (2015) Control of Geometrically Nonlinear Vibrations of Functionally Graded Magneto-Electro-Elastic Plates. International Journal of Mechanical Sciences, 99, 154-167. https://doi.org/10.1016/j.ijmecsci.2015.05.012

Notation Index

1. a: radius of the shell

2. Thickness of the shell

3. $\phi $ : the cone Angle

4. ${E}_{\theta}$ : θ 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. $\alpha $ : Coefficient of linear thermal expansion

20. ${\sigma}_{r}$ : r direction of the stress

21. ${\sigma}_{\theta}$ : θ direction of the stress

22. A: Tensile rigidity

23. D: Flexural rigidity