Chaos caused by fatigue crack growth

The nonlinear dynamic responses including chaotic oscillations caused by a fatigue crack growth are presented. Fatigue tests have been conducted on a novel fatigue-testing rig, where the loading is generated from inertial forces. The nonlinearity is in the form of discontinuous stiffness caused by the opening and closing of a growing crack. Nonlinear dynamic tools such as Poincaré maps and bifurcation diagrams are used to unveil the global dynamics of the system. The results obtained indicate that fatigue crack growth strongly influences the dynamic response of the system leading to chaos.     

Response of a tethered aerostat to simulated turbulence

Aerostats are lighter-than-air vehicles tethered to the ground by a cable and used for broadcasting, communications, surveillance, and drug interdiction. The dynamic response of tethered aerostats subject to extreme atmospheric turbulence often dictates survivability. This paper develops a theoretical model that predicts the planar response of a tethered aerostat subject to atmospheric turbulence and simulates the response to 1000 simulated hurricane scale turbulent time histories. The aerostat dynamic model assumes the aerostat hull to be a rigid body with non-linear fluid loading, instantaneous weathervaning for planar response, and a continuous tether. Galerkin’s method discretizes the coupled aerostat and tether partial differential equations to produce a non-linear initial value problem that is integrated numerically given initial conditions and wind inputs. The proper orthogonal decomposition theorem generates, based on Hurricane Georges wind data, turbulent time histories that possess the sequential behavior of actual turbulence, are spectrally accurate, and have non-Gaussian density functions. The generated turbulent time histories are simulated to predict the aerostat response to severe turbulence. The resulting probability distributions for the aerostat position, pitch angle, and confluence point tension predict the aerostat behavior in high gust environments. The dynamic results can be up to twice as large as a static analysis indicating the importance of dynamics in aerostat modeling. The results uncover a worst case wind input consisting of a two-pulse vertical gust.     

Dynamic solid mechanics using finite volume methods

A procedure for evaluating the dynamic structural response of elastic solid domains is presented. A prerequisite for the analysis of dynamic fluid–structure interaction is the use of a consistent set of finite volume (FV) methods on a single unstructured mesh. This paper describes a three-dimensional (3D) FV, vertex-based method for dynamic solid mechanics. A novel Newmark predictor–corrector implicit scheme was developed to provide time accurate solutions and the scheme was evaluated on a 3D cantilever problem. By employing a small amount of viscous damping, very accurate predictions of the fundamental natural frequency were obtained with respect to both the amplitude and period of oscillation. This scheme has been implemented into the multi-physics modelling software framework, Physica, for later application to full dynamic fluid structure interaction.     

Multidimensional turbulence spectra – Statistical analysis of turbulent vortices

Strong nonlinear or very fast phenomena such as mixing, coalescence and breakup in chemical engineering processes, are not correctly described using average turbulence properties. Since these phenomena are modeled by the interaction of fluid particles with single or paired vortices, distribution of the properties of individual turbulent vortices should be studied and understood. In this paper, statistical analysis of turbulent vortices was performed using a novel vortex tracking algorithm. The vortices were identified using the normalized Q-criterion with extended volumes calculated using the Biot–Savart law in order to capture most of the coherent structure related to each vortex. This new and fast algorithm makes it possible to estimate the volume of all resolved vortices. Turbulence was modeled using large-eddy simulation with the dynamic Smagorinsky–Lilly subgrid scale model for different Reynolds numbers. Number density of turbulent vortices were quantified and compared with different models. It is concluded that the calculated number densities for vortices in the inertial subrange and also for the larger scales are in very good agreement with the models proposed by Batchelor and Martinez-Bazán. Moreover, the associated enstrophy within the same size of coherent structures is quantified and its distribution is compared to models for distribution of turbulent kinetic energy. The associated enstrophy within the same size of coherent structures has a wide distribution that is normal distributed in the logarithmic scale.     

动载下缝端应力强度因子计算的扩展有限元法

在扩展有限元法(extended finite element methods, XFEM)的理论框架下,重点研究了动荷载作用下稳定裂纹尖端动态应力强度因子(dynamic stress intensity factors, DSIFs)的求解方法．根据XFEM的位移模式,推导了动力XFEM的支配方程,采用Newmark隐式算法进行时间积分同时,提出一种XFEM质量矩阵的集中策略,给出了求解DSIFs的相互作用积分方法,与静态问题的相互作用积分方法相比,增加了惯性项的贡献．最后,若干典型算例的计算结果表明:XFEM可以准确评价稳定裂纹尖端的DSIFs,建议的质量矩阵集中策略是有效的,为得到正确的DSIFs,惯性项的贡献不可忽略.     

