CN108694296A - 一种针对大型复杂模型电场数值计算的有限元建模方法 - Google Patents

一种针对大型复杂模型电场数值计算的有限元建模方法 Download PDF

Info

Publication number
CN108694296A
CN108694296A CN201810649276.0A CN201810649276A CN108694296A CN 108694296 A CN108694296 A CN 108694296A CN 201810649276 A CN201810649276 A CN 201810649276A CN 108694296 A CN108694296 A CN 108694296A
Authority
CN
China
Prior art keywords
model
entity
finite element
layer
outsourcing
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.)
Granted
Application number
CN201810649276.0A
Other languages
English (en)
Other versions
CN108694296B (zh
Inventor
杜志叶
田雨
金硕
谢鸣
谢一鸣
金颀
王建
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN201810649276.0A priority Critical patent/CN108694296B/zh
Publication of CN108694296A publication Critical patent/CN108694296A/zh
Application granted granted Critical
Publication of CN108694296B publication Critical patent/CN108694296B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

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)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种针对大型复杂模型电场数值计算的有限元建模方法。当应用有限元法求解大型复杂模型的表面电场时,必须针对全模型展开建模和计算。同时为了保证计算精度,往往需要对关注设备表面网格进行精细控制。但是当求解域内实体较多且结构复杂时,实现整体有限元模型的精细建模较为困难,且易造成大量网格堆积,难以保证计算效率。对此,公开了一种针对大型复杂模型电场数值计算的有限元建模方法。该方法在各关注设备表面生成多层相同几何结构外包实体,完成全模型区域划分后,基于各层外包实体参数化实现模型表面多层网格加密。在源头上解决整体模型中关注设备表面网格难以精细控制的问题,提高了大型复杂有限元模型的建模效率及质量。

Description

