• Advances in nano research
• /
• v.12 no.2
• /
• pp.117-137
• /
• 2022

Effect of thickness stretching on free vibration, bending and buckling behavior of carbon nanotubes reinforced composite (CNTRC) laminated nanoplates rested on new variable elastic foundation is investigated in this paper using a developed four-unknown quasi-3D higher-order shear deformation theory (HSDT). The key feature of this theoretical formulation is that, in addition to considering the thickness stretching effect, the number of unknowns of the displacement field is reduced to four, and which is more than five in the other models. Two new forms of CNTs reinforcement distribution are proposed and analyzed based on cosine functions. By considering the higher-order nonlocal strain gradient theory, microstructure and length scale influences are included. Variational method is developed to derive the governing equation and Galerkin method is employed to derive an analytical solution of governing equilibrium equations. Two-dimensional variable Winkler elastic foundation is suggested in this study for the first time. A parametric study is executed to determine the impact of the reinforcement patterns, nonlocal parameter, length scale parameter, side-t-thickness ratio and aspect ratio, elastic foundation and various boundary conditions on bending, buckling and free vibration responses of the CNTRC plate.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.1A
• /
• pp.9-18
• /
• 2011

Generally the connection of structural members is assumed as hinge, rigid and semi-rigid connections. The exact tangent stiffness matrix of a semi-rigid frame element is newly derived using the stability functions considering shear deformations. Also, linearized elastic- and geometric-stiffness matrices of shear deformable semi-rigid frame are newly proposed. For the exact stiffness matrix, an accurate displacement field is introduced by equilibrium equation for beam-column under the bending and the axial forces. Also, stability functions considering sway deformation and force-displacement relations with elastic rotational spring on ends are defined. In order to illustrate the accuracy of this study, various numerical examples are presented and compared with other researcher's results. Lastly, shear deformation and semi-rigid effects on buckling behaviors of structure are parametrically investigated.

• Journal of Korean Society of Steel Construction
• /
• v.20 no.1
• /
• pp.81-92
• /
• 2008

Generally the connection of members is defined as hinge or rigid. But, real joints on structure have to be considered semi-rigid connections because this permits relative rotation for members on joints. The purpose of this study is to derive a generalized tangential stiffness matrix of frames with semi-rigid connections and to develop a buckling analysis program. For the exact stiffness matrix, an accurate displacement field is introduced using an equilibrium equation for beam-columns under the bending and axial forces. Also, stability functions that consider sway deformation and force-displacement relations with rotational spring on ends were defined. In order to illustrate the accuracy of this study and the characteristics of semi-rigid for system buckling load, samples of angle-, portal- and 3-story frames with semi-rigid connections are presented, where the proposed approach is found to be in excellent agreement with other research results. Meanwhile, the application of codes such as Eurocode 3 and LRFD led to significant inaccuracies.

• Journal of the Korean Geotechnical Society
• /
• v.27 no.3
• /
• pp.75-83
• /
• 2011

Emerging issues related with fully coupled Thermo-Hydro-Mechanical (THM) behavior of unsaturated soil demand the development of a numerical tool in diverse geo-mechanical and geo-environmental areas. This paper presents general governing equations for coupled THM processes in unsaturated porous media. Coupled partial differential equations are derived from three mass balances equations (solid, water, and air), energy balance equation, and force equilibrium equation. With Galerkin formulation and time integration of these governing equations, finite element code is developed to find nonlinear solution of four main variables (displacement-u, gas pressure-$P_g$), liquid pressure-$P_1$), and temperature-T) using Newton's iterative scheme. Three cases of numerical simulations are conducted and discussed: one-dimensional drainage experiments (u-$P_g-P_1$), thermal consolidation (u-$P_1$-T), and effect of pile on surrounding soil due to surface temperature variation (u-$P_1$-T).

• Journal of the Korean Geotechnical Society
• /
• v.30 no.9
• /
• pp.67-75
• /
• 2014

Methane gas hydrate which is considered energy source for the next generation has an urgent need to develop reliable numerical simulator for coupled THM phenomena in the porous media, to minimize problems arising during the production and optimize production procedures. International collaborations to improve previous numerical codes are in progress, but they still have mismatch in the predicted value and unstable convergence. In this paper, FEM code for fully coupled THM phenomena is developed to analyze methane hydrate dissociation in the porous media. Coupled partial differential equations are derived from four mass balance equations (methane hydrate, soil, water, and hydrate gas), energy balance equation, and force equilibrium equation. Five main variables (displacement, gas saturation, fluid pressure, temperature, and hydrate saturation) are chosen to give higher numerical convergence through trial combinations of variables, and they can analyze the whole region of a phase change in hydrate bearing porous media. The kinetic model is used to predict dissociation of methane hydrate. Developed THM FEM code is applied to the comparative study on a Masuda's laboratory experiment for the hydrate production, and verified for the stability and convergence.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.23 no.1
• /
• pp.95-110
• /
• 2010

