Nonlinear dynamic response and active vibration control of the viscoelastic cable with small sag

IntroductionCablesareveryefficientstructuralmembersandhencehavebeenwidelyusedinmanylong_spanstructures,includingcable_supportbridges,guyedtowersandcable_supportroofs.Sincecablesarelight,veryflexibleandlightlydamped ,structuresutilizingcables,i.e .,cable_structuresystems,usuallyhavevariousdynamicproblems.Theirmodelsarethereforeverimportantinpredictingandcontrollingtheirresponses.Inthelastdecade,thenonlineardynamicvibrationandstabilitybehaviorofcablesandcable_structureshavedrawntheattentionofman…     

Homogenization of Wall-Slip Gas Flow Through Porous Media

The permeability of reservoir rocks is most commonly measured with an atmospheric gas. Permeability is greater for a gas than for a liquid. The Klinkenberg equation gives a semi-empirical relation between the liquid and gas permeabilities. In this paper, the wall-slip gas flow problem is homogenized. This problem is described by the steady state, low velocity Navier–Stokes equations for a compressible gas with a small Knudsen number. Darcy's law with a permeability tensor equal to that of liquid flow is shown to be valid to the lowest order. The lowest order wall-slip correction is a local tensorial form of the Klinkenberg equation. The Klinkenberg permeability is a positive tensor. It is in general not symmetric, but may under some conditions, which we specify, be symmetric. Our result reduces to the Klinkenberg equation for constant viscosity gas flow in isotropic media.     

Approximate kinetic equations in rarefied gas theory

The fundamental kinetic equation of gas theory, the Boltzniann equation, is a complex integrodiffcrential equation. The difficulties associated with its solution are the result not only of the large number of independent variables, seven in the general case, but also of the very complicated structure of the collision integral. However, for the mechanics of rarefied gases the primary interest lies not in the distribution function itself, which satisfies the Boltzmann equation, but rather in its first few moments, i.e., the averaged characteristics. This circumstance suggests the possibility of obtaining the averaged quantities by a simpler way than the direct method of direct solution of the Boltzmann equation with subsequent calculation of the integrals.It is well known that if a distribution function satisfies the Boltzmann equation, then its moments satisfy an infinite system of moment equations. Consequently, if we wish to obtain with satisfactory accuracy some number of first moments, then we must require that these moments satisfy the exact system of moment equations. However, this does not mean that to determine the moments of interest to us we must solve this system, particularly since the system of moment equations is not closed. The closure of the system by specifying the form of the distribution function (method of moments) can be considered only as a rough approximate method of solving problems. First, in this case it is not possible to satisfy all the equations and we must limit ourselves to certain of the equations; second, generally speaking, we do not know which equation the selected distribution function satisfies, and, consequently, we do not know to what degree it has the properties of the distribution function which satisfies the Boltzmann equation.A more reliable technique for solving the problems of rarefied gasdynamics is that based on the approximation of the Boltzmann equation, more precisely, the approximation of the collision integral. The idea of replacing the collision integral by a simpler expression is not new [1–4]. The kinetic equations obtained as a result of this replacement are usually termed model equations, since their derivation is usually based on physical arguments and not on the direct use of the properties of the Boltzmann collision integral. In this connection we do not know to what degree the solutions of the Boltzmann equation and the model equations are close, particularly since the latter do not yield the possibility of refining the solution. Exceptions are the kinetic model for the linearized Boltzmann equation [5] and the sequence of model equations of [6], constructed by a method which is to some degree analogous with that of [5].In the present paper we suggest for the simplification of the solution of rarefied gas mechanics problems a technique for constructing a sequence of approximate kinetic equations which is based on an approximation of the collision integral. For each approximate equation (i.e., equation with an approximate collision operator) the first few moment equations coincide with the exact moment equations. It is assumed that the accuracy of the approximate equation increases with increase of the number of exact moment equations. Concretely, the approximation for the collision integral consists of a suitable approximation of the reverse collision integral and the collision frequency. The reverse collision integral is represented in the form of the product of the collision frequency and a function which characterizes the molecular velocity distribution resulting from the collisions, where the latter is selected in the form of a locally Maxwellian function multiplied by a polynomial in terms of the components of the molecular proper velocities. The collision frequency is approximated by a suitable expression which depends on the problem conditions. For the majority of problems it may obviously be taken equal to the collision frequency calculated from the locally Maxwellian distribution function; if necessary the error resulting from the inexact calculation of the collision frequency may be reduced by iterations.To illustrate the method, we solve the simplest problem of rarefied gas theory-the problem on the relaxation of an initially homogeneous and isotropic distribution in an unbounded space to an equilibrium distribution.The author wishes to thank A. A. Nikol'skii for discussions of the study and V. A. Rykov for the numerical results presented for the exact solution.     

