• Title/Summary/Keyword: Galerkin integration

Search Result 86, Processing Time 0.021 seconds

Multiscale Wavelet-Galerkin Method in General Two-Dimensional Problems (일반 형상의 2차원 영역에서의 멀티스케일 웨이블렛-갤러킨 기법)

  • Kim, Yun-Yeong;Jang, Gang-Won;Kim, Jae-Eun
    • Transactions of the Korean Society of Mechanical Engineers A
    • /
    • v.26 no.5
    • /
    • pp.939-951
    • /
    • 2002
  • We propose a new multiscale Galerkin method based on interpolation wavelets for two-dimensional Poisson's and plane elasticity problems. The major contributions of the present work are: 1) full multiresolution numerical analysis is carried out, 2) general boundaries are handled by a fictitious domain method without using a penalty term or the Lagrange multiplier, 3) no special integration rule is necessary unlike in the (bi-)orthogonal wavelet-based methods, and 4) an efficient adaptive scheme is easy to incorporate. Several benchmark-type problems are considered to show the effectiveness and the potentials of the present approach. is 1-2m/s and impact deformation of the electrode depends on the strain rate at that velocity, the dynamic behavior of the sinter-forged Cu-Cr is a key to investigate the impact characteristics of the electrodes. The dynamic response of the material at the high strain rate is obtained from the split Hopkinson pressure bar test using disc-type specimens. Experimental results from both quasi-static and dynamic compressive tests are Interpolated to construct the Johnson-Cook model as the constitutive relation that should be applied to simulation of the dynamic behavior of the electrodes. The impact characteristics of a vacuum interrupter are investigated with computer simulations by changing the value of five parameters such as the initial velocity of a movable electrode, the added mass of a movable electrode, the wipe spring constant, initial offset of a wipe spring and the virtual fixed spring constant.

FEM Numerical Formulation for Debris Flow (토석류 유동해석을 위한 유한요소 수식화)

  • Shin, Hosung
    • Journal of the Korean Geotechnical Society
    • /
    • v.30 no.10
    • /
    • pp.55-65
    • /
    • 2014
  • Recent researches on debris flow is focused on understanding its movement mechanism and building a numerical simulator to predict its behavior. However, previous simulators emulating fluid-like debris flow have limitations in numerical stability, geometric modeling and application of various boundary conditions. In this study, depth integration is applied to continuity equation and force equilibrium for debris flow. Thickness of sediment, and average velocities in x and y flow direction are chosen for main variables in the analysis, which improve numerical stability in the area with zero thickness. Petrov-Galerkin formulation uses a discontinuous test function of the weighted matrix from DG scheme. Presented mechanical constitutive model combines fluid and granular behaviors for debris flow. Effects on slope angle, inducing debris height, and bottom friction resistance are investigated for a simple slope. Numerical results also show the effect of embankment at the bottom of the slope. Developed numerical simulator can assess various risk factors for the expected area of debris flow, and facilitate embankment design in order to minimize damage.

On the local stability condition in the planar beam finite element

  • Planinc, Igor;Saje, Miran;Cas, Bojan
    • Structural Engineering and Mechanics
    • /
    • v.12 no.5
    • /
    • pp.507-526
    • /
    • 2001
  • In standard finite element algorithms, the local stability conditions are not accounted for in the formulation of the tangent stiffness matrix. As a result, the loss of the local stability is not adequately related to the onset of the global instability. The phenomenon typically arises with material-type localizations, such as shear bands and plastic hinges. This paper addresses the problem in the context of the planar, finite-strain, rate-independent, materially non-linear beam theory, although the proposed technology is in principle not limited to beam structures. A weak formulation of Reissner's finite-strain beam theory is first presented, where the pseudocurvature of the deformed axis is the only unknown function. We further derive the local stability conditions for the large deformation case, and suggest various possible combinations of the interpolation and numerical integration schemes that trigger the simultaneous loss of the local and global instabilities of a statically determined beam. For practical applications, we advice on a procedure that uses a special numerical integration rule, where interpolation nodes and integration points are equal in number, but not in locations, except for the point of the local instability, where the interpolation node and the integration point coalesce. Provided that the point of instability is an end-point of the beam-a condition often met in engineering practice-the procedure simplifies substantially; one of such algorithms uses the combination of the Lagrangian interpolation and Lobatto's integration. The present paper uses the Galerkin finite element discretization, but a conceptually similar technology could be extended to other discretization methods.

