• Title/Summary/Keyword: the Newmark method

Search Result 347, Processing Time 0.022 seconds

Flow-induced Vibration Analysis for Cascades with Stator-rotor Interaction and Viscosity Effect (스테이터-로터 상호간섭 및 점성효과를 고려한 케스케이드의 유체유발 진동해석)

  • Oh, Se-Won;Park, Oung;Kim, Dong-Hyun
    • Transactions of the Korean Society for Noise and Vibration Engineering
    • /
    • v.16 no.10 s.115
    • /
    • pp.1082-1089
    • /
    • 2006
  • In this study, advanced computational analysis system has been developed in order to investigate flow-induced vibration(FIV) phenomenon for general stator-rotor cascade configurations. Relative movement of the rotor with respect to stator is reflected by modeling Independent two computational domains. Fluid domains are modeled using the unstructured grid system with dynamic moving and local deforming methods. Unsteady, Reynolds-averaged Wavier-stokes equations with one equation Spalart-Allmaras and two-equation SST ${\kappa}-{\varepsilon}$ turbulence models are solved for unsteady flow problems and also relative moving and vibration effects of the rotor cascade are fully considered. A coupled implicit time marching scheme based on the Newmark integration method is used for computing the governing equations of fluid-structure interaction problems. Detailed vibration responses for different flow conditions are presented and then vibration characteristics are physically investigated in the time domain as computational virtual tests.

Aerodynamic and Structural Design of 6kW Class Vertical-Axis Wind Turbine (공탄성 변형효과를 고려한 5MW급 풍력발전 블레이드의 피치각에 따른 성능해석)

  • Kim, Yo-Han;Kim, Dong-Hyun;Hwang, Mi-Hyun;Kim, Kyung-Hee;Hwang, Byung-Sun;Hong, Un-Sung
    • The KSFM Journal of Fluid Machinery
    • /
    • v.14 no.3
    • /
    • pp.39-44
    • /
    • 2011
  • In this study, performance analyses have been conducted for a 5MW class wind turbine blade model. Advanced computational analysis system based on computational fluid dynamics(CFD) and computational structural dynamics(CSD) has been developed in order to investigate detailed dynamic responsed of wind turbine blade. Reynolds-averaged Navier-Stokes (RANS) equations with K-${\epsilon}$ turbulence model are solved for unsteady flow problems of the rotating turbine blade model. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of the 3D turbine blade for fluid-structure interaction (FSI) problems. Predicted aerodynamic performance considering structural deformation effect of the blade show different results compared to the case of rigid blade model.

Inelastic displacement ratios for evaluation of stiffness degrading structures with soil structure interaction built on soft soil sites

  • Aydemir, Muberra Eser
    • Structural Engineering and Mechanics
    • /
    • v.45 no.6
    • /
    • pp.741-758
    • /
    • 2013
  • In this study, inelastic displacement ratios are investigated for existing systems with known lateral strength considering soil structure interaction. For this purpose, SDOF systems for period range of 0.1-3.0 s with different hysteretic behaviors are considered for a number of 18 earthquake motions recorded on soft soil. The effect of stiffness degradation on inelastic displacement ratios is investigated. The Modified Clough model is used to represent structures that exhibit significant stiffness degradation when subjected to reverse cyclic loading and the elastoplastic model is used to represent non-degrading structures. Soil structure interaction analyses are conducted by means of equivalent fixed base model effective period, effective damping and effective ductility values differing from fixed-base case. For inelastic time history analyses, Newmark method for step by step time integration was adapted in an in-house computer program. A new equation is proposed for inelastic displacement ratio of system with SSI with elastoplastic or degrading behavior as a function of structural period ($\tilde{T}$), strength reduction factor (R) and period lengthening ratio ($\tilde{T}$/T). The proposed equation for $\tilde{C}_R$ which takes the soil-structure interaction into account should be useful in estimating the inelastic deformation of existing structures with known lateral strength.

Transonic Flutter Characteristics of Supercritical Airfoils Considering Shockwave and Flow Separation Effects (충격파 및 유동박리 효과를 고려한 초임계 에어포일의 천음속 플러터 특성)

  • Kim, Dong-Hyun;Kim, Yu-Sung;Kim, Yo-Han;Kim, Seok-Soo
    • Proceedings of the Korean Society for Noise and Vibration Engineering Conference
    • /
    • 2008.11a
    • /
    • pp.167-174
    • /
    • 2008
  • In this study, flutter analyses for supercritical airfoil have been conducted in transonic region. Advanced computational analysis system based on computational fluid dynamics (CFD) and computational structural dynamics (CSD) has been developed in order to investigate detailed static and dynamic responses of supercritical airfoil. Reynolds-averaged Navier-Stokes equations with Spalart-Allmaras (S-A) and SST ${\kappa}-{\omega}$ turbulence models are solved for unsteady flow problems. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of cascades for fluid-structure interaction (FSI) problems. Also, flow-induced vibration (FIV) analyses for various supercritical airfoil models have been conducted. Detailed flutter responses for supercritical are presented to show the physical performance and vibration characteristics in various angle of attack.

  • PDF

