Next Article in Journal
Exact Two-Dimensional Analytical Calculations for Magnetic Field, Electromagnetic Torque, UMF, Back-EMF, and Inductance of Outer Rotor Surface Inset Permanent Magnet Machines
Next Article in Special Issue
On the Mixed Dirichlet–Steklov-Type and Steklov-Type Biharmonic Problems in Weighted Spaces
Previous Article in Journal
Mathematical Analysis of a Prey–Predator System: An Adaptive Back-Stepping Control and Stochastic Approach
Previous Article in Special Issue
On a Problem Arising in Application of the Re-Quantization Method to Construct Asymptotics of Solutions to Linear Differential Equations with Holomorphic Coefficients at Infinity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Elastodiffusive Orthotropic Euler–Bernoulli Beam Considering Diffusion Flux Relaxation

by
Dmitry Tarlakovskii
1,2 and
Andrei Zemskov
1,2,*
1
Research Institute of Mechanics, Lomonosov Moscow State University, 119192 Moscow, Russia
2
Moscow Aviation Institute (National Research University), 125993 Moscow, Russia
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2019, 24(1), 23; https://doi.org/10.3390/mca24010023
Submission received: 29 December 2018 / Revised: 30 January 2019 / Accepted: 30 January 2019 / Published: 6 February 2019
(This article belongs to the Special Issue Mathematical Modeling in Physical Sciences)

Abstract

:
This article considers an unsteady elastic diffusion model of Euler–Bernoulli beam oscillations in the presence of diffusion flux relaxation. We used the model of coupled elastic diffusion for a homogeneous orthotropic multicomponent continuum to formulate the problem. A model of unsteady bending for the elastic diffusive Euler–Bernoulli beam was obtained using Hamilton’s variational principle. The Laplace transform on time and the Fourier series expansion by the spatial coordinate were used to solve the obtained problem.

1. Introduction

In order to solve various types of fundamental, applied, and technical problems, scientists and engineers often have to consider the unsteady interaction between various physical fields. An example of this is the interaction between mechanical and diffusion fields. Among the recent publications, it can be noted that [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19] were devoted to this problem. In particular, [1,2,3,4,5,6,7,8,9,10,11] considered the thermal effects, while the electromagnetic effects were studied in [12,13,14,15]. In addition, the analysis of fast unsteady processes in relatively short time periods requires the relaxation of thermal and diffusion perturbations to be taken into account [1,2,4,5,6,8,9,11,12,13].
Based on the reviewed publications, it can be concluded that numerical–analytical methods constructed on the Laplace and Fourier integral transformations are used to solve unsteady mechanodiffusion problems. In this case, the Durbin method is mainly used for the Laplace transform inversion. This method allows the Mellin integral to be expressed through the Fourier transform. Special numerical algorithms are used to do the inverse Fourier transform. Descriptions of these algorithms can be found in [20,21,22].
However, there is no universal algorithm for the inversion of Laplace transform, as noted in [23]. Each specific algorithm fits a certain class of functions. At the same time, the specificity of Laplace images influences the choice of suitable systems of functions with which one can approximate the unknown original values. In addition, Laplace images are very cumbersome to obtain by solving the coupled problems. It is not always possible in practice to verify the capability of using one or the other algorithm to find Laplace transform originals. Thus, it can be concluded that the Laplace transform inversion is the main mathematical complexity that arises when solving unsteady problems (mechanodiffusion problems in particular).
This article considers an unsteady elastic diffusion problem of Euler–Bernoulli beam oscillations in the presence of diffusion flux relaxation. It also proposes a method to construct a solution based on the use of the integral Laplace transform and expansion into series by eigenfunctions. The expansion coefficients are represented as rational functions, which simplifies the issue of inverse Laplace transform. Thus, the problem considered in the article is solved analytically.
It should be noted that analogues of the considered problem are absent in well-known scientific publications.

2. Problem Formulation

