• Journal of Mechanical Science and Technology
• /
• v.14 no.1
• /
• pp.19-29
• /
• 2000

A viscoelastic finite element analysis is presented to investigate residual stresses occurred in a laminated cylindrical shell during cure. An incremental viscoelastic constitutive equation that can describe stress relaxation during the cure is derived as a recursive formula which can be used conveniently for a numerical analysis. The finite element analysis program is developed on the basis of a 3-D degenerated shell element and the first order shear deformation theory, and is verified by comparing with an one dimensional exact solution. Viscoelastic effect on the residual stresses in the laminated shell during the cure is investigated by performing both the viscoelastic and linear elastic analyses considering thermal deformation and chemical shrinkage simultaneously. The results show that there is big difference between viscoelastic stresses and linear elastic stresses. The effect of cooling rates and cooling paths on the residual stresses is also examined.

• International Journal of Highway Engineering
• /
• v.14 no.5
• /
• pp.1-9
• /
• 2012

PURPOSES: A viscoelastic axisymmetric finite element analysis code has been developed for stress analysis of asphalt pavement structures. METHODS: Generalized Maxwell Model (GMM) and 4-node isoparametric element were employed for finite element formulation. The code was developed using $C^{+}^{+}$ computer program language and named as KICTPAVE. For the verification of the developed code, a structural model of a pavement system was constructed. The structural model was composed of three layers: asphalt layer, crushed stone layer, and soil subgrade. Two types of analysis were considered for the verification: (1)elastic static analysis, (2)viscoelastic time-dependent analysis. For the elastic static analysis, linear elastic material model was assigned to all the layers, and a static load was applied to the structural model. For the viscoelastic time-dependent analysis, GMM and linear elastic material model were assigned to the asphalt layer and all the other layers respectively, and a cyclic loading condition was applied to the structural model. RESULTS: The stresses and deformations from KICTPAVE were compared with those from ABAQUS. The analysis results obtained from the two codes showed good agreement in time-dependent response of the element under the loading area as well as the surface deformation of asphalt layer, and horizontal and vertical stresses along the axisymmetric axis. CONCLUSIONS: The validity of KICTPAVE was confirmed by showing the agreement of the analysis results from the two codes.

• Structural Engineering and Mechanics
• /
• v.15 no.6
• /
• pp.735-752
• /
• 2003

The quasi-static and dynamic responses of a linear viscoelastic circular beam on Winkler foundation are studied numerically by using the mixed finite element method in transformed Laplace-Carson space. This element VCR12 has 12 independent variables. The solution is obtained in transformed space and Schapery, Dubner, Durbin and Maximum Degree of Precision (MDOP) transform techniques are employed for numerical inversion. The performance of the method is presented by several quasi-static and dynamic example problems.

• Journal of the Society of Naval Architects of Korea
• /
• v.46 no.2
• /
• pp.155-164
• /
• 2009

In this paper, a spectral finite element method for a rectangular sandwich plate with viscoelastic core having the Levy-type boundary conditions has been plated. The sandwich plate consists of two isotropic and elastic face plates with a surfaced-bonded viscoelastic core. For the analysis, the in-plane and transverse energy in the face plates and only shear energy in the core are considered, respectively. To account for the frequency dependent complex shear modulus of the viscoelastic core, the Golla-Hughes-McTavish model is adopted. To evaluate the validity and accuracy of the proposed method, the frequency response function and dynamic responses of the sandwich plate with all edges simply supported subject to an impact load are calculated and compared with those calculated by a finite element method. Though these calculations, it is confirmed that the proposed method is very reliable and efficient one for vibration analysis of a rectangular sandwich plate with viscoelastic core having the Levy-type boundary conditions.

• Computational Structural Engineering
• /
• v.9 no.1
• /
• pp.101-113
• /
• 1996

In this paper, the finite element formulation for the thermal analysis of quasi-static, uncoupled, homogeneous, isotropic and linear viscoelastic problems is presented based on the principle of virtual work. The viscoelastic material is assumed to be thermorheologically simple, which is well known material property in a large class of high polymers. The variational formulation and the finite element equation in matrix from are derived. Effective generation and storage of the hereditary stiffness matrices are given in detail especially for the case of the steady state temperature distribution T=T(x). Some numerical examples are given and compared with published results to show the versatility of the derived finite element formulations.

• Computational Structural Engineering
• /
• v.6 no.2
• /
• pp.103-114
• /
• 1993

