• Transactions of the Korean Society of Mechanical Engineers
• /
• v.16 no.8
• /
• pp.1429-1437
• /
• 1992

A new finite element formulation based on the relaxation type hereditary integral is presented for a time-domain analysis of isotropic, linear viscoelastic problems. The semi-discrete variational approximation and elastic-viscoelastic correspondence principle are used in the theoretical development of the proposed method. In a time-stepping procedure of final, linear algebraic system equations, only a small additional computation for past history is required since the equivalent stiffness matrix is constant. The viscoelasticity matrices are derived and the stress computation algorithm is given in matrix form. The effect of time increment and Gauss point numbers on the numerical accuracy is examined. Two dimensional numerical examples of plane strain and plane stress are solved and compared with the analytical solutions to demonstrate the versatility and accuracy of the present method.

• Transactions of the Korean Society for Noise and Vibration Engineering
• /
• v.20 no.10
• /
• pp.923-930
• /
• 2010

Finite element models of dynamic systems can be updated in two stages. In the first stage, mass and stiffness matrices are updated neglecting damping, and in the second stage, damping matrices are estimated with the mass and stiffness matrices fixed. Three methods to estimate damping matrices for this purpose are proposed in this paper. The methods include one for proportional damping systems and two for non-proportional damping systems. Method 1 utilizes orthogonality of normal modes and estimates damping matrices using the modal parameters extracted from the measured responses. Method 2 estimates damping matrices from impedance matrices which are the inverse of FRF matrices. Method 3 estimates damping using the equation which relates a damping matrix to the difference between the analytical and measured FRFs. The characteristics of the three methods are investigated by applying them to simulated discrete system data and experimental cantilever beam data.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.23 no.1
• /
• pp.19-29
• /
• 2019

Recent earthquakes in Korea caused some damages to stone pagodas and thereby awakened the importance of earthquake preparedness. Korean stone pagodas which have been built with very creative style of material use and construction method are worthy of world heritage. Each stone pagoda consists of three parts: top; body; and base. However each tower is uniquely defined by its own features, which makes it more difficult to generalize the seismic assessment method for stone pagodas. This study has focused on qualitative preliminary evaluation of stone pagodas that enables us to compare the relative seismic performance across major aspects among many various Korean pagodas. Specifically an analytical model for multi-block stone pagodas is to be proposed upon the investigation of structural characteristics of stone pagoda and their dynamic behavior. A strategy for seismic evaluation of heritage stone pagodas is to be established and major evaluation factors appropriate for the qualitative evaluation are identified. The evaluation factors for overall seismic resisting behavior of stone pagodas are selected based on the dynamic motions of a rigid block and its limit state. Numerical simulation analysis using discrete element method is performed to analyze the sensitivity of each factor to earthquake and discuss some effects on seismic performance.

• Geomechanics and Engineering
• /
• v.5 no.5
• /
• pp.379-399
• /
• 2013

In practice, analysis of laterally loaded piles is carried out using beams on non-linear Winkler springs model (often known as p-y method) due to its simplicity, low computational cost and the ability to model layered soils. In this approach, soil-pile interaction along the depth is characterized by a set of discrete non-linear springs represented by p-y curves where p is the pressure on the soil that causes a relative deformation of y. p-y curves are usually constructed based on semi-empirical correlations. In order to construct API/DNV proposed p-y curve for clay, one needs two values from the monotonic stress-strain test results i.e., undrained strength ($s_u$) and the strain at 50% yield stress (${\varepsilon}_{50}$). This approach may ignore various features for a particular soil which may lead to un-conservative or over-conservative design as not all the data points in the stress-strain relation are used. However, with the increasing ability to simulate soil-structure interaction problems using highly developed computers, the trend has shifted towards a more theoretically sound basis. In this paper, principles of Mobilized Strength Design (MSD) concept is used to construct a continuous p-y curves from experimentally obtained stress-strain relationship of the soil. In the method, the stress-strain graph is scaled by two coefficient $N_C$ (for stress) and $M_C$ (for strain) to obtain the p-y curves. $M_C$ and $N_C$ are derived based on Semi-Analytical Finite Element approach exploiting the axial symmetry where a pile is modelled as a series of embedded discs. An example is considered to show the application of the methodology.

• Journal of Korean Society of Steel Construction
• /
• v.24 no.3
• /
• pp.289-299
• /
• 2012

The goal of this paper is to apply the multi-step Taylor method to a space truss, a non-linear discrete dynamic system, and analyze the non-linear dynamic response and unstable behavior of the structures. The accurate solution based on an analytical approach is needed to deal with the inverse problem, or the dynamic instability of a space truss, because the governing equation has geometrical non-linearity. Therefore, the governing motion equations of the space truss were formulated by considering non-linearity, where an accurate analytical solution could be obtained using the Taylor method. To verify the accuracy of the applied method, an SDOF model was adopted, and the analysis using the Taylor method was compared with the result of the 4th order Runge-Kutta method. Moreover, the dynamic instability and buckling characteristics of the adopted model under step excitation was investigated. The result of the comparison between the two methods of analysis was well matched, and the investigation shows that the dynamic response and the attractors in the phase space can also delineate dynamic snapping under step excitation, and damping affects the displacement of the truss. The analysis shows that dynamic buckling occurs at approximately 77% and 83% of the static buckling in the undamped and damped systems, respectively.

• Geomechanics and Engineering
• /
• v.23 no.5
• /
• pp.455-466
• /
• 2020