Transonic Flutter Characteristics of Supercritical Airfoils Considering Shockwave and Flow Separation Effects (충격파 및 유동박리 효과를 고려한 초임계 에어포일의 천음속 플러터 특성)

  • Lin, Han;Kim, Dong-Hyun;Kim, Yu-Sung;Kim, Yo-Han;Kim, Seok-Soo
    • Journal of the Korean Society for Aviation and Aeronautics
    • /
    • v.17 no.2
    • /
    • pp.8-17
    • /
    • 2009
  • In this study, flutter analyses for supercritical airfoil have been conducted in transonic region. Advanced computational analysis system based on computational fluid dynamics (CFD) and computational structural dynamics (CSD) has been developed in order to investigate detailed static and dynamic responses of supercritical airfoil. Reynolds-averaged Navier-Stokes equations with Spalart-Allmaras (S-A) and SST ${\kappa}-{\omega}$ turbulence models are solved for unsteady flow problems. A fully implicit time marching scheme based on the Newmark direct integration method is used for computing the coupled aeroelastic governing equations of cascades for fluid-structure interaction (FSI) problems. Also, flow-induced vibration (FIV) analyses for various supercritical airfoil models have been conducted. Detailed flutter responses for supercritical are presented to show the physical performance and vibration characteristics in various angle of attack.

  • PDF

Earthquake response of nanocomposite concrete pipes conveying and immersing in fluid using numerical methods

  • Maleki, Mostafa;Bidgoli, Mahmood Rabani;Kolahchi, Reza
    • Computers and Concrete
    • /
    • v.24 no.2
    • /
    • pp.125-135
    • /
    • 2019
  • Concrete pipelines are the most efficient and safe means for gas and oil transportation over a long distance. The use of nano materials and nono-engineering can be considered for enhancing concrete pipelines properties. the tests show that SiO2 nanoparticles can improve the mechanical behavior of concrete. Moreover, severe hazard for pipelines is seismic ground motion. Over the years, scientists have attempted to understand pipe behavior against earthquake most frequently via numerical modeling and simulation. Therefore, in this paper, the dynamic response of underwater nanocomposite submerged pipeline conveying fluid is studied. The structure is subjected to the dynamic loads caused by earthquake and the governing equations of the system are derived using mathematical model via Classic shell theory and Hamilton's principle. Navier-Stokes equation is employed to calculate the force due to the fluid in the pipe. As well, the effect of external fluid is modeled with an external force. Mori-Tanaka approach is used to estimate the equivalent material properties of the nanocomposite. 1978 Tabas earthquake in Iran is considered for modelling seismic load. The dynamic displacement of the structure is extracted using differential quadrature method (DQM) and Newmark method. The effects of different parameters such as SiO2 nanoparticles volume percent, boundary conditions, thickness to radius ratios, length to radius ratios, internal and external fluid pressure and earthquake intensity are discussed on the seismic response of the structure. From results obtained in this paper, it can be found that the dynamic response of the pipe is increased in the presence of internal and external fluid. Furthermore, the use of SiO2 nanoparticles in concrete pipeline reduces the displacement of the structure during an earthquake.

Combined multi-predict-correct iterative method for interaction between pulsatile flow and large deformation structure

  • Wang, Wenquan;Zhang, Li-Xiang;Yan, Yan;Guo, Yakun
    • Coupled systems mechanics
    • /
    • v.1 no.4
    • /
    • pp.361-379
    • /
    • 2012
  • This paper presents a fully coupled three-dimensional solver for the analysis of interaction between pulsatile flow and large deformation structure. A partitioned time marching algorithm is employed for the solution of the time dependent coupled discretised problem, enabling the use of highly developed, robust and well-tested solvers for each field. Conservative transfer of information at the fluid-structure interface is combined with an effective multi-predict-correct iterative scheme to enable implicit coupling of the interacting fields at each time increment. The three-dimensional unsteady incompressible fluid is solved using a powerful implicit time stepping technique and an ALE formulation for moving boundaries with second-order time accurate is used. A full spectrum of total variational diminishing (TVD) schemes in unstructured grids is allowed implementation for the advection terms and finite element shape functions are used to evaluate the solution and its variation within mesh elements. A finite element dynamic analysis of the highly deformable structure is carried out with a numerical strategy combining the implicit Newmark time integration algorithm with a Newton-Raphson second-order optimisation method. The proposed model is used to predict the wave flow fields of a particular flow-induced vibrational phenomenon, and comparison of the numerical results with available experimental data validates the methodology and assesses its accuracy. Another test case about three-dimensional biomedical model with pulsatile inflow is presented to benchmark the algorithm and to demonstrate the potential applications of this method.

