• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.4
• /
• pp.737-752
• /
• 1992

Analysis of elasto-dynamic deformation of flexible linkage mechanism is conducted using the finite element method. The equations of motion of the system are derived from the static structural problem in which dynamic inertia, gravitational and driving forces are treated as external loads. Linear spring model is included in the formulation of equation of motions to represent the effects of deformation of elastic bearings of revolute joints on the system behavior. A computer program is constructed and applied to analyze a specific crank-lever 4-bar mechanism. The algorithm of the program is as follows. First, the natural frequencies and the mode shapes of the system are calculated by solving the eigenproblem of the mechanism system which can be considered as a static structure by assuming the input shaft (crank shaft) to be fixed at any given configuration of mechanism. And finally, the elasto-dynamic deformation of the whole system is obtained using mode superposition method for the case of constant input speed. The effect of geometric stiffness on the mechamism is included in the program with the axial forces of links obtained through the quasi-static displacement analysis. It is found that the geometric stiffness exerts an important effect upon the elasto-dynamic behavior of the flexible linkage mechanism. Elastic deformation of bearing lowers the natural frequencies of the system, resulting smaller elastic displacement at the mid-point of the links and bigger elestic displacement at the ends of the links than rigid bearing. The above investigation of flexible linkage mechanism shows that the effects of the elastic deformation of bearing on the mechanism should be considered to design the mechanism which satisfies more preciously the purpose and the condition of design.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.1
• /
• pp.1-8
• /
• 2014

Using the level set and the meshfree methods, we develop a topological shape optimization method applied to linear elasticity problems. Design gradients are computed using an efficient adjoint design sensitivity analysis(DSA) method. The boundaries are represented by an implicit moving boundary(IMB) embedded in the level set function obtainable from the "Hamilton-Jacobi type" equation with the "Up-wind scheme". Then, using the implicit function, explicit boundaries are generated to obtain the response and sensitivity of the structures. Global nodal shape function derived on a basis of the reproducing kernel(RK) method is employed to discretize the displacement field in the governing continuum equation. Thus, the material points can be located everywhere in the continuum domain, which enables to generate the explicit boundaries and leads to a precise design result. The developed method defines a Lagrangian functional for the constrained optimization. It minimizes the compliance, satisfying the constraint of allowable volume through the variations of boundary. During the optimization, the velocity to integrate the Hamilton-Jacobi equation is obtained from the optimality condition for the Lagrangian functional. Compared with the conventional shape optimization method, the developed one can easily represent the topological shape variations.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.26 no.5
• /
• pp.385-392
• /
• 2013

The finite element method has become the most widely used method of structural analysis and recently, the method has often been applied to complex dynamic and nonlinear structural analyses problems. Even for these complex problems, where the responses are hard to predict, finite element analyses yield reliable results if appropriate element types and meshes are used. However, the dynamic and nonlinear behaviors of a structure often include large deformations in various portions of the structure and if the same mesh is used throughout the analysis, some elements may deform to shapes beyond the reliable limits; thus dynamically adapting finite element meshes are needed in order for the finite element analyses to be accurate. In addition, to satisfy the users requirement of quick real run time of finite element programs, the algorithms must be computationally efficient. This paper presents an adaptive finite element mesh generation scheme for dynamic analyses of structures that may adapt at each time step. Representative strain values are used for error estimates and combinations of the h-method(node movement) and the r-method(element division) are used for mesh refinements. A coefficient that depends on the shape of an element is used to limit overly distorted elements. A simple frame example shows the accuracy and computational efficiency of the scheme. The aim of the study is to outline the adaptive scheme and to demonstrate the potential use in general finite element analyses of dynamic and nonlinear structural problems commonly encountered.

• Geophysics and Geophysical Exploration
• /
• v.15 no.1
• /
• pp.33-38
• /
• 2012

A resistivity method has been applied to wide range of engineering and environmental problems with the help of automatic and precise data acquisition. Thus, more accurate modeling and inversion of time-lapse monitoring data are required since resistivity monitoring has been introduced to quantitatively find out subsurface changes With respect to time. Here, we used the finite element method (FEM) for 3D resistivity modeling since the method is easy to realize complex topography and arbitrary shaped anomalous bodies. In the FEM, the linear elements, also referred to as first order elements, have certain advantages of simple formulation and narrow bandwidth of system equation. However, the linear elements show the poor accuracy and slow convergence of the solution with respect to the number of elements or nodes. To achieve the higher accuracy of finite element solution, high order elements are generally used. In this study, we developed a 3D resistivity modeling program using high order Serendipity elements. Comparing the Serendipity element solutions for a cube model with the linear element solutions, we assured that the Serendipity element solutions are more accurate than the linear element solutions in the 3D resistivity modeling.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.37 no.1
• /
• pp.49-56
• /
• 2024

