• Proceedings of the KSME Conference
• /
• 2000.11a
• /
• pp.529-534
• /
• 2000

A finite element analysis for a rotating cantilever beam is presented in this study. Based on a dynamic modelling method using the stretch deformation instead of the conventional axial deformation, three linear partial differential equations are derived from Hamilton's principle. Two of the linear differential equations show the coupling effect between stretch and chordwise deformations. The other equation is an uncoupled one for the flapwise deformation. From these partial differential equations and the associated boundary conditions, are derived two weak forms: one is for the chordwise motion and the other is for the flapwise motion. The weak forms are spatially discretized with newly defined two-node beam elements. With the discretized equations or the matrix-vector equations, the behaviours of the natural frequencies are investigated for the variation of the rotating speed.

• Structural Engineering and Mechanics
• /
• v.53 no.2
• /
• pp.205-226
• /
• 2015

Finite element method (FEM) is an effective quantitative method to solve complex engineering problems. The basic idea of FEM for a complex problem is to be able to find a solution by reducing the problem made simple. If mathematical tools are inadequate to obtain precise result, even approximate result, FEM is the only method that can be used for structural analyses. In FEM, the domain is divided into a large number of simple, small and interconnected sub-regions called finite elements. FEM has been used commonly for linear and nonlinear analyses of different types of structures to give us accurate results of plane stress and plane strain problems in civil engineering area. In this paper, FEM is used to investigate stress analysis of a shear wall which is subjected to concentrated loads and fundamental principles of stress analysis of the shear wall are presented by using matrix displacement method in this paper. This study is consisting of two parts. In the first part, the shear wall is discretized with constant strain triangular finite elements and stiffness matrix and load vector which is attained from external effects are calculated for each of finite elements using matrix displacement method. As to second part of the study, finite element analysis of the shear wall is made by ANSYS software program. Results obtained in the second part are presented with tables and graphics, also results of each part is compared with each other, so the performance of the matrix displacement method is demonstrated. The solutions obtained by using the proposed method show excellent agreements with the results of ANSYS. The results show that this method is effective and preferable for the stress analysis of shell structures. Further studies should be carried out to be able to prove the efficiency of the matrix displacement method on the solution of plane stress problems using different types of structures.

• Advances in Computational Design
• /
• v.5 no.3
• /
• pp.305-321
• /
• 2020

In this paper is presented the solution method for three-dimensional problem of transversely isotropic body's elastoplastic deformation by the finite element method (FEM). The process of problem solution consists of: determining the effective parameters of a transversely isotropic medium; construction of the finite element mesh of the body configuration, including the determination of the local minimum value of the tape width of non-zero coefficients of equation systems by using of front method; constructing of the stiffness matrix coefficients and load vector node components of the equation for an individual finite element's state according to the theory of small elastoplastic deformations for a transversely isotropic medium; the formation of a resolving symmetric-tape system of equations by summing of all state equations coefficients summing of all finite elements; solution of the system of symmetric-tape equations systems by means of the square root method; calculation of the body's elastoplastic stress-strain state by performing the iterative process of the initial stress method. For each problem solution stage, effective computational algorithms have been developed that reduce computational operations number by modifying existing solution methods and taking into account the matrix coefficients structure. As an example it is given, the problem solution of fibrous composite straining in the form of a rectangle with a system of circular holes.

• Journal of KSNVE
• /
• v.8 no.5
• /
• pp.949-956
• /
• 1998

Complex and large lattice type structures are frequently used in design of bridge, tower, crane and aerospace structures. In general, in order to analyze these structures we have used the finite element method(FEM). This method is the most widely used and powerful method for structural analysis lately. However, it is necessary to use a large amount of computer memory and computational time because the FEM requires many degrees of freedom for solving dynamic problems exactly for these complex and large structures. For analyzing these structures on a personal computer, the authors developed the transfer stiffness coefficient method(TSCM). This method is based on the concept of the transfer of the nodal dynamic stiffness coefficient matrix which is related to force and displacement vector at each node. And we suggested TSCM for free vibration analysis of complex and large lattice type structures in the previous report. In this paper, we formulate forced vibration analysis algorithm for complex and large lattice type structures using extened TSCM. And we confirmed the validity of TSCM through computational results by the FEM and TSCM, and experimental results for lattice type structures with harmonic excitation.

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

• Proceedings of the KIEE Conference
• /
• 1995.07a
• /
• pp.3-5
• /
• 1995

For a linear induction motor(LIM), the constants of each phase are different due to the structure. In this paper, a vector control analysis method of a LIM taking into consideration its asymmetrical constants are proposed. And, in order to prove the propriety of proposed vector control method and to analyze the dynamic characteristics of LIM's vector control, FEM taking into account of movement and using stator tapped winding is used in the analysis region. So, It is confirmed that the proposed asymetrical constants vector control theory and simulation method of mixing with FEM is appropriate to dynamic characteristics analysis of LIM.

• 한국전산유체공학회:학술대회논문집
• /
• 2006.10a
• /
• pp.61-65
• /
• 2006

This paper describes a recent development on the divergence free basis function based on a hermite stream function. The well-known cavity problem has been used to compare the accuracy and the convergence of the present method with those of a modified residual method known as one of the stabilized finite element methods. The comparison showed the present method performs better in the accuracy and convergence.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.33 no.5
• /
• pp.167-172
• /
• 1984

This paper describes the finite element analysis of magnetostatic field problems with permanent magnets. Two kinds of algorithms, one using the magnetic vector potential and the other using the magnetic scalar potential, are introduced. The magnetization of the pemanent magnet is used as the source instead of the magnetic equivalent current in both of the formulations using the magnetic vector potential and the magnetic scalar potential. A simple functional, which has only the region integral instead of the region integral and boundary integral, is derived in the formulation using the magnetic scalar potential. These make the formulation of the system equations simpler and more convenient than the conventional methods. The numerical results by the two proposed algorithms for a C-type permanent magnet model are compared with the analytic solutions respectively. The numerical results are in good agreement with the analytic solutions.

• The Transactions of the Korean Institute of Electrical Engineers B
• /
• v.54 no.10
• /
• pp.483-491
• /
• 2005

To solve the 3D eddy current problems by using FE(finite element) method with MVP(magnetic vector potential) and electric scalar potential, Coulomb gauge condition and current continuity condition have to be considered. Coulomb gauge condition enforced on existing FE formulations to insure the uniqueness of MVP looks unnatural and current continuity condition which can be driven from Ampere's law looks unnecessary. So in this paper the effect of two conditions on FE formulations are investigated in order to help to obtain accurate numerical simulation results.

• Proceedings of the KIEE Conference
• /
• 1994.07b
• /
• pp.1213-1215
• /
• 1994

The shape of lossy material inside a waveguide matched load is optimally designed to give low reflection over a given frequency range. The 3 dimensional vector finite element mettled is used as an analysis tool which does not generate spurious mode. The optimizing process used in this parer is the Powell technique. The designed load gives the low reflection about -30 dB around 10GHz with 1.5 wavelength of the load length.