• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.29 no.6 s.237
• /
• pp.913-920
• /
• 2005

Compared to previous in-vitro test, FE model showed reliable motion patterns. A finite element model of a 50th percentile male neck was developed to study the mechanics of whiplash injury while the rear impacts. The model was consisted of the whole cervical vertebrae including part of occipital, intervertebral discs. which were modeled using linear viscoelastic materials and posterior elements. The sliding interfaces were defined to simulate contact phenomena in facet joints and in odontoid process. All ligaments and atlanto-occipital membrane were modeled as nonlinear bar elements. Only muscle elements were not considered. Motion of each cervical vertebra was obtained from the dynamic simulation with a MADYMO model for 15 km/h $40\%$ rear end offset impacts. Soft tissue neck injury(STNI) was investigated with a developed FE model. In FE model analysis, the high stress was appeared at C3/C4 disc in offset impact. Further research is still needed in order to improve the developed neck FE model for many different crash patterns.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2005.10a
• /
• pp.63-66
• /
• 2005

Nanoimprint lithography (NIL) process at room temperature has been newly proposed in recent years to overcome the shape accuracy and sticking problem induced in a conventional NIL process. Success of the room temperature NIL relies on the accurate understand of the mechanical behavior of the polymer. Since a conventional NIL process has to heat a polymer above the glass transition temperature to deform the physical shape of the polymer with a mold pattern, viscoelastic property of polymer have major effect on the NIL process. However, rate dependent behavior of polymer is important in the room temperature NIL process because a mold with engraved patterns is rapidly pressed onto a substrate coated with the polymer by the hydraulic equipment. In this paper, finite element analysis of the room temperature NIL process is performed with considering the strain rate dependent behavior of the polymer. The analyses with the variation of imprinting speed and imprinting pattern are carried out in order to investigate the effect of such process parameters on the room temperature NIL process. The analyses results show that the deformed shape and imprint force is quite different with the variation of punch speed because the dynamic behavior of the polymer is considered with the rate dependent plasticity model. The results provide a guideline for the determination of process conditions in the room temperature NIL process.

• Structural Engineering and Mechanics
• /
• v.85 no.2
• /
• pp.163-177
• /
• 2023

The design of segmental bridges, a structure that typically employs precast prestressed concrete elements and the balanced cantilever construction method for the deck, may demand a highly complex structural analysis for increased precision of the results. This work presents a comprehensive numerical analysis of a 3D finite element model using the software ANSYS, version 21.2, to simulate the constructive deck stages of the New Guaiba Bridge, a structure located in Porto Alegre city, southern Brazil. The materials concrete and steel were considered viscoelastic. The concrete used a Generalized Kelvin model, with subroutines written in FORTRAN and added to the main model through the customization tool UPF (User Programmable Features). The steel prestressing tendons used a Generalized Maxwell model available in ANSYS. The balanced cantilever constructive steps of a span of the New Guaiba Bridge were then numerically simulated to follow the actual constructive sequence of the bridge. A comparison between the results obtained with the numerical model and the actual vertical displacement data monitored during the bridge's construction was carried out, showing a good correlation.

• Journal of KSNVE
• /
• v.9 no.1
• /
• pp.174-183
• /
• 1999

An analytical model for the lateral vibration of beams with a bonded lap joint and partial layered dampers has been proposed in this paper. Both shear and normal forces acting along the interface between the elastic and viscoelastic layers were considered in the vibration analysis. Analytical results were compared with those obtained by a finite element method. Effects of the size and location of layers in partial dampers on system loss factor($\eta_s$) and resonant frequency($\omega_r$) were studied. which showed that partial dampers adhered to the site exhibiting the maximum amplitude of vibration were most influential in the increase of $\eta_s$ and the decrease of $\omega_r$. Specific system loss factor( $\eta_s$ divided by total mass of system) was also evaluated in the analysis.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.05a
• /
• pp.311-315
• /
• 2002

The coupled finite/boundary element method is used in numerical analysis for acoustic radiation from the vibration of rectangular composite plate which is simply supported. This analysis is validated using the Wallace equation for an isotropic plate. Active control of sound fields has been carried out using 3 pairs of piezoelectric sensor/actuator and a pair of viscoelastic material by passive constrained layer damping treatment. The results show that the optimal placement of piezoelectric sensor/actuator and VE patch is required to control the sound fields from a vibrating composite plate.

• Journal of Mechanical Science and Technology
• /
• v.15 no.2
• /
• pp.199-209
• /
• 2001

Modal analysis method (MAM) is introduced for the fully coupled structural dynamic problems. In this paper, the beam with active constrained layered damping (ACLD) treatment is considered as a representative problem. The ACLD beam consists of a viscoelastic layer that is sandwiched between the base beam structure and an active piezoelectric layer. The exact damped natural modes are spectrally formulated from a set of fully coupled dynamic equations of motion. The orthogonality property of the exact damped natural modes is then derived in a closed form to complete the modal analysis method. The accuracy of the present MAM is evaluated through some illustrative examples: the dynamic characteristics obtained by the present MAM are compared with the results by spectral element method (SEM) and finite element method (FEM). It is numerically proved that MAM solutions become identical to the accurate SEM solutions as the number of exact natural used in MAM is increased.