Application of a Convolution Method for the Fast Prediction of Wind-Induced Surface Current in the Yellow Sea and the East China Sea (표층해류 신속예측을 위한 회선적분법의 적용)

  • 강관수;정경태
    • Journal of Korean Society of Coastal and Ocean Engineers
    • /
    • v.7 no.3
    • /
    • pp.265-276
    • /
    • 1995
  • In this Paper, the Performance of the convolution method has been investigated as an effort to develop a simple system of predicting wind-driven surface current on a real time basis. In this approach wind stress is assumed to be spatially uniform and the effect of atmospheric pressure is neglected. The discrete convolution weights are determined in advance at each point using a linear three-dimensional Galerkin model with linear shape functions(Galerkin-FEM model). Four directions of wind stress(e.g. NE, SW, NW, SE) with unit magnitude are imposed in the model calculation for the construction of data base for convolution weights. Given the time history of wind stress, it is then possible to predict with-driven currents promptly using the convolution product of finite length. An unsteady wind stress of arbitrary form can be approximated by a series of wind pulses with magnitude of 6 hour averaged value. A total of 12 pulses are involved in the convolution product To examine the accuracy of the convolution method a series of numerical experiments has been carried out in the idealized basin representing the scale of the Yellow Sea and the East China Sea. The wind stress imposed varies sinusoidally in time. It was found that the predicted surface currents and elevation fields were in good agreement with the results computed by the direct integration of the Galerkin model. A model with grid 1/8$^{\circ}$ in latitude, l/6$^{\circ}$ in longitude was established which covers the entire region of the Yellow Sea and the East China Sea. The numerical prediction in terms of the convolution product has been carried out with particular attention on the formation of upwind flow in the middle of the Yellow Sea by northerly wind.

  • PDF

Dynamic Response and Vibration Characteristics of an Isolation Rail Track under a Traveling Mass (주행질량하의 방진 궤도레일의 동적응답 및 진동특성)

  • Oh, B.J.;Ryu, B.J.;Kim, J.H.;Lee, Y.S.
    • Transactions of the Korean Society for Noise and Vibration Engineering
    • /
    • v.21 no.4
    • /
    • pp.365-373
    • /
    • 2011
  • This paper presents the dynamic response and the vibration characteristics for a rail-track supported by discrete springs and dampers. Recently, automatic conveyer system, rail-track, rack-master system demand the soundproof facilities and vibration suppression measures in order to satisfy the strict environmental standards. The equations of motions of the dynamic characteristics for a vibration suppression rail-track under a traveling mass were derived by Galerkin's mode summation method considering gravity, centrifugal force, Coriolis force, inertia force of the moving mass, transverse inertia of the rail-track. Also, numerical results were calculated by Runge-Kutta integration method. In order to investigate vibration characteristics and dynamic responses, modal testing and measurement of the responses of the rail-track were performed. Through the experiment and numerical simulations, numerical results have a good agreement with experimental ones.

Analysis of 2-Dimensional Elasto-Plastic Stress by a Time-Discontinuous Variational Integrator of Hamiltonian (해밀토니안의 시간 불연속 변분적분기를 이용한 2차원 탄소성 응력파 해석)

  • Chol, S.S.;Huh, H.;Park, K.C.
    • Proceedings of the Korean Society for Technology of Plasticity Conference
    • /
    • 2008.10a
    • /
    • pp.263-266
    • /
    • 2008
  • This paper is concerned with the analysis of elasto-plastic stress waves in a mode I semi-infinite cracked solid subjected to Heaviside pulse load. This study adopts a time-discontinuous variational integrator based on Hamiltonian in order to reduce the numerical dispersive and dissipative errors. This also utilizes an integration scheme of the constitutive model with 2nd-order accuracy which is formulated on the strain space for a rate and temperature dependent material model. Finite element analyses of elasto-plastic stress waves are carried out in order to compare the accuracy between a conventional Galerkin method and the time- discontinuous variational integrator.

  • PDF

