导电圆柱薄壳的热磁弹性效应分析

研究了磁场环境中受机械载荷作用的导电圆柱薄壳的热磁弹性问题．首先,根据电动力学方程和广义Ohm定律,得到了导电薄壳电流密度的分布,考虑到Joule热效应及热平衡方程,得到了导电薄壳的温度分布．其次,通过几何方程、物理方程、运动方程和电动力学方程导出了导电薄壳在机械场、电磁场以及温度场作用下的基本方程．最后,采用差分法及准线性化方法,得到了可以应用离散正交法求解的准线性微分方程组．对于导电圆柱薄壳,得到了Lorentz力表达式,并且推导了温度场积分特征值．讨论了导电圆柱薄壳应力、温度及变形随外加电磁参量的变化规律,并通过实例证实了可以通过改变电、磁、力场的参数来实现对薄壳的应力、应变、温度的控制．     

Some general properties of Eshelby’s tensor fields in transport phenomena and anti-plane elasticity

Consider an infinite thermally conductive medium characterized by Fourier’s law, in which a subdomain, called an inclusion, is subjected to a prescribed uniform heat flux-free temperature gradient. The second-order tensor field relating the gradient of the resulting temperature field over the medium to the uniform heat flux-free temperature gradient is referred to as Eshelby’s tensor field for conduction. The present work aims at deriving the general properties of Eshelby’s tensor field for conduction. It is found that: (i) the trace of Eshelby’s tensor field is equal to the characteristic function of the inclusion, independently of the latter’s shape; (ii) the isotropic part of Eshelby’s tensor field over the inclusion of arbitrary shape is identical to Eshelby’s tensor field over a 2D circular or 3D spherical inclusion; (iii) when the medium is made of an isotropic material and when the inclusion has some specific rotational symmetries, the value of the Eshelby’s tensor field evaluated at the inclusion gravity center and the symmetric average of Eshelby’s tensor fields are both equal to Eshelby’s tensor field for a 2D circular or 3D spherical inclusion. These results are then extended, with the help of a linear transformation, to the general case where the medium consists of an anisotropic conductive material. The method elaborated and results obtained by the present work are directly transposable to the physically analogous transport phenomena of electric conduction, dielectrics, magnetism, diffusion and flow in porous media and to the mathematically identical phenomenon of anti-plane elasticity.     

Magnetic and Gravitational Convection of Air in a Porous Cubic Enclosure with a Coil Inclined Around the Y Axis

The natural convection heat transfer of air in a porous media can be controlled by gradient magnetic field. Thermomagnetic convection of air in a porous cubic enclosure with an electric coil inclined around the $Y$ axis was numerically investigated. The Biot–Savart law was used to calculate the magnetic field. The governing equations in primitive variables were discretized by the finite-volume method and solved by the SIMPLE algorithm. The flow and temperature fields for the air natural convection were presented and the mean Nusselt number on the hot wall was calculated and compared. The results show that both the magnetic force and coil inclination have significant effect on the flow field and heat transfer in a porous cubic enclosure, the natural convection heat transfer of air can be enhanced or controlled by applying gradient magnetic field.     

磁-电-弹性半空间在轴对称热载荷作用下的三维问题研究

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.      

THREE-DIMENSIONAL ANALYSIS OF A MAGNETOELECTROELASTIC HALF-SPACE UNDER AXISYMMETRIC TEMPERATURE LOADING 1)

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.     

MAGNETO-THERMO-ELASTIC COUPLED DYNAMICS EQUATIONS OF AXIALLY MOVING CARRY CURRENT PLATE IN MAGNETIC FIELD

针对磁场环境中轴向运动载电流导电板磁热弹性耦合动力学建模问题进行研究. 考虑几何非线性和热效应条件下, 给出薄板运动的动能、应变能以及外力虚功的表达式.应用哈密顿变分原 理, 推得力、运动、电、磁和热效应相互作用下轴向运动导电板的非线性磁热弹性耦合振动方程.基于麦克斯韦电磁场方程, 考虑相应的电磁本构关系和电磁边界条件, 推得磁场环境中轴向运动载电流导电板满足的电动力学方程和所受电磁力表达式, 并给出焦耳热作用下耦合形式的热传导方程. 算例表明, 磁场等参量对动力学系统分岔特性有明显影响.所得结果可为此类问题的进一步求解和分 析提供理论参考.     

A microstructure continuum approach to electromagneto-elastic conductors

A micromorphic continuum model of a deformable electromagnetic conductor is established introducing microdensities of bound and free charges. The conductive part of electric current consists of contributions due to free charges and microdeformation. Beside the conservation of charge, we derive suitable evolution equations for electric multipoles which are exploited to obtain the macroscopic form of Maxwell’s equations. A constitutive model for electromagneto-elastic conductors is considered which allows for a natural characterization of perfect conductors independently on the form of the constitutive equation for the conduction current. A generalized Ohm’s law is also derived for not ideal conductors which accounts for relaxation effects. The consequences of the linearized Ohm’s law on the classic magnetic transport equation are shown.     