Thermo-Gravitational Convection in a Near-Critical Fluid in a Side-Heated Enclosed Cavity

Convection near a thermodynamic critical point in a square cavity with lateral heating is investigated numerically on the basis of the Navier-Stokes equations for a compressible gas with a Van-der-Waals equation of state. Comparison of a near-critical fluid and a perfect gas with parameters equal to those of the real fluid near the critical point shows that, with the development of convection, the dynamics of these two media are qualitatively different; however, a certain similarity is observed for the steady-state regime. The dependence of the steady-state flow and heat transfer characteristics on the nondimensional governing parameters is investigated.     

A method for studying waves with spatially localized envelopes in a class of nonlinear partial differential equations

A methodology for investigating stationary and travelling waves with spatially localized envelopes is presented. The nonlinear governing partial differential equations considered possess a constant first integral of motion, and are separable in space and time when the small parameter of the problem is set to zero. To study stationary waves, a coordinate transformation on the governing nonlinear partial differential equation is imposed which eliminates the time dependence from the problem. An amplitude modulation function is then introduced to express the response of the system at an arbitrary point as a nonlinear function of a reference response. Analytic approximations to the amplitude modulation function are developed by expressing it in power series, and asymptotically solving sets of singular functional equations at the various orders of approximation. Travelling solutions may be computed from stationary ones, by imposing appropriate Lorentz transformations. As an application of the methodology, stationary and travelling breathers of a nonlinear partial differential equation are analytically computed.     

A shape optimisation method of a body located in adiabatic flows

The purpose of this study is to derive an optimal shape of a body located in adiabatic flow. In this study, we use the equation of motion, the equation of continuity and the pressure–density relation derived from the Poisson’s law as the governing equation. The formulation is based on an optimal control theory in which a performance function of fluid force is taken into consideration. The performance function should be minimised satisfying the governing equations. This problem can be solved without constraints by using the adjoint equation with adjoint variables corresponding to the state equation. The performance function is defined by the drag and lift forces acting on the body. The weighted gradient method is applied as a minimisation technique, the Galerkin finite element method is used as a spatial discretisation and the implicit scheme is used as a temporal discretisation to solve the state equations. The mixed interpolation, the bubble function for velocity and the linear function for density, is employed as the interpolation. The optimal shape is obtained for a body in adiabatic flows.     

The effect of radiation on transient natural convection past a doubly infinite plate

The present research note is concerned with the transient (short time) simultaneous free convection and radiation analysis of a viscous fluid along a doubly infinite vertical isothermal flat plate. To simplify a very complicated problem, an incompressible flow field is used in the analysis. Generally, the exact numerical solution of this problem is quite lengthy. However, by considering an optically thick radiating gas, expressed by the Rosseland diffusion approximation, the solution is much simpler. Moreover, this case leads to a complete similarity transformation of the governing partial differential equations into a set of ordinary differential equations. An exact numerical solution is obtained of the resulting ordinary differential equations for a Prandtl number equal to 0.733 and for a wide range of involved parameters.     

Applications of the group of equations of motion of a polytropic gas

The machinery of Lie theory (groups and algebras) is applied to the system of equations governing the unsteady flow of a polytropic gas. The action on solutions of transformation groups which leave the equations invariant is considered. Using the invariants of the transformation groups, various symmetry reductions are achieved in both the steady state and the unsteady cases. These reduce the system of partial differential equations to systems of ordinary differential equations for which some closed-form solutions are obtained. It is then illustrated how each solution in the steady case gives rise to time-dependent solutions.     

