• 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.

• Journal of the Korean Magnetics Society
• /
• v.7 no.2
• /
• pp.109-116
• /
• 1997

This paper describes an implementation of the rotational hysteresis analysis using the finite element method and a vector hysteresis model. The vector hysteresis model is extended from the magnetization-dependent Preisach model. The magnetization of each element is calculated for a given history of the rotating field. The lagging angle and the magnitude of the magnetization for an applied field are evaluated. The experimental results show that the magnetization characteristics of the magnetic materials under the rotating field can be analyzed accurately by the proposed method.

• Journal of Power Electronics
• /
• v.11 no.4
• /
• pp.401-407
• /
• 2011

In order to enhance torque density of five-phase permanent magnetic synchronous motor with third harmonic injection for electric vehicles (EVs) applications, optimum seeking method for injection ratio of third harmonic was proposed adopting theoretical derivation and finite element analysis method, under the constraint of same amplitude for current and air-gap flux. By five-dimension space vector decomposition, the mathematic model in two orthogonal space plane, $d_1-q_1$ and $d_3-q_3$, was deduced. And the corresponding dual-plane vector control method was accomplished to independently control fundamental and third harmonic currents in each vector plane. A five-phase PMSM prototype with quasi-trapezoidal flux pattern and its fivephase voltage source inverter were designed. Also, the dual-plane vector control was digitized in a single XC3S1200E FPGA. Simulation and experimental results prove that using the proposed optimum seeking method, the torque density of five-phase PMSM is enhanced by 20%, without any increase of power converter capacity, machine size and iron core saturation.

• Journal of applied mathematics & informatics
• /
• v.15 no.1_2
• /
• pp.29-51
• /
• 2004

In this paper, we propose a new mixed finite element method, called the characteristics-mixed method, for approximating the solution to Burgers' equation. This method is based upon a space-time variational form of Burgers' equation. The hyperbolic part of the equation is approximated along the characteristics in time and the diffusion part is approximated by a mixed finite element method of lowest order. The scheme is locally conservative since fluid is transported along the approximate characteristics on the discrete level and the test function can be piecewise constant. Our analysis show the new method approximate the scalar unknown and the vector flux optimally and simultaneously. We also show this scheme has much smaller time-truncation errors than those of standard methods. Numerical example is presented to show that the new scheme is easily implemented, shocks and boundary layers are handled with almost no oscillations. One of the contributions of the paper is to show how the optimal error estimates in $L^2(\Omega)$ are obtained which are much more difficult than in the standard finite element methods. These results seem to be new in the literature of finite element methods.

• Transactions of the Korean Society of Machine Tool Engineers
• /
• v.10 no.1
• /
• pp.16-22
• /
• 2001

In this study for reduction degree of freedom of dynamic model, a mixed method to combined finite element method and transfer matrix method is presented. This offers the advantages of an automatic reduction in the size of the eigenvalues problem and of a straightforward means of dynamic substructuring. The analytical procedure in this method for dynamic analysis of 3-dimensional cantilevered box beam are described. the result of numerical example is shown to demonstate the efficiency and accuracy of this method. The result form this example agree well those obtained by ANSYS, By using this technique, the number of nodes required in the regular finite element method is reduced and therefore a smaller com-puter can be used.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.10a
• /
• pp.69-76
• /
• 1998

Due to the drastic improvement of computer hardware and operating system, it is easy to break through the main defects of limited computer memory and processing time, etc. To keep up with this situation, this paper is focused on developing the preprocessor program with the input method based on vector graphic editor and the preprocessing technique including automatic node generation algorithm for the $\rho$-version finite element program. To develop this preprocessor program, the special data structure and the OOP(Object Oriented Programming) have been used by the Visual Basic 4.0. The Special data structure is proposed to describe the geometric data of node numberings and coordinates suitable for the $\rho$-version finite element program, which are quite different from the comvential h-version finite element program.

• Proceedings of the KIEE Conference
• /
• 1987.07a
• /
• pp.589-592
• /
• 1987

• Structural Engineering and Mechanics
• /
• v.40 no.2
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.17 no.9
• /
• pp.2242-2251
• /
• 1993

The mesh-based frictional contact model has been developed which does not rely on the spatial derivatives of the tool surface. Only points on the surface are evaluated from the description. which can then be simplified because of the relaxed demands placed on it. The surface tangents, normals, and corresponding derivatives at each finite-element node are evaluated directly from the finite-element mesh, in terms of the connecting nodal positions. The advantages accrue because there is no longer a need for a smooth tool surface to assure reasonable normals and derivatives. Furthermore, it can be shown that the equilibrium equations can only be properly written with a special normal derived from the mesh itself. The validity, accuracy, computation time, and stability of mesh-based contact model were discussed with the numerical examples of rounded flat-top and rough, flat-top rounded punch forming operations. Also, the forming process of a automobile inner panel section was simulated for testing the robustness of new contact model. In the discussion, the superiority of new model was examined, comparing with tool-based contact one.

• Structural Engineering and Mechanics
• /
• v.75 no.5
• /
• pp.585-594
• /
• 2020

Detecting the damage of indeterminate trusses is of major importance in the literature. This paper proposes a quick approach in this regard, utilizing a precise mathematical approach based on Finite Element Method. Different to a general two-step method defined in the literature essentially based on optimization approach, this method consists of three steps including Damage-Suspected Element Identification step, Imminent Damaged Element Identification step, and finally, Damage Severity Detection step and does not need any optimizing algorithm. The first step focuses on the identification of damage-suspected elements using an index based on modal residual force vector. In the second step, imminent damage elements are identified among the damage-suspected elements detected in the previous step using a specific technique. Ultimately, in the third step, a novel relation is derived to calculate the damage severity of each imminent damaged element. To show the efficiency and quick function of the proposed method, three examples including a 25-bar planar truss, a 31-bar planar truss, and a 52-bar space truss are studied; results of which indicate that the method is innovatively capable of suitably detecting, for indeterminate trusses, not only damaged elements but also their individual damage severity by carrying out solely one analysis.