• Journal of Mechanical Science and Technology
• /
• v.16 no.10
• /
• pp.1201-1210
• /
• 2002

A penalty frame method is proposed for the coupled analysis of finite elements with independently modeled substructures. Although previously reported hybrid interface method by Aminpour et al (IJNME, Vol 38, 1995) is accurate and reliable, it requires non-conventional special solution algorithm such as multifrontal solver. In present study, an alternative method has been developed using penalty frame constraints, which results in positive symmetric global stiffness matrices. Thus the conventional skyline solver or band solver can be utilized in the solution routine, which makes the present method applicable in the environment of conventional finite element commercial software. Numerical examples show applicability of the present method.

• Structural Engineering and Mechanics
• /
• v.4 no.4
• /
• pp.437-450
• /
• 1996

A study on the behaviour of fibre reinforced concrete deep beams with and without web openings is carried out using nonlinear finite element analysis. Eight node isoparametric plane stress elements are employed to model the fibre reinforced concrete materials. Steel bars are treated using a compatible three node truss elements. The constitutive equations for fibre reinforced concrete materials take into account the softening effect of co-existing shear strains. Element stiffness at each step is formulated based on the tangent modulus at the current level of principal strains. Transformation between principal directions and global coordinate system is imposed. Comparison of analytical results with experimental values indicates reasonably good agreement. The proposed numerical model can be used to study the behaviour of this composite structures of practically any geometries.

• International Journal of CAD/CAM
• /
• v.6 no.1
• /
• pp.73-79
• /
• 2006

In this paper, we propose a new triangular finite element mesh generation method based on simplification of high-density mesh and adaptive Level-of-Detail (LOD) methods for efficient CAE. In our method, mesh simplification is used to control the mesh properties required for FE mesh, such as the number of triangular elements, element shape quality and size while keeping the specified approximation tolerance. Adaptive LOD methods based on vertex hierarchy according to curvature and region of interest, and global LOD method preserving density distributions are also proposed in order to construct a mesh more appropriate for CAE purpose. These methods enable efficient generation of FE meshes with properties appropriate for analysis purpose from a high-density mesh. Finally, the effectiveness of our approach is shown through evaluations of the FE meshes for practical use.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.5 no.6
• /
• pp.44-52
• /
• 1997

As part of the effort to reduce the weight of powertrain, a hollow crankshaft has been designed. The mass reduction of the crankshaft changes the dynamic properties of the crankshaft such as moment of inertia, and torsional, bending stiffness. The purpose of this paper is to compare the dynamic behavior of the hollow crankshaft with that of the original, solid crankshaft. Global dynamic behavior of the crankshaft is analyzed bgy the transfer matrix method(TMM). The crankshaft has been modeled by 38 lumped mass and stiffness elements. The dynamic patameters of each lumped element are provided by Finite Element Method(FEM). The responses of the crankshaft from TMM are fed back as loading conditions to the Finite Element model to obtain dynamic stresses for critical areas of the crankshaft.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.859-865
• /
• 2007

This paper presents a method to analyze the unbalance response of a high speed polygon mirror scanner motor supported by sintered bearing and flexible supporting structures by using the finite element method and the mode superposition method. The appropriate finite element equations for polygon mirror are described by rotating annular sector element using Kirchhoff plate theory and von Karman non-linear strain, and its rigid body motion is also considered. The rotating components except for the polygon mirror are modeled by Timoshenko beam element including the gyroscopic effect. The flexible supporting structures are modeled by using a 4-node tetrahedron element and 4-node shell element with rotational degrees of freedom. Finite element equations of each component of the polygon mirror scanner motor and the flexible supporting structures are consistently derived by satisfying the geometric compatibility in the internal boundary between each component. The rigid link constraints are also imposed at the interface area between sleeve and sintered bearing to describe the physical motion at this interface. 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 solving the associated eigenvalue problem by using the restarted Arnoldi iteration method. Unbalance responses in time and frequency domain are performed by superposing the eigenvalues and eigenvectors from the free vibration analysis. The validity of the proposed method is verified by comparing the simulated unbalance response with the experimental results. This research also shows that the flexibility of supporting structures plays an important role in determining the unbalance response of the polygon mirror scanner motor.

• Computers and Concrete
• /
• v.13 no.2
• /
• pp.209-220
• /
• 2014

Vibration based damage detection is very popular in the civil engineering area. Especially, special structures like dams, long-span bridges and high-rise buildings, need continues monitoring in terms of mechanical properties of material, static and dynamic behavior. It has been stated in the International Commission on Large Dams that more than half of the large concrete dams were constructed more than 50 years ago and the old dams have subjected to repeating loads such as earthquake, overflow, blast, etc.,. So, some unexpected failures may occur and catastrophic damages may be taken place because of theloss of strength, stiffness and other physical properties of concrete. Therefore, these dams need repairs provided with global damage evaluation in order to preserve structural integrity. The paper aims to show the effectiveness of the model updating method for global damage detection on a laboratory arch dam model. Ambient vibration test is used in order to determine the experimental dynamic characteristics. The initial finite element model is updated according to the experimentally determined natural frequencies and mode shapes. The web thickness is selected as updating parameter in the damage evaluation. It is observed from the study that the damage case is revealed with high accuracy and a good match is attained between the estimated and the real damage cases by model updating method.

