带铁磁薄膜悬臂板的磁场微感应器磁弹性特征研究

对于在可变形非磁材料悬臂梁式板单表面粘贴可磁化材料薄膜所构成的磁场微传感器件结构,研究了其处在磁场环境中的磁弹性弯曲变形的磁场-力学特征。为此,建立了由有限元方法分析磁场与有限差分法计算挠曲变形相结合来计算其结构在磁场中产生磁弹性变形的定量分析程序。在此基础上,对于这一微传感结构的算例给出了其结构变形随外加磁场环境变化的磁场-挠度特征关系等定量结果。结果表明：微传感器件不仅可以测量出磁场的大小,而且给出了测量磁场矢量方向的可能性。     

A theoretical model of magnetoelastic bucking for soft ferromagnetic thin plates

As an essential model of magnetoelastic interaction between magnetic field and mechanical deformation, the study on magnetoelastic buckling phenomenon of soft ferromagnetic plates in a magnetic environment has been conducted. One of the key steps for the theoretical prediction of the critical magnetic field is how to formulate magnetic force exerted on the magnetized medium. Till today, the theoretical predictions, from theoretical models in publications, of the magnetoelastic buckling of ferromagnetic cantilevered beam-plate in transverse magnetic field are all higher than their experimental data. Sometimes, the discrepancy between them is as high as 100%. In this paper, the macroscope formulation of the magnetic forces is strictly obtained from the microscope Amperion current model. After that, a new theoretical model is established to describe the magnetoelastic buckling phenomenon of ferromagnetic thin plates with geometrically nonlinear deformation in a nonuniform transverse magnetic field. The numerical method for quantitative analysis is employed by combining the finite elemental method for magnetic fields and the finite difference method for deformation of plates. The numerical results obtained from this new theoretical model show that the theoretical predictions of critical values of the buckling magnetic field for the ferromagnetic cantilevered beam-plate are in excellent agreement with their experimental data. By the way, the region of applicability to the Moon-Pao's model, or the couple model, is checked by quantitative results. This project was supported in part by the National Natural Science Foundation of China and the Foundation of the SEdC of China for Returned Chinese Scholars from Abroad.     

Magnetoelastic bending and stability of soft ferromagnetic rectangular plates

I.IntroductionSofar,twokindsofphenomenaofmagnetoelasticinteractiontoferromagneticplatesinappliedmagneticfieldshavebeenfound.Oneisthemagnetoelasticbuchlingorillstabllitytoacantileveredbeam-plateinatransversemagneticfieldt'],thatisthephen1Onlenollot'negiltivemagneticstiffness,whiletheotheristheincreaseofnaturalfrequencyoftheplatewhichisfreelyviratinginanin-planemagneticfieldl=],thatisthephenomenonofpositivemagneticstiffness.Inordertosimulatethesetwokindsofexperimentalphenomena,sometheoreticalmo…     

软铁磁薄板磁弹性屈曲的理论模型

铁磁弹性薄板的磁弹性屈曲问题一直作为电磁——弹性力学相互作用的一个基本模型进行研究，而作用在其磁介质上的磁力计算则是定量理论预测准确与否的关键．到目前为止，文献上已有的理论模型对悬臂铁磁梁式悬臂板在横向磁场中磁弹性屈曲的理论预测值始终高于实验值，有的甚至相差１００％左右．本文基于电磁力计算的微观安培电流模型，严格给出了软铁磁薄板等效横向磁力的宏观计算表达式．在此基础上，建立了电磁——力学相互耦合作用的非线性理论模型．该模型能描述铁磁薄板结构在非均匀横向磁场环境中的磁弹性失稳（或屈曲）．其定量分析采用了有限元法和有限差分法相结合．数值结果显示：本模型给出的磁弹性屈曲的临界磁场值与实验值符合良好．与此同时，文中还对文献中认为较成功的Ｍｏｏｎ－Ｐａｏ模型的基本假设进行了分析．定量结果发现：Ｍｏｏｎ－Ｐａｏ理论模型的基本假设仅在梁式板的长厚比Ｌ／ｈ比较大时（约在２００左右），是可以接受的，而当Ｌ／ｈ较小时，该假设将导致理论值与实验值的较大误差．Ｌ／ｈ比值越小，理论值与实验值的误差越大     

Magnetoelastic bending and stability of current-carrdying coil structures

