1.Introduction

Tubular moving-magnet linear oscillating motors(TMMLOMs)are specific linear actuators that produce a high-frequency reciprocating motion.Their outstanding characteristics,such as efficiency,dynamic performance,and simple structure,give their broad application prospect,but their model has not been clearly described and their structure has not been optimized adequately.TMMLOMs work at a designed frequency with the stroke of mover oscillating no more than one pole pitch,which are widely applied in many equipment,such as artificial hearts,compressors,refrigerators,etc.1–5As short-travel actuations are propelled directly without conventional cranks,TMMLOMs enjoy a 20%-30%better efficiency compared with rotary motors in the application of household refrigerator compressors.5By enlarging the power whilst reducing acoustic noise and vibrations accordingly,TMMLOMs are tested for air-conditioner compressors with an efficiency of over 92%at a rated condition.6Due to a relatively higher frequency response benefited from the absence of rotor inertia in the rotary mechanism,TMMLOMs are developed to drive the servo pump in a hydraulic system.7Meanwhile,TMMLOMs are utilized in linear electrohydrostatic actuators in more-electric aircraft systems with a high power density.8In addition,a cooperative con figuration is introduced where two movers of dual TMMLOMs work with a phase angle of π/2.9A reliable coordinated control of two TMMLOMs requires a regulation based on an accurate model.Furthermore,a powerful output with limited space and restricted mass in aircraft also needs more optimized analysis of the motor structure and principle based on a precise calculation of the air-gap magnetic field.A comprehensive model,therefore,is indispensable for TMMLOMs.

Some well-known modeling optimization approaches,such as changing the direction of magnetization on a mover,con figuring the groups of PM arrays on movers,decreasing the slots of the stator,and regulating the split ratio of E-Core PMs,have been con firmed to improve the characteristics of TMMLOMs.However,current research pays little attention to the effects of mover gaps(air gaps among PMs or between PMs and materials)which exist inherently in the segments of PMs.Despite the fact that the effects of discontinuity on PMs has not been proven to be very prominent,an introduction of their effects on the distribution of a magnetic field through an analytical mathematics and simulation model is also a necessity when aiming to optimize a TMMLOM for a high-performance aeronautical facility and create a dual motor collaborative control strategy on this basis.

Many studies about the TMMLOMs have been published over the last three decades,focusing on improving the output performance analytically and experimentally through reducing the mover weight and increase the frequency.3,4

引导基金按照政府引导、市场运作、科学决策、防范风险的原则运作，主要采取阶段参股模式，即与投资机构、企业集团等社会资本共同组建子基金，由子基金对企业进行股权投资。引导基金出资一般不超过子基金规模的25%，天使基金不超过50%，子基金由专业投资机构管理，以市场化方式来运作。引导基金实现约定目标从子基金退出时，可将其出资部分的预期收益让渡给社会资本，其中天使基金收益返还比例不超过40%，其他基金不超过30%。

The topology of PMs is an important factor that influences the output performance.Kim et al.proposed a linear oscillatory actuator whose magnetization direction of PMs was parallel to the stroke axis,10and the machine shared a high-power density and a low-cogging force.A transient model of this motor was validated experimentally.11In the contrary,Kim et al.discussed a TMMLOM with its magnetization direction perpendicular to the motion direction.12,13The magnetic field was analyzed by an analytical model,and the prototype gave capabilities of high force density and low detent force.By the ways of magnetic vector potential and cylindrical coordinate formulation,Wang et al.presented a general framework for TMMLOMs.His solutions presented the analytical field distribution for both axial polarized and radial polarized tubular linear PM motors.14The air-gap flux densities of both topologies were compared,15and a finite element model showed that the radial attractive force was higher in the axial polarized structures when considering the eccentricity effect.

In a publication series presented by Wang et al.,a quasi-Halbach structure of TMMLOMs which was a combination of radial polarized and axial polarized PMs arrays was deeply analyzed,14,16–19including an analytical flux expression of the quasi-Halbach topology,14,16,18comparative analysis to conventional magnetization methods,17and parameters optimizations on the size of the structure.19Because of some selfshielding property,the quasi-Halbach PM arrays would reduce the moving mass to improve the dynamic capability and result in a higher flux density.

The coupling between windings and PMs is another area that researchers have focused on.Considering the con figuration of stator windings in TMMLOMs,Wang et al.deduced a set of formulations to describe the distribution of magnetic flux theoretically based on Maxwell’s equations.16Zhu et al.studied a variety of winding arrangements with a PMs assembly.20–24A rule was outlined that the stator tooth number was supposed to differ the mover pole number by 1 in TMMLOMs.20Subsequently,simulation calculations and experiment validations were finished by Zhu et al.,which demonstrated that the E-core wind con figuration yielded less flux leakage and performed conducive to the oscillation.20,22Eventually,design optimization and prototype validation was implemented for the topology of a TMMLOM with E-type windings and quasi-Halbach PMs.23,24This type of motor was further extended into E-type series TMMLOMs,25which consisted of multi-pair E-cores with a quasi-Halbach topology to enhance the thrust.Furthermore,Jiao et al.presented compound Halbach PMs in a TMMLOM,whose topology was a dual-layer integration of quasi-Halbach PMs and axial polarized PMs.26In Jiao’s structure,the conventional back-iron was not needed any more,hence the mover mass could be decreased further with a better dynamic performance.

(1).Region I is the space of the air-gap.

