Overall stability of pin-jointed frameworks after the onset of elastic buckling

Conclusions The total potential energy surfaces (V-surfaces) in the adjacent neighborhood of the initial position of the two perfect pin-jointed frames studied in examples 1 and 2 are confined to a three-dimensional space since only two of the members were assumed to be on the point of buckling. If however more than two members of a frame can contract flexurally, the representation of the V-surfaces is not so simple. The negative regions of the V-surface are of primary interest in the study of the post-buckling hehaviour of a frame since then an accelerated motion of the system away from the initial position of the frame may ensue. Therefore in the case of a multi-dimensional surface the location of the negative regions of the surface and the evaluation of the corresponding unstable modes of the given frame may be carried out more conveniently analytically using the expressions given by (2.8) and (2.9).In order to prove that the system is in unstable equilibrium, when the frame is on the point of buckling, it suffices to find at least one adjacent position for which V is negative. This may be a much easier task in the case of a complex frame than representing the entire shape of the surface in a multi-dimensional adjacent space so that the choice of a particular post-buckling mode or the possibility of a snap-through from one mode into another may be fully understood. Similarly, in the case of stable equilibrium of the system when a given frame is on the point of buckling, the result may be obtained quickly if the test for positive definiteness of the quadratic from in the expression for V is carried out.The effect of initial imperfections in the members becomes apparent on comparing the loaddisplacement characteristic of the frame in example 1 with the slope to the equilibrium path for the perfect frame. It is observed that adequate agreement between the two is established when the buckling mode is more fully developed.Then clearly, the initial imperfections in the members can be expected to affect the initial regions of the total potential energy surfaces, but agreement between such surfaces for perfect as well as in perfect frames may be adequate after the initial distortions of the imperfect frame have been overcome.Consequently the V surfaces as well as the stability criterion formulated for perfect frames by expressions 2.10 and 2.11 will in general give an indication of the post-buckling behavior of pin-jointed elastic frames when their members possess initial imperfections.This paper represents a part of the Investigation into the Post Buckling Behavior of Frames sponsored by the Aluminium Development Association of Great Britain and carried out by the author at Cambridge University from the years 1957 to 1960.     

Structure of a collisionless boundary layer and turbulent damping of ions

The flow of a low-pressure plasma in a MHD channel is unstable in a number of cases. The instability can be caused by a current flowing across the magnetic field. In this study we investigate an unstable, turbulent flow of a rarefied plasma near the magnetized electrodes, representing plane magnetic dipoles. Owing to the growth of microscopic turbulence near the electrodes, the maximum density of the current that is induced in the plasma is localized and turbulent damping of the incoming flow occurs. The energy of damping goes into the turbulent heating of the plasma. Under these conditions a structure of the boundary layer is found for a stationary flow. The characteristic transverse dimension of the boundary layer is considerably less than the particle mean free path; therefore, such a boundary layer can be called collisionless.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 5, pp. 17–28, September–October, 1972.     

Dynamic Buckling of Autonomous Systems Having Potential Energy Universal Unfoldings of Cuspoid Catastrophe

This work deals with dynamic buckling universal solutions of discrete nondissipative systems under step loading of infinite duration. Attention is focused on total potential energy functions associated with universal unfoldings of cuspoid type catastrophes with one active coordinate. The fold, dual cusp and tilted cusp catastrophes under statically applied loading occurring via limit points, asymmetric/symmetric bifurcations and nondegenerate hysteresis points are extended to the case of dynamic loading. Catastrophe manifolds of these types showing imperfection sensitivity under both types of loading are fully assessed. Important findings regarding dynamic buckling of imperfect systems generated by perfect systems associated with imperfect bifurcations are explored. The analysis is supplemented by a numerical application of a system exhibiting imperfect bifurcation when it is perfect as well as a hysteresis point associated with a tilted cusp catastrophe, when it becomes imperfect.     

On the nonlinear dynamic buckling mechanism of autonomous dissipative/nondissipative discrete structural systems

Summary Nonlinear dynamic buckling of nonlinearly elastic dissipative/nondissipative multi-mass systems, mainly under step load of infinite duration, is studied in detail. These systems, under the same loading applied statically, experience a limit point instability. The analysis can be readily extended to the case of dynamic buckling under impact loading. Energy, topological and geometrical aspects for the total potential energyV, which is constrained to lie in a region of phase-space whereV0, allow conclusions to be drawn directly regarding dynamic buckling. Criteria leading to very good, approximate and lower/upper bound dynamic buckling estimates are readily established without solving the highly nonlinear set of equations of motion. The theory is illustrated with several analyses of a two-degree-of-freedom model.     