As a results of magnetoelastic interaction, the mechanical behavior of current-carrying coil structures, such as deformation and instability, is a key problem in the design of strong field magnets. In this paper, a nonlinear mathematical model is presented to described the deformation and buckling of D-type current-carrying coil, based on the Biot-Savart law and the bending theory of curved beams. The bending deformation, the critical value of current for the magnetoelastic buckling of the current-carrying coil, and the effect of the type and number of supports at middle part of the bendling coil on the critical value are quantitatively investigated by a semi-analytical and semi-numerical method. The numerical results are shown to be in good agreement with the experimental data. The project supported by the National Natural Science Foundation of China, the Science Foundation of the State Education Committee of China for Outstanding Teacher in Universities, and the Natural Science Foundation of Gansu Province of China     

METHOD OF GREEN＇S FUNCTION OF CORRUGATED SHELLS

By using the fundamental equations of axisymmetric shallow shells of revolution, the nonlinear bending of a shallow corrugated shell with taper under arbitrary load has been investigated. The nonlinear boundary value problem of the corrugated shell was reduced to the nonlinear integral equations by using the method of Green＇s function. To solve the integral equations, expansion method was used to obtain Green＇s function. Then the integral equations were reduced to the form with degenerate core by expanding Green＇s function as series of characteristic function. Therefore, the integral equations become nonlinear algebraic equations. Newton＇ s iterative method was utilized to solve the nonlinear algebraic equations. To guarantee the convergence of the iterative method, deflection at center was taken as control parameter. Corresponding loads were obtained by increasing deflection one by one. As a numerical example,elastic characteristic of shallow corrugated shells with spherical taper was studied.Calculation results show that characteristic of corrugated shells changes remarkably. The snapping instability which is analogous to shallow spherical shells occurs with increasing load if the taper is relatively large. The solution is close to the experimental results.     

Study on magneto–elastic–plastic deformation characteristics of ferromagnetic rectangular plate with simple supports

In this paper, the magnetic-elastic-plastic deformation behavior is studied for a ferromagnetic plate with simple supports. The perturbation formula of magnetic force is first derived based on the perturbation technique, and is then applied to the analysis of deformation characteristics with emphasis laid on the analyses of modes, symmetry of deformation and influences of incident angle of applied magnetic field on the plate deformation. The theoretical analyses offer explanations why the configuration offer- romagnetic rectangular plate with simple supports under an oblique magnetic field is in-wavy type along the x-direction, and why the largest deformation of the ferromagnetic plate occurs at the incident angle of 45°for the magnetic field. A numerical code based on the finite element method is developed to simulate quantitatively behaviors of the nonlinearly coupled multi-field problem. Some characteristic curves are plotted to illustrate the magneto--elastic-plastic deflections, and to reveal how the deflections can be influenced by the incident angle of applied magnetic field. The deformation characteristics obtained from the numerical simulations are found in good agreement with the theoretical analyses.     

A 1D time-domain method for in-plane wave motions in a layered half-space

A 1D finite element method in time domain is developed in this paper and applied to calculate in-plane wave motions of free field exited by SV or P wave oblique incidence in an elastic layered half-space. First, the layered half-space is discretized on the basis of the propagation characteristic of elastic wave according to the Snell law. Then, the finite element method with lumped mass and the central difference method are incorporated to establish 2D wave motion equations, which can be transformed into 1D equations by discretization principle and explicit finite element method. By solving the 1D equations, the displacements of nodes in any vertical line can be obtained, and the wave motions in layered half-space are finally determined based on the characteristic of traveling wave. Both the theoretical analysis and the numerical results demonstrate that the proposed method has high accuracy and good stability. The project supported by the National Natural Science Foundation of China (50478014), the National 973 Program (2007CB714200) and the Beijing Natural Science Foundation (8061003). The English text was polished by Yunming Chen.     

On multiple circular inclusions in plane magnetoelasticity

A general series solution to the magnetoelastic problem of interacting circular inclusions in plane magnetoelasticity is provided in this paper. By the use of complex variable theory and Laurent series expansion method, the general expression of the magnetic and the magnetoelastic complex potentials for the circular inclusion problem is derived. Expanding the definition of the Airy’s stress function of pure elastic field into the magnetoelastic field and applying the superposition method, the general expression then can be reduced to a set of linear algebraic equations and solved in a series form. An approximate closed form solution for the case of two arbitrarily located inclusions is also provided. For illustrating the effect of the pertinent parameters, the numerical results of the interfacial magnetoelastic stresses are displayed in graphic form.     