By introducing a Carter coefficient,Wang et al.modified the analytical model of tubular linear PM motors with a slotted stator14and a semi-closed slotted stator.27They also applied this correction to calculate the distribution of flux density of a TMMLOM with quasi-Halbach PMs.This method was also followed by Wang et al.to model the E-type series TMMLOMs.28

The detailed formulas and the expression of Ln(x)are presented in Appendix B.According to Eqs.(1)and(7),the distribution of the magnet field is given as follows.

Experimental,numerical,and analytical methods are the three normal research approaches on LMMLOM research.Some studies have been implemented entirely experimentally,6,10,29,30whose results could exhibit the performance of a motor directly.However,they didn’t reveal the principle of LMMLOMs.Numerical approaches have been widely applied in LMMLOM designs,7,17,31–35which allow to solve a complex design without much approximation.However,it didn’t reflect the relation between topological parameters and the motor performance,which results in a difficulty in optimization.

where μ0is the permeability of vacuum with a value of 4π× 10-7H/m,μrIIthe relative permeability of PMs,and M0the residual magnetization vector in A/m.

An accurate model depends on the numbers of factors which should be taken into account during modeling.Most of the literature studied TMMLOMs,and took the con figurations of mover PMs and stator topology into consideration.The flux density in the air-gap determines the output capacity of the motor.A quasi-Halbach PMs array proves its better performance in the flux density due to the self-shielding property,and its model has been analyzed in detail.16–24,28It produces a larger flux distribution with an equal mover mass,which can increase properties of both thrust and dynamic.Reduction of the stator slot,in addition,is also an effective approach to improve the flux distribution,whose model has also drawn a lot of focus from researchers,7,14,15,17,23,26,27because permeability is quite sensitive to air medium.

Air medium also exists in the mover of a TMMLOM,the surface of which is assembled by several individual PM modules.The manufacturing method results in some gaps among PM arrays as well as between PMs and the mover frame.Bianchi analyzed that the magnetic flux could move through the small air-gap in an axial magnetized mover causing a radial force.15However,he just focused on the radial force in a TMMLOM by the method of finite element and didn’t give a model to explain how the gaps affected the flux distribution.Via analytical analysis,39Xia et al.presented the negative influence caused by gaps between Halbach PM segments in rotary motors,discussing the relationship between the effect of gaps and the parameters of PMs through an analytical model.In Xia’s paper,however,the gaps were only investigated in a circumferential dimension.Similar to rotary motors,mover gaps also exist in linear motors as well as TMMLOMs,whose gaps are more complex than those of rotary machines because of a more complex structure in the mover.The surface of the mover is mounted by PM segments including radial polarized ones and axial polarized ones.Due to the installation error,the stuffed glue,and the interaction force among PMs,gaps are inevitable between and beneath segments in quasi-Halbach arrays.The existence of gaps would affect the magnetic flux distribution and output thrust.Therefore,the thrust model of TMMLOMs considering mover gaps is a necessity for designing a high-precision control strategy.

In this paper,Section 1 introduces the problem statement of TMMLOMs.Section 2 establishes a modeling method considering the mover gaps.Section 3 gives model validation and performance analysis.Finally,conclusions are presented in Section 4.

2.Schematic structure

In this section,the structure of a TMMLOM is clarified for further understanding.Subsequently,the existence of gaps in the mover is presented in detail.

主要仪器与试剂：立式高压蒸汽灭菌锅为YXQ-LS，购自上海博讯实业有限公司；生物安全柜为ESCO Class II，购自珠海市珠峰仪器仪表有限公司；气质联用仪为6890N-5973MSD，购自安捷伦科技有限公司；二氯甲烷为ACS级，购自北京百灵威科技有限公司；4-辛醇为GC级，购自阿拉丁试剂有限公司。

2.1.Structure of a TMMLOM

Fig.1 Structure of a TMMLOM.

The structure of a TMMLOM consisting of 4 pairs of E-type windings with a quasi-Halbach PMs con figuration is shown in Fig.1.The Matierial of PMs is NdFeB N38SH.Silicon steel is laminated circumferentially to make up the stator core,and coils are connected in parallel through the stator,while the hollow mover is made of electrical iron with PMs covered on the surface.Including radial polarized PMs and axial polarized PMs,a quasi-Halbach array can enhance the magnetic flux over the outer surface of the mover resulting from its selfshielding effect.The stroke of this machine is designed to be 5 mm in maximum.

The hollow core lightens the mover mass,which enables the TMMLOM to work at a high frequency of 100 Hz.Moreover,this design also improves the dynamic performance.High power and good dynamic response make it possible to meet the requirements of an aircraft actuator.

2.2.Gap topology

As shown in Fig.2 and Fig.3,all the individual PM segments assembled on the mover core can be classi fied into two types.An axial polarized PM segment is a torus whose magnetization direction is along the z-axis.A radial polarized PM segment is a radial magnetized one,which was fabricated with a tile-like shape as well as a central angle of 60°

