直线电机COMSOL仿真案例

本文介绍了一种改进的COMSOL官方二维模型,用于仿真有限长度直线电机的电磁特性。通过与ANSYSMaxwell模型的对比,证实了这种方法的可行性。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

0 Preface

COMSOL官方直线电机二维模型案例应用的是“周期性条件”特征,无法准确评估有限长度直线电机的电磁特性。因此本例结合官方案例探索出一种有限长度直线电机的电磁仿真方法。

特别地,本案例方法中参数不一定合理,设置不一定完全正确,方法也不是最优方法。但是本案例已通过ANSYS Maxwell建立相同模型进行了仿真验证,两者电磁性能基本一致,因此可以保证的是本案例方法是可行的,可以提供一种参考思路。至于是不是完全准确,那就再说吧…… 💨

该案例仅供学习使用,请勿传播。除此以外,还可以参考视频教程

(有需要模型的联系我:narcos@foxmail.com)

链接: https://caiyun.139.com/m/i?125Cmp6gLasAl  提取码:XUGi

00 准备

1 建模操作说明

建议配合视频步骤

1.1 新建

新建窗口中, 单击模型向导

1 二维。

2 AC/DC> 电磁场 > 磁场 (mf)。

3 稳态。

1.2 几何1

单位 mm

(此部分采用DXF导入建模,所以涉及到提取、并集和装配操作,十分繁琐。下述步骤较为简单,建议参考视频步骤一步一步来)

1.2.1 导入 (imp1)

1 右键几何选择导入。

2 导入 linearmotor.dxf 。

1.2.2 stator1 (extract1)

1 右键选择提取 (extract)。

2 几何实体层选择域。

3 选择上部分铁心。

4 输入对象处理选择创建其余对象。

5 构建选定对象。

(下述相同方法)

1.2.3 stator2 (extract2)

1.2.4 pms (extract3)

1.2.5 coils1 (extract4)

1.2.6 coils2 (extract5)

1.2.7 statorair1 (extract6)

1.2.8 statorair2 (extract7)

1.2.9 mapair1 (extract8)

1.2.10 mapair2 (extract9)

1.2.11 moverair1 (extract10)

1.2.12 moverair2 (extract11)

1.2.13 moverair (extract12)

1.2.14 Stator1 (uni1)

1.2.15 Stator2 (uni2)

1.2.16 Mover (uni3)

1.2.17 形成装配 (fin)

1.3 定义

此部分可以方便后续选择域,也可以忽略。

1.3.1 选择(显式)

  • PM_N

域 55, 57, 59, 61, 63

  • PM_S

域 56, 58, 60, 62, 64

  • Iron

域 1, 67

  • U1

域 68–71, 80–83, 92–95, 104–107

  • V1

域 72–75, 84–87, 96–99, 108–111

  • W1

域 76–79, 88–91, 100–103, 112–115

  • U2

域 3–6, 15–18, 27–30, 39–42

  • V2

域 7–10, 19–22, 31–34, 43–46

  • W2

域 11–14, 23–26, 35–38, 47–50

  • mapedge1

边界 255–256

  • mapedge2

边界 305, 308

1.4 材料

不同域给定不同材料

1.4.1 Air (mat1)

几何实体层
选择几何 geom1: 维度 2: 所有域

1.4.2 Copper (mat2)

几何实体层
选择几何 geom1: 维度 2: 域 3–50, 68–115

1.4.3 Soft Iron (Without Losses) (mat3)

几何实体层
选择几何 geom1: 维度 2: 域 1, 67

1.4.4 BMN-48SH (mat4)

几何实体层
选择几何 geom1: 维度 2: 域 55–64

1.5 定义

给定仿真参数

1.5.1 参数

名称表达式描述
v2[m/s]2 m/s
t0[s]0 s
airgap1[mm]0.001 m
Lef55[mm]0.055 m
ipp11[A]11 A
tau10[mm]0.01 m
fsv/2/tau100 1/s

1.5.2 变量1

名称表达式单位描述
iaipp*sin(2*pi*fs*t)A
ibipp*sin(2*pi*fs*t - 2*pi/3)A
icipp*sin(2*pi*fs*t + 2*pi/3)A

1.5.3 组件1变量2(可省略)

名称表达式单位描述
F_phimf.Bx*mf.By/mu0_constPa
Force1intop1(F_phi)/(airgap/4)*LefN
Force2intop2(F_phi)/(airgap/4)*LefN

1.5.4 非局部耦合(可省略)

  • 积分1
几何实体层
选择几何 geom1: 维度 2: 域 54

  • 积分2
几何实体层
选择几何 geom1: 维度 2: 域 52

1.6 动网格1

此部分参考电磁柱塞

1.6.1 指定变形1

几何实体层
选择几何 geom1: 维度 2: 域 52–64

指定变形

描述
指定变形{v*t, 0, 0}

