Post-buckling analysis of viscoelastic plates with fractional derivative models

The post-buckling response of thin plates made of linear viscoelastic materials is investigated. The employed viscoelastic material is described with fractional order time derivatives. The governing equations, which are derived by considering the equilibrium of the plate element, are three coupled nonlinear fractional partial evolution type differential equations in terms of three displacements. The nonlinearity is due to nonlinear kinematic relations based on the von Kármán assumption. The solution is achieved using the analog equation method (AEM), which transforms the original equations into three uncoupled linear equations, namely a linear plate (biharmonic) equation for the transverse deflection and two linear membrane (Poisson’s) equations for the inplane deformation under fictitious loads. The resulting initial value problem for the fictitious sources is a system of nonlinear fractional ordinary differential equations, which is solved using the numerical method developed recently by Katsikadelis for multi-term nonlinear fractional differential equations. The numerical examples not only demonstrate the efficiency and validate the accuracy of the solution procedure, but also give a better insight into this complicated but very interesting engineering plate problem     

A nonlinear fractional viscoelastic material model for polymers

In this contribution a test scheme based on tensile tests at different velocities, relaxation experiments and deformation controlled loading and unloading processes with intermediate relaxations has been used to experimentally characterize the nonlinear, inelastic material behavior. Based on the experimental observations a small strain nonlinear fractional viscoelastic material model is derived. In order to use the model within a finite element analysis, the constitutive equations have been generalized for the multiaxial case. The experimental test scheme and the fractional viscoelastic material model are subsequently applied to characterize and compute the mechanical behavior of the thermoplastic Polypropylene. After the identification of the material parameters several uniaxial and multiaxial simulations have been carried out and compared with experimental results.     

Application of fractional derivative models in linear viscoelastic problems

Appropriate knowledge of viscoelastic properties of polymers and elastomers is of fundamental importance for a correct modelization and analysis of structures where such materials are present, especially when dealing with dynamic and vibration problems. In this paper experimental results of a series of compression and tension tests on specimens of styrene-butadiene rubber and polypropylene plastic are presented; tests consist of creep and relaxation tests, as well as cyclic loading at different frequencies. Experimental data are then used to calibrate some linear viscoelastic models; besides the classical approach based on a combination in series or parallel of standard mechanical elements as springs and dashpots, particular emphasis is given to the application of models whose constitutive equations are based on differential equations of fractional order (Fractional Derivative Model). The two approaches are compared analyzing their capability to reproduce all the experimental data for given materials; also, the main computational issues related with these models are addressed, and the advantage of using a limited number of parameters is demonstrated.     

Universal linear viscoelastic approximation property of fractional viscoelastic models with application to asphalt concrete

The paper presents a comprehensive linear viscoelastic characterization of asphalt concrete using fractional viscoelastic models. For this purpose, it is shown that fractional viscoelastic models are universal approximators of relaxation and retardation spectra. This essentially means that any spectrum can be mathematically represented by fractional viscoelastic models. Characterization of asphalt concrete is performed by constructing the dynamic modulus master curve and determining the parameters of the generalized fractional Maxwell model (GFMM). This procedure is similar to the widely used one of determining the master curve of asphalt concrete using a statistical function such as the sigmoidal model. However, from the GFMM, the relaxation modulus, creep compliance, continuous relaxation spectrum, and Prony series parameters can be determined analytically. A further advantage of the GFMM is that unlike the sigmoidal model, which only gives a representation of either the dynamic modulus or the storage modulus, the GFMM gives a representation of both the storage modulus and loss modulus (and therefore also the dynamic modulus and phase angle). The procedure was successfully applied to ten different mixes used in the State of Virginia.     

On the dynamic stochastic response of FE models

A numerical procedure to compute the mean and covariance matrix of the response of nonlinear structures modeled by large FE models is presented. Non-white, non-zero mean, non-stationary Gaussian distributed excitation is represented by the well known Karhunen–Loéve expansion, which allows to describe any type of non-white Gaussian excitation in contrast to filtered white noise which might not be easily adjusted to available statistical data. The solution procedure differs considerably from standard methodologies using a state space representation. In the proposed approach, step-by-step integration procedures developed for deterministic FE analysis are applied to compute the first two moments of the stochastic response.     

An analytical and numerical study of a dynamically accelerating semi-infinite crack in a linear viscoelastic material

In a previous paper, a general method was presented for constructing the solution to the problem of a semi-infinite, mode III crack propagating dynamically through an infinite, general linear viscoelastic body. The only restrictions placed upon the crack tip speed were that it have constant sign and in magnitude not exceed the glassy shear wave speed. In the present contribution, those previous analytical results are applied to a study of dynamic unsteady crack growth in a linear viscoelastic body. In particular, a numerical algorithm for computing the stress intensity factor is given along with example simulations of running cracks using the Achenbach-Chao viscoelastic model and a stress intensity factor (SIF) fracture criterion. We also compare the transient SIF with the dynamic steady state SIF, and examine the transition to constant crack speed for a dynamically accelerating crack in a viscoelastic material.     

On the thermodynamically consistent fractional wave equation for viscoelastic solids

Since the well-known methods for the computation of harmonically induced LAMB waves in elastic plates cannot be applied directly to viscoelastic material, an enhanced material model is developed. It is based on fractional time derivatives and incorporates a fractional KELVIN?CVOIGT model. A proof of the thermodynamic consistency of this model is given. On the basis of the developed material model, the fractional wave equation is derived which allows for the calculation of LAMB waves in viscoelastic solids.     

Marginal analysis guided design justification: a material handling example