Characteristic-value analysis for plastic dynamic buckling of columns under elastoplastic compression waves

The characteristic-value analysis of plastic dynamic buckling is presented for columns under the action of elastoplastic compression wave caused by an axial-step load. Two critical conditions constituting a dynamic instability criterion are derived on the basis of transformation and conservation of energy. The governing equations, the boundary conditions and the continuity conditions derived by the use of the first critical condition are the same as those given by the adjacent-equilibrium criterion and are insufficient for determining two characteristic parameters involved in the governing equations. A supplementary restraint equation for buckling deformations at the plastic-wave front and the elastic-wave front is derived by the use of the second critical condition. Then, a couple of characteristic equations for two characteristic parameters are derived on the condition that the governing equations have non-trivial solutions satisfying the boundary conditions, the continuity conditions and the supplementary restraint equation. The critical-load parameters, dynamic characteristic parameter (exponent parameter of inertia term) and dynamic buckling modes are calculated from the solutions of the characteristic equations.     

Chemical relaxation in a viscous shock

We consider the flow of a nonequilibrium dissociating diatomic gas in a normal compression shock with account for viscosity and heat conductivity. The distribution of gas parameters in the flow is found by numerically solving the Navier-Stokes and chemical kinetics equations. The greatest difficulty in numerical integration comes from the singular points of this system at which the initial conditions are given. These points lead to instability of the numerical results when the problem is solved by standard numerical methods. An integration method is proposed that yields stable numerical results-continuous profiles of the distribution of the basic gas parameters in the shock are obtained.We consider steady one-dimensional flow in which the gas passes from equilibrium state 1 to another equilibrium state 2, which has higher values for temperature, density, and pressure. Such a flow is termed a normal compression shock.The parameter distribution in normal shock for nonequilibrium chemical processes has usually been calculated [1–3] without account for the transport phenomena (viscosity, heat conduction, and diffusion). The presence of an infinitely thin shock front perpendicular to the flow velocity direction was postulated. It was assumed that the flow is undisturbed ahead of the shock front. The gas parameters (velocity, density, and temperature) change discontinuously across the shock front, but the gas composition does not change. The composition change due to reactions takes place behind the shock front. The gas parameter distribution behind the front was calculated by solving the system of gasdynamic and chemical kinetics equations using the initial values determined from the Hugoniot conditions at the front to state 2 far downstream.Several studies (for example, [4, 5]) do account for transport phenomena in calculating parameter distribution in a compression shock, but not for nonequilibrium chemical reactions. These problems are solved by integrating the Navier-Stokes equations continuously from state 1 in the oncoming flow to state 2 downstream.We present a solution to the problem of normal compression shock in nonequilibrium dissociating oxygen with account for viscosity and heat conduction using the Navier-Stokes equations.     

Delta Shock Wave Solution of the Riemann Problem for the Non-homogeneous Modified Chaplygin Gasdynamics

The motivation of the present study is to derive the solution of the Riemann problem for modified Chaplygin gas equations in the presence of constant external force. The analysis leads to the fact that in some special circumstances delta shock appears in the solution of the Riemann problem. Also, the Rankine–Hugoniot relations for delta shock wave which are utilized to determine the strength, position and propagation speed of the delta shocks have been derived. Delta shock wave solution to the Riemann problem for the modified Chaplygin gas equation is obtained. It is found that the external force term, appearing in the governing equations, influences the Riemann solution for the modified Chaplygin gas equation.

     

Quasi‐One‐Dimensional Riemann Problems and Their Role in Self‐Similar Two‐Dimensional Problems

We study two‐dimensional Riemann problems with piecewise constant data. We identify a class of two‐dimensional systems, including many standard equations of compressible flow, which are simplified by a transformation to similarity variables. For equations in this class, a two‐dimensional Riemann problem with sectorially constant data becomes a boundary‐value problem in the finite plane. For data leading to shock interactions, this problem separates into two parts: a quasi‐one‐dimensional problem in supersonic regions, and an equation of mixed type in subsonic regions. We prove a theorem on local existence of solutions of quasi‐one‐dimensional Riemann problems. For 2 × 2 systems, we generalize a theorem of Courant & Friedrichs, that any hyperbolic state adjacent to a constant state must be a simple wave. In the subsonic regions, where the governing equation is of mixed hyperbolic‐elliptic type, we show that the elliptic part is degenerate at the boundary, with a nonlinear variant of a degeneracy first described by Keldysh. (Accepted December 4, 1997)     

