• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2003.04a
• /
• pp.235-242
• /
• 2003

The solids and the fluids in porous media have a relative velocity to each other. Due to physically and chemically different material properties and their relative velocity, the behavior of saturated porous media is extremely complicated. Thus, in order to describe and clarify the deformation behavior of saturated porous media, constitutive models for deformation of porous media coupling several effects such as flow of the fluids or thermodynanical change need to be developed in frame of Arbitrary Lagrangian Eulerian (ALE) description. The aim of ALE formulations is to maximize the advantages of Lagrangian and Eulerian elements, and to minimize the disadvantages. Therefore, this method is appropriate for the analysis of porous media that are considered for the behavior of the solids and the fluids. In this work, governing equations of porous media based on ALE description are obtained from governing equations in frame of updated Lagrangian description. Then, weak forms of these equations are derived using arbitrary weighting functions.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.04a
• /
• pp.68-74
• /
• 1998

A finite element analysis for physical phenomenon which are governed by parabolic equation, has some inefficiencies caused by much computational time and large storage space. In this paper, a comparative study is performed to suggest the best efficient transient analysis algorithms for parabolic equations. First, the general finite element analysis techniques are summarized in views of formulation procedures, treatments of convection terms. and time stepping methods. Results of several combinations applied to one dimensional convection-diffusion equation and Burger equation are represented and compared using some criteria such as accuracy, stability, and computational time. Through the results, some guidelines to select a algorithm for solving parabolic equations are proposed for diffusion dominant and convection dominant cases. Finally applicability of two dimensional extension of the result is also discussed.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1998.10a
• /
• pp.109-116
• /
• 1998

The method of vibration analysis used is the one developed by the senior author. He developed and reported, in 1974, a simple but exact method of calculating the natural frequency of beam and tower structures with irregular cross-sections and attached mass/masses. Since 1989, this method has been extended to two-dimensional problems with several types of given conditions and has been reported at several international conferences. This method uses the deflection influence surfaces. The finite difference method is used for this purpose, in this paper. In order to reduce the pivotal points required, the three simultaneous partial differential equations of equilibrium with three dependent variables, w, M$_{x}$, and $M_{y}$, are used instead of the one forth order partial differential equation. By neglecting the M$_{x}$ terms, the size of the matrices needed to solve the resulting linear equations are reduced to two thirds of the "non-modified" equations.tions.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 1999.10a
• /
• pp.421-428
• /
• 1999

The elastic critical load of a slender compression member plays an important role when the proper design of that member is required. For the tapered compression members, however, there are cases when the conventional neutral equililbrium or energy method can't be applied to the determination of critical loads of those members. In this paper, finite element method is applied to the approximate determination of the symmetrically tapered bars. Here in this paper, the bars are assumed to take sinusoidally changing shapes along their axes. The parameters considered in this study are taper parameter, $\alpha$ and the sectional property parameter, m. The computed results by finite element method are represented in the forms of algebraic equations. Regression technique is employed to determine the coefficients of algebraic equations. The critical loads estimated by the proposed algebraic equations coincide fairly well with those of finite element method.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2000.10a
• /
• pp.299-306
• /
• 2000

The elastic critical load of a slender compression member plays an important role when the proper design of that member is required. For tapered compression members, however, there are cases when the conventional neutral equilibrium or energy method can't be applied to the determination of critical loads. In this paper, the finite element method is applied to the approximate determination of the linearly tapered members. In this paper, the bars are assumed to be tapered linearly along their axes. The parameters considered in this study are taper parameter, α and the sectional property parameter, m. The member ends are either hinged or fixed. The computed results using the finite element method are represented in the forms of algebraic equations. The regression technique is employed to determine the coefficients of the algebraic equations. Critical loads estimated by the proposed algebraic equations coincide flirty well with those employing the finite element method.

• Structural Engineering and Mechanics
• /
• v.29 no.4
• /
• pp.355-380
• /
• 2008

