The stability of a free rotating layer in an electrically conducting fluid in an axial field

The subject considered is a homogeneous electrically conducting incompressible medium with a current in a homogeneous external magnetic field and bounded by parallel insulating planes normal to the induction vector. When the current is fed by means of a system of coaxial electrodes located on one or both of the insulating planes, regions arise in which the medium is in rotational motion. If the lateral wall is at a sufficient distance from the electrodes, the rotating layer which forms as a result of the interaction of the axial magnetic field and the radial component of the electric current has free lateral boundaries. A study is made of the way in which the Reynolds number for the loss of stability in such a layer depends on the Hartmann number and on the geometric parameter for high values of the Hartmann number and low values of the magnetic Reynolds number.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 166–173, September–October, 1984.     

End effects in magnetohydrodynamic channels at finite magnetic reynolds numbers

The effects of the magnetic Reynolds number have been examined via the distribution of the magnetic fields induced by the motion of a medium in a rectangular channel with conducting walls in the presence of an inhomogeneous magnetic field; the effects of wall conductivity and geometry of the external field are also examined as regards the distribution of the induced currents, the Joule loss, and the electric and magnetic fields over the cross section. The problem has previously been considered for a channel with insulating walls [1].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 19–27, May–June, 1971.We are indebted to A. B. Vatazhin for his interest.     

Simply approximate estimation of local failure position of a free-free slender shell

Dynamic plastic failure characteristics of a space free-free slender shell subjected to intense dynamic loading of suddenly applied pressure unsymmetrical triangle distributed along its span was studied. Both rigid perfectly plastic (r-p-p) analytical method and finite element method based elastic perfectly plastic (e-p-p) material idealization and shell element model were adopted to predict the local failure position in the structure. It was shown that both r-p-p and e-p-p model could estimate a plastic “kink” taking place in the slender shell, which reflects the strain localization of deformation. The comparison for the position of “kink” predicted by using r-p-p and e-p-p methods is found to be reasonable good.     

Outflow boundary conditions for one-dimensional finite element modeling of blood flow and pressure waves in arteries

Flow and pressure waves, originating due to the contraction of the heart, propagate along the deformable vessels and reflect due to tapering, branching, and other discontinuities. The size and complexity of the cardiovascular system necessitate a “multiscale” approach, with “upstream” regions of interest (large arteries) coupled to reduced-order models of “downstream” vessels. Previous efforts to couple upstream and downstream domains have included specifying resistance and impedance outflow boundary conditions for the nonlinear one-dimensional wave propagation equations and iterative coupling between three-dimensional and one-dimensional numerical methods. We have developed a new approach to solve the one-dimensional nonlinear equations of blood flow in elastic vessels utilizing a space-time finite element method with GLS-stabilization for the upstream domain, and a boundary term to couple to the downstream domain. The outflow boundary conditions are derived following an approach analogous to the Dirichlet-to-Neumann (DtN) method. In the downstream domain, we solve simplified zero/one-dimensional equations to derive relationships between pressure and flow accommodating periodic and transient phenomena with a consistent formulation for different boundary condition types. In this paper, we also present a new boundary condition that accommodates transient phenomena based on a Green’s function solution of the linear, damped wave equation in the downstream domain.     

Passive effects of rare-earth permanent magnets on flexible conductive structures

The paper presents a theoretical and experimental study of vibrating structures where paramagnetic or diamagnetic systems interact with rare-earth passive magnets.The theoretical model of the system is focused on the damping properties of permanent magnets and on their interactions with the dynamic behaviour of an Euler–Bernoulli beam. In particular, the magnetic model is based on the analogy of the equivalent currents method in a quasi-static open-circuit-type configuration and it is used to determine the influence of eddy currents on the dynamic behaviour of conducting material structures. The magnetic effects are characterised by a viscous-type damping and by a stiffening dynamic effect of the structure, called “phantom effect”.The authors present the experimental outcomes for uniform cantilever clamped-free beams of different kinds of paramagnetic or diamagnetic conducting materials. It appears that the system frequency response can be modified by the presence of a pair of concordant or discordant permanent magnets of high residual induction settled at the free end.Through the comparison between theoretical and experimental results, the paper demonstrates the validity of the model, that is able to describe both the above mentioned effect of dynamic stiffening of the structure and the considerable localised damping properties in paramagnetic or diamagnetic materials having low electric resistivity.     

