• Structural Engineering and Mechanics
• /
• v.42 no.4
• /
• pp.471-488
• /
• 2012

An efficient one dimensional finite element model has been presented for the dynamic analysis of composite laminated beams, using the efficient layerwise zigzag theory. To meet the convergence requirements for the weak integral formulation, cubic Hermite interpolation is used for the transverse displacement ($w_0$), and linear interpolation is used for the axial displacement ($u_0$) and shear rotation (${\psi}_0$). Each node of an element has four degrees of freedom. The expressions of variationally consistent inertia, stiffness matrices and the load vector are derived in closed form using exact integration. The formulation is validated by comparing the results with the 2D-FE results for composite symmetric and sandwich beams with various end conditions. The employed finite element model is free of shear locking. The present zigzag finite element results for natural frequencies, mode shapes of cantilever and clamped-clamped beams are obtained with a one-dimensional finite element codes developed in MATLAB. These 1D-FE results for cantilever and clamped beams are compared with the 2D-FE results obtained using ABAQUS to show the accuracy of the developed MATLAB code, for zigzag theory for these boundary conditions. This comparison establishes the accuracy of zigzag finite element analysis for dynamic response under given boundary conditions.

• Proceedings of the Korean Society for Technology of Plasticity Conference
• /
• 2004.08a
• /
• pp.72-82
• /
• 2004

The deformation mechanism of hexagonal close-packed materials is quite complicate including slips and twins. A deformation mechanism, which accounts for both slip and twinning, was investigated for polycrystalline hop materials. The model was developed in a finite element polycrystal model formulated with initial strain method where the stiffness matrix in FEM is based on the elastic modulus. We predicted numerically the texture of Mg alloy(AZ31B) sheet by using FEM based on crystal plasticity theory. Also, we introduced the recrystallized texture employed the maximum energy release theory after rolling. From the numerical study, it was clarified that the shrink twin could not be the main mechanism for shortening of c-axis, because the lattice rotation due to twin rejects fur c-axis to become parallel to ND(normal direction of plate). It was showed that the deformation texture with the pyramidal slip gives the ring type pole figure having hole in the center.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.25 no.11
• /
• pp.1697-1704
• /
• 2001

A finite element model for polarization switching in electro-mechanically coupled materials is proposed and applied to predict the switching behavior of a two-dimensional ferroelectric ceramic. A complicated micro-structure existing in the material is modeled as il continuum body and a simple 3 node triangle finite element with nodal displacement and voltage degrees of freedom is used for a finite element analysis. The elements use nonlinear constitutive equations, switching criterion and kinetic relation, fur representation of material response at strong electric and stress fields. The polarization state of the material is represented by internal variables in each element, which are updated at each simulation step based on the proposed constitutive equations. The model reproduces strain and electric displacement hysteresis loops observed in the material.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.12 no.2
• /
• pp.141-149
• /
• 2002

Although the finite element method has become an indispensible tool for the dynamic analysis of structures, difficulty remains to quantify the errors associated with discretization. To improve the modeling accuracy, this paper proposes a method to make a combined use of finite elements and exact dynamic elements. Exact interpolation functions for the Timoshenko beam element are derived using the exact dynamic element modeling (EDEM) and compared with interpolation functions of the finite element method (FEM). The exact interpolation functions are tested with the Laplace variable varied. A combined use of finite element method and exact interpolation functions is presented to gain more accurate mode shape functions. This paper also presents a combined use of finite elements and exact dynamic elements in design/reanalysis problems. Timoshenko flames with tapered sections are tested to demonstrate the design procedure with the proposed method. The numerical study shows that the combined use of finite element model and exact dynamic element model is very useful.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2004.05a
• /
• pp.645-650
• /
• 2004

There have been lots of efforts to analyze the dynamic characteristics of mechanical systems. However, it is difficult to know the dynamic characteristics of mechanical systems composed of many parts with joints. Specially, in case of a bolted joint structure, no effective modeling method has been defined to acquire dynamic characteristics of the structure, using the finite element (FE) analysis. In this research, a linear dynamic model is developed for bolted joints and large interfaces using con frusta method and linear spring elements, respectively. The developed modeling method for bolted joints is verified based on the experimental result.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2002.11b
• /
• pp.1049-1053
• /
• 2002

An unconstrained tuning fork with a 3-D model has been numerically analyzed by Finite Element Method (FEM) and Boundary Element Method (BEM). The first three natural frequencies were calculated by the FEM modal analysis. Then the change of the modal frequencies was examined with the variation of the tuning fork length and width. Analytical model equations were derived from the numerically relating results of the modal frequency-tuning fork length by approximating minimization. Finally the BEM was used for the sound pressure field calculation from the structural displacement data.

