• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.28 no.6A
• /
• pp.789-798
• /
• 2008

The stiffness method provides a framework to calculate the structural deformations directly from solving the equilibrium state. However, to use the displacement shape functions leads to approximate estimation of stiffness matrix and resisting forces, and accordingly results in a low accuracy. The conventional flexibility method uses the relation between sectional forces and nodal forces in which the equilibrium is always satisfied over all sections along the element. However, the determination of the element resisting forces is not so straightforward. In this study, a new fiber finite element mixed method has been developed for nonlinear anaysis of steel-concrete composite structures in the context of a standard finite element analysis program. The proposed method applies the Newton method based on the load control and uses the incremental secant stiffness method which is computationally efficient and stable. Also, the method is employed to analyze the steel-concrete composite structures, and the analysis results are compared with those obtained by ABAQUS. The comparison shows that the proposed method consistently well predicts the nonlinear behavior of the composite structures, and gives good efficiency.

• Structural Engineering and Mechanics
• /
• v.68 no.4
• /
• pp.385-398
• /
• 2018

In the present paper, we offer a new flat shell finite element. It is the result of the combination of a membrane element and a bending element, both based on the strain-based formulation. It is known that $C^{\circ}$ plane membrane elements provide poor deflection and stress for problems where bending is dominant. In addition, they encounter continuity and compliance problems when they connect to C1 class plate elements. The reach of the present work is to surmount these problems when a membrane element is coupled with a thin plate element in order to construct a shell element. The membrane element used is a triangular element with four nodes, three nodes at the vertices of the triangle and the fourth one at its barycenter. Each node has three degrees of freedom, two translations and one rotation around the normal. The coefficients related to the degrees of freedom at the internal node are subsequently removed from the element stiffness matrix by using the static condensation technique. The interpolation functions of strain, displacements and stresses fields are developed from equilibrium conditions. The plate element used for the construction of the present shell element is a triangular four-node thin plate element based on Kirchhoff plate theory, the strain approach, the four fictitious node, the static condensation and the analytic integration. The shell element result of this combination is robust, competitive and efficient.

• Transactions of the Korean Society of Mechanical Engineers
• /
• v.15 no.1
• /
• pp.1-10
• /
• 1991

Two-dimensional explicit finite element code was developed. The transient dynamics code can analyse large deformations of non-linear materials subjected to extremely high strain rates. The Lagrangian finite element program uses an explicit time integration operator to integrate the equations of motion, thus the stiffness matrix is not introduced. Cylinder upsetting and ring compression problems are simulated to check the effects of friction and inertia force. It is shown that (1) calculated results agree very well with experimental results, (2) constant shear friction method overestimates the decrease of inner ring radius and then underestimates after on in comparison with the Coulomb friction method, and (3) the effect of the increase in initial strain rate is similar to the effect of higher frictional coefficient.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.7 no.5
• /
• pp.27-35
• /
• 2007

This study presents a non -prismatic beam element for modeling the elastic and inelastic behavior of the steel beam, which has the post-Northridge connections in steel moment frames that are subjected to earthquake ground motions. The elastic stiffness matrix for non-prismatic members with reduced beam section (RES) connection is in the closed-form. The plasticity model is of a discrete type and is composed of a series of nonlinear hinges connected by rigid links. The hardening rules can model the inelastic behavior for monotonic and random cyclic loading, and the effects of local buckling. Verification and calibration of the model are presented in a companion paper.

• Structural Engineering and Mechanics
• /
• v.26 no.6
• /
• pp.725-739
• /
• 2007

In recent years there are many plate bending elements that emerged for solving both thin and thick plates. The main features of these elements are that they are based on mix formulation interpolation with discrete collocation constraints. These elements passed the patch test for mix formulation and performed well for linear analysis of thin and thick plates. In this paper a member of this family of elements, namely, the Discrete Reissner-Mindlin (DRM) is further extended and developed to analyze both thin and thick plates with geometric nonlinearity. The Von K$\acute{a}$rm$\acute{a}$n's large displacement plate theory based on Lagrangian coordinate system is used. The Hu-Washizu variational principle is employed to formulate the stiffness matrix of the geometrically Nonlinear Discrete Reissner-Mindlin (NDRM). An iterative-incremental procedure is implemented to solve the nonlinear equations. The element is then tested for plates with simply supported and clamped edges under uniformly distributed transverse loads. The results obtained using the geometrically NDRM element is then compared with the results of available analytical solutions. It has been observed that the NDRM results agreed well with the analytical solutions results. Therefore, it is concluded that the NDRM element is both reliable and efficient in analyzing thin and thick plates with geometric non-linearity.

• Advances in Automotive Engineering
• /
• v.1 no.1
• /
• pp.41-59
• /
• 2018

