• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.13 no.3
• /
• pp.225-231
• /
• 2003

A new three-dimensional thin shell element for a structure containing an integrated distributed piezoelectric sensor and actuator is Proposed. The assumed strain formulation and the bubble function are introduced to improve the performance of the shell element. A finite element formulation gives a general tool that can predict the static or dynamic responses of the shell with piezoelectric sensor/actuator. The verification through the calculation of the static response for the piezoelectric bimorph beam shows that the results agree with those from the theoretical analysis very well. Dynamic response of a shell shows that the reduction of vibration is possible with the introduction of the piezoelectric shell sensor and actuator. However. the curvature of sensor/actuator is an obstacle for application, since the flexible PVDF is not strong enough and the PZT with curvature should be made specially.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2001.11a
• /
• pp.225-228
• /
• 2001

A new three-dimensional thin shell element for the structure containing an integrated distributed piezoelectric sensor and actuator is proposed. A finite element formulation for the static response of the shell with piezoelectric sensor／actuator is derived. The assumed strain formulation and the bubble function improves the performance of the shell element. The verification through the calculation of the static response for the piezoelectic bimorph beam shows that the results agree with those from the theoretical analysis very well.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.9
• /
• pp.1579-1588
• /
• 2003

In this paper, the finite element equations for the transient linear viscoelastic stress analysis are presented in time domain, whose variational formulation is derived by using the Galerkin's method based on the equations of motion in time integral. Since the inertia terms are not included in the variational formulation, the time integration schemes such as the Newmark's method widely used in the classical dynamic analysis based on the equations of motion in time differential are not required in the development of that formulation, resulting in a computationally simple and stable numerical algorithm. The viscoelastic material is assumed to behave as a standard linear solid in shear and an elastic solid in dilatation. To show the validity of the presented method, two numerical examples are solved nuder plane strain and plane stress conditions and good results are obtained.

• Smart Structures and Systems
• /
• v.14 no.4
• /
• pp.699-718
• /
• 2014

The indirect displacement estimation using acceleration and strain (IDEAS) method is extended to various types of beam structures beyond the previous validation on the prismatic or near-prismatic beams. By fusing different types of responses, the IDEAS method is able to estimate displacements containing pseudo-static components with high frequency noise to be significantly reduced. However, the concerns to the IDEAS method come from possible disagreement of the assumed sinusoidal mode shapes to the actual mode shapes, which allows the IDEAS method to be valid only for simply-supported prismatic beams and limits its applicability to real world problems. In this paper, the extension of the IDEAS method to the general types of beams is investigated by the mathematical formulation of the modal mapping matrix only for the monitored substructure, so-called monitoring span. The formulation particularly considers continuous and wide beams to extend the IDEAS method to general beam structures that reflect many real bridges. Numerical simulations using four types of beams with various irregularities are presented to show the effectiveness and accuracy of the IDEAS method in estimating displacements.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 1996.10a
• /
• pp.38-47
• /
• 1996

Sectional analysis program for plane strain or axisymmetric geometry of aluminum alloy sheet metals was developed. For modeling the anomalous behavior of aluminum alloy, Barlat's strain rate potential and Hill's 1990 non-quadratic yield theory arranged under the plane stress assumption were employed. 2-D rigid-viscoplastic FEM formulation based on the bending-augmented membrane theory was derived, solving simultaneously force equilibrium as well as non-penetration condition. Isotropic hardening law was also assumed for yielding behavior. To verify the validity and availability of the developed program, 2-D stretch/draw forming process for plane strain geometry and cylindrical cup deep drawing process for axisymmetric geometry were simulated.

• Structural Engineering and Mechanics
• /
• v.49 no.6
• /
• pp.793-810
• /
• 2014

This research work deals with the development of a new Triangular finite element for the linear analysis of plate bending with transverse shear effect. It is developed in perspective to building shell elements. The displacements field of the element has been developed by the use of the strain-based approach and it is based on the assumed independent functions for the various components of strain insofar as it is allowed by the compatibility equations. Its formulation uses also concepts related to the fourth fictitious node, the static condensation and analytic integration. It is based on the assumptions of tick plate.s theory (Reissner-Mindlin theory). The element possesses three essential external degrees of freedom at each of the four nodes and satisfies the exact representation of the rigid body modes of displacements. As a result of this approach, a new bending plate finite element (Pep43) which is competitive, robust and efficient.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.19 no.2
• /
• pp.486-494
• /
• 1995