We consider the unsteady oscillation problem of elastic diffusion Euler–Bernoulli beam in the presence of diffusion flux relaxation. The beam in general formulation is under the action of longitudinal and transverse forces. Also, bending moments are at its ends. The load scheme is presented in Figure 1.
For the mathematical problem formulation, we use the coupled elastic diffusion continuum model in a rectangular Cartesian coordinate system, which has the following form [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18]:
u ¨ i = σ i j x j + F i , η ˙ ( q ) = J i ( q ) x i + Y ( q ) ( q = 1 , N ¯ ) .
where the point denotes a time derivative; σ i j and J i ( q ) are the stress tensor components and the diffusion flux vector, respectively. They are defined as follows (the beam material is perfect solid solution):
σ i j = C i j k l u k x l q = 1 N α i j ( q ) η ( q ) , J i ( q ) + τ ( q ) J ˙ i ( q ) = D i j ( q ) η ( q ) x j + Λ i j k l ( q ) 2 u k x j x l ( q = 1 , N ¯ ) .
Here, the dots denote the time derivative. All quantities in (1) and (2) are dimensionless. For them, the following notations are used:
x i = x i l , u i = u i l , τ = C t l , C i j k l = C i j k l C 1111 , C 2 = C 1111 ρ , α i j ( q ) = α i j ( q ) C 1111 , τ q = C τ ( q ) L , D i j ( q ) = D i j ( q ) C l , Λ i j k l ( q ) = m ( q ) D i j ( q ) α k l ( q ) n 0 ( q ) ρ R T 0 C l , F i = ρ l F i С 1111 , Y ( q ) = l Y ( q ) C ,
where t is time; x i * are rectangular Cartesian coordinates; ρ is the medium density; u i * are displacement vector components; C i j k l are elastic constant tensor components; T 0 is the initial temperature; D i j * ( q ) are the self-diffusion coefficients; R is the universal gas constant; m ( q ) is the molar mass; η ( q ) = n ( q ) n 0 ( q ) is the concentration increment of q -th component in the N -component medium; n 0 ( q ) and n ( q ) are the initial and actual concentrations (mass fractions); α i j * ( q ) are coefficients characterizing the medium volumetric changes due to diffusion; l is beam length; F i and Y ( q ) are mechanical and diffusive bulk perturbations; and τ ( q ) is the relaxation time of diffusion perturbations.
The formulation of the problem is completed by the initial and boundary conditions.
Initial conditions:
u i | τ = t 0 = u i 0 , u ˙ i | τ = t 0 = v i 0 , η ( q ) | τ = t 0 = η 0 ( q ) ( q = 1 , N ¯ ) .
Here, u i 0 , v i 0 , η 0 ( q ) are the given functions of spatial coordinates. Further in the paper, we assume that t 0 = 0 , u i 0 = 0 , v i 0 = 0 , η 0 ( q ) = 0 .
Boundary conditions (domain G is bounded; n i are components of the outer normal unit vector to G , G = Π u Π σ = Π η Π J ):
u i | Π u = U i , σ i j n j | Π σ = P i , η ( q ) | Π η = N ( q ) , ( J i ( q ) + τ q J ˙ i ( q ) ) | Π J = I i ( q ) ( τ > 0 , q = 1 , N ¯ ) .
The quantities on the boundary conditions on right sides are surface kinematic U i , N ( q ) and dynamic P i , I i ( q ) perturbations.
To construct the beam bending equations, a transition to the variational formulation of the problem (1)–(5) is used. According to Hamilton’s variational principle, the relations (1)–(5) can be regarded as a condition for the stationarity of a certain functional H ( u i , η ( q ) ) , whose variation is written thus:
δ H = t 1 t 2 d τ G ( u ¨ i σ i j x j F i ) δ u i d G + q = 1 N t 1 t 2 d τ G ( 1 + τ q τ ) ( η ˙ ( q ) + J i ( q ) x i Y ( q ) ) δ η ( q ) d G + + t 1 t 2 Π σ ( σ i j n j P i ) δ u i d S d τ + q = 1 N t 1 t 2 Π J ( J i ( q ) + τ q J ˙ i ( q ) I i ( q ) ) n i δ η ( q ) d S d τ .
Further, we will assume that
  • The axis O x 3 is the central axis of the cross section. In this case,
    D x 2 d x 2 d x 3 = 0 .
  • The side surface Π b is free from mechanical loads, i.e.,
    σ i j n j | Π b = 0 .
  • We also assume that there is no mass transfer through the side surface,
    ( J i ( q ) + τ q J ˙ i ( q ) ) | Π b = 0 .
  • The beam material is a homogeneous ortotropic continuum.
    The bending of beam is considered in plane x 1 O x 2 . Then u k = u k ( x 1 , x 2 , τ ) , k = 1 , 2 ; u 3 = 0 ; ε i 3 = 0 . Mass transfer occurs also in the plane x 1 O x 2 , i.e., η(q) = η(q)(x1,x2,τ).
  • Transverse deflections are considered small. Then, the linearization of the unknown quantities with respect to the variable x 2 will have the following form:
    u 1 ( x 1 , x 2 , τ ) = u ( x 1 , τ ) x 2 χ ( x 1 , τ ) ,   u 2 ( x 1 , x 2 , τ ) = v ( x 1 , τ ) + x 2 ψ ( x 1 , τ ) , η ( q ) ( x 1 , x 2 , τ ) = N q ( x 1 , τ ) + x 2 H q ( x 1 , τ ) .
  • The cross sections after deformation remain normal to the neutral line of the beam (Euler–Bernoulli ’s beam theory). Also, we assume that there are no deformations along the O x 2 axis [24,25].
    ε 22 = u 2 x 2 = ψ = 0 ψ = 0 , χ ( x 1 , τ ) = v ( x 1 , τ ) .
