• Proceedings of the KSME Conference
• /
• 2003.11a
• /
• pp.1346-1351
• /
• 2003

There is few research about contact problem in SPH because it is primarily suitable to analyze the large deformation problem. However, an elasto-plastic problem with small deformation need to be considered about contact characteristics. The numerical formulating methods for SPH is induced to be able to obtain solutions based on a variational method in contact problem. The contact algorithm presented is applied to the elastic impact problem in 1D and 2D. The results show thai an imaginary tension and a numerical instability which happen in impacting between different materials can be removed and contact forces which could not have been calculated are able to obtain.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.450-455
• /
• 2004

There is few research about contact problem for a rigid surface with an arbitrary shape in SPH. The variational equation based on the virtual work principle is derived and its solution is obtained by the penalty method. It is proposed a new method that can determine the parameters for a penetration and a penetration rate used in the penalty method. The reproducing condition is adopted to correct the deficiency of kernel on the boundary. In order to calculate a penetration of particles, after checking boundary particles for deformable body boundary normal vectors were determined on the rigid surface. Numerical simulations for models which have rigid surface with an arbitrary shape were conducted to validate the proposed method in 2D. The results of those analysis represent that the contact algorithm proposed in this study works properly.

• International Journal of Aeronautical and Space Sciences
• /
• v.18 no.3
• /
• pp.467-473
• /
• 2017

This paper focuses on the investigation of three-dimensional (3D) warping effect on the stiffness constants of composite beams with closed section profiles. A finite element (FE) cross-sectional analysis is developed based on the Reissner's multifield variational principle. The 3D in-plane and out-of-plane warping displacements, and sectional stresses are approximated as linear functions of generalized sectional stress resultants at the global level and as FE shape functions at the local sectional level. The classical elastic couplings are taken into account which include transverse shear and Poisson deformation effects. A generalized Timoshenko level $6{\times}6$ stiffness matrix is computed for closed section composite beams with and without warping. The effect of neglecting the 3D warping on stiffness constants is shown to be significant indicating large errors as high as 93.3%.

• Nonlinear Functional Analysis and Applications
• /
• v.27 no.2
• /
• pp.405-432
• /
• 2022

Cracks in beams and shallow arches are modeled by massless rotational springs. First, we introduce a specially designed linear operator that "absorbs" the boundary conditions at the cracks. Then the equations of motion are derived from the first principles using the Extended Hamilton's Principle, accounting for non-conservative forces. The variational formulation of the equations is stated in terms of the subdifferentials of the bending and axial potential energies. The equations are given in their abstract (weak), as well as in classical forms.

• Journal of Korean Society of Transportation
• /
• v.24 no.2 s.88
• /
• pp.169-178
• /
• 2006

Traffic assignment has been used both for predicting travel demands and for evaluating the tools for alleviating congestion on road network in advance. Some assignment models have been proposed such as equivalent mathematical minimization method, variational inequality problem, nonlinear complementary problem and fixed point method, in following the principle of Wardrop (1952) that no driver can not Improve his travel cost by unilaterally changing his route. Recently Jin(2005a) presented a traffic assignment model based on dynamic process. This paper proposes a solution algorithm for the model of Jin and assesses the performances. Compared to the Frank-Wolfe method, which has been wildly used for solving the existing assignment models, the proposed algorithm is expected to be more efficient because it does not need to evaluate the objective function. Two numerical examples are used for assessing the algorithm, and they show that the algorithm converges to user equilibrium of Wardrop.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.27 no.3
• /
• pp.173-182
• /
• 2014

The mixed convolved action principle provides a new rigorous weak variational formalism for a broad range of initial boundary value problems in mathematical physics and mechanics in terms of mixed formulation, convolution, and fractional calculus. In this paper, its potential in the development of numerical methods for transient problems in various dynamical systems when adopting temporally second order approximation is investigated. For this, the classical single-degree-of-freedom linear elastic dynamical systems are primarily considered to investigate computational characteristics of the developed algorithms. For the undamped system, all the developed algorithms are symplectic with respect to the time step. For the damped system, they are shown to be accurate with good convergence characteristics.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.21 no.2
• /
• pp.153-160
• /
• 2008

