电动机发电机的参数灵敏度分析 步骤一基于有限元法采集数据 步骤二基于Matlab程序进行参数灵敏度分析 步骤三分析结果绘图第一步有限元暗房操作用ANSYS Maxwell搭个永磁同步电机模型重点盯着磁钢厚度、气隙长度、绕组匝数三个刺头参数。脚本控制参数批量跑仿真偷摸记录下空载反电势和铁损数据# ANSYS Maxwell脚本片段 oDesign.ChangeProperty( [ NAME:AllTabs, [ NAME:LocalVariableTab, [NAME:PropServers, LocalVariables], [NAME:ChangedProps,[Name:$磁钢厚度, Value:, f{thk}mm]] ] ] ) # 导出CSV时顺手给文件名打标签 with open(fdata_{thk}_{gap}_{turn}.csv, w) as f: f.write(output_data)第二步Matlab猎杀时刻数据灌进Matlab上sobol法全局灵敏度分析。注意标准化操作别翻车% 敏感度刺客在线教学 data readtable(妖艳数据合集.csv); X normalize([data.磁钢厚度, data.气隙, data.匝数]); Y normalize(data.反电势); [sobol_indices, st_error] sobol_sensitivity(X, Y, order,2); % 交互作用检测器 interaction_effect sobol_indices(:,2) - sum(sobol_indices(:,1),2);这里sobol_indices第一列是各参数单干的影响力第二列是狼狈为奸的合作效果。突然发现气隙和绕组匝数居然有负交互作用——这俩参数同时增大会互相拖后腿电动机发电机的参数灵敏度分析 步骤一基于有限元法采集数据 步骤二基于Matlab程序进行参数灵敏度分析 步骤三分析结果绘图第三步可视化摊牌用热力图暴露参数间的爱恨纠葛% 热力修罗场 heatmap({磁钢,气隙,匝数}, {磁钢,气隙,匝数}, interaction_effect); colormap(jet); title(参数相爱相杀关系图); % 单参数影响力排行榜 bar(sobol_indices(:,1),FaceColor,[0.2 0.6 0.8]); set(gca,XTickLabel,{磁钢Thk,气隙Gap,匝数Turns});结果图显示气隙长度以0.78的敏感度称王但它的交互作用全是负值——提醒我们调参时得盯着其他参数同步配平。磁钢厚度看似存在感不强却在铁损分析中突然黑化这告诉我们敏感度这玩意得看具体优化目标跟找对象一样不能只看一面。小贴士有限元数据别傻乎乎跑全参数组合Latin超立方采样能让计算量腰斩。Matlab并行计算工具箱开满咖啡还没凉结果就出来了。