带人工雨线的拉索在风激励下的响应

通过在风洞中对一具有可调动力特性、雨线位置和风向角的带人工雨线的拉索模型进行了试验,研究了拉索的风雨激振特性.试验结果和其他研究人员的结果进行了仔细对比,得到了一些新的结论,澄清了过去的一些模糊认识.结果表明,带人工雨线的水平索在风向为零时的响应可以用Den Hartog驰振机制来解释,而风向不为零时,拉索表现为限速振动或限速和驰振的混合型振动.     

Stability of a steady flow guided by flexible walls

Flexible channels conveying fluid may loose stability by divergence or flutter, thus leading to undesirable dynamic behaviour. A lot of investigations have already been done on this subject, each involving structural models expressing a specific type of stiffness. The aim of this contribution is to investigate and compare systematically the influence of various types of structural stiffness on the stability behaviour of the coupled problem. Therefore, a simple, yet general model is built. The fluid is bounded by a rigid and an elastical supported, flexible wall, which exhibits bending and extensional stiffness and considers effects of pre-load. The influence of the various characteristics on the stability of the steady state is discussed and a perturbation approach is used to investigate the influence of small nonlinearities. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The effect of the elasticity of an orbital tether system on the oscillations of a satellite

An orbital tether system, including a satellite (a rigid body), an elastic ponderable tether and a terminal load, is investigated. A mathematical model is obtained using Lagrange's equation of the second kind, which enables the plane translational motion of the centres of mass of the elements of the system and the rotational motion of the satellite and the tether to be investigated. It is shown that the equations of motion for the new independent variable, that is, the true anomaly angle, obtained on the assumption that the motion of the centre of mass of the system is independent of the relative motion of its elements, are an extension of the known mathematical models. The effect of the elasticity of the tether on the angular oscillations of the tether and the satellite is investigated. The model constructed can be used both to analyse of the deployment of a tether system as well as to investigate of the combined behaviour of a satellite and a tether about the natural centres of mass.     

CFD simulation of methane dispersion and innovative methane management in underground mining faces

This study addresses methane dispersion in a mine tunnel with discrete methane sources and various methods to handle it. Air flow behavior and methane dispersion in the tunnel are simulated utilizing the computational fluid dynamics (CFD) approach. Various possible conditions which may occur in a mine tunnel are investigated. Simulation results indicate that methane dispersion inside the mine tunnel is influenced significantly by the number as well as location of the sources and quantity of methane released from each source. Furthermore, application of an innovative flow divider which comprises volumetric flow control and flow director, is investigated. It is found that by properly directing the ventilation flow to the location where methane is accumulating can reduce methane concentration below the safe level. In addition, it is noted that focusing the ventilation flow at a point is more effective as compared to dispersing it at several points. This study provides some new ideas for designing an “intelligent” underground mine ventilation system which can cost-effectively maintain methane concentration below the critical value.     

Inertial Particle Motion in Recirculating Stokes Flow

Recirculation occurs in many cavity flows. In particular, alveolar flow models have been shown to exhibit recirculation patterns. However, many particles that are inhaled by the lungs do not follow this flow. Instead, they may diffuse into the surrounding flow or possess enough inertia to propel them from fluid particle paths. In this study, we construct a minimal model to observe the behavior of inertial particles caught within a recirculating Stokes flow. We find that, given favorable conditions, inertial particles can be cleared from the cavity or deposited on walls. This depends on the strength of inertia in zero gravity, but can be enhanced when gravity and the orientation of the cavity are taken into account. It is also possible for these effects to balance one another, producing a skewed limit cycle. These combined effects may play a significant part in the retention, deposition, and clearance of aerosols and particulates from alveolar cavities.     

Vibration analysis of truncated conical shells subjected to flowing fluid

A method to predict the dynamic behaviour of anisotropic truncated conical shells conveying fluid is presented in this paper. It is a combination of finite element method and classical shell theory. The displacement functions are derived from exact solutions of Sanders’ shell equilibrium equations of conical shells. The velocity potential, Bernoulli’s equation and impermeability condition have been applied to the shell–fluid interface to obtain an explicit expression for fluid pressure which yields three forces (inertial, centrifugal, Coriolis) of the moving fluid. To the best of the authors’ knowledge, this paper reports the first comparison made between two works which deal with conical shells subjected to internal flowing fluid effects. The results obtained by this method for conical shells with various boundary condition and geometries, in vacuum, fully-filled and when subjected to flowing fluid were compared with those of other experimental and numerical investigations and good agreement was obtained.     