Elastic wave scattering by a three-dimensional inhomogeneity in an elastic half space

In this paper we will consider scattering of elastic waves in a half space. The half space is an isotropic, linear and homogeneous medium except for a finite inhomogeneity. The T-matrix method (also called the “extended boundary condition method” or “null field approach”) is extended to derive expressions for the elastic field inside the half space and the surface field on the interface. The assumptions on the source that excites the half space are fairly weak. In the numerical applications found in this paper we assume a Rayleigh surface wave to be the incoming field, and we only compute the surface displacements. We make illustrations on some simple types of scatterers (spheres and spheroids; the latter ones can be arbitrarily oriented).     

Finite deformation of incompressible fiber-reinforced elastomers: A computational micromechanics approach

The in-plane finite deformation of incompressible fiber-reinforced elastomers was studied using computational micromechanics. Composite microstructure was made up of a random and homogeneous dispersion of aligned rigid fibers within a hyperelastic matrix. Different matrices (Neo-Hookean and Gent), fibers (monodisperse or polydisperse, circular or elliptical section) and reinforcement volume fractions (10–40%) were analyzed through the finite element simulation of a representative volume element of the microstructure. A successive remeshing strategy was employed when necessary to reach the large deformation regime in which the evolution of the microstructure influences the effective properties. The simulations provided for the first time “quasi-exact” results of the in-plane finite deformation for this class of composites, which were used to assess the accuracy of the available homogenization estimates for incompressible hyperelastic composites.     

Joule loss due to compressibility of gas in a channel of variable section

It is known that closed electric currents arise in a conducting medium moving in a non-uniform magnetic field. These currents lead to additional energy loss and adversely affect the characteristics of magnetohydrodynamic channels. (The numerous investigations of these effects are dealt with in the review [2, 3].) Eddy electric currents are also formed, however, when a medium flows in a uniform magnetic field perpendicular to the to the plane of motion if the channel has a variable cross section and the medium is compressible [1], This paper is devoted to an investigation of some features of these flows. It is assumed in the analysis that the gas flows in channels whose geometry varies slightly.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, Vol. 9, No. 4, pp. 3–9, July–August, 1968.     

Geomechatronics – Interaction between ground and construction machinery and its application to construction robotics

“Geomechatronics” is a technical field in which “Geotechniques” is fused with “Mechatronics” that is the technical field to promote the automatic control of machines by using the electronics. In the field of “Geomechatronics”, a construction machine, which treats geotechnical materials such as soil and rock, automatically evaluates the properties and conditions of the ground and determines the optimum controlling method of itself for the ground with the base of the machine–ground interaction. Some researches for practical use in the field of geomechatronics are introduced, and then the progressing view of this research and technical filed is explained in this paper.     

Effect of ion slip on the time-varying Hartmann flow of a non-Newtonian viscoelastic fluid with heat transfer

Ion slip in a time-varying Hartmann flow of a conducting incompressible non-Newtonian viscoelastic fluid between two parallel horizontal insulating porous plates is studied with allowance for heat transfer. A uniform and constant pressure gradient is applied in the axial direction. An external uniform magnetic field and uniform suction and injection through the surface of the plates are applied in the normal direction. The two plates are maintained at different but constant temperatures; the Joule and viscous dissipations are taken into consideration. Numerical solutions for the governing momentum and energy equations are obtained with the use of finite differences, and the effect of various physical parameters on both the velocity and temperature fields is discussed.     

On constitutive modelling of porous neo-Hookean composites