In this study, the new three-dimensional finite element analysis model of guideway structures considering ultra high-speed magnetic levitation train-bridge interaction, in which the various improved finite elements are used to model structural members, is proposed. The box-type bridge deck of guideway structures is modeled by Nonconforming Flat Shell finite elements with six DOF (degrees of freedom). The sidewalls on a bridge deck are idealized by using beam finite elements and spring connecting elements. The vehicle model devised for an ultra high-speed Maglev train is employed, which is composed of rigid bodies with concentrated mass. The characteristics of levitation and guidance force, which exist between the super-conducting magnet and guideway, are modeled with the equivalent spring model. By Lagrange's equations of motion, the equations of motion of Maglev train are formulated. Finally, by deriving the equations of the force acting on the guideway considering Maglev train-bridge interaction, the complete system matrices of Maglev train-guideway structure system are composed.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.6
• /
• pp.627-633
• /
• 2010

This paper presents a method for identifying the parameter set of inelastic constitutive equations, which is based on an Evolutionary Algorithm. The advantage of the method is that appropriate parameters can be identified even when the measured data are subject to considerable errors and the model equations are inaccurate. The design of experiments suited for the parameter identification of a material model by Chaboche under the uniaxial loading and stationary temperature conditions was first considered. Then the parameter set of the model was identified by the proposed method from a set of experimental data. In comparison to those by other methods, the resultant stress-strain curves by the proposed method correlated better to the actual material behaviors.

• Structural Engineering and Mechanics
• /
• v.36 no.5
• /
• pp.529-544
• /
• 2010

Nonlinear equations of structures are generally solved numerically by the iterative solution of linear equations. However, this iterative procedure diverges when the tangent stiffness is ill-conditioned which occurs near limit points. In other words, a major challenge with simple iterative methods is failure caused by a singular or near singular Jacobian matrix. In this paper, using the Newton-Raphson algorithm based on Davidenko's equations, the iterations can traverse the limit point without difficulty. It is argued that the propose algorithm may be both more computationally efficient and more robust compared to the other algorithm when tracing path through severe nonlinearities such as those associated with structural collapse. Two frames are analyzed using the proposed algorithm and the results are compared with the previous methods. The ability of the proposed method, particularly for tracing the limit points, is demonstrated by those numerical examples.

• Structural Engineering and Mechanics
• /
• v.73 no.1
• /
• pp.83-95
• /
• 2020

In this work, a mathematical model of beam-column system carrying a double eccentric end mass system is investigated, and solved analytically based on the exact solution analysis. The model considers the case in which the double eccentric end mass is a rigid storage tank containing fluid. Both Timoshenko and Bernoulli-Euler beam bending theories are considered. Equation of motion, general solution and boundary conditions for the present system model are developed and presented in dimensional and non-dimensional format. Several important non-dimensional design parameters are introduced. Symbolic and/or explicit formulae of the frequency and mode shape equations are formulated. To the authors knowledge, the present reduced closed form symbolic and explicit frequency equations have not appeared in literature. For different applications, the results are validated using commercial finite element package, namely ANSYS. The beam-column system investigated in this paper is significant for many engineering applications, especially, in mechanical and structural systems.

• Korean Management Science Review
• /
• v.18 no.1
• /
• pp.119-133
• /
• 2001

It is important to find the equilibrium level of real interest rate for it affects real and financial sector of economy. However, it is difficult to find the equilibrium level because like the most macroeconomic model the real interest model has parameter instability problem caused by structural change and it is supported by various theories and definitions. Hence, in order to cover these problems structural change detection model of real interest rate is developed to combine the real interest rate equilibrium model and the procedure to detect structural change points. 3 equations are established to find various effects of other interest-related macroeconomic variables and from each equation, structural changes are found. Those structural change points are consistent with common expectation. Oil Crisis (December, 1987), the starting point of Economic Stabilization Policy (January, 1982), the starting point of capital liberalization (January, 1988), the starting and finishing points of Interest deregulation (January, 1992 and December, 1994), Foreign Exchange Crisis (December, 1977) are detected as important points. From the equation of fisher and real effects, real interest rate level is estimated as 4.09% (October, 1988) and dependent on the underlying model, it is estimated as 0%∼13.56% (October, 1988), so it varies so much. It is expected that this result is connected to the large scale simultaneous equations to detect the parameter instability in real time, so induces the flexible economic policies.