一种针对大型复杂模型电场数值计算的有限元建模方法
技术领域
本发明涉及一种有限元建模技术,尤其涉及一种专门针对大型复杂模型电场数值计算的有限元高效、精确建模方法。
背景技术
有限元法是电场计算中应用较为广泛的一种数值方法。当应用有限元法求解诸如特高压直流阀厅等大型复杂模型的表面电场时,必须针对全模型展开建模和计算。同时为了保证有限元法的计算精度,往往需要对关注设备表面网格进行精细控制。但是当求解域内实体较多且结构复杂时,实现整体有限元模型的精细建模较为困难,且易造成大量网格堆积,难以保证计算效率。
目前,针对大型复杂模型应用较为普遍的有限元建模方法主要是子模型法。子模型法将子模型与整体较粗糙的模型分割开,单独进行网格控制,提高了各关注区域电场的计算精度。但子模型法需要额外进行子模型的剖分及求解,计算过程较为繁杂;同时子模型法要求子模型切割边界远离载荷集中区域,其边界插值存在一定的误差传递。
实际工程问题的有限元电场求解过程中,大型复杂模型一般需采用专业建模软件进行外部建模。采用专业建模软件Solidworks与有限元仿真软件ANSYS进行联合仿真,在模型建立时既在各重点关注设备表面批量生成多层相同几何形态的外包实体,后续剖分时基于各层外包实体参数化实现模型表面多层网格加密,从而在源头上解决整体模型中关注设备表面网格难以精细控制的问题,提高大型复杂有限元模型的建模效率及质量。
发明内容
本发明的目的是提供一种针对大型复杂模型电场数值计算的有限元建模方法。采用本发明的方法可实现大型复杂有限元模型的高效、精确建模,有效实现关注设备表面网格的精细控制。
本发明提供的方法主要由关注设备表面多层外包实体生成、全模型区域划分、模型表面网格加密等几大部分组成。
本发明所提供的一种针对大型复杂模型电场数值计算的有限元建模方法,其特征在于,
关注设备表面多层外包实体生成步骤:从原模型基本几何形态出发,基于Solidworks软件在关注设备表面沿法向向外建立多层厚度可控的外包曲面,与ANSYS软件搭接后,使用参数化命令将外包曲面批量转化为实体。
全模型区域划分步骤:求解域存在多个设备时,所生成外包实体间及其与原实体模型间易发生交叠。对外包实体与原实体模型间位置关系进行分析,在保证原实体模型完整的前提下,借助布尔操作完成全模型中各区域的划分。并以模型各实体体积及形心坐标作为特征量,进行模型组件定义,完成实体属性归类。
模型表面网格加密步骤:结合有限元电场计算网格划分原则,基于所建立多层外包实体在模型表面逐层生成加密网格。由于导体设备表面电场主要沿其法向分布,可控制各层网格沿设备表面切向尺寸大于法向尺寸,同时网格尺寸由内向外逐层增大,保证网格具有良好的划分过渡性。
在上述的一种针对大型复杂模型电场数值计算的有限元建模方法,所述关注设备表面多层外包实体生成步骤的具体方法包括:
步骤1、在Solidworks软件中对关注设备各个表面进行“等距曲面”操作,在设备表面形成一层面包围。控制等距参数,初步生成厚度由内向外逐层增大的多层外包曲面。
步骤2、将模型保存为.x_t格式,导入ANSYS软件中实现两个软件的搭接,在ANSYS软件中设置模型几何类型为所有实体,保证外包曲面成功导入。
步骤3、由于每层外包曲面都是各自独立的面包围,基于APDL参数化命令,将各层外包曲面分别选取,使用“VA”命令,将各层外包曲面所包围区域实体化,最终在结构上形成多层外包实体。
在上述的一种针对大型复杂模型电场数值计算的有限元建模方法,所述整体模型区域划分步骤的具体方法包括:
步骤1、在外包曲面实体化步骤之前,采用参数化命令“*GET”得到原模型中各实体的体积及形心位置坐标信息,将其存储于矩阵M中。
步骤2、对整体模型进行布尔操作,首先使用“VSBV”命令将外包实体中与原模型交叠的区域删除,确保原模型完整性,然后对整体模型进行“VOVLAP”操作,实现全模型的搭接。
步骤2、调取矩阵M中所储存各实体的体积及形心位置坐标信息,以其为特征量,在整体模型中快速对比选取原模型各实体,进而完成原模型及外包实体的组件定义。对各区域属性归类后,最终实现整体模型的区域划分。
在上述的一种针对大型复杂模型电场数值计算的有限元建模方法,所述模型表面网格加密步骤的具体方法包括:
步骤1、外包实体层次有序地分布于关注设备表面,通过选取与所定义原模型相接触的实体,由内向外逐层定义各层外包实体。
步骤2、每层外包实体单独剖分为一层加密网格,采用自由剖分方式,控制关注设备表面网格沿设备表面切向尺寸大一些,沿法向尺寸小一些,形成由内向外厚度逐层增加的多层加密网格,最后完成全场域的剖分。
因此,本发明具有如下优点:(1)本发明对全模型展开建模和计算,确保了模型的完整性,保证了大型复杂模型有限元电场数值计算的准确性。(2)本发明实现关注设备表面网格的精细控制的同时,有效减小了整体模型网格划分的单元数量,有效提高了计算精度与计算效率。(3)本发明利用Solidworks软件与ANSYS软件进行联合仿真,采用APDL命令可参数化实现有限元建模过程,操作便捷,实用性强。
附图说明
图1是本方法实现的流程示意图。
图2是±800kV特高压直流换流站高端阀厅模型。
图3是阀塔中性管母避雷器金具外包曲面。
图4是阀塔中性管母避雷器金具表面网格加密。
图5是±800kV特高压直流换流站高端阀厅表面电场分布。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实时方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
以某±800kV特高压直流换流站高端阀厅为例,阐述其金具表面电场计算时的有限元建模过程。如图2所示为所搭建阀厅整体模型,为了降低模型剖分难度并减少不必要的计算,对套管、阀塔等部位进行合理简化。
由于阀厅内屏蔽金具种类多样,不同金具结构大小差异较大,为了对网格进行合理控制,将避雷器均压环、OCT均压环等内径较小金具与出线套管均压球等尺寸较大的屏蔽金具进行分类,分别生成不同厚度外包实体进行网格控制。在各关注设备表面生成3层外包实体,第i层外包实体厚度为Δhi,具体参数如下表1所示。
表1阀厅金具外包实体参数
根据上述外包实体参数,基于Solidworks软件中“等距曲面”操作在阀厅模型中重点关注金具表面生成相应外包曲面。如图3所示为阀厅模型中阀塔中性管母避雷器金具所生成外包曲面。外包曲面与原实体模型基本几何结构相同,尺寸可控。将模型导入有限元分析软件ANSYS中,由于每个金具的每层外包曲面都是一个独立的面包围,使用参数化命令“VA”即可将每层外包曲面所包围区域转化为实体,程序化实现所有外包实体的生成。
在对阀厅全模型进行区域划分之前,首先采用参数化命令记录阀厅原模型各实体的体积及形心坐标信息。进行布尔操作时,使用“VSBV”命令将外包实体中与原模型交叠的区域删除,确保原模型完整性,然后对整体模型进行“VOVLAP”操作,实现整体模型的搭接。布尔操作完成后,利用所记录原模型各实体的体积及形心坐标信息作为特征量快速选取求解域中阀厅原模型各实体,完成原模型及外包实体的组件定义,完成各区域属性归类。
对模型进行网格划分时,采用自由剖分方式,通过其表面多层外包实体进行网格控制,每层外包实体单独进行剖分控制生成一层加密网格。由于阀厅内各金具自身表面电位恒等,故其表面电场沿其切向上分量基本为0,主要分布在其表面法向上。为保证网格具有良好的划分过渡性,同时有效控制全场域网格数量,外包实体生成时已控制其厚度由内向外逐层增大,网格划分时可控制关注设备表面网格沿设备表面切向尺寸大一些,沿法向尺寸小一些,最终形成由内向外厚度逐层增加的多层加密网格,然后再完成全场域的剖分。阀塔中性管母避雷器表面网格加密如图4所示。由于外包实体层次有序地分布于模型表面,通过所定义原模型可将与其相接触的第一层外包实体快速选中,进而逐层选取各层外包实体,使用参数化命令可批量实现其剖分控制,操作便捷,效率高。
按照上述步骤完成阀厅全模型剖分,网格数量得到有效控制,总数约1800万。使用瞬时电位加载法,在额定功率运行状态下,计算得到阀厅设备表面电场分布如图5所示。可以看出,模型表面电场过渡平滑,电场强度最大值出现在D侧换流变C相出线下部套管圆柱形均压球表面,最大值为13.1kV/cm,其余部位的电场强度相对较小。证明了本发明的可行性及有效性。
上述实施例所述是用以具体说明本专利,文中虽通过特定的术语进行说明,但不能以此限定本专利的保护范围,熟悉此技术领域的人士可在了解本专利的精神与原则后对其进行变更或修改而达到等效目的,而此等效变更和修改,皆应涵盖于权利要求范围所界定范畴内。