• Korea-Australia Rheology Journal
• /
• v.16 no.4
• /
• pp.183-191
• /
• 2004

High Deborah or Weissenberg number problems in viscoelastic flow modeling have been known formidably difficult even in the inertialess limit. There exists almost no result that shows satisfactory accuracy and proper mesh convergence at the same time. However recently, quite a breakthrough seems to have been made in this field of computational rheology. So called matrix-logarithm (here we name it tensor-logarithm) formulation of the viscoelastic constitutive equations originally written in terms of the conformation tensor has been suggested by Fattal and Kupferman (2004) and its finite element implementation has been first presented by Hulsen (2004). Both the works have reported almost unbounded convergence limit in solving two benchmark problems. This new formulation incorporates proper polynomial interpolations of the log­arithm for the variables that exhibit steep exponential dependence near stagnation points, and it also strictly preserves the positive definiteness of the conformation tensor. In this study, we present an alternative pro­cedure for deriving the tensor-logarithmic representation of the differential constitutive equations and pro­vide a numerical example with the Leonov model in 4:1 planar contraction flows. Dramatic improvement of the computational algorithm with stable convergence has been demonstrated and it seems that there exists appropriate mesh convergence even though this conclusion requires further study. It is thought that this new formalism will work only for a few differential constitutive equations proven globally stable. Thus the math­ematical stability criteria perhaps play an important role on the choice and development of the suitable con­stitutive equations. In this respect, the Leonov viscoelastic model is quite feasible and becomes more essential since it has been proven globally stable and it offers the simplest form in the tensor-logarithmic formulation.

• Korea-Australia Rheology Journal
• /
• v.21 no.1
• /
• pp.47-58
• /
• 2009

We have investigated the transient behavior of 1D fully developed Poiseuille viscoelastic flow under finite pressure gradient described by the Oldroyd-B and Leonov constitutive equations. For analysis we employ a simple $2^{nd}$ order discretization scheme such as central difference for space and the Crank-Nicolson for time approximation. For the analysis of the Oldroyd-B model, we also apply the analytical solution, which is obtained again in this work in terms of elementary solution procedure simpler than the previous one (Waters and King, 1970). Both models demonstrate qualitatively similar solutions, but their eventual steady flowrate exhibits noticeable difference due to the absence or presence of shear thinning behavior. In the inertialess flow, the flowrate instantaneously attains a large value corresponding to the Newtonian creeping flow and then decreases to its steady value when the applied pressure gradient is low. However with finite liquid density the flow field shows severe fluctuation even accompanying reversals of flow directions. As the assigned pressure gradient increases, the flowrate achieves its steady value significantly higher than its value during oscillations after quite long period of time. We have also illustrated comparison between 1D and 2D results and possible mechanism of complex 2D flow rearrangement employing a previous solution of [mite element computation. In addition, we discuss some mathematical points regarding missing boundary conditions in 2D modeling due to the change of the type of differential equations when varying from inertialess to inertial flow.

• Smart Structures and Systems
• /
• v.18 no.6
• /
• pp.1125-1143
• /
• 2016

Forced vibration analysis of a simple supported viscoelastic nanobeam is studied based on modified couple stress theory (MCST). The nanobeam is excited by a transverse triangular force impulse modulated by a harmonic motion. The elastic medium is considered as Winkler-Pasternak elastic foundation.The damping effect is considered by using the Kelvin-Voigt viscoelastic model. The inclusion of an additional material parameter enables the new beam model to capture the size effect. The new non-classical beam model reduces to the classical beam model when the length scale parameter is set to zero. The considered problem is investigated within the Timoshenko beam theory by using finite element method. The effects of the transverse shear deformation and rotary inertia are included according to the Timoshenko beam theory. The obtained system of differential equations is reduced to a linear algebraic equation system and solved in the time domain by using Newmark average acceleration method. Numerical results are presented to investigate the influences the material length scale parameter, the parameter of the elastic medium and aspect ratio on the dynamic response of the nanobeam. Also, the difference between the classical beam theory (CBT) and modified couple stress theory is investigated for forced vibration responses of nanobeams.

• Structural Engineering and Mechanics
• /
• v.41 no.1
• /
• pp.113-137
• /
• 2012

Frame structures with viscoelastic (VE) dampers mounted on them are considered in this paper. It is the aim of this paper to compare the dynamic characteristics of frame structures with VE dampers when the dampers are modelled by means of different models. The classical rheological models, the model with the fractional order derivative, and the complex modulus model are used. A relatively large structure with VE dampers is considered in order to make the results of comparison more representative. The formulae for dissipation energy are derived. The finite element method is used to derive the equations of motion of the structure with dampers and such equations are written in terms of both physical and state-space variables. The solution to motion equations in the frequency domain is given and the dynamic properties of the structure with VE dampers are determined as a solution to the appropriately defined eigenvalue problem. Several conclusions concerning the applicability of a family of models of VE dampers are formulated on the basis of results of an extensive numerical analysis.