Model experiments of the transient response to flooding of the box shaped barge

Coupling of the flooded water and ship motions was studied experimentally. Roll decay tests for one flooded compartment and transient abrupt flooding tests were performed for the box shaped barge model. The tests were conducted to obtain information on the flooding process for the development of numerical tools and to provide validation data. Quantitative values on the effect of flooded water on the roll damping were obtained. Flooded water behaves in a different manner in undivided and divided compartments. Flooded water in divided compartment increases roll damping significantly. In undivided compartment roll damping was high at low amount of flooded water. For higher amounts damping was of the same order as for the intact model. Initial flooding is a complex process where the ship and flooded water motions are coupled. Propagation of the flooding water inside the compartment, at a dam-break type abrupt flooding, was studied by tracking the surface of the flooded water. An image processing algorithm was used to obtain the tracked surface. Flooded water volume and its center of gravity were estimated from the tracked surface. Different internal layouts of the flooded compartment can lead to a totally different roll response. The inflooding jet plays an important role in the response in case of the undivided compartment. While, for a divided compartment, asymmetric flooding due to the obstructions causes high heel angle on the damage side.     

Evolution of governing mass and momentum balances following an abrupt pressure impact in a porous medium

A mathematical model is developed of an abrupt pressure impact applied to a compressible fluid flowing through a porous medium domain. Nondimensional forms of the macroscopic fluid mass and momentum balance equations yield two new scalar numbers relating storage change to pressure rise. A sequence of four reduced forms of mass and momentum balance equations are shown to be associated with a sequence of four time periods following the onset of a pressure change. At the very first time period, pressure is proven to be distributed uniformly within the affected domain. During the second time interval, the momentum balance equation conforms to a wave form. The behavior during the third time period is governed by the averaged Navier-Stokes equation. After a long time, the fourth period is dominated by a momentum balance similar to Brinkman's equation which may convert to Darcy's equation when friction at the solid-fluid interface dominates.     

The effect of rigid-body motions on the whipping response of a ship hull subjected to an underwater bubble

The dynamic elastic response of a floating ship hull girder to an underwater explosion bubble is normally composed of two parts: rigid-body motion and elastic deformation. However, the effects of rigid-body motion have consistently been neglected in the current literature based on the assumption that they are small. In this paper, our focus is on the study of rigid-body motion effects on the hull girder's elastic deformation, also known as the ‘whipping response’. A theory of interaction between a gas bubble and a hull girder is presented. The bubble dynamic equations combined with the bubble migration, free surface effect and drag force considerations are solved numerically using the Runge–Kutta method. The rigid-body and elastic responses of the hull that are induced by the impulsive pressure of a bubble are calculated using the methods presented herein. Two different examples of real ships are given to demonstrate the effect of rigid-body motion on whipping responses. The numerical results show that rigid-body motions reduce the amplitudes and vibration natural periods of the bending moments of the hull girder. These effects can be ignored for slender hulls, but must be taken into consideration for shorter/wider hulls so as not to underestimate the longitudinal strength.     

水中爆炸冲击波载荷作用下舰船结构动态响应的数值模拟

采用大型通用有限元软件MSC.DYTRAN,对某型水面舰船全船结构在水下爆炸冲击波载荷作用下的动态响应进行了数值模拟。分别从结构的变形损伤形式、能量吸收和冲击环境等几个方面研究了舰船结构在水下爆炸载荷作用下的破坏机理和响应特征。     

On the convergence and causality of a frequency domain method for dynamic structural analysis

Venanico-Filho et al. developed an elegant matrix formulation for dynamic analysis by frequency domain (FD), but the convergence, causality and extended period need further refining. In the present paper, it was argued that: (1) under reasonable assumptions (approximating the frequency response function by the discrete Fourier transform of the discretized unitary impulse response function), the matrix formulation by FD is equivalent to a circular convolution; (2) to avoid the wraparound interference, the excitation vector and impulse response must be padded with enough zeros; (3) provided that the zero padding requirement satisfied, the convergence and accuracy of direct time domain analysis, which is equivalent to that by FD, are guaranteed by the numerical integration scheme; (4) the imaginary part of the computational response approaching zero is due to the continuity of the impulse response functions. The English text was polished by Yunming Chen     

On the stability of the solution to a gonorrhea discrete mathematical model