This paper presents analysis of the interaction between tunnel and Qanat with a particular interest for the optimization of Qanat shape using the discrete element code, PFC2D, and the results will be compared with the FEM results of PLAXIS2D. For these concerns, using software PFC2D based on Discrete Element Method (DEM), a model with dimension of 100m ^⁎ 100 m was prepared. A circular tunnel with dimension of 9 m was situated 20 m below the ground surface. Also one Qanat was situated perpendicularly above the tunnel roof. Distance between Qanat center and ground surface was 8 m. Five different shapes for Qanat were selected i.e., square, semi-circular, vertical ellipse, circular and horizontal ellipse. Confining pressure of 5 MPa was applied to the model. The vertical displacement of balls situated in ground surface was picked up to measure the ground subsidence. Also two measuring circles were situated at the tunnel roof and at the Qanat roof to check the vertical displacements. The properties of the alluvial soil of Tehran city are: γ_dry=19 (KN/㎥), E= 750 (kg/㎠), ν=0.35, c=0.3(kg/㎠), φ=34°. In order to validate the DEM results, a comparison between the numerical results (obtained in this study) and analytical and field monitoring have been done. The PFC2D results are compared with the FEM results. The results shows that when Qanat has rectangular shape, the tensile stress concentration at the Qanat corners has maximum value while it has minimum value for vertical ellipse shape. The ground subsidence for Qanat rectangular shape has maximum value while it has minimum value for ellipse shape of Qanat. The vertical displacements at the tunnel roof for Qanat rectangular shape has maximum value while it has minimum value for ellipse shape of Qanat. Historical shape of Qante approved the finding of this research.

• The KIPS Transactions:PartC
• /
• v.12C no.3 s.99
• /
• pp.409-418
• /
• 2005

This paper presents an efficient multicast method based on a locality in mobile IPv6 networks. We exploit the repetitive movement pattern of mobile node to reduce the total number of experience of graft and join procedure. We defined the locality scope by a movement pattern. While the network is included in the locality scope, the network should maintain a multicast tree even when the mobile node moves to the other network. In this way, the mobile host can receive a multicast service without a delay when it moves to the network in the locality scope later. We compare our scheme with existing schemes under the total signaling cost and the service delay time by using a discrete analytical model for cost analysis. Analytical results demonstrated that the total signaling cost and service delay time was significantly reduced through our proposed scheme.

• Geomechanics and Engineering
• /
• v.11 no.1
• /
• pp.1-39
• /
• 2016

The finite element method (FEM), discrete element method (DEM), and Discontinuous deformation analysis (DDA) are among the standard numerical techniques applied in computational geo-mechanics. However, in some cases there no possibility for modelling by traditional finite analytical techniques or other mesh-based techniques. The solution presented in the current study as a completely Lagrangian and mesh-free technique is smoothed particle hydrodynamics (SPH). This method was basically applied for simulation of fluid flow by dividing the fluid into several particles. However, several researchers attempted to simulate soil-water interaction, landslides, and failure of soil by SPH method. In fact, this method is able to deal with behavior and interaction of different states of materials (liquid and solid) and multiphase soil models and their large deformations. Soil indicates different behaviors when interacting with water, structure, instrumentations, or different layers. Thus, study into these interactions using the mesh based grids has been facilitated by mesh-less SPH technique in this work. It has been revealed that the fast development, computational sophistication, and emerge of mesh-less particle modeling techniques offer solutions for problems which are not modeled by the traditional mesh-based techniques. Also it has been found that the smoothed particle hydrodynamic provides advanced techniques for simulation of soil materials as compared to the current traditional numerical methods. Besides, findings indicate that the advantages of applying this method are its high power, simplicity of concept, relative simplicity in combination of modern physics, and particularly its potential in study of large deformations and failures.

• Structural Engineering and Mechanics
• /
• v.11 no.6
• /
• pp.591-604
• /
• 2001

In this paper, a finite element with embedded displacement discontinuity which eliminates the need for remeshing of elements in the discrete crack approach is applied for the progressive fracture analysis of concrete structures. A finite element formulation is implemented with the extension of the principle of virtual work to a continuum which contains internal displacement discontinuity. By introducing a discontinuous displacement shape function into the finite element formulation, the displacement discontinuity is obtained within an element. By applying either a nonlinear or an idealized linear softening curve representing the fracture process zone (FPZ) of concrete as a constitutive equation to the displacement discontinuity, progressive fracture analysis of concrete structures is performed. In this analysis, localized progressive fracture simultaneous with crack closure in concrete structures under mixed mode loading is simulated by adopting the unloading path in the softening curve. Several examples demonstrate the capability of the analytical technique for the progressive fracture analysis of concrete structures.

• Structural Engineering and Mechanics
• /
• v.55 no.3
• /
• pp.655-674
• /
• 2015

In this paper, a new and simplified method is presented in which the natural frequencies of the uniform and non-uniform beams are calculated through simple mathematical relationships. The various vibration problems such as: Rayleigh beam under variable axial force, axial vibration of a bar with and without end discrete spring, torsional vibration of a bar with an attached mass moment of inertia, flexural vibration of the beam with laterally distributed elastic springs and also flexural vibration of the beam with effects of viscose damping are investigated. The governing differential equations are first obtained and then; according to a harmonic vibration, are converted into single variable equations in terms of location. Through repetitive integrations, the governing equations are converted into weak form integral equations. The mode shape functions of the vibration are approximated using a power series. Substitution of the power series into the integral equations results in a system of linear algebraic equations. The natural frequencies are determined by calculation of a non-trivial solution for system of equations. The efficiency and convergence rate of the current approach are investigated through comparison of the numerical results obtained with those obtained from other published references and results of available finite element software.