In this study, we propose a new three-node hybrid-mixed curved beam element with the relaxed-equiribrium stress functions for static analysis. The proposed element considering shear deformation is based on the Hellinger-Reissner variational principle. The stress functions are carefully chosen from three important considerations: (i) all the kinematic deformation modes must be suppressed, and (ii) the spurious constraints must be removed in the limiting behaviors via the field-consistency, and (iii) the relaxed equilibrium conditions could be incorporated because it might be impossible to select the stress functions and parameters to fully satisfy both the equiribrium conditions and the suppression of kinematic deformation modes in the three-node curved beam hybrid-mixed formulation. Numerical examples confirm the superior and stable behavior of the proposed element regardless of slenderness ratio and curvature. Besides, the proposed element shows the outstanding performance in predicting the stress resultant distributions.

• Structural Engineering and Mechanics
• /
• v.30 no.3
• /
• pp.263-278
• /
• 2008

Structural pipe-in-pipe cross-sections have significant potential for application in offshore oil and gas production systems because of their property that combines insulation performance with structural strength in an integrated way. Such cross-sections comprise inner and outer thin walled pipes with the annulus between them fully filled by a selectable thick filler material to impart an appropriate combination of properties. Structural pipe-in-pipe cross-sections can exhibit several different collapse mechanisms and the basis of the preferential occurrence of one over others is of interest. This paper presents an elastic analyses of a structural pipe-in-pipe cross-section when subjected to external hydrostatic pressure. It formulates and solves the static and elastic buckling problem using the variational principle of minimum potential energy. The paper also investigates a simplified formulation of the problem where the outer pipe and its contact with the filler material is considered as a 'pipe on an elastic foundation'. Results are presented to show the variation of elastic buckling pressure with the relative elastic modulus of the filler and pipe materials, the filler thickness and the thicknesses of the inner and outer pipes. The range of applicability of the simplified 'pipe on an elastic foundation' analysis is also presented. A brief review of the types of materials that could be used as the filler is combined with the results of the analysis to draw conclusions about elastic buckling behaviour of structural pipe-in-pipe cross-sections.

• Advances in nano research
• /
• v.6 no.2
• /
• pp.163-182
• /
• 2018

Based on the Reissner mixed variational theorem (RMVT), the authors present a nonlocal Timoshenko beam theory (TBT) for the nonlinear free vibration analysis of multi-walled carbon nanotubes (MWCNT) embedded in an elastic medium. In this formulation, four different edge conditions of the embedded MWCNT are considered, two different models with regard to the van der Waals interaction between each pair of walls constituting the MWCNT are considered, and the interaction between the MWCNT and its surrounding medium is simulated using the Pasternak-type foundation. The motion equations of an individual wall and the associated boundary conditions are derived using Hamilton's principle, in which the von $K{\acute{a}}rm{\acute{a}}n$ geometrical nonlinearity is considered. Eringen's nonlocal elasticity theory is used to account for the effects of the small length scale. Variations of the lowest frequency parameters with the maximum modal deflection of the embedded MWCNT are obtained using the differential quadrature method in conjunction with a direct iterative approach.

• Computers and Concrete
• /
• v.2 no.2
• /
• pp.125-146
• /
• 2005

The objective of this paper is setting out, for a cable-stayed bridge with a curtain suspension, a method to determine the modes of vibration of the structure. The system of differential equations governing the vibrations of the bridge, derived by means of a variational formulation in a nonlinear field, is reported in Appendix C. The whole analysis results from the application of Hamilton's principle, using the expressions of potential and kinetic energies and of the virtual work made by viscous damping forces of the various parts of the bridge (Monaco and Fiore 2003). This paper focuses on the equation concerning the transversal motion of the girder of the cable-stayed bridge and in particular on its final form obtained, restrictedly to the linear case, neglecting some quantities affecting the solution in a non-remarkable way. In the hypotheses of normal mode of vibration and of steady-state, we propose the resolution of this equation by a particular method based on a numerical approach. Respecting the boundary conditions, we derive, for each mode of vibration, the corresponding frequency, both natural and damped, the shape-function of the girder axis and the exponential function governing the variability of motion amplitude in time. Finally the results so obtained are compared with those deriving from the dynamic analysis performed by a finite elements calculation program.