• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1999.04a
• /
• pp.42-49
• /
• 1999

Because the linear elastic tincture analysis has been proved to be insufficient in predicting the failure of cracked bodies, in recent years, a number of fracture concepts have been studied which remain applicable in the presence of large-scale plasticity near a crack tip. This work thereby presents a new finite element model, as accurate as possible, to analyze plane problems of ductile fracture under large-scale yielding conditions. Based on the incremental theory of plasticity, the p-version finite element analysis is employed to account for the values of J-integral, the most dominant fracture parameter, and the shape of plastic zone near a crack tip by using the J-integral method and equivalent domain integral method. The numerical results by the proposed model are compared with the theoretical solutions in literatures and the numerical solutions by the i,-version of F.E.M.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1999.04a
• /
• pp.169-176
• /
• 1999

An improved formulation for spatial stability md free vibration of thin-walled curved beams with variable curvature and non-symmetric cross sections are presented based on the displacement field considering the second order terms of finite semitangential rotations. By introducing Vlasov's assumptions, the total potential energy is derived from the principle of linearized virtual work for a continuum. In this formulation, all displacement parameters and the warping function are defined at the centroid axis so that the coupled terms of bending and torsion are added to the elastic strain energy. Also, the potential energy due to initial stress resultants is consistently derived corresponding to the semitangential rotation and moment. The cubic Hermitian polynomials are utilized as shape functions for development of the curved thin-walled beam element having eight degrees of freedom. In order to illustrate the accuracy and practical usefulness of this study, . numerical solutions for free vibration of arches are presented and compared with resells of other researchers and solutions analyzed by the ABAQUS's shell element.

• Journal of computational fluids engineering
• /
• v.9 no.4
• /
• pp.1-6
• /
• 2004

This study investigates the transition of flows in a natural convection problem with periodic wall temperatures of the form, T/sub L/=T₁＋δ Tsinκχ and T/sub L/=T₂＋δ Tsinκχ .The fluid considered is air with P/sub γ/=0.7. In the conduction-dominated regime with a small Rayleigh number, two large cells are formed over one wave length, for all wave numbers. When k≤1.8, the flow becomes unstable with increase of the Rayleigh number, and multicellular convection occurs above a critical Rayleigh number. The flow patterns are classified by the number of eddies over one wave length, and several kinds of transition phenomena, such as 2→3→4, 4→3→2, and 2→4 eddy flow, occur with increase( or decrease) of the Rayleigh number. Dual solutions are found above a critical Rayleigh number, and an anomalous bifurcation is observed.

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

A finite element procedure to estimate ultimate strength of space frames considering spread of plasticity is presented. The improved displacement field is introduced based on inclusion of second order terms of finite rotations. All the nonlinear terms due to bending and torsional moment as well as axial force are precisely considered. The concept of plastic hinge is introduced and the incremental load/displacement method is applied for the elasto-plastic analysis. The initial yield surface is defined based on the residual stress and the full plastification surface is considered under the combined action of axial force, bending and torsional moments. The elasto-plastic stiffness matrices are derived using the flow rule and the normality condition of the limit function. Finite element solutions for ultimate strength of space frames are compared with available solutions and experimental results.

• Computational Structural Engineering
• /
• v.1 no.2
• /
• pp.83-90
• /
• 1988

The present work is concerned with the improvement of finite element for the analysis of plate bending structures. The element formulation is based upon Mindlin plate concept. The displacement field of this element is formed by adding nonconforming modes to two rotational displacement components of a 'heterosis plate element. The element has the requisite numbers of zero eigenvalues associated with rigid body modes to avoid the spurious zero energy mode. It is shown that the results obtained by the element converged to the exact solutions very rapidly as the mesh is refined and exhibited reliable solutions through numerical studies for standard benchmark problems. This element is shown to overcome the shear locking problem completely in very thin plate situation even for irregular meshes.

• Journal of computational fluids engineering
• /
• v.4 no.2
• /
• pp.1-8
• /
• 1999

Natural convection in a horizontal annulus with the inner cylinder heated by the application of a constant heat flux and the isothermally cooled outer cylinder is considered, and the transition of flows and the bifurcation phenomenon are numerically investigated for air with Pr=0.7. The zero initial condition always induces a crescent-sheped eddy flow. A bicellular flow in which the fluid descends along the vertical central plane of the annulus can be obtained at high Rayleigh number by introducing artificial numerical disturbances. Dual solutions are found above a certain critical Rayleigh number. Hysteresis phenomena have not been observed.

• Journal of the Korea Society for Simulation
• /
• v.11 no.3
• /
• pp.13-22
• /
• 2002

This paper describes an automated evaluation of electrostatic field for micro motors whose sizes range 10 to 103um. Electric field modeling in micro motors has been generally restricted to in-plane two-dimensional finite element analysis (FEA). In this paper, the actual three-dimensional geometry of the micro motor is considered. An automatic FE mesh generation technique, which is based on the fuzzy knowledge processing and computational geometry techniques, is incorporated in the system, together with one of commercial FE analysis codes and one of commercial solid modelers. The system allows a geometry model of concern to be automatically converted to different FE models, depending on physical phenomena to be analyzed, electrostatic analysis and stress analysis and so on. The FE models are then exported to the FE analysis code, and then analyses are peformed. Then, analytical analysis and FE analysis about the torque generated by electrostatic micro motor are performed. The starting torque is proportional to $V^2$, the calculated starting torque from the two-dimensional analytical solutions are three times larger than those from the three-dimensional FE solutions.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.449-456
• /
• 1998

A consistent finite element formation and analytic solutions are presented for spatial stability of thin-walled circular arch. The total potential energy is derived by applying the principle of linearized virtual work and including second order terms of finite semitangential rotations. As a result the energy functional corresponding to the semitangential rotation is obtained, in which the elastic strain energy terms are considered restrained warping effects. We have obtained analytic solution for the lateral buckling of monosymmetric thin-walled curved beam subjected to pure bending or uniform compression and it's boundary conditions are simply supported. For finite element analysis, the two node cubic Hermitian polynomials are utilized as shape Auctions. In order to illustrate the accuracy of this study, parameter studies for lateral buckling problems of circular arch are presented and compared with available solutions and numerical results analyzed by the FEM using straight beam element.

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

A stability problems of isotropic shells under pure bending is investigated based on the classical shells theory. The governing equations of stability problem presented by Donnell and Love, are developed and the solutions for the cylindrical shells are obtained by using Galerkin method. Bending moment is applied at the ends of the cylindrical shell as a from of distributed load in the shape of sine curve. For the isotropic materials, the result of the general purpose structural analysis program based on the finite element method are compared with the critical moment obtained from the classical shell theories. The critical loads for the cylindrical shells with various geometry can not be evaluated with a simple equation. However, accurate solutions for the stability problems of cylindrical shells can be obtained through the equilibrium equation developed in the study.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1997.10a
• /
• pp.239-246
• /
• 1997

The general formulation of free vibration and stability analysis of unsymmetric thin-walled space frames and beam-columns is presented. The kinetic and total potential energy is derived by applying the extended virtual work principle, introducing displacement parameters defined at the arbitrarily chosen axis and including second order terms of finite semitangential rotations. In formulating the finite element procedure, cubic Hermitian polynomials are utilized as shape functions of the two node space frame element. Mass, elastic stiffness, and geometric stiffness matrices for the unsymmetric thin-walled section are evaluated. In order to illustrate the accuracy and practical usefulness of this formulation, finite element solutions for the free vibration and stability problems of thin-walled beam-columns and space frames are presented and compared with available solutions.