• Transactions of the Korean Society of Automotive Engineers
• /
• v.6 no.3
• /
• pp.54-62
• /
• 1998

In this study, collapse behavior of frame composed of thin-walled rectangular tube is investigated. Considering the collapse of frame, the bending and compression members undergo large deformation. The stiffness of the compound element is obtained from analytical moment-rotation relationship and approximated load-deflection relationsh- ip of thin-walled rectangular tube. A computer program is developed for the large deformation analysis of frame. An incremental displacement method is used in the program and at each incremental stage, the stiffness matrix of the total structure is checked with the state of each element for bending and compression.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1994.04a
• /
• pp.57-64
• /
• 1994

The purpose of this paper is to find minimum surface shape of pneumatic structure using the finite element method. The pneumatic membrane structure is a kind of large deformation problem and very flexible composite material, which mean geomatric nonlinearity. It is not to resist for compression and resultant moment. As the displacement due to internal pressure is getting bigger, it should be considered the direction of forces. It becomes non-linear problem with the non-conservative force. The follower-force depends on the deformation and the direction of force is normal to each element. The solution process is obtained the new stiffness matrix (load correction matrix) depending on deformation through each iterated step. However, the stiffness matrix have not the symmetry and influence on the time of covergence. So in this paper Newton-Rhapson method for solving non-linear problem and for using symmetic matrix, the load direction is changed in each iterated step using the transformation matrix.

• Journal of the Korean Society of Manufacturing Technology Engineers
• /
• v.6 no.3
• /
• pp.58-64
• /
• 1997

A micromechanics model to describe the elastic behavior of fiber or whisker reinforced metal matrix composites was developed and the stress concentrations between reinforcements were investigated using the modified shear lag model with the comparison of finite element analysis (FEA). The rationale is based on the replacement of the matrix between fiber ends with the fictitious fiber to maintain the compatibility of displacement and traction. It was found that the new model gives a good agreement with FEA results in the small fiber aspect ratio regime as well as that in the large fiber aspect ratio regime. By the calculation of the present model, stress concentration factor in the matrix and the composite elastic modulus were predicted accurately. Some important factors affecting stress concentrations, such as fiber volume fraction, fiber aspect ratio, end gap size, and modulus ratio, were also discussed.

• Proceedings of the Korean Society of Machine Tool Engineers Conference
• /
• 1998.10a
• /
• pp.261-266
• /
• 1998

In this study, crush characteristics of thin-walled rectangular tube is investigated. The stiffness of the element is obtained from analytical moment-rotation relationship and approximated load-deflection relationship of thin-walled rectangular tube. A computer program is developed for the large deformation analysis of frame. An incremental displacement method is used in the program and at each incremental stage, the stiffness matrix of the total structure is checked with the state each element for bending and compression.

• 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.

• 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.

• Structural Engineering and Mechanics
• /
• v.37 no.6
• /
• pp.633-648
• /
• 2011

This paper investigates the aspect-ratio locking of the isoparametric 8-node quadrilateral (QUAD8) element. An important finding is that, if finite element solution is carried out with in exact arithmetic (i.e., with no truncation and round off errors), the locking tendency of the element is completely avoided even for aspect-ratios as high as 100000. The current finite element codes mostly use floating point arithmetic. Thus, they can only avoid this locking for aspect-ratios up to 100 or 1000. A novel method is proposed in the paper to avoid aspect-ratio locking in floating point computations. In this method, the offending terms of the strain-displacement matrix (i.e., $\mathbf{B}$-matrix) are multiplied by suitable scaling factors to avoid ill-conditioning of stiffness matrix. Numerical examples are presented to demonstrate the efficacy of the method. The examples reveal that aspect-ratio locking is avoided even for aspect-ratios as high as 100000.

• Journal of the Korea Institute of Military Science and Technology
• /
• v.12 no.4
• /
• pp.540-547
• /
• 2009

In the near future, military robots are likely to be substituted for military personnel in the field of battle. The power system of a legged robot is considerably more complex than the one used for a land vehicle because of the coordination and stability issues due to the large number of degree of freedom. In this paper, a servovalve-piston combination system for a straight-line motion of robot leg is modeled as three degree of freedom based on double inputs and single output transfer function. The output is the displacement of piston from neutral. The inputs are valve displacement from neutral and arbitrary load force in this system. LQR(Linear Quadratic Regulator) technique is applied in order to achieve robust stability and fast responses of the system. The Kalman filter loop, rejection of disturbance and noise, riccati equation, filter gain matrix, and frequency domain equality are analyzed and designed.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.05a
• /
• pp.117-119
• /
• 2005

For large structural eigen-analysis, the whole structure is divided into some partitioned structures and through synthesis of partitioned structural model the eigen-data of structure can be obtained. In that case, eigenvalue problem consists of semidefinite mass matrix form because of displacement constraint condition. In this work the eigenvalue problem is considered by means of several method, determinant search and null space reduction method.

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

Stepped bean with immovable ends for large amplitude of vibration including effects of longitudinal displacement, shear deformation and rotary inertia is investigated for free and forced vibration using finite element method. Modified harmonic force matrix is introduced for analysis of vibration with finite amplitude of the stepped beam under uniform hamonic loading and beam with nonuniform harmonic loading. Numerical examples of stepped beam with various support conditions are analysed for deflections and natural frequencies. Results show that the proposed method is valid and efficient.