Shear-induced aggregation and break up of fibril clusters close to the percolation concentration

To probe the behaviour of fibrillar assemblies of ovalbumin under oscillatory shear, close to the percolation concentration, c_p (7.5%), rheo-optical measurements and Fourier transform rheology were performed. Different results were found close to c_p (7.3%), compared to slightly further away from c_p (6.9 and 7.1%). For 6.9 and 7.1%, a decrease in complex viscosity, and a linear increase in birefringence, n, with increasing strain was observed, indicating deformation and orientation of the fibril clusters. For 7.3%, a decrease in complex viscosity was followed by an increase in complex viscosity with increasing strain, which coincided with a strong increase in n, dichroism, n, and the intensity of the normalized third harmonic (I₃/I₁). This regime was followed by a second decrease in complex viscosity, where n,n and I₃/I₁ decreased. In the first regime where the viscosity was decreasing with increasing strain, deformation and orientation of existing clusters takes place. At higher oscillatory shear, a larger deformation occurs and larger structures are formed, which is most likely aggregation of the clusters. Finally, at even higher strains, the clusters break up again. An increase in complex viscosity, n, n and I₃/I₁ was observed when a second strain sweep was performed 30 min after the first. This indicates that the shear-induced cluster formation and break up are not completely reversible, and the initial cluster size distribution is not recovered after cessation of flow.     

Stability of the tumbled state in shear flow of a nematic liquid crystal

In shear flow of a nematic liquid crystal with ₃ 0 flow alignment cannot occur. The stability of the stationary in-plane solution, the tumbled state, is investigated using abstract techniques. Employing the existence of an elastic energy a sufficient criterion for stability is formulated. This criterion depends on the in-plane solution which is obtainable as a quadrature that is non-elementary except in special cases. It is shown that the tumbled state is stable and asymptotically stable for some physical configurations. The criterion presented is not a necessary condition for stability and thus only gives a lower bound.     

Nonlinear stabilization of unstable acoustical vibrations in a bounded heat-evolving medium

The working process in many power units, taking place with the evolution of energy, under determined conditions is acoustically unstable in a linear approximation. The problem of determination of the amplitudes of the unstable waves, set up as a result of the pumping of energy from the unstable mode to the damping mode with their linear interaction is of practical interest. In this paper equations are derived for the fully established amplitudes of plane acoustical vibrations in a three-wave approximation, taking account of boundary impedances, breaking down the internal resonance of the acoustical overtones excited. The discussion regards a high-temperature heat-evolving gas, whose general stability conditions were formulated in [1] and discussed also in [2, 3].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 34–41, November–December, 1978.     

Heat transfer for laminar flow in internally finned pipes with different fin heights and uniform wall temperature

An analysis is presented for fully developed laminar convective heat transfer in a pipe provided with internal longitudinal fins, and with uniform outside wall temperature. The fins are arranged in two groups of different heights. The governing equations have been solved numerically to obtain the velocity and temperature distributions. The results obtained for different pipe-fins geometries show that the fin heights affect greatly flow and heat transfer characteristics. Reducing the height of one fin group decreases the friction coefficient significantly. At the same time Nusselt number decreases inappreciably so that such reduction is justified. Thus, the use of different fin heights in internally finned pipes enables the enhancement of heat transfer at reasonably low friction coefficient.Nomenclature A_f dimensionless flow area of the finned pipe, Eq. (8) - a_f flow area of the finned pipe - C_p specific heat at constant pressure - f coefficient of friction, Eq. (12) - H₁, H₂ dimensionless fin height h₁/r_o h₂/r_o - h₁, h₂ fin heights - average heat transfer coefficient at solid-fluid interface - KR fin conductance parameter, k_s/k_f - k_f thermal conductivity of fluid - k_s thermal conductivity of fin - l pipe length - mass flow rate - N number of fins - Nu Nusselt number, Eqs. (15) and (16) - P pressure - Q total heat transfer rate at solid fluid interface - Q_f1, Q_f2 heat transfer rate at fin surface - q_w average heat flux at pipe-wall, Q/(2 r_ol) - R dimensionless radial coordinate r/r_o - Re Reynolds Number, Eq. (13) - r radial coordinate - r_o radius of pipe - r₁, r₂ radii of fin tips - T temperature - T_b bulk temperature - U dimensionless velocity, Eq. (2) - U_b dimensionless bulk velocity - u_z axial velocity - z axial coordinate - angle between the flanks of two adjacent fins - half the angle subtended by a fin - angle between the center-lines of two adjacent fins - angular coordinate - dynamic viscosity - density - dimensionless temperature, Eq. (6) - _b dimensionless bulk temperature     