In this paper a hyperelastic constitutive model is developed for neo-Hookean composites with aligned continuous cylindrical pores in the finite elasticity regime. Although the matrix is incompressible, the composite itself is compressible because of the existence of voids. For this compressible transversely isotropic material, the deformation gradient can be decomposed multiplicatively into three parts: an isochoric uniaxial deformation along the preferred direction of the material (which is identical to the direction of the cylindrical pores here); an equi-biaxial deformation on the transverse plane (the plane perpendicular to the preferred direction); and subsequent shear deformation (which includes “along-fibre” shear and transverse shear). Compared to the multiplicative decomposition used in our previous model for incompressible fibre reinforced composites [Guo, Z., Peng, X.Q., Moran, B., 2006, A composites-based hyperelastic constitutive model for soft tissue with application to the human annulus fibrosus. J. Mech. Phys. Solids 54(9), 1952–1971], the equi-biaxial deformation is introduced to achieve the desired volume change. To estimate the strain energy function for this composite, a cylindrical composite element model is developed. Analytically exact strain distributions in the composite element model are derived for the isochoric uniaxial deformation along the preferred direction, the equi-biaxial deformation on the transverse plane, as well as the “along-fibre” shear deformation. The effective shear modulus from conventional composites theory based on the infinitesimal strain linear elasticity is extended to the present finite deformation regime to estimate the strain energy related to the transverse shear deformation, which leads to an explicit formula for the strain energy function of the composite under a general finite deformation state.     

The emergence of shear bands in plane strain

In this paper, we examine the effect of a slight material imperfection on the deformation field in an otherwise homogeneous body subjected to a plane bi-axial stretch at infinity. Both hyper-elastic and hypo-elastic materials are considered, with the constitutive equations assumed to be such that the governing equilibrium equations lose ellipticity at some strain level. A straightforward regular perturbation analysis is performed and attention is focussed on the features of the first order terms. It is found that the effect of the imperfection is negligible at small values of the applied load. As the load increases, the imperfection more or less abruptly gets “activated” and causes a rapid concentration of strain within certain narrow bands—shear bands—passing through the imperfection. In order to estimate the accuracy of the linearized analysis, a finite element solution of the nonlinear problem is also carried out and results are compared.     

Orientation effect on the Boussinesq indentation of a transversely isotropic piezoelectric material

Three-dimensional finite element analysis was used to study the effect of the angle between the loading direction and the axisymmetric direction on the indentation behavior of a transversely isotropic piezoelectric half-space by a cylindrical indenter of flat end. Two cases were considered in the analysis, which included (a) the indentation by an insulating indenter, and (b) the indentation by a conducting indenter. Both the indentation load and the indentation-induced potential were found to be proportional to the indentation depth. Using the simulation results and the analytical relationship for the indentation by a rigid, insulating indenter, semi-analytical relationships were developed between the indentation load and the indentation depth and between the indentation-induced potential on the indenter and the indentation depth for the conducting indenter, respectively. The proportionality between the indentation-induced potential and the indentation depth is only a function of the angle between the loading direction and the poling direction, independent of the type of indenters, which may be used to measure the relative direction of the loading axis to the axisymmetric axis of transversely piezoelectric materials from the indentation test.     

压电薄板板边导电裂纹尖端的力电场分析

采用有限元方法，分析了压电薄板板边不同长度导电裂纹尖端的力电场分布规律，发现导电裂纹尖端的应力场和电场强度存在明显的集中和奇异现象，集中和奇异的程度与裂纹长度有关。而且，在裂纹延长线上分别存在两点，这里的应力和电场对裂纹长度不太敏感，总等于无裂纹时薄板的均匀应力和均匀电场强度；同时，还研究了导电裂纹尖端的应力强度因子和电场强度因子对裂纹长度的依赖关系，发现在线性本构的前提下，导电裂纹尖端的应力强度因子与电场强度因子之间具有近似的线性关系。     

Electromechanical behaviour of a finite piezoelectric layer under a flat punch

This paper solves the problem of a smooth and frictionless punch on a piezoelectric ceramic layer. Different electrical boundary conditions that employ conducting or insulating punches are presented. The stress and electric displacement intensity factors are used to characterize the electromechanical fields at the punch tip. The field intensity factors are obtained numerically for finite layer thickness. Effects of the thickness of the piezoelectric layer on the stress and electric displacement, and the stress and electric displacement intensity factors at the punch tip are discussed. Solution technique for two identical and collinear surface punches on the piezoelectric layer is also provided and the effect of relative distance between the punches is investigated. Numerical results for some interesting special cases, such as conducting punch and insulating punch, and infinite piezoelectric layer thickness, are presented.     

