• Journal of Advanced Marine Engineering and Technology
• /
• v.14 no.4
• /
• pp.46-56
• /
• 1990

Stress analysis of axi-symetric body with non-symetric loading and non-symetric given displacements is investigated in this paper using the finite element method. As the non-symetric load and non-symetric given displacements of axi-symetric body are generally periodic functions of angle .theta., the nodal forces and nodal displacements can be expanded in cosine and sine series, that is, Fourier series. Furthermore, using Euler's formula, the cosine and sine series can be converted into exponential series and it is prooved that the related calculus become more clear. Substituting the nodal displacements expanded in Fourier series into the strain components of cylindrical coordinates system, the element strains are expressed in series form and by the principal of virtual work, the element stiffness martix and element load vector are obtained for each order. It is also showed that if the non-symetric loads are even or odd functions of angle ${\theta}$ the stiffness matrix and load vector of the system are composed with only real numbers and relatively small capacity fo computer memory is enough for calculation.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2007.04a
• /
• pp.577-582
• /
• 2007

A variational formulation for plane elasticity problems is derived based on an isogeometric approach. The isogeometric analysis is an emerging methodology such that the basis functions in analysis domain arc generated directly from NURBS (Non-Uniform Rational B-Splines) geometry. Thus. the solution space can be represented in terms of the same functions to represent the geometry. The coefficients of basis functions or the control variables play the role of degrees-of-freedom. Furthermore, due to h-. p-, and k-refinement schemes, the high order geometric features can be described exactly and easily without tedious re-meshing process. The isogeometric sensitivity analysis method enables us to analyze arbitrarily shaped structures without re-meshing. Also, it provides a precise construction method of finite element model to exactly represent geometry using B-spline base functions in CAD geometric modeling. To obtain precise shape sensitivity, the normal and curvature of boundary should be taken into account in the shape sensitivity expressions. However, in conventional finite element methods, the normal information is inaccurate and the curvature is generally missing due to the use of linear interpolation functions. A continuum-based adjoint sensitivity analysis method using the isogeometric approach is derived for the plane elasticity problems. The conventional shape optimization using the finite element method has some difficulties in the parameterization of boundary. In isogeometric analysis, however, the geometric properties arc already embedded in the B-spline shape functions and control points. The perturbation of control points in isogeometric analysis automatically results in shape changes. Using the conventional finite clement method, the inter-element continuity of the design space is not guaranteed so that the normal vector and curvature arc not accurate enough. On tile other hand, in isogeometric analysis, these values arc continuous over the whole design space so that accurate shape sensitivity can be obtained. Through numerical examples, the developed isogeometric sensitivity analysis method is verified to show excellent agreement with finite difference sensitivity.

• Journal of the Korean Society for Precision Engineering
• /
• v.16 no.6
• /
• pp.148-159
• /
• 1999

A numerical method is presented for the dynamic analysis of spur gears rotating with very high angular speeds. For an efficient computation each gear is assumed to consist of a rotating rigid disk and an elastic tooth having mass, and finite element formulations are used for the equations of motion of the tooth. The geometric constraint is imposed between the rigid disk and the elastic tooth to fix them, and contact condition is imposed between the meshing teeth of the gears. At each iteration of each time step the Lagrange multiplier and contact force are revised by using the constraint error vector, and then the whole equations of motion are time integrated with the given Lagrange multiplier and contact force. For the accurate solution the velocity and acceleration constraints as well as the displacement constraint are satisfied by the monotone reductions of the constraint error vectors. Computing procedures associated with the iterative schemes are explained and numerical simulations are conducted with the spur gears.

• Proceedings of the Korean Institute of Electrical and Electronic Material Engineers Conference
• /
• 1998.06a
• /
• pp.33-36
• /
• 1998

This paper is a study on a linear ultrasonic motor with a first longitudinal$(L_1)$ and fourth bending $(B_4)$ double-mode rectangular plate. The stator vibrator is composed of an elastic material plate and of a piezo-ceramic element having a motion by electrical excitation. Each strain vector differs by $90^{\circ}$ generate travelling wave with the elliptical displacement motion of a point on the surface. To magnify displacement of longitudinal direction in elliptical displacement motion, the motor has a mechanism of the.displacement enlargement. In this paper, the vibration shape of the stator is simulated using the finite element method. A detailed model considered of the piezoelectric effect and of the exact geometry of the stator is used to calculate the displacement. The position of displacement mechanism is decided by the maximum displacement.

• Proceedings of the KIEE Conference
• /
• 2000.11b
• /
• pp.262-264
• /
• 2000

A numerical method for the analysis of 3D eddy current in conductors due to applied time varying field is suggested using the finite element method. In the approximation of the field quantifies, the edge element is used, because it reduce the required computer memory and the computing time compared with the nodal elements. With edge elements, furthermore, the field governing equations become simple because the electric scalar potential ${\phi}$ can be set to zero. The modified magnetic vector potential($A^*$) is used as a state variable. The analysed results are compared with the experimentally measured ones for the TEAM workshop problem3.

• The Transactions of the Korean Institute of Electrical Engineers
• /
• v.41 no.7
• /
• pp.723-734
• /
• 1992

