• Structural Engineering and Mechanics
• /
• v.51 no.4
• /
• pp.601-625
• /
• 2014

Based on continuum mechanics and the principle of virtual displacements, incremental total Lagrangian formulation (T.L.) and incremental updated Lagrangian formulation (U.L.) were presented. Both T.L. and U.L. considered the large displacement stiffness matrix, which was modified to be symmetrical matrix. According to the incremental updated Lagrangian formulation, small strain, large displacement, finite rotation of three dimensional Timoshenko fiber beam element tangent stiffness matrix was developed. Considering large displacement and finite rotation, a new type of tangent stiffness matrix of the beam element was developed. According to the basic assumption of plane section, the displacement field of an arbitrary fiber was presented in terms of nodal displacement of centroid of cross-area. In addition, shear deformation effect was taken account. Furthermore, a nonlinear finite element method program has been developed and several examples were tested to demonstrate the accuracy and generality of the three dimensional beam element.

• Smart Structures and Systems
• /
• v.15 no.3
• /
• pp.645-663
• /
• 2015

Recently, an indirect displacement estimation method using data fusion of acceleration and strain (i.e., acceleration-strain-based method) has been developed. Though the method showed good performance on beam-like structures, it has inherent limitation in applying to more general types of bridges that may have complex shapes, because it uses assumed analytical (sinusoidal) mode shapes to map the measured strain into displacement. This paper proposes an improved displacement estimation method that can be applied to more general types of bridges by building the mapping using the finite element model of the structure rather than using the assumed sinusoidal mode shapes. The performance of the proposed method is evaluated by numerical simulations on a deck arch bridge model and a three-span truss bridge model whose mode shapes are difficult to express as analytical functions. The displacements are estimated by acceleration-based method, strain-based method, acceleration-strain-based method, and the improved method. Then the results are compared with the exact displacement. An experimental validation is also carried out on a prestressed concrete girder bridge. The proposed method is found to provide the best estimate for dynamic displacements in the comparison, showing good agreement with the measurements as well.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.22 no.2
• /
• pp.278-288
• /
• 1998

An efficient and useful method to improve the local stress field in a displacement-formulated finite element solution has been proposed using the theory of conjugate approximations for a stress field and the Loubignac's iterative method for a displacement field. Validity of the proposed method has been tested through three test examples, to improve the stress field and displacement field in the whole domain and the local regions. As a result of analysis on the test examples, it is found that the stress field in the local regions are approximated to those in the whole domain within a few iterations which have satisfied the original finite element equilibrium equation. In addition, it is found that the local stress field are by far better approximated to the exact stress field than the displacement-based stress field with the reduction of the finite-element mesh-size.

• Structural Engineering and Mechanics
• /
• v.26 no.1
• /
• pp.43-68
• /
• 2007

Using the Mindlin-Reissner plate theory, many quadrilateral plate bending elements have been developed so far to analyze thin and moderately thick plate problems via displacement based finite element method. Here new formulation has been made to analyze thin and moderately thick plate problems using force based finite element method called Integrated Force Method (IFM). The IFM is a novel matrix formulation developed in recent years for analyzing civil, mechanical and aerospace engineering structures. In this method all independent/internal forces are treated as unknown variables which are calculated by simultaneously imposing equations of equilibrium and compatibility conditions. In this paper the force based new bilinear quadrilateral plate bending element (MQP4) is proposed to analyze the thin and moderately thick plate bending problems using Integrated Force Method. The Mindlin-Reissner plate theory has been used in the formulation of this element which accounts the effect of shear deformation. Standard plate bending benchmark problems are analyzed using the proposed element MQP4 via Integrated Force Method to study its performance with respect to accuracy and convergence, and results are compared with those of displacement based 4-node quadrilateral plate bending finite elements available in the literature. The results are also compared with the exact solutions. The proposed element MQP4 is free from shear locking and works satisfactorily in both thin and moderately thick plate bending situations.

• Journal of the Korean Society of Hazard Mitigation
• /
• v.9 no.1
• /
• pp.9-13
• /
• 2009

A simple outline for teaching adaptive scheme based finite element method for planar problems as a part of computer aided teaching of structural engineering curriculum is presented. Displacement based finite element formulation for planar problems and representative strain value based adaptive scheme for mesh generation are considered. As examples, a cantilever beam with a concentrated load treated as a planar problem and stretching of a plate with a circular hole are analyzed with displacement based finite element method with adaptive meshes. The examples and outlines show how adaptive based finite element method may become an essential part of computer aided teaching of structural engineering.

• Proceedings of the Korean Society For Composite Materials Conference
• /
• 2003.04a
• /
• pp.143-146
• /
• 2003