Then, Equation (10) will take the following form:
u 1 = u x 2 v , u 2 = v , u 3 = 0 , u = u ( x 1 , τ ) , v = v ( x 1 , τ ) .
Here, the prime denotes the derivative with respect to the variable x 1 .
From (8) and (9), it follows that,
σ 22 n 2 | Π b = σ 22 n 2 | Γ = σ 22 | γ 2 ( x 3 ) γ 1 ( x 3 ) = 0 ;
G σ 22 x 2 d G = 0 1 d x 1 D σ 22 x 2 d x 2 d x 3 = 0 1 d x 1 a b γ 2 ( x 3 ) γ 1 ( x 3 ) σ 22 x 2 d x 2 d x 3 = 0 1 d x 1 a b σ 22 | γ 2 ( x 3 ) γ 1 ( x 3 ) d x 3 = 0 .
The components of the stress tensor and the diffusion flow vector will have the following form:
σ 11 = C 11 u 1 x 1 + C 12 u 2 x 2 q = 1 N α 1 ( q ) η ( q ) = ( u x 2 v ) q = 1 N α 1 ( q ) ( N q + x 2 H q ) , σ 12 = C 66 ( u 1 x 2 + u 2 x 1 ) = 0 , σ 22 = C 12 u 1 x 1 + C 22 u 2 x 2 q = 1 N α 2 ( q ) η ( q ) = C 12 ( u x 2 v ) q = 1 N α 2 ( q ) ( N q + x 2 H q ) , J 1 ( q ) + τ q J ˙ 1 ( q ) = D 1 ( q ) η ( q ) x 1 + Λ 11 ( q ) 2 u 1 x 1 2 + Λ 12 ( q ) 2 u 2 x 1 x 2 = D 1 ( q ) ( N q + x 2 H q ) + Λ 11 ( q ) ( u x 2 v ) , J 2 ( q ) + τ q J ˙ 2 ( q ) = D 2 ( q ) η ( q ) x 2 + Λ 22 ( q ) 2 u 2 x 2 2 + Λ 21 ( q ) 2 u 1 x 1 x 2 = D 2 ( q ) H q Λ 21 ( q ) v ( q = 1 , N ¯ ) ,
where
C 11 = C 1111 , C 12 = C 1122 , C 22 = C 2222 , C 66 = C 1212 , α 1 ( q ) = α 11 ( q ) , α 2 ( q ) = α 22 ( q ) , D 1 ( q ) = D 11 ( q ) , D 2 ( q ) = D 22 ( q ) , Λ 11 ( q ) = Λ 1111 ( q ) , Λ 12 ( q ) = Λ 1122 ( q ) , Λ 21 ( q ) = Λ 2211 ( q ) , Λ 22 ( q ) = Λ 2222 ( q ) .
Substituting the relations (7), (10)–(14) into (6), we obtain the following:
δ H = t 1 t 2 d τ 0 1 [ F ( u ¨ u + q = 1 N α 1 ( q ) N ) n ] δ u d x 1 t 1 t 2 d τ 0 1 [ J 3 ( v ¨ v I V q = 1 N α 1 ( q ) H q ) + m ] δ v d x 1 + + t 1 t 2 d τ 0 1 ( F v ¨ q ) δ v d x 1 + q = 1 N t 1 t 2 d τ 0 1 [ F ( N ˙ q + τ ( q ) N ¨ q D 1 ( q ) N q + Λ 11 ( q ) u ) y ( q ) ] δ N q d x 1 + + q = 1 N t 1 t 2 d τ 0 1 [ J 3 ( H ˙ q + τ ( q ) H ¨ q D 1 ( q ) H q Λ 11 ( q ) v I V ) z ( q ) ] δ H q d x 1 + + t 1 t 2 [ J 3 ( v ¨ v q = 1 N α 1 ( q ) H q ) + m ] δ v | x 1 = 0 x 1 = 1 d τ + + t 1 t 2 { [ F ( u q = 1 N α 1 ( q ) N q ) N ] δ u + [ J 3 ( v + q = 1 N α 1 ( q ) H q ) + M ] δ v Q δ v } | x 1 = 0 x 1 = 1 d τ + + q = 1 N t 1 t 2 { [ F ( D 1 ( q ) N q + Λ 11 ( q ) u ) Γ ( q ) ] δ N q [ J 3 ( D 1 ( q ) H q + Λ 11 ( q ) v ) + Ω ( q ) ] δ H q } | x 1 = 0 x 1 = 1 d τ .
We introduce the following notation:
  • D d x 2 d x 3 = F is the cross-sectional area,
  • D x 2 2 d x 2 d x 3 = J 3 is moment of inertia of the beam section relative to the axis O x 3 ,
  • D F 1 d x 2 d x 3 = n is the linearly distributed axial load,
  • D F 1 x 2 d x 2 d x 3 = m is the linearly distributed moment,
  • D F 2 d x 2 d x 3 = q is the linearly distributed transverse load,
  • D ( Y ( q ) + τ ( q ) Y ˙ ( q ) ) d x 2 d x 3 = y ( q ) is the linear density of bulk mass transfer sources,
    D ( Y ( q ) + τ ( q ) Y ˙ ( q ) ) x 2 d x 2 d x 3 = z ( q ) , N ( x 1 ) = D P 1 ( x 1 , x 2 , x 3 , τ ) d x 2 d x 3 ,     M ( x 1 ) = D P 1 ( x 1 , x 2 , x 3 , τ ) x 2 d x 2 d x 3 ,     Q ( x 1 ) = D P 2 ( x 1 , x 2 , x 3 , τ ) d x 2 d x 3 , Γ ( q ) ( x 1 ) = D I 1 ( q ) ( x 1 , x 2 , x 3 , τ ) d x 2 d x 3 ,     Ω ( q ) ( x 1 ) = D x 2 I 1 ( q ) ( x 1 , x 2 , x 3 , τ ) d x 2 d x 3 .
Thus, the following boundary-value problems are the necessary conditions for the stationarity of the Hamilton functional considering (7)–(11):
u ¨ = u q = 1 N α 1 ( q ) N q + n F , N ˙ q + τ q N ¨ q = D 1 ( q ) N q Λ 11 ( q ) u + y ( q ) F ;
( u q = 1 N α 1 ( q ) N q ) | x 1 = 0 = N 0 F , ( u q = 1 N α 1 ( q ) N q ) | x 1 = 1 = N 1 F ;
( D 1 ( q ) N q + Λ 11 ( q ) u ) | x 1 = 0 = Γ 0 ( q ) F , ( D 1 ( q ) N q + Λ 11 ( q ) u ) | x 1 = 1 = Γ 1 ( q ) F ;
v ¨ F J 3 v ¨ = v I V + j = 1 N α 1 ( q ) H j q J 3 m J 3 , H ˙ q + τ q H ¨ q = D 1 ( q ) H q + Λ 11 ( q ) v I V + z ( q ) J 3 ;
( v + j = 1 N α 1 ( q ) H j ) | x 1 = 0 = M 0 J 3 , ( v + j = 1 N α 1 ( q ) H j ) | x 1 = 1 = M 1 J 3 ;
( v + j = 1 N α 1 ( q ) H j v ¨ ) | x 1 = 0 = Q 0 + m | x 1 = 0 J 3 , ( v + j = 1 N α 1 ( q ) H j v ¨ ) | x 1 = 1 = Q 1 + m | x 1 = 1 J 3 ;
( D 1 ( q ) H q + Λ 11 ( q ) v ) | x 1 = 0 = Ω 0 ( q ) J 3 , ( D 1 ( q ) H q + Λ 11 ( q ) v ) | x 1 = 1 = Ω 1 ( q ) J 3 ,
where
N 0 = N ( 0 ) , N 1 = N ( 1 ) , M 0 = M ( 0 ) , M 1 = M ( 1 ) , Q 0 = Q ( 0 ) , Q 1 = Q ( 1 ) , Γ 0 ( q ) = Γ ( q ) ( 0 ) , Γ 1 ( q ) = Γ ( q ) ( 1 ) , Ω 0 ( q ) = Ω ( q ) ( 0 ) , Ω 1 ( q ) = Ω ( q ) ( 1 ) .
In accordance with the variational principle of Lagrange, the boundary conditions (15)–(21) combined with the kinematic boundary conditions are as follows:
u | x 1 = 0 = U 0 , u | x 1 = 1 = U 1 ;
N q | x 1 = 0 = N q 0 , N q | x 1 = 1 = N q 1 ;
v | x 1 = 0 = V 0 , v | x 1 = 1 = V 1 ;
v | x 1 = 0 = V 0 , v | x 1 = 1 = V 1 ;
H q | x 1 = 0 = H q 0 , H q | x 1 = 1 = H q 1 ;
H q | x 1 = 0 = H q 0 , H q | x 1 = 1 = H q 1 .
Problems (15), (16), and (23) as well as (15), (17), and (22) are equivalent to one-dimensional problems of elastic diffusion for a layer. Their solutions are constructed using Laplace transform and decomposition into trigonometric Fourier series [26,27,28]. The papers [29,30] also considered problems with other combinations of boundary conditions that are solved by the method of equivalent boundary conditions.
Thus, the object of further consideration in this article will be initial boundary value problems for Equation (18) with a combination of boundary conditions (19), (20), and (24)–(27). Initial conditions will be assumed to be zero.