In the study, investigation of fiber- reinforced composite materials that can be an alternative to conventional steel was performed by finite element analysis with the help of software. Steel and composite materials have been studied on a four axle truck chassis model. Three-dimensional finite element model was created with software, and then analyzes were performed. The analyses were performed for static and dynamic/fatigue cases. Fatigue cases are formed with the help of design spectra model and fatigue analyses were performed as static analyses with this design spectra. First, analyses were performed for steel and after that optimization analyses were made for the AS4-PEEK carbon fiber composite and Eglass-Epoxy fiber composite materials. Optimization of composite material analyzes include determining the total laminate thickness, thickness of each ply, orientation of each ply and ply stacking sequence. Analyzes were made according to macro mechanical properties of composite, micromechanics case has not been considered. Improvements in weight reduction up to %50 provided at the end of the composite optimization analyzes with satisfying stiffness performance of chassis. Fatigue strength of the composite structure depends on various factors such as, fiber orientation, ply thickness, ply stack sequence, fiber ductility, ductility of the matrix, loading angle. Therefore, the accuracy of theoretical calculations and analyzes should be correlated by testing.

• Structural Engineering and Mechanics
• /
• v.23 no.1
• /
• pp.15-27
• /
• 2006

A locally varying temperature field or a mixture of two or more different materials can cause local variation of elasticity properties of a beam. In this paper, a new Euler-Bernoulli beam element with varying Young's modulus along its longitudinal axis is presented. The influence of axial forces according to the linearized 2nd order beam theory is considered, as well. The stiffness matrix of this element contains the transfer constants which depend on Young's modulus variation and on axial forces. Occurrence of the polynomial variation of Young's modulus has been assumed. Such approach can be also used for smooth local variation of Young's modulus. The critical loads of the straight slender columns were studied using the new beam element. The influence of position of the local Young's modulus variation and its type (such as linear, quadratic, etc.) on the critical load value and rate of convergence was investigated. The obtained results based on the new beam element were compared with ANSYS solutions, where the number of elements gradually increased. Our results show significant influence of the locally varying Young's modulus on the critical load value and the convergence rate.

• Journal of the Korean Society for Precision Engineering
• /
• v.12 no.3
• /
• pp.101-109
• /
• 1995

Axisymmetric and quasi-static finite element analysis of an inflated tire rotating with constant angular velocity and contact to road has been performed. Centrifugal force effect was added to load stiffness matrix and equation of effective material properties were calculated by the Halpin-Tsai formulation. In this report, radial truck/bus tire was analyzed. It was inflated and rotated at speeds up to 140 km/h. Then, contact problem was performed to calculate stress-strain field of tire wiht flat rigid road under the load due to the self-weight of a vehicle. Significant changes of stress-strain field of tire were observed in the finite element analysis. Shear stress, strain and strain energy density were rapidly increased at the dege of #2 belt at freely rotating state. This concentrated stress and strain made belt edge sparation. Under the condition of flat riged road contact, strain energy density of #2 belt, carcass turn-up part were concentrated and bigger values than only freely rotation state. Therefore, dynamic behaivor of tire has to considered as design factors which are affected to belt edge separation and bead breakage.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.19 no.3
• /
• pp.36-41
• /
• 2020

Mechanical systems installed in transport devices, such as vehicles, airplanes, and ships, are mostly subject to translational accelerations at the joints during operations. This base acceleration excitation has a large influence on the performance of the system, therefore, its response must be well analyzed. However, the existing methods for dynamic analysis of structures have some limitations in use. This study presents a new numerical method using relative acceleration to solve these limitations. If the governing equation of motion is linear and the mass matrix, the damping matrix, and the stiffness matrix are constant over time in the finite element analysis, the proposed method can be applied to the transient behavior analysis and the harmonic response analysis of the structure. Because it is not necessary to introduce a virtual mass and the rigid body motions are removed from the analysis, it is possible to use not only the direct integration method in the time domain but also the mode superposition method to obtain the dynamic responses. This paper demonstrates with three examples how the present method is suitable for the dynamic analysis of a structure with multi-degree of freedom.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.1
• /
• pp.11-20
• /
• 2004

The main goal of this study is to develop MATLAB programming for an analysis of distortional deformations and stresses of the straight box girder. For this purpose, a distortional deformation theory is firstly summarized and then a BEF (Beam on Elastic Foundation) theory is presented using analogy of the corresponding variables. Finally, with governing equations of the beam-column element on elastic foundation, an exact element stiffness matrix of the beam element and nodal forces equivalent to concentrated and distributed loads are evaluated via a generalized linear eigenvalue problem. In order to verify the efficiency and accuracy of this method, distortional stresses of box girders with multiple diaphragms are presented and compared with results by FEA.