• Magazine of the Korean Society of Agricultural Engineers
• /
• v.27 no.2
• /
• pp.38-46
• /
• 1985

The reinforced concrete farm silos on the elastic foundatin are widely used in agricultural engineering because of their superior structural performance, economy and attractive appearance. Various methods for the analysis and design of farm silo, such as the analytical method, the finite difference method, and the finite element methods, can be used. But the analytical procedure can not be applied for the intricate conditions in practice. Therefore lately the finite element method has been become in the structural mechanics. In this paper, a method of finite element analysis for the cylindrical farm silo on ffness matrix for the elastic foundation governed by winkler's assumption. A complete computer programs have been developed in this paper can be applicable not only to the shell structures on elastic foundation but also to the arbitrary three dimensional structures. Assuming the small deflection theory, the membrane and plate bending behaviours of flat plate element can be assumed mutually uncoupled. In this case, the element has 5 degrees of freedom per node when defined in the local coordinate system. However, when the element properties are transformed to the global coordinates for assembly, the 6th degree of freedom should be considered. A problem arises in this procedure the resultant stiffness in the 6th degree of freedom at this node will be zero. But this singularity of the stiffness matrix can be eliminated easily by merely replacing the zero diagonal by dummy stiffness.

• Structural Engineering and Mechanics
• /
• v.83 no.3
• /
• pp.353-361
• /
• 2022

Linear Elastic Fracture Mechanics (LEFM) has been developed by applying stress analysis to determine the stress intensity factor (SIF, K). The finite element method (FEM) is widely used as a standard tool for evaluating the SIF for various crack configurations. The prediction accuracy can be achieved by applying an adaptive Delaunay triangulation combined with a FEM. The solution can be solved using either direct or iterative solvers. This work adopts the element-by-element preconditioned conjugate gradient (EBE-PCG) iterative solver into an adaptive FEM to solve the solution to heal problem size constraints that exist when direct solution techniques are applied. It can avoid the formation of a global stiffness matrix of a finite element model. Several numerical experiments reveal that the present method is simple, fast, and efficient compared to conventional sparse direct solvers. The optimum convergence criterion for two-dimensional LEFM analysis is studied. In this paper, four sample problems of a two-edge cracked plate, a center cracked plate, a single-edge cracked plate, and a compact tension specimen is used to evaluate the accuracy of the prediction of the SIF values. Finally, the efficiency of the present iterative solver is summarized by comparing the computational time for all cases.

• Proceedings of the Korean Society of Precision Engineering Conference
• /
• 1996.04a
• /
• pp.842-846
• /
• 1996

Leaf springs are widely used as a major suspension component in many commercial vehicles, such as buses, trucks, etc. They have a complex dynamic behavior due to the geometric nonlinear and the contact mechanism between the leaves. The interface conditions between the leaves play a significant role in the global behavior of the comfort and ride of the vehicle system. The paper concentrates on modeling leaf springs and contact frictions between the leaves using a nonlinear finite element approach. A nonlinear load-displacement hysteresis curve for the leaf spring is simulated and its results are compared with test results.

• Structural Engineering and Mechanics
• /
• v.24 no.6
• /
• pp.647-664
• /
• 2006

In this study, an improved finite element formulation with a scheme of solution for the large deflection analysis of inextensible prismatic and nonprismatic slender beams is developed. For this purpose, a three-noded Lagrangian beam-element with two dependent degrees of freedom per node (i.e., the vertical displacement, y, and the actual slope, $dy/ds=sin{\theta}$, where s is the curved coordinate along the deflected beam) is used to derive the element stiffness matrix. The element stiffness matrix in the global xy-coordinate system is achieved by means of coordinate transformation of a highly nonlinear ($6{\times}6$) element matrix in the local sy-coordinate. Because of bending with large curvature, highly nonlinear expressions are developed within the global stiffness matrix. To achieve the solution after specifying the proper loading and boundary conditions, an iterative quasi-linearization technique with successive corrections are employed considering these nonlinear expressions to remain constant during all iterations of the solution. In order to verify the validity and the accuracy of this study, the vertical and the horizontal displacements of prismatic and nonprismatic beams subjected to various cases of loading and boundary conditions are evaluated and compared with analytic solutions and numerical results by available references and the results by ADINA, and excellent agreements were achieved. The main advantage of the present technique is that the solution is directly obtained, i.e., non-incremental approach, using few iterations (3 to 6 iterations) and without the need to split the stiffness matrix into elastic and geometric matrices.