• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.15 no.7 s.100
• /
• pp.836-842
• /
• 2005

Nonlinear dynamic characteristics of active piezolaminated plates are investigated with respect to the thermopiezoelastic behaviors. For largely deformed structures with small strain, the incremental total Lagrangian formulation is presented based on the virtual work principles. A multi-field layer-wise finite shell element is proposed for assuring high accuracy and non-linearity of displacement, electric and thermal fields. For dynamic consideration of thermopiezoelastic snap-through phenomena, the implicit Newmark's scheme with the Newton-Raphson iteration is implemented for the transient response of various piezolaminated models with symmetric or eccentric active layers. The bifurcate thermal buckling of symmetric structural models is first investigated and the characteristics of piezoelectric active responses are studied for finding snap-through piezoelectric potentials and the load-path tracking map. The thermoelastic stable and unstable postbuckling, thermopiezoelastic snap-through phenomena with several attractors are proved using the nonlinear time responses for various initial conditions and damping loss factors. Present results show that thermopiezoelastic snap-through phenomena can result in the difficulty of buckling and postbuckling control of intelligent structures.

• Journal of Theoretical and Applied Mechanics
• /
• v.2 no.1
• /
• pp.31-50
• /
• 1996

In this paper, a finite element analysis based on the local approach concept to fracture in the continuum damage mechanics is performed to analyze ductile fracture in two dimensional quasi-static state. First an isotropic damage model based on the generalized concept of effective stress is proposed for structural materials in the context of large deformation. In this model, the stiffness degradation is taken as a measure of damage and so, the fracture phenomenon can be explained as the critical deterioration of stiffness at a material point. The modified Riks' continuation technique is used to solve incremental iterative equations. Crack propagation is achieved by removing critically damaged elements. The mesh size sensitivity analysis and the simulation of the well known shearing mode failure in plane strain state are carried out to verify the present formulation. As numerical examples, an edge cracked plate and the specimen with a circular hole under plane stress are taken. Load-displacement curves and successively fractured shapes are shown. From the results, it can be concluded that the proposed model based on the local approach concept in the continuum damage mechanics may be stated as a reasonable tool to explain ductile fracture initiation and crack propagation.

• Structural Engineering and Mechanics
• /
• v.56 no.6
• /
• pp.917-938
• /
• 2015

The nonlinear numerical analysis of the impact response of reinforced concrete/mortar beam incorporated with the updated Lagrangian method, namely the Smoothed Particle Hydrodynamics (SPH) is carried out in this study. The analysis includes the simulation of the effects of high mass low velocity impact load falling on beam structures. Three material models to describe the localized failure of structural elements are: (1) linear pressure-sensitive yield criteria (Drucker-Prager type) in the pre-peak regime for the concrete/mortar meanwhile, the shear strain energy criterion (Von Mises) is applied for the steel reinforcement (2) nonlinear hardening law by means of modified linear Drucker-Prager envelope by employing the plane cap surface to simulate the irreversible plastic behavior of concrete/mortar (3) implementation of linear and nonlinear softening in tension and compression regions, respectively, to express the complex behavior of concrete material during short time loading condition. Validation upon existing experimental test results is conducted, from which the impact behavior of concrete beams are best described using the SPH model adopting an average velocity and erosion algorithm, where instability in terms of numerical fragmentation is reduced considerably.

• Journal of the Korean Society for Precision Engineering
• /
• v.15 no.9
• /
• pp.117-126
• /
• 1998

For investigating rubber pad sheet metal forming process, the rubber pad deformation characteristics as well as the contact problem of rubber pad-sheet metal has been analyzed. In this paper, the behavior of the rubber deformation is represented by hyper-elastic constitutive relations based on a generalized Mooney-Rivlin model. Finite element procedures for the two-dimensional responses, employing total Lagrangian formulations are implemented in an implicit form. The volumetric incompressibility condition of the rubber deformation is included in the formulation by using penalty method. The sheet metal is characterized by elasto-plastic material with strain hardening effect and analyzed by a commercial code. The contact procedure and interface program between rubber pad and sheet metal are implemented. Inflation experiment of circular rubber pad identifies the behaviour of the rubber pad deformation during the process. The various form dies and scaled down apparatus of the rubber-pad forming process are fabricated for simulating realistic forming process. The obtaining experimental data and FEM solutions were compared. The numerical solutions illustrate fair agreement with experimental results. The forming pressure distribution according to the dimensions of sheet metal and rubber pads, various rubber models and rubber material are also compared and discussed.

• Structural Engineering and Mechanics
• /
• v.17 no.6
• /
• pp.751-766
• /
• 2004

This paper extends the use of the hierarchic degenerated shell element to geometric non-linear analysis of composite laminated skew plates by the p-version of the finite element method. For the geometric non-linear analysis, the total Lagrangian formulation is adopted with moderately large displacement and small strain being accounted for in the sense of von Karman hypothesis. The present model is based on equivalent-single layer laminate theory with the first order shear deformation including a shear correction factor of 5/6. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. A wide variety of linear and non-linear results obtained by the p-version finite element model are presented for the laminated skew plates as well as laminated square plates. A numerical analysis is made to illustrate the influence of the geometric non-linear effect on the transverse deflections and the stresses with respect to width/depth ratio (a/h), skew angle (${\beta}$), and stacking sequence of layers. The present results are in good agreement with the results in literatures.

• Smart Structures and Systems
• /
• v.5 no.1
• /
• pp.23-42
• /
• 2009

