ARTICLE DETAIL

资讯详情

深耕网站建设与运营推广的一线实战洞察。

ANSYS Maxwell轴向磁通电机3D建模:RMxprt参数映射与实战

ANSYS Maxwell轴向磁通电机3D建模:RMxprt参数映射与实战 1. 轴向磁通电机建模为什么值得单独拎出来讲轴向磁通电机这几年在新能源车、无人机、工业伺服领域的热度一直往上走核心原因就一个同样体积下它的功率密度和转矩密度比传统径向磁通电机高出一截。磁路走的是轴向气隙是平面定转子像两个盘子面对面扣在一起铁芯利用率天然就高。但问题也恰恰出在这个“盘子”结构上——二维模型根本描述不了它的端部效应和三维磁路分布你拿Maxwell 2D去算结果和实测能差出百分之二三十这不是仿真精度的问题是模型维度选错了。所以做轴向磁通电机绕不开3D建模。而ANSYS Maxwell这套工具链里RMxprt负责的是“快速路”你输入额定功率、转速、极槽配合、绕组形式这些顶层参数它帮你把定转子尺寸、绕组匝数、线径这些几何量自动算出来生成一个参数化的电机模型。但RMxprt原生输出的是径向磁通电机模板轴向磁通的结构它不直接支持。这就带来一个很现实的工程问题怎么把RMxprt算出来的电磁参数高效地迁移到3D轴向磁通模型里而不是从零开始画几何。我自己的做法是把RMxprt当成“电磁参数计算器”来用而不是“几何生成器”。具体来说先用RMxprt的通用电机模板比如三相感应电机或永磁同步电机输入目标性能指标让它算出每极磁通、电负荷、磁负荷、绕组匝数这些关键量。然后拿这些量去驱动3D轴向磁通模型的参数化建模。这样做的好处是电磁方案有理论依据不是拍脑袋定尺寸同时3D模型又能真实反映轴向磁通的三维磁路两边各取所长。这个思路适合谁如果你已经会用Maxwell 2D做径向磁通电机仿真现在要转轴向磁通那这套流程能帮你省掉大量试错时间。如果你是完全新手建议先把Maxwell 2D的瞬态场求解流程跑通一遍再来看3D否则很容易被三维网格和求解时间劝退。2. RMxprt参数提取与轴向磁通模型的映射逻辑2.1 RMxprt里到底该看哪些输出量打开RMxprt选一个最接近的电机类型模板。轴向磁通永磁同步电机通常用“Adjust-Speed Synchronous Machine”或者“Brushless Permanent-Magnet DC Machine”模板来近似。输入额定功率、额定转速、直流母线电压、极对数、定子外径这些顶层参数求解之后重点看这几组输出每极磁通在“Winding”或“Magnetic”结果里找这个值直接决定你3D模型里永磁体的厚度和面积。电负荷A单位安培导体数反映绕组设计的电流密度水平。磁负荷B气隙磁密幅值轴向磁通电机一般取0.6~0.8T太高了铁耗压不住。绕组匝数与线径RMxprt会给出每相串联匝数和并绕根数这是你3D模型里绕组几何的直接输入。槽满率这个值决定了你3D模型里槽形尺寸能不能塞下绕组。注意RMxprt算出来的尺寸是基于径向磁通结构的定子轭部厚度、齿宽这些几何量不能直接照搬到轴向磁通模型里。你要提取的是电磁参数不是几何参数。2.2 从电磁参数反推轴向磁通几何尺寸拿到每极磁通Φ和磁负荷B之后永磁体的有效面积A_m可以用A_m Φ / (B * k_leak)来估算k_leak是漏磁系数轴向磁通电机端部漏磁比径向磁通大一般取1.15~1.3。永磁体厚度h_m则根据去磁曲线和工作点来定简单估算可以用h_m (B_r - B) * l_g / (μ_0 * H_c)其中B_r是剩磁l_g是气隙长度H_c是矫顽力。定子齿宽和轭部厚度用磁通连续性来算齿宽b_t Φ / (B_t * k_fe * L_stk)B_t是齿部磁密一般取1.6~1.8Tk_fe是叠压系数L_stk是轴向有效长度。轭部厚度类似但轭部磁通是齿部磁通的一半因为磁通在轭部分流所以b_y Φ / (2 * B_y * k_fe * L_stk)。这些公式在RMxprt的输出里都能找到对应的输入量你不需要手算每一个但理解这个映射关系很重要——它让你知道3D模型里改一个尺寸电磁性能会往哪个方向变。2.3 参数化建模的变量表设计在Maxwell 3D里建轴向磁通模型强烈建议用参数化变量来驱动几何。我一般会定义这么一组变量变量名含义初始值来源D_out定子外径RMxprt定子外径D_in定子内径按D_out的0.5~0.6倍初选h_m永磁体厚度上述公式估算g气隙长度0.5~1.0mm按加工精度定b_t定子齿宽磁通连续性计算b_y定子轭厚磁通连续性计算N_t每相匝数RMxprt输出I_rms相电流有效值按功率和电压反推把这些变量在Maxwell的“Design Properties”里定义好后面画几何的时候直接引用变量名改一个值整个模型自动更新。这个习惯能让你在优化阶段省掉大量重复劳动。3. 3D轴向磁通模型的具体搭建步骤3.1 几何建模从盘式结构到参数化草图轴向磁通电机的基本结构是中间一个定子盘两边各一个转子盘或者反过来中间转子两边定子。以双定子单转子结构为例建模顺序如下画定子盘在XY平面上画一个圆环外径D_out内径D_in拉伸厚度L_stk。这个圆环就是定子铁芯的轭部。开槽在定子盘的一个端面上画槽形。轴向磁通电机的槽通常是梯形或平行槽槽底宽度和槽口宽度不一样。用参数化草图槽宽用b_s表示槽深用h_s表示槽数用N_s表示。然后用“Array”功能沿圆周阵列。放绕组在槽内画绕组截面截面面积由线径和并绕根数决定。Maxwell里可以用“Winding”功能自动生成绕组但3D模型里我更喜欢手动画导体因为端部形状对电感影响很大。画永磁体在转子盘表面画扇形永磁体扇形角度等于极距角乘以极弧系数。永磁体厚度h_m充磁方向沿轴向。画转子轭永磁体背面加一层转子轭材料用低碳钢厚度按磁通连续性算。设置气隙定子盘和永磁体之间留出气隙g这个气隙在3D模型里是空气域网格要加密。实操心得画扇形永磁体的时候用“Draw Rectangle”然后“Sweep Along Arc”比直接画扇形更可控因为你可以精确控制永磁体的内外半径和角度。阵列的时候用“Duplicate Around Axis”角度步长设为360/N_pN_p是极数。3.2 材料定义与边界条件设置材料这块有几个坑要避开定子铁芯用“steel_1008”或者自定义的硅钢片材料BH曲线要输入完整。轴向磁通电机频率高铁耗系数要设对否则涡流损耗算不准。永磁体钕铁硼材料剩磁B_r和矫顽力H_c按牌号输入。充磁方向用“Cylindrical”坐标系方向沿Z轴。绕组铜电导率5.8e7 S/m。如果考虑涡流效应要勾选“Eddy Effects”。气隙空气相对磁导率1。边界条件方面轴向磁通电机通常用“Master/Slave”边界来模拟周期性。如果模型是1/8周期比如8极电机取45度在切面上设Master和Slave相位差设0因为磁场是周期的。外边界用“Balloon”或者“Zero Tangential H”都行Balloon更省网格。3.3 网格划分策略与求解设置3D瞬态场的网格量是最大的瓶颈。我的经验是气隙网格沿轴向至少3层每层厚度不超过0.3mm。用“Inside Selection”加“Length Based”控制。永磁体网格表面网格尺寸不超过2mm因为涡流主要在表面。定子齿部齿宽方向至少2个单元否则磁密峰值会被平均掉。绕组如果只算电磁绕组可以用粗网格如果要算温升绕组网格要细。求解设置里瞬态场的时间步长很关键。步长太大转子位置更新不及时转矩波形会失真。一般取一个电周期至少100步。比如转速3000rpm8极电机电频率200Hz电周期5ms步长取50μs。总求解时间至少3个电周期让瞬态响应衰减掉。注意Maxwell 3D瞬态场求解时间很长一个模型跑几个小时很正常。建议先用2D近似模型验证绕组连接和激励设置确认无误再跑3D。4. 常见报错与排查技巧实录4.1 RMxprt导出模型失败怎么办RMxprt导出到Maxwell 2D/3D的时候最常见的报错是“failover feature ansys electronics_desktop is not available”。这个通常不是license问题而是RMxprt的设计文件损坏了。解决办法新建一个RMxprt项目把原来的输入参数重新填一遍不要直接复制旧项目。如果还不行检查一下安装目录下的“syslib”文件夹里有没有对应的材料库文件。另一个常见问题是“maxwell中rm如何设置输出几分之几的模型”。RMxprt默认输出全模型要输出部分模型在“Machine”选项卡里找到“Fraction of Model”填1/2、1/4或者1/8。但注意这个设置只影响几何输出不影响电磁计算——RMxprt内部还是按全模型算的。4.2 3D模型求解不收敛的排查顺序3D瞬态场不收敛按这个顺序查检查绕组激励电流方向对不对三相相位差是不是120度用“Winding”功能自动生成的话检查一下“Polarity”有没有设反。检查永磁体充磁方向相邻永磁体极性要交替用“Cylindrical”坐标系的时候角度偏移要设对。检查网格质量用“Mesh Statistics”看有没有负体积单元。有的话把几何里的尖角倒圆或者把网格尺寸调小。检查时间步长步长太大导致转子“跳过”了某些位置磁场突变引起不收敛。把步长减半试试。检查材料BH曲线如果BH曲线在低磁场区斜率太大迭代容易震荡。把BH曲线在0~0.5T区间加密。4.3 电感计算偏差大的修正方法Maxwell算出来的电感经常和实测对不上尤其是轴向磁通电机端部漏磁大2D模型根本算不准。修正方法用3D模型算电感3D能捕捉端部漏磁但网格要够细。加端部绕组在3D模型里把绕组端部画出来不要用“End Connection”简化。用能量法验证在“Parameters”里勾选“Matrix”算完看电感矩阵。如果自感和互感比例不对检查绕组相序。实操心得我一般会先用RMxprt算一个电感初值然后3D模型算一个值两个值差在15%以内就认为合理。如果差太多优先检查气隙长度和永磁体厚度这两个量对电感影响最大。5. 从仿真到实测的闭环验证思路5.1 空载反电动势对比空载反电动势是验证模型最直接的指标。仿真里把绕组开路转子按额定转速转看相电压波形。实测用示波器测反电动势对比峰值和波形畸变率。如果仿真峰值比实测高10%以上大概率是永磁体工作点设高了把B_r调低5%再试。如果波形畸变率差很多检查极弧系数和槽口宽度。5.2 负载转矩与效率Map负载转矩对比要在相同电流下进行。仿真里给定额定电流看平均转矩。实测用测功机加载记录转矩-转速曲线。效率Map的话仿真可以扫电流角和超前角实测要逐点测工作量很大。我的做法是仿真扫一个粗Map找出高效区实测只测高效区里的几个点验证趋势就行。5.3 温升对电磁性能的修正轴向磁通电机散热条件比径向磁通差温升对永磁体性能影响很大。仿真里如果只算冷态转矩会偏高。修正方法用Maxwell和Icepak做耦合或者简单点把永磁体温度设到实际工作温度比如80°CB_r按温度系数修正。钕铁硼的温度系数大约是-0.11%/°C80°C时B_r比20°C时低6.6%左右。6. 参数化优化与脚本自动化6.1 用Optimetrics做多目标优化Maxwell自带的Optimetrics模块可以做参数扫描和优化。我一般会定义这几个优化变量永磁体厚度h_m、气隙g、定子齿宽b_t、绕组匝数N_t。目标函数设两个平均转矩最大化和转矩脉动最小化。约束条件槽满率不超过0.75磁密不超过1.8T。优化算法选“Quasi Newton”或者“Genetic Algorithm”。Quasi Newton收敛快但容易陷局部最优Genetic Algorithm全局搜索能力强但计算量大。我的做法是先用Genetic Algorithm跑20代找到大致区域再用Quasi Newton精细搜索。6.2 用IronPython脚本批量处理Maxwell支持IronPython脚本可以批量改参数、跑仿真、导出结果。比如你要扫10个气隙值手动改太慢写个脚本import ScriptEnv ScriptEnv.Initialize(Ansoft.ElectronicsDesktop) oProject oDesktop.GetActiveProject() oDesign oProject.GetActiveDesign() for g in [0.5, 0.6, 0.7, 0.8, 0.9, 1.0]: oDesign.ChangeProperty([NAME:AllTabs, [NAME:LocalVariableTab, [NAME:PropServers, LocalVariables], [NAME:ChangedProps, [NAME:g, Value:, str(g) mm]]]]) oDesign.Analyze(Setup1) oModule oDesign.GetModule(ReportSetup) oModule.CreateReport(Torque_ str(g), Transient, Rectangular Plot, Setup1 : Transient, [Domain:, Sweep], [Time:, [All]], [Torque:, [Moving1.Torque]])这个脚本会把每个气隙值下的转矩波形导出成单独的报告。跑完一轮你就有完整的数据了。6.3 脚本自动化的注意事项变量名要统一脚本里引用的变量名必须和Design Properties里定义的一模一样大小写敏感。求解器要提前设好脚本只改参数不建设置。求解步长、收敛精度这些要在脚本运行前手动设好。结果导出路径默认导出到项目目录下如果要指定路径在CreateReport里加“OutputPath”参数。错误处理脚本跑批量的时候某个参数组合不收敛会中断整个循环。加个try-except把不收敛的组合跳过记录到日志里。7. 我踩过的坑和最后分享的几个技巧第一个坑RMxprt里设了1/8模型导出到3D之后发现几何只有45度但边界条件没自动加Master/Slave结果磁场全乱了。后来才知道RMxprt的分数模型只影响几何边界条件要手动加。第二个坑3D模型跑了一晚上结果出来转矩是负的。查了半天发现是绕组相序设反了A相和C相对调了一下就好了。所以跑长仿真之前一定先用短时间步长跑一个电周期看转矩方向对不对。第三个坑永磁体用“Cylindrical”坐标系充磁角度偏移没设对相邻磁极极性一样气隙磁密直接掉了一半。后来改成用“Cartesian”坐标系每个永磁体单独设充磁方向虽然麻烦但不会错。最后分享一个技巧Maxwell 3D的网格划分很吃时间如果只是做参数扫描可以把模型切成1/4或者1/8用Master/Slave边界网格量能降到1/4以下求解时间大幅缩短。但切模型的时候注意切面要选在磁通对称的位置否则边界条件加不对。还有一个技巧如果只需要看电磁性能不需要看温升可以把绕组用“Stranded”导体代替“Solid”导体求解速度快很多。Stranded导体不考虑涡流但轴向磁通电机的绕组涡流损耗占比不大忽略掉对转矩影响很小。这个流程我前后跑了大概十几个轴向磁通模型从最初的RMxprt导出失败到后来能稳定跑通3D瞬态场中间花了不少时间在排查各种报错上。现在回头看最关键的还是把RMxprt的电磁参数和3D几何的映射关系理清楚这个逻辑通了后面就是调参数和等结果的事了。
返回列表