Influence of Vibrational Relaxation on the Pulsation Activity in Flows of an Excited Diatomic Gas

The influence of vibrational relaxation on the nonlinear evolution of a large vortex structure in a shear flow of a highly nonequilibrium diatomic gas is studied. Calculations are performed using the equations of twotemperature gas dynamics for a viscous heatconducting gas. Relaxation of the temperature of vibrational levels of gas molecules to equilibrium is described by the Landau–Teller equation. The contribution of the relaxation of rotational levels is taken into account by the bulk viscosity in the stress tensor. It is shown that in the presence of only the relaxation process with no viscous dissipation, the damping of the kinetic energy of perturbations and Reynolds stresses increases by up to 10 % compared to the case of thermal equilibrium. For high (actually attainable) degrees of excitation of the vibrational mode, moderate dynamic and bulk viscosities, and a typical relaxation time comparable to flow time, the relative effect of perturbation damping reaches 15%.     

An application of non-singularity BEM to plate bending, buckling and natural vibration problems

In this paper the fundamental solution of the singular governing equation of plate static bending is taken as the Green's function, which can satisfy the governing equation precisely in the plate region. Based on the principle of superposition, let the function values on the plate boundary, induced by a set of the Green's function sources (including the known sources in the plate region and the unknown sources in the fictitious region), satisfy the prescribed conditions on specially chosen boundary matching points, and the corresponding semi-analytical and semi-numerical solution can be obtained, which is free from the restraint of boundary forms and boundary conditions. The more matching points there are on the boundary, the better the accuracy of results is. Finally, in static bending problems a set of linear algebraic equations has to be computed; in buckling problems the minimum value of buckling eigenvalue equation has to be found; in natural vibration problems the eigenvalues of the frequency equation have to be calculated. Numerical examples are given and the results are compared with those by the analytical method and other methods. It can be seen that they are very close to each other.     

Stability of elastic bodies under omnilateral compression

Conclusions We have analyzed here the stability of the equilibrium of a simply connected isotropic compressible body with the elastic potential of arbitrary form and under uniform omnilateral deformation. A survey has been given here of earlier results obtained by other authors. The basic celations have been stated in a general form covering the theory of finite subcritical strains and two variants of the theory of small subcritical strains. For the latter theory new relations have been rigorously derived from which perturbations of tracking surface loads can be calculated, on the basis of corresponding expressions in the theory of finite subcritical strains. It has been proven that the sufficient conditions for the applicability of the static method of analysis are satisfied when the same boundary conditions are given over the entire body surface, as well as in several cases of different boundary conditions given at different segments of the boundary surface. It has been shown in a general form, for the theory of finite subcritical strains and for two variants of the theory of small subcritical strains, that the equilibrium of an elastic body under omnilateral deformation is stable, if a tracking load, is given over the entire boundary surface. As an example of problems with different boundary conditions at different segments of the boundary surface, we have considered the conventional problem concerning the stability of a bar on hinge supports and under uniform omnilateral deformation. It has been rigorously proven that in this case the equilibrium is stable when tracking loads are given at the lateral surfaces and is unstable when dead loads are given at the lateral surfaces. These conclusions apply to the theory of finite subcritical strains as well as to the theory of small subcritical strains, and they represent the complete version pertaining to compressible bodies.Institute of Mechanics, Academy of Sciences of the Ukrainian SSR, Kiev. Translated from Prikladnaya Mekhanika, Vol. 12, No. 6, pp. 3–27, June, 1976.     

Validation of Thermal Equilibrium Assumption in Transient Forced Convection Flow in Porous Channel

The validity of the local thermal equilibrium assumption in the transient forced convection channel flow is investigated analytically. Closed form expressions are presented for the temperatures of the fluid and solid domains and for the criterion which insures the validity of the local thermal equilibrium assumption. It is found that four dimensionless parameters control the local thermal equilibrium assumption. These parameters are the porosity , the volumetric Biot number Bi, the dimensionless channel length _max and the solid to fluid total thermal capacity ratio C _R. The qualitative and quantitative aspects of the effects of these four parameters on the channel thermal equilibrium relaxation time are investigated.     

Numerical analysis of axisymmetric buckling of conical shells

