Determination of the rheological properties of thin plate under transient vibration

Determination of the rheological properties of thin plate under transient vibration

Jan Sko č ilas^I; Blanka Sko č ilasová^II; Josef Soukup^III

^IJan Skočilas, M.Sc., Ph.D.: Czech Technical University in Prague, Faculty of Mechanical Engineering, Department of Process Engineering, Technická 4, 166 07 Praha, Czech Republic ([e-mail])

^IIBlanka Skočilasová, M.SC. (Ph.D. student): Jan Evangelista Purkyně University in Ústí nad Labem, Faculty of Production Technology and Management, Department of Machines and Mechanics, Na Okraji 1001, 400 01 Ústí nad Labem, Czech Republic ([e-mail])

^IIIJosef Soukup, assoc. prof., ing., CSc.: Jan Evangelista Purkyně University in Ústí nad Labem, Faculty of Production Technology and Management, Department of Machines and Mechanics, Na Okraji 1001, 400 01 Ústí nad Labem, Czech Republic ([e-mail])

ABSTRACT

The article deals with systematic analysis of the transient vibration of rectangular viscoelastic orthotropic thin 2D plate. The analysis is focused on specific deformation models of plate and for specific linear models of rheological properties. The viscoelastic isotropic and anisotropic (orthotropic) materials were investigated. The plate is loaded by general transient impulse. The time and coordinate dependencies of the fundamental quantities - displacements, rotations, stress and deformations in arbitrary points of plate under transient vibration, e. g. the analysis of stress and deformation waves in plates were investigated. The selected models are defined by constitutive equations for stress and deformation dependence. The isotropic and orthotropic model is considered. The analysis results depend on quality of the determined mechanical properties of the materials. The standard methods commonly used are time-consuming and not accurate enough. The new methodology of determining the material parameters directly from investigated plate is proposed and proven.

Keywords: rheologic properties, thin plate, transient vibration

1 INTRODUCTION

The theory of the plate vibration belongs to the most explored realm of mechanics of deformable plates. Authors have dealt with this realm more than 200 years. The hundred studies and monographs have been published till present time. The results of vibration investigation were obtained by many analytical, numerical and experimental methods. The most sophisticated analytical solution is focused on transient vibrations of thin elastic isotropic plate. The anisotropic plate e.g. orthotropic plate is less commonly discussed, e.g. [2]. The investigation of transient vibration of viscoelastic plate appears in the first half of the last century and this problem has been examined by significantly less articles [7]. The essential evolution of transient vibration of the viscoelastic plate begun after technology development of polymers (macromolecular materials) [6] and composites until second half of last century, first for Kirchhoff's model of plate [3, 4], later for Rayleigh model. The investigation of the viscoelastic anisotropic plate is progressing realm of mechanics now.

The article deals with the investigation of transient vibrations and stress wave propagation in viscoelastic orthotropic thin plate. The phenomenological models are commonly applied for description of material mechanical properties under transient states of viscoelastic bodies. The constitutive equations define relations of stress and deformation of investigation continuum by superposition of elastic and viscoelastic components (represented by spring and viscous damper respectively). The different models of viscoelastic heredity material were developed by various combinations (parallel or serial) of these elements, e.g. Voight-Kelvin, Maxwell, Zener etc [5].

The general problem is determination of the stress modulus of the viscoelastic orthotropic material (e.g. composite). The contemporary methods of stress modulus determination are based on static measurement of the modulus of individual elements. The final modulus of model is determined from conversion of measured modulus by volume fractions. This method is not satisfactory accurate and the results obtained for dynamic loading of the model are disputable.

These materials are wide used, e.g. in automotive industry, rail vehicles, planes, ships, space research, civil engineering etc. The composites and laminated materials have indisputable preferences (lower weight, higher strength, time stability etc.). It is obvious that the high requirements will be demanded to computational methods for these materials applications in the future. The presented contribution deals with one of the important assumptions for accuracy improvement of these calculations and determination of elasticity modulus of these materials under dynamic loading. The proposed method would be applicable for composite and sandwich construction and laminates.

2 SOLUTION

Most of the commonly used works for investigation of stress analysis in linear viscoelastic bodies solve constitutive equations in differential form. The features of these linear models for one axis tension are described in Alfrey [1} by equation

where p ≥ 0, q ≥ 0 are material parameters, ε- deformation, σ- stress, n - number of material parameters p, m - number of material parameters q.

If m= n then it is trivial case