1.6.2 变形域1

几何实体层
选择几何 geom1: 维度 2: 域 51, 65

平滑处理

描述
网格平滑类型拉普拉斯

1.6.3 固定边界1

几何实体层边界
选择几何 geom1: 维度 1: 边界 254, 309

1.6.4 指定法向网格位移1

几何实体层边界
选择几何 geom1: 维度 1: 边界 255–256, 305, 308

1.7 磁场(mf)

1.7.1 安培定律1(无需设置)

1.7.2 磁绝缘1(无需设置)

1.7.3 初始值1(无需设置)

1.7.4 Iron

几何实体层
选择几何 geom1: 维度 2: 域 1, 67

本构关系 B-H Settings

DescriptionValue
磁化模型B - H 曲线
磁场模来自材料
磁共能密度来自材料

1.7.5 PMs

几何实体层
选择几何 geom1: 维度 2: 域 55–64

磁体Settings

DescriptionValue
方向法指定南北边界
图案类型线性图案(或基于域索引)
周期性类型交变

本构关系 B-H Settings

DescriptionValue
磁化模型剩余磁通密度
回复磁导率来自材料
剩余磁通密度模来自材料

  • 北1
几何实体层边界
选择几何 geom1: 维度 1: 边界 265

  • 南1
几何实体层边界
选择几何 geom1: 维度 1: 边界 266

1.7.6 U1

几何实体层
名称U1
选择命名 sel4: 几何 geom1: 维度 2: 域 68–71, 80–83, 92–95, 104–107

线圈Settings

DescriptionValueUnit
线圈名称U1
导线模型均匀多匝
线圈组
线圈激励电流
线圈电流iaA

均匀多匝导线Settings

DescriptionValueUnit
匝数10
线圈导线电导率6E7S/m
线圈导线截面积用户定义
线圈导线截面积1E-6

  • 反向电流方向1
几何实体层
选择几何 geom1: 维度 2: 域 69–70, 80, 83, 93–94, 104, 107

1.7.7 V1

几何实体层
名称V1
选择命名 sel5: 几何 geom1: 维度 2: 域 72–75, 84–87, 96–99, 108–111

线圈Settings

DescriptionValueUnit
线圈名称V1
导线模型均匀多匝
线圈组
线圈激励电流
线圈电流ibA

均匀多匝导线Settings

DescriptionValueUnit
匝数10
线圈导线电导率6E7S/m
线圈导线截面积用户定义
线圈导线截面积1E-6

  • 反向电流方向1
几何实体层
选择几何 geom1: 维度 2: 域 72, 75, 85–86, 96, 99, 109–110

1.7.8 W1

几何实体层
名称W1
选择命名 sel6: 几何 geom1: 维度 2: 域 76–79, 88–91, 100–103, 112–115

线圈Settings

DescriptionValueUnit
线圈名称W1
导线模型均匀多匝
线圈组
线圈激励电流
线圈电流icA

均匀多匝导线Settings

DescriptionValueUnit
匝数10
线圈导线电导率6E7S/m
线圈导线截面积用户定义
线圈导线截面积1E-6

  • 反向电流方向1
几何实体层
选择几何 geom1: 维度 2: 域 77–78, 88, 91, 101–102, 112, 115

1.7.9 U2

几何实体层
名称U2
选择命名 sel7: 几何 geom1: 维度 2: 域 3–6, 15–18, 27–30, 39–42

线圈Settings

DescriptionValueUnit
线圈名称U2
导线模型均匀多匝
线圈组
线圈激励电流
线圈电流iaA

均匀多匝导线Settings

DescriptionValueUnit
匝数10
线圈导线电导率6E7S/m
线圈导线截面积用户定义
线圈导线截面积1E-6

  • 反向电流方向1
几何实体层
选择几何 geom1: 维度 2: 域 4–5, 15, 18, 28–29, 39, 42

1.7.10 V2

几何实体层
名称V2
选择命名 sel8: 几何 geom1: 维度 2: 域 7–10, 19–22, 31–34, 43–46

线圈Settings

DescriptionValueUnit
线圈名称V2
导线模型均匀多匝
线圈组
线圈激励电流
线圈电流ibA

均匀多匝导线Settings

DescriptionValueUnit
匝数10
线圈导线电导率6E7S/m
线圈导线截面积用户定义
线圈导线截面积1E-6

  • 反向电流方向1
几何实体层
选择几何 geom1: 维度 2: 域 7, 10, 20–21, 31, 34, 44–45

1.7.11 W2

几何实体层
名称W2
选择命名 sel9: 几何 geom1: 维度 2: 域 11–14, 23–26, 35–38, 47–50

线圈Settings

DescriptionValueUnit
线圈名称W2
导线模型均匀多匝
线圈组
线圈激励电流
线圈电流icA

均匀多匝导线Settings