弹性压应力波下直杆动力失稳的机理的判据

基于应力波理论和失稳瞬间能量的转换和守恒,导出了一个直杆动力分岔失稳的准则：（1）直杆在发生分岔失稳的瞬间所释放出的压缩变形能等于屈曲所需变形能与屈曲动能之和;（2）在上述能量转换过程中,能量对时间的变化率服从守恒定律。应用临界条件（1）推导出的直杆动力失稳的控制方程和杆端边界条件以及连续条件,与应用哈密顿原理推导的结果完全相同,但不足以构成求解直杆动力失稳问题的完备定解条件,导出包含两个特征参数的一对特征方程。从而建立了求解直杆动力失稳模态和两个特征参数（临界力参数和失稳惯性项指数参数即动力特征参数）的较严密理论方法。     

Finite Oscillations of a Cylinder in a Coaxial Duct Subjected to Annular Compressible Flow

As a generalization considering small fluid-structural vibrations, the present paper examines the finite magnitude oscillatory motion of an elastically supported rigid cylinder in a cylindrical rigid duct conveying a compressible flow. The fluid is assumed to be inviscid and irrotational and free purely transverse vibrations of the body are dealt with. The governing equations of motion are the fully nonlinear Euler equations together with the continuity equation and a state equation (here for an ideal gas), the ordinary differential equation for the vibrating cylinder, and the kinematical transition and boundary conditions at the moving contact interface between fluid and body and the outside fluid border, respectively. A pertubation analysis is performed to calculate not only the dynamic characteristics for small coupled oscillations but also the corrections due to the inherent nonlinearities of the vibroacoustic problem. To make the calculation steps more transparent, the simpler problem of a two-dimensional channel flow between a rigid wall and an elastically supported rigid plate is also included in the present study. As an outlook, the influence of flexibility of the cylinder (or the plate) is addressed and the problem of forced vibrations is touched. This revised version was published online in July 2006 with corrections to the Cover Date.     

Multicomponent fluid flow analysis using a new set of conservation equations

In this work hydrodynamics of multicomponent ideal gas mixtures have been studied. Starting from the kinetic equations, the Eulerian approach is used to derive a new set of conservation equations for the multicomponent system where each component may have different velocity and kinetic temperature. The equations are based on the Grad's method of moment derived from the kinetic model in a relaxation time approximation (RTA). Based on this model which contains separate equation sets for each component of the system, a computer code has been developed for numerical computation of compressible flows of binary gas mixture in generalized curvilinear boundary conforming coordinates. Since these equations are similar to the Navier-Stokes equations for the single fluid systems, the same numerical methods are applied to these new equations. The Roe's numerical scheme is used to discretize the convective terms of governing fluid flow equations. The prepared algorithm and the computer code are capable of computing and presenting flow fields of each component of the system separately as well as the average flow field of the multicomponent gas system as a whole. Comparison of the present code results with those of a more common algorithm based on the mixture theory in a supersonic converging-diverging nozzle provides the validation of the present formulation. Afterwards, a more involved nozzle cooling problem with a binary ideal gas (helium-xenon) is chosen to compare the present results with those of the ordinary mixture theory. The present model provides the details of the flow fields of each component separately which is not available otherwise. It is also shown that the separate fluids treatment, such as the present study, is crucial when considering time scales on the order of (or shorter than) the intercollisions relaxation times.     

Solution of the Rayleigh problem for a powerlaw non-Newtonian conducting fluid via group method