Transverse Vibration Analysis of the Deploying Beam by Simulation and Experiment (시뮬레이션과 실험을 통한 전개하는 보의 횡 방향 진동 분석)

  • Kim, Jaewon;Zhu, Kefei;Chung, Jintai
    • Transactions of the Korean Society for Noise and Vibration Engineering
    • /
    • v.25 no.12
    • /
    • pp.866-873
    • /
    • 2015
  • The transverse vibration of the deploying beam from rigid hub was analyzed by simulation and experiment. The linear governing equation of the deploying beam was obtained using the Euler-Bernoulli beam theory. To discretize the governing equation, the Galerkin method was used. After transforming the governing equation into the weak form, the weak form was discretized. The discretized equation was expressed by the matrix-vector form, and then the Newmark method was applied to simulate. To consider the damping effect of the beam, we conducted the modal test with various beam length. The mass proportional damping was selected by the relation of the first and second damping ratio. The proportional damping coefficient was calculated using the acquired natural frequency and damping ratio through the modal test. The experiment was set up to measure the transverse vibration of the deploying beam. The fixed beam at the carriage of the linear actuator was moved by moving the carriage. The transverse vibration of the deploying beam was observed by the Eulerian description near the hub. The deploying or retraction motion of the beam had the constant velocity and the velocity profile with acceleration and deceleration. We compared the transverse vibration results by the simulation and experiment. The observed response by the Eulerian description were analyzed.

An algorithm for quantifying dynamic buckling and post-buckling behavior of delaminated FRP plates with a rectangular hole stiffened by smart (SMA) stitches

  • Soltanieh, Ghazaleh;Yam, Michael C.H.
    • Smart Structures and Systems
    • /
    • v.28 no.6
    • /
    • pp.745-760
    • /
    • 2021
  • Dynamic buckling of structure is one of the failure modes that needs to be considered since it may result in catastrophic failure of the structure in a short period of time. For a thin fiber-reinforced polymer (FRP) plate under compression, buckling is an inherent hazard which will be intensified by the existence of defects like holes, cracks, and delamination. On the other hand, the growth of the delamination is another prime concern for thin FRP plates. In the current paper, reinforcing the plates against buckling is realized by using SMA wires in the form of stitches. A numerical framework is proposed to simulate the dynamic instability emphasizing the effect of the SMA stitches in suppressing delamination growth. The suggested algorithm is more accurate than the other methods when considering the transformation point of the SMA wires and the modeling of the cohesive zone using simple and yet reliable technique. The computational design of the method by producing the line by line orders leads to a simple algorithm for simulating the super-elastic behavior. The Lagoudas constitutive model of the SMA material is implemented in the form of user material subroutines (VUMAT). The normal bilinear spring model is used to reproduce the cohesive zone behavior. The nonlinear finite element formulation is programmed into FORTRAN using the Newmark-beta numerical time-integration approach. The obtained results are compared with the results obtained by the finite element method using ABAQUS/Explicit solver. The obtained results by the proposed algorithm and those by ABAQUS are in good agreement.

Auto-parametric resonance of framed structures under periodic excitations

  • Li, Yuchun;Gou, Hongliang;Zhang, Long;Chang, Chenyu
    • Structural Engineering and Mechanics
    • /
    • v.61 no.4
    • /
    • pp.497-510
    • /
    • 2017
  • A framed structure may be composed of two sub-structures, which are linked by a hinged joint. One sub-structure is the primary system and the other is the secondary system. The primary system, which is subjected to the periodic external load, can give rise to an auto-parametric resonance of the second system. Considering the geometric-stiffness effect produced by the axially internal force, the element equation of motion is derived by the extended Hamilton's principle. The element equations are then assembled into the global non-homogeneous Mathieu-Hill equations. The Newmark's method is introduced to solve the time-history responses of the non-homogeneous Mathieu-Hill equations. The energy-growth exponent/coefficient (EGE/EGC) and a finite-time Lyapunov exponent (FLE) are proposed for determining the auto-parametric instability boundaries of the structural system. The auto-parametric instabilities are numerically analyzed for the two frames. The influence of relative stiffness between the primary and secondary systems on the auto-parametric instability boundaries is investigated. A phenomenon of the "auto-parametric internal resonance" (the auto-parametric resonance of the second system induced by a normal resonance of the primary system) is predicted through the two numerical examples. The risk of auto-parametric internal resonance is emphasized. An auto-parametric resonance experiment of a ${\Gamma}$-shaped frame is conducted for verifying the theoretical predictions and present calculation method.