• Proceedings of the KIEE Conference
• /
• 2000.11b
• /
• pp.315-317
• /
• 2000

Linear induction motor(LIM) have static and dynamic end effects due to its finite core length, so that per-phase impedances are asymmetric and the air gap flux distribution is distorted. So, this paper propose the d-q axis equivalent circuit and vector control method considering both static and dynamic end effects of the LIM. This vector control method consists of the slip frequency control, the time-invariant control and decoupling control. As a result, it is shown that the results of equivalent circuit method have a goof agreement with the results of finite element method.

• Proceedings of the Korea Water Resources Association Conference
• /
• 2005.05b
• /
• pp.409-413
• /
• 2005

Even though the relative importance of length scale of flow system allow us to simplify three dimensional flow problem to one or two dimensional representation, many systems still require three dimensional analysis. The objective of this study is to develop an efficient and accurate finite element model for analyzing and predicting three dimensional flow features in natural rivers and to offend to model spreading of pollutants and transport of sediments in the future. Firstly, three dimensional Reynolds averaged Navier-Stokes equations with the hydrostatic pressure assumption in generalized curvilinear coordinates were combined with the kinematic free-surface condition. Secondly. to simulate realistic high Reynolds number flow, the model employed the Streamline Upwind/Petrov-Galerkin(SU/PG) scheme as a weighting function for the finite element method in conjunction with an appropriate turbulence model(Smagorinsky scheme for the horizontal plain and Mellor-Yamada scheme for the vertical direction). Several tests is performed for the purpose of validation and verification of the developed model. A simple rectangular channel, 5-shaped and U-shaped channel are used for tests and comparisons are made with RMA-10 model. Runs for each case is converged stably without a oscillation and calculated water-surface deformation, longitudinal and transversal velocities, and velocity vector fields are in good agreement with the results of RMA-10 model.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.11a
• /
• pp.330-336
• /
• 2005

This paper presents a method to analyze the vibration of a flexible spinning disk-spindle system with FDBs, flexible base structure and an actuator in a HDD by using the FEM. Finite element equations of each component of a HDD spindle system from the spinning flexible disk to the flexible base plate are consistently derived by satisfying the geometric compatibility in the internal boundary between each component. A global matrix equation obtained by assembling the finite element equations of each substructure is transformed to a state-space matrix-vector equation, and both damped natural frequencies and modal damping ratios are calculated by using the restarted Arnoldi iteration method. The validity of the proposed method is verified by comparing the simulated natural frequencies, mode shapes with the experimental results.

• Structural Engineering and Mechanics
• /
• v.13 no.2
• /
• pp.135-154
• /
• 2002

This paper presents the convected material frame approach to study the nonlinear behavior of inelastic frame structures. The convected material frame approach is a modification of the co-rotational approximation by incorporating an adaptive convected material frame in the basic definition of the displacement vector and strain tensor. In the formulation, each discrete element is associated with a local coordinate system that rotates and translates with the element. For each load increment, the corresponding strain-displacement and nodal force-stress relationships are defined in the updated local coordinates, and based on the updated element geometry. The rigid body motion and deformation displacements are decoupled for each increment. This modified approach incorporates the geometrical nonlinearities through the continuous updating of the material frame geometry. A generalized nonlinear function is used to derive the inelastic constitutive relation and the kinematic hardening is considered. The equation of motion is integrated by an explicit procedure and it involves only vector assemblage and vector storage in the analysis by assuming a lumped mass matrix of diagonal form. Several numerical examples are demonstrated in close agreement with the solutions obtained by the ANSYS code. Numerical studies show that the proposed approach is capable of investigating large deflection of inelastic planar structures and providing an excellent numerical performance.

• 한국전산유체공학회:학술대회논문집
• /
• 1997.10a
• /
• pp.49-54
• /
• 1997

Rice's monotone streamline upwind finite element method, which was proposed to treat convection-dominated flows, is applied to the linear triangular element. An alignment technique of unstructured grids with given velocity fields is used to prevent the interpolation error produced in evaluating the convection term in the upwind method. The alignment of grids is accomplished by optimizing a target function defined with the inner-product of a properly chosen side vector in the element with the velocity field. Two pure advection problems are considered to demonstrate the superiorities of the present approach in solving the convection-dominated flow on the unstructured grid. Solutions obtained with aligned grids are much closer to the exact solutions than those with initial regular grids. The capability of the present approach in predicting the appearance of the secondary vortex in the laminar confined jet impingement is shown by comparing streamlines to those produced by SIMPLE on a highly stretched grid toward the impingement plate.

• Restorative Dentistry and Endodontics
• /
• v.19 no.1
• /
• pp.231-254
• /
• 1994

