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

• Journal of Mechanical Science and Technology
• /
• v.17 no.1
• /
• pp.77-84
• /
• 2003

Virtual work in multibody systems is frequently expressed as the inner product of the virtual displacement and the resultant force at the centroid. But provided that the resultant force is converted into the equipollent forces there is no restriction on where the analysis reference point is placed. There are basically three candidate points : the centroid, joint point and the instant global origin. The traditional fully recursive formulation uses the centroid, but the present work verifies that the instant global origin always shows better efficiency (e.g. 86% CPU time of the centroid for quarter car model) and joint point shows the efficiency between that of the centroid and the instant global origin. A discussion on how important it is to define the analysis reference point properly in a fully recursive formulation is also presented.

• Steel and Composite Structures
• /
• v.24 no.5
• /
• pp.523-536
• /
• 2017

A layerwise (LW) formulation based on the Galerkin method is presented to investigate the three-dimensional stress state in long sandwich plate which is subjected to tension force and pure bending moment. Based on the Galerkin method and the LW discretization approach, the equilibrium equations of elasticity for the long plate are written in the weak form and discretized through the thickness of the plate. The discretized equations are written in terms of displacement components of the numerical layers. The governing equations of the plate are solved analytically for the free edge boundary conditions. The distribution of stress state especially the 3D stress state in the vicinity of the edges of the sandwich plate which is subjected to tension and pure bending is studied. In order to increase the accuracy, the out of plane stresses are obtained by integrating the equilibrium equations of elasticity. The convergence and accuracy of the predictions are studied and various numerical results are presented for distribution of the in-plane and out of plane stresses in symmetric and un-symmetric sandwich plates.

• Proceedings of the KSME Conference
• /
• 2001.06a
• /
• pp.674-679
• /
• 2001

For a large deformation nonlinear dynamic analysis of two-dimensional frictional contact, the linear complementarity formulation combined with a linearization is used. The solution procedure is based on the total Lagrangian formulation with a predictor and corrector scheme. For contact searching, a hierarchical scheme with a circular territory is used. A second-order approximation of displacements is used to detect impact time and position. The formulation is illustrated by means of numerical examples.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1996.04a
• /
• pp.182-189
• /
• 1996

The focus of this study is on the problem of the design of structure of undetermined topology. This problem has been regarded as being the most challenging of structural optimization problems, because of the difficulty of allowing topology to change. Conventional approaches break down when element sizes approach to zero, due to stiffness matrix singularity. In this study, a novel nonlinear Programming formulation of the topology Problem is developed and examined. Its main feature is the ability to account for topology variation through zero element sizes. Stiffness matrix singularity is avoided by embedding the equilibrium equations as equality constraints in the optimization problem. Although the formulation is general, two dimensional plane elasticity examples are presented. The design problem is to find minimum weight of a plane structure of fixed geometry but variable topology, subject to constraints on stress and displacement. Variables are thicknesses of finite elements, and are permitted to assume zero sizes. The examples demonstrate that the formulation is effective for finding at least a locally minimal weight.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1994.10a
• /
• pp.121-128
• /
• 1994

A layerwise theory for the dynamic response of a laminated composite plate with integrated piezoelectric actuators and sensors subjected to both mechanical and electrical loadings is proposed. The formulation is derived form the variational principle with consideration for both total potential energy of the structures and the electrical potential energy of the piezoceramics. The governing equations of the present theory account for direct and converse effects of piezoelectrics, and layerwise variation of displacement field through the thickness of a laminate.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1991.10a
• /
• pp.33-39
• /
• 1991

This Paper presents a geometrically nonlinear behaviors of shell problems by using the three-dimensional curved shell element, which includs large displacements and large rotations. The standard formulation of the geometrically nonlinearity is restricted to the assumption of infinitesmal rotation increments. This standard formulation for the displacement function is numerically improved by considering the second order expansions of Tayler series. The nonlinear behaviors of the single and double curved shells are compared wi th the other results.

• Computational Structural Engineering
• /
• v.7 no.4
• /
• pp.103-111
• /
• 1994

Based on a variational principle of the consistent shear deformable discrete laminate theory derived in the companion paper Part I, a finite element procedure for the vibration analysis of laminated composite plates is presented. The present formulation takes the in-plane displacements of an arbitrary layer, the rotations of the cross section of each layer and transverse displacement of the plate as the state variables at a nodal point of finite element, resulting in total nodal degree of freedom of 2(n+l) +1 for the n-layered laminate. Thus, it allows to specify displacement boundary conditions of layer stretching and/or rotation of layer cross sections around the plate edge and/or lateral displacement. The developed procedure is applied to the free vibration problem for sandwich-type hybrid laminates composed of layers with drastically different material properties whose elasticity solutions are known. Comparison of analysis results with other FEM solutions showed that the present formulation yields better accuracy.

• Computational Structural Engineering
• /
• v.9 no.3
• /
• pp.179-186
• /
• 1996

A simple nonlinear programming(NLP) formulation for the optimal topology problem of structures is developed and examined. The NLP formulation is general, and can handle arbitrary objective functions and arbitrary stress, displacement constraints under multiple loading conditions. The formulation is based on simultaneous analysis and design approach to avoid stiffness matrix singularity resulting from zero sizing variables. By embedding the equilibrium equations as equality constraints in the nonlinear programming problem, we avoid constructing and factoring a system stiffness matrix, and hence avoid its singularity. The examples demonstrate that the formulation is effective for finding an optimal solution, and shown to be robust under a variety of constraints.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.3 no.1
• /
• pp.55-66
• /
• 1999

This paper presents an analytical prediction of the hysteresis behavior of reinforced concrete long column with rectangular section under the cyclic loading state. The mechanical characteristic of cracked concrete and reinforcing bar in concrete has been modeled, considering the bond effect between reinforcing bar and concrete, the effect of aggregate interlocking at crack surface and the stiffness degradation after the crack has taken place. The strength increase of concrete due to the lateral confining reinforcement has been also taken into account to model the confined concrete. The formulation of these models for concrete and reinforcing bar has been based on the smeared crack concept that the stress-strain relationship of reinforced concrete element would be defined using the average values. In addition to the material nonlinear properties, the algorithm for large displacement problem that may give an additional deformation has been formulated using total Lagrangian formulation. The analytically predicted behavior was compared with test result and they showed good agreement in overall behavior.