An automated dynamic structural analysis module stands as a crucial element within a structural integrated mitigation system. This module must deliver prompt real-time responses to enable timely actions, such as evacuation or warnings, in response to the severity posed by the structural system. The finite element method, a widely adopted approximate structural analysis approach globally, owes its popularity in part to its user-friendly nature. However, the computational efficiency and accuracy of results depend on the user-provided finite element mesh, with the number of elements and their quality playing pivotal roles. This paper introduces a computationally efficient adaptive mesh generation scheme that optimally combines the h-method of node movement and the r-method of element division for mesh refinement. Adaptive mesh generation schemes automatically create finite element meshes, and in this case, representative strain values for a given mesh are employed for error estimates. When applied to dynamic problems analyzed in the time domain, meshes need to be modified at each time step, considering a few hundred or thousand steps. The algorithm's specifics are demonstrated through a standard cantilever beam example subjected to a concentrated load at the free end. Additionally, a portal frame example showcases the generation of various robust meshes. These examples illustrate the adaptive algorithm's capability to produce robust meshes, ensuring reasonable accuracy and efficient computing time. Moreover, the study highlights the potential for the scheme's effective application in complex structural dynamic problems, such as those subjected to seismic or erratic wind loads. It also emphasizes its suitability for general nonlinear analysis problems, establishing the versatility and reliability of the proposed adaptive mesh generation scheme.

• Journal of Korean Society of Steel Construction
• /
• v.21 no.6
• /
• pp.597-608
• /
• 2009

In this paper, large-deformation elasto-plastic analysis of space frames that considersjoint connection properties is presented. This method is based on the large-deformation formula with finite rotation, which was developed initially for elastic systems, and is extended herein to include the elasto-plastic effect and the member joint connection properties of semi-rigid what?. The analytical method was derived from the Eulerian concept, which takes into consideration the effects of large joint translations and rotations. The localmember force-deformation relationships were obtained from the beam-column approach, and the change caused by the axial strain in the member chord lengths and flexural bowing were taken into account. The effect of the axial force of the member on bending and torsional stiffness, and on the plastic moment capacity, is included in the analysis. The material is assumed to be ideally elasto-plastic, and yielding is considered concentrated at the member ends in the form of plastic hinges. The semi-rigid properties of the member joint connection are considered based on the power or linear model. The arc length method is usedto trace the post-buckling range of the elastic and elasto-plastic problems with the semi-rigid connection. A sample non-linear buckling analysis was carried out with the proposed space frame formulations to demonstrate the potential of the developed method in terms of its accuracy and efficiency.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.4
• /
• pp.375-387
• /
• 2004

A new finite element model will be presented to analyze the nonlinear behavior of RC beams and slabs strengthened by a patch repair. The numerical approach is based on the p-version degenerate shell element including theory of anisotropic laminated composites, theory of materially and geometrically nonlinear plates. In the nonlinear formulation of this model, the total Lagrangian formulation is adopted with large deflections and moderate rotations being accounted for in the sense of von Karman hypothesis. The material model is based on hardening rule, crushing condition, plate-end debonding strength model and so on. The Gauss-Lobatto numerical quadrature is applied to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed p-version nonlinear finite element model is demonstrated through the load-deflection curves, the ultimate loads, and the failure modes of RC beams or slabs bonded with steel plates or FRP plates compared with available result of experiment and other numerical methods.

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

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.1
• /
• pp.115-130
• /
• 2013

The sectional methods of current design codes have been broadly used for the design of various kinds of reinforced concrete pile caps. Lately, the strut-tie model approach of current design codes also became one of the attracting methods for pile caps. However, since the sectional methods and the strut-tie model approach of current design codes have been established by considering the behaviors of structural concrete without D-regions and two-dimensional concrete structures with D-regions, respectively, it is inappropriate to apply the methods to the pile caps dominated by 3-dimensional structural behavior with disturbed stress regions. In this study, the refined 3-dimensional strut-tie models, which consider the strength characteristics of 3-dimensional concrete struts and nodal zones and the load-carrying capacity of concrete ties in tension regions, are proposed for the rational analysis and design of pile caps. To examine the validity of the proposed models and to verify the necessity of appropriate constituent elements for describing 3-dimensional structural behavior and load-transfer mechanism of pile caps, the ultimate strength of 78 reinforced concrete pile caps tested to failure was examined by the proposed models along with the sectional and strut-tie model methods of current design codes.

• Journal of the Korea Concrete Institute
• /
• v.16 no.6 s.84
• /
• pp.841-850
• /
• 2004

The modeling of uplift pressure within dam, on the foundation on which it was constructed, and on the interface between the dam and foundation is a critical aspect in the analysis of concrete gravity dams, i.e. crack stability in concrete dam can correctly be predicted when uplift pressures are accurately modelled. Current models consider a uniform uplift distribution, but recent experimental results show that it varies along the crack faces and the procedures for modeling uplift pressures are well established for the traditional hand-calculation methods, but this is not the case for finite element (FE) analysis. In large structures, such as dams, because of smaller size of the fracture process zone with respect to the structure size, limited errors should occur under the assumptions of linear elastic fracture mechanics (LEFM). In this paper, the fracture behaviour of concrete gravity dams mainly subjected to uplift Pressure at the crack face was studied. Triangular type, trapezoidal type and parabolic type distribution of the uplift pressure including uniform type were considered in case of evaluating stress intensity factor by surface integral method. The effects of body forces, overtopping pressures are also considered and a parametric study of gravity dams under the assumption of LEFM is performed.