Finite element formulation of viscoelastic sandwich beams using fractional derivative operators

This paper presents a finite element formulation for transient dynamic analysis of sandwich beams with embedded viscoelastic material using fractional derivative constitutive equations. The sandwich configuration is composed of a viscoelastic core (based on Timoshenko theory) sandwiched between elastic faces (based on Euler–Bernoulli assumptions). The viscoelastic model used to describe the behavior of the core is a four-parameter fractional derivative model. Concerning the parameter identification, a strategy to estimate the fractional order of the time derivative and the relaxation time is outlined. Curve-fitting aspects are focused, showing a good agreement with experimental data. In order to implement the viscoelastic model into the finite element formulation, the Grünwald definition of the fractional operator is employed. To solve the equation of motion, a direct time integration method based on the implicit Newmark scheme is used. One of the particularities of the proposed algorithm lies in the storage of displacement history only, reducing considerably the numerical efforts related to the non-locality of fractional operators. After validations, numerical applications are presented in order to analyze truncation effects (fading memory phenomena) and solution convergence aspects.     

Transient dynamic analysis of tapered FRP composite transmission poles using finite element method

Power transmission poles are subjected to dynamic cantilever bending due to wind gusts and cable unilateral failure, or may also be subjected to vehicle impacts. In this paper, transient dynamic analysis of tapered fiber-reinforced polymer (FRP) composite transmission poles with circular thin-walled cross-section subjected to dynamic cable tension and vehicle impacts is investigated by combination of tapered beam finite element and precise time integration method. It is assumed that the material behavior is linearly elastic and the laminate of the cross-section of the wall is symmetric or antisymmetric angle-ply. The effect of fiber type and orientation, the pole geometry, and the concentrated mass at the pole tip are evaluated by performing the dynamic analysis of FRP poles under step, triangular and sine pulses. There is a good agreement between the results of the present method and those obtained from the poles modeled by ANSYS commercial finite element software and existed literatures. Also, there is a significant shorter run-time in the present method. It is concluded that beyond 10 layers for the laminate with constant thickness of the wall, the pole tip deflection does not change.     

The effect of loading time on flexible pavement dynamic response: a finite element analysis

Dynamic response of asphalt concrete (AC) pavements under moving load is a key component for accurate prediction of flexible pavement performance. The time and temperature dependency of AC materials calls for utilizing advanced material characterization and mechanistic theories, such as viscoelasticity and stress/strain analysis. In layered elastic analysis, as implemented in the new Mechanistic-Empirical Pavement Design Guide (MEPDG), the time dependency is accounted for by calculating the loading times at different AC layer depths. In this study, the time effect on pavement response was evaluated by means of the concept of “pseudo temperature.” With the pavement temperature measured from instrumented thermocouples, the time and temperature dependency of AC materials was integrated into one single factor, termed “effective temperature.” Via this effective temperature, pavement responses under a transient load were predicted through finite element analysis. In the finite element model, viscoelastic behavior of AC materials was characterized through relaxation moduli, while the layers with unbound granular material were assumed to be in an elastic mode. The analysis was conducted for two different AC mixtures in a simplified flexible pavement structure at two different seasons. Finite element analysis results reveal that the loading time has a more pronounced impact on pavement response in the summer for both asphalt types. The results indicate that for reasonable prediction of dynamic response in flexible pavements, the effect of the depth-dependent loading time on pavement temperature should be considered.     

Formulation of reference surface element and its applications in dynamic analysis of delaminated composite beams

This paper presents a new finite element formulation, referred to as reference surface element (RSE) model, for numerical prediction of dynamic behaviour of delaminated composite beams and plates using the finite element method. The RSE formulation can be readily incorporated into all elements based on the Timoshenko beam theory and the Reissner–Mindlin plate theory taking into account the transverse shear deformations. The ‘free model' and ‘constrained model' for dynamic analysis of delaminated composite beams and/or plates have been unified in this RSE formulation. The RSE formulation has been applied to an existing 2-node Timoshenko beam element taking into account the transverse shear deformations and the bending–extension coupling. Frequencies and vibration mode shapes are determined through solving an eigenvalue problem. Numerical results show that the present RSE model is reliable and practical when used to predict frequencies and mode shapes of delaminated composite beams. The RSE formulation has also been used to investigate the effects of the number, size and interfacial loci of delaminations on frequencies and mode shapes of composite beams.     