DescriptionValueUnit
匝数10
线圈导线电导率6E7S/m
线圈导线截面积用户定义
线圈导线截面积1E-6

  • 反向电流方向1
几何实体层
选择几何 geom1: 维度 2: 域 12–13, 23, 26, 36–37, 47, 50

1.7.12 连续性1

几何实体层边界
选择几何 geom1: 维度 1: 所有边界

对选择Settings

DescriptionValue
一致边界对 1 (ap1)

1.7.13 连续性2

几何实体层边界
选择几何 geom1: 维度 1: 所有边界

对选择Settings

DescriptionValue
一致边界对 2 (ap2)

1.7.14 力计算1

几何实体层
选择几何 geom1: 维度 2: 域 55–64

力计算Settings

DescriptionValueUnit
力名称Fe
扭矩轴,x 分量0
扭矩轴,y 分量0
扭矩轴,z 分量1
扭矩旋转点,x 分量0m
扭矩旋转点,y 分量0m
扭矩旋转点,z 分量0m

1.8 网格1

1.8.1 大小 (size)

描述
最大单元大小13.7
最小单元大小0.0612
曲率因子0.3
最大单元增长率1.3

1.8.2 Iron (size1)

几何实体层
选择几何 geom1: 维度 2: 域 1, 67

描述
最小单元大小0.0612
最小单元大小
曲率因子0.3
曲率因子
狭窄区域分辨率
最大单元增长率1.3
最大单元增长率
定制单元大小定制

1.8.3 Coil (size2)

几何实体层
选择几何 geom1: 维度 2: 域 3–50, 68–115

描述
最小单元大小0.0612
最小单元大小
曲率因子0.3
曲率因子
狭窄区域分辨率
最大单元增长率1.3
最大单元增长率
定制单元大小定制

1.8.4 statorair (size3)

几何实体层
选择几何 geom1: 维度 2: 域 2, 66

描述
最大单元大小0.5
最小单元大小0.0612
最小单元大小
曲率因子0.3
曲率因子
狭窄区域分辨率
最大单元增长率1.3
最大单元增长率
定制单元大小定制

1.8.5 moverair (size4)

几何实体层
选择几何 geom1: 维度 2: 域 52–54

描述
最小单元大小0.0612
最小单元大小
曲率因子0.3
曲率因子
狭窄区域分辨率
最大单元增长率1.3
最大单元增长率
定制单元大小定制

1.8.6 PMs (size5)

几何实体层
选择几何 geom1: 维度 2: 域 55–64

描述
最大单元大小0.6
最小单元大小0.0612
最小单元大小
曲率因子0.3
曲率因子
狭窄区域分辨率
最大单元增长率1.3
最大单元增长率
定制单元大小定制

1.8.7 映射1 (map1)

这一步比较关键

几何实体层
选择几何 geom1: 维度 2: 域 51, 65

描述
迭代次数4
要处理的最大单元深度4
  • 分布1 (dis1)
几何实体层边界
选择几何 geom1: 维度 1: 边界 254, 309

描述
单元数20
  • 分布2 (dis2)
几何实体层边界
名称mapedge1
选择命名 sel10: 几何 geom1: 维度 1: 边界 255–256

描述
分布类型预定义
单元数50
单元大小比5
  • 分布3 (dis3)
几何实体层边界
名称mapedge2
选择命名 sel11: 几何 geom1: 维度 1: 边界 305, 308

描述
分布类型预定义
单元数50
单元大小比5
反向

1.8.8 自由三角形网格1 (ftri1)

几何实体层
选择剩余部分

描述
迭代次数4
要处理的最大单元深度4

2 研究2

添加一个瞬态的研究

2.1 瞬态

时步单位
range(0,1/fs/100,2/fs)s

2.2 求解器配置

瞬态求解器1>全耦合1

方法和终止

描述
雅可比矩阵更新在每次迭代中
最大迭代次数25
容差因子0.2

开始计算

3 结果

3.1 磁通密度模 (mf)

3.2 Thrust

1 右键结果选择一维绘图组

2 命名为Thrust

3 设置轴坐标、图例等

4 右键Thrust选择全局

5 添加表达式mf.Forcex_Fe

6 绘图

3.3 VCoil1

3.4 VCoil2

3.5 Fluxlinkage1

3.6 Fluxlinkage2

3.7 Current1

3.8 Current2

3.9 Force_cal(可忽略)

该方法为取气隙中的电磁力波,然后通过麦克斯韦张力法计算得到电磁力

3.10 Force_cal_contrast(可忽略)

两种电磁力计算方法对比

4 ANSYS Maxwell对比

基于相同的模型和激励设置

4.1 模型

仿真步骤不一一给出,可参考开源仿真文件

4.2 结果

推力波形周期一致,值差不多,因此COMSOL按照此案例方法可行。


comsol模型: Download - MEGA

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值