Self-similar magnetogasdynamic flows accompanied by detonation and combustion waves

Magnetogasdynamic (MGD) flows with detonation waves and combustion fronts have attracted more and more attention in recent years. Intensive heat supply assures such a significant increase in the temperature and pressure behind the heat liberation fronts that the gaseous combustion products become conductive so that the flow map in the electric and magnetic fields can vary substantially as compared with ordinary gasdynamics. In the case of finite gas conductivity, when the magnetic Reynolds numbers R_m are low, the asymptotic laws of detonation wave propagation which either go over into the Chapman-Jouguet (CJ) mode (in a number of cases at a finite distance from the initiation source) or remain overcompressed, have been studied [1]. Stationary flow modes behind detonation waves have been investigated in [2] and the problem of the detonation wave originating at the closed end of the tube emerging in the stationary mode in crossed homogeneous magnetic and electric fields has been examined. Results are presented in this paper of an investigation of one-dimensional self-similar flows caused by piston motion in a hot gas mixture in which a detonation wave or combustion front is propagated. The motion is realized in external electric and magnetic fields which exert a substantial effect on the flow of the conductive combustion products. Domains of application of the governing parameters in which the various flow modes are realized are found by using a qualitative and numerical analysis. The results obtained are used to solve problems about the hypersonic gas flow around a thin wedge in an axial magnetic field.     

压电、压磁材料球对称问题的通解

研究压电、压磁材料在球坐标系下,不计体力、体电荷和体电流的情况下,由平衡方程、梯度方程、压电和压磁的本构方程导出应力、应变、位移、电位移、电场强度、电位势、磁感强度和磁位势各未知量的通解．考虑不同的边界条件,将其通解应用到应力、电学短路以及位移、电学开路和磁场分布的边界条件中,得到不同边界条件下问题的解。     

孔隙地层中声波诱导的电磁场测井响应模拟

声电效应测井过程中,当充流体井中发射的声波传播至井外时,在孔隙地层中引起含带电离 子的孔隙流体相对于固体骨架流动,导致电磁场. 在忽略转换电磁场对声场影响的条件下,先采用Biot-Rosenbaum模型计算孔隙地层中的 声场,再计算声波诱导的电场和磁场. 计算结果表明,有伴随纵波波群和横波波群的电磁场、有伴随斯通利波的电磁场,还有微弱的 临界折射电磁波. 计算还表明: 接收器偏离井轴越远,电场和磁场的强度越大; 在井壁处,磁场强度的波形与轴向电场强度的波形相似.     

柔性微粒介电泳分离过程的多尺度模拟

介电泳分离是一种高效的微细颗粒分离技术,利用非均匀电场极化并操纵分离微流道中的颗粒. 柔性微粒在介电泳分离过程中同时受多种物理场、多相流和微粒变形等复杂因素的影响,仅用单一的计算方法对其进行模拟存在一定的难度,本文采用有限单元——格子玻尔兹曼耦合计算的方法处理这一难题.介观尺度的格子玻尔兹曼方法将流体看成由大量微小粒子组成,在离散格子上求解玻尔兹曼输运方程,易于处理多相流及大变形问题,特别适合模拟柔性颗粒在介电泳分离过程中的变形情况.另一方面,介电泳分离过程的模拟需求解流体、电场和微粒运动方程,计算量相当庞大,通过有限单元法求解介电泳力,提高计算效率.利用这种多尺度耦合计算方法,对一款现有的介电泳芯片分离过程进行了模拟.分析了微粒在电场作用下产生的介电泳力,揭示了介电泳力与电场变化率等因素之间的关系.对微粒运动轨迹及其变形的情况进行了研究,发现微粒的变形主要与流体剪切作用有关.这种多尺度耦合计算方法,为复杂微流体的计算提供了一种有效的解决方案.      

Multi-frequency excitation of magnetorheological elastomer-based sandwich beam with conductive skins

The present work deals with the dynamic stability of a symmetric sandwich beam with magnetorheological elastomer (MRE) embedded viscoelastic core and conductive skins subjected to time varying axial force and magnetic field. The conductive skins induce magnetic loads and moments under the application of magnetic field during vibration. The MRE part works in shear mode and hence the dynamic properties of the sandwich beam can be controlled by magnetic fields due to the field dependent shear modulus of MRE material. Considering the core to be incompressible in transverse direction, classical sandwich beam theory has been used along with extended Hamilton's principle and Galarkin's method to derive the governing equation of motion. The resulting equation reduces to that of a multi-frequency parametrically excited system. Second order method of multiple scales has been used to study the stability of the system for simply supported and clamped free sandwich beams. Here the experimentally obtained properties of magnetorheological elastomers based on natural rubber have been considered in the numerical simulation. The results suggest that the stability of the MRE embedded sandwich beam can be improved by using magnetic field.     