A SOLUTION PROCEDURE FOR LOWER BOUND LIMIT AND SHAKEDOWN ANALYSIS BY SGBEM

The symmetric Galerkin boundary element method (SGBEM) instead of the finite element method is used to perform lower bound limit and shakedown analysis of structures. The self-equilibrium stress fields are constructed by a linear combination of several basic self-equilibrium stress fields with parameters to be determined. These basic self-equilibrium stress fields are expressed as elastic responses of the body to imposed permanent strains and obtained through elastic-plastic incremental analysis. The complex method is used to solve nonlinear programming and determine the maximal load amplifier. The limit analysis is treated as a special case of shakedown analysis in which only the proportional loading is considered. The numerical results show that SGBEM is efficient and accurate for solving limit and shakedown analysis problems. Project supported by the National Natural Science Foundation of China (No. 19902007), the National Foundation for Excellent Doctorial Dissertation of China (No. 200025) and the Basic Research Foundation of Tsinghua University.     

Incremental Magnetoelastic Deformations,with Application to Surface Instability

In this paper the equations governing the deformations of infinitesimal (incremental) disturbances superimposed on finite static deformation fields involving magnetic and elastic interactions are presented. The coupling between the equations of mechanical equilibrium and Maxwell’s equations complicates the incremental formulation and particular attention is therefore paid to the derivation of the incremental equations, of the tensors of magnetoelastic moduli and of the incremental boundary conditions at a magnetoelastic/vacuum interface. The problem of surface stability for a solid half-space under plane strain with a magnetic field normal to its surface is used to illustrate the general results. The analysis involved leads to the simultaneous resolution of a bicubic and vanishing of a 7×7 determinant. In order to provide specific demonstration of the effect of the magnetic field, the material model is specialized to that of a “magnetoelastic Mooney–Rivlin solid”. Depending on the magnitudes of the magnetic field and the magnetoelastic coupling parameters, this shows that the half-space may become either more stable or less stable than in the absence of a magnetic field.      

Changes in the natural frequency of a ferromagnetic rod in a magnetic field due to magnetoelastic interaction

Based on the magnetoelastic generalized variational principle and Hamilton＇s principle, a dynamic theoretical model characterizing the magnetoelastic interaction of a soft ferromagnetic medium in an applied magnetic field is developed in this paper. From the variational manipulation of magnetic scale potential and elastic displacement, all the fundamental equations for the magnetic field and mechanical deformation, as well as the magnetic body force and magnetic traction for describing magnetoelastic interaction are derived. The theoretical model is applied to a ferromagnetic rod vibrating in an applied magnetic field using a perturbation technique and the Galerkin method. The results show that the magnetic field will change the natural frequencies of the ferromagnetic rod by causing a decrease with the bending motion along the applied magnetic field where the magnetoelastic buckling will take place, and by causing an increase when the bending motion of the rod is perpendicular to the field. The prediction by the mode presented in this paper qualitatively agrees with the natural frequency changes of the ferromagnetic rod observed in the experiment.     

Interval Arithmetic and Static Interval Finite Element Method

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｔｈｅａｎａｌｙｓｉｓａｎｄｄｅｓｉｇｎｏｆｓｔｒｕｃｔｕｒｅｓ,ｓｏｍｅｕｎａｖｏｉｄａｂｌｅｕｎｃｅｒｔａｉｎｔｉｅｓ ,ｓｕｃｈａｓｔｈａｔｏｆｍａｔｅｒｉａｌａｎｄｇｅｏｍｅｔｒｉｃａｌｐｒｏｐｅｒｔｉｅｓ,ｌｏａｄｓ ,ａｎｄｓｏｏｎ ,ｓｈｏｕｌｄｂｅｒｅａｓｏｎａｂｌｙｔａｋｅｎｉｎｔｏａｃｃｏｕｎｔ.Ｉｎｔｈｅｐａｓｔｄｅｃａｄｅｓ,ｔｈｅｓｅｕｎｃｅｒｔａｉｎｔｉｅｓｗｅｒｅｍｏｓｔｌｙｔｒｅａｔｅｄｗｉｔｈｐｒｏｂａｂｉｌｉｔｙｔｈｅｏｒｙｏｒｒａｎｄｏｍｐ…     

