• Structural Engineering and Mechanics
• /
• 제42권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.

• 한국자기학회지
• /
• 제7권2호
• /
• pp.109-116
• /
• 1997

교번 히스테리시스와 회전 히스테리시스를 동시에 고려할 수 있는 벡터 히스테리시스 모델을 이용하여, 회전 히스테리시스 특성을 갖는 전자계의 유한 요소 해석 알고리즘을 제안한다. 벡터 히스테리시스 모델은 자화의존 프라이자흐 모델에서 확장된 모델을 이용하였다. 회전 자계가 자성체에 가해진 경우, 각 요소에서의 자화값이 벡터 모델로 부터 구해지며, 이 값은 다시 유한요소 해석의 입력이 되어서 인가자계와 자화와의 지연각을 구할 수 있다. 실험 결과와의 비교를 통해, 제안된 방법이 회저자계의 영향을 받는 전자계 시스템의 자화 특성을 정확히 해석 할 수 있음을 확인하였다.

• Journal of Power Electronics
• /
• 제11권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
• /
• 제15권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.

• 한국공작기계학회논문집
• /
• 제10권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.

• 한국전산구조공학회:학술대회논문집
• /
• 한국전산구조공학회 1998년도 가을 학술발표회 논문집
• /
• 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.

• 대한전기학회:학술대회논문집
• /
• 대한전기학회 1987년도 전기.전자공학 학술대회 논문집(I)
• /
• pp.589-592
• /
• 1987

• Structural Engineering and Mechanics
• /
• 제40권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.

• 대한기계학회논문집
• /
• 제17권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
• /
• 제75권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.