Claims (4)

1.一种针对大型复杂模型电场数值计算的有限元建模方法,其特征在于,
关注设备表面多层外包实体生成步骤:从原模型基本几何形态出发,基于Solidworks软件在关注设备表面沿法向向外建立多层厚度可控的外包曲面,与ANSYS软件搭接后,使用参数化命令将外包曲面批量转化为实体;
全模型区域划分步骤:求解域存在多个设备时,所生成外包实体间及其与原实体模型间易发生交叠;对外包实体与原实体模型间位置关系进行分析,在保证原实体模型完整的前提下,借助布尔操作完成全模型中各区域的划分;并以模型各实体体积及形心坐标作为特征量,进行模型组件定义,完成实体属性归类;
模型表面网格加密步骤:结合有限元电场计算网格划分原则,基于所建立多层外包实体在模型表面逐层生成加密网格,控制各层网格沿模型表面切向尺寸大于法向尺寸,同时网格尺寸由内向外逐层增大。
2.根据权利要求1所述的一种针对大型复杂模型电场数值计算的有限元建模方法,其特征在于,所述关注设备表面多层外包实体生成步骤的具体方法包括:
步骤1、在Solidworks软件中对关注设备各个表面进行“等距曲面”操作,在设备表面形成一层面包围;控制等距参数,初步生成厚度由内向外逐层增大的多层外包曲面;
步骤2、将模型保存为.x_t格式,导入ANSYS软件中实现两个软件的搭接,在ANSYS软件中设置模型几何类型为所有实体,保证外包曲面成功导入;
步骤3、由于每层外包曲面都是各自独立的面包围,基于APDL参数化命令,将各层外包曲面分别选取,使用“VA”命令,将各层外包曲面所包围区域实体化,最终在结构上形成多层外包实体。
3.根据权利要求1所述的一种针对大型复杂模型电场数值计算的有限元建模方法,其特征在于,所述整体模型区域划分步骤的具体方法包括:
步骤1、在外包曲面实体化步骤之前,采用参数化命令“*GET”得到原模型中各实体的体积及形心位置坐标信息,将其存储于矩阵M中;
步骤2、对整体模型进行布尔操作,首先使用“VSBV”命令将外包实体中与原模型交叠的区域删除,确保原模型完整性,然后对整体模型进行“VOVLAP”操作,实现全模型的搭接;
步骤2、调取矩阵M中所储存各实体的体积及形心位置坐标信息,以其为特征量,在整体模型中快速对比选取原模型各实体,进而完成原模型及外包实体的组件定义;对各区域属性归类后,最终实现整体模型的区域划分。
4.根据权利要求1所述的一种针对大型复杂模型电场数值计算的有限元建模方法,其特征在于,所述模型表面网格加密步骤的具体方法包括:
步骤1、外包实体层次有序地分布于关注设备表面,通过选取与所定义原模型相接触的实体,由内向外逐层定义各层外包实体;
步骤2、每层外包实体单独剖分为一层加密网格,采用自由剖分方式,控制关注设备表面网格沿设备表面切向尺寸大一些,沿法向尺寸小一些,形成由内向外厚度逐层增加的多层加密网格,最后完成全场域的剖分。
CN201810649276.0A 2018-06-22 2018-06-22 一种针对大型复杂模型电场数值计算的有限元建模方法 Active CN108694296B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810649276.0A CN108694296B (zh) 2018-06-22 2018-06-22 一种针对大型复杂模型电场数值计算的有限元建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810649276.0A CN108694296B (zh) 2018-06-22 2018-06-22 一种针对大型复杂模型电场数值计算的有限元建模方法