Thomson effect with hyperbolic two-temperature on magneto-thermo-visco-elasticity

The study considers a homogeneous isotropic thermo-visco-elastic solid with hyperbolic two-temperature to cope up with its two-dimensional(2 D) deformations. The heat conduction equation is influenced by the Thomson coefficient. Lord-Shulman's theory is used to modify the basic governing equations. A method called "normal mode analysis" is utilized to attain the magnetic field, stress, conductive and thermodynamic temperature, and displacement components. Also, a number of numerical calculations are performed and discussed to understand the impact of hyperbolic two-temperatures,Thomson parameter, and viscosity on the material mentioned above.     

基于Hamilton体系的辛半解析法在各向异性电磁波导中的应用

力学中的Hamilton体系使用对偶变量来描述问题,而电磁场正好有电场和磁场这一对对偶变量。本文将力学中的Hamilton体系应用到电磁波导问题。根据电磁波导的Hamilton体系理论,辛几何可用于任意各向异性材料。将横向的电场和磁场构成对偶向量,基于Hamilton变分原理做半解析横向离散,并保持结构辛体系。本文以各向异性材料电磁波导为例,求解了问题的辛本征值,得到了镜像线的色散曲线。     

Heating of a medium as a result of Joule energy dissipation

It is well known that when a magnetic field is present and electric currents flow through a gas, terms over and above those present in the case of ordinary thermal conductivity appear in the heat flux density vector. If the gas is dense enough and the magnetic field not over large, then the anisotropy caused by the magnetic field may be neglected. However, for a sufficiently large electric current a term proportional to the temperature and to the current density vector remains in the heat flux density vector. This effect explains, for example, the asymmetry of heat fluxes in the electrodes of a continuously operated electromagnetic accelerator (1).We shall consider this situation in relation to the example of a fully ionized quasineutral gas with identical electron and ion temperatures.     

Modelling turbulent and dispersion heat fluxes in turbulent porous medium flows using the resolved LES data

Resolved large eddy simulations (LESs) of turbulent conjugate heat transfer in porous media are performed by the lattice Boltzmann method (LBM) for modelling turbulent and dispersion heat flux terms of the double-averaged energy equation. The considered porous structures are square rod arrays, staggered cube arrays and body centred cubic foam. In the LBM, the double-distribution function method which solves the distribution functions for the velocity and the internal energy is used. For the velocity and thermal fields, the D3Q27 multiple-relaxation-time method and the regularized D3Q19 single-relaxation-time method are applied, respectively. A priori tests using the LES data suggest that the trends of the sum of the dispersion and volume-averaged turbulent heat fluxes can be well captured by the second order gradient diffusion model.     

Fracture mechanics for a mode III crack in a magnetoelectroelastic solid

In this paper, we developed a Stroh-type formalism for anti-plane deformation and then investigated the fracture mechanics for an elliptical cavity in a magnetoelectroelastic solid under remotely uniform in-plane electromagnetic and/or anti-plane mechanical loading, which allowed us to take the electromagnetic field inside the cavity into account. Reducing the cavity into a crack, we had explicit solutions in closed forms for a mode III crack, which included the extreme cases for an impermeable crack and a permeable crack. The results were illustrated with plots, showing that in the absence of mechanical loads, an applied electric or magnetic field, positive or negative, always tended to close the crack. On the other hand, in the presence of a mechanical load, a negative electric or magnetic field retarded crack growth, while a positive field could either enhance or retard crack propagation, depending on the strengths of the applied electric/magnetic fields and the level of the mechanical load as well. In other words, the effect of electric/magnetic fields on the fracture behavior is mechanical load-dependent.     

Exact Electromagnetothermoelastic Solution for a Transversely Isotropic Piezoelectric Hollow Sphere Subjected to Arbitrary Thermal Shock

This paper presents analytical study for electromagnetothermoelastic transient behavior of a transversely isotropic hollow sphere, placed in a uniform magnetic field, subjected to arbitrary thermal shock. Exact solutions for the transient responses of stresses, perturbation of magnetic field vector, electric displacement and electric potential in the transversely isotropic piezoelectric hollow sphere are obtained by means of the Hankel transform, the Laplace transform and their inverse transforms. An interpolation method is used to solve the Volterra integral equation of the second kind caused by interactions among electric, magnetic, thermal and elastic fields. From the sample numerical calculations, it is seen that the present method is suitable for the transversely isotropic hollow sphere, placed in a uniform magnetic field, subjected to arbitrary thermal shock. Finally, the result can be used as a reference to solve other transient coupling problems of electromagnetothermoelasticity.