Modeling and dynamics of a two-line kite

A mathematical model of a kite connected to the ground by two straight tethers of varying lengths is presented and used to study the traction force generated by kites flying in cross-wind conditions. The equations of motion are obtained by using a Lagrangian formulation, which yields a low-order system of ordinary differential equations free of constraint forces. Two parameters are chosen for the analysis. The first parameter is the wind velocity. The second parameter is one of the stability derivatives of the aerodynamic model: the roll response to the sideslip angle, known also as effective dihedral. This parameter affects significantly the lateral dynamics of the kite. It has been found that when the effective dihedral is below a certain threshold, the kite follows stable periodic trajectories, and naturally flies in cross-wind conditions while generating a high tension along both tethers. This result indicates that kite-based propulsion systems could operate without controlling tether lengths if kite design, including the dihedral and sweep angles, is done appropriately. If both tether lengths are varied out-of-phase and periodically, then kite dynamics can be very complex. The trajectories are chaotic and intermittent for values of the effective dihedral below a certain negative threshold. It is found that tether tensions can be very similar with and without tether length modulation if the parameters of the model are well-chosen. The use of the model for pure traction applications of kites is discussed.     

Conservation laws of turbulence models

The conservation of mass, momentum, energy, helicity, and enstrophy in fluid flow are important because these quantities organize a flow, and characterize change in the flow's structure over time. In turbulent flow, conservation laws remain important in the inertial range of wave numbers, where viscous effects are negligible. It is in the inertial range where energy, helicity (3d), and enstrophy (2d) must be accurately cascaded for a turbulence model to be qualitatively correct. A first and necessary step for an accurate cascade is conservation; however, many turbulent flow simulations are based on turbulence models whose conservation properties are little explored and might be very different from those of the Navier-Stokes equations.We explore conservation laws and approximate conservation laws satisfied by LES turbulence models. For the Leray, Leray deconvolution, Bardina, and Nth order deconvolution models, we give exact or approximate laws for a model mass, momentum, energy, enstrophy and helicity. The possibility of cascades for model quantities is also discussed.     

Numerical simulation of scour around a submarine pipeline using computational fluid dynamics and discrete element method

Scour under a submarine pipeline can lead to structural failure; hence, a good understanding of the scour mechanism is paramount. Various numerical methods have been proposed to simulate scour, such as potential flow theory and single-phase and two-phase turbulent models. However, these numerical methods have limitations such as their reliance on calibrated empirical parameters and inability to provide detailed information. This paper investigates the use of a coupled computational fluid dynamics-discrete element method (CFD-DEM) model to simulate scour around a pipeline. The novelty of this work is to use CFD-DEM to extract detailed information, leading to new findings that enhance the current understanding of the underlying mechanisms of the scour process. The simulated scour evolution and bed profile are found to be in good agreement with published experimental results. Detailed results include the contours of the fluid velocity and fluid pressure, particle motion and velocity, fluid forces on the particles, and inter-particle forces. The sediment transport rate is calculated using the velocity of each single particle. The quantitative analysis of the bed load layer is also presented. The numerical results reveal three scour stages: onset of scour, tunnel erosion, and lee-wake erosion. Particle velocity and force distributions show that during the tunnel erosion stage, the particle motion and particle–particle interactive forces are particularly intense, suggesting that single-phase models, which are unable to account for inter-particle interactions, may be inadequate. The fluid pressure contours show a distinct pressure gradient. The pressure gradient force is calculated and found to be comparable with the drag force for the onset of scour and the tunnel erosion. However, for the lee-wake erosion, the drag force is shown to be the dominant mechanism for particle movements.     

输送流体管道的固——液耦合动力学研究

根据Ｈａｍｉｌｔｏｎ原理推导输送流体管道固—液耦合振动方程，得到反对称的固—液耦合阻尼矩阵和对称的固—液耦合刚度矩阵；用ＱＲ法计算管道固有频率，给出了管道前４阶固有频率—流速曲线；讨论了流体的流速、压强变化以及固—液耦合阻尼和固—液耦合刚度对管道固有频率的影响；用Ｎｅｗｍａｒｋ法计算不同流速时管道对阶跃载荷的动力响应；发现了各阶固有频率都有随流速的提高而降低、再提高、再降低的周而复始现象·     

Steady-state response of a geared rotor system with slant cracked shaft and time-varying mesh stiffness