Gaps are very clear between these segments as marked in Fig.2.In the installation process of Quasi-Halbach arrays,the adhesive would take up less room inevitably between PMs.In addition,the design tolerance as well as the interaction force between segments will also account for the assembly gaps.For the permeability in non-ferromagnetic gaps is different from that in ferromagnetic ones,the existence of gaps will affect the distribution of the magnetic flux.Allowing for the structure of Quasi-Halbach segments,the gaps mentioned above can be divided into three categories,i.e.,radial gaps,axial gaps,and circumferential gaps,which are illustrated in Fig.3.x,y,z are three axes in Cartesian coordinate system and r,θ,z are axes in cylindrical coordinate system,where era,eci,eaxare the unit vectors in radial,circumferential and axial directions respectively.Radial gaps exist between PMs and the hollow core,axial gaps between axial polarized PMs and radial polarized PMs,and circumferential gaps between tile-like PMs.The three kinds of gaps have individual effects on the magnet flux,and subsequently,the characteristics of radial and axial gaps are about to be discussed in the following section.

Fig.2 Mover prototype.

Fig.3 Mover structure.

3.Mathematical model considering mover gaps

In this section,an analytical model is yielded in cylindrical coordinates by the method of magnetization vector.Based on Laplace’s and Poisson’s equations,radial and axial gaps are taken into consideration.Moreover,boundary conditions are presented and the distribution of the magnetic flux is deduced in the air-gap.On this basis,an output thrust model is given.

3.1.Model assumption and magnetization vector

To research the effect of gaps in the mover on the air-gap magnetic field,it’s necessary to build an analytical model considering the gaps,and in order to simplify the derivation of the analytical model,some assumptions are stated here.

(1).The axial length of the theoretical model is in finite.It comprises a series of repetitive permanent magnet armatures in the TMMLOM extending to infinity,which facilitates a Fourier series representation,and the cycle length is much longer than the length of the motor.

在我国最新发布的护理事业发展规划（2016-2020）中指出，要将护理服务内涵与群众健康需求密切对接起来，提高护理专业的素质以及服务能力。这就要求：当前护理专业教学中应以医院临床的护理需求为导向，提高护理专业人员的岗位胜任力，切实提高我国护理行业的能力和发展水平。但是，在当前我国高职院校护理专业教学中还存在着护理教学理念严重落后于时代，护理教师信息化操作能力不能满足教育信息化发展的需求以及护理专业教学与医院护理岗位工作需求脱节等严重问题。

(2).Each individual permanent magnet shares the same size and magnetic property with a smooth surface,and the gaps of the same type have the same thickness.

(3).The analytical model ignores the saturation phenomenon.

(4).A Carter coefficient will be introduced to compensate the impact of the stator slot.

(5).The stator and the hollow core are in finitely permeable while the relative permeability of PMs is approximated as 1.

In formulation of the magnetic field,the room to be studied is divided into the following three regions based on magnetic characteristics,as shown in Fig.4.

To improve the properties of a TMMLOM,a stator slot was taken into account by researchers.Bianchi et al.analyzed the factor of a stator slot that affected the force.15His research showed that slotless motors shared a higher mean force than slotted ones.Kim et al.designed a slotted TMMLOM13and a slotless TMMLOM12for an eco-pedal system,and the later presented a smaller force ripple experimentally.An analytical model with a quasi-Halbach magnetized armature and a semi-closed slot stator was given by Chen et al.,23and a closed slot structure was employed into a TMMLOM to reduce the cogging force in Liang’s paper,7of which the leakage permeability was much larger than that of the former.Besides,Kim et al.proposed a novel method for stator lamination to laminate the teeth and yoke of a stator respectively,which decreases the gaps in the stator and multiplies the flux density.

(2).Region II denotes the permanent magnets arrays filled with a rare-earth magnetic material with axial gaps.

(3).Region III represents the radial gaps between the permanent magnets and the hollow core.

graand gaxare radial and axial gap in Fig.4,τwthe width of stator winding,btthe width of stator teeth,g the length of air-gap,τralength of radially magnetized PMs, τaxlength of axially magnetized PMs,τpthe pole pitch,Rsthe outer radius of mover,Rhthe outer radius of mover hollow core and Rshthe inner radius of mover hollow core.The magnetic field properties in the three regions are characterized by an expression between the magnetic field intensity HI，HII,HIII(in A/m)and the magnetic flux density BI,BII,BIII(in Tesla)as

Fig.4 Magnets topology and leading parameters of TMMLOM.

The key step to model a TMMLOM is to describe the distribution of magnetic flux,resulting from which that analytical methods can be grouped into three categories.One is lumped magnetic circuit models,7,24,36–38which consider a motor made of several parts with individual magnetic parameters,and they can only model the motor roughly due to the lumped approximation method.Another choice is the equivalent surface current method,which introduces several layers of surface current to take the place of PMs.37Its approximate structure can only outline the regularities of PMs with a simple geometric form.The last one is the analytical solution from Maxwell’s equations,16,18,23,25,28that shares a better approximation of the real magnetic flux.It is restricted that indispensable approximate conditions need be set to ensure the model solvable.

The magnetization vector M of the PMs array can be expressed in the cylindrical coordinates as

where Mrais the radial components of the magnetization vector,Maxthe axial components of the magnetization vector.Their spatial distributions are shown in Fig.5 where M0is the absolute value of magnetization vector which is given by

DICSSAC系统水和空气通过间壁式换热器强制，利用水热容大的特点及其蓄冷作用，使得冷却系统稳定，可以抵抗大风引起的背压波动。

where Breis the magnet remanence of permanent magnets.Equations of Mraand Maxderived from the con figuration of PMs as shown in Fig.5,where Rseis the effective stator inner radius,Rmis the effective PMs inner radius,and Rpis the effective hollow core outer radius.

They can be resolved into Fourier series,which are given by

The magnetic vector potential in cylindrical coordinates can be expressed as