3. Method of Solution

We consider the problems (18), (19), (24), (26). The solutions of the problem are represented in the form ( q = 1 , N + 1 ¯ ):
v ( x , τ ) = k = 1 N + 2 0 τ [ G 1 k ( x , τ t ) f k 1 ( t ) + G 1 k ( 1 x , τ t ) f k 2 ( t ) ] d t + k = 1 N + 1 0 τ 0 1 G ˜ 1 k ( x , ξ , τ t ) F k ( ξ , t ) d ξ d t , H q ( x , τ ) = k = 1 N + 2 0 τ [ G q + 1 , k ( x , τ t ) f k 1 ( t ) + G q + 1 , k ( 1 x , τ t ) f k 2 ( t ) ] d t + + k = 1 N + 1 0 τ 0 1 G ˜ q + 1 , k ( x , ξ , τ t ) F k ( ξ , t ) d ξ d t .
Here, x = x 1 ; G m k and G ˜ m k are surface and bulk Green’s functions; F 1 ( x , τ ) = ( q + m ) / J 3 and F q + 1 ( x , τ ) = z ( q ) / J 3 are body forces entering into Equation (18); and f k l ( t ) are surface disturbances that have the following form:
f 11 ( τ ) = M 0 ( τ ) J 3 , f 12 ( τ ) = M 1 ( τ ) J 3 , f 21 ( τ ) = V 0 , f 22 ( τ ) = V 1 , f q + 2 , 1 ( τ ) = H q 0 ( τ ) , f q + 2 , 2 ( τ ) = H q 1 ( τ ) .
The surface Green’s functions G m k satisfy
G ¨ 1 k F J 3 G ¨ 1 k = G 1 k I V + j = 1 N α 1 ( j ) G j + 1 , k , G ˙ q + 1 , k + τ q G ¨ q + 1 , k = D 1 ( q ) G q + 1 , k + Λ 11 ( q ) G 1 k I V ,
( G 1 k + j = 1 N α 1 ( j ) G j + 1 , k ) | x 1 = 0 = δ 1 k δ ( τ ) , ( G 1 k + j = 1 N α 1 ( j ) G j + 1 , k ) | x 1 = 1 = 0 , G 1 k | x 1 = 0 = δ 2 k δ ( τ ) , G 1 k | x 1 = 1 = 0 , G q + 1 , k | x 1 = 0 = δ q + 2 , k δ ( τ ) , G q + 1 , k | x 1 = 1 = 0 .
The bulk Green’s functions G ˜ m k satisfy
G ˜ ¨ 1 k F J 3 G ˜ ¨ 1 k = G ˜ 1 k I V + j = 1 N α 1 ( j ) G ˜ j + 1 , k δ 1 k δ ( x ξ ) δ ( τ ) , G ˜ ˙ q + 1 , k + τ q G ˜ ¨ q + 1 , k = D 1 ( q ) G ˜ q + 1 , k + Λ 11 ( q ) G ˜ 1 k I V + δ q + 1 , k δ ( x ξ ) δ ( τ ) ,
and homogeneous boundary conditions corresponding to (30).
We consider the problem of finding the surface Green’s functions G m k . Applying the Laplace transformation with respect to time to (29) and (30), we get the following ( s is the Laplace transform parameter):
s 2 G 1 k L ( x , s ) F J 3 s 2 G 1 k L ( x , s ) = [ G 1 k I V ( x , s ) ] L + j = 1 N α 1 ( j ) G j + 1 , k L ( x , s ) , ( s + τ q s 2 ) G q + 1 , k L ( x , s ) = D 1 ( q ) G q + 1 , k L ( x , s ) + Λ 11 ( q ) [ G 1 k I V ( x , s ) ] L , ( G 1 k L + j = 1 N α j G j + 1 , k L ) | x 1 = 0 = δ 1 k , ( G 1 k L + j = 1 N α j G j + 1 , k L ) | x 1 = 1 = 0 , G 1 k L | x 1 = 0 = δ 2 k , G 1 k L | x 1 = 1 = 0 , G q + 1 , k L | x 1 = 0 = δ q + 1 , k , G q + 1 , k L | x 1 = 1 = 0 .
Next, we multiply each equation by sin λ n x and integrate in the interval [ 0 , 1 ] . The result is as follows:
k 1 ( λ n , s ) G 1 k L s ( λ n , s ) λ n 2 j = 1 N α 1 ( j ) G j + 1 , k L s ( λ n , s ) = F 1 k ( λ n ) , λ n 4 Λ 11 ( q ) G 1 k L s ( λ n , s ) + k q + 1 ( λ n , s ) G q + 1 , k L s ( λ n , s ) = F q + 1 , k ( λ n ) ,
where
k 1 ( λ n , s ) = ( λ n 2 + F J 3 ) s 2 + λ n 4 , F q + 1 , k ( λ n ) = 2 λ n ( Λ 11 ( q ) δ 1 k Λ 11 ( q ) λ n 2 δ 2 k + D 1 ( q ) δ q + 2 , k Λ 11 ( q ) j = 1 N α 1 ( j ) δ j + 2 , k ) , k q + 1 ( λ n , s ) = s + τ q s 2 + D 1 ( q ) λ n 2 , F 1 k ( λ n ) = 2 λ n δ 1 k + 2 λ n ( s 2 + λ n 2 ) δ 2 k , G m k L ( x , s ) = n = 1 G m k L s ( λ n , s ) sin λ n x , G m k L s ( λ n , s ) = 2 0 1 G m k L ( x , s ) sin λ n x d x , λ n = π n .
The solution of system (31) has the following form:
G 1 k L s ( λ n , s ) = P 1 k ( λ n , s ) P ( λ n , s ) , G q + 1 , 1 L s ( λ n , s ) = 2 λ n Λ 11 ( q ) k q + 1 ( λ n , s ) + P q + 1 , 1 ( λ n , s ) Q q ( λ n , s ) , G q + 1 , 2 L s ( λ n , s ) = 2 Λ 11 ( q ) λ n 3 k q + 1 ( λ n , s ) + P q + 1 , 2 ( λ n , s ) Q q ( λ n , s ) , G q + 1 , p + 2 L s ( λ n , s ) = 2 λ n ( D 1 ( q ) δ q p Λ 11 ( q ) α p ) k q + 1 ( λ n , s ) + P q + 1 , p + 2 ( λ n , s ) Q q ( λ n , s ) ( q , p = 1 , N ¯ , k = 1 , N + 1 ¯ ) .
Here,
P ( λ n , s ) = k 1 ( λ n , s ) Π ( λ n , s ) λ n 6 j = 1 N Λ 11 ( j ) α 1 ( j ) Π j ( λ n , s ) , Q q ( λ n , s ) = k q + 1 ( λ n , s ) P ( λ n , s ) , P 11 ( λ n , s ) = 2 λ n [ Π ( λ n , s ) λ n 2 j = 1 N α 1 ( j ) Λ 11 ( j ) Π j ( λ n , s ) ] , P 12 ( λ n , s ) = 2 λ n [ ( s 2 + λ n 2 ) Π ( λ n , s ) λ n 4 j = 1 N α 1 ( j ) Λ 11 ( j ) Π j ( λ n , s ) ] , P 1 , q + 2 ( λ n , s ) = 2 α 1 ( q ) λ n 3 [ D 1 ( q ) Π q ( λ n , s ) j = 1 N α 1 ( j ) Λ 11 ( j ) Π j ( λ n , s ) ] , P q + 1 , 1 ( λ n , s ) = 2 λ n 5 Λ 11 ( q ) [ Π ( λ n , s ) λ n 2 j = 1 N α 1 ( j ) Λ 11 ( j ) Π j ( λ n , s ) ] , P q + 1 , 2 ( λ n , s ) = 2 λ n 5 Λ 11 ( q ) [ ( s 2 + λ n 2 ) Π ( λ n , s ) λ n 4 j = 1 N α 1 ( j ) Λ 11 ( j ) Π j ( λ n , s ) ] , P q + 1 , p + 2 ( λ n , s ) = 2 α 1 ( p ) λ n 7 Λ 11 ( q ) [ D 1 ( p ) Π p ( λ n , s ) j = 1 N α 1 ( j ) Λ 11 ( j ) Π j ( λ n , s ) ] , Π ( λ n , s ) = j = 1 N k j + 1 ( λ n , s ) , Π j ( λ n , s ) = r = 1 , r j N k r + 1 ( λ n , s ) .
The Laplace originals of functions in (32) have the following form [28]:
G 1 k s ( λ n , τ ) = j = 1 2 N + 2 A 1 k ( j ) ( λ n ) e s j ( λ n ) τ , ξ 1 , 2 ( λ n ) = 1 ± 1 4 τ q D 1 ( q ) λ n 2 2 τ q , G q + 1 , 1 s ( λ n , τ ) = l = 1 2 [ 2 Λ 11 ( q ) λ n 1 + 2 τ q ξ l ( λ n ) + A q + 1 , 1 ( 2 N + 2 + l ) ( λ n ) ] e ξ l ( λ n ) τ + j = 1 2 N + 2 A q + 1 , 1 ( j ) ( λ n ) e s j ( λ n ) τ , G q + 1 , 2 s ( λ n , τ ) = l = 1 2 [ 2 Λ 11 ( q ) λ n 3 1 + 2 τ q ξ l ( λ n ) + A q + 1 , 2 ( 2 N + 2 + l ) ( λ n ) ] e ξ l ( λ n ) τ + j = 1 2 N + 2 A q + 1 , 2 ( j ) ( λ n ) e s j ( λ n ) τ , G q + 1 , p + 2 s ( λ n , τ ) = l = 1 2 [ 2 ( D 1 ( q ) δ q p Λ 11 ( q ) α 1 ( p ) ) λ n 1 + 2 τ q ξ l ( λ n ) + A q + 1 , p + 2 ( 2 N + 2 + l ) ( λ n ) ] e ξ l ( λ n ) τ + j = 1 2 N + 2 A q + 1 , p + 2 ( j ) ( λ n ) e s j ( λ n ) τ , A 1 k ( j ) ( λ n ) = P 1 k ( λ n , s j ( λ n ) ) P ( λ n , s j ( λ n ) ) ( j = 1 , N + 3 ¯ , k = 1 , N + 2 ¯ , r = 1 , N + 4 ¯ , q , p = 1 , N ¯ ) , A q + 1 , k ( r ) ( λ n ) = P q + 1 , k ( λ n , s r ( λ n ) ) Q q ( λ n , s r ( λ n ) ) , A q + 1 , k ( 2 N + 2 + l ) = P q + 1 , k ( λ n , ξ l ( λ n ) ) Q q ( λ n , ξ l ( λ n ) ) ( l = 1 , 2 ) ,
where s j ( λ n ) , j = 1 , 2 N + 2 ¯ are zeros of polynomial P ( λ n , s ) .