ＩｎｔｒｏｄｕｃｔｉｏｎＦｌｕｉｄｓｔｈａｔｏｂｅｙＮｅｗｔｏｎ’ｓｌａｗｏｆｖｉｓｃｏｓｉｔｙａｒｅｃａｌｌｅｄＮｅｗｔｏｎｉａｎｆｌｕｉｄｓ.Ｎｅｗｔｏｎ’ｓｌａｗｏｆｖｉｓｃｏｓｉｔｙｉｓτ=μｄｕ/ｄｔ,ｗｈｅｒｅτｉｓｔｈｅｓｈｅａｒｓｔｒｅｓｓａｎｄμｉｓｔｈｅｖｉｓｃｏｓｉｔｙ .ＮｏｔａｌｌｆｌｕｉｄｓｆｏｌｌｏｗｔｈｅＮｅｗｔｏｎｉａｎｓｔｒｅｓｓ＿ｓｔｒａｉｎｒｅｌａｔｉｏｎ .Ｓｏｍｅｆｌｕｉｄｓ ,ｓｕｃｈａｓｋｅｔｃｈｕｐ ,ａｒｅｓｈｅａｒｔｈｉｎｎｉｎｇ ,ｔｈａｔｉｓ,ｔｈｅｃ…     

Fracture analysis of a functionally graded strip with arbitrary distributed material properties

In this paper, the plane elasticity problem for a crack in a functionally graded strip with material properties varying arbitrarily is studied. The governing equation in terms of Airy stress function is formulated and exact solutions are obtained for several special variations of material properties in Fourier transformation domain. A multi-layered model is employed to model arbitrary variations of material properties based on two linear-distributed material softness parameters. The mixed boundary problem is reduced to a system of singular integral equations that are solved numerically. Comparisons with other two existing multi-layered models have been made. Some numerical examples are given to demonstrate the accuracy, efficiency and versatility of the model. Numerical results show that fracture toughness of materials can be greatly improved by graded variation of elastic modulus and the influence of the specific form of elastic modulus on the fracture behavior of FGM is limited.     

Non-axisymmetrical vibration of elastic circular plate on layered transversely isotropic saturated ground

The non-axisymmetrical vibration of elastic circular plate resting on a layered transversely isotropic saturated ground was studied.First,the 3-d dynamic equations in cylindrical coordinate for transversely isotropic saturated soils were transformed into a group of governing differential equations with 1-order by the technique of Fourier ex- panding with respect to azimuth,and the state equation is established by Hankel integral transform method,furthermore the transfer matrixes within layered media are derived based on the solutions of the state equation.Secondly,by the transfer matrixes,the general solutions of dynamic response for layered transversely isotropic saturated ground excited by an arbitrary harmonic force were established under the boundary conditions, drainage conditions on the surface of.ground as well as the contact conditions.Thirdly, the problem was led to a pair of dual integral equations describing the mixed boundary- value problem which can be reduced to the Fredholm integral equations of the second kind solved by numerical procedure easily.At the end of this paper,a numerical result concerning vertical and radical displacements both the surface of saturated ground and plate is evaluated.     

Analytical Solution of Compression,Free Swelling and Electrical Loading of Saturated Charged Porous Media

Analytical solutions are derived for one-dimensional consolidation, free swelling and electrical loading of a saturated charged porous medium. The governing equations describe infinitesimal deformations of linear elastic isotropic charged porous media saturated with a mono-valent ionic solution. From the governing equations a coupled diffusion equation in state space notation is derived for the electro-chemical potentials, which is decoupled introducing a set of normal parameters, being a linear combination of the eigenvectors of the diffusivity matrix. The magnitude of the eigenvalues of the diffusivity matrix correspond to the time scales for Darcy flow, diffusion of ionic constituents and diffusion of electrical potential.     

Kinetic model for the motion of compressible bubbles in a perfect fluid

Collective behavior of compressible gas bubbles moving in an inviscid incompressible fluid is studied. A kinetic approach is employed, based on an approximate calculation of the fluid flow potential and formulation of Hamilton's equations for generalized coordinates and momenta of bubbles. Kinetic equations governing the evolution of a distribution function of bubbles are derived. These equations are similar to Vlasov's equations. Conservation laws which are direct consequences of the kinetic system are found. It is shown that for a narrowly peaked distribution function they form a closed system of hydrodynamical equations for the mean flow parameters. The system yields the analogue of Rayleigh–Lamb's equation governing oscillations of bubbles. A variational principle for the hydrodynamical system is established and the linear stability analysis is performed.