In this paper, the characteristics of IPMSM(Interior-type Permanent Magnet Synchronous Motor) are simulated using 2-D. finite element method. This paper deals with the following characteristics : air gap flux density considering skew, back e.m.f., torque and inductance. Back e.m.f. is calculated using the flux obtained from the vector potential of FEM solution. Torque is calculated using improved Maxwell stress tensor method and current angle which is obtained from the controller. Direct axis inductance and quadrature axis inductance are also calculated using energy perturbation method. Computed results are found in satisfactory agreement with experimental ones. This method also can be applied for the computation and analysis of the characteristics of SPMSM, current-excited synchronous motor and reluctance motor.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.04a
• /
• pp.288-295
• /
• 2003

An adaptive procedure in finite element analysis is presented by p-refinement of meshes in conjunction with a posteriori error estimator that is based on the recovery technique. In case of the recovery technique, the SPR(superconvergent patch recovery) approach has been modified for p-adaptive mesh refinement. The strategy of finding a nearly optimal distribution of polynomial degrees on a fixed finite element mesh is discussed such that a particular element has to be refined automatically to obtain an acceptable level of accuracy by increasing p-levels non-uniformly. To verify the proposed algorithm, the limit value approach is proposed which utilizes the exact strain energy computed from the extrapolation equation. A new pre-processor is developed for the p-version finite element program in which the vector graphic editor is used for the automatic generation of node connection and coordinate by halfedge solid data structure according to uniform or nonuniform p-distribution. The general 2-D algorithm is also developed to generate face modes and internal modes in accordance with different mesh types. The quality of the error estimator is investigated with the help of two mumerical examples. The results show that the sequences of p-distributions obtained by the proposed error indicator closely follow the optimal trajectory.

• Structural Engineering and Mechanics
• /
• v.81 no.4
• /
• pp.411-428
• /
• 2022

This paper presents the Campbell diagram analysis of the rotordynamic system using the full order model (FOM) and the reduced order model (ROM) techniques to determine the critical speeds, identify the stability and reduce the computational time. Due to the spin-speed-dependent matrices (e.g., centrifugal stiffening matrix), several model order reduction (MOR) techniques may be considered, such as the modal superposition (MS) method and the Krylov subspace-based MOR techniques (e.g., Ritz vector (RV), quasi-static Ritz vector (QSRV), multifrequency quasi-static Ritz vector (MQSRV), multifrequency/ multi-spin-speed quasi-static Ritz vector (MMQSRV) and the combined Ritz vector & modal superposition (RV+MS) methods). The proposed MMQSRV method in this study is extended from the MQSRV method by incorporating the rotational-speed-dependent stiffness matrices into the Krylov subspace during the MOR process. Thus, the objective of this note is to respond to the question of whether to use the MS method or the Krylov subspace-based MOR technique in establishing the Campbell diagram of the shaft-disc-blade assembly systems in three-dimensional (3D) finite element analysis (FEA). The Campbell diagrams produced by the FOM and various MOR methods are presented and discussed thoroughly by computing the norm of relative errors (ER). It is found that the RV and the MS methods are dominant at low and high rotating speeds, respectively. More precisely, as the spinning velocity becomes large, the calculated ER produced by the RV method is significantly increased; in contrast, the ER produced by the MS method is smaller and more consistent. From a computational point of view, the MORs have substantially reduced the time computing considerably compared to the FOM. Additionally, the verification of the 3D FE rotordynamic model is also provided and found to be in close agreement with the existing solutions.

• Earthquakes and Structures
• /
• v.17 no.1
• /
• pp.63-73
• /
• 2019

This paper studies damage detection as an optimization problem. A new objective function based on changes in natural frequencies, and Natural Frequency Vector Assurance Criterion (NFVAC) was developed. Due to their easy and fast acquisition, natural frequencies were utilized to detect structural damages. Moreover, they are sensitive to stiffness reduction. The method presented here consists of two stages. Firstly, Finite Element Model (FEM) is updated. Secondly, damage severities and locations are determined. To minimize the proposed objective function, a new bio-inspired optimization algorithm called salp swarm was employed. Efficiency of the method presented here is validated by three experimental examples. The first example relates to three-story shear frame with two single damage cases in the first story. The second relates to a five-story shear frame with single and multiple damage cases in the first and third stories. The last one relates to a large-scale eight-story shear frame with minor damage case in the first and third stories. Moreover, the performance of Salp Swarm Algorithm (SSA) was compared with Particle Swarm Optimization (PSO). The results show that better accuracy is obtained using SSA than using PSO. The obtained results clearly indicate that the proposed method can be used to determine accurately and efficiently both damage location and severity in multi-story shear frames.

• Journal of the Korean Magnetics Society
• /
• v.7 no.3
• /
• pp.159-165
• /
• 1997

Recently, many researchers are interested in the analysis of the magnetic fields considering the vector relation between the magnetic flux density and the magnetic field intensity. In this paper, Finite Element Method with two dimensional tensor permeability is adopted to calculate the magnetic field quantities in the system having the rotational magnetic field. In addition, the method of analyzing the magnetic field considering the harmonic components is introduced. The validity of the proposed method is proved compared with the experimental results.