A numerical model to simulate bond-slip behavior of composite beam bridges is introduced in this paper. Assuming a linear bond stress-slip relation along the interface between the slab and girder, the slip behavior is implemented into a finite element formulation. Adopting the introduced model, the slip behavior can be taken account even in a beam element which is composed of both end nodes only. Governing equation of the slip behavior, based on the linear partial interaction theory, can be determined from the force equilibrium and a constant curvature distribution across the section of a composite beam. Since the governing equation for the slip behavior requires the moment values at both end nodes, the piecewise linear distribution of the constant bending moment in an element is assumed. Analysis results by the model are compared with numerical results and experimental values, and load-displacement relations of composite beams were then evaluated to verify the validity of the proposed model.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.1
• /
• pp.25-35
• /
• 2000

The superstructure of general bridge like slab bridge and slab on girder bridge is composed of elastically supported isotropic plate. The objective of this study is to develop the new analysis method for elastically supported plate with general edge beam or girder(boundaries) under arbitrary out of plane loading. The displacement solutions for the macro-element of plate and beam are obtained by solving for the unknown interactive forces and moments at the beam or nodal line locations after satisfying equilibrium equation along the nodal line. The displacement functions for macro-elements ate proposed in single Fourier series using harmonic analysis, and the equilibrium equations of nodal line are composed by using slope-deflection method. The proposed analysis method is programmed by MS-Fortran and can be applied to all types of isotropic decks with bridge-type boundaries. Numerical examples involving elastically supported plates with various aspect ratio, loading cases, and bridge-type boundary conditions are presented to demonstrate the accuracy of this program. The major advantage of this new analysis method is the development of a simple solution algorithm, leads to obtain rapidly responses of bridge deck system. This proposed method can be used in parametric study of behavior of bridge decks.

• Journal of Korean Society of Coastal and Ocean Engineers
• /
• v.7 no.2
• /
• pp.156-162
• /
• 1995

A numerical method for a 2-D oil boom model considering the flexibility of skirt has been developed The neater is assumed rigid and the skirt is tensioned membrane having a point mass at its end The fluid motion is potential. The kinematic condition which demands the continuity of the displacement is imposed at the joint between the floater and the skirt. The dynamic condition for the point mass is imposed at the bottom end of the skirt. The numerical method is based on the Green's function method in the frame of linear potential theory. It finds it's solution simultaneously from the total system of three equations, integral equation, the equation of motion of the floater and the equilibrium equation of the deformation of the skirt. Integral equation is derived by applying the Green's theorem to radiation potential and Green's function. Proper descretization of those three equations leads to the system of a linear algebraic equation. Due to the flexibility of skirt the motion of floater can be diminished in some range of wave frequency and furthermore the mechanism of resonance of the oil boom can be changed. The motion responses of various oil booms have been compared varying the length of the skirt and the point mass. The numerical method has been validated indirectly from the good correspondence between the motion responses of the flexible skirt model and the rigid skirt model in low frequency limit.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.12 no.4
• /
• pp.71-79
• /
• 1992

In general, the accelerogram of earthquake ground motion or the accelerogram obtained from dynamic tests contain various errors. In these errors of the accelerograms, there are instrumental errors(magnitude and phase distortion) due to the response characteristics of accelerometer and the digitizing error concentrated in low and high frequency components and random errors. Then, these errors may be detrimental to the results of data processing and dynamic analysis. An efficient method which can correct the errors of the accelerogram is proposed in this study. The correction of errors can be accomplished through four steps as followes ; 1) using an interpolation method a data form appropriate to the error correction is prepared, 2) low and high frequency errors of the accelerogram are removed by band-pass filter between prescribed frequency limits, 3) instrumental errors are corrected using dynamic equilibrium equation of the accelerometer, 4) velocity and displacement are obtained by integrating corrected accelerogram. Presently, infinite impulse response(IIR) filter and finite impulse response (FIR) filter are generally used as band-pass filter. In the proposed error correction procedure, the deficiencies of FIR filter and IIR filter are reduced and, using the properties of the differentiation and the integration of Fourier transform, the accuracy of instrument correction and integration is improved.

• Korean Journal of Soil Science and Fertilizer
• /
• v.44 no.1
• /
• pp.15-21
• /
• 2011

Transport of heavy metals such as Cd is affected by several rate-limiting processes including adsorption and desorption by exchange reactions in soils. In this study, column transport and batch kinetic experiments were performed to assess Cd mobility in a double-layered soil with a reclaimed saline and sodic soil (SSS) as top soil and macroporous granule (MPG) as a bottom layer. For individual soil layer having different physical and chemical properties, Cd was considered to be nonlinear reactivity with the soil matrix in layered soils. The dispersive equation for reactive solutes was solved with three types of boundary conditions for the interface between soil layers. The adsorption of Cd with respect to the saline-sodic sandy loam and the MPG indicated that the nature of the sites or the mechanisms involved in the sorption process of Cd was different and the amounts of Cd for both of samples increases with increasing amounts of equilibrium concentration whereas the amount of Cd adsorbed in saline-sodic sandy loam soil was higher than that in MPG. The results of breakthrough curve indicating relative Cd retardation accompanied by layer material and sequence during leaching showed that the number of pore volumes to reach the maximum relative concentration of 1 increased in the order of MPG, SSS, and double layer of SSS-MPG. Breakthrough curves (BTCs) from column experiments were well predicted with our double-layered model where independently derived solute physical and retention parameters were implemented.