The design and justification of an integrated manufacturing system and all of the supporting sub-systems (e.g. material handling system) is a difficult task. One of the most difficult aspects in the design process is determining how and where to proceed. This paper proposes the use of a marginal analysis-directed branch and bound approach as the basis for directing the design process so that a combined economic and performance justification path results. The background and rationale for a marginal analysis guided design approach is given and an example of a material handling system design is presented to illustrate the development of a design justification decision path.     

Experimental and numerical analysis of a bus component in composite material

The aim of this work is the fatigue design of a structural component manufactured in composite materials. The use of composite materials in structural applications implies that the reliability and safety requirements are met, in particular from a point of view of the fatigue life. A bus component is considered: the composite material is obtained by means of the pultrusion technique: the longitudinal reinforcement is glass fibre, while the matrix is polyester resin. At first, a static and fatigue characterization of the pultruded composite is performed to identify the mechanical behaviour along the fibres and in the normal direction, following an extensive testing programme on specimens. During the fatigue tests, the variation of the elastic modulus and the residual strength are monitored to characterize the damage; further experimental tests are performed to fit parameters required for models of fatigue life prevision. Fatigue bending tests are performed on the component, showing crack propagation through the section: a check at the SEM is performed to identify the main internal type of damage. A final comparison with a numerical simulation is proposed: this numerical model will be useful for the optimal fatigue design of the section.     

On the numerical integration of a class of pressure-dependent plasticity models

An unconditionally stable algorithm for the numerical integration of elastoplastic pressure-dependent constitutive relations is analysed in detail in this paper. The application of the method to plane stress problems, in which the out-of-plane strain component is not defined kinematically, is discussed. The tangent moduli resulting from this integration algorithm are obtained by consistent linearization of the elastoplastic constitutive equations. The algorithm is applied to Gurson's constitutive model, some one-dimensional problems are solved, and comparisons with exact solutions are made. The paper closes with a numerical study of the necking of an axi-symmetric specimen using Gurson's plasticity model to describe the constitutive behaviour of the material.     

Experiments on crack propagation in a viscoelastic material

The problem of plane stress crack growth in a rectangular linearly viscoelastic plate under a pair of equal and opposite forces applied to the crack surface is studied experimentally. A fatigue crack is generated at the end of a slot cut on the specimen. After application of two opposite forces to the specimen with a fatigue crack, the crack grows gradually. It is found that under sufficiently small load, the speed of crack growth decreases as time increases, and crack arrest occurs finally. However when the applied load is sufficiently large, the crack grows with increasing speed until a catastrophic failure of the specimen takes place. Between these two extreme cases, the crack grows with an oscillatory speed. A theory based on the models of an infinite plate and double cantilevers is proposed for the explanation of the behaviour of crack growth under various loads. Our experimental results check fairly well with theoretical predictions.     

Finite element analysis of viscoelastic composite structures based on a micromechanical material model

The stress and creep analysis of structures made of micro-heterogeneous composite materials is treated as a two-scale problem, defined as a mechanical investigation on different length scales. Reinforced composites show by definition a heterogeneous texture on the microlevel, determined by the constitutive behaviour of the matrix material and the embedded fibres as well as the characteristics of the bonding properties in the interphase. All these heterogeneities are neglected by the finite element analysis of structural elements on the macroscale, since a ficticious and homogeneous continuum with averaged properties is assumed. Therefore, the constitutive equations of the substitute material should well reflect the mechanical behaviour of the existing micro-heterogeneous composite in an average sense.The paper at hand starts with the brief outline of a micromechanical model, named generalized method of cells (GMC), which provides the macrostress responses due to macrostrain processes as well as the homogenised constitutive tensor of the substitute material. The macroscopic stresses and strains are obtained as volume averages of the corresponding microfields within a representative volume element. The effective material tensor constitutes the mapping between the macro-strains and the macro-stresses. The cells method is used for the homogenisation of the unidirectionally reinforced single layers of laminates made of viscoelastic resins and flexibly embedded elastic fibres. The algorithm for the homogenisation of the constitutive properties runs simultaneously to the finite element analysis at each point of numerical integration and provides the macro-stresses and the homogenised constitutive properties. The validity of the proposed two-scale simulation is investigated by solving boundary value problems and comparing the numerical results for the structures to the experimental data of creep and relaxation tests or analytical solutions.     

Determination of the relaxation modulus of a linearly viscoelastic material

In this paper, a numerical method for computing the relaxation modulus of a linearly viscoelastic material is presented. The method is valid for relaxation tests where a constant strain rate is followed by a constant strain. The method is similar to the procedure suggested by Zapas and Phillips. Unlike Zapas-Phillips approach, this new method can be also applied for times shorter than t ₁/2, where t ₁ denotes time when the maximum strain is achieved. Therefore this method is very suitable for materials that experiences fast relaxation. The method is verified with numerical simulations. Results from the simulations are compared with analytical solution and Zapas-Phillips method. Results indicate that the presented approach is suitable for estimating the relaxation modulus.     

Application of fractional derivatives to the analysis of damped vibrations of viscoelastic single mass systems

Summary Free damped vibrations of an oscillator, whose viscoelastic properties are described in terms of the fractional calculus Kelvin-Voight model, Maxwell model, and standard linear solid model are determined. The problem is solved by the Laplace transform method. When passing from image to pre-image one is led to find the roots of an algebraic equation with fractional exponents. The method for solving such equations is proposed which allows one to investigate the roots behaviour in a wide range of single-mass system parameters. A comparison between the results obtained on the basis of the three models has been carried out. It has been shown that for all models the characteristic equations do not possess real roots, but have one pair of complex conjugates, i.e. the test single-mass systems subjected to the impulse excitation do not pass into an aperiodic regime in none of magnitudes of the relaxation and creep times. Main characteristics of vibratory motions of the single-mass system as functions of the relaxation time or creep time, which are equivalent to the temperature dependencies, are constructed and analyzed for all three models.