• Interaction and multiscale mechanics
• /
• v.1 no.1
• /
• pp.1-31
• /
• 2008

The present paper presents multiscale modelling via coupling of the discrete and finite element methods. Theoretical formulation of the discrete element method using spherical or cylindrical particles has been briefly reviewed. Basic equations of the finite element method using the explicit time integration have been given. The micr-macro transition for the discrete element method has been discussed. Theoretical formulations for macroscopic stress and strain tensors have been given. Determination of macroscopic constitutive properties using dimensionless micro-macro relationships has been proposed. The formulation of the multiscale DEM/FEM model employing the DEM and FEM in different subdomains of the same body has been presented. The coupling allows the use of partially overlapping DEM and FEM subdomains. The overlap zone in the two coupling algorithms is introduced in order to provide a smooth transition from one discretization method to the other. Coupling between the DEM and FEM subdomains is provided by additional kinematic constraints imposed by means of either the Lagrange multipliers or penalty function method. The coupled DEM/FEM formulation has been implemented in the authors' own numerical program. Good performance of the numerical algorithms has been demonstrated in a number of examples.

• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2005.04a
• /
• pp.295-300
• /
• 2005

The freeplay, one of the concentrated structural nonlinearities, is inevitable for control surfaces of a real air vehicle due to normal wear of components and manufacturing mismatches. Also aerodynamic nonlinearities caused by a shock wave occur in transonic region. In practice, these nonlinearities induce the limit cycle oscillation (LCO) and decrease the transonic flutter speed. In this study, the fictitious mass method is used to apply a modal approach to nonlinear structural models due to freeplay. The transonic small-disturbance (TSD) equation is used to calculate unsteady aerodynamic forces in transonic region. Nonlinear aeroelastic time responses are predicted by the coupled time integration method (CTIM). This method was also applied to a 3D all-movable control wing to investigate its nonlinear aeroelastic responses. The angle of attack effect on the LCO characteristics has been found to be closely related with the initial pitching moment.

• 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.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2001.05a
• /
• pp.675-680
• /
• 2001

In this study, a fluid/structure coupled analysis system using computational fluid dynamics and computational structural dynamics has been developed. The unsteady flow fields are predicted using unstructured Euler code. Coupled time-integration method (CTIM) was applied to computer simulation of the flow-induced vibration phenomena. To investigate the interaction effect of shock motions, 2-DOF airfoil systems have been studied in the subsonic and transonic flow region. Also, aeroelastic analyses for the airfoil with an arbitrary object are performed to show the analysis capability and interference effects for the complex geometries. The present results show the flutter stabilities and characteristics of aeroelastic responses with moving shock effects.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2007.05a
• /
• pp.275-282
• /
• 2007

In this study, a fluid-structure interaction (FSI) analysis system has been developed in order to evaluate the turbine cascade performance with blade structural deformation effect. Relative movement of the rotor with respect to stator is reflected by modeling independent two computational domains. To consider the deformed position of rotor airfoil, dynamic moving grid method is applied. Reynolds-averaged Navier-Stokes equations with one equation Spalart-Allmaras and two-equation SST $k-{\varepsilon}$ turbulence models are solved to predict unsteady fluid dynamic loads. A fully implicit time marching scheme based on the Newmark direct integration method with high artificial damping is used to compute the fluid-structure interaction problem. Cascade performance evaluations for different elastic axis positions are presented and compared each other. It is importantly shown that the predicted aerodynamic performance considering structural deformation effect of blade can show some deviations compared to the data generally computed from rigid blade configurations and the position of elastic axis also tend to give sensitive effect.

• Ocean Systems Engineering
• /
• v.4 no.2
• /
• pp.83-97
• /
• 2014

Wave propagation in a three-dimensional (3D) fully nonlinear numerical wave tank (NWT) is studied based on velocity potential theory. The governing Laplace equation with fully nonlinear boundary conditions on the moving free surface is solved using the indirect desingularized boundary integral equation method (DBIEM). The fourth-order predictor-corrector Adams-Bashforth-Moulton scheme (ABM4) and mixed Eulerian-Lagrangian (MEL) method are used for the time-stepping integration of the free surface boundary conditions. A smoothing algorithm, B-spline, is applied to eliminate the possible saw-tooth instabilities. The artificial wave speed employed in MTF (multi-transmitting formula) approach is investigated for fully nonlinear wave problem. The numerical results from incorporating the damping zone (DZ), MTF and MTF coupled DZ (MTF+DZ) methods as radiation condition are compared with analytical solution. An effective MTF+DZ method is finally adopted to simulate the 3D linear wave, second-order wave and irregular wave propagation. It is shown that the MTF+DZ method can be used for simulating fully nonlinear wave propagation very efficiently.