When m ≠ n, the important cases are for m= (n - 1)or m= (n+1). These models correspond to fundamental phenomenological models of viscous continuum. It is necessary to experimentally obtain the values of (m+n)material parameters, when these models are applied. The differential equations of constitutive equations are more complicated under multi-axis tension.

For two-axis tension (normal stress σ_x, σ_y, in axis direction, and shear stress σ_xy), the constitutive equations for hereditary model of viscoelastic material could be expressed by equation that has analogous form to Alfrey's form for one-axis tension.

Normal stress

for i, j ε {x, y} and i≠j.

Shear stress

The constitutive equation could be expressed in form by technical parameters (considered basic simplest models).

Voigt-Kelvin model Normal stress for isotropic material

Shear stress

where E, G are the Young and Kirchhoff modules, λ, ή are the coefficients of normal and tangential (shear) viscosity, μ, ν are the elastic and viscous Poisson ratios.

Maxwell model Normal stress for isotropic material

Shear stress

It is possible to obtain the equations for stress components in integral form by solution of differential equation (5a) for σ_i, or (5c) for σ_xyrespectively.

Normal stress for isotropic material

After integration by parts, we got

Normal stress for orthotropic material

After integration by parts

Shear stress

Zener model Normal stress for isotropic material

Normal stress for orthotropic material

Shear stress

We get the equations for stress components in integral form by solution of differential equation (7b) for σ_i, for orthotropic material and zero initial conditions.

The specific material parameters E, G, μ, λ, ή are necessary to determine by experiments. More often it is determined from simple bodies with static experiments and then it is applied to bodies with more complicated shapes and complicated loading process. The number of parameters arises from particular cases mentioned above.

The relationship between stress and deformation for hereditary viscoelastic materials is possible to express by Duhamel integral. Integral operator for one-axis tension expresses the fundamental constitutive equation.

The relationship of both expressions in the differential (Alfrey) and integral (Volterra) form is possible to formulate after Laplace transformation (s - parameter of transformation)

where R(s).J(s) = s^-2, and R(t-ז).J(ז) = t.

The Boltzmann principle of superposition allows to write by adjustment of (1) for p₀ = 1 substituting into (8)

The relationship of both solutions is obvious from comparison of (5a), (7a), (8) and (10).

Conclusion

It is more comfortable and less time-consuming to use constitutive equations in proposed form (approximation of time function, relaxation module K(t - ז)), than to determine of many parameters.

From equation (10) arise the necessity of determination of parameter q₀ and function K(t)=Ae^-ςt, which approximate function of viscous component relaxation (module), parameters A and ς. The approximation of this function is less difficult than determination of material parameters (E, G, μ, λ, η).

In case of transient states investigation for thin viscoelastic plate vibration is possible to derive equations for stress components or displacement (time dependent) by integral expression of constitutive equations (8) and (10) respectively. The required parameters could be obtained from comparison of stress or displacement components find out by experimental and analytical methods.

The relaxation module function could be experimentally determined directly on investigated body. The methodology for experimental investigation of this function was proposed and it is under patent application.

Acknowledgement

This work has been supported by grant IGA - Internal Agency of University of Jan Evangelista Purkyně in Ústí nad Labem and grand GA ČR, number of project 101/07/0946.

• [1] Alfrey, T.: Non-homogeneous stresses in visco-elastic media Q. Appl. Math. 2, pp. 113-119, 1944
• [2] Ghugal, Y., M., Sayyad, A., S.: Free vibration of thick orthotropic plates using trigonometric shear deformation theory. pp. 229-244, Vol. 8, No. 3/2009, Latin American Journal of Solids and Structures.
• [3] Kolsky, H.: The propagatin of stress waves in viscoelastic solids. J. Appl. Mechanics Revies, Vol. 11, No 9, pp. 465-468, 1958.
• [4] Leitman, M., J.: The linear theory of viscoelasticity. Encyklopedia of Physics, ed. S. Flügge, Vol. Via/3 Mechanics of solids III., ed. C. Trusdell, Spinger, Berlin, 1973.
• [5] Macosko, Ch., W.: Rheology. Principles, Measurements, and Applications. Wiley-VCH, New York, 1994.
• [6] Riande, E., et al: Polymer Viscoelasticity Marcel Dekker, N.Y.- Basel, 2000
• [7] Tschoegel, N., W.: The Phenomenological Tudory of Linear Viscoelastic Behavior Springer, Berlin, 1989

Publication Dates

• Publication in this collection
  25 Feb 2013
• Date of issue
  Jan 2013