4. Examples

For the calculation example, we take the one-component (N = 1) aluminum medium with the following characteristics [31]:
C 1111 = 12.05 10 10 N m 2 , T 0 = 800 K , ρ = 2700 k g m 3 , α 11 ( 1 ) = 4.2 10 5 J m o l , D 11 ( 1 ) = 7.73 10 14 m 2 s , L = 0.1 m , τ ( q ) = 200 s .
We assume that the beam has a rectangular section with height h = 0.1 m and width b = 0.05 m . The geometrical characteristics of the section is as follows:
F = b h = 5.00 10 5 m 2 , J 3 = b h 3 12 = 4.17 10 10 m 4 .
Example 1. We assume in the boundary conditions (19), (24), and (26):
f 11 ( τ ) = M 0 ( τ ) J 3 = H ( τ ) , f 12 ( τ ) = M 1 ( τ ) J 3 = H ( τ ) , f 21 ( τ ) = V 0 ( τ ) = 0 , f 22 ( τ ) = V 1 ( τ ) = 0 , f 31 ( τ ) = H 20 ( τ ) = 0 , f 32 ( τ ) = H 21 ( τ ) = 0 .
Then, according to (28), in the absence case of volume perturbations, we have the following:
v ( x , τ ) = 0 τ [ G 11 ( x , τ t ) + G 11 ( 1 x , τ t ) ] H ( t ) d t = = 2 n = 1 sin λ n 2 cos λ n ( 1 2 x ) j = 1 4 A 11 ( j ) ( λ n ) e s j ( λ n ) τ 1 s j ( λ n ) , H 1 ( x , τ ) = 0 τ [ G 21 ( x , τ t ) + G 21 ( 1 x , τ t ) ] H ( t ) d t = = 2 n = 1 sin λ n 2 cos λ n ( 1 2 x ) l = 1 2 ( 2 Λ 11 ( 1 ) λ n 1 + 2 τ q ξ l ( λ n ) + A 21 ( 4 + l ) ( λ n ) ) e ξ l ( λ n ) τ 1 ξ l ( λ n ) + + 2 n = 1 sin λ n 2 cos λ n ( 1 2 x ) j = 1 4 A 21 ( j ) ( λ n ) e s j ( λ n ) τ 1 s j ( λ n ) .
The calculation results are shown in Figure 2. It should be noted that the effects associated with the diffusion flux relaxation for given perturbations (33) do not appear at all. The graphs for the beam deflection v ( x , τ ) and the concentration increment H 1 ( x , τ ) without considering the relaxation have the same view as in Figure 2.
Example 2. Now, we assume in the boundary conditions (19), (24), and (26) that
f 11 ( τ ) = M 0 ( τ ) J 3 = 0 , f 12 ( τ ) = M 1 ( τ ) J 3 = 0 , f 21 ( τ ) = V 0 ( τ ) = 0 , f 22 ( τ ) = V 1 ( τ ) = 0 , f 31 ( τ ) = H 20 ( τ ) = H ( τ ) , f 32 ( τ ) = H 21 ( τ ) = H ( τ ) .
That (28)
v ( x , τ ) = 0 τ [ G 13 ( x , τ t ) + G 13 ( 1 x , τ t ) ] H ( t ) d t = = 2 n = 1 sin λ n 2 cos λ n ( 1 2 x ) j = 1 4 A 13 ( j ) ( λ n ) e s j ( λ n ) τ 1 s j ( λ n ) , H 1 ( x , τ ) = 0 τ [ G 22 ( x , τ t ) + G 22 ( 1 x , τ t ) ] H ( t ) d t = = 2 n = 1 sin λ n 2 cos λ n ( 1 2 x ) l = 1 2 ( 2 ( D 1 ( 1 ) α 1 ( 1 ) Λ 11 ( 1 ) ) λ n 1 + 2 τ q ξ l ( λ n ) + A 23 ( 4 + l ) ( λ n ) ) e ξ l ( λ n ) τ 1 ξ l ( λ n ) + + 2 n = 1 sin λ n 2 cos λ n ( 1 2 x ) j = 1 4 A 23 ( j ) ( λ n ) e s j ( λ n ) τ 1 s j ( λ n ) .
H 1 ( x , τ ) = 0 τ [ G 22 ( x , τ t ) + G 22 ( 1 x , τ t ) ] H ( t ) d t = = 2 n = 1 sin λ n 2 cos λ n ( 1 2 x ) l = 1 2 ( 2 ( D 1 ( 1 ) α 1 ( 1 ) Λ 11 ( 1 ) ) λ n 1 + 2 τ q ξ l ( λ n ) + A 23 ( 4 + l ) ( λ n ) ) e ξ l ( λ n ) τ 1 ξ l ( λ n ) + + 2 n = 1 sin λ n 2 cos λ n ( 1 2 x ) j = 1 4 A 23 ( j ) ( λ n ) e s j ( λ n ) τ 1 s j ( λ n ) .  
The calculation results are shown in Figure 3.
Relaxation effects are manifested under the condition τ τ q (dimensionless time 10 7 is corresponds to 150 s), as can be seen in Figure 3. The relaxation effects fade out with increasing time, which is demonstrated in Figure 4. A similar situation is observed for the concentration increment H q ( x , τ ) . Based on the completed study, it can be concluded that the relaxation effects manifestation essentially depends on both the material properties and the type of specified perturbations. The results obtained in this article are consistent with the results obtained by other researchers, for example, [6,12,32,33].