Space adaptive finite element methods for dynamic Signorini problems

Space adaptive techniques for dynamic Signorini problems are discussed. For discretisation, the Newmark method in time and low order finite elements in space are used. For the global discretisation error in space, an a posteriori error estimate is derived on the basis of the semi-discrete problem in mixed form. This approach relies on an auxiliary problem, which takes the form of a variational equation. An adaptive method based on the estimate is applied to improve the finite element approximation. Numerical results illustrate the performance of the presented method. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Finite element/finite volume approaches with adaptive time stepping strategies for transient thermal problems

Transient thermal analysis of engineering materials and structures by space discretization techniques such as the finite element method (FEM) or finite volume method (FVM) lead to a system of parabolic ordinary differential equations in time. These semidiscrete equations are traditionally solved using the generalized trapezoidal family of time integration algorithms which uses a constant single time step. This single time step is normally selected based on the stability and accuracy criteria of the time integration method employed. For long duration transient analysis and/or when severe time step restrictions as in nonlinear problems prohibit the use of taking a larger time step, a single time stepping strategy for the thermal analysis may not be optimal during the entire temporal analysis. As a consequence, an adaptive time stepping strategy which computes the time step based on the local truncation error with a good global error control may be used to obtain optimal time steps for use during the entire analysis. Such an adaptive time stepping approach is described here. Also proposed is an approach for employing combinedFEM/FVM mesh partitionings to achieve numerically improved physical representations. Adaptive time stepping is employed thoughout to practical linear/nonlinear transient engineering problems for studying their effectiveness in finite element and finite volume thermal analysis simulations. Additional support and computing times were furnished by Minnesota Supercomputer Institute at the University of Minnesota.     

Rapid implementation of material models within finite element analysis

This paper presents a simple procedure for obtaining a numerical approximation to the consistent tangent matrix, together with a straightforward implicit (Euler backward) integration algorithm. The combined algorithm is used to incorporate four models into the commercial finite element package ABAQUS/Standard; illustrating how it can be used to rapidly implement material models within finite element analysis. The models have been chosen, not only because they help to illuminate the structure of the algorithm, but also because they illustrate its wide ranging applicability and permit the procedure to be tested against analytical results and an existing, well established, model.     

Hybrid finite element for analysis of functionally graded beams

A hybrid finite element model is presented, where stiffness and mass distributions over a beam with functionally graded material (FGM) are accurately modeled for both elastic and inelastic material responses. Von Mises and Drucker-Prager plasticity models are implemented for metallic and ceramic parts of FGM, respectively. Three-dimensional stress-strain relations are solved by a general closest point projection algorithm, and then condensed to the dimensions of the beam element. Numerical examples and verification studies on a proposed element demonstrate accuracy and robustness under inelastic material response as well as capturing fundamental, higher, and mix modes of vibration frequencies and shapes.     

动态镦粗过程的动力显式有限元分析

为了研究金属的三维动态锤锻成形过程,基于连续介质力学及有限变形理论,建立了一种有限元模型.采用动力分析方法,在运动方程中加入惯性力项考虑锤锻中显著的惯性效应;根据设备的工作原理按照能量守恒定律计算变形期间的锤头速度;同时,将变形视为一个绝热过程计算变形期间试样内部的温度升高.基于建立的模型开发了动力显式有限元分析程序,模拟了铅块试样在落锤打击下的动态镦粗过程,给出了试样内部的位移、等效应变、等效应力和温度分布规律.将变形后试样几何形状、成形载荷-时间曲线和锤头速度-时间曲线的计算结果与实验结果相对比,表明了开发程序计算结果的准确性.     

Elasto-plastic finite element analysis of a crack in an infinite plate