An advanced time-domain finite element formulation is presented for the displacement and stress analysis of isotropic, linear viscoelastic problems of a hereditary-type constitutive law. The semidiscrete finite element method with linear time-stepping scheme and an elastic-viscoelastic correspondence principle are used in the theoretical development. An efficient treatment of hereditary integral is introduced to improve the numerical accuracy, to reduce the computation time, and to avoid the use of large memory storage. Two-dimensional numerical examples of plane strain and plane stress are solved under the assumption on the material property of being elastic in dilatation and like three-element Voigt model in distorsion, and compared with the analytical solutions and the past numerical results to show the versatility and efficiency of the proposed method.

• Journal of Theoretical and Applied Mechanics
• /
• v.1 no.1
• /
• pp.97-109
• /
• 1995

In this paper, a finite viscoelastic continuum model for rubber and its finite element analysis are presented. This finite viscoelatic model based on continuum mechanics is an extended model of Johnson and Wuigley's 1-D model. In this extended model, continuum based kinematic measures are rigorously defied and by using this kinematic measures, elastic stage law and flow rule are introduced. In kinematics, three configuration are introduced. In kinematics, three configuration are introduced. They are reference, current and virtual visco configurations. In elastic state law, it is assumed that at a certain time, there exists an elastic potential which describes the recoverable elastic energy. From this elastic potential, elastic state law is derived. The proposed flow rule is based on phenomenological observation. The flow rule gives precise relaxation response. In finite element approximation, mixed Lagrangian description is used, where total and similar method of updated Lagrangian descriptions are used together. This approach reduces the numerical job and gives simple nonlinear syatem of equations. To satisfy the incompressible condition, penalty-type modified Mooney-Rivlin energy function is adopted. By this method nearly incompressible condition is obtain the virtual visco configuration. For verification, uniaxial stretch tests are simulated for various stretch rates. The simulated results show good agreement with experiments. As a practical experiments. As a preactical example, pressurized rubber plate is simulated. The result shows finite viscoelastic effects clearly.

• Geomechanics and Engineering
• /
• v.24 no.1
• /
• pp.91-103
• /
• 2021

This paper concerns with forced dynamic response of thick functionally graded (FG) beam resting on viscoelastic foundation including porosity impacts. The dynamic point load is proposed to be triangle point loads in time domain. In current analysis the beam is assumed to be thick, therefore, the two-dimensional plane stress constitutive equation is proposed to govern the stress-strain relationship through the thickness. The porosity and void included in constituent is described by three different distribution models through the beam thickness. The governing equations are obtained by using Lagrange's equations and solved by finite element method. In frame of finite element analysis, twelve-node 2D plane element is exploited to discretize the space domain of beam. In the solution of the dynamic problem, Newmark average acceleration method is used. In the numerical results, effects of porosity coefficient, porosity distribution and foundation parameters on the dynamic responses of functionally graded viscoelastic beam are presented and discussed. The current model is efficient in many applications used porous FGM, such as aerospace, nuclear, power plane sheller, and marine structures.

• Journal of Mechanical Science and Technology
• /
• v.17 no.2
• /
• pp.171-180
• /
• 2003

One of the most significant problems in the processing of composite materials is residual stress. The high residual stress may cause cracking in the matrix without external loads and degrade the integrity of composite structures. In this study, thermo-viscoelastic residual stresses occurred in an aluminum liner-inserted polymer composite cylinder are investigated. This type of the structure is used for rocket fuselage due to the convenience to attach payloads and equipment to the metal liner by machining. The time and degree of cure dependent thermo-viscoelastic constitutive equations are developed and coupled with a thermo-chemical process model. These equations are solved with the finite element method to predict the residual stresses in the composite cylinder and also in the interface between the liner and the composite during cure.

• Structural Engineering and Mechanics
• /
• v.47 no.1
• /
• pp.27-44
• /
• 2013

A finite element computational procedure for the accurate analysis of quasistatic thermorheological complex structures response is developed. The geometrical nonlinearity, arising from large displacements and rotations (but small strains), is accounted for by the total Lagrangian description of motion. The Schapery's nonlinear single-integral viscoelastic constitutive model is modified for a time-stress-temperature-dependent behavior. The nonlinear thermo-viscoelastic constitutive equations are incrementalized leading to a recursive relationship and thereby the resulting finite element equations necessitate data storage from the previous time step only, and not the entire deformation history. The Newton-Raphson iterative scheme is employed to obtain a converged solution for the non-linear finite element equations. The developed numerical model is verified with the previously published works and a good agreement with them is found. The applicability of the developed model is demonstrated by analyzing two examples with different thermal/mechanical loading histories.