In order to reduce vibration or to control shape of structures made of metal or composites, piezoelectric materials have been extensively used since their discovery in 1880's. A recent trend is also seen to apply piezoelectric materials to flexible structures made of rubber-like materials. In this paper a non-linear finite element model using updated Lagrangian (UL) approach has been developed for static analysis of rubber-elastic material with surface-bonded piezoelectric patches. A compressible stain energy function has been used for modeling the rubber as hyperelastic material. For formulation of the nonlinear finite element model a twenty-node brick element is used. Four degrees of freedom u, v and w and electrical potential ${\varphi}$ per node are considered as the field variables. PVDF (polyvinylidene fluoride) patches are applied as sensors/actuators or sensors and actuators. The present model has been applied to bimorph PVDF cantilever beam to validate the formulation. It is then applied to study the smart rubber components under different boundary and loading conditions. The results predicted by the present formulation are compared with the analytical solutions as well as the available published results. Some results are given as new ones as no published solutions available in the literatures to the best of the authors' knowledge.

• Earthquakes and Structures
• /
• v.26 no.1
• /
• pp.59-75
• /
• 2024

This paper presents a dam-reservoir interaction model that includes, water compressibility, sloshing of surface water, and radiation damping at the far-end reservoir, to investigate the influence of concrete deterioration on seismic behavior along with seismic performance of gravity dams. Investigations on seismic performance of the dam body have been conducted using the linear time-history responses obtained under six real and 0.3 g normalized earthquake records with time durations from 10 sec to 80 sec. The deterioration of concrete is assumed to develop due to mechanical and chemical actions over the dam lifespan. Several computer programs have been developed in FORTRAN 90 and MATLAB programming languages to analyze the coupled problem considering two-dimensional (2D) plane-strain condition. According to the results obtained from this study, the dam structure shows critical responses at the later ages (75 years) that could cause disastrous consequences; the critical effects of some earthquake loads such as Chi-Chi with 36.5% damage and Loma with 56.2% damage at the later ages of the selected dam body cannot be negligible; and therefore, the deterioration of concrete along with its effects on the dam response should be considered in analysis and design.

• Computers and Concrete
• /
• v.16 no.6
• /
• pp.933-961
• /
• 2015

Damage by high-speed impact fracture is a dominant mode of failure in several applications of concrete structures. Numerical modelling can play a crucial role in understanding and predicting complex fracture processes. The commonly used mesh-based Finite Element Method has difficulties in accurately modelling the high deformation and disintegration associated with fracture, as this often distorts the mesh. Even with careful re-meshing FEM often fails to handle extreme deformations and results in poor accuracy. Moreover, simulating the mechanism of fragmentation requires detachment of elements along their boundaries, and this needs a fine mesh to allow the natural propagation of damage/cracks. Smoothed Particle Hydrodynamics (SPH) is an alternative particle based (mesh-less) Lagrangian method that is particularly suitable for analysing fracture because of its capability to model large deformation and to track free surfaces generated due to fracturing. Here we demonstrate the capabilities of SPH for predicting brittle fracture by studying a slender concrete structure (column) under the impact of a high-speed projectile. To explore the effect of the projectile material behaviour on the fracture process, the projectile is assumed to be either perfectly-elastic or elastoplastic in two separate cases. The transient stress field and the resulting evolution of damage under impact are investigated. The nature of the collision and the constitutive behaviour are found to considerably affect the fracture process for the structure including the crack propagation rates, and the size and motion of the fragments. The progress of fracture is tracked by measuring the average damage level of the structure and the extent of energy dissipation, which depend strongly on the type of collision. The effect of fracture property (failure strain) of the concrete due to its various compositions is found to have a profound effect on the damage and fragmentation pattern of the structure.

• Composites Research
• /
• v.14 no.6
• /
• pp.16-23
• /
• 2001

Analysis of mechanical behavior for a pultruded-wound hollow rod is presented. For this purpose, the pultruded-wound hollow rod is manufactured by the new winder attached to the conventional pultrusion system. And the conventional pultrusion process is newly altered to manufacture pultruded-wound specimens. A computer program, POST II, is modified to perform this study, In the nonlinear finite element formulation, the updated Lagrangian description method based on the second Piolar-Kirchhoff stress tensor and the Green strain tensor are used. For the finite element modeling of the composite hollow rod, the eight-node degenerated shell element is utilized. In order to estimate the failure, the maximum stress criterion is adopted to the averaged stress in the each layer of the finite elements. As numerical examples, the behavior of glass/up composite hollow rod is investigated from the initial loading to the final collapse. Present finite element results considering stiffness degradation and stress unload due to failure shows excellent agreement with experiments in the ultimate load, failure and deformations.

• Journal of Korean Association for Spatial Structures
• /
• v.6 no.4 s.22
• /
• pp.53-61
• /
• 2006

In th is paper, to predicting the large deformation and cyclic plastic behavior of steel members under loading, 3-Dimensional elastic-plastic FE analysis method is developed by using finite deformation theory and proposed cyclic plasticity model. finite deformation theory, described the large deformation, is formulated by using Updated-lagrangian formulation and Green's strain tensor, Jaumann's derivative of Kirchoff stress. Also, cyclic plasticity model proposed by author is applied to developed analysis method. To verification of developed analysis method, analysis result of steel plate specimen compare to the analysis result using infinitesimal deformation theory and test result. Also, load-displacement and deflection shape, analysis result of pipe-section steel column, compare to test result. The good agreement between analysis result and experiment result shown that developed 3-dimensional finite element analysis can be predict the large deformation and cyclic plastic behavior of steel members.