ＯＮＴＨＥＳＴＡＢＩＬＩＴＹＯＦＴＨＥＳＯＬＵＴＩＯＮＴＯＡＧＯＮＯＲＲＨＥＡＤＩＳＣＲＥＴＥＭＡＴＨＥＭＡＴＩＣＡＬＭＯＤＥＬＪｉｎＪｕｎ（金均）（ＳｈａｎｇｈａｉＴｅａｃｈｅｒｓＵｎｉｖｅｒｓｉｔｙ,Ｓｈａｎｇｈａｉ）（ＲｅｃｅｉｖｅｄＭａｒｃｈ...     

Disturbances to a rotating fluid with a sphere moving along the axis of rotation

I.IntroductionWhenabodyprojectedinthewaterisintranslationalmotionthroughwateritwillcertainlycausesdisturbances.Howtodetectthedisturbancesf'arawayfromthebodyandhowtodetermillebytheintbrlnationobtainedthepositionandvelocityofthebodyandthesiteofprojectionareimportant.Astheprojectedbodyisrotating,thefluidaroundisinauniformlyrotationrelativetothereferenceframerotatinginsynchronismwiththebody,andthisisaquestionofthemotionofabodyinarotatingfluid.SuchproblemswerefirststudiedbyProudman(1916)l']andTay…     

Transient response of a finite crack in an orthotropic strip under the in-plane shear impact

The transient response of a central crack in an orthotropic strip under the in-plane shear impact loading is studied by using the dual integral equation method proposed by Copson and Sih. The general formula for the shear stress intensity factor near the crack tip is derived. Numerical results of with in various cases are obtained by solving the second kind Fredholm integral equation and by performing the inverse Laplace transform.     

Boundary Integral Analysis of the Symmetric Dynamic Problem for an Infinite Bimaterial Solid with an Embedded Crack

The symmetric frequency domain problem for two ideally bonded elastic half-spaces with a perpendicular plane crack is considered. It is reduced to the boundary integral equation (BIE) with integration over the limited crack region. The contact conditions on the bimaterial interface are satisfied identically in the initial stage of obtaining the equation. After boundary element solution of the equation, the stress concentration in the vicinity of a penny-shaped crack under time-harmonic loading of constant amplitude is studied. The mode I stress intensity factors as functions of angular coordinate of a crack front point and wave number for various relations between the material parameters are computed. The crack depth relative to the bimaterial interface is determined, when the effect of the material dissimilarity on the crack can be neglected.     

On the distributed Lagrange multiplier/fictitious domain method for rigid‐particle‐laden flows: a proposition for an alternative formulation of the Lagrange multipliers

The distributed Lagrange multiplier/fictitious domain method proposed for the direct numerical simulation of particle‐laden flows is considered in this work. First, it is demonstrated that improved accuracy is obtained with a coupled numerical scheme, whereby the pressure and the Lagrange multiplier fields enforcing incompressibility and rigid body motion, respectively, are calculated and applied together. However, the convergence characteristics of the iterative solution of the coupled scheme are poor because symmetric but indefinite and poorly conditioned matrices are produced. An analysis is then presented, which suggests that the cause for the matrix pathologies lies in the interaction of the respective matrix operators enforcing incompressibility and rigid body motion. On the basis of this analysis, an alternative formulation is developed for the Lagrange multipliers, being now composed of a set of forces distributed only on the particle boundary together with a set of couples distributed within the particle core. The new formulation is tested with several types of flows with stationary or moving particles under creeping or finite Reynolds number conditions and it is demonstrated that it produces correct results and better conditioned matrices, thus enabling faster and more reliable convergence of the conjugate gradient method. The analysis and tests, therefore, support the expectation that the proposed formulation is promising and worthy of further study and improvement. Copyright © 2011 John Wiley & Sons, Ltd.     

Measurement modes of the response time of a magneto-rheological fluid (MRF) for changing magnetic flux density