Fracture of cusp, on posterior teeth, especially those carious or restored, is major cause of tooth loss. Inappropriate treatments, such as unnecessarily wide cavity preparations, increase the potential of further trauma and possible fracture of the remaining tooth structures. Fracture potential may be directly related to the stresses exerted upon the tooth during masticatory function. The purpose of this study is to evaluate the fracture resistance of tooth, restored with composite resin inlay. In this study, MOD inlay cavity prepared on maxillary first premolar and restored with composite resin inlay. Three dimensional finite element models with eight nodes isoparametric solid element, developed by serial grinding-photographing technique. These models have various occlusal isthmus and depth of cavity, 1/2, 1/3 and 1/4 of isthmus width and 0.7, 0.85 and 1.0 of depth of cavity. The magnitude of load was 474 N and 172 N as presented to maximal biting force and normal chewing force. These loads applied onto ridges of buccal and lingual cusp. These models analyzed with three dimensional finite element method. The results of this study were as follows : 1. There is no difference of displacement between width of occlusal isthmus and depth of cavity. 2. The stress concentrated at bucco-mesial comer, bucco-disal comer, pulpal line angle and the interface area between internal slopes of cusp and resin inlay. 3. The vector of stress direct to buccal and lingual side from center of cavity, to tooth surface going on to enamel. The magnitude of vector increase from occlusal surface to cervix. 4. The crack of tooth start interface area, between internal slop of buccal cusp and resin inlay. It progresses through buccopulpal line angle to cervix at buccomesial and buccodistal comer. 5. The influence with depth of cavity to fracture of tooth was more than width of isthmus. 6. It would be favorable to make the isthmus width narrower than a third of the intercuspal distance and depth of cavity is below 1 : 0.7.

• Structural Engineering and Mechanics
• /
• v.74 no.1
• /
• pp.1-18
• /
• 2020

For the large deformation of shear walls under vertical and horizontal loads, there are difficulties in obtaining accurate simulation results using the response analysis method, even with fine mesh elements. Furthermore, concrete material nonlinearity, stiffness degradation, concrete cracking and crushing, and steel bar damage may occur during the large deformation of reinforced concrete (RC) shear walls. Matrix operations that are involved in nonlinear analysis using the traditional finite-element method (FEM) may also result in flaws, and may thus lead to serious errors. To solve these problems, a planar four-node element was developed based on vector mechanics. Owing to particle-based formulation along the path element, the method does not require repeated constructions of a global stiffness matrix for the nonlinear behavior of the structure. The nonlinear concrete constitutive model and bilinear steel material model are integrated with the developed element, to ensure that large deformation and damage behavior can be addressed. For verification, simulation analyses were performed to obtain experimental results on an RC shear wall subjected to a monotonically increasing lateral load with a constant vertical load. To appropriately evaluate the parameters, investigations were conducted on the loading speed, meshing dimension, and the damping factor, because vector mechanics is based on the equation of motion. The static problem was then verified to obtain a stable solution by employing a balanced equation of motion. Using the parameters obtained, the simulated pushover response, including the bearing capacity, deformation ability, curvature development, and energy dissipation, were found to be in accordance with the experimental observation. This study demonstrated the potential of the developed planar element for simulating the entire process of large deformation and damage behavior in RC shear walls.

• Computers and Concrete
• /
• v.2 no.5
• /
• pp.345-366
• /
• 2005

This work presents an assessment of the computational performance of a vector-parallel implementation of probabilistic model for concrete cracking in 3D. This paper shows the continuing efforts towards code optimization as reported in earlier works Paz, et al. (2002a,b and 2003). The probabilistic crack approach is based on the direct Monte Carlo method. Cracking is accounted by means of 3D interface elements. This approach considers that all nonlinearities are restricted to interface elements modeling cracks. The heterogeneity governs the overall cracking behavior and related size effects on concrete fracture. Computational kernels in the implementation are the inexact Newton iterative driver to solve the non-linear problem and a preconditioned conjugate gradient (PCG) driver to solve linearized equations, using an element by element (EBE) strategy to compute matrix-vector products. In particular the paper analyzes code behavior using OpenMP directives in parallel vector processors (PVP), such as the CRAY SV1 and CRAY T94. The impact of the memory architecture on code performance, and also some strategies devised to circumvent this issue are addressed by numerical experiment.

• Journal of KSNVE
• /
• v.3 no.1
• /
• pp.77-82
• /
• 1993

In general, the dynamic analysis with FEM(Finite Element Method) of large structures requires large computer memory space and long computational time. For the purpose of economical dynamic analysis of large structures, most of engineers want to use an efficient solution algorithm. This paper reports the modified CMS(Component Mode Synthesis) method which uses more efficient algorithm than the classical CMS method. In this paper, it is shown that Ritz vector sets can play the role of normal mode vector sets of substurctures in the CMS algorithm. The modified CMS method has good convergence performance compared with that of the classical CMS method.

• Architectural research
• /
• v.13 no.1
• /
• pp.17-24
• /
• 2011

Isogeometric solution of Poisson equation is provided. NURBS (NonUniform B-spline Surface) is introduced to express both geometry of structure and unknown field of governing equation. The terms of stiffness matrix and load vector are consistently derived with very accurate geometric definition. The validity of the isogeometric formulation is demonstrated by using two numerical examples such as square plate and L-shape plate. From numerical results, the present solutions have a good agreement with analytical and finite element (FE) solutions with the use of a few cells in isogeometric analysis.