“亦如汉睢阳王收拷谈生也”，出自魏曹丕志怪小说《列异传》中的《谈生》篇：谈生年四十仍无妇，夜半读《诗经》，有十五六岁女子来，与生结为夫妇。自言与人不同，三年之后始可以火相照。生一子，已二岁。谈生不能忍，偷以火照之，女腰以下枯，遂不能复生肉，赠一珠袍而别。此袍后为睢阳王家买去，发现是已故王女之袍，于是收拷谈生。谈生以实相告，王乃认谈生为婿，并表其子为侍中。

3.2.Distribution of flux density

According to the fundamental theory of an electromagnetic field,the magnetic field in the TMMLOM is a solenoid field or source-free field,namely,

原文：And since everyone is here, this is where we meet each other.

where i=I，II，III that represent the three regions.It can be derived that for any vector,the divergence of its curl is always equal to zero.Thus,a magnetic vector potential Aiin region i can be introduced,so that

总之，恩格斯为唯物矛盾观的深化发展做了大量研究工作，为乐曲添上了精彩的一笔，这表现在：其一，用大量事实说明了矛盾的客观性和普遍性；其二，从侧面补充和发展了唯物矛盾观的观点，如：“辩证的否定”“运动本身就是矛盾”等；其三，论述自然、历史和思维中的辩证法，并深入探讨了辩证法的基本规律、对象、性质及作用。但是，恩格斯并没有着手构建一个系统的矛盾体系，即使他提出了很多深刻的矛盾观点，遗憾的是其本身还是不成系统和体系的。

Furthermore,deducing from Appendix A,the governing equations in the three regions are obtained as

wherezthe axial position,and τlpthe cycle length utilized in the Fourier series.De fine thatand，and the coefficients are determined as

where Ai，rais the radial component of Ai,Ai，cithe circumferential component of Aiand Ai，axthe axial component of Ai.For the magnetic field distribution of a tubular linear oscillating motor is axially symmetric,Aihas only one component Ai，ci.Therefore,the governing equations can be transformed into

Fig.5 Distribution of magnetization vector.

whereare constants and r is the radial position.In(x)and Kn(x)are the modified Bessel functions of the first and second kinds of order n,respectively,which are defined as

再生细骨料的研究与应用，对于建筑垃圾的充分资源化和环境保护意义重大，符合我国国情和可持续发展的理念，同时也是建筑材料技术的一种拓展和革新。但目前我国在再生砂浆的大规模应用上还有许多问题亟待解决。

《卜算子慢》最早见《全唐诗·附词》中所录五代钟辐词。词调中有“慢”字的作品，唐五代仅仅有此一首。慢词创作，宋人可能直接取名《卜算子》，例如柳永的《卜算子·江枫渐老》，“调下原无‘慢’字，据毛本、《词律》、《钦定词谱》补”。柳词和张先的《卜算子慢·溪山别意》比较，张词较柳词，在字数上多添四字，句式上也异于柳词的各片均八句，上阕四仄韵，下阕五仄韵。句数上，张词双阕各添一句，用韵上，上阕五仄韵，下阕六仄韵。

Bi，ra,Bi，axare the radial and axial component of Biin region i,Hi，ra,Hi，axare the radial and axial component of Hiin region i.In regions I and III,

where i=I，III.In region II,

The compensation of the stator slot in the TMMLOM is accounted for by introducing a Carter coefficient Kc,40and the equations is given as follows:

Therefore,the effective air-gap ge,the effective stator inner radius Rse,the effective PMs inner radius Rmand the effective hollow core outer radius Rpare deduced as follows:

Utilizing the expressions of aI，k，bI，k，aII，k，bII，k，aIII，k，bIII，kin Appendix C,the distribution of the magnetic flux can be determined uniquely.

In the TMMLOM,the magnetic flux density is continuous in the normal direction of the interface between two mediums,while the magnetic field intensity continuous in the tangential direction.Accordingly,boundary conditions are given as

黄淮海地区为我国优势小麦产区，本文结合多年实地深入的调研和一线的工作经验，遵循小麦的生长发育规律，并在关键物候期给予科学化的管理，以期实现小麦的增产提质的目标，帮助广大种植户节本增效。

3.3.Models of electromotive force and thrust

Fig.6 Topology of an armature.

As shown in Fig.6,zdis the travel distance of the mover,and the con figuration of coils is composed of 4 pairs of windings.In terms of the positions of windings,the coils region can be divided into two categories.One is the normal coils winding whose arrangement is next to two other windings.The other one is the end coils winding which is located at the end of the TMMLOM.Moreover,the flux-linkage of a single winding Ψccan be deduced by integration as

where Ncis the number of coils turns in a single winding region.

Applying the flux distribution of a normal coils region calculated before to Eq.(17)yields

where

KWkreflects the influence of the windings con figuration,and the more turns that coils are arranged in the width,the better output the TMMLOM performs,while KSkis a parameter that contains the information of the structure of an armature and the topology of PMs.When neglecting the magnetic saturation effects,KSkkeeps constant.Hence,the electromotive force of a single winding can be given by

学生社团等第二课堂活动为培养学科竞赛人才提供了重要保障。外国语学院建立的学生国学社致力于在全校范围内广招对中华传统文化和人文知识有浓厚兴趣的青年人才，以书法、国画、民乐、诗词、相声等为活动载体，培育了一批对我国文化精神高度认同、具备一定国学素养、掌握一定传统艺术基础的学生骨干，而这个群体中的佼佼者大多成为了学科竞赛中的骨干力量。