• Proceedings of the Acoustical Society of Korea Conference
• /
• spring
• /
• pp.465-468
• /
• 2002

An unconstrained tuning fork with a 3-D model has been numerically analyzed by Finite Element Method (FEM) and Boundary Element Method (BEM). The first three natural frequencies were calculated by the FEM modal analysis. Then the change of the modal frequencies was examined with the variation of the tuning fork length and width. Analytical model equations were derived from the numerically relating results of the modal frequency-tuning fork length by approximating minimization. Finally the BEM was used for the sound pressure field calculation from the structural displacement data.

• Coupled systems mechanics
• /
• v.4 no.3
• /
• pp.263-277
• /
• 2015

This paper presents a coupled FEM-BEM strategy for the numerical analysis of elastodynamic problems where infinite-domain models and complex heterogeneous media are involved, rendering a configuration in which neither the Finite Element Method (FEM) nor the Boundary Element Method (BEM) is most appropriate for the numerical analysis. In this case, the coupling of these methodologies is recommended, allowing exploring their respective advantages. Here, frequency domain analyses are focused and an iterative FEM-BEM coupling technique is considered. In this iterative coupling, each sub-domain of the model is solved separately, and the variables at the common interfaces are iteratively updated, until convergence is achieved. A relaxation parameter is introduced into the coupling algorithm and an expression for its optimal value is deduced. The iterative FEM-BEM coupling technique allows independent discretizations to be efficiently employed for both finite and boundary element methods, without any requirement of matching nodes at the common interfaces. In addition, it leads to smaller and better-conditioned systems of equations (different solvers, suitable for each sub-domain, may be employed), which do not need to be treated (inverted, triangularized etc.) at each iterative step, providing an accurate and efficient methodology.

• Journal of The Korean Society of Agricultural Engineers
• /
• v.57 no.6
• /
• pp.1-7
• /
• 2015

Most of chloride diffusion models based on finite difference method (FDM) could not express the diffusion in horizontal direction at each elevation. To overcome these weakness, two dimensional chloride ion penetration model based on finite element method (FEM) to be able to combine various multi-physics simultaneously was suggested by introducing moving mesh technique. To avoid the generation of mesh being able to be distorted depending on the relative movement of water level to static concrete, a rectangular type of mesh was intentionally adopted and the total number of meshes was empirically selected. The simulated results showed that the contents of surface chloride decreased following to the increase of elevation in the top part of low sea level, whereas there were no changes in the bottom part of low level. In the DuraCrete model, the diffusion coefficient of splashed zone is generally smaller than submerged zone, whereas the trend of Life365 model is reverse. Therefore, it could be understood that the developed model using moving mesh technique effectively reflects $DuraCrete^{TM}$ model rather than $Life365^{TM}$ model. In the future, the model will be easily expanded to be combined with various multi-physics models considering water evaporation, heat of hydration, irradiation effect of sun and so on because it is based on FEM.

• Smart Structures and Systems
• /
• v.13 no.4
• /
• pp.587-614
• /
• 2014

In recent years the interest in online monitoring of lightweight structures with ultrasonic guided waves is steadily growing. Especially the aircraft industry is a driving force in the development of structural health monitoring (SHM) systems. In order to optimally design SHM systems powerful and efficient numerical simulation tools to predict the behaviour of ultrasonic elastic waves in thin-walled structures are required. It has been shown that in real industrial applications, such as airplane wings or fuselages, conventional linear and quadratic pure displacement finite elements commonly used to model ultrasonic elastic waves quickly reach their limits. The required mesh density, to obtain good quality solutions, results in enormous computational costs when solving the wave propagation problem in the time domain. To resolve this problem different possibilities are available. Analytical methods and higher order finite element method approaches (HO-FEM), like p-FEM, spectral elements, spectral analysis and isogeometric analysis, are among them. Although analytical approaches offer fast and accurate results, they are limited to rather simple geometries. On the other hand, the application of higher order finite element schemes is a computationally demanding task. The drawbacks of both methods can be circumvented if regions of complex geometry are modelled using a HO-FEM approach while the response of the remaining structure is computed utilizing an analytical approach. The objective of the paper is to present an efficient method to couple different HO-FEM schemes with an analytical description of an undisturbed region. Using this hybrid formulation the numerical effort can be drastically reduced. The functionality of the proposed scheme is demonstrated by studying the propagation of ultrasonic guided waves in plates, excited by a piezoelectric patch actuator. The actuator is modelled utilizing higher order coupled field finite elements, whereas the homogenous, isotropic plate is described analytically. The results of this "semi-analytical" approach highlight the opportunities to reduce the numerical effort if closed-form solutions are partially available.