A finite element based on the efficient higher order zig-zag theory with multiple delaminations Is developed to refine the predictions of frequency and mode shapes. Displacement field through the thickness are constructed by superimposing linear zig-zag field to the smooth globally cubic varying field. The layer-dependent degrees of freedom of displacement fields are expressed in terms of reference primary degrees of freedom by applying interface continuity conditions including delaminated interfaces as well as free hounding surface conditions of transverse shear stresses. Thus the proposed theory is not only accurate but also efficient. This displacement field can systematically handle the number, shape, size, and locations of delaminations. Throught the dynamic version of variational approach, the dynamic equilibrium equations and variationally consistent boundary conditions are obtained. Through the natural frequency analysis and time response analysis of composite plate with multiple delaminations, the accuracy and efficiency of the present finite element are demonstrated. The present finite element is suitable in the predictions of the dynamic response of the thick composite plate with multiple delaminations.

• Coupled systems mechanics
• /
• v.8 no.5
• /
• pp.393-414
• /
• 2019

This paper presents seismic analysis of concrete gravity dams considering soil-structure-fluid interaction. Displacement based plane strain finite element formulation is considered for the dam and foundation domain whereas pressure based finite element formulation is considered for the reservoir domain. A direct coupling method has been adopted to obtain the interaction effects among the dam, foundation and reservoir domain to obtain the dynamic responses of the dam. An efficient absorbing boundary condition has been implemented at the truncation surfaces of the foundation and reservoir domains. A parametric study has been carried out considering each domain separately and collectively based on natural frequencies, crest displacement and stress at the neck level of the dam body. The combined frequency of the entire coupled system is very less than that of the each individual sub-system. The crest displacement and neck level stresses of the dam shows prominent enhancement when coupling effect is taken into consideration. These outcomes suggest that a complete coupled analysis is necessary to obtain the actual responses of the concrete gravity dam. The developed methodology can easily be implemented in finite element code for analyzing the coupled problem to obtain the desired responses of the individual subdomains.

• Structural Engineering and Mechanics
• /
• v.2 no.4
• /
• pp.395-402
• /
• 1994

An improved finite element-transfer matrix method is applied to the transient analysis of plates with large displacement under various excitations. In the present method, the transfer of state vectors from left to right in a combined finite element-transfer matrix method is changed into the transfer of generally incremental stiffness equations of every section from left to right. Furthermore, in this method, the propagation of round-off errors occurring in recursive multiplications of transfer and point matrices is avoided. The Newmark-${\beta}$ method is employed for time integration and the modified Newton-Raphson method for equilibrium iteration in each time step. An ITNONDL-W program based on this method using the IBM-PC/AT microcomputer is developed. Finally numerical examples are presented to demonstrate the accuracy as well as the potential of the proposed method for dynamic large deflection analysis of plates with random boundaries under various excitations.

• Structural Engineering and Mechanics
• /
• v.7 no.2
• /
• pp.203-216
• /
• 1999

A previously proposed finite element formulation method is refined and modified to generate a new type of elements. The method is based on selecting a set of general solution modes for element formulation. The constant strain modes and higher order modes are selected and the formulation method is designed to ensure that the element will pass the basic single element test, which in turn ensures the passage of the basic patch test. If the element is to pass the higher order patch test also, the element stiffness matrix is in general asymmetric. The element stiffness matrix depends only on a nodal displacement matrix and a nodal force matrix. A symmetric stiffness matrix can be obtained by either modifying the nodal displacement matrix or the nodal force matrix. It is shown that both modifications lead to the same new element, which is demonstrated through numerical examples to be more robust than an assumed stress hybrid element in plane stress application. The method of formulation can also be used to arrive at the conforming displacement and hybrid stress formulations. The convergence of the latter two is explained from the point of view of the proposed method.

• Advanced Composite Materials
• /
• v.18 no.1
• /
• pp.21-42
• /
• 2009

The vibration characteristics of post-buckled laminated composite doubly curved shells are investigated. The finite element method is used for the analysis of post-buckling and free vibration of post-buckled laminated shells. The geometric non-linear finite element model includes the general non-linear terms in the strain-displacement relationships. The shell geometry used in the present formulation is derived using an orthogonal curvilinear coordinate system. Based on the principle of virtual work the non-linear finite element equations are derived. Arc-length method is implemented to capture the load-displacement equilibrium curve. The vibration characteristics of post-buckled shell are performed using tangent stiffness obtained from the converged deflection. The code is first validated and then employed to generate numerical results. Parametric studies are performed to analyze the snapping and vibration characteristics. The relationship between loads and fundamental frequencies and between loads and the corresponding displacements are determined for various parameters such as thickness ratio and shallowness.