The dynamic behavior of geared rotor system with defects is helpful for the failure diagnosis and state detecting of the system. Extensive efforts have been devoted to study the dynamic behaviors of geared systems with tooth root cracks. When surface cracks (especially for slant cracks) appear on the transmission shaft, the dynamic characteristics of the system have not gained sufficient attentions. Due to the parametric excitations induced by slant crack breathing and time-varying mesh stiffness, the steady-state response of the cracked geared rotor system differs distinctly from that of the uncracked system. Thus, utilizing the direct spectral method (DSM), the forced response spectra of a geared rotor system with slant cracked shaft and time-varying mesh stiffness under transmission error, unbalance force and torsional excitations are, respectively, obtained and discussed in detail. The effects of crack types (straight or slant crack) and crack depth on the forced response spectra of the system without and with torsional excitation are considered in the analysis. In addition, how the frequency response characteristics change after considering the crack is also investigated. It is shown that the torsional excitations have significant influence on the forced response spectra of slant cracked system. Sub-critical resonances are also found in the frequency response curves. The results could be used for shaft crack detection in geared rotor system.     

Mathieu equation with application to analysis of dynamic characteristics of resonant inertial sensors

In this paper, Mathieu equation is applied to analyze the dynamic characteristics of resonant inertial sensors. Unlike previous work, Mathieu equation is not just a differential equation and analyzes the stability of the transition curves, but become an important method in analyzing parametric resonant characteristics and approximate output of resonant inertial sensors. It is demonstrated that the mathematical model of resonant inertial sensors is described by Mathieu equation. The relevant Mathieu equation theory and dynamic characteristics analysis methods were proposed, which include both stability and dynamic linear output. Finally, theoretical and experimental analysis show that the correlation of the theoretical curve and the experimental result coincide so perfectly, which means proposed analysis methods for Mathieu equation could be used to analyze the dynamic output characteristic of resonant inertial sensors. The theoretical analyzing approach of Mathieu equation and experimental results of resonant inertial sensors are obtained, which provide an application area for Mathieu equation and a reference for the robust design for resonant inertial sensors.     

Magnetohydrodynamics flow of a Burgers’ fluid in an orthogonal rheometer

The magnetohydrodynamics flow of an electrically conducting, incompressible Burgers’ fluid in an orthogonal rheometer is investigated. An exact solution is obtained. The effects of various dimensionless parameters existing in the model on the velocity field, vorticity and traction are studied graphically. It is noted that boundary layers form for a variety of reasons. It form as the Reynolds number increases. Also, as the Weissenberg number increases a distinct boundary layer formation is observed. It can develop at low Reynolds number provided the Weissenberg number is sufficiently high, however, it is not possible in the case of a Newtonian fluid. It is shown that no torque is exerted by the fluid on one of the disks. Results are compared with Oldroyd-B fluid.     

计入气穴的水下爆炸作用下结构流固耦合三维数值模拟

水下爆炸在结构物面附近产生的气穴现象,严重影响水下爆炸作用下的流固耦合动响应,是舰船水下爆炸领域的难点,传统的边界元方法、有限元方法(FEM)难以解决水下爆炸气穴现象这类强非线性问题．针对此问题,计及流体中的气穴现象,考虑流体的可压缩型,忽略流体粘性,建立水下爆炸瞬态强非线性流固耦合三维数值模型,采用流体谱单元方法(SEM)和结构有限元方法求解该模型．计算结果表明:相对有限元法,谱单元法具有更高的计算精度,且谱单元解与解析解、试验值吻合良好．在此基础上,结合ABAQUS软件,分别探讨三维球壳、船体板架在水下爆炸作用下的瞬态流固耦合机理,给出气穴区域及其对水中结构物动响应的影响特征,旨在为舰船水下爆炸瞬态流固耦合问题的相关研究提供参考．     

Pitfalls in the frequency response represented onto Polynomial Chaos for random SDOF mechanical systems

Uncertainties are present in the modeling of dynamical systems and they must be taken into account to improve the prediction of the models. It is very important to understand how they propagate and how random systems behave. This study aims at pointing out the somehow complex behavior of the structural response of stochastic dynamical systems and consequently the difficulty to represent this behavior using spectral approaches. The main objective is to find numerically the probability density function (PDF) of the response of a random linear mechanical systems. Since it is found that difficulties can occur even for a single-degree-of-freedom system when only the stiffness is random, this work focuses on this application to test several methods. Polynomial Chaos performance is first investigated for the propagation of uncertainties in several situations of stiffness variances for a damped single-degree-of-freedom system. For some specific conditions of damping and stiffness variances, it is found that numerical difficulties occur for the standard polynomial bases near the resonant frequency, where it is generally observed that the shape of the system response PDFs presents multimodality. Strategies to build enhanced bases are then proposed and investigated with varying degrees of success. Finally, a multi-element approach is used in order to gain robustness.