where v is the instantaneous speed of the mover.Thus,the thrust of a normal coil Fcis obtained as

in which i is the instantaneous current in a single winding.In terms of the end coils,the thrust Fecis expressed by

As the TMMLOM has 6 normal coils windings and 2 end coils windings,the total thrust F is derived finally as follows:

Fig.12 are the thrust outputs affected by gaps in an end coils winding.The thrust curves of end coils are no more asymmetrical like those in normal coils.The reason is the unsymmetrical distribution of the flux density generated by end PMs.Both of the two end coils have an opposite force trend,which results in the approximate symmetry in dual end coils.

4.Simulation and validation

In this section,analytical calculations and FEM experiments are implemented based on the major parameters of the TMMLOM given in Table 1,and the model is validated by measured results.

4.1.Validation of the analytical model

Finite element model(FEM)analysis achieves a good accuracy in the magnetic flux distribution of a tubular linear motor.41To verify the analytical model proposed in the last section,an FEM is established in the ANSYS Maxwell software.The open-circuit magnetic field whose coils are not energized is simulated,and the flux density distributed in the air-gap is generated by a quasi-Halbach PM arrays form mover.The radial flux density at the radius of r=18.3 mm is illustrated in Fig.7(a)and the axial flux density in Fig.7(b).The calculated results from the analytical model are also depicted in both figures as well.As shown in Fig.7(a),the analytical results matches well with the red curves,which is precise enough to describe the flux regularity and calculate the output thrust.

Besides,the axial flux density obtained from the analytical results also shares a similar trend to the FEM data,and the practicability of the proposed model has been demonstrated.

The output force of the TMMLOM is related to the distribution of the radial density.According to the axial information of flux calculated above,the thrust can be derived analytically.The analytical relations and finite element results are exhibited in Fig.8.

Table 1 Leading parameters of the TMMLOM.

Symbol bt(mm) τw(mm) g(mm) τra(mm) τax(mm) τp(mm) Rsh(mm) Rh(mm) Rs(mm) Bre(T)Value 12 3 1 16 8 24 10 14 18 1.24

Fig.7 Validation of flux density.

Fig.8 Validation of thrust.

It can be validated that the normal coils and end coils have different regulations of output force.The normal coils,which produce the thrust primarily, fit very well with the FEM result.It is demonstrated that the thrust of the normal coils is a linear function to the current,and the small error in condition of high current is caused by the magnetic saturation.The end coils model,however,presents a larger error.Because the end effect of the liner motor,the magnetic flux is unsymmetrical with respect to the end coils.Due to the fact that the TMMLOM has 6 normal coils windings and only 2 end coils windings,what’s more,the force generated by the end coils is only half of that by the normal coils.Hence,the error is acceptable,and the thrust model is suitable in gap analysis.

4.2.Effect of the radial and axial gaps

4.2.1.Effect of the radial and axial gaps on the flux distribution

我国在电子商务相关立法方面还存在立法空白，相对于西方国家，我国对客户的信息安全、交易安全以及知识产权保护还是不到位。宿迁市司法机关和相关政府部门应当严格落实既有的法律和行政法规，维护从事电子商务双方主体的合法权益。通过调研，对于实践中立法的缺失，宿迁市人大常委会可以报请省级人大及其常委会在不抵触上位法的前提下，及时制定地方性法规，宿迁市政府也可以制定不与法律、法规相抵触的行政规章，从法律层面规范电子商务产业的发展。同时，在农村地区，宿迁市政府应当加强对相关法律的宣传工作，使得村民在从事电子商务过程中能够知法守法，并通过法律解决纠纷，维护交易的安全性和稳定性[7-8]。

After validation of the proposed model considering the mover gaps,the effect of gaps is discussed.For the radial magnetic flux is a key factor that affects the output,the distribution of the magnetic flux in the radial direction is analyzed considering radial gaps in Fig.9(a)and axial gaps in Fig.9(b),respectively.

Fig.9(a)shows the distribution of the radial flux in the airgap with different sizes of radial gaps.It demonstrates that radial gaps would decrease the amplitude comprehensively.Even though they do not change the waveform of the flux distribution,the value of the whole magnetic field declines heavily with increasing radial gaps.When it comes to axial gaps,the effect is different.As depicted in Fig.9(b),only the peak value of the flux field falls with these gaps.In the contrary,the central area of the flux distribution is not affected by gax.Eventually,the average value of the magnetic flux under the in fluence of axial gaps is much higher than that of radial gaps with the same gap thickness.Similar simulation results could also be drawn by the FEM results as shown in Fig.10,which have verified the conclusion of the proposed model.

Fig.9 Effect of gaps on radial flux density analytically.

Fig.10 Effect of gaps on radial flux density by FEM.

Comparing radial gaps with axial gaps,the former have a more significant impact on the radial flux density,and the negative effect could decrease the energy density of the magnet flied in the air-gap,as well as the power density in the TMMLOM.

4.2.2.Effect of the radial and axial gaps on the thrust

Using the validated model,the relationship between the thrust and mover gaps is further investigated.Fig.11 reflect the influence of gaps in a normal coils winding.