The response of a magneto-rheological fluid (MRF) to a change of magnetic flux density is investigated by using a commercial plate–plate magneto-rheometer MCR501 (Anton Paar GmbH) at constant shear rate. The instrument was modified to allow an online determination of the transient flux density in the MRF. Both current and voltage imposition to the magneto-cell were applied by using a power operational amplifier to drive the electromagnet. Assuming a Maxwell behavior with switching time λ and a linear increase in shear stress with flux density, analytic relations for the transient shear stress are derived for sinusoidal and single exponential flux densities vs time. True switching times of a few milliseconds are only obtained if the low pass filter in the original MCR501 torque signal is surpassed by a firmware allowing a sampling rate of 0.1 ms. For a sinusoidal flux density, the switching time is derived from the modulation depth of the shear stress. An upper bound of λ < 3 ms for a flux density of 0.8 T was found. For step coil current imposition of 1 T magnitude, switching times of 2.8 ms (start-up) and 1.8 ms (shutdown) allowed to fit the transient torque signal more than 2/3 of the total change. Finally, the effect of a sigmoidal characteristic on the switching time determination is addressed. This paper was presented at Annual European Rheology Conference (AERC) held in Hersonisos, Crete, Greece, April 27-29, 2006.     

Application of Nonlinear Normal Mode Analysis to the Nonlinear and Coupled Dynamics of a Floating Offshore Platform with Damping

The nonlinear dynamics of ships and floating offshore platforms hasattracted much attention over the last several years. However the topicof multiple-degrees-of-freedom systems has almost been completely ignoredwith very few exceptions. This is probably due to the complexity ofanalyzing strongly nonlinear and coupled systems. It turns out thatcoupling may be particularly important for certain critical dynamicssuch as the dynamics of a floating offshore platform about its diagonalaxis. In a previous work, Kota et al. [1] applied the recently developed nonlinearnormal mode technique to analyze the coupled nonlinear dynamics of afloating offshore platform. Although this previous work was restrictedto unforced and undamped systems, in this work a comparison of the twoalternative nonlinear normal mode analysis techniques was completed.Considering the relative practical importance of damping versus externalforcing for this system, in the present work, we utilize just one of thetwo major techniques available [2] to analyze damped multiple-degrees-of-freedom nonlinear dynamics. Specifically, we investigate the effect ofnonlinearity, and non-proportionate damping. Our results show that thistechnique allows one to simply consider the effect of nonlinearity andgeneral damping on the resulting normal modes. This technique isparticularly powerful because it allows one to visualize the modes in ageometric fashion using the invariant manifold concept from dynamicalsystems.     

Effect of the viscosity properties of ICE on the deflection of an ICE sheet subjected to a moving load

Steady-state rectilinear motion of a load on an ice sheet modeled by a viscoelastic plate is considered. The viscoelastic properties of ice are described using the linear Maxwell and Kelvin-Voigt models and a generalized Maxwell-Kelvin model. Calculated vertical displacements and strains of the ice plate are compared with available experimental data. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 3, pp. 147–157, May–June, 2009.     

纳米铜自修复性能评价标准研究

为了建立纳米抗磨添加剂自修复性能的评价标准,提出了以磨损试验后试块质量增重为衡量依据的实验设想。以商品纳米Cu抗磨自修复剂作为考察对象,在高速环块磨损试验机上,发现了在不同试验力、转速、试验时间和添加剂浓度等试验条件下,试块摩擦增重的变化规律,从而提出了对添加剂自修复性能进行评定的试验方法和条件,并进行了重复性和再现性测试。通过扫描电子显微镜、X射线能量色散谱仪等对磨痕表面形貌和元素的分析,证实了添加剂的抗磨自修复性能和评价方法的可行性。     

Symmetric-Galerkin boundary element analysis of the dynamic T-stress for the interaction of a crack with an auxetic inclusion

In this work, the symmetric-Galerkin boundary element method (SGBEM) for 2-D elastodynamics in the Laplace-space frequency domain (Laplace domain) is employed to study the dynamic stress intensity factors (DSIFs) and the dynamic T-stress (DTS) during the interaction between a crack and an auxetic inclusion under impact loading conditions. It is found that, while the auxeticity has virtually no effect on the DSIFs, its influence on the DTS is noticeable. This finding is particularly important as it implies the imperative need of fracture criteria based on both the DSIFs and DTS for predicting crack growth in composite materials with auxetic phases.     

变质量、变阻尼、变刚度结构系统的动力响应

本文首次提出了变质量、变阻尼、变刚度结构系统的动力响应问题,并根据抓斗卸船机的工作过程导出了结构系统的动力方程及解法。最后给出了工程实际算例。     

