电机振动噪声（NVH）气隙磁场推导

最新推荐文章于 2024-05-02 21:39:04 发布

夭夭耀

最新推荐文章于 2024-05-02 21:39:04 发布

阅读量2.7k

点赞数 6

分类专栏：电机震动噪声文章标签：经验分享 1024程序员节

本文链接：https://blog.csdn.net/weixin_43210097/article/details/120911573

版权

电机震动噪声专栏收录该内容

6 篇文章 90 订阅

订阅专栏

前言

正文

最后的话

前言

其实我之前振动噪声系列中，涉及气隙磁密的分析的，都是给一个推导的逻辑链，说什么知道磁动势了，知道气隙磁导了，就可以推出气隙磁密了。这话吧，不能说是错的，但是仅仅凭着这个逻辑链，是无法从量上给出确切的气隙磁场的分布情况的。言归正传，之前的笔记旨在给大家搭建一个整体的理论框架，建立起这么一个印象就足够了。因为你拿着之前的笔记去细究，仔细一看其实哪都没说明白，但是振动噪声关心的问题好像又说了很多了。

接下来，我想从源头上，给大家梳理梳理，电机振动噪声那么关心的气隙磁场究竟是怎么推出来的，参考文献肯定还是要的（我拍脑袋也推不出来呀），这篇《Instantaneous Magnetic Field Distribution in Brushless Permanent Magnet dc Motors, Part I: Open-circuit Field》就是参考的主要依据了。以下正文。

推导过程对电磁场的相关知识使用得比较多，建议大家先自己复习一下。

正文

先明确一下讨论的对象，目前只讨论永磁体产生的磁场（Open-circuit Field），而且我们为了简化分析，先讨论表贴式径向充磁且定子无槽的情况下产生的磁场。

其实大家不要把这个问题想太复杂了，本质上跟解应用题是一样的，都是为了对已有的问题列方程去刻画它们，然后解方程，思路就这么简单。很多同学看到整篇整篇的推导，就不想看了，这样是不行的。每一个方程我都会把它的目的写下来，这样比大家自己硬看已经友好很多了。

下面开始列方程了哈，首先我们把视线放在两个地方，一个是气隙，另一个是永磁体。我们的目标很清晰，为了解出气隙磁场，那为啥还要去看永磁体呢，因为它是源，是我们问题的组成部分。列出两个很基础但是很重要.

忽略了空气的相对磁导率，它太接近于1了，第二个式子大家可能看上去比较陌生，M是磁化强度，是描述宏观磁性体磁性强弱程度的物理量。在经典电磁学中，磁化强度(magnetization)或磁性极化( magnetic polarization)是表示磁性物质永久的或者诱发的偶极磁矩的矢量场。通常用符号M表示。由于铁磁物质的M和H不是简单的线性关系，无法进行合并导致的。另外，本文讨论的场量均为矢量。

接着，我们需要用磁位去对方程进行刻画，当然选择标量位了（矢量怎么解微分方程，想想就可怕）。磁位是什么？我带着大家一起推导一遍。首先是磁场环路定律：

$\oint H\cdot dl=\int_{s}^{}J\cdot ds$

如果回路里面没有电流，等式右边就是0了，在这样的条件下，我们的磁场是满足下列方程的：

$rotH=0$

$divB=0$

旋度和散度的符号在线编辑打不出来，大家凑合着看看。同时呢在任何介质中，都有（2）式的存在，所以我们可以得到下面的关系：

$divH=-divM$

然后，我们假想一个磁荷密度出来（主要是为了跟电荷密度对应）公式编辑器确实不太方便，我用图片的方式给出来：

有了这个式子，我们可以分别知道H的散度和旋度了。我把散度用磁荷密度来表达一下，因为之后会用到这个式子对气隙磁场进行求解：

然后，引入磁标势 $\varphi$ （标量位）来刻画磁场了，定义如下：

$H=-gra\varphi$

接着，再把话题切回来，继续求解我们的气隙磁场。用（3）式分别刻画气隙磁场和永磁体磁场，在气隙中，有：

$gra^{2}\varphi_{air} =0$

$gra^{2}\varphi =\frac{divM}{\mu _{0}}$

其实得到这两个表达式，我们的方程列写环节就已经结束了。但是，只给这两个式子，这个方程组压根解不了。下面做的事情就是选取合适的数学工具进行推导了。还记得我们的永磁体是怎样的充磁方式吗？极坐标就是为此而生的，在极坐标下，有如下关系：

$\dot{M}=M1\dot{R}+M2\dot{\theta }$

R和 $\theta$ 分别是极坐标下的单位矢量，接着我们再结合永磁体径向充磁的特点，有如下的关系：

$M1=\sum M_{n}\cos (np\theta)$

$M2=0$

在极坐标下对M求旋度，直接给结果了哈：

$divM=\frac{M1}{r}$

这样一顿操作之后，前面的两个关于磁标量位的表达式就变成了极坐标下的相应表达（极坐标的表达式其实就是柱坐标系去掉Z轴），如下所示：

上方分别是在气隙中和在永磁体中的磁标量关于径向长度和偏角的偏微分方程，再加上边界条件进行约束，就可以把两个磁标量位给求解出来，进而求得H，再进一步将磁场的B求出来。下图是我们对求解磁场给定的边界条件：

真的是不得不加一个水平线分割一下，看到这大家好多都晕了。但是目前为止，真的只是将这个问题转化到数学的维度上，是的，我们还没开始对问题进行求解。只是用数学去表述了这个问题。至于求解的过程，我目前没有能力对这个偏微分方程组进行手算。的确是能力有限了，之后我会跟我的同伴讨论一下这个问题，期待着后文吧。

最后的话

我记录跟推导无关的一些题外话，电机振动噪声的仿真特点是跨多个领域进行，从电磁到振动再到声音的传播，仿真链长，时间长，传递误差大。最让人困惑的是，仿真结果的准确性如何验证，在实际运行的过程中，电机连接件产生的机械噪声与电磁噪声难以分离，这一问题使得电机的振动噪声没有一个量化的标准，更别说通过耦合仿真来准确地得到有意义的振动噪声仿真数据了。

也是基于上方描述的一些现象和无奈，目前学界对电磁环节电磁力波的研究力度比较大，因为单单是研究这个点，似乎还处于我们能够去思考琢磨的范围，从电磁到振动，是真的很复杂很复杂。我并不是对这个问题感到悲观，只是有些时候面对这些问题，经常会产生一种无力感。

所以呢，大家不要讨厌这满篇的推导过程，它很可能是面对这复杂问题时，我们用自己的智慧和时间能够做到的一些努力。

夭夭耀

关注

6
点赞
踩
19

收藏

觉得还不错? 一键收藏
4
评论
电机振动噪声（NVH）气隙磁场推导

其实我之前振动噪声系列中，涉及气隙磁密的分析的，都是给一个推导的逻辑链，说什么知道磁动势了，知道气隙磁导了，就可以推出气隙磁密了。这话吧，不能说是错的，但是就跟大海一样，全是水。言归正传，之前的笔记旨在给大家搭建一个整体的理论框架，建立起这么一个印象就足够了。因为你拿着之前的笔记去细究，仔细一看其实哪都没说明白，但是振动噪声关心的问题好像又说了很多了。接下来，我想从源头上，给大家梳理梳理，电机振动噪声那么关心的气隙磁场究竟是怎么推出来的
复制链接

扫一扫