This paper introduces two kinds of new four-node quadrilateral axisymmetric elements with independently-assumed strains. They are formulated by the modified Hellinger-Reissner principle so as to employ incompatible displacements and assumed strains. In one of the present elements, the strains from incompatible displacements are corrected to pass the constant strain patch test. The other contains incompatible functions that are obtained from the element boundary condition. The elements are evaluated on the several problems of bending and material incompressibility with regular and distorted meshes. The results show that the new element performs excellently in deformation and stress calculation.

• Smart Structures and Systems
• /
• v.30 no.3
• /
• pp.221-237
• /
• 2022

This paper is concerned with nonlinear modeling and efficient numerical simulation of electrostrictive materials and structures. Two types of such materials are considered: relaxor ferroelectric ceramics and electrostrictive polymers. For ceramics, a geometrically linear formulation is developed, whereas polymers are studied in a geometrically nonlinear regime. In the paper, we focus on constitutive modeling first. For the reversible constitutive response under consideration, we introduce the augmented Helmholtz free energy, which is composed of a purely elastic part, a dielectric part and an augmentation term. For the elastic part, we involve an additive decomposition of the strain tensor into an elastic strain and an electrostrictive eigenstrain, which depends on the polarization of the material. In the geometrically nonlinear case, a corresponding multiplicative decomposition of the deformation gradient tensor replaces the additive strain decomposition used in the geometrically linear formulation. For the dielectric part, we first introduce the internal energy, to which a Legendre transformation is applied to compute the free energy. The augmentation term accounts for the contribution from vacuum to the energy. In our formulation, the augmented free energy depends not only on the strain and the electric field, but also on the polarization and an internal polarization; the latter two are internal variables. With the constitutive framework established, a Finite Element implementation is briefly discussed. We use high-order elements for the discretization of the independent variables, which include also the internal variables and, in case the material is assumed incompressible, the hydrostatic pressure, which is introduced as a Lagrange multiplier. The elements are implemented in the open source code Netgen/NGSolve. Finally, example problems are solved for both, relaxor ferroelectric ceramics and electrostrictive polymers. We focus on thin plate-type structures to show the efficiency of the numerical scheme and its applicability to thin electrostrictive structures.

• Journal of Ocean Engineering and Technology
• /
• v.16 no.6
• /
• pp.55-59
• /
• 2002

Delamination means that cracking occurs on the interface layer between composite laminates. In this paper, to predict the delamination growth in composite laminates subjected to low-velocity impact, the unit load method was introduced, and an eighteen-node 3-D finite element analysis, based on assumed strain mixed formulation, was conducted. Strain energy release rate, necessary to determine the delamination growth, was calculated by using the virtual crack closure technique. The unit load method saves the computation time more than the re-meshing method. The virtual crack closure technique enables the strain energy release rate to be easily calculated, because information of the singular stress field near the crack tip is not required. Hence, the delamination growth in composite laminates that are subjected to low-velocity impact can be efficiently predicted using the above-mentioned methods.

• Land and Housing Review
• /
• v.1 no.1
• /
• pp.59-65
• /
• 2010

The paper presents a comparison between a semi-implicit time integration linear finite element implementation and fully-implicit nonlinear Newton-Raphson finite element implementation of a triphasic small strain mixture formulation of an elastic partially saturated porous medium. The pore air phase pressure pa is assumed atmospheric, i.e., $p_a$ = 0, although the formulation and implementation are general to handle increase in pore air pressure as a result of loading, if needed. The solid skeleton phase is assumed linear isotropic elastic and partially saturated 'consolidation' in the presence of surface infiltration and traction is simulated. The verification of the implementation against an analytical solution for partially saturated pore water flow (no deformation) and comparison between the two implementations is presented and the important of the porosity-dependent nature of the partially saturated permeability is assessed on comparison with a commercial code for the partially saturated flow with deformation. As a result, the response of partially saturated permeability subjected to the porosity influences on the saturation of a soil, and the different behaviors of the partially saturated soil between staggered and monolithic coupled programs is worth of attention because the negative pore water pressure in the partially saturated soil depends on the difference.