Assign the current as a constant of 2 A,and the result shows that the thrust varies with the mover position.When the stoke distance is less than 2 mm,the thrust is approximate to a stable value.Nevertheless,if the mover moves further,the thrust would decrease gradually.Because the width of radial magnetized PMs is 16 mm,the width of coils in a winding is 12 mm.The mover moves only in a range of plus and minus 2 mm,so the thrust would not decline,and the stable stroke could be defined as

有限温存，无限辛酸——道出了诸多人情的真相。犹如月盈则亏，犹如盛极而衰，犹如一切生命的终点无非死亡，万事万物阴暗或者说令人悲哀的一面，亦是人生谜底中令人无法回避、逃避的答案之一。

The radial and axial gaps are another aspect of causes,which lead to a thrust loss.As illustrated in Fig.11,both types of gaps lead to an overall drop of the thrust curve,and it declines faster in the effect of radial gaps compared with that of axial gaps.This result agrees well with the distribution of the flux density,for the output force is directly proportional to magnet remanence of permanent magnets Bre.

Fig.11 Effect of gaps on thrust in normal coils winding.

After the design finalization of the TMMLOM,Eq.(24)indicates that the thrust is a function of the current and the mover position.

The effect of mover gaps in end coils is unsymmetrical,either.The attenuation is obvious in a positive stoke and unobvious in a negative stoke.As the end coils are not in a strong magnetic area,the thrust as well as the decrement of the thrust drops heavily in the minus mover position.The other effect is the same as that in normal coils.

The total thrust of the TMMLOM is shown in Fig.13.For the reason that the total force depends mostly on the normal coils,the thrust property is similar to that with the normal coils.Hence,the total thrust of the motor is much more sensitive to the radial gaps than the axial gaps.

5.Comparison with experimental results

Fig.12 Effect of gaps on thrust in end coils winding.

He et al.measured the radial flux density of a TMMLOM prototype.37She also built an FEM and lumped magnetic circuit model to describe the radial flux density in the air-gap.The error of her analytical model is much greater than that of the FEM results.Her experiment results and FEM calculations are presented in Fig.14.The radial flux density marked black was measured by a Gauss meter with the mover assembled inside the stator.From the figure,the magnetic field intensity of the a mar an thine curve is obviously higher than the experimental points.Using the model proposed in this paper to express the magnet field,the analytical results without gaps are close to the FEM curve,which didn’t take the mover gaps into account.When taking the gaps into consideration,the mean value is close to the measured data,as shown in the red curve whose axial gap is 0.2 mm and radial gap is 0.5 mm.There still exist some errors because of the stator slot.The comparison con firms the rationality and accuracy of the proposed model.

6.Conclusions

This paper studies the effects of mover gaps in a TMMLOM on the distribution of the magnetic flux density and the output thrust.An analytical model is proposed considering radial and axial gaps,and its performance is validated through FEM experiments.

Based on this approach,both radial and axial gaps would weaken the performance of the TMMLOM.Radial gaps would lead to a comprehensive damping in the radial flux density in the air-gap,while the axial mover only decreases that at the end of the radial polarized magnet.Besides,85%of the thrust in the motor is produced by normal coils,so the effect of gaps on the total force is similar to that of normal coils.In addition,the existence of radial gaps may cause a thrust decay 5 times as much as that in fluenced by axial gaps.Thus,in the process of designing and manufacturing TMMLOMs,it is noticeable to pay more attention to limit the thickness of radial gaps than to axial ones.

Fig.13 Effect of gaps on total thrust.

Fig.14 Comparison with experimental results.

Acknowledgements

The authors would like to appreciate the supports of National Basic Research Program of China(973 Program)(No.2014CB046402),National Natural Science Foundation of China(Nos.51620105010,51575019,51675019,51505015),and 111 Program of China.

Appendix A

Because the curl of any scalar function’s(f)gradient is always equal to zero,we could have

which indicates that Aimay have multiple solutions.To uniquely determine its value,Coulomb gauge,∇·Ai=0,is applied as a constraint.Hence,we could have

The combination of Maxwell’s equations and Eq.(1)gives

Substituting Eq.(A2)into Eq.(A3)yields

where J(A/m2)is the current density in the field.In the air–gap of a linear oscillating motor,J=0.Therefore,the Laplace’s equation for RegionⅠis obtained as

Similarly,the Laplace equation for Region III is derived as

In Region II,the combination of Maxwell’s equations and Eq.(1)gives

For J=0,Eqs.(A2)and(A7)obtain the Poisson equation analogously for Region II as

Appendix B

Considering that expressionis any order of Fourier series of Ai，ci,since Ai，ciis only a function of r and θ,we could assume that

where r is the unique variable of Ri(r)and z the unique variable of Zi(z)according to the method of variable separation.

General solution to Laplace equation

The Laplace Eqs.(A5)and(A6)in cylindrical coordinates can be simplified as

For the magnetic field distribution of a tubular linear oscillating motor is axially symmetric,the expression becomes

namely,

where i=I，III andis any order of Fourier series of Aci.Substituting Eq.(B1)into Eq.(B4)yields

where r and z are independent variables,and the last term as a function of z must be a constant.Therefore,the following formula is established:

Then Eq.(B5)can be written as

Eq.(B6)becomes

So Eq.(B5)equals to

and there are three possible solutions according to the variation of λ2.

When λ2=0,the following equations are obtained:

The general solution to Eq.(B11)is

where C0，C1，C2，C3are constants.Thus

Forshould be a periodic function of z,Eq.(B13)is not a valid solution of

When λ2＞ 0,the following equations are obtained:

The general solution to Eq.(B11)is