Publications (2)

Publication Number Publication Date
CN108694296A true CN108694296A (zh) 2018-10-23
CN108694296B CN108694296B (zh) 2021-09-24

Family

ID=63849898

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810649276.0A Active CN108694296B (zh) 2018-06-22 2018-06-22 一种针对大型复杂模型电场数值计算的有限元建模方法

Country Status (1)

Country Link
CN (1) CN108694296B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110472318A (zh) * 2019-07-31 2019-11-19 中国南方电网有限责任公司超高压输电公司 一种基于apdl语言的换流站阀厅金具数值仿真批量化建模仿真方法
CN110750916A (zh) * 2019-10-31 2020-02-04 南方电网科学研究院有限责任公司 用于柔性直流换流站阀厅电场计算的电位加载方法及装置
CN111125949A (zh) * 2019-12-06 2020-05-08 北京科技大学 一种有限元分析的大规模并行网格划分系统及方法
CN111191390A (zh) * 2018-10-26 2020-05-22 中国航发商用航空发动机有限责任公司 表面带有内凹部的零件建模方法、设备以及电子设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020120430A1 (en) * 2001-02-23 2002-08-29 Korea Advanced Institute Of Science And Technology Hexahedral finite element modeling method for controlling element size and storage medium therefor
CN101615214A (zh) * 2009-07-23 2009-12-30 清华大学 一种对磁共振成像rf线圈进行性能分析的方法
CN102201669A (zh) * 2010-04-30 2011-09-28 北京网联直流工程技术有限公司 一种±660kV阀厅直流金具选型方法
CN106294914A (zh) * 2015-06-04 2017-01-04 长沙理工大学 一种考虑非线性各向异性的换流变压器阀侧三维电场仿真方法
CN107944130A (zh) * 2017-11-22 2018-04-20 国网江西省电力有限公司电力科学研究院 一种基于有限元电场仿真的金具起晕电压确定方法
CN108133089A (zh) * 2017-12-11 2018-06-08 三峡大学 一种输电线路屏蔽环静电场有限元数值计算的荷载施加方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020120430A1 (en) * 2001-02-23 2002-08-29 Korea Advanced Institute Of Science And Technology Hexahedral finite element modeling method for controlling element size and storage medium therefor
CN101615214A (zh) * 2009-07-23 2009-12-30 清华大学 一种对磁共振成像rf线圈进行性能分析的方法
CN102201669A (zh) * 2010-04-30 2011-09-28 北京网联直流工程技术有限公司 一种±660kV阀厅直流金具选型方法
CN106294914A (zh) * 2015-06-04 2017-01-04 长沙理工大学 一种考虑非线性各向异性的换流变压器阀侧三维电场仿真方法
CN107944130A (zh) * 2017-11-22 2018-04-20 国网江西省电力有限公司电力科学研究院 一种基于有限元电场仿真的金具起晕电压确定方法
CN108133089A (zh) * 2017-12-11 2018-06-08 三峡大学 一种输电线路屏蔽环静电场有限元数值计算的荷载施加方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
金硕 等: "特高压直流输电系统阀厅金具表面电场精细求解", 《高电压技术》 *
阮江军 等: "±800kV特高压直流换流站阀厅金具表面电场计算", 《高电压技术》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111191390A (zh) * 2018-10-26 2020-05-22 中国航发商用航空发动机有限责任公司 表面带有内凹部的零件建模方法、设备以及电子设备
CN111191390B (zh) * 2018-10-26 2023-09-01 中国航发商用航空发动机有限责任公司 表面带有内凹部的零件建模方法、设备以及电子设备
CN110472318A (zh) * 2019-07-31 2019-11-19 中国南方电网有限责任公司超高压输电公司 一种基于apdl语言的换流站阀厅金具数值仿真批量化建模仿真方法
CN110472318B (zh) * 2019-07-31 2022-11-22 中国南方电网有限责任公司超高压输电公司 一种基于apdl语言的换流站阀厅金具数值仿真批量化建模仿真方法
CN110750916A (zh) * 2019-10-31 2020-02-04 南方电网科学研究院有限责任公司 用于柔性直流换流站阀厅电场计算的电位加载方法及装置
CN111125949A (zh) * 2019-12-06 2020-05-08 北京科技大学 一种有限元分析的大规模并行网格划分系统及方法