Analytical study of the dynamic response of an embedded railway track to a moving load

Summary In this paper, the steady-state dynamic response of an embedded railway track to a moving train is investigated theoretically. The model for the track consists of a flexible plate performing vertical vibrations, two beams that are connected to the plate by continuous visco-elastic elements and an elastic foundation that supports the plate. Two harmonic loads that move uniformly along the beams describe the train load. The plate, the beams and the elastic foundation are employed to model a concrete slab of an embedded track, the rails and the ground reaction, respectively. The problem is studied by employing the Fourier integral transforms in the following way. Firstly, the dispersion analysis of waves that may propagate along the system is accomplished in the frequency-wavenumber domain. On the basis of this analysis, critical velocities of the loads are found both for the in-phase and anti-phase vibrations of the loads. Secondly, the vertical displacement of the rails and the slab, along with the stresses in the slab, are investigated as functions of the velocity and frequency of the loads. Finally, the response of the two-dimensional model is compared to that of a simplified one-dimensional model.     

Stability of an axisymmetric swirled flow in a supersonic cocurrent stream for volume energy supply in the viscous vortex core

The results of numerical calculations of the stability of axisymmetric swirled flows in a viscous vortex embedded in a supersonic cocurrent stream with a constant circulation of the azimuthal velocity component are presented. The stability characteristics of the swirled three-dimensional viscous flow in the streamwise vortex are determined on the basis of the linearized system of Navier-Stokes equations for a viscous heat-conducting gas under the assumption that the basic undisturbed flow is locally plane-parallel. The disturbed flow stability is studied in the temporal formulation with respect to both symmetric and asymmetric three-dimensional waves traveling along the vortex axis and corresponding to both positive and negative values of the azimuthal wavenumber. It is shown that at external inviscid flow Mach numbers M = 2 and 3 thermal energy supply in a small region near the vortex axis leads to considerable restructuring of the basic undisturbed flow in the vicinity of the vortex core, the growth of the adverse pressure gradient along the vortex axis, and a significant change in the small perturbation stability and behavior.__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, 2005, pp. 71–80. Original Russian Text Copyright © 2005 by Kazakov.     

A fracture damage model for Jointed rock masses and its application to the stability analysis of dam abutments

Based on continuum damage mechanics, for jointed rock masses, a fracture damage model is presented in this paper. First, the damage tensors are defined through the elastic-flexibility of intact rock and the equivalent elastic-damage flexibility for rock mass. Then, by the self-consistent principle of solid mechanics, the equivalent elastic-damage flexibility tensors involving the interaction between multicracks are deduced. The damage evolution law is proposed involving the mechanism of crack propagation process: frictional sliding, crack kinking, growing of branched tension cracks, interlinking of the microcracks near branched crack tips leading to the breakthrough of macro-cracks and finally the failure of rock mass. Thus the evolution of damage variables reasonably unified with the process of crack propagation is given. Finally, a plastic-brittle damage constitutive relation including brittle coupled strain rate, developed and applied to the stability analysis of complicated rock foundation of a dam in China, is described in this paper.     

An iterative method for the computation of the response of linearised Navier–Stokes equations to harmonic forcing and application to forced cylinder wakes

The central aim of this paper is the development and application of an efficient, iterative methodology for the computation of the perturbation fields induced by harmonic forcing of the linearised Navier–Stokes equations. The problem is formulated directly in the frequency domain, and the resulting system of equations is solved iteratively until convergence. The method is easily implemented to any implicit code that can solve iteratively the steady‐state Navier–Stokes equations. In this paper, it is applied to investigate the flow around a static cylinder with pulsating approaching flow and a cylinder undergoing forced stream‐wise oscillations. All terms of the perturbation kinetic energy equation are computed, and it is shown that perturbations grow by extracting energy from two sources: the underlying base flow field and the externally provided energy that maintains the imposed oscillation. The periodic drag force acting on the cylinder is also computed, and it is demonstrated that Morrison's equation is a simple model that can estimate with good accuracy the amplitude and phase of this force with respect to the approaching flow. The perturbation fields induced by periodic inlet flow (static cylinder) and forced stream‐wise cylinder oscillation are closely related: the velocity fields are identical in the appropriate reference frames, and a simple expression is derived, which links the pressures in the two flow cases. Copyright © 2014 John Wiley & Sons, Ltd.