A finite deformation analysis of polycrystal by considering the influence of grain boundary

By combining grain boundary (GB) and its influence zone, a micromechanic model for polycrystal is established for considering the influence of GB. By using the crystal plasticity theory and the finite element method for finite deformation, numerical simulation is carried out by the model. Calculated results display the microscopic characteristic of deformation fields of grains and are in qualitative agreement with experimental results.The project is supported by National Natural Science Foundation of China.     

悬臂铁磁板磁弹性耦合作用的力学分析

对有限板宽的悬臂铁磁板在外磁场中的磁弹性相互耦合作用的力学行为，建立了描述板弯曲和稳定性的理论模型及有限元定量分析程序，研究了外磁场倾斜角对板磁弹性失稳的临界磁场值的影响。     

考虑尺寸效应板梁的磁弹性屈曲

研究表明,有限尺寸效应是影响板——梁磁弹性稳定性分析精度的主要原因。本文用边界元法分析了板边沿的磁感强度集中,然后用里兹法计算了软铁磁材料悬臂板——梁的临界屈曲磁场。通过和Moon-Pao的理论和实验结果以及Miga的有限元结果比较,表明本方法具有良好的精度,且计算方法简便。     

ANALYSIS ON THE MAGNETO-ELASTIC-PLASTIC BUCKLING/SNAPPING OF CANTILEVER RECTANGULAR FERROMAGNETIC PLATES

An analysis of buckling/snapping and bending behaviors of magneto-elastic-plastic interaction and coupling for cantilever rectangular soft ferromagnetic plates is presented.Based on the expression of magnetic force from the variational principle of ferromagnetic plates,the buckling and bending theory of thin plates,the Mises yield criterion and the increment theory for plastic deformation,we establish a numerical code to quantitatively simulate the behaviors of the nonlinearly multi-fields coupling problems by the finite element method.Along with the phenom- ena of buckling/snapping and bending,or the characteristic curve of deflection versus magnitude of applied magnetic fields being numerically displayed,the critical loads of buckling/snapping, and the influences of plastic deformation and the width of plate on these critical loads,the plastic regions expanding with the magnitude of applied magnetic field,as well as the evolvement of deflection configuration of the plate are numerically obtained in a case study.     

Experimental and theoretical study of the nonlinear response of agiant magnetostrictive rod

In this paper, a magnetomechanical coupling constitutive relation of the giant magnetostrictive material was investigated experimentally and theoretically. A grain-oriented magnetostrictive rod of iron and rare earth was tested under a combined magnetomechanical loading. Two types of experimental curves were obtained, i.e., the magnetostrictive curve of the extensional strain vs the magnetic field, and the curve of the magnetic polarization intensity vs the pre-stress. A new theoretical constitutive model, based on the density of domain switching, is developed. Comparison of the theoretical predictions with the experimental results indicates that this model can capture the main characteristics of the magnetoelastic coupling deformation of a giant magnetostrictive rod. The project supported by the National Natural Science Foundation of China (10025209, 10132010, 10102007)     

An analysis for imperfect sensitivity of a superconducting rod

On the basis of the interaction between the electromagnetic field and a solid, the imperfect sensitivity of a superconducting rod is studied in this paper. The influences of the initial deflection of the rod and the initial magnetic field on the stability are discussed. One can see that the bifurcation response of the rod with imperfections is similar to the universal unfolding of the fork bifurcation from the singularity theory.Project supported by the Foundation of the State Education Commission of China and the Natural Science Foundation of Gansu Province.     

横向磁场中轴向运动铁磁梁的磁弹性主共振

研究在外激励力与磁场作用下轴向运动铁磁梁的磁弹性非线性主共振问题．基于弹性理论和电磁理论，给出梁的动能和弹性势能表达式，根据哈密顿原理，推导出磁场中轴向运动铁磁梁的磁弹性双向耦合非线性振动方程．通过伽辽金积分法进行离散，得出两端简支边界条件下铁磁梁磁弹性非线性强迫振动方程．应用多尺度法对方程进行求解，得出幅频响应方程．最后通过算例，给出铁磁梁的幅频特性曲线、振幅-磁感应强度和振幅-外激励力曲线并进行分析．结果显示，在幅频响应曲线中铁磁梁的轴向运动速度、外激励力、轴向拉力越大，共振振幅越大；而磁感应强度越大，振幅越小．