where C0，C1，C2，C3are constants,Jn(x)is the Bessel function of the first kind of order n,Yn(x)the Bessel function of the second kind of order n,and Γ(x)the Gamma function,which are defined as

Thus

Forishould be a periodic function of z,Eq.(B17)is not a valid solution of

When λ2＜ 0,let λ =jmk,where j is the imaginary unit.The following equations are obtained:

The general solution to Eq.(B18)is

Thus

where C4，C5，C6，C7are constants.In this case,is a periodic function of z.Hence,Eq.(B20)is probably a solution of Laplace’s equation.

Eq.(7)in cylindrical coordinates is expressed as

Considering the constraint of symmetry at z=0,we have the radial component of the flux density Bra|z=0=0.Moreover,Substituting Eq.(B20)into Eq.(B21)gives

Thus we draw the conclusion that

Therefore,the general solution to the Laplace’s Eq.(B6)is

where i=I，III and ai，k，bi，kare constants.

General solution to Poisson equation

The Poisson Eq.(A8)in cylindrical coordinates can be simplified as

Substituting Eq.(4)into Eq.(B25)yields

Similarly,the general solution to the corresponding homogeneous equation of the Poisson Eq.(B26)is

Assuming thatis a particular solution of the Poisson equation,the general solution to the Poisson equation is obtained as

Substitutinginto Eq.(B26)gives

Substituting Eq.(B30)into Eq.(B29)yields

whereThe general solution to the modi fied Struve differential Eq.(B31)is

where Hn(x)are the Struve functions of order n and Ln(x)the modi fied Struve functions of order n.

Hence,the particular solution of the Poisson equation is

Substituting Eq.(B34)into Eq.(B28)derives the general solution to the Poisson equation as

Appendix C

Submitting Eqs.(12)and(13)into Eq.(14),we can yield the following equation:

11.Tutelea LN,Kim MC,Topor M,Ju L,Boldea I.Linear permanent magnet oscillatory machine:Comprehensive modeling for transients with validation by experiments.IEEE Trans Industr Electron 2008;55(2):492–500.

12.Kim YK,Gu BG,Rhyu SH,Jung IS.Analysis and design of slotless tubular linear actuator for high performance on the ecopedal system of vehicles.Int J Appl Electromagnet Mech 2012;39(1):817–23.

13.Kim YK,Gu BG,Jung IS,Won SH,Hur J.Analysis and design of slotted tubular linear actuator for the eco-pedal system of a vehicle.IEEE Trans Magn 2012;48(2):939–42.

14.Wang JB,Jewell GW,Howe D.General framework for the analysis and design of tubular linear permanent magnet machines.IEEE Trans Magn 1999;35(3):1986–2000.

The solutions of Eq.(C1) fix the undetermined coefficients in Eqs.(12)and(13).

References

1.Ebihara D,Watada M.Development of a single-winding linear oscillatory actuator.IEEE Trans Magn 1992;28(5):3030–2.

2.Watada M,Yanashima K,Oishi Y,Ebihara D,Dohmeki H.Improvement on characteristics of linear oscillatory actuator for artificial hearts.IEEE Trans Magn 1993;29(6):3361–3.

3.Redlich R.A summary of twenty years experience with linear motors and alternators.Proceedings of linear divers for industry applications;1995 May 31–June 2 Nagasaki,Japan.1995.p.1–7.

4.Boldea I,Nasar SA.Linear electric actuators and generatorsElectric machines and drives conference.Cambridge,UK.Piscataway,NJ:IEEE Press;1997.p.5.

5.Lee HK,Song GY,Park JS,Hong EP,Jung WH.Development of the linear compressor for a household refrigerator.Proceedings of fifteenth international compressor Engineering conference;2000 July 25–28;West Lafayette,USA;2000.p.31–38.

6.Lee H,Jeong SS,Lee CW,Lee HK.Linear compressor for airconditioner.Proceedings of the international compressor engineering conference;2004 July 12–15;West Lafayette,USA.2004.p.C047.

7.Liang HS,Jiao ZX,Yan L,Zhao LF,Wu S,Li Y.Design and analysis of a tubular linear oscillating motor for directly-driven eha pump.Sens Actuators A –Phys 2014;210:107–18.

8.Wang TY,Liang HS,Jiao ZX,He P,Yan L.Dynamics modeling and load analysis of linear motor for leha system.Proceedings of lf uid power and mechatronics(FPM);2015 Aug 5–7;Harbin,China.2015.p.1128–33.

9.Wang ZM,Jiao ZX,Li Y.Analysis,design and simulation of a collaborative rectification hydraulic pump.Proceedings of guidance,navigation and control conference(CGNCC),2014 Aug 8–10;Beijing,China.2014.p.2065–70.

10.Kim TH,Lee HW,Kim YH,Lee J.Development of a flux concentration-type linear oscillatory actuator.IEEE Trans Magn 2004;40(4):2092–4.

15.Bianchi N,Bolognani S,Corte DD,Tonel F.Tubular linear permanent magnet motors:An overall comparison.IEEE Trans Ind Appl 2003;39(2):466–75.

16.Wang JB,Howe D,Lin ZY.Comparative study of winding con figurations of short-stroke,single phase tubular permanent magnet motor for refrigeration applications.Proceedings of the 42nd IEEE IAS annual meeting;2007 Sept 23–27;New Orleans,USA.2007.p.31123–8.