• Coupled systems mechanics
• /
• v.9 no.2
• /
• pp.163-182
• /
• 2020

Coupled thermo-mechanical analysis of reinforced concrete slab at elevated temperatures from a fire accounting for nonlinear thermal parameters is carried out. The main focus of the paper is put on a one-way continuous reinforced concrete slab exposed to fire from the single (bottom) side as the most typical working condition under fire loading. Although contemporary techniques alongside the fire protection measures are in constant development, in most cases it is not possible to avoid the material deterioration particularly nearby the exposed surface from a fire. Thereby the structural fire resistance of reinforced concrete slabs is mostly influenced by a relative distance between reinforcement and the exposed surface. A parametric study with variable concrete cover ranging from 15 mm to 35 mm is performed. As the first part of a one-way coupled thermo-mechanical analysis, transient nonlinear heat transfer analysis is performed by applying the net heat flux on the exposed surface. The solution of proposed heat analysis is obtained at certain time steps of interest by α-method using the explicit Euler time-integration scheme. Spatial discretization is done by the finite element method using a 1D 2-noded truss element with the temperature nodal values as unknowns. The obtained results in terms of temperature field inside the element are compared with available numerical and experimental results. A high level of agreement can be observed, implying the proposed model capable of describing the temperature field during a fire. Accompanying thermal analysis, mechanical analysis is performed in two ways. Firstly, using the guidelines given in Eurocode 2 - Part 1-2 resulting in the fire resistance rating for the aforementioned concrete cover values. The second way is a fully numerical coupled analysis carried out in general-purpose finite element software DIANA FEA. Both approaches indicate structural fire behavior similar to those observed in large-scale fire tests.

• Journal of computational fluids engineering
• /
• v.13 no.4
• /
• pp.33-38
• /
• 2008

For study on the unsteady blockage effect, flows around a rotationally oscillating circular cylinder with relatively low forcing frequency in closed test-section wind tunnels have been numerically investigated by solving compressible Navier-Stokes equations. The numerical scheme is based on a node-based finite-volume method with the Roe's flux-difference splitting and an implicit time-integration method coupled with dual time-step sub-iteration. The computed results of the oscillating cylinder in the test section showed that the fluctuations of lift and drag are augmented by the blockage effects. The drag further increases because of low base pressure. The pressure on the test section wall shows the harmonics having the oscillating and the shedding frequencies contained in the blockage effect.

• 한국전산유체공학회:학술대회논문집
• /
• 2000.05a
• /
• pp.66-72
• /
• 2000

Aerodynamic sensitivity analysis is performed for the Navier-Stokes equations coupled with two-equation turbulence models using a discrete adjoint method and a direct differentiation method respectively. Like the mean flow equations, the turbulence model equations are also hand-differentiated to accurately calculate the sensitivity derivatives of flow quantities with respect to design variables in turbulent viscous flows. Both the direct differentiation code and the adjoint variable code adopt the same time integration scheme with the flow solver to efficiently solve the differentiated equations. The sensitivity codes are then compared with the flow solver in terms of solution accuracy, computing time and computer memory requirements. The sensitivity derivatives obtained from the sensitivity codes with different turbulence models are compared with each other. Using two-equation turbulence models, it is observed that a usual assumption of constant turbulent eddy viscosity in adjoint methods may lead to seriously inaccurate results in highly turbulent flows.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.26 no.12
• /
• pp.2511-2518
• /
• 2002

Dynamic behaviors and stability of an optical disk drive coupled with an automatic ball balancer (ABB) are analyzed by a theoretical approach. The feeding system is modeled a rigid body with six degree-of-freedom. Using Lagrange's equation, we derive the nonlinear equations of motion for a non -autonomous system with respect to the rectangular coordinate. To investigate the dynamic stability of the system in the neighborhood of the equilibrium positions, the monodromy matrix technique is applied to the perturbed equations. On the other hand, time responses are computed by the Runge -Kutta method. We also investigate the effects of the damping coefficient and the position of ABB on the dynamic behaviors of the system.