Nonlinear boundaryvalue problems of axisymmetric buckling of conical shells under a uniformly distributed normal pressure are solved by the shooting method. The problems are formulated for a system of six firstorder ordinary differential equations with independent rotation and displacement fields. Simply supported and clamped cases are considered. Branching solutions of the boundaryvalue problems are studied for different pressures and geometrical parameters of the shells. The nonmonotonic and discontinuous curves of equilibrium states obtained show that collapse, i.e., snapthrough instability is possible. For a simply supported shell, multivalued solutions are obtained for both external and internal pressure. For a clamped thinwalled shell, theoretical results are compared with experimental data.     

Numerical Analysis of Axisymmetric Buckling of Plates under Radial Compression

Nonlinear boundaryvalue problems of axisymmetric buckling of simply supported and clamped plates under radial compression are formulated for a system of six firstorder ordinary differential equations with independent fields of finite displacements and rotations. Multivalued solutions are obtained by the shooting method with specified accuracy. Bifurcation of the solutions of the problem is studied, and a parametric bifurcation diagram is constructed for various values of the loading parameter. Curves of buckling modes are given for three branches of the solution. The numerical results agree with available theoretical data.     

Turbulent mixing,induced by the Richtmyer-Meshkov instability

Turbulent mixing development, induced by passing of shock waves through the interface, has been studied on the basis of two semi-empirical models. The process of transition from random initial perturbations to turbulence is a rather complicated one, but it proceeds in a finite time. Here it is assumed that passing of shock waves through the initial zone of roughness leads to the instant development of mixing. The expression for the turbulent kinetic energy, transfered by the arriving shock wave, has been derived. The value of this energy is determined by the value of the non-dimensional parameter which is equivalent to the Richardson number. The analytical solution describing the turbulent kinetic energy decay and the mixing zone width change in time, has been obtained.The analytical form of constructed solutions allows to perform comparison between the models and to make the choice of constants on the basis of comparison with the known experiments. Comparison of the obtained theoretical results with Zaitzev's experimental data is given for the impulsive acceleration.This article was processed using Springer-Verlag T_EX Shock Waves macro package 1.0 and the AMS fonts, developed by the American Mathematical Society.     

Supersonic air flow past a sphere with account for vibrational relaxation

A numerical solution is obtained for the problem of air flow past a sphere under conditions when nonequilibrium excitation of the vibrational degrees of freedom of the molecular components takes place in the shock layer. The problem is solved using the method of [1]. In calculating the relaxation rates account was taken of two processes: 1) transition of the molecular translational energy into vibrational energy during collision; 2) exchange of vibrational energy between the air components. Expressions for the relaxation rates were computed in [2]. The solution indicates that in the state far from equilibrium a relaxation layer is formed near the sphere surface. A comparison is made of the calculated values of the shock standoff with the experimental data of [3].Notation uV_max, vV_max velocity components normal and tangential to the sphere surface - V_max maximal velocity - P V _max ² pressure - density - TT temperature - e_viRT vibrational energy of the i-th component per mole (i=–O₂, N₂) - =rb–1 shock wave shape - a _f the frozen speed of sound - HRT/m gas total enthalpy     

Slow-motion manifolds,dormant instability,and singular perturbations

If the coexistence of two phases at the transition temperature is kept under observation for a long time, then one observes that the system is not exactly in equilibrium and a very slow evolution driven by surface tension is taking place. Theoretically, one should eventually see a spatially homogeneous state, but the time for settling down is so long that what one actually observes is motion towards a stable state. The complexity of the spatial distribution of the two phases keeps decreasing but appears to be stable for very long periods of time with intermittent periods of fast motion when there are small inclusions of one of the two regions embedded in the other phase. For a simple reaction diffusion model, it is shown that this phenomenon can be explained by investigating the flow on the attractor and the unstable manifolds of equilibria.     

Free convection above a horizontal circular disk in a saturated porous medium

The boundary-layer flow above a horizontal impermeable circular disk embedded in a saturated porous medium is considered in the cases both when the disk is held at a constant temperature above ambient and when heat is supplied to the convective fluid by the disk at a constant rate. Series solutions are obtained based on the flat plate solution, which holds at the edge of the disk, as the leading order terms. These series solutions can then be used to describe the flow nearly all the way across the disk. A simple approximate solution, based on an integrated form of the energy equation, is also obtained and is shown, for the constant wall temperature case, to be useful in indicating how the solution behaves near the centre of the disk. The solutions asr0, wherer measures distance from the centre is discussed in both cases, and it is shown that the solution develops a singularity with the boundary layer having a thickness of 0[(–logr)^1/2].