5. Conclusions

The coupled unsteady oscillations model of an elastic diffusion Euler–Bernoulli beam was constructed using Hamilton’s variational principle. An algorithm for constructing surface Green’s functions was proposed. The use of the unknown function expansion into series by eigenfunctions allowed us to solve the problem associated with the Laplace transform inversion. Such an approach made it possible to find an analytical solution to the oscillation problem of an elastic diffusion Euler–Bernoulli beam.
On the basis of the developed model, the interaction of mechanical and diffusion fields was investigated. The influence of diffusion relaxation effects on the beam was analyzed. A number of test calculations showed that relaxation effects appeared only at the initial moments of time and under certain boundary conditions. The results of the calculations are presented in analytical and graphical forms.

Author Contributions

D.T. and A.Z. contributed equally to the writing of the main manuscript and preparation of the figures. Both authors reviewed the manuscript. Both authors contributed equally to the paper.

Funding

This work was funded by a subsidy from RFBR (Project 17-08-00663 А).

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Afram, A.Y.; Khader, S.E. 2D Problem for a Half-Space under the Theory of Fractional Thermoelastic Diffusion. Am. J. Sci. Ind. Res. 2014, 6, 47–57. [Google Scholar]
  2. Atwa, S.Y.; Egypt, Z. Generalized Thermoelastic Diffusion with Effect of Fractional Parameter on Plane Waves Temperature-Dependent Elastic Medium. J. Mater. Chem. Eng. 2013, 1, 55–74. [Google Scholar]
  3. Belova, I.V.; Murch, G.E. Thermal and diffusion-induced stresses in crystalline solids. J. Appl. Phys. 1995, 77, 127–134. [Google Scholar] [CrossRef]
  4. Choudhary, S.; Deswal, S. Mechanical Loads on a Generalized Thermoelastic Medium with Diffusion. Meccanica 2010, 45, 401–413. [Google Scholar] [CrossRef]
  5. Elhagary, M.A. Generalized thermoelastic diffusion problem for an infinitely long hollow cylinder for short times. Acta Mech. 2011, 218, 205–215. [Google Scholar] [CrossRef]
  6. El-Sayed, A.M. A two-dimensional generalized thermoelastic diffusion problem for a half-space. Math. Mech. Solids 2016, 21, 1045–1060. [Google Scholar] [CrossRef]
  7. Knyazeva, A.G. Model of medium with diffusion and internal surfaces and some applied problems. Mater. Phys. Mech. 2004, 7, 29–36. [Google Scholar]
  8. Kumar, R.; Chawla, V. Green’s Functions in Orthotropic Thermoelastic Diffusion Media. Eng. Anal. Bound. Elem. 2012, 36, 1272–1277. [Google Scholar] [CrossRef]
  9. Olesiak, Z.S.; Pyryev, Yu.A. A coupled quasi-stationary problem of thermodiffusion for an elastic cylinder. Int. J. Eng. Sci. 1995, 33, 773–780. [Google Scholar] [CrossRef]
  10. Pidstryhach, Y.S. Differential equations of the problem of thermodiffusion in a solid deformable isotropic body. Dop. Akad. Nauk USSR 1961, 2, 169–172. [Google Scholar]
  11. Sherief, H.H.; El-Maghraby, N.M. A Thick Plate Problem in the Theory of Generalized Thermoelastic Diffusion. Int. J. Thermophys. 2009, 30, 2044–2057. [Google Scholar] [CrossRef]
  12. Aouadi, M. Variable electrical and thermal conductivity in the theory of generalized thermoelastic diffusion. Z. Angew. Math. Phys. 2005, 57, 350–366. [Google Scholar] [CrossRef]
  13. Deswal, S.; Kalkal, K. A two-dimensional generalized electro-magneto-thermoviscoelastic problem for a half-space with diffusion. Int. J. Therm. Sci. 2011, 50, 749–759. [Google Scholar] [CrossRef]
  14. Tarlakovskii, D.V.; Vestyak, V.A.; Zemskov, A.V. Dynamic Processes in Thermoelectromagnetoelastic and Thermoelastodiffusive Media. In Encyclopedia of Thermal Stress; Springer: London, UK, 2014; Volume 2, pp. 1064–1071. [Google Scholar]
  15. Zhang, J.; Li, Y. A Two-Dimensional Generalized Electromagnetothermoelastic Diffusion Problem for a Rotating Half-Space. Hindawi Publ. Corp. Math. Probl. Eng. 2014, 2014, 1–12. [Google Scholar] [CrossRef]
  16. Chu, J.L.; Lee, S. Diffusion-induced stresses in a long bar of square cross section. J. Appl. Phys. 1993, 73, 3211–3219. [Google Scholar] [CrossRef]
  17. Freidin, A.B.; Korolev, I.K.; Aleshchenko, S.P.; Vilchevskaya, E.N. Chemical affinity tensor and chemical reaction front propagation: Theory and FE-simulations. J. Fract. 2016, 202, 245–259. [Google Scholar] [CrossRef]
  18. Hwang, C.C.; Chen, K.M.; Hsieh, J.Y. Diffusion-induced stresses in a long bar under an electric field. J. Phys. D Appl. Phys. 1994, 27, 2155–2162. [Google Scholar] [CrossRef]
  19. Indeitsev, D.A.; Semenov, B.N.; Sterlin, M.D. The Phenomenon of Localization of Diffusion Process in a Dynamically Deformed Solid. Dokl. Phys. 2012, 57, 171–173. [Google Scholar] [CrossRef]
  20. Crump, K. Numerical Inversion of Laplace Transforms Using a Fourier Series Approximation. J. Assoc. Comput. Mach. 1976, 23, 89–96. [Google Scholar] [CrossRef]
  21. Durbin, F. Numerical inversion of Laplace transforms: An efficient improvement to Dubner and Abate’s method. Comput. J. 1974, 17, 371–376. [Google Scholar] [CrossRef]
  22. Honig, G.; Hirdes, U. A method for the numerical inversion of Laplace transforms. J. Comput. Appl. Math. 1984, 10, 113–132. [Google Scholar] [CrossRef]
  23. Poroshina, N.I.; Ryabov, V.M. Methods for Laplace Transform Inversion; Vestnik of the St. Petersburg 239 University: Cambridge, UK, 2011; Volume 44, pp. 214–222. [Google Scholar]
  24. Mikhailova, E.Yu.; Tarlakovskii, D.V.; Fedotenkov, G.V. Obshchaya Teoriya Uprugikh Obolochek; MAI: Moscow, Russia, 2018. [Google Scholar]
  25. Gorshkov, A.G.; Medvedsky, A.L.; Rabinsky, L.N.; Tarlakovsky, D.V. Volny v Sploshnykh Sredakh; Fizmatlit: Moscow, Russia, 2004. [Google Scholar]
  26. Davydov, S.A.; Zemskov, A.V.; Tarlakovskii, D.V. An Elastic Half-Space under the Action of One-Dimensional Time-Dependent Diffusion Perturbations. Lobachevskii J. Math. 2015, 36, 503–509. [Google Scholar] [CrossRef]
  27. Davydov, S.A.; Vestyak, A.V.; Zemskov, A.V.; Tarlakovskii, D.V. Unsteady one-dimensional problem of thermoelastic diffusion for homogeneous multicomponent continuum with a plane boundaries. Uchenye Zapiski Kazanskogo Univ. Ser. Fiziko-Matematicheskie Nauki 2018, 160, 183–196. [Google Scholar]
  28. Davydov, S.A.; Zemskov, A.V. Unsteady One-dimensional Perturbations in Multicomponent Thermoelastic Layer with Cross-diffusion Effect. J. Phys. Conf. Ser. 2018, 1129, 012009. [Google Scholar] [CrossRef]
  29. Zemskov, A.V.; Tarlakovskiy, D.V. Method of the equivalent boundary conditions in the unsteady problem for elastic diffusion layer. Mater. Phys. Mech. 2015, 23, 36–41. [Google Scholar]
  30. Davydov, S.A.; Zemskov, A.V.; Igumnov, L.A.; Tarlakovskiy, D.V. Non-stationary model of mechanical diffusion for half-space with arbitrary boundary conditions. Mater. Phys. Mech. 2016, 28, 72–76. [Google Scholar]
  31. Grigoriev, I.S.; Meylikhov, I.Z. Fizicheskiye Velichiny: Sprovochnik; Energoatomizdat: Moscow, Russia, 1991. [Google Scholar]
  32. Polyanin, A.D.; Vyazmin, A.V. differential-difference heat-conduction and diffusion models and equations with a finite relaxation time. Theor. Found. Chem. Eng. 2013, 47, 217–224. [Google Scholar] [CrossRef]
  33. Szekeres, A.; Fekete, B. Continuummechanics—Heat Conduction—Cognition. Period. Polytech. Mech. Eng. 2015, 59, 8–15. [Google Scholar] [CrossRef]