The elastic support method was recently developed to simulate the effects of unbounded solids in the finite element analysis of stresses and displacements. The method eliminates all the computational disadvantages encountered in the use of `infinite' elements or coupled finite element boundary element methods while retaining all the computational advantages of the finite element method. In this paper, the method is extended to the elasto-plastic analysis of fracture in infinite solids by using the load increment approach and including the effects of strain hardening. Numerical tests and parametric study are conducted by analysing a straight crack in an infinite plate. Present results for J integrals and plastified zones are compared, respectively, with analytical solutions and available results obtained by using the body force method. The agreement between the results is found to be very good even if the truncation boundary of the finite element model is located very close to the crack tip or the plastified zone.     

Stress intensity factor error index for finite element analysis with singular elements

An error index for the stress intensity factor (SIF) obtained from the finite element analysis results using singular elements is proposed. The index was developed by considering the facts that the analytical function shape of the crack tip displacement is known and that the SIF can be evaluated from the displacements only. The advantage of the error index is that it has the dimension of the SIF and converges to zero when the actual error of the SIF by displacement correlation technique converges to zero. Numerical examples for some typical crack problems, including a mixed mode crack, whose analytical solutions are known, indicated the validity of the index. The degree of actual SIF error seems to be approximated by the value of the proposed index.     

Iterative mesh partitioning optimization for parallel nonlinear dynamic finite element analysis with direct substructuring

 This work presents a novel iterative approach for mesh partitioning optimization to promote the efficiency of parallel nonlinear dynamic finite element analysis with the direct substructure method, which involves static condensation of substructures' internal degrees of freedom. The proposed approach includes four major phases – initial partitioning, substructure workload prediction, element weights tuning, and partitioning results adjustment. The final three phases are performed iteratively until the workloads among the substructures are balanced reasonably. A substructure workload predictor that considers the sparsity and ordering of the substructure matrix is used in the proposed approach. Several numerical experiments conducted herein reveal that the proposed iterative mesh partitioning optimization often results in a superior workload balance among substructures and reduces the total elapsed time of the corresponding parallel nonlinear dynamic finite element analysis. Received 22 August 2001 / Accepted 20 January 2002     

A frequency-domain approach for modelling transient elastodynamics using scaled boundary finite element method

This study develops a frequency-domain method for modelling general transient linear-elastic dynamic problems using the semi-analytical scaled boundary finite element method (SBFEM). This approach first uses the newly-developed analytical Frobenius solution to the governing equilibrium equation system in the frequency domain to calculate complex frequency-response functions (CFRFs). This is followed by a fast Fourier transform (FFT) of the transient load and a subsequent inverse FFT of the CFRFs to obtain time histories of structural responses. A set of wave propagation and structural dynamics problems, subjected to various load forms such as Heaviside step load, triangular blast load and ramped wind load, are modelled using the new approach. Due to the semi-analytical nature of the SBFEM, each problem is successfully modelled using a very small number of degrees of freedom. The numerical results agree very well with the analytical solutions and the results from detailed finite element analyses.     

Boundary element method analyses of transient heat conduction in an unbounded solid layer containing inclusions

The boundary element method (BEM) is used to compute the three-dimensional transient heat conduction through an unbounded solid layer that may contain heterogeneities, when a pointwise heat source placed at some point in the media is excited. Analytical solutions for the steady-state response of this solid layer when subjected to a spatially sinusoidal harmonic heat line source are presented when the solid layer has no inclusions. These solutions are incorporated into a BEM formulation as Greens functions to avoid the discretization of flat media interfaces. The solution is obtained in the frequency domain, and time responses are computed by applying inverse (Fast) Fourier Transforms. Complex frequencies are used to prevent the aliasing phenomena. The results provided by the proposed Greens functions and BEM formulation are implemented and compared with those computed by a BEM code that uses the Greens functions for an unbounded media which requires the discretization of all solid interfaces with boundary elements. The proposed BEM model is then used to evaluate the temperature field evolution through an unbounded solid layer that contains cylindrical inclusions with different thermal properties, when illuminated by a plane heat source. In this model zero initial conditions are assumed. Different simulation analyses using this model are then performed to evaluate the importance of the thermal properties of the inclusions on transient heat conduction through the solid layer.     

Multiscale analysis of stiffness and fracture of nanoparticle-reinforced composites using micromechanics and global–local finite element models

A combined micromechanics analysis and global–local finite element method is proposed to study the interaction of particles and matrix at the nano-scale near a crack tip. An analytical model is used to obtain the effective elastic modulus of nanoparticle-reinforced composites, then a global–local multi-scale finite element model with effective homogeneous material properties is used to study the fracture of a compact tension sample. For SiO₂ particle-reinforced epoxy composites with various volume fractions, the simulation results for effective elastic modulus, fracture toughness, and critical strain energy release rate show good agreement with previously published experimental data. It is demonstrated that the proposed parametric multi-scale model can be used to efficiently study the toughness mechanisms at both the macro and nano-scale.     

2D finite element analysis of thermally bonded nonwoven materials: Continuous and discontinuous models

Due to a random structure of nonwoven materials, their non-uniform local material properties and nonlinear properties of single fibres, it is difficult to develop a numerical model that adequately accounts for these features and properly describes their performance. Two different finite element (FE) models – continuous and discontinuous – are developed here to describe the tensile behaviour of nonwoven materials. A macro-level continuum finite element model is developed based on the classic composite theory by treating the fibrous network as orthotropic material. This model is used to analyse the effect of thermally bonding points on the deformational behaviour and deformation mechanisms of thermally bonded nonwoven materials at macro-scale. To describe the effects of discontinuous microstructure of the fabric and implement the properties of polypropylene fibres, a micro-level discontinuous finite element model is developed. Applicability of both models to describe various deformational features observed in experiments with a real thermally bonded nonwoven is discussed.     

A finite element model for contact analysis of multiple Cosserat bodies

The objective of this paper is to develop a finite element model for multi-body contact analysis of Cosserat materials. Based on the parametric virtual work principle, a quadratic programming method is developed for finite element analysis of contact problems. The contact problem with friction between two Cosserat bodies is treated in the same way as in plastic analysis. The penalty factors, that are normally introduced into the algorithm for contact analysis, have a direct influence on accuracy of solution. There is no available rule for choosing a reasonable value of these factors for simulation of contact problems of Cosserat materials, and they are therefore cancelled through a special technique so that the numerical results can be of high accuracy. Compared with the conventional work on Cosserat elasticity, the newly developed model is on the contact analysis of the Cosserat materials and is seldom found in the existing literatures. Four examples are computed to illustrate the validity and importance of the model developed.     

A finite element analysis of quasistatic crack growth in a pressure sensitive constrained ductile layer

Polymeric adhesive layers are employed for bonding two components in a wide variety of technological applications. It has been observed that, unlike in metals, the yield behavior of polymers is affected by the state of hydrostatic stress. In this work, the effect of pressure sensitivity of yielding and layer thickness on quasistatic interfacial crack growth in a ductile adhesive layer is investigated. To this end, finite deformation, finite element analyses of a cracked sandwiched layer are carried out under plane strain, small-scale yielding conditions for a wide range of mode mixities. The Drucker–Prager constitutive equations are employed to represent the behavior of the layer. Crack propagation is simulated through a cohesive zone model, in which the interface is assumed to follow a prescribed traction–separation law. The results show that for a given mode mixity, the steady state fracture toughness |K|_ss is enhanced as the degree of pressure sensitivity increases. Further, for a given level of pressure sensitivity, |K|_ss increases steeply as mode II loading is approached.     

Combination of finite element and reliability methods in nonlinear fracture mechanics

This paper presents a probabilistic methodology for nonlinear fracture analysis in order to get decisive help for the reparation and functioning optimization of general cracked structures. It involves nonlinear finite element analysis. Two methods are studied for the coupling of finite element with reliability software: the direct method and the quadratic response surface method. To ensure the response surface efficiency, we introduce new quality measures in the convergence scheme. An example of a cracked pipe is presented to illustrate the proposed methodology. The results show that the methodology is able to give accurate probabilistic characterization of the J-integral in elastic–plastic fracture mechanics without obvious time consumption. By introducing an “analysis re-using” technique, we show how the response surface method becomes cost attractive in case of incremental finite element analysis.