---

Freie Konvektion über einer waagerechten Kreisscheibe in einem gesättigten porösen Medium

Zusammenfassung Es wird die Grenzschicht über einer waagerechten undurchlässigen Kreisscheibe, die in einem gesättigten porösen Medium eingebettet ist, untersucht. Zwei Möglichkeiten werden betrachtet: bei der einen ist die Oberflächentemperatur der Scheibe konstant, aber größer als die der Umgebung, bei der anderen ist die Scheibe gleichförmig beheizt. Es werden Lösungen in Form von Reihenentwicklungen erzielt, wobei die ersten Glieder den Lösungen einer ebenen Platte entsprechen, die an dem Rand der Scheibe gültig sind. Danach werden diese Lösungen zur Beschreibung der Bewegung der Flüssigkeit über einen großen Bereich der Scheibe benutzt. Eine einfache Lösung für den Fall der konstanten Wandtemperatur wird durch Integration der Energiegleichung erhalten und zur Beschreibung des Verhaltens der Lösung in der Nähe der Scheibenmitte verwendet. Man untersucht die Lösung fürr0 in beiden betrachteten Fällen, wobeir der Abstand von der Scheibenmitte ist. Bei einer Dicke von 0[(–logr)^1/2] der Grenzschicht weist die Lösung eine Singularität auf.

---

Nomenclature a radius of disk - g acceleration due to gravity - K permeability of the porous medium - k thermal conductivity - q prescribed wall heat flux - Q (non-dimensional) heat transfer coefficient - r, coordinate measuring distance from the centre of the disk - R _a,R _a ^* Rayleigh number - T temperature of the convective fluid - T ₀ ambient temperature - T ₁ prescribed wall temperature - u Darcy's law velocity in ther-direction - V _w (non-dimensional) wall velocity - w Darcy's law velocity in thez-direction - x coordinate measuring distance from the edge of the disk - z, coordinate measuring distance normal to the disk - equivalent thermal diffusivity - coefficient of thermal expansion - non-dimensional temperature - _w (non-dimensional) wall temperature - viscosity of the convective fluid - stream-function - stream function at the edge of the boundary layer     

Pseudolinear vibroimpact systems: Non-white random excitation

Response analyses of vibroimpact systems to random excitation are greatly facilitated by using certain piecewise-linear transformations of state variables, which reduce the impact-type nonlinearities (with velocity jumps) to nonlinearities of the common type — without velocity jumps. This reduction permitted to obtain certain exact and approximate asymptotic solutions for stationary probability densities of the response for random vibration problems with white-noise excitation. Moreover, if a linear system with a single barrier has its static equilibrium position exactly at the barrier, then the transformed equation of free vibration is found to be perfectly linear in case of the elastic impact. The transformed excitation term contains a signature-type nonlinearity, which is found to be of no importance in case of a white-noise random excitation. Thus, an exact solution for the response spectral density had been obtained previously for such a vibroimpact system, which may be called pseudolinear, for the case of a white-noise excitation. This paper presents analysis of a lightly damped pseudolinear SDOF vibroimpact system under a non-white random excitation. Solution is based on Fourier series expansion of a signum function for narrow-band response. Formulae for mean square response are obtained for resonant case, where the (narrow-band) response is predominantly with frequencies, close to the system's natural frequency; and for non-resonant case, where frequencies of the narrow-band excitation dominate the response. The results obtained may be applied directly for studying response of moored bodies to ocean wave loading, and may also be used for establishing and verifying procedures for approximate analysis of general vibroimpact systems.     

Nonlinear dynamic buckling of pinned–fixed shallow arches under a sudden central concentrated load

The structural behavior of a shallow arch is highly nonlinear, and so when the amplitude of the oscillation of the arch produced by a suddenly-applied load is sufficiently large, the oscillation of the arch may reach a position on its unstable equilibrium paths that leads the arch to buckle dynamically. This paper uses an energy method to investigate the nonlinear elastic dynamic in-plane buckling of a pinned–fixed shallow circular arch under a central concentrated load that is applied suddenly and with an infinite duration. The principle of conservation of energy is used to establish the criterion for dynamic buckling of the arch, and the analytical solution for the dynamic buckling load is derived. Two methods are proposed to determine the dynamic buckling load. It is shown that under a suddenly-applied central load, a shallow pinned–fixed arch with a high modified slenderness (which is defined in the paper) has a lower dynamic buckling load and an upper dynamic buckling load, while an arch with a low modified slenderness has a unique dynamic buckling load.