时谐磁场下频率对螺线管线圈电磁力场的影响研究毕业论文.pdf

时谐磁场下频率对螺线管线圈电磁力场的影响研究时谐磁场下频率对螺线管线圈电磁力场的影响研究时谐磁场下频率对螺线管线圈电磁力场的影响研究摘要电磁铸造技术是一种新兴的连铸技术，由于它是一种无模型连铸技术，铸件表面十分光滑平整，同时，其内部组织均匀细小，无偏析现象。因此电磁铸造技术可以说是生产铸锭的较好的方法之一。而对于电磁铸造技术，其工艺核心就是励磁线圈在金属中产生的时变电磁力场。针对这一时变电磁力场的研究，迄今为止国内外尚无十分明确的数值分析研究以探求其变化规律。本文将通过对谐波磁场下交变电磁力场中洛仑兹力的大小、方向的讨论来研究其变化规律，并通过一些数据及图表来进行详细分析。在电磁铸造工艺参数中，许多因素都会影响时变电磁力场，而本文特别针对励磁电流的频率这一特定因素，借助ANSYS软件进行二维模型感应电磁场的初步计算并输出结果，而后使用指导教师开发的FORTRAN数值研究程序对ANSYS导出的数据文件进行处理，得出感应电磁场的幅值和初相位数据并据此计算得出时变电磁力场的特征数据，而后利用ORIGIN软件对这些力场特征数据进行绘图并分析，ANSYS计算分析中涉及的其他因素如金属及线圈的电阻率和相对磁导率，以及环境温度等均设为常量。对190Hz频率下的力场分析研究结果表明：时变电磁力场中金属中的每一个点上的时变力的端部都会以椭圆的形式进行大小及方向的变化，而励磁电流的频率会影响力椭圆的大小、旋向及其分布区域。随着励磁电流频率的增大，力椭圆的长轴以及短长轴比都会随之增大，而其原本有旋向特征的区域，其旋向保持不变，没有旋向的区域则不断减小，转化为有旋向的区域。关键词：电磁铸造，时变电磁力场，数值分析，励磁电流的频率时谐磁场下频率对螺线管线圈电磁力场的影响研究ASTUDYOFTHEINFLUENCEOFFREQUENCYONTHETME-VARYINGLORENTZFORCEFIELDUNDERHARMONICMAGNETICFIELDABSTRACTElectromagneticcastingtechnologyisanewcontinuouscastingtechnologybecauseitisanon-modelcastingtechnologycastingsurfaceisverysmoothwhileitsinternalorganizationunilyfinenosegregation.Thereforeelectromagneticcastingtechnologycanbesaidtobeoneofthebestwaytoteachtheproductionofingots.Forelectromagneticcastingtechnologyitscoretechnologyistheexcitationcoilsinthemetalalternatingelectromagneticforcefield.NumericalAnalysisofElectromagneticFieldvaryingstudyonthiswhenthereisnoveryclearsofarathomeandabroadtoseektheirvariation.ThisarticlewilldiscussthealternatingmagneticfieldharmonicelectromagneticforcefieldLorentzforcesizedirectiontostudythevariationandtocarryoutadetailedanalysisofdataandthroughsomecharts.Inelectromagneticcastingprocessparametersmanyfactorsaffectthetime-varyingelectromagneticforcefieldandthefrequencyoftheexcitationcurrentarticleespeciallyforthisparticularfactorusingANSYSsoftwaretocalculatetheinitialtwo-dimensionalmodelofelectromagneticfieldandtheoutputandthenusetheinstructordevelopmentFORTRANNumericalStudyprogramANSYSexporteddatafileprocessingdataobtainedamplitudeandinitialphaseinductionelectromagneticfieldandcalculatesthealternatingelectromagneticforcefieldwhencharacteristicdataobtainedthentheuseoftheseforcefieldsORIGINsoftwarefeaturesdataplottingandanalysisANSYStocalculateotherfactorsinvolvedsuchastheresistivityofmetalandthecoilandrelativepermeabilityandissettoaconstantambienttemperatureetc.OnField190Hzfrequencyanalysisresultsshow:时谐磁场下频率对螺线管线圈电磁力场的影响研究time-varyingelectromagneticforcefieldinthemetalofeachpointonthevariableforceendsareintheofanellipsevarythesizeanddirectionandtheexcitationcurrentfrequencywillaffectthesizeoftheforceoftheellipserotationanddistributionarea.Withtheincreasingfrequencyoftheexcitationcurrentpowerandshort-majoraxisoftheellipsemajoraxisratioincreaseswillbeanditsoriginalfeatureshaveregionalrotationwhichrotationremainsunchangedthereisnorotationoftheregioncontinuetoSavesmallconvertedtohaveaspintothearea.Keywords:ElectromagneticCastingTime-varyingelectromagneticforcefieldNumericalAnalysisFrequencyofexcitationcurrent时谐磁场下频率对螺线管线圈电磁力场的影响研究目录1绪论.11.1引言.11.2研究现状.21.2.1技术研究现状.21.2.2理论研究现状.31.3研究内容及方法.41.3.1研究内容.41.3.2具体研究方法.52计算原理及模型条件.62.1时谐电磁场.62.1.1时谐电磁场.62.1.2感应电磁场的特点.62.2二维简化模型的建立.72.2.1三维模型.72.2.2二维模型简化.72.2.3模型属性的设置.82.3力场分析原理.93数值模拟分析.103.1ANSYS二维时谐磁场分析.103.1.1ANSYS介绍.103.1.2ANSYS分析过程.103.2FORTRAN数据处理.123.2.1FORTRAN介绍.123.2.2FORTRAN处理过程.123.3ORIGIN绘图.12时谐磁场下频率对螺线管线圈电磁力场的影响研究3.3.1ORIGIN介绍.123.3.2ORIGIN绘图过程.134力场分析.144.1f=60Hz时电磁力场的分析结果研究.144.1.1幅值与相位的变化规律.144.1.2电磁力.164.2幅值及相位随f的变化.214.2.1幅值随f的变化.214.2.2相位随f的变化.234.3电磁力随f的变化.254.3.1力轨迹随f的变化.254.3.2力轨迹旋向随f的变化.264.3.3焦耳热随f的变化.275结论.30参考文献.31致谢.32译文及原文.34时谐磁场下频率对螺线管线圈电磁力场的影响研究11绪论绪论1.1引言电磁铸造（EMC）技术1始于20世纪60年代，由前苏联的Getselev等人发明。这是一种利用电磁感应原理进行铸造的无模连铸技术，对于铸造领域而言，这是一种飞跃性的突破。长期以来困扰铸造领域的几大问题如材料的利用率，铸件内部组织结构问题，铸件表面处理问题等在电磁铸造技术出现以后都有望得到解决。与传统的直接冷却（DC）法相比，电磁铸造的铸模是一个电磁线圈。在感应线圈中通上中频电流后，液态金属内部就会产生感应电流，感应电流又既而产生感应磁场，在感应磁场与感应电场相互作用下，又会产生感应时变电磁力场。这种时变电磁力场会使得液态金属不流散而形成稳定的液柱。与此同时，铸锭下方的冷却水喷口朝上喷水，使液态金属逐渐凝固，接着铸机拖动底模向下，并供给一定量的液态金属，使其得到补充。这就是电磁铸造的大体过程。在电磁铸造过程中，液态金属没有铸模，因而是在自由状态下凝固成型，所以，铸造出来的金属坯件的表面十分光滑平整，完全不存在DC法中出现的金属表面粘附，偏析瘤等问题。作为坯件，它甚至可以不用精加工，这使得金属的利用率大大的提高。再者，由于液态金属始终处于电磁力场的搅拌状态，其内部组织均匀细小，因而极大地减少了偏析现象。所以说，电磁铸造技术为连铸的质量和经济性提供了更多的可能。因此自其问世以来，就一直被人们广泛关注。对于这样一种极具时代意义的铸造技术，人们越发迫切的想要去了解和掌握它的核心。在电磁铸造的整个过程中，时谐磁场2一直扮演着重要的角色。理论上来说，只要液态的铸造材料能导电，时谐磁场就能通过改变其感应电磁力场的变化，既而控制其温度场、流场、成分场以及其耦合过程，从而达到控制凝固成型的目的。因此可以说电磁铸造的核心就在于对时谐磁场所产生的电磁力场的研究。在时谐磁场的影响下，感应电磁力场中洛仑兹力的大小及方向既会随着时间时谐磁场下频率对螺线管线圈电磁力场的影响研究2的变化而变化，也会受到其空间位置的影响。其变化规律非常复杂。现阶段的数学分析表明，时变洛仑兹力的变化周期是有规律可循的，其周期为励磁电流周期的一半。因此时变电磁力场是一种有固定周期的矢量力场，这为我们的研究提供了思路和切入点。在数学分析的基础上，进而对时变洛仑兹力场的特征进行数学推导，便可以将其进行定量化的描述。本文的研究意义就在于借助软件分析和数据处理技术，来定量地研究时变电磁力场的变化规律，如变化轨迹、轨迹旋向、力场结构分布等。这将十分有助于对材料的电磁处理研究，从而为进一步的电磁铸造技术发展提供理论依据和数据基础。1.2研究现状1.2.1技术研究现状电磁铸造技术发展至今已有50余年，对于电磁铸造技术的应用也越来越广泛和深入。1960年，苏联的Getselev等人首先提出了电磁铸造法(EMC)。而在6年后人们首次通过EMC技术制得了第一个铸锭。1969年其开始应用于工业生产，并成功地铸造出了直径为200-500mm的铸锭。后来这种技术在捷克斯洛伐克、匈牙利、东德等东欧各国迅速发展。到了1973年美国Kaiser公司申请了专利技术，同年Alusuisse公司也开始研究改进这项技术。另外许多大型的制铝工厂也同时期引入了这项技术，如Alcoa、Reynolds、Pechiney、日本三菱化成等公司。从苏联引进的EMC技术只有最基本的单块铸锭生产技术这还远远满足不了西欧等国对大批量、高质量铸锭的工业和军工需求。所以在该项技术被引进之后每家公司都在不断地改进和发展这项技术。1975年，Alusuisse公司美国分公司Comalco开始生产和研发较为成熟的整套设备其生产的连铸设备能进行4块尺寸为500 x13002mm的铸锭的连续铸造。1981年，Kaiser公司则研发出了5块尺寸为500 x13002mm的铸锭的连续铸造设备。在EMC技术的发展过程中能实现工业生产的公司只有Kaiser及Alusuisse它们进行了数十年的研究才有了这样的进步。EMC技术可以是一种非常有技术含量的铸造方法。目前来看，其技术发展现状主要体现在不同形状铸坯的连铸、不同材料铸坯的连铸、不同规格铸坯的连铸以及不同铸造方式的连铸34。时谐磁场下频率对螺线管线圈电磁力场的影响研究31.2.2理论研究现状从1966年生产的第一块铸锭到今天的批量化生产铝合金，电磁铸造技术一直在不断地发展着。其应用领域从单一的铝合金的连铸延伸至铜、钢等常用金属的连铸。我国在这一方面也取得了一定的成绩，发展至今已铸造出多种优质铝合金铸锭并对EMC温度场、应力场进行了一定数量的实验与数值模拟的研究。然而由于电磁铸造过程中产生的时变电磁力场十分复杂，对于其定量特征的研究一直止步不前。目前来看，对于时变电磁力场的研究主要有以下几种处理方式：第一种方式是在简化模型的基础上，通过特征长度、电导率和空载磁感应强度来计算出洛仑兹力5，从而进行各种感应器性能的研究。然而这种方法既忽略了时变电磁力场的时变特征，亦忽略了其随空间位置变化的特性。所以不能很好的对时变电磁力场进行定量描述。第二种方式是在空载条件下，仅仅分析感应器产生的感应磁场，从而估算洛仑兹力的大小和方向；这样就完全忽略了感应器产生的感应电场和材料本身电阻率及相对磁导率的影响。所以这种方式也不能较好的定量描述时变电磁力场的特征。第三种方式是通过数学处理，将所研究的时变洛仑兹力拆分为两部分。其一被称为时均洛仑兹力，其变化是一个仅与力场位置相关的三角函数，与时间变化并无关系，这一部分的力较好研究，可以通过ANSYS等软件直接进行分析处理从而得出具体的数据。其二则称之为脉动洛仑兹力，其变化是一个与时间和力场位置都相关的三角函数，变化周期为所施加的励磁电流周期的一半。对于这一部分很难直接用商用软件进行分析处理，所以较难研究。在进行电磁铸造过程中流场、温度场、电磁场的耦合计算6时，许多学者会选择忽略脉动洛仑兹力而只研究和运用时均部分，并以此来代替整个时变洛仑兹力的变化特征。有研究表明在励磁电流频率较高时，脉动部分6的影响相对较小，忽略这一部分进行研究分析是合理的。但当励磁电流频率较低时，脉动洛仑兹力的影响是不能忽略的。所以这种方法并不是十分完整的，其在进行低频励磁电流产生的时变洛仑兹力的研究时捉襟见肘。因此，可以说虽然学者们对于电磁铸造过程中的物理场都有所深入研究，如时谐磁场下频率对螺线管线圈电磁力场的影响研究4流场（大多数以时均部分作为源项）、温度场、成分场等，但对于电磁力场中的脉动洛仑兹力的周期性特征，研究者们却一直缺乏详尽的了解，其周期性的空间时变规律一直难以描述，如随着时间的变化，时变洛仑兹力以何种形式发生改变，其大小、方向、变化轨迹等都没有明确而细致的描述。而这其中的难点就在于对脉动洛仑兹力的变化特点未能把握。因为时变洛仑兹力的脉动部分是时变的主体部分，其变化会随时间和空间的改变而改变。这是一直困扰着当今学者们的一大问题。1.3研究内容及方法1.3.1研究内容（1）了解电磁铸造的发展，从而明确课题的背景和意义电磁铸造技术发展至今已逾50年，从起初的铝锭连铸已发展至如今的铜、钢的铸造。其铸造方式和铸型都有所增加。随着电磁铸造被越来越多的关注和应用，其理论分析的支持理应跟上技术的应用需求。本次课题正是从当前较为难解决的含脉动洛仑兹力的电磁力场入手，着重分析脉动部分的时变规律，通过软件的仿真分析和辅助计算对其变化轨迹、轨迹旋向等特征给出详尽的定量描述，并绘制出直观的图表使其变化特征更清楚的表现出来。这将有利于我们精准的掌握感应电磁场的周期性变化特征，加深我们对各种因素对周期性时变电磁力场的影响的理解，为之后电磁铸造感应器等机器设备的设计和改进提供可靠的理论依据，从而推动电磁铸造技术更深入的发展。（2）了解时谐磁场的概念及其感应电磁力场的特点了解时谐磁场是一种以正弦规律变化的磁场，其感应电磁力场是具有时间和空间分布特征的场。（3）了解模型的计算原理和方法了解通过简化模型和数值仿真来进行理论分析的方法。（4）了解时变电磁力场的分析原理了解其分析原理是通过数学处理将力场中的力拆分成两部分后进行分析。并依据其幅值及相位等的基础数据得出具体的力场特征数据。（5）进行时变电磁力场的分析，研究励磁电流频率对其产生的影响归纳不同的特征进行电磁力场的定量描述，并着重研究随着励磁电流频率时谐磁场下频率对螺线管线圈电磁力场的影响研究5的改变，其特征有何种变化。1.3.2具体研究方法（1）将实际三维模型简化为研究所需的二维模型，使数据量减少，以便于实际研究操作和简化分析。再通过ANSYS仿真分析软件对简化模型进行二维时谐磁场分析，得出原始感应电磁场数据。（2）运用FORTRAN语言，通过数学模型和实际物理特征将原始感应电磁力场数据转化为时变洛仑兹力的特征数据。（3）通过ORIGIN绘图软件，将生成的时变洛仑兹力特征数据绘制成所需的二维矢量图、三维散点图，更直观的去观察感应电磁力场随励磁电流频率的变化规律。（4）从力的变化轨迹、轨迹旋向、力场结构分布等方面进行力场分析，给出时变电磁力场的定量描述，然后进行不同励磁电流频率下，其电磁力场特征的比较分析，研究其变化规律并得出结论。时谐磁场下频率对螺线管线圈电磁力场的影响研究62计算原理及模型条件计算原理及模型条件2.1时谐电磁场通过对螺线管线圈施加励磁电流，就会产生时谐磁场。时谐电磁场下频率的变化会影响到其产生的电磁力场。2.1.1时谐电磁场在时变电磁场中，如果其场源随着时间的变化，以一定的角频率做时谐（正弦或余弦）变化，则其所产生的电磁场也以相同的角频率呈时谐变化。这种以一定的角频率作时谐变化的电磁场，我们称之为时谐电磁场或正弦电磁场。所以时谐电磁场是一种感应电磁场，其角频率与所施加的励磁电流频率是一致的，并且其变化规律是正弦变化规律，其有固定的周期、幅值及相位，只要通过改变励磁电流的频率就能改变时谐电磁场的频率，这是本次研究的重要条件之一。2.1.2感应电磁场的特点本次研究中涉及的感应电磁场是一个时变电磁场，而任意的时变场都可以表示为不同幅值及相位的时变场的叠加。我们讨论的是二维时谐磁场问题，所以不妨让励磁电流zJ和感应涡流GzJ（xy）与Z轴平行，感应磁场xB（xy）和yB（xy）均在XOY平面上，它们的频率和角频率分别设为f和，则励磁电流可以写作：)2sin()sin(ftJtJJzmzmz（2-1）而涡流区中某一点的感应电磁场则可以表示为：)(sin)()()(sin)()()(sin)()(321yxtyxJyxJyxtyxByxByxtyxByxBGzmGzymyxmx（2-2）上式中，下标m表示幅值，t为时间，)321(ii表示各分量相对励磁电流的相位。其计算方法见参考文献7此处不做赘述。由此，我们就可以量化涡流区中的感应电磁场了。时谐磁场下频率对螺线管线圈电磁力场的影响研究72.2二维简化模型的建立通过对模型的简化能使得研究的计算量及难度降低，同时简化的模型能更直观的反映其变化特征，本节将介绍本次研究中对于三维模型向二维模型的简化过程。2.2.1三维模型本次研究的三维模型包括一个一定匝数的螺线管线圈和与其同轴的金属圆柱。螺线管由实心铜金属棒绕制而成。金属圆柱的直径为80mm，高度为100mm，线圈铜棒的直径为10mm。其三维图如下图所示：图2-1螺线管线圈的三维图图中，黄色部分为螺线管线圈，灰色部分为金属圆柱。线圈和金属圆柱之间存在间隙。2.2.2二维模型简化在二维模型中，将通过圆柱中心轴截取任意平面作为简化的二维平面。由于其为对称模型，故在二维模型的建立中，只需截取截面的一半进行分析即可。另外为了简化模型，将线圈数设置为3个，每个线圈的圆心距为15mm，圆心与界面边线的距离为10mm。简化后的二维模型如下图所示：时谐磁场下频率对螺线管线圈电磁力场的影响研究8图2-2简化二维模型图中矩形为过轴心的金属圆柱截面，尺寸为（40100）2mm，三个圆形为线圈截面。通过这种方式简化的二维模型更便于我们进行数值模拟分析。2.2.3模型属性的设置在本次分析案例中，金属圆柱所选取的材料为金属钛（Ti），线圈所选取的材料为金属铜(Cu)。另外在ANSYS谐波电磁场分析中，需确定近场区和远场区特性，此处将其均设置为空气，近场区直径为500mm远场区直径为1000mm。具体材料特性如下表所示：表2-1材料特性模型材料电阻率m相对磁导率金属圆柱线圈近场区及远场区钛铜空气61071.1810673.1111注：由于本次数值分析中所选取的频率较小，由其对材料产生的热效应等不做考虑，环境温度稳定20，且材料特性也不会随温度变化，为恒定值。因为励磁电流的频率较低时，脉动洛伦兹力的影响较为明显，故本次数值分析选取了8个较低的励磁电流频率，其具体数值为1Hz、10Hz、20Hz、30Hz、45Hz、60Hz、75Hz、90Hz。线圈上施加的励磁电流密度均为26101mAJs。时谐磁场下频率对螺线管线圈电磁力场的影响研究92.3力场分析原理由式（2-1），我们可以计算出电磁体积力F（3mN）。设i、j分别为X、Y轴的单位矢量，则可以得到以下公式：pulsatingmeanxmGzmxmGzmymGzmymGzmGFFtbtabajtBJjBJitBJiBJBJF)()()()2cos()cos()2cos()cos(210013132323（2-3）式中，电磁体积力由0a、0b、)(ta、)(tb四项组成，0a和0b分别表示i、j方向上与时间变量不相关的分量，)(ta和)(tb则表示i、j方向上与时间变量相关的分量。由此可将电磁体积力分为两部分8，其一为时均洛仑兹力meanF，其二为脉动洛仑兹力pulsatingF。jBiBJFxmymGzmmean)cos()cos(211323（2-4）jtBitBJFxmymGzmpulsating)2cos()2cos(211323（2-5）由式（2-4）、（2-5）不难看出，meanF由0a和0b组成，其为与时间变量不相关的函数，仅与该点在涡流区中的位置相关。而pulsatingF则由)(ta和)(tb组成，其为时间变量相关函数，不仅受到该点在涡流区中的位置影响，还会随时间变化而变化。通过上述数学处理，我们可以清晰的对感应电磁力场进行量化描述，将时变洛伦兹力拆分为了时均部分和脉动部分。通过两个部分的数据处理，我们就能计算出电磁体积力的大小和方向等特征信息。进而研究和分析出每一点上的时变洛仑兹力的变化规律。这就是本文的力场分析原理9。时谐磁场下频率对螺线管线圈电磁力场的影响研究103数值模拟分析数值模拟分析通过对三维模型的简化，我们可以借用软件对二维模型进行数值模拟分析，这里我们用到的软件主要有ANSYS、FORTRAN、ORIGIN。本章将简要介绍其进行数值分析及数据处理的过程。3.1ANSYS二维时谐磁场分析3.1.1ANSYS介绍自问世以来，数值模拟技术在工程中就有着十分广泛的应用。其通过计算机软件程序能实现非常强大的仿真分析功能。80年代初期，世界上应用于工程的通用有限元程序就已达到几百种，其中为人们所熟知的有：ANSYS，NASTRAN，ASKA，ADINA，SAP等。它们大多是用FORTRAN语言来编写的，编写量有几万条甚至几十万条源代码，所以其作用也越来越大。其不仅拥有能处理多种情况的有限元分析程序，还有功能强大的前处理和后处理程序。因为有限元通用程序使用方便且计算精度高，其数值仿真结果已成为了各种工业设计和生产过程的可靠的理论数据依据。而在众多的通用有限元软件中，ANSYS10是最为人们所广泛接受的一款软件。ANSYS软件是由美国宾夕法尼亚州匹兹堡的ANSYS公司研发的有限元分析程序。其可以通过有限元的划分，实现结构、力场、电场、磁场和流体的仿真分析，是一款功能强大的大型通用CAE软件。自其问世以来，在各种科学研究和工程设计领域做出了巨大的贡献，为诸多研究者提供了便利可靠的仿真分析功能。3.1.2ANSYS分析过程（1）定义首选项ANSYS中有诸多分析类型，选择正确的分析类型能简化分析过程中不必要的特征数据，使得分析结果更加准确。所以首先要确定数值分析类型，本例中要进行二维时谐磁场分析，所以选择Magnetic-nodal(电磁节点法)。（2）定义单元类型时谐磁场下频率对螺线管线圈电磁力场的影响研究11第二部是要选择单元类型，不同的单元类型代表不同的材料区和自由度。本例中模型的材料为钛、铜和空气，且为二维平面分析，所以选取PLANE53作为单元类型一，另外，在ANSYS分析中需要设置近场区和远场区，近场区表示理论上受到实验影响的区域，远场区则表示无限远处的区域（这是有限元的一种定义无限远处的方法）。因为近场区空气的实验特性，所以PLANE53同样适用于近场区。而远场单元选取INFIN110作为单元类型二。两种单元的设置中，自由度均设置为AZ方向，其他设置不变。（3）定义材料特性不同的实验材料具备不同的物理特性，在ANSYS中要进行其特性的区分。所以分别设置钛、铜和空气的电阻率及相对磁导率。其详细数值见表2-1。（4）创建平面二维模型材料的单元类型及材料特性设置完成后，就可以开始创建分析模型了。在前文中已提到，为方便仿真分析，对实际的三维模型进行了二维简化。如图3-1所示为ANSYS中绘制的简化二维模型，其中较小的半圆表示近场区域，直径为500mm，较大的半圆表示远场区域，直径为1000mm。其他尺寸值参见图2-2。另外，需要注意的是，为了正确区分不同的区域，需对所有区域进行布尔搭接运算，使重叠区域消除。图3-1ANSYS绘制的二维模型图（5）划分网格划分网格是有限元分析的重要步骤。也是成功进行仿真分析的关键。在划分网格过程中，首先要为每个区域赋予材料特性及单元类型，这在之前的步骤中已经设置完毕，只需一一对应即可。接着进行网格划分。需要注意的是靠近线圈的时谐磁场下频率对螺线管线圈电磁力场的影响研究12位置，磁感应线较为密集，划分网格时此处网格应适当密集一些。另外在能选择Mapped划分的情况下应选择Mapped划分，因为其划分的网格具备更高的规整度，可使最终的仿真结果更加精确。（6）求解进行最后的求解前，首先要定义分析类型为Harmonic（谐波分析），然后为线圈施加电流密度并输入电流频率。接着进行Solve（求解）操作，并选择从当前载荷步开始求解。需要注意的是，在求解之前要进行SelectEverything的操作，否则系统默认选取上一步选中的对象进行求解。这将使得正常数据结果无法显示。（7）输出数据在本例中，为了得到下一步数据，我们需要输出的数据有节点坐标（Nodes-Coordinatesonly）、单元节点特性（Element-Nodes+Attributes）、节点磁通密度、矢量磁势以及单元磁通密度、电流密度和焦耳热等。3.2FORTRAN数据处理3.2.1FORTRAN介绍FORTRAN作为计算机高级程序设计语言已有很长的历史，它最强大的功能就在于可以直接对矩阵和复数进行复杂的运算，其最大特性就是接近数学公式的自然描述，并且在计算机里具有很高的执行效率。深受理论研究人员的喜爱。3.2.2FORTRAN处理过程将ANSYS导出的数据文件直接导入已编好的FORTRAN程序中，就能直接生成所需要的时变电磁力场特征数据。这些数据包括感应电磁力场的磁通量密度数据、焦耳热数据，还有电磁力场特征数据，如时均力数据、电磁体积力数据以及每一节点洛仑兹力变化的最大值和变化轨迹等。这些数据将对我们的分析研究有着至关重要的作用。3.3ORIGIN绘图3.3.1ORIGIN介绍ORIGIN11是由美国OriginLab公司研发的专业绘图软件，其主要具有两大功能：绘图和数据分析。绘图功能包含几十种二维和三维绘图模板，而数据分析时谐磁场下频率对螺线管线圈电磁力场的影响研究13功能则可以进行包括统计、信号处理、图像处理、峰值分析和曲线拟合等各种完善的数学分析。3.3.2ORIGIN绘图过程（1）导入数据ORIGIN可以直接读取TXT格式的文件，所以我们只需将所需文件拖入ORIGIN的Worksheet中就能直接导入数据了。（2）设置坐标轴在ORIGIN绘图前，我们需要先设置坐标轴。根据需要我们可以将不同列的数据分别设置为X轴、Y轴和Z轴。（3）绘图ORIGIN绘图功能中一共有几十种二维和三维模板，我们需要绘制的是三维散点图（3DScatter）、三维曲面图（3DSurface）以及二维矢量图（VectorXYAM）。前两者需要在数据中设置X、Y、Z三个轴，而后者需要在数据中设置X、Y、Y、Y四个轴。需要注意的是ORIGIN是根据固定模板绘图的，所以我们需要对图片进行适当修改，以符合我们的实际要求。（4）导出图片在ORIGIN中有十几种图片格式可供选择，从图片处理方式及清晰度的角度考虑，我们选择导出PNG格式的图片。根据分析需要，我们最终需要导出的图片有各个频率下的磁通密度B的幅值变化三维散点图xmB、ymB以及相位变化三维散点图1、2，感应电流密度GJ的幅值变化三维散点图zmJ以及相位变化三维散点图3，焦耳热三维曲面图JouleHeat，涡流区节点力最大值图maxF，力变化轨迹主轴图major-axis，力变化轨迹旋向图sn（顺逆）以及力变化轨迹主轴方向图。至此，前期的模型仿真分析，以及数据处理已全部完成。通过ANSYS、FORTRAN以及ORIGIN的帮助，我们成功地将抽象的三维时谐磁场分析结果转化为直观的数值数据及图片，使我们能够更清晰的去观察力场的变化特征。另外还有一些辅助软件如Solidworks、PS等未及详尽介绍，通过这些软件我们能建立三维模型和二维动态图来辅助分析，此处不做赘述。时谐磁场下频率对螺线管线圈电磁力场的影响研究144力场分析力场分析4.1f=60Hz时电磁力场的分析结果研究感应电磁力场会随着励磁电流的频率变化而产生大小及方向的改变，为了对力场的描述有更清晰的认识，本节将先从某一特定频率下的感应电磁力场特征进行分析，这里选取频率的中间值60Hz作为分析对象。4.1.1幅值与相位的变化规律当f=60Hz时，磁通量密度B的幅值xmB、ymB，感应电流密度GJ的幅值zmJ以及各自相对应的相位1、2、3分布如图4-1所示。图中水平横坐标与水平纵坐标表示每个节点在感应电磁场中的位置，纵坐标表示该点对应的数值。需要注意的是由于数据处理上会存在误差，所以图中会存在一些异常点，不影响我们对整体情况的研究。a)b)时谐磁场下频率对螺线管线圈电磁力场的影响研究15c)d)e)f)图4-1感应电磁场的幅值及相位分布（f=60Hz）a)xmB，b)1，c)ymB，d)2，e)zmJ，f)3由图4-1a)可以看出，xmB在y=0处（三个线圈竖直方向的中点位置）取得最时谐磁场下频率对螺线管线圈电磁力场的影响研究16小值，而ymB、zmJ则在此处取得最大值，见图4-1c)、e)。再看图4-1b），1在y=0处发生跃变8，当y0时，其值约为1.5（将图中数据除以后得到此值），当y0（金属圆柱上半区）时，1与2相近，xmB与ymB几乎同相位。而当y01.51.52y0部分即可。当3t时靠近线圈位置的电磁力方向约与X轴呈-45，当t时靠近线圈位置的电磁力方向约与X轴呈75，当35t时靠近线圈位置的电磁力方向约与X轴呈160。而电磁力的大小也随之发生变化。其变化轨迹近似为椭圆，旋向为逆时针旋向（y0部分和y0部分，力椭圆的主轴为2、4象限方向，在y0部分，力椭圆的主轴为1、3象限方向。对比图4-7，其表现的变化规律也是一致的。如图4-7所示为力椭圆的短长轴比值图，从图中我们发现，短长轴比值的最大值约为0.07，出现在（4050）、（40，-50）处（金属圆柱表面外侧的两端）。这个值非常小，所以我们绘制出的图4-4所显示的椭圆轨迹为一个扁长的椭圆。对比图4-5，我们发现图4-7底部凸起部分的边缘值为0，这一边缘所表示的区域与图4-5的区域边界线相一致。另外，在y=0处，其值也为0。说明在y=0处电磁力的变化轨迹也有椭圆退化成为一条线段。这在图4-5中并未表现出来。如图4-8所示为力椭圆的主轴模长图，可以看出，在靠近（400）处的主轴较长，可知，该点的力椭圆较为扁长，而该点在一个变化周期内的电磁力最大值可能也较大。对比图4-9我们发现这个猜测是正确的，且最大电磁力的方向与主轴方向相近。由上述分析，我们发现，时变电磁力场中的力是规律性的变化的，其变化轨迹为一椭圆（或退化成一条线段）。而其轨迹的旋向及分布也有一定的规律。4.2幅值及相位随f的变化从前文分析中已知感应电磁力场中磁通量密度B的幅值及相位的变化，那么其随励磁电流频率f的变化规律是怎样的呢？4.2.1幅值随f的变化如图4-10所示为f=1Hz、f=45Hz、f=90Hz时的幅值图。a)b)时谐磁场下频率对螺线管线圈电磁力场的影响研究22c)d)e)f)图4-10感应电磁场幅值图：a)、b)f=1Hzc)、d)f=45Hze)、f)f=90Hz由图4-10我们发现随着频率的变化，感应电磁场的幅值并无明显变化，为了确定其具体变化情况，绘制了图4-11。图4-11为xB的幅值变化率图。其最大值约为-1.2%，其值非常小。所以随着f的改变，电磁场的幅值变化是非常小的。图4-11磁通量密度幅值变化率图时谐磁场下频率对螺线管线圈电磁力场的影响研究23如图4-12所示则为f=1Hz、f=45Hz、f=90Hz时的感应涡流幅值图，通过对比不同频率的感应涡流幅值图我们发现，zmJ的图像表面越来越陡峭，这表示随着f的增加，zmJ的变化率越来越大。a)b)c)图4-12感应涡流幅值图：a)f=1Hzb)f=45Hzc)f=90Hz4.2.2相位随f的变化如图4-13所示为f=1Hz、f=45Hz、f=90Hz时的相位图。从图4-13中可以看出1上的点做线性变化，随着f的变大，其线