Nonlinear stability and bifurcations of an axially accelerating beam with an intermediate spring-support

  • Ghayesh, Mergen H.;Amabili, Marco
    • Coupled systems mechanics
    • /
    • v.2 no.2
    • /
    • pp.159-174
    • /
    • 2013
  • The present work aims at investigating the nonlinear dynamics, bifurcations, and stability of an axially accelerating beam with an intermediate spring-support. The problem of a parametrically excited system is addressed for the gyroscopic system. A geometric nonlinearity due to mid-plane stretching is considered and Hamilton's principle is employed to derive the nonlinear equation of motion. The equation is then reduced into a set of nonlinear ordinary differential equations with coupled terms via Galerkin's method. For the system in the sub-critical speed regime, the pseudo-arclength continuation technique is employed to plot the frequency-response curves. The results are presented for the system with and without a three-to-one internal resonance between the first two transverse modes. Also, the global dynamics of the system is investigated using direct time integration of the discretized equations. The mean axial speed and the amplitude of speed variations are varied as the bifurcation parameters and the bifurcation diagrams of Poincare maps are constructed.

Nonlinear Dynamic Response of Cantilevered Carbon Nanotube Resonator by Electrostatic Excitation (정전기력 가진에 의한 외팔보형 탄소나노튜브 공진기의 비선형 동적 응답)

  • Kim, Il-Kwang;Lee, Soo-Il
    • Transactions of the Korean Society for Noise and Vibration Engineering
    • /
    • v.21 no.9
    • /
    • pp.813-819
    • /
    • 2011
  • This paper predicted nonlinear dynamic responses of a cantilevered carbon nanotube(CNT) resonator incorporating the electrostatic forces and van der Waals interactions between the CNT cantilever and ground plane. The structural model of CNT includes geometric and inertial nonlinearities to investigate various phenomena of nonlinear responses of the CNT due to the electrostatic excitation. In order to solve this problem, we used Galerkin's approximation and the numerical integration techniques. As a result, the CNT nano-resonator shows the softening effect through saddle-node bifurcation near primary resonance frequency with increasing the applied AC and DC voltages. Also we can predict nonlinear secondary resonances such as superharmonic and subharmonic resonances. The superharmonic resonance of the nano-resonator is influenced by applied AC voltage. The period-doubling bifurcation leads to the subharmonic resonance which occurs when the nano-resonator is actuated by electrostatic forces as parametric excitation.

Formulation of fully coupled THM behavior in unsaturated soil (불포화지반에 대한 열-수리-역학 거동의 수식화)

  • Shin, Ho-Sung
    • Proceedings of the Korean Geotechical Society Conference
    • /
    • 2010.09a
    • /
    • pp.808-812
    • /
    • 2010
  • A great deal of attention is focused on coupled Thermo-Hydro-Mechanical (THM) behavior of multiphase porous media 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 3 mass balances equations (solid, water, and air), energy balance equation, and force equilibrium equation. Finite element code is developed from the Galerkin formulation and time integration of these governing equations for 4 main variables (displacement $\underline{u}$, gas pressure $P_g$, liquid pressure $P_l$), and temperature T). The code is validated with theoretical solutions for linear material with simple boundary conditions.

  • PDF

Electro-Magnetic Field Computation Using the Natural Element Method (Natural Element Method를 이용한 전자장 해석)

  • Kim, H.K.;Jung, J.K.;Oh, Y.H.;Park, K.Y.
    • Proceedings of the KIEE Conference
    • /
    • 2003.04a
    • /
    • pp.59-61
    • /
    • 2003
  • The natural element method is a kind of meshless Galerkin method. The shape function is derived from the natural neighbor coordinates interpolation scheme. Natural neighbor shape functions are $C^0$ everywhere, except the nodes where they are $C^0$. The numerical integration is carried out using the Delaunay triangles as the background cells. The method is applied to the test problems and simulation results show that the natural element method can give accurate solutions for the electromagnetic field problems.

  • PDF