CN110210116B - 一种适用于带集磁器管件电磁成形的二维轴对称仿真模型及建模方法 - Google Patents
一种适用于带集磁器管件电磁成形的二维轴对称仿真模型及建模方法 Download PDFInfo
- Publication number
- CN110210116B CN110210116B CN201910465106.1A CN201910465106A CN110210116B CN 110210116 B CN110210116 B CN 110210116B CN 201910465106 A CN201910465106 A CN 201910465106A CN 110210116 B CN110210116 B CN 110210116B
- Authority
- CN
- China
- Prior art keywords
- conductor ring
- conductor
- workpiece
- close
- magnetic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000004088 simulation Methods 0.000 title claims abstract description 25
- 239000004020 conductor Substances 0.000 claims abstract description 103
- 239000000463 material Substances 0.000 claims abstract description 29
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 239000011810 insulating material Substances 0.000 claims abstract description 6
- 238000009826 distribution Methods 0.000 claims description 21
- 230000008569 process Effects 0.000 claims description 8
- 230000004907 flux Effects 0.000 claims description 2
- 230000035699 permeability Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 12
- 238000005516 engineering process Methods 0.000 description 5
- 230000002708 enhancing effect Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 229910052782 aluminium Inorganic materials 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000017525 heat dissipation Effects 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 238000003913 materials processing Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000007769 metal material Substances 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000003466 welding Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Shaping Metal By Deep-Drawing, Or The Like (AREA)
Abstract
一种适用于带集磁器管件电磁成形的二维轴对称仿真模型及建模方法,采用靠近驱动线圈一侧的导体环、靠近待加工工件一侧的导体环、中间导体环这三个导体环代替集磁器;靠近驱动线圈一侧的导体环和靠近待加工工件一侧的导体环设置为集磁器导磁材料,中间导体环设置为理想绝缘材料;设置靠近驱动线圈一侧的导体环中的总电流与靠近待加工工件一侧的导体环中的总电流大小相等、方向相反。本发明一种适用于带集磁器管件电磁成形的二维轴对称仿真模型及建模方法,与现有三维模型相比,本发明提供的仿真模型与方法计算速度更快;与现有二维轴对称模型对比,本发明提供的仿真模型与方法更贴合物理实际,计算精度更高。
Description
技术领域
本发明属于金属材料加工技术领域,特别是一种适用于带集磁器管件电磁成形的二维轴对称仿真模型及建模方法。
背景技术
电磁成形技术作为一种绿色、柔性的高速成形技术,完美契合时代需求并具有显著优势。电磁成形过程中,当待加工工件的电阻率过高时,通常采用带集磁器电磁成形技术;即在驱动线圈与待加工工件之间增加一集磁器,达到增强待加工工件电磁力的目的。此时,因为模型不具有轴对称性,模拟这一过程就显得难度很大。
中国专利“一种新型电磁脉冲焊接用集磁器(CN 105965184B)”公开了一种新的集磁器结构,设置冷却通道,散热性能好,集磁器和线圈使用寿命长;然而该专利并未涉及如何对其进行仿真分析。文献“3D simulation of magnetic field distribution inelectromagnetic forming systems with field-shaper(Journal of materialsprocessing technology,2009(209):2295-2301)”建立了带集磁器管件电磁成形的三维仿真模型,这一模型能够真实的反映带集磁器管件电磁成形工件的电磁力分布与变形过程,但是模型复杂,不利于分析优化。文献“板料电磁成形集磁器工作原理的模拟(塑性工程学报,2011,18(1):36-42)”采用ANSYS二维轴对称模型,研究集磁器的相对直径、厚度、中心圆孔半径等结构参数对工件电磁力分布和强度的影响,并提出提高集磁器效率和改善板料成形效果的途径;然而,这一模型仅单纯的将集磁器看做一个导体环,与物理现实不符,导致其计算精度偏差很大。
电磁成形过程中,当待加工工件的电阻率过高时,通常采用带集磁器电磁成形技术;即在驱动线圈与待加工工件之间增加一集磁器,达到增强待加工工件电磁力的目的。集磁器是一带微小间隙的圆环状导体,因为其整体不具有轴对称性,导致目前带集磁器电磁成形的仿真计算只能采用三维模型;建立三维模型时,因为集磁器存在微小间隙,导致划分网格的时候异常艰难,计算难度非常大,计算精度得不到保障。
发明内容
为解决上述技术问题,本发明提供一种适用于带集磁器管件电磁成形的二维轴对称仿真模型及建模方法,与现有三维模型相比,本发明提供的仿真模型与方法计算速度更快;与现有二维轴对称模型对比,本发明提供的仿真模型与方法更贴合物理实际,计算精度更高。
本发明采取的技术方案为:
一种适用于带集磁器管件电磁成形的二维轴对称仿真模型,采用靠近驱动线圈一侧的导体环、靠近待加工工件一侧的导体环、中间导体环这三个导体环代替集磁器;靠近驱动线圈一侧的导体环和靠近待加工工件一侧的导体环设置为集磁器导磁材料,中间导体环设置为理想绝缘材料;设置靠近驱动线圈一侧的导体环中的总电流与靠近待加工工件一侧的导体环中的总电流大小相等、方向相反。
一种适用于带集磁器管件电磁成形的二维轴对称仿真模型建模方法,包括以下步骤:
步骤1:用与集磁器相同的导体材料,一般为纯铜或纯铝材料。替换集磁器的微小间隙,使集磁器变成一整块具有轴对称型的导体环。
步骤2:将所述导体环分割为三个导体环,分别为:靠近驱动线圈一侧的导体环、靠近待加工工件一侧的导体环、中间导体环。
步骤3:建立二维轴对称模型,包括具有轴对称型的驱动线圈、待加工工件、三个导体环,以及相应的空气域。其中,相应的空气域指的是驱动线圈、待加工工件、三个导体环周围的空气,是驱动线圈的几何尺寸5-10倍。
步骤4:将靠近驱动线圈一侧的导体环内电流密度面积分记为I1,将靠近待加工工件一侧的导体环内电流密度面积分记为I2,设置约束条件I1=-I2。
步骤4中的约束条件设置是为了使集磁器的电流分布能够实现圆环等效。若不做约束条件,模型会默认三个导体环是一个整体,这与集磁器模型不符;设置了约束条件之后,电流会自动在靠近驱动线圈一侧的导体环、靠近待加工工件一侧的导体环中流动,且形成大小相等、方向相反的电流,与集磁器的电流分布一致。
步骤5:将靠近驱动线圈一侧的导体环、靠近待加工工件一侧的导体环的材料设置为集磁器导体材料,将中间导体环的材料设置为理想绝缘材料。其中,理想绝缘材料指的是电导率为0的材料。
步骤6:设置二维轴对称模型中的源、边界及其他材料属性,进行带集磁器电磁成形的二维轴对称计算,得到待加工工件的电磁力分布。
其中,二维轴对称模型中的源是指驱动线圈中的脉冲电流;
二维轴对称模型中的边界指在模型外边界,磁通密度或矢量磁位满足的条件。
二维轴对称模型中的其他材料属性指的是空气域的材料属性,包括空气域的电阻率和磁导率。
图4所示为待加工工件的电磁力分布示意图。
所述靠近驱动线圈一侧的导体环的厚度等于电磁成形过程中的趋肤深度靠近待加工工件一侧的导体环的厚度等于电磁成形过程中的趋肤深度其中,f为电磁成形放电等效频率,μ为集磁器材料磁导率,ε为集磁器材料介电常数。
上述是为了描述靠近驱动线圈一侧的导体环、靠近待加工工件一侧的导体环的厚度,即三个导体环的分割原则。
所述靠近驱动线圈一侧的导体环内电流密度面积分靠近待加工工件一侧的导体环内电流密度面积分其中,J1为靠近驱动线圈一侧的导体环内电流密度,J2为靠近待加工工件的导体环内电流密度,S1为靠近驱动线圈一侧的导体环对应的二维轴对称区域;S2为靠近待加工工件的导体环对应的二维轴对称区域。
上述是为了实现集磁器的电流分布,对导体环内的电流做了约束。保证了集磁器的模型能够简化成二维轴对称。
本发明一种适用于带集磁器管件电磁成形的二维轴对称仿真模型及建模方法,将微小间隙采用与集磁器相同的导体代替,使其变成一整块具有轴对称型的导体环。将带集磁器电磁成形模型简化为二维轴对称模型,达到简化计算带集磁器电磁成形模型的目的。
与现有三维模型相比,本发明提供的仿真模型与方法维数降低,计算量大大减小,计算速度快。
现有的二维轴对称模型没有做物理意义上的等效,不科学,故精度很差。与现有二维轴对称模型对比,本发明提供的仿真模型与方法更贴合物理实际,计算精度更高。
附图说明
图1(a)为带集磁器管件电磁成形的原始三维几何示意图;
图1(b)为采用一具有轴对称型的导体环代替集磁器的电磁成形三维几何示意图;
图1(c)为将导体环分割为三个导体环的电磁成形三维几何示意图;
图1(d)为图1(c)对应的电磁成形二维轴对称模型示意图。
图2(a)为集磁器原始状态时的电流分布示意图。
图2(b)为采用一具有轴对称型的导体环代替集磁器后的电流分布示意图。
图2(c)为将导体环分割为三个导体环后的电流分布示意图。
图2(d)为图2(c)三维几何模型对应的二维轴对称模型电流分布示意图。
图3为适用于带集磁器管件电磁成形的二维轴对称仿真方法流程示意图。
图4为待加工工件的电磁力分布示意图。
其中:1-驱动线圈;2-集磁器;21-集磁器的微小气隙;22-具有轴对称型的导体环;23-靠近驱动线圈一侧的导体环;24-中间导体环;25-靠近待加工工件一侧的导体环;3-待加工工件;4-空气域;51-内层电流;52-外层电流;53-气隙过渡区域电流。
其中:
51-内层电流,指的是集磁器内侧流过的电流。
52-外层电流,指的是集磁器外侧流过的电流。
53-气隙过渡区域电流,指的是集磁器间隙侧流过的电流。
具体实施方式
图1(a)~图1(d)为适用于带集磁器管件电磁成形的二维轴对称仿真模型简化过程示意图。
图2(a)~图2(d)为集磁器电流分布示意图。
图3为适用于带集磁器管件电磁成形的二维轴对称仿真方法流程示意图。
模型简化时,首先,以图1(a)所示的带集磁器管件电磁成形的原始三维几何示为原型,得到集磁器原始状态时的电流分布,如图2(a)所示。
按照图3执行第一步S1:用与集磁器相同的导体材料替换集磁器的微小间隙,使其变成一整块具有轴对称型的导体环,此时采用一具有轴对称型的导体环代替集磁器的电磁成形三维几何示意图如图1(b)所示,采用一具有轴对称型的导体环代替集磁器后的电流分布示意图如图2(b)所示。
执行第二步S2:将导体环分割为三个导体环,分别为靠近驱动线圈一侧的导体环、靠近待加工工件一侧的导体环、中间导体环,此时将导体环分割为三个导体环的电磁成形三维几何示意图如图1(c)所示,将导体环分割为三个导体环后的电流分布示意图如图2(c)所示。
执行第三步S3:采用商业有限元软件,如ansys、comsol,建立二维轴对称模型,包括具有轴对称型的驱动线圈、待加工工件与三个导体环,以及相应的空气域,此时与图1(c)对应的电磁成形二维轴对称模型示意图如图1(d)所示,与图2(c)三维几何模型对应的二维轴对称模型电流分布示意图如图2(d)所示。
执行第四步S4:将靠近驱动线圈一侧的导体环内电流密度面积分记为I1,将靠近待加工工件的导体环内电流密度面积分记为I2,设置约束条件I1=-I2。
执行第五步S5:将靠近驱动线圈一侧的导体环和靠近待加工工件的导体环的材料设置为集磁器导体材料,将中间导体环的材料设置为理想绝缘材料。
执行第六步S6:设置二维轴对称模型中的源、边界及其他材料属性,进行带集磁器电磁成形的二维轴对称计算,得到待加工工件的电磁力分布。具体是设置二维轴对称模型中驱动线圈的脉冲电流为实验测量得到的驱动线圈电流,在模型的外边缘设置磁绝缘边界条件,划分有限元网格,利用商用有限元软件计算得到待加工工件的电磁力分布。
按照这一方法,即可将不具有对称性的带集磁器电磁成形模型简化为二维轴对称模型,可提高模型计算速度。
Claims (3)
1.一种适用于带集磁器管件电磁成形的二维轴对称仿真模型建模方法,其特征在于包括以下步骤:
步骤1:用与集磁器(2)相同的导体材料,替换集磁器的微小间隙(21),使其变成一整块具有轴对称型的导体环(22);
步骤2:将所述导体环分割为三个导体环,分别为:靠近驱动线圈(1)一侧的导体环(23)、靠近待加工工件一侧的导体环(25)、中间导体环(24);
步骤3:建立二维轴对称模型,包括具有轴对称型的驱动线圈、待加工工件(3)、三个导体环,以及相应的空气域(4);
步骤4:将靠近驱动线圈(1)一侧的导体环(23)内电流密度面积分记为I1,将靠近待加工工件一侧的导体环(25)内电流密度面积分记为I2,设置约束条件I1=-I2;
步骤5:将靠近驱动线圈(1)一侧的导体环(23)、靠近待加工工件一侧的导体环(25)的材料设置为集磁器导体材料,将中间导体环(24)的材料设置为理想绝缘材料;
步骤6:设置二维轴对称模型中的源、边界及其他材料属性,进行带集磁器电磁成形的二维轴对称计算,得到待加工工件(3)的电磁力分布;其中:
二维轴对称模型中的源是指驱动线圈中的脉冲电流;
二维轴对称模型中的边界指在模型外边界,磁通密度或矢量磁位满足的条件;
二维轴对称模型中的其他材料属性指的是空气域的材料属性,包括空气域的电阻率和磁导率。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910465106.1A CN110210116B (zh) | 2019-05-30 | 2019-05-30 | 一种适用于带集磁器管件电磁成形的二维轴对称仿真模型及建模方法 |
CN202310124543.3A CN116257994A (zh) | 2019-05-30 | 2019-05-30 | 一种适用于带集磁器管件电磁成形的二维轴对称仿真模型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910465106.1A CN110210116B (zh) | 2019-05-30 | 2019-05-30 | 一种适用于带集磁器管件电磁成形的二维轴对称仿真模型及建模方法 |
Related Child Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310124543.3A Division CN116257994A (zh) | 2019-05-30 | 2019-05-30 | 一种适用于带集磁器管件电磁成形的二维轴对称仿真模型 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110210116A CN110210116A (zh) | 2019-09-06 |
CN110210116B true CN110210116B (zh) | 2023-04-07 |
Family
ID=67789781
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310124543.3A Pending CN116257994A (zh) | 2019-05-30 | 2019-05-30 | 一种适用于带集磁器管件电磁成形的二维轴对称仿真模型 |
CN201910465106.1A Active CN110210116B (zh) | 2019-05-30 | 2019-05-30 | 一种适用于带集磁器管件电磁成形的二维轴对称仿真模型及建模方法 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310124543.3A Pending CN116257994A (zh) | 2019-05-30 | 2019-05-30 | 一种适用于带集磁器管件电磁成形的二维轴对称仿真模型 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN116257994A (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112906254A (zh) * | 2019-11-19 | 2021-06-04 | 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 | 一种变压器绕组漏磁场仿真建模方法 |
CN112387843B (zh) * | 2020-10-27 | 2021-11-16 | 三峡大学 | 圆环电磁胀形分析、控制方法及胀形装置 |
CN114871326B (zh) * | 2022-04-02 | 2024-05-28 | 三峡大学 | 采用金属涂层集磁器的管件胀形方法及装置 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AUPR179500A0 (en) * | 2000-11-30 | 2000-12-21 | Saintech Pty Limited | Ion source |
JP4744339B2 (ja) * | 2006-04-04 | 2011-08-10 | 株式会社神戸製鋼所 | 電磁成形コイル |
CN106660101B (zh) * | 2014-05-04 | 2019-05-14 | 贝瓦克生产机械有限公司 | 容器的电磁成形系统和方法 |
CN105127284B (zh) * | 2015-09-29 | 2017-04-26 | 华中科技大学 | 一种分层控制的电磁渐进成形方法 |
CN106529100B (zh) * | 2016-12-23 | 2019-09-27 | 大工(青岛)新能源材料技术研究院有限公司 | 一种基于有限元分析螺管超导磁体线圈磁场仿真方法 |
CN208195342U (zh) * | 2018-04-02 | 2018-12-07 | 三峡大学 | 一种铝合金管件交错变形电磁连接装置 |
-
2019
- 2019-05-30 CN CN202310124543.3A patent/CN116257994A/zh active Pending
- 2019-05-30 CN CN201910465106.1A patent/CN110210116B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN116257994A (zh) | 2023-06-13 |
CN110210116A (zh) | 2019-09-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110210116B (zh) | 一种适用于带集磁器管件电磁成形的二维轴对称仿真模型及建模方法 | |
CN106513679B (zh) | 一种适用于大型金属零件的激光增材制造方法 | |
Pantokratoras et al. | EMHD free-convection boundary-layer flow from a Riga-plate | |
CN104174842B (zh) | 一种基于交变磁场的金属丝材增材设备及增材方法 | |
CN112191849A (zh) | 基于温度分布的梯度多孔散热装置设计及增材制造方法 | |
Zhu et al. | Electrochemical machining on blisk cascade passage with dynamic additional electrolyte flow | |
CN105138772B (zh) | 变截面构件电子束焊的有限元模拟方法 | |
CN104268317B (zh) | 机械零件圆角结构形状优化方法 | |
CN111187884B (zh) | 一种用于焊管中频感应加热热源优化方法 | |
CN104132086A (zh) | 平面不对称电磁阻尼器 | |
Zhang et al. | Study on debris evacuation of EDM small hole processing on titanium alloy | |
CN108363858A (zh) | 基于极坐标转换的圆筒直线振荡电机磁场确定方法 | |
CN104361227A (zh) | 一种定子铁心背部支持筋感应电压计算方法 | |
Ai et al. | Cathodic modification, numerical simulation and experimental investigation on electrochemical machining for the small inner-walled ring groove | |
CN105871257B (zh) | 一种动圈式磁悬浮永磁平面电机的转矩电流分配方法 | |
CN108509681A (zh) | 一种基于光伏隔离变压器的涡流分析方法 | |
CN113255191A (zh) | 一种感应加热模型的等效电路参数辨识方法 | |
CN204639401U (zh) | 转炉炉壳水冷管安装用工装 | |
CN103737164A (zh) | 等离子弧柱压缩器 | |
CN218387979U (zh) | 一种基于金属增材制造用铜感应加热器 | |
Wu et al. | On the fluid behavior and stability of Ti-6Al-4V titanium alloy GMAW molten pool: Effect of the longitudinal magnetic field | |
CN204004126U (zh) | 平面不对称电磁阻尼器 | |
CN112016228A (zh) | 一种基于水冷压缩系数的水下焊接热源模型的建模方法 | |
CN117272762B (zh) | 水冷磁体线圈对流换热系数确定方法及系统 | |
CN104091023A (zh) | 一种在电磁场棱边单元法中施加规范的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |