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

This paper deals with the development of a variable-node element and its application to the adaptive h-version mesh refinement-recovery for the incompressible viscous flow analysis. The element which has variable mid-side nodes can be used in generating the transition zone between the refined and unrefined elements and efficiently used for construction of a refined mesh without generating distorted elements. A modified Gaussian quadrature is needed to evaluate the element matrices due to the discontinuity of derivatives of the shape functions used for the element. The penalty function method which can reduce the number of independent variables is adopted for the purpose of computational efficiency and the selective reduced integration is carried out for the convection and pressure terms to preserve the stability of solution. For the economical analysis of transient problems, not only the mesh refinement but also the mesh recovery is needed. The numerical examples show that the optimal mesh for the finite element analysis of a wind around the structures can be obtained automatically by the proposed scheme.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1990.04a
• /
• pp.18-23
• /
• 1990

The paper is concerned with the elasto-plastic and geometrically nonlinear analysis of shell structures using an improved degenerated shell element. In the formulation of the improved degenerated shell element, an enhanced interpolation of transverse shear strains in the natural coordinate system is used to overcome the shear locking problems; the reduced integration technique in in-plane strains is applied to avoid membrane locking behavior; selective addition the nonconforming displacement modes improve the element performances. This element is free of serious locking problems and undesirable compatible or commutable spurious kinematic deformation modes and passes the patch tests. An incremental total Lagrangian formulation is presented which allows the calculation of arbitrarily large displacements and rotations. The resulting nonlinear equations are solved by the Newton-Raphson solution scheme. The versatility and accuracy of this improved degenerated shell element are demonstrated by solving several numerical examples.

• Journal of computational fluids engineering
• /
• v.2 no.2
• /
• pp.26-34
• /
• 1997

An implicit viscous turbulent flow solver is developed for two-dimensional geon unstructured triangular meshes. The flux terms are discretized based on a cell-centered formulation with the Roe's flux-difference splitting. The solution is advanced in time us backward-Euler time-stepping scheme. At each time step, the linear system of equation approximately solved wi th the Gauss-Seidel relaxation scheme. The effect of turbulence is with a standard k-ε two-equation model which is solved separately from the mean flow equation the same backward-Euler time integration scheme. The triangular meshes are generated advancing-front/layer technique. Validations are made for flows over the NACA 0012 airfoil. Douglas 3-element airfoil. Good agreements are obtained between the numerical result experiment.

• Computers and Concrete
• /
• v.25 no.2
• /
• pp.181-192
• /
• 2020

Nowadays, steel fiber reinforced concretes (SFRCs) are widely used in practical applications. Significant experimental research has thus been carried out to determine the constitutive equations that represent the behavior of SFRCs under multiaxial loadings. However, numerical modelling of SFRCs via FEM has been challenging due to the complexities of the implementation of these constitutive equations. In this study, following the literature, a plasticity model is constructed for the behavior of SFRCs that involves the Willam-Warnke failure surface with the relevant evolution laws and a non-associated flow rule for determining the plastic deformations. For the precise (yet rapid) integration of the constitutive equations, an explicit substepping scheme consisting of yield intersection and drift correction algorithms is employed and thus implemented in ABAQUS via UMAT. The FEM model includes various material parameters that are determined from the experimental data. Three sets of parameters are used in the numerical simulations. While the first set is from the experiments that are conducted in this study on SFRC specimens with various contents of steel fibers, the other two sets are from the experiments reported in the literature. The response of SFRCs under multiaxial compression obtained from various numerical simulations are compared with the experimental data. The good agreement between numerical results and the experimental data indicates that not only the adopted plasticity model represents the behavior of SFRCs very well but also the implemented integration scheme can be employed in practical applications of SFRCs.

• Journal of KSNVE
• /
• v.10 no.1
• /
• pp.97-106
• /
• 2000

The dynamic behavior of crankshaft-bearing system in scroll compressor has been investigated using the combined methodologies of finite elements and transfer matrices. The finite element formulation is proposed including the field element for a shaft section and the point element at balancer weight locations, bearing locations, etc., whereas the conventional method is used with the elements. The Houbolt method is used to consider the time march for the integration of the system equations. The linear stiffness and damping coefficients are calculated for a finite cylindrical fluid-film bearing by solving the Reynolds equation, using finite difference method. The orbital response of crankshaft supported on the linear bearing model is obtained, considering balancer weights of motor rotor. And, the steady state displacement of crankshaft are compared with a variation in balancer weight. The loci of crankshaft at bearing locations are composed of the synchronous whirl component and the non-synchronous whirl component.

• Journal of Ocean Engineering and Technology
• /
• v.12 no.1
• /
• pp.10-22
• /
• 1998

In the present work the elastic-plastic FE formulations using dynamic explicit time integration schemes are used for numerical analysis of a large auto-body panel stamping processes. For analyses of more complex cases with larger and more refined meshes, the explicit method is more time effective than implicit method, and has no convergency problem and has the robust nature of contact and friction algorithms while implicit method is widely used because of excellent accuracy and reliability. The elastic-plastic scheme is more reliable and rigorous while the rigid-plastic scheme require small computation time. In finite element simulation of auto-body panel stamping processes, the roobustness and stability of computation are important requirements since the computation time and convergency become major points of consideration besides the solution accuracy due to the complexity of geometry conditions. The performnce of the dynamic explicit algorithms are investigated by comparing the simulation results of formaing of complicate shaped autobody parts, such as a fuel tank and a rear hinge, with the experimental results. It has been shown that the proposed dynamic explicit elastic-plastic finite element method enables an effective computation for complicated auto-body panel stamping processes.

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

An implicit viscous turbulent flow solver is developed for two-dimensional geometries on unstructured triangular meshes. The flux terms are discretized based on a cell-centered finite-volume formulation with the Roe's flux-difference splitting. The solution is advanced in time using an implicit backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with the Gauss-Seidel relaxation scheme. The effect of turbulence effects is approximated with a standard $k-{\varepsilon}$ two-equation model which is solved separately from the mean flow equations using the same backward-Euler time integration scheme. The triangular meshes are generated using an advancing-front/layer technique. Validations are made for flows over the NACA0012 airfoil and the Douglas 3-element airfoil. Good agreements are obtained between the numerical results and the experiment.

• International Journal of Aeronautical and Space Sciences
• /
• v.9 no.1
• /
• pp.85-99
• /
• 2008

The dynamic response analysis for large structures using finite element method requires a large amount of computational resources. This paper presents an efficient vibration analysis procedure by combining node-based substructuring reduction method with a response analysis scheme for structures with undamped, proportional or nonproportional damping. The iterative form of substructuring reduction scheme is derived to reduce the full eigenproblem and to calculate the dynamic responses. In calculating the time response, direct integration scheme is used because it can be applied directly to the reduced model. Especially for the non proportional damping matrix, the transformation matrices defined in the displacement space are used to reduce the system. The efficiency and the effectiveness of the present method are demonstrated through the numerical examples.

• Structural Engineering and Mechanics
• /
• v.42 no.2
• /
• pp.211-228
• /
• 2012

Present study deals with the development of finite element based solution methodology to investigate active control of dynamic response of delaminated composite shells with piezoelectric sensors and actuators. The formulation is based on first order shear deformation theory and an eight-noded isoparametric element is used. A coupled piezoelectric-mechanical formulation is used in the development of the constitutive equations. For modeling the delamination, multipoint constraint algorithm is incorporated in the finite element code. A simple negative feedback control algorithm coupling the direct and converse piezoelectric effects is used to actively control the dynamic response of delaminated composite shells in a closed loop employing Newmark's time integration scheme. The validity of the numerical model is demonstrated by comparing the present results with those available in the literature. A number of parametric studies such as the locations of sensor/actuator patches, delamination size and its location, radius of curvature to width ratio, shell types and loading conditions are carried out to understand their effect on the transient response of piezoceramic delaminated composite shells.

• Structural Engineering and Mechanics
• /
• v.20 no.2
• /
• pp.241-257
• /
• 2005

A numerical approach combining the finite element method with two different stability criteria namely the Budiansky and the phase-plane buckling criteria is used to study the dynamic buckling phenomena of plate and shell structures subjected to sudden applied loading. In the finite element analysis an explicit time integration scheme is used and the two criteria are implemented in the Finite Element analysis. The dynamic responses of the plate and shell structures have been investigated for different values of the plate and shell imperfection factors. The results indicate that the dynamic buckling time, which is normally considered in predicting elasto-plastic buckling behavior, should be taken into consideration with the buckling criteria for elastic buckling analysis of plate and shell structures. By selecting proper control variables and incorporating them with two dynamic buckling criteria, the unique dynamic buckling load can be obtained and the problems of ambiguity and contradiction of dynamic buckling load of plate and shell structure can be resolved.