Stabilized finite element formulation applied to the kinematic Ponomarenko dynamo problem

A stabilized finite element (B, q) formulation is developed to solve the kinematic dynamo problem. As a test case, we solve the induction equation for a given solid body helical flow, embedded in a cylindrical conducting shell. This problem corresponds to the well-known Ponomarenko dynamo. It has the interesting property to have an exact dispersion relation giving the magnetic growth rate as a function of the flow properties. Therefore, it is a good benchmark to test our kinematic dynamo code. We calculated the dynamo threshold and plotted the geometry of the generated magnetic field. We also evaluated the residual error due to our stabilized formulation.     

Steady shear banding in complex fluids

Shear banding (SB) is manifested by the abrupt “demixing” of the flow into regions of high and low shear rate. In this paper, we first relate analytically the rheological parameters of the fluid with the range of shear rates and stresses of SB occurrence. For this, we accept that the origin of shear banding is constitutive, and adopt a non-linear viscoelastic expression able to accommodate the double-valuedness of the stress with flow intensity, under certain conditions. We then implement the model for the case of pressure driven flow through a cylindrical pipe; we derive approximate expressions for the velocity profile in the two-banded regions (core and outer annular), the overall throughput in the presence or absence of “spurt”, and the radial location limits of the shear rate discontinuity.     

Bounds for the effective properties of heterogeneous plates

This paper presents new bounds for heterogeneous plates which are similar to the well-known Hashin–Shtrikman bounds, but take into account plate boundary conditions. The Hashin–Shtrikman variational principle is used with a self-adjoint Green-operator with traction-free boundary conditions proposed by the authors. This variational formulation enables to derive lower and upper bounds for the effective in-plane and out-of-plane elastic properties of the plate. Two applications of the general theory are considered: first, in-plane invariant polarization fields are used to recover the “first-order” bounds proposed by Kolpakov [Kolpakov, A.G., 1999. Variational principles for stiffnesses of a non-homogeneous plate. J. Meth. Phys. Solids 47, 2075–2092] for general heterogeneous plates; next, “second-order bounds” for n-phase plates whose constituents are statistically homogeneous in the in-plane directions are obtained. The results related to a two-phase material made of elastic isotropic materials are shown. The “second-order” bounds for the plate elastic properties are compared with the plate properties of homogeneous plates made of materials having an elasticity tensor computed from “second-order” Hashin–Shtrikman bounds in an infinite domain.     

Legendre integrators, post-processing and quasiequilibrium

A toolbox for the development and reduction of the dynamical models of nonequilibrium systems is presented. The main components of this toolbox are: Legendre integrators, dynamical post-processing, and the thermodynamic projector. The thermodynamic projector is the tool to transform almost any anzatz to a thermodynamically consistent model. The post-processing is the cheapest way to improve the solution obtained by the Legendre integrators. Legendre integrators give the opportunity to solve linear equations instead of nonlinear ones for quasiequilibrium (“maximum entropy”, MaxEnt) approximations. The essentially new element of this toolbox, the method of thermodynamic projector, is demonstrated on application to the FENE-P model of polymer kinetic theory. The multi-peak model of polymer dynamics is developed.     

MHD Couette two-fluid flow and heat transfer in presence of uniform inclined magnetic field

The MHD Couette flow of two immiscible fluids in a parallel plate channel in the presence of an applied electric and inclined magnetic field is investigated in the paper. One of the fluids is assumed to be electrically conducting, while the other fluid and the channel plates are assumed to be electrically insulating. Separate solutions with appropriate boundary conditions for each fluid are obtained and these solutions are matched at the interface using suitable matching conditions. The partial differential equations governing the flow and heat transfer are transformed to ordinary differential equations and closed-form solutions are obtained in both fluid regions of the channel. The results for various values of the Hartmann number, the angle of magnetic field inclination, the loading parameter and the ratio of the heights of the fluids are presented graphically to show their effect on the flow and heat transfer characteristics.