^{1}

^{1}

Novel exact solutions of one-dimensional transient dynamic piezoelectric problems for thickness polarized layers and disks, or length polarized rods, are obtained. The solutions are derived using a time-domain Green’s function method that leads to an exact analytical recursive procedure which is applicable for a wide variety of boundary conditions including nonlinear cases. A nonlinear damper boundary condition is considered in more detail. The corresponding nonlinear relationship between stresses and velocities at a current time moment is used in the recursive procedure. In addition to the exact recursive procedure that is effective for calculations, some new practically important explicit exact solutions are presented. Several examples of the time behavior of the output electric potential difference are given to illustrate the effectiveness of the proposed exact approach.

Piezoelectric materials and devices have been widely used in many technical applications. Nowadays, the coupling between electrical and mechanical beha- viors is used in different devices based both on the so-called “direct piezoelectric effect” or the “converse piezoelectric effect” [

Some newer relevant applications include (among others) the high voltage generation from transient dynamic impact processes in vehicles [

Analysis of operating electrical and mechanical parameters of such processes can be done by using various analytical and numerical methods. Although ana- lytical approaches are limited to rather simple geometries and other restrictions (homogeneous or piecewise-homogeneous bodies, linear governing equations, etc.), they often provide exact solutions.

Analytical methods have been successfully used for many transient dynamic one-dimensional piezoelectric problems [

In this paper, a time-domain Green’s function method is implemented for so- lution of one-dimensional transient dynamic piezoelectric problems for thick- ness polarized disks or length polarized rods. This method stems from a time- domain representation formulas approach for transient dynamic piezoe-lectric problems described in [

Consider a transversely isotropic homogeneous piezoelectric material (piezoe- lectric element) with the x 3 -axis as the poling direction and the x 1 − x 2 plane as the isotropic plane. Let this piezoelectric material occupy a disk (or a cylinder) Ω bounded in x 3 -direction by planes x 3 = 0 and x 3 = h where h > 0 is the thickness of the disk (or the length of the cylinder). Consider a uniaxial strain state or a stress stress state in x 3 direction when there is only one non- zero component of strain γ 33 or stress σ 33 the others being zero. We assume that the non-zero stress and strain components, and also the displacement u 3 and electric displacement D 3 in the x 3 -direction, and the electric potential ϕ , depend only on x 3 and t which is usually the case for a longitudinal mode piezoelectric element [

σ 33 = σ 33 ( x 3 , t ) , D 3 = D 3 ( x 3 , t ) , u 3 = u 3 ( x 3 , t ) , ϕ = ϕ ( x 3 , t ) in Ω (1)

Under conditions (1), we can use the following one-dimensional constitutive equations (both for the uniaxial strain state and for the uniaxial stress state) that relate the mechanical and electrical fields in (1):

σ 33 = C u 3 , 3 + e ϕ , 3 ; D 3 = e u 3 , 3 − ϵ ϕ , 3 (2)

where coefficients are different for the uniaxial strain and uniaxial stress cases.

Then the corresponding equations of motions can be written as

ρ u ¨ 3 − C u 3 , 33 − e ϕ , 33 = b 3 ; − D 3 , 3 = ϵ ϕ , 33 − e u 3 , 33 = − q in Ω (3)

where b 3 = b 3 ( x 3 , t ) and q = q ( x 3 , t ) denote the body force in x 3 -direction and electric charge.

To simplify further notations we will denote x 3 and derivatives with respect to x 3 by x and the prime, respectively, and will skip subindex 3 for the elastic displacement, electric displacement and body force components presented in (3). Then system (3) becomes

ρ u ¨ − C u ″ − e ϕ ″ = b ; − D ′ = ϵ ϕ ″ − e u ″ = − q in Ω . (4)

The Green’s functions for vector { u , ϕ } can be obtained using concentrated impulses instead of b or q in (4) when Ω is substituted by the infinite media.

Since ϕ ″ can be expressed through u ″ due to the second equation in (4), then the first equation in (4) can be presented as the one-dimensional wave equation for displacement u :

ρ u ¨ − C D u ″ = b − e ϵ q (5)

where

C D = C + e 2 ϵ

is the Young’s modulus measured at constant D .

The wave speed corresponding to Equation (5) is denoted below by

c = C D ρ

The Green’s function for u corresponding to load { b = δ ( x ) δ ( t ) , q = 0 } is the well-known Green’s function for the one-dimensional wave Equation (5):

U ( x , t ) = 1 2 ρ c H ( t − | x | / c ) (6)

where H ( t ) is the Heaviside step function (right-continuous), i.e. H ( t ) = 0 for t < 0 and H ( t ) = 1 for t ≥ 0 .

The corresponding Green’s function for ϕ is calculated using the second equation in (4):

U ϕ ( x , t ) = e ϵ U ( x , t ) (7)

Based on (6) and (7), the representation formula for the displacement vector in 3-D case described in [

u ( x , t ) = 1 2 [ u ( 0 , t − x / c ) + u ( h , t − ( h − x ) / c ) ] − 1 2 ρ c ∫ − ∞ t − x / c [ σ ( 0 , τ ) + e ϵ D ( 0 , τ ) ] d τ + 1 2 ρ c ∫ − ∞ t − ( h − x ) / c [ σ ( h , τ ) + e ϵ D ( h , τ ) ] d τ + 1 2 ρ c ∫ 0 h ∫ − ∞ t − | ξ − x | / c [ b ( ξ , τ ) − e ϵ q ( ξ , τ ) ] d τ d ξ in Ω , (8)

where σ 33 is denoted by σ and it is taken into account that the outward normals to the lower and upper boundaries of the layer 0 ≤ x ≤ h have opposite directions.

In many practical applications, the electric volume charges are absent. There- fore, we consider henceforth only the case when q = 0 . Then the terms related to D in the above expression can be simplified since, based on Equation (4) in this case, D ( x , t ) is spatially uniform:

D ( x , t ) = D ( t ) . (9)

Due to the property (9) the representation formula (8) can be rewritten as

u ( x , t ) = 1 2 [ u ( 0 , t − x / c ) + u ( h , t − ( h − x ) / c ) ] − 1 2 ρ c ∫ − ∞ t − x / c [ σ ( 0 , τ ) + e ϵ D ( τ ) ] d τ + 1 2 ρ c ∫ − ∞ t − ( h − x ) / c [ σ ( h , τ ) + e ϵ D ( τ ) ] d τ + 1 2 ρ c ∫ 0 h ∫ − ∞ t − | ξ − x | / c b ( ξ , τ ) d τ d ξ in Ω . (10)

To obtain a representation formula for ϕ ( x , t ) , let us consider an auxiliary function

ψ ( x , t ) = ϕ ( x , t ) − e ϵ u ( x , t ) (11)

that has the following connection to the electric displacement:

D = − ϵ ψ ′ .

According to (3), ψ ″ = − D ′ / ϵ = 0 . Then, using the corresponding Green’s function | x | / 2 and Equation (11), we get a representation formula for ψ ( x , t ) involving only boundary value of function ψ ( x , t ) and a spatially uniform electric displacement:

ψ ( x , t ) = 1 2 [ ψ ( 0 , t ) + ψ ( h , t ) ] + h − 2 x 2 ϵ D ( t ) in Ω . (12)

Formulas (12) and (11) lead to the following expression for ϕ ( x , t ) :

ϕ ( x , t ) = 1 2 [ ϕ ( 0 , t ) + ϕ ( h , t ) ] − e 2 ϵ [ u ( 0 , t ) + u ( h , t ) ] + e ϵ u ( x , t ) + h − 2 x 2 ϵ D ( t ) in Ω . (13)

After u ( x , t ) is calculated, ϕ ( x , t ) can be determined using this calculated value, D ( t ) and boundary values of ϕ ( x , t ) .

The representation formula (10) allows us to get representation formulas for the velocity v ( x , t ) = u ˙ ( x , t ) and stress σ ( x , t ) . Differentiating (10) with respect to time provides the following representation formula for the velocity:

v ( x , t ) = 1 2 [ v ( 0 , t − x / c ) + v ( h , t − ( h − x ) / c ) ] − 1 2 ρ c [ σ ( 0 , t − x / c ) + e ϵ D ( t − x / c ) ] + 1 2 ρ c [ σ ( h , t − ( h − x ) / c ) + e ϵ D ( t − ( h − x ) / c ) ] + 1 2 ρ c ∫ 0 h b ( ξ , t − | ξ − x | / c ) d ξ in Ω . (14)

To get a representation formula for the stress we need to use the first contitu- tive equation from (2) (in the new notations introduced after equations (3)) and expression (13) which gives the following expression for the stress:

σ ( x , t ) = C u ′ ( x , t ) + e [ e ϵ u ′ ( x , t ) − 1 ϵ D ( t ) ] = C D u ′ ( x , t ) − e ϵ D ( t ) . (15)

After differentiating (10) with respect to x and substituting the result into (15) we get

σ ( x , t ) = ρ c 2 [ v ( 0 , t − x / c ) + v ( h , t − ( h − x ) / c ) ] + 1 2 [ σ ( 0 , t − x / c ) + σ ( h , t − ( h − x ) / c ) ] + e 2 ϵ [ D ( t − x / c ) + D ( t − ( h − x ) / c ) − 2 D ( t ) ] + 1 2 ∫ 0 h b ( ξ , t − | ξ − x | / c ) sgn ( x − ξ ) d ξ in Ω . (16)

A representation formula for D ( x , t ) is not needed under assumption that q = 0 since the electric displacement is uniform in space in this case and deter- mined solely by the electric boundary conditions.

The velocity representation formula (14) generates two boundary equations when x tends to the upper and lower boundaries of the piezoelectric element, that is, when x tends to h or 0:

v ( h , t ) = v ( 0 , t − θ ) + 1 ρ c [ σ ( h , t ) − σ ( 0 , t − θ ) ] + e ρ c ϵ [ D ( t ) − D ( t − θ ) ] + 1 ρ c ∫ 0 h b ( ξ , t − θ + ξ / c ) d ξ , (17)

v ( 0, t ) = v ( h , t − θ ) + 1 ρ c [ σ ( h , t − θ ) − σ ( 0, t ) ] + e ρ c ϵ [ D ( t − θ ) − D ( t ) ] + 1 ρ c ∫ 0 h b ( ξ , t − ξ / c ) d ξ (18)

where θ denotes the time taken by the elastic wave to travel the thickness of the piezoelectric layer:

θ = h c .

Similarly, the stress representation formula (16) generates the following boun- dary equations:

σ ( h , t ) = σ ( 0, t − θ ) + ρ c [ v ( h , t ) − v ( 0, t − θ ) ] + e ϵ [ D ( t − θ ) − D ( t ) ] − ∫ 0 h b ( ξ , t − θ + ξ / c ) d ξ , (19)

σ ( 0 , t ) = σ ( h , t − θ ) + ρ c [ v ( h , t − θ ) − v ( 0 , t ) ] + e ϵ [ D ( t − θ ) − D ( t ) ] + ∫ 0 h b ( ξ , t − ξ / c ) d ξ . (20)

It is easy to verify that Equations (17) and (19), though presented in different forms, are equivalent. The same is true for the pair of Equations (18) and (20). Therefore, we shall use the equations in these pairs interchangeably.

We will not work with boundary equations that can be obtained directly from the displacement representation formula (10), since it is computationally more effective to determine at first unknown boundary values of the velocity v ( x , t ) , and then calculate unknown boundary values of the displacement u ( x , t ) by integrating the boundary velocity over time (using also an initial condition for u ( x , t ) ).

We also need to consider boundary values of the expression (13) for the electric potential. It is important to emphasize that two equations obtained from (13) when x tends to h or to 0 are equivalent and, therefore, they are pre- sented below as one equation:

ϕ ( h , t ) − ϕ ( 0, t ) = e ϵ [ u ( h , t ) − u ( 0, t ) ] − h ϵ D ( t ) (21)

The boundary equations presented above will be used in the next section to create an exact time domain solution procedure in the case when nonlinear damper boundary conditions are sprecified.

Suppose that the lower end face of the piezoelectric element is fixed to a non- linear damper. Let F be a damping force acting on the lower end face which is defined by the following nonlinear relationship [

F = − k α | v ( 0 , t ) | α sgn ( v ( 0 , t ) ) . (22)

where k α > 0 is the damping constant, α > 0 is the damping exponent, and sgn ( . ) is the signum function defined for all real numbers (including 0 where its value is also 0). If v ( 0, t ) ≠ 0 , the direction of F is opposite to v ( 0, t ) . The exponent α has a value 1 for a linear damper, but may vary in practice in the interval ( 0,2 ] [

σ ( 0 , t ) = k α A | v ( 0 , t ) | α sgn ( v ( 0 , t ) ) (23)

where A is the lower end face area.

Consider additional assumptions that will be used to get exact solutions for the damper boundary conditions based on the results of the previous section. We suppose that the values of u , σ , b , ϕ , D are defined for − ∞ < t < ∞ . In addition, let us assume henceforth that

u ( x , t ) = 0 , ϕ ( x , t ) = 0 if 0 < x < h , t < 0 (24)

which means, based on (2) and (3), that σ , D and b are also zero inside the piezoelectric body at negative times. The next additional assumption is that

b ( x , t ) = 0 (25)

inside the piezoelectric body at any time in the sense of generalized functions. This also includes the assumption that the initial conditions for the elastic displacement u ( x , t ) are zero, as discussed in [

Regarding the design of the piezoelectric element, we assume that it is a cylinder (or a rod) with two coated electrodes at x = 0 and x = h . The elec- trodes are considered to be of negligible thickness (from the mechanical point of view) and their deformation is neglected. The output voltage, which is defined as the electric potential difference between the lower and upper electrodes Δ ϕ = ϕ ( 0 , t ) − ϕ ( h , t ) , is of primary interest below.

The electric boundary condition at x = 0 corresponds to the grounded electrode:

ϕ ( 0 , t ) = 0 if t ≥ 0 (26)

At the upper end face, the following mechanical boundary condition is used:

σ ( h , t ) = p ( t ) if t ≥ 0 (27)

where p ( t ) is an applied normal stress load which is assumed to be known and negative.

The electric boundary condition at the upper end face x = h is formulated as follows:

D ( h , t ) = 0 if t ≥ 0. (28)

So, based on (9), D ( t ) = 0 .

Using the above assumptions the representation formulas (14) and (16) for the velocity and stress take the following simplified forms:

v ( x , t ) = 1 2 [ v ( 0 , t − x / c ) + v ( h , t − ( h − x ) / c ) ] + 1 2 ρ c [ p ( t − ( h − x ) / c ) − σ ( 0 , t − x / c ) ] in Ω , (29)

σ ( x , t ) = ρ c 2 [ v ( 0 , t − x / c ) + v ( h , t − ( h − x ) / c ) ] + 1 2 [ σ ( 0 , t − x / c ) + p ( t − ( h − x ) / c ) ] in Ω (30)

where all the time dependent functions are equal to zero for negative times.

In the representation formulas (29) and (30), there are three unknown boundary functions v ( 0, t ) , σ ( 0, t ) and v ( h , t ) first two of which are related by Equation (23). Two additional equations needed for determination of these three functions will be derived below based on (17) and (18).

After the the velocity v ( x , t ) is determined for any particular x over time, the corresponding displacement u ( x , t ) can be obtained (due to the zero initial conditions) as

u ( x , t ) = ∫ 0 t v ( x , τ ) d τ . (31)

Boundary values of the displacement provide (according to (21) and (26)) the electric potential value at x = h :

ϕ ( h , t ) = e ϵ [ u ( h , t ) − u ( 0, t ) ] (32)

The solution of the above problem will be obtained by using an exact recursive procedure based on the following equations obtained from (17) and (18) under the boundary conditions (23) (26) (27) (28):

v ( h , t ) = 2 v ( 0 , t − θ ) − v ( h , t − 2 θ ) + 1 ρ c [ p ( t ) − p ( t − 2 θ ) ] , (33)

v ( 0 , t ) + k α ρ c A | v ( 0 , t ) | α sgn ( v ( 0 , t ) ) = v ( h , t − θ ) + 1 ρ c p ( t − θ ) . (34)

There are two unknowns v ( h , t ) and v ( 0, t ) at each time moment t in these equations. The right-hand sides of the equations are known at each time point since they contain either p ( t ) or time-dalayed function values at t − θ that had to be determined at a previous step of the recursive process.

In order to simplify deriving next results, we need to introduce some addi- tional notations:

γ = k α ρ c A , ξ = v ( 0 , t ) , r = v ( h , t − θ ) + 1 ρ c p ( t − θ ) . (35)

Then, Equation (34) reads as

ξ + γ | ξ | α sgn ( ξ ) = r . (36)

Let the left-hand side of Equation (36) be denoted by f ( ξ ) . Since α > 0 , f ( ξ ) is a continuous strictly monotonically increasing function on ( − ∞ , ∞ ) ranging from − ∞ to ∞ . Therefore, for any real r , there exists one and only one solution of Equation (36) in ( − ∞ , ∞ ) .

Denote by Q α the operator that tranforms r into this solution of equation (36). Thus, Q α is the left inverse operator of the nonlinear operator acting on ξ in the left-hand side of Equation (36). If α = 2 , 1, 1 / 2 or 1 / 3 , the corres- ponding expressions of Q α r are very simple for computations:

{ Q 2 r = 1 2 γ ( − 1 + 1 + 4 | r | γ ) sgn ( r ) , Q 1 r = r 1 + γ , Q 1 / 2 r = 1 4 ( − γ + γ 2 + 4 | r | ) 2 sgn ( r ) , Q 1 / 3 r = − γ [ 1 6 ( 108 r + 12 12 γ 3 + 81 r 2 ) 1 / 3 − 2 γ ( 108 r + 12 12 γ 3 + 81 r 2 ) 1 / 3 ] + r . (37)

The calculation of Q α r for other values of α can effectively be imple- mented using a symbolic computation software like Maple [

With help of the inverse operator Q α Equation (34) can be rewritten in the following explicit form for calculating v ( 0, t ) :

v ( 0 , t ) = Q α [ v ( h , t − θ ) + 1 ρ c p ( t − θ ) ] . (38)

Equation (38) combined with (33) creates the recursive procedure that can be used directly for calculations or can lead to building explicit exact solution for vector { v ( 0, t ) , v ( h , t ) } step by step over consecutive time intervals j θ ≤ t < ( j + 1 ) θ ( j = 0 , 1 , 2 , ⋯ ) . In doing so, it is helpful to substitute v ( 0, t − θ ) in (33) by its expression obtained from (38) which provides the following recursive equation for v ( h , t ) :

v ( h , t ) = 2 Q α [ v ( h , t − 2 θ ) + 1 ρ c p ( t − 2 θ ) ] − [ v ( h , t − 2 θ ) + 1 ρ c p ( t − 2 θ ) ] + 1 ρ c p ( t ) ,

or, using the identity operator I (that leaves unchanged the element on which it operates),

v ( h , t ) = ( 2 Q α − I ) [ v ( h , t − 2 θ ) + 1 ρ c p ( t − 2 θ ) ] + 1 ρ c p ( t ) . (39)

Now we derive some explicit exact solutions for v ( h , t ) and v ( 0, t ) corres- ponding to three practically important ranges of the duration t 1 of the stress load at x = h . Our goal is to present the boundary velocities directly through the transient stress load at x = h that generates the dynamic process in the piezoelectric body.

So, p ( t ) = 0 if t ∉ [ 0 , 2 θ ) . Using the recursive Equation (39) under this condition for consecutive intervals [ 2 k θ , 2 ( k + 1 ) θ ) , k = 0 , 1 , 2 , ⋯ , we finally obtain the following explicit expression for v ( h , t ) :

v ( h , t ) = { 1 ρ c p ( t ) if 0 ≤ t < 2 θ , ( 2 Q α − I ) k [ 2 ρ c p ( t − 2 k θ ) ] if 2 k θ ≤ t < 2 ( k + 1 ) θ , k = 1 , 2 , ⋯ . (40)

Substituting (40) into (38) we get the corresponding explicit expression for v ( 0, t ) :

v ( 0 , t ) = { 0 if 0 ≤ t < θ , Q α ( 2 Q α − I ) k − 1 [ 2 ρ c p ( t − 2 k θ ) ] if ( 2 k − 1 ) θ ≤ t < ( 2 k + 1 ) θ , k = 1 , 2 , ⋯ . (41)

In this case, p ( t ) = 0 if t ∉ [ 0 , 4 θ ) . Acting similarly to section 4 we derive the following explicit exact solutions:

v ( h , t ) = { 1 ρ c p ( t ) if 0 ≤ t < 2 θ , ( 2 Q α − I ) k − 1 [ 1 ρ c p ( t − 2 ( k − 1 ) θ ) + ( 2 Q α − I ) ( 2 ρ c p ( t − 2 k θ ) ) ] if 2 k θ ≤ t < 2 ( k + 1 ) θ , k = 1 , 2 , ⋯ ; (42)

v ( 0 , t ) = { 0 if 0 ≤ t < θ , Q α [ 2 ρ c p ( t − θ ) ] if θ ≤ t < 3 θ , Q α ( 2 Q α − I ) k − 1 [ 1 ρ c p ( t − 2 k θ + θ ) + ( 2 Q α − I ) ( 2 ρ c p ( t − 2 k θ − θ ) ) ] if ( 2 k − 1 ) θ ≤ t < ( 2 k + 1 ) θ , k = 2 , 3 , ⋯ . (43)

So, p ( t ) = 0 if t ∉ [ 0 , 6 θ ) , and the corresponding exact solutions have the following closed form:

v ( h , t ) = { 1 ρ c p ( t ) if 0 ≤ t < 2 θ , 1 ρ c p ( t ) + ( 2 Q α − I ) [ 2 ρ c p ( t − 2 θ ) ] if 2 θ ≤ t < 4 θ , ( 2 Q α − I ) k − 1 [ 1 ρ c p ( t − 2 ( k − 1 ) θ ) + ( 2 Q α − I ) [ 2 ρ c p ( t − 2 k θ ) + ( 2 Q α − I ) ( 2 ρ c p ( t − 2 ( k + 1 ) θ ) ) ] ] if 2 ( k + 1 ) θ ≤ t < 2 ( k + 2 ) θ , k = 1 , 2 , ⋯ ; (44)

v ( 0 , t ) = { 0 if 0 ≤ t < θ , Q α [ 2 ρ c p ( t − θ ) ] if θ ≤ t < 3 θ , Q α [ 1 ρ c p ( t − θ ) + ( 2 Q α − I ) ( 2 ρ c p ( t − 3 θ ) ) ] if 3 θ ≤ t < 5 θ , Q α ( 2 Q α − I ) k − 1 [ 1 ρ c p ( t − 2 k θ + θ ) + ( 2 Q α − I ) [ 2 ρ c p ( t − 2 k θ − θ ) + ( 2 Q α − I ) ( 2 ρ c p ( t − 2 k θ − 3 θ ) ) ] ] if ( 2 k + 3 ) θ ≤ t < ( 2 k + 5 ) θ , k = 1 , 2 , ⋯ . (45)

Similar explicit formulas for t 1 ≥ 6 θ are excessively cumbersome. In this case, it is easier to directly use the recursive procedure based on (33) and (38) which has the same simple form regardless of the transient load duration and also provides exact results.

Consider some examples of using the results of the previous section for mathematical modeling of piezoelectric cylindrical devices installed in a car as proposed in [

C = 5.32 * 10 10 N / m 2 ; e = 19.89 N / ( V ⋅ m ) ; ϵ = 76.12 * 10 − 10 farad / m ; ρ = 7750.0 kg / m 3

Then the elastic wave speed c in the piezoelectric material is equal to 3684.06 m / s . Next, we take into account that the total force instantaneously applied to the top of a piston in an internal combustion engine is around 6300 pounds, which corresponds to approximately 28,640 N [

p ( t ) = − p a [ H ( t ) − H ( t − t 1 ) ] (46)

where t 1 is the duration of the pressure impulse.

We assume first that α = 0.5 and k α = 1000 N ⋅ s 1 / 2 ⋅ m − 1 / 2 (see, e.g., [

Comparing the graphs we can see that the maximum or peak value does not depend on the pressure impulse duration t 1 . However, the number of peaks in each figure depends on the t 1 . The time distance between two neighboring peaks is approximately equal to 2 θ . After the pressure load is removed, there is an attenuation of the output voltage vibrations.

Now let us consider another set of the damper parameters: α = 2.0 and k α = 250 N ⋅ s 2 ⋅ m − 2 . The parameters of the material and the impulse durations are the same as above. The operator Q α is also calculated according to (37).

The calculated results for the output voltage V = Δ ϕ (in kV) are presented for t ≤ 50 μs in Figures 4-6.

Comparison of these graphs shows that the maximum value of the output voltage does not depend on the pressure impulse duration t 1 which similar to the case when α = 0.5 . The difference is that now there is only one peak but its width depends on the t 1 . After the pressure load is removed, the attenuation of the output voltage vibrations is very pronounced: the amplitude of vibrations after the load removal is almost negligible.

One-dimensional transient dynamic piezoelectric problems for thickness polarized layers and disks, or length polarized rods, are considered here in the framework of a time-domain Green’s function method. As the result, a novel exact analytical recursive procedure is derived which is applicable for a wide variety of boundary conditions including the nonlinear damper case. Some new practically important explicit exact solutions are presented. The effectiveness of the proposed exact approach is demonstrated by examples of the time behavior of the output electric potential difference between two electrodes coated at the

end faces of a piezoelectric cylinder fixed to a nonlinear damper at one end, and subjected to impulsive loading at the other.

The authors would like to thank the anonymous reviewers for their comments.

Khutoryansky, N.M. and Genis, V. (2017) Exact Time Do- main Solutions of 1-D Transient Dynamic Piezoelectric Problems with Nonlinear Dam- per Boundary Conditions. Journal of Applied Mathematics and Physics, 5, 873-888. https://doi.org/10.4236/jamp.2017.54077