Figure 1. The formulation of the problem.
Figure 1. The formulation of the problem.
Mca 24 00023 g001
Figure 2. The results of calculations using Equation (34): (a) The beam deflections v ( x , τ ) ; (b) The concentration increment H 1 ( x , τ ) .
Figure 2. The results of calculations using Equation (34): (a) The beam deflections v ( x , τ ) ; (b) The concentration increment H 1 ( x , τ ) .
Mca 24 00023 g002
Figure 3. The beam deflections v ( x , τ ) , x = 0.5 . The solid line is the solution with relaxation of diffusion fluxes, and the dashed line is the solution without relaxation of diffusion fluxes: (a) Effect of diffusion processes relaxation on the beam deflections at τ ~ 10 7 ; (b) Effect of diffusion processes relaxation on the beam deflections at τ ~ 10 8 .
Figure 3. The beam deflections v ( x , τ ) , x = 0.5 . The solid line is the solution with relaxation of diffusion fluxes, and the dashed line is the solution without relaxation of diffusion fluxes: (a) Effect of diffusion processes relaxation on the beam deflections at τ ~ 10 7 ; (b) Effect of diffusion processes relaxation on the beam deflections at τ ~ 10 8 .
Mca 24 00023 g003
Figure 4. The beam deflections v ( x , τ ) , x = 0.5 . The solid line is the solution with relaxation of diffusion fluxes, and the dashed line is the solution without relaxation of diffusion fluxes: (a) Effect of diffusion processes relaxation on the beam deflections at τ ~ 10 9 ; (b) Effect of diffusion processes relaxation on the beam deflections at τ ~ 10 12 .
Figure 4. The beam deflections v ( x , τ ) , x = 0.5 . The solid line is the solution with relaxation of diffusion fluxes, and the dashed line is the solution without relaxation of diffusion fluxes: (a) Effect of diffusion processes relaxation on the beam deflections at τ ~ 10 9 ; (b) Effect of diffusion processes relaxation on the beam deflections at τ ~ 10 12 .
Mca 24 00023 g004

Share and Cite

MDPI and ACS Style

Tarlakovskii, D.; Zemskov, A. An Elastodiffusive Orthotropic Euler–Bernoulli Beam Considering Diffusion Flux Relaxation. Math. Comput. Appl. 2019, 24, 23. https://doi.org/10.3390/mca24010023

AMA Style

Tarlakovskii D, Zemskov A. An Elastodiffusive Orthotropic Euler–Bernoulli Beam Considering Diffusion Flux Relaxation. Mathematical and Computational Applications. 2019; 24(1):23. https://doi.org/10.3390/mca24010023

Chicago/Turabian Style

Tarlakovskii, Dmitry, and Andrei Zemskov. 2019. "An Elastodiffusive Orthotropic Euler–Bernoulli Beam Considering Diffusion Flux Relaxation" Mathematical and Computational Applications 24, no. 1: 23. https://doi.org/10.3390/mca24010023

APA Style

Tarlakovskii, D., & Zemskov, A. (2019). An Elastodiffusive Orthotropic Euler–Bernoulli Beam Considering Diffusion Flux Relaxation. Mathematical and Computational Applications, 24(1), 23. https://doi.org/10.3390/mca24010023

Article Metrics

Back to TopTop