17.Wang JB,Howe D,Lin ZY.Comparative studies on linear motor topologies for reciprocating vapor compressors.2007 IEEE international electric machines&drives conference;2007 May 3–5;Antalya,Turkey.2007.p.364–9.

18.Wang JB,Lin ZY,Howe D.Analysis of a short-stroke,singlephase,quasi-halbach magnetised tubular permanent magnet motor for linear compressor applications.IET Electr Power Appl 2008;2(3):193–200.

19.Wang JB,Howe D,Lin ZY.Design optimization of short-stroke single-phase tubular permanent-magnet motor for refrigeration applications.IEEE Trans Industr Electron 2010;57(1):327–34.

20.Zhu ZQ,Chen X,Howe D,Iwasaki S.Electromagnetic modeling of a novel linear oscillating actuator.IEEE Trans Magn 2008;44(11):3855–8.

21.Chen X,Zhu ZQ,Howe D,Dai JS.Comparative study of alternative permanent magnet linear oscillating actuators.Proceedings of 2008 international conference on electrical machines and systems.2008 Oct 17–20;Sapporo,Japan.2008.p.2826–31.

22.Zhu ZQ,Chen X.Analysis of an E-core interior permanent magnet linear oscillating actuator.IEEE Trans Magn 2009;45(10):4384–7.

23.Chen X,Zhu ZQ,Howe D.Modeling and analysis of a tubular oscillating permanent-magnet actuator.IEEE Trans Ind Appl 2009;45(6):1961–70.

24.Chen X,Zhu ZQ.Analytical determination of optimal split ratio of e-core permanent magnet linear oscillating actuators.IEEE Trans Ind Appl 2011;47(1):25–33.

25.Wang TY,He P,Yan L,Jiao ZX.Modeling,simulation and experiment study of electromagnetic performance for E-type series linear oscillating motor.Proceedings of IEEE international confer-ence on advanced intelligent mechatronics;2015 Jul 7–11;Busan,Korea.2015.p.607–12.

26.Jiao ZX,Wang TY,Yan L.Design of a tubular linear oscillating motor with novel compound halbach magnet array.IEEE/ASME Trans Mechatron 2016;22(1):498–508.

27.Wang JB,Howe D.Tubular modular permanent-magnet machines equipped with quasi-halbach magnetized magnets-part I:Magnetic field distribution,emf,and thrust force.IEEE Trans Magn 2005;41(9):2470–8.

28.Wang TY,Yan L,Jiao ZX,He P.Analytical modeling of linear oscillating motor with a mixed method considering saturation effect.Sens Actuators A – Phys 2015;234:375–83.

29.Jang S,Choi J,You D.Dynamic performance of tubular linear actuator with halbach array and mechanical spring driven by pwm inverter.IEEE Trans Magn 2006;42(10):3518–20.

30.Karunanithi R,Jacob S,Singh GAS,Damu C,Das M.Development of moving magnet type linear motor for dual piston compressor for pulse tube cryocooler.Adv Cryogenic Eng:Trans 2012;1434(1):525–31.

31.Nowak L.Dynamic FE analysis of quasi-axisymmetrical electromechanical converters.IEEE Trans Magn 1994;30(5):3268–71.

32.Sadowski N,Carlson R,Beckert AM,Bastos JPA.Dynamic modeling of a newly designed linear actuator using 3d edge elements analysis.IEEE Trans Magn 1996;32(3):1633–6.

33.Tomczuk B,Sobol M.A field-network model of a linear oscillating motor and its dynamics characteristics.IEEE Trans Magn 2005;41(8):2362–7.

34.Ji JH,Yan SJ,Zhao WX,Liu GH,Zhu XY.Minimization of cogging force in a novel linear permanent-magnet motor for artificial hearts.IEEE Trans Magn 2013;49(7):3901–4.

35.Kim KH,Park HI,Jeong SS,Jang SM,Choi JY.Comparison of characteristics of permanent-magnet linear oscillating actuator according to laminated method of stator core.IEEE Trans Appl Supercond 2016;26(4):1–4.

36.Wang TY,Jiao ZX,Yan L,Chen CY,Chen IM.Design and analysis of an improved halbach tubular linear motor with nonferromagnetic mover tube for direct-driven EHA.Proceedings of IEEE chinese guidance,navigation and control conference;2014 Oct 8–10;Yantai,China.2014.p.797–802.

37.He P,Jiao ZX,Yan L,Liang HS.Analysis of magnet layout in circumferential and axial direction for halbach PM arrays.Proceedings of guidance,navigation and control conference;2014 Oct 8–10;Yantai,China.2014.p.2013–8.

38.Wang TY,Jiao ZX,Yan L,He P.Design of novel double-layer compound stator for tubular linear oscillating motor.Proceedings of 2015 IEEE 10th conference on industrial electronics and applications(ICIEA);2015 Jun 15–17;Auckland,New Zealand.2015.p.1910–5.

39.Xia CL,Guo L,Wang H.Modeling and analyzing of magnetic field of segmented halbach array permanent magnet machine considering gap between segments.IEEE Trans Magn 2014;50(12):1–9.

40.Gu QS,Gao HZ.Effect of slotting in PM electric machines.Electr Mach Power Syst 1985;10(4):273–84.

41.Tomczuk B,Schroder G,Waindok A.Finite-element analysis of the magnetic field and electromechanical parameters calculation for a slotted permanent-magnet tubular linear motor.IEEE Trans Magn 2007;43(7):3229–36.