Also Published As

Publication number Publication date
CN108694296B (zh) 2021-09-24

Similar Documents

Publication Publication Date Title
CN108694296A (zh) 一种针对大型复杂模型电场数值计算的有限元建模方法
Kataoka et al. Voltage stability limit of electric power systems with generator reactive power constraints considered
CN103093007B (zh) 输电铁塔三维虚拟装配方法
CN106600700B (zh) 一种三维模型数据处理系统
Zheng et al. Interpolation over arbitrary topology meshes using a two-phase subdivision scheme
CN106970999B (zh) 一种bpa数据库建模与维护方法
CN102270261A (zh) 基于cad图模一体化的智能变电站模型设计方法
CN104866613A (zh) 一种电网设备设施三维组件化模型数据库的建立方法
CN102118026B (zh) 基于etap的大规模电网的谐波负序潮流方法
CN109284536A (zh) 基于bim的配电网模型族库搭建方法
CN111724479A (zh) 一种变电站电气一次设备三维建模方法和系统
Fu et al. Generalized Debye sources-based EFIE solver on subdivision surfaces
CN105808881A (zh) 高压电缆中间接头的有限元优化设计方法
CN113742534B (zh) 一种多能源系统控制和保护系统的自动建模及可视化方法
Wu et al. Hybridizing FD-TD analysis with unconditionally stable FEM for objects of curved boundary
Kapur et al. Modeling of integrated RF passive devices
WO2023206281A1 (zh) 一种模拟电路的快速仿真优化方法及其系统
WO2023226113A1 (zh) 直流电弧仿真方法、装置、计算机设备和存储介质
Shafieipour et al. On new triangle quadrature rules for the locally corrected Nyström method formulated on NURBS-generated Bézier surfaces in 3-D
Du et al. An efficient modeling method for large-scale complex finite-element electric field numerical calculation
CN113609556B (zh) 一种基础c#的gim格式与xdb格式无损转换方法
Guarin et al. Construction and simulation of a planar transformer prototype
Wang et al. Gauss-Seidel Progressive Iterative Approximation (GS-PIA) for Loop Surface Interpolation.
Liu et al. Optimization of electrode shape using the boundary element method
Slavkova et al. Computer-aided design of current transformer for power electronics converters with voltage feedback

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