CN108052697B - 一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法 - Google Patents
一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法 Download PDFInfo
- Publication number
- CN108052697B CN108052697B CN201711149404.7A CN201711149404A CN108052697B CN 108052697 B CN108052697 B CN 108052697B CN 201711149404 A CN201711149404 A CN 201711149404A CN 108052697 B CN108052697 B CN 108052697B
- Authority
- CN
- China
- Prior art keywords
- conductor
- shell
- temperature
- calculation
- workbench
- 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.)
- Expired - Fee Related
Links
- 238000004364 calculation method Methods 0.000 title claims description 92
- 238000000034 method Methods 0.000 claims abstract description 37
- 230000008878 coupling Effects 0.000 claims abstract description 13
- 238000010168 coupling process Methods 0.000 claims abstract description 13
- 238000005859 coupling reaction Methods 0.000 claims abstract description 13
- 230000008602 contraction Effects 0.000 claims abstract description 12
- 239000004020 conductor Substances 0.000 claims description 91
- 239000003570 air Substances 0.000 claims description 57
- 230000005855 radiation Effects 0.000 claims description 50
- 239000007787 solid Substances 0.000 claims description 37
- 238000002076 thermal analysis method Methods 0.000 claims description 32
- 239000012530 fluid Substances 0.000 claims description 22
- 239000000463 material Substances 0.000 claims description 22
- 238000004458 analytical method Methods 0.000 claims description 18
- 230000005540 biological transmission Effects 0.000 claims description 18
- 230000002500 effect on skin Effects 0.000 claims description 16
- 230000001052 transient effect Effects 0.000 claims description 16
- 238000010438 heat treatment Methods 0.000 claims description 14
- 230000003068 static effect Effects 0.000 claims description 8
- 238000012916 structural analysis Methods 0.000 claims description 8
- 230000008646 thermal stress Effects 0.000 claims description 8
- 239000012080 ambient air Substances 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 238000012544 monitoring process Methods 0.000 claims description 6
- 230000020169 heat generation Effects 0.000 claims description 5
- 230000004323 axial length Effects 0.000 claims description 4
- 239000011343 solid material Substances 0.000 claims description 4
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 230000017525 heat dissipation Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 4
- 238000011156 evaluation Methods 0.000 abstract 1
- 238000012216 screening Methods 0.000 abstract 1
- 238000002474 experimental method Methods 0.000 description 12
- 238000004088 simulation Methods 0.000 description 12
- 238000012546 transfer Methods 0.000 description 9
- 238000013461 design Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000006698 induction Effects 0.000 description 3
- 238000001816 cooling Methods 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 238000012805 post-processing Methods 0.000 description 2
- 230000035882 stress Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 239000012811 non-conductive material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000009423 ventilation 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
- 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
- G06F2113/00—Details relating to the application field
- G06F2113/16—Cables, cable trees or wire harnesses
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
一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法，根据廊管GIL布置的结构特性，在workbench中建立GIL整体环境下的电气热多物理场耦合计算模型，根据相应运行工况，计算获得GIL三维温度场，再利用有限元方法在workbench中获得相应的GIL胀缩形变数据。该计算方法可为廊管GIL元件选择及筛查提供基础参考，并基于运行基本参量，计算出GIL运行过程中的温度及胀缩形变的理论数据，为系统运行状态评估提供技术支撑；本发明能较准确的计算气体绝缘输电线路的温度场，并且可以通过改变施加的载荷及边界条件等可以对多种不同运行条件下的气体绝缘输电线路温度场进行仿真计算，同时可以计算GIL管道的热变形，对气体绝缘输电线路的结构设计及运行条件的选择有参考意义。
Description
技术领域
本发明涉及特高压气体绝缘输电技术领域，特别涉及一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法。
背景技术
气体绝缘输电线路(GIL)具有电压等级高、输送容量大、电能损耗低、电磁干扰小等优点，且在传统架空线路受约束的区域，如需跨越海峡、大江大海等条件下具有巨大优势，因此GIL将得到广泛应用。
在运行过程中，由于焦耳热损耗，输电线路一直处于发热状态，导致线路温度升高。温升不仅降低了线路的最大载流能力，而且温度超过一定限度时将影响内部非导电材料的绝缘性能，降低线路使用寿命，甚至发生严重的放电事故；同时温度的变化会使输电管道发生热变形，由热变形产生的应力及疲劳损伤对管路寿命及运行情况有很大影响。因此在设计阶段对线路的运行温度进行合理的预测，并计算管路热应力及热变形对输电线路的结构设计及运行条件的控制有着重要意义。
目前，GIL温度计算方法主要有两种：解析法和有限元法。解析法通过建立焦耳热损耗与对流换热、辐射换热的热平衡关系迭代求解导体与外壳的温度，该方法计算速度快，可以根据运行条件很快的计算出输电线路导体及外壳的温度，但不能反映GIL整体温度的分布。有限元法主要有两种：一种是吴晓文等人建立的包含外部空气域在内的二维多物理场耦合有限元模型，避免了在固体表面施加对流换热边界条件；一种是王健等人建立的不包含外部空气域的三维多物理场耦合模型，该模型可以同时得到输电线路径向与轴向的温度分布，但是需要在固体表面施加对流及辐射换热边界条件。而且以上两种有限元方法都无法考虑在有风速的情况下的GIL管道温度分布。
在计算出管路温度场的情况下对管路的热变形进行计算，可以对管路的设计，如伸缩节伸缩量的选择，连接、固定装置的设计等有一定的指导意义及参考价值。
发明内容
本发明的目的在于克服现有方法的不足，提供一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法，建立了三维仿真模型，通过设置空气域出入口边界条件，来计算在考虑廊管内由风机通风的情况下GIL管道的温度分布，通过设置气体域的辐射项，并对气体边界设置合理的辐射系数，辐射换热自动迭代求解；简化了边界条件的设置过程，仿真结果与实验结果的误差在5％以内，并且在计算出管路温度场的温度后，可以计算管路热变形，对管路的设计有参考价值。
为了达到上述目的，本发明的目的通过以下技术方案实现：
一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法，包括以下步骤：
A：根据结构特性，对模型进行简化，使用SolidWorks软件建立包含绝缘气体及环境空气的计算模型；
B：把物理模型的几何文件导入网格划分软件ANSYS ICEM CFD中，进行流体和固体模块的网格划分；
C：在workbench中建立CFX模块，将网格导入CFX中，建立流体、固体域，并对各个计算域进行设置；
D：在固体计算域上施加考虑集肤效应的焦耳热载荷；
E：对流、固计算域进行耦合，并施加边界条件，然后在CFX模块中进行温度场计算；
F：将固体域网格文件导入workbench稳态或瞬态热分析模块，并将CFX的计算结果导入到热分析模块中的模型的对应边界上，再次对导体及外壳进行温度场计算；
G：将热分析模块中的计算结果导入workbench的静态或瞬态结构分析模块，计算热变形。
所述步骤A包括以下步骤：
A1：对模型进行简化，由于每一段管道两端的触头处的接触良好，导电性能优良，所以可以忽略触头处的热损耗对管道温度的影响，所以模型只截取部分轴向长度的管道来计算温度场的分布；
A2：使用SolidWorks进行建模，建立的模型包括空气域、圆筒状导体、圆筒状外壳、导体内部的圆柱形绝缘气体区域即内层绝缘气体以及导体和外壳之间的圆筒形绝缘气体区域即外层绝缘气体。
所述步骤B网格划分具体为：
将建立好的各个区域模型分别导入ICEM CFD软件中，采用分块划分网格的方式得到O型及H型结构化网格，并对气体域定义边界层，对边界层网格进行加密。
所述步骤C包括以下步骤：
C1：在workbench工作界面中放入CFX模块，将划分好的网格导入CFX中；
C2：创建计算中所需的新材料，固体材料的热物性视为定值，通过直接赋值的方式直接输入相应的导热率、比热容参数；气体为不可压缩流体，定压比热容视为定值，密度、动力粘度、导热系数考虑温度的影响，通过CEL语言定义气体属性：
密度计算公式为：
导热系数计算公式为：
动力粘度计算公式为：
其中：ρ(T)为气体密度，单位为kg/m3；ρ0为0℃时气体的密度，单位为kg/m3；λ(T)为导热系数，λ0为0℃时气体的导热系数；η(T)为气体动力粘度，单位为Pa·s；η0为0℃时气体的动力粘度；T0为参考温度，273K；S为Sutherland常数，单位为K；；T为该网格温度，单位为K，在迭代时由求解器实时给出；
当考虑空气湿度的影响时，可以采用相同方法建立湿空气的密度、导热系数、动力粘度公式，然后赋予新材料；
C3：创建流、固计算域，并为各个计算域指定相应材料；
C4：对于流体计算域，需要通过添加浮力项考虑气体由于受热产生的密度变化发生的自然对流，其中浮力项参考密度计算公式为：
其中ρc为浮力项参考密度，单位为kg/m3；Tc为环境温度，单位为K；
同时通过打开气体域辐射项，设置辐射面的辐射率考虑结构的辐射散热，只考虑气体界面上的辐射，不考虑气体内部对辐射的吸收。
所述步骤D包括以下步骤：
D1：考虑集肤效应的影响，导体和外壳的电阻计算公式为：
其中：Ri为导体或外壳电阻，单位为Ω/m；Kf为集肤效应系数；ρ20为20℃时导体或外壳的直流电阻率，单位为Ω·m；α20为电阻温度系数；Si为导体或外壳的横截面积，单位为mm2；Ti为导体或外壳的运行温度，单位为℃，在迭代时由求解器实时给出；
集肤效应系数计算公式为：
其中：C为导体或外壳壁厚，单位为mm；D为导体或外壳的外径，单位为mm；
D2：导体或外壳每米损耗功率为：
P＝I2Ri
其中：P为导体或外壳的每米损耗，即每米的焦耳热功率，单位为W/m；I为在导体中为GIL额定电流，单位为A，在外壳中为感应电流，当GIL长度小于20m时，感应电流为额定的0.95，当GIL长度大于20m时，感应电流等于GIL额定电流；
导体或外壳单位体积的焦耳热功率为：
其中：PV，为导体或外壳单位体积的焦耳热功率，单位为W/m3；
D3：通过CEL语言定义导体和外壳单位体积的焦耳热功率，并将其作为生热率载荷施加在导体及外壳上。
所述步骤E包括以下步骤：
E1：将流体域与固体域交界面进行耦合，包含热量交换；
E2：考虑空气由于受热发生自然对流时，空气域选用层流模型，空气两端面设置为开口边界条件；考虑风速时，根据廊管尺寸及空气流速综合计算空气流动的雷诺数，选择合适的流动模型，然后设置空气入口流速、温度，出口压强；空气外边界面设置为壁面，温度为环境温度，考虑环境吸收外壳的所有辐射，所以空气外边界面辐射系数设置为1；空气与输电管道交界面辐射系数根据输电管道外壳外表面的具体情况确定；
E3：绝缘气体由于只发生自然对流，所以选用层流模型；端面设置为对称面；内层绝缘气体不考虑辐射；外层绝缘气体边界的辐射系数根据输电管道外壳内表面及导体外表面的具体情况确定；
E4：设置求解器：包括迭代步数、时间步、收敛条件；设置温度监测点；然后进行求解，得到温度场的分布结果。
所述步骤F包括以下步骤：
F1：在workbench工作界面中建立有限元建模器Finite Element Modeler模块，并在该模块中输入由ICEM CFD生成的.uns格式的固体域网格文件；
F2：在workbench工作界面中建立稳态或瞬态热分析模块，将上述有限元建模器中的网格导入热分析模块中；
F3：在热分析模块中输入材料导热率；
F4：将CFX模块中的固体边界的温度计算结果导入到热分析模块中的模型的对应边界上，然后重新进行导体或外壳的温度场计算。
所述步骤G包括以下步骤：
所述步骤G包括以下步骤：
G1：在workbench工作界面建立静态或瞬态结构分析模块，将有限元建模器中的网格导入结构分析模块中；
G2：在结构分析模块中输入材料的密度、热膨胀系数、零热应力的参考温度；
G3：施加的约束条件为在导体或外壳的一端截面施加位移约束，仅约束轴向位移，并将热分析模块中计算的温度导入结构分析模块；同时在分析选项中打开弱弹簧选项，然后进行热变形计算。
本发明公开的一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法，建立了三维仿真模型。由于包含了气体域，对流换热自动迭代求解；通过设置气体域的辐射项，并对气体边界设置合理的辐射系数，辐射换热自动迭代求解；如此一来将简化边界条件的设置过程。同时由于是三维模型，可以通过设置空气域出入口边界条件，来计算在考虑廊管内由风机通风的情况下GIL管道的温度分布。采用此方法，对Minaguchi等人的实验进行仿真计算，仿真结果与实验结果的误差在5％以内。并且在计算出管路温度场的温度后，可以计算管路热变形，对管路的设计有参考价值。
附图说明
图1为本发明的流程图。
图2为实验模型图；
图3为有限元计算模型图；
图4为第三组实验条件下外壳温度场分布图：其中图4A为外壳径向温度分布图，图4B为外壳轴向温度分布图；
图5为第三组实验条件下外壳形变图；其中图5A为外壳总形变图，图5B为外壳轴向形变图。
具体实施方式
下面结合附图对本发明做详细叙述。
如图1所示，本发明一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法，依次包括以下步骤：
A：根据结构特性，对模型进行简化，使用SolidWorks软件建立包含绝缘气体及环境空气的计算模型，具体步骤为：
以廊道内三相两回特高压GIL输电线路为研究对象，对模型进行适当简化：(1)由于每一段管道两端的触头处的接触良好，导电性能优良，所以可以忽略触头处的热损耗对管道整体温度的影响，所以模型可以只截取部分轴向长度的管道来计算温度场的分布，所得温度场计算结果依然较准确。(2)忽略GIL管道的支架部分，只考虑廊道内的空气及GIL管道。
根据以上简化，利用SolidWorks软件建立的模型包括：空气域、圆筒状导体、圆筒状外壳、导体内部的圆柱形绝缘气体区域(内层绝缘气体)以及导体和外壳之间的圆筒形绝缘气体区域(外层绝缘气体)。
B：把物理模型的几何文件导入网格划分软件ANSYS ICEM CFD中，进行流体和固体模块的网格划分：
将建立好的各个区域模型分别导入ICEM CFD软件中，采用分块划分网格的方式得到O型及H型混合结构化网格，并对气体域定义边界层，对边界层网格进行加密。
C：在workbench中建立CFX模块，将网格导入CFX中，建立流体、固体域，并对各个计算域进行设置，具体步骤为：
C1：在workbench工作界面中放入CFX模块，将划分好的网格导入CFX中；
C2：创建计算中所需的新材料。固体材料的热物性视为定值，通过直接赋值的方式直接输入相应的密度、导热率、比热容等参数，常用的固体域材料及参数见下表。
气体为不可压缩流体，定压比热容视为定值，密度、动力粘度、导热系数考虑温度的影响，通过CEL语言定义气体属性：
密度计算公式为：
导热系数计算公式为：
动力粘度计算公式为：
其中：ρ(T)(kg/m3)为气体密度，ρ0(kg/m3)为0℃时气体的密度；λ(T)为导热系数，λ0为0℃时气体的导热系数；η(T)(Pa·s)为气体动力粘度，η0为0℃时气体的动力粘度；T0为参考温度，273K；S(K)为Sutherland常数；T(K)为该网格温度，在迭代时由求解器实时给出。
当考虑空气湿度的影响时，可以采用相同方法建立湿空气的密度、导热系数、动力粘度公式，然后赋予新材料。
C3：创建流、固计算域，并为各个计算域指定相应材料。
C4：对于流体计算域，需要通过添加浮力项考虑气体由于受热产生的密度变化发生的自然对流。其中浮力项参考密度计算公式为：
其中ρc(kg/m3)为浮力项参考密度；Tc(K)为环境温度。
同时通过打开气体域辐射项，设置辐射面的辐射率考虑结构的辐射散热，只考虑气体界面上的辐射，不考虑气体内部对辐射的吸收。
D：在固体计算域上施加考虑集肤效应的焦耳热载荷，具体步骤为：
D1：考虑集肤效应的影响，导体和外壳的电阻计算公式为：
其中：Ri(Ω/m)为导体或外壳电阻；Kf为集肤效应系数；ρ20(Ω·m)为20℃时导体或外壳的直流电阻率；α20为电阻温度系数；Si(mm2)为导体或外壳的横截面积；Ti(℃)为导体或外壳的运行温度，在迭代时由求解器实时给出。
集肤效应系数计算公式为：
其中：C(mm)为导体或外壳壁厚；D(mm)为导体或外壳的外径。
D2：导体或外壳每米损耗功率为：
P＝I2Ri
其中：P(W/m)，为导体或外壳的每米损耗，即每米的焦耳热功率；I(A)，在导体中为GIL额定电流，在外壳中为感应电流，当GIL长度小于20m时，感应电流为额定的0.95，当GIL长度大于20m时，感应电流等于GIL额定电流。
导体或外壳单位体积的焦耳热功率为：
其中：PV(W/m3)，为导体或外壳单位体积的焦耳热功率。
D3：通过CEL语言定义导体和外壳单位体积的焦耳热功率，并将其作为生热率载荷施加在导体及外壳上。
E：对流、固计算域进行耦合，并施加边界条件，然后在CFX模块中进行温度场计算，具体步骤为：
E1：将流体域与固体域交界面进行耦合，包含热量交换。
E2：考虑空气由于受热发生自然对流时，空气域选用层流模型，空气两端面设置为开口边界条件；考虑风速时，根据廊管尺寸及空气流速综合计算雷诺数，根据雷诺数判断空气流动状态，然后选择合适的流动计算模型，设置空气入口流速、温度，出口压强；空气外边界面设置为壁面，温度为环境温度，考虑环境吸收外壳的所有辐射，所以空气外边界面辐射系数设置为1；空气与输电管道交界面辐射系数根据输电管道外壳外表面的具体情况确定。
E3：绝缘气体由于只发生自然对流，所以选用层流模型；端面设置为对称面；内层绝缘气体不考虑辐射；外层绝缘气体边界的辐射系数根据输电管道外壳内表面及导体外表面的具体情况确定。
E4：设置求解器：包括迭代步数、时间步、收敛条件等；设置温度监测点；然后进行求解。当温度监测点的温度达到稳定时，即视为求解收敛，最后在后处理器中进行后处理，得到温度场的分布等其他结果。
F：将固体域网格文件导入workbench稳态或瞬态热分析模块，并将CFX的计算结果导入到热分析模块中的模型的对应边界上，再次对导体及外壳进行温度场计算，具体步骤为：
F1：在workbench工作界面中建立Finite Element Modeler(有限元建模器)模块，并在该模块中输入由ICEM CFD生成的.uns格式的固体域网格文件。
F2：在workbench工作界面中建立稳态或瞬态热分析模块，将上述有限元建模器中的网格导入热分析模块中。
F3：在热分析模块中输入材料导热率。
F4：将CFX模块中的固体边界的温度计算结果导入到热分析模块中的模型的对应边界上，然后重新进行导体或外壳的温度场计算。
G：将热分析模块中的计算结果导入workbench的静态或瞬态结构分析模块，计算热变形，具体步骤为：
G1：在workbench工作界面建立静态或瞬态结构分析模块，将有限元建模器中的网格导入结构分析模块中。
G2：在结构分析模块中输入材料的密度、热膨胀系数、零热应力的参考温度。
G3：将热分析模块中计算的温度导入结构分析模块；由于计算模型为部分管道模型，所以在管道端面不施加位移约束，只需在分析选项中打开弱弹簧选项即可；最后进行热应力、热变形计算。
实施例
按本发明提出的方法，对Minaguchi等人的实验进行仿真计算，实验模型如图2所示：
导体及外壳尺寸参数如下：
| 内径/mm | 外径/mm | |
| 导体 | 140 | 160 |
| 外壳 | 480 | 490 |
实验环境为：室内，无太阳辐射；环境空气为自然对流状态。根据实验结果得到，在此环境下，管道整体的温度场分布在轴向上是均匀的。实验在三组运行条件下进行：
| 第一组 | 第二组 | 第三组 | |
| 运行电流/(A) | 2600 | 4400 | 5600 |
| 环境温度/(℃) | 26 | 16 | 24 |
根据以上描述，对于此实验的模拟本方法依次包括以下步骤：
A：根据结构特性，对模型进行简化，使用SolidWorks软件建立包含绝缘气体及环境空气的计算模型，具体步骤为：
对模型进行适当简化：由于管道两端的触头处的接触良好，导电性能优良，所以可以忽略触头处的热损耗对管道整体温度的影响，模型可以只截取部分轴向长度的管道来计算温度场的分布。
因此，根据导体及外壳尺寸，利用SolidWorks软件建立的模型如图3所示，自外向内包括：空气域、圆筒状外壳，导体和外壳之间的圆筒形绝缘气体区域(外层绝缘气体)，圆筒状导体，导体内部的圆柱形绝缘气体区域(内层绝缘气体)。其中空气域半径为750mm，建立的模型的轴向(Z轴方向)长度为100mm。
B：把物理模型的几何文件导入网格划分软件ANSYS ICEM CFD中，进行流体和固体模块的网格划分：
将建立好的各个区域模型分别导入ICEM CFD软件中，采用分块划分网格的方式得到O型结构化网格，并对气体域定义边界层，对边界层网格进行加密。
C：在workbench中建立CFX模块，将网格导入CFX中，建立流体、固体域，并对各个计算域进行设置，具体步骤为：
C1：在workbench工作界面中建立CFX模块，将划分好的网格导入CFX模块中；
C2：创建计算中所需的新材料。固体材料的热物性视为定值，通过直接赋值的方式直接输入相应的密度、导热率、比热容等参数，实验模型中导体及外壳的材料热物性参数如下：
空气及绝缘气体视为不可压缩流体，定压比热容视为定值，密度、动力粘度、导热系数考虑温度的影响，通过CEL语言定义气体属性：
密度计算公式为：
导热系数计算公式为：
动力粘度计算公式为：
其中：ρ(T)为气体密度，ρ0(kg/m3)为0℃时气体的密度；λ(T)为导热系数，λ0(W/(m·K))为0℃时气体的导热系数；η(T)为气体动力粘度，η0(Pa·s)为0℃时气体的动力粘度；T0(K)为参考温度；S(K)为Sutherland常数；T(K)为该网格温度，在迭代时由求解器实时给出。
0℃时气体的物性参数如下：
热物性参数计算公式中的常数值为：
C3：创建流、固计算域，并为各个计算域指定相应材料。
C4：对于流体计算域，由于不同位置气体的温度不同，因此会导致不同位置气体的密度不同，从而在重力作用下气体发生自然对流。此时需要打开流体域浮力项，考虑气体的自然对流传热。其中浮力项的重力方向为Y轴方向，大小为-9.8m/s2，浮力项参考密度计算公式为：
其中ρc(kg/m3)为浮力项参考密度；Tc(K)为预期的气体域平均温度。
通过打开气体域辐射项，考虑模型的辐射传热。考虑到气体发射和吸收辐射的能力，将气体视为热辐射的透明体，因此辐射模型选择Discrete Transfer模型，辐射模式为Surface to Surface，通过在边界上设置辐射面的辐射率计算模型的辐射传热。通过这种方法考虑模型的辐射散热，将大大简化辐射散热边界条件的施加。
D：在固体计算域上施加考虑集肤效应的焦耳热载荷，作为发热源，具体步骤为：
D1：考虑集肤效应的影响，通过CEL语言定义导体和外壳的电阻计算公式为：
其中：Ri(Ω/m)为导体或外壳电阻；Kf为集肤效应系数；ρ20(Ω·m)为20℃时导体或外壳的直流电阻率；α20为电阻温度系数；Si(mm2)为导体或外壳的横截面积；Ti(℃)为导体或外壳的运行温度，在迭代时由求解器实时给出。公式中的参数如下：
| ρ<sub>20</sub>/(Ω·m) | α<sub>20</sub> | |
| 导体 | 0.029 | 0.004 |
| 外壳 | 0.033 | 0.004 |
通过CEL语言定义集肤效应系数计算公式为：
其中：C(mm)为导体或外壳壁厚；D(mm)为导体或外壳的外径。
D2：通过CEL语言定义导体或外壳每米损耗功率为：
P＝I2Ri
其中：P(W/m)，为导体或外壳的每米损耗，即每米的焦耳热功率；I(A)，在导体中为GIL额定电流，在外壳中为感应电流，当GIL长度小于20m时，感应电流为额定电流的0.95，当GIL长度大于20m时，感应电流等于GIL额定电流。在实验中，导体和外壳是直接连通的，所以流过两者的电流大小相同。
通过CEL语言定义导体或外壳单位体积的焦耳热功率为：
其中：PV(W/m3)，为导体或外壳单位体积的焦耳热功率。
D3：通过CEL语言建立好体积焦耳热功率计算公式后，在导体及外壳计算域中设置子域，将导体及外壳单位体积的焦耳热功率施加在子域上作为生热源。
E：对流、固计算域进行耦合，并施加边界条件，然后在CFX模块中进行温度场计算，具体步骤为：
E1：将流体域与固体域交界面进行耦合，考虑流固体域之间的换热，需要在耦合面上打开传热选项；同时由于流、固体域网格划分的密度不同，所以耦合面上网格的连接方式选择GGI连接方式。
E2：在此实验中，由于实验环境为室内无风速条件，所以环境空气的流动状态为由于受热而发生自然的对流。此时空气域选用层流计算模型，空气两端面设置为开口边界条件。
空气外边界面设置为无滑移等温壁面，温度为环境温度。空气与外壳交界面辐射系数根据外壳外表面的具体情况确定，在此实验中，外壳外表面的辐射系数为0.8；考虑环境吸收外壳发出的所有辐射，所以空气外边界面辐射系数设置为1。
E3：绝缘气体处于封闭空间内，只因受热不均发生自然对流，所以气体流动状态为层流，选用层流计算模型。由于只对部分管道进行了建模仿真，所以两部分绝缘气体域的端面都设置为对称面。由于导体内壁面温差不大，所以不考虑其辐射传热，所以不启用内层绝缘气体域的热辐射选项。由于导体与外壳间的温差较大，需要考虑两者之间的辐射传热，所以需要使用外层绝缘气体域的热辐射选项，具体设置如C3中所述，外层绝缘气体边界的辐射系数根据外壳内表面及导体外表面的具体情况确定，在此实验中，外壳内表面辐射系数为0.9，导体外表面辐射系数为0.9。
E4：设置求解器：包括迭代步数、时间步、收敛条件等；设置温度监测点；然后进行求解。当温度监测点的温度达到稳定时，即视为求解收敛，最后在后处理器中进行后处理，得到温度场、流场等结果。其中对第三组实验进行仿真计算得到外壳的温度场结果如图4所示：
F：将固体域网格文件导入workbench稳态或瞬态热分析模块，并将CFX的计算结果导入到热分析模块中的模型的对应边界上，再次对导体及外壳进行温度场计算，具体步骤为：
F1：在workbench工作界面中建立Finite Element Modeler(有限元建模器)模块，并在该模块中导入由ICEM CFD生成的.uns格式的固体域网格文件。
F2：在workbench工作界面中建立稳态或瞬态热分析模块，将上述有限元建模器中的网格导入热分析模块中。
F3：在热分析模块中输入材料导热率。
F4：将CFX模块中的固体边界的温度计算结果导入到热分析模块中的模型的对应边界上，然后重新进行导体或外壳的温度场计算。
G：将热分析模块中的计算结果导入workbench的静态或瞬态结构分析模块，计算热变形，具体步骤为：
G1：在workbench工作界面建立静态或瞬态结构分析模块，将有限元建模器中的网格导入结构分析模块中。
G2：在结构分析模块中输入材料的密度、热膨胀系数、零热应力的参考温度。对于此实验，零热应力参考温度为每组实验条件下的环境温度，材料的热膨胀系数如下：
| 材料 | 热膨胀系数/(1/℃) | |
| 导体 | A6063 | 2.34*10<sup>-5</sup> |
| 外壳 | A5005 | 2.38*10<sup>-5</sup> |
G3：将热分析模块中计算的温度导入结构分析模块；由于计算模型为部分管道模型，所以在管道端面不施加位移约束，只需在分析选项中打开弱弹簧选项即可；最后进行热应力、热变形计算。其中对第三组实验进行仿真计算得到外壳的形变场结果如图5所示。
以上即为本方法的完整过程，
使用本方法对Minaguchi等人的三组实验分别进行仿真计算，得到的温度结果对比如下：
可以看到，使用本方法进行仿真计算，误差最大为2.13％，说明本方法能有效求解GIL温度场的分布，并以温度场计算结果为依据进行应力、应变计算是可行的。
Claims (7)
1.一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法，其特征在于，包括以下步骤：
A：根据结构特性，对模型进行简化，使用SolidWorks软件建立包含绝缘气体及环境空气的计算模型；
B：把物理模型的几何文件导入网格划分软件ANSYS ICEM CFD中，进行流体和固体模块的网格划分；
C：在workbench中建立CFX模块，将网格导入CFX中，建立流体、固体域，并对各个计算域进行设置；
D：在固体计算域上施加考虑集肤效应的焦耳热载荷；
E：对流、固计算域进行耦合，并施加边界条件，然后在CFX模块中进行温度场计算；
F：将固体域网格文件导入workbench稳态或瞬态热分析模块，并将CFX的计算结果导入到热分析模块中的模型的对应边界上，再次对导体及外壳进行温度场计算；
G：将热分析模块中的计算结果导入workbench的静态或瞬态结构分析模块，计算热变形；
所述步骤C包括以下步骤：
C1：在workbench工作界面中放入CFX模块，将划分好的网格导入CFX中；
C2：创建计算中所需的新材料，固体材料的热物性视为定值，通过直接赋值的方式直接输入相应的导热率、比热容参数；气体为不可压缩流体，定压比热容视为定值，密度、动力粘度、导热系数考虑温度的影响，通过CEL语言定义气体属性：
密度计算公式为：
导热系数计算公式为：
动力粘度计算公式为：
其中：ρ(T)为气体密度，单位为kg/m3；ρ0为0℃时气体的密度，单位为kg/m3；λ(T)为导热系数，λ0为0℃时气体的导热系数；η(T)为气体动力粘度，单位为Pa·s；η0为0℃时气体的动力粘度；T0为参考温度，273K；S为Sutherland常数，单位为K；T为该网格温度，单位为K，在迭代时由求解器实时给出；
当考虑空气湿度的影响时，可以采用相同方法建立湿空气的密度、导热系数、动力粘度公式，然后赋予新材料；
C3：创建流、固计算域，并为各个计算域指定相应材料；
C4：对于流体计算域，需要通过添加浮力项考虑气体由于受热产生的密度变化发生的自然对流，其中浮力项参考密度计算公式为：
其中ρc为浮力项参考密度，单位为kg/m3；Tc为环境温度，单位为K；
同时通过打开气体域辐射项，设置辐射面的辐射率考虑结构的辐射散热，只考虑气体界面上的辐射，不考虑气体内部对辐射的吸收。
2.根据权利要求1所述的一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法，其特征在于，所述步骤A包括以下步骤：
A1：对模型进行简化，由于每一段管道两端的触头处的接触良好，导电性能优良，所以可以忽略触头处的热损耗对管道温度的影响，所以模型只截取部分轴向长度的管道来计算温度场的分布；
A2：使用SolidWorks进行建模，建立的模型包括空气域、圆筒状导体、圆筒状外壳、导体内部的圆柱形绝缘气体区域即内层绝缘气体以及导体和外壳之间的圆筒形绝缘气体区域即外层绝缘气体。
3.根据权利要求1所述的一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法，其特征在于，所述步骤B网格划分具体为：
将建立好的各个区域模型分别导入ICEM CFD软件中，采用分块划分网格的方式得到O型及H型结构化网格，并对气体域定义边界层，对边界层网格进行加密。
4.根据权利要求1所述的一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法，其特征在于，所述步骤D包括以下步骤：
D1：考虑集肤效应的影响，导体和外壳的电阻计算公式为：
其中：Ri为导体或外壳电阻，单位为Ω/m；Kf为集肤效应系数；ρ20为20℃时导体或外壳的直流电阻率，单位为Ω·m；α20为电阻温度系数；Si为导体或外壳的横截面积，单位为mm2；Ti为导体或外壳的运行温度，单位为℃，在迭代时由求解器实时给出；
集肤效应系数计算公式为：
其中：C为导体或外壳壁厚，单位为mm；D为导体或外壳的外径，单位为mm；
D2：导体或外壳每米损耗功率为：
P＝I2Ri
其中：P为导体或外壳的每米损耗，即每米的焦耳热功率，单位为W/m；I为在导体中为GIL额定电流，单位为A，在外壳中为感应电流，当GIL长度小于20m时，感应电流为额定的0.95，当GIL长度大于20m时，感应电流等于GIL额定电流；
导体或外壳单位体积的焦耳热功率为：
其中：PV，为导体或外壳单位体积的焦耳热功率，单位为W/m3；
D3：通过CEL语言定义导体和外壳单位体积的焦耳热功率，并将其作为生热率载荷施加在导体及外壳上。
5.根据权利要求1所述的一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法，其特征在于，所述步骤E包括以下步骤：
E1：将流体域与固体域交界面进行耦合，包含热量交换；
E2：考虑空气由于受热发生自然对流时，空气域选用层流模型，空气两端面设置为开口边界条件；考虑风速时，根据廊管尺寸及空气流速综合计算空气流动的雷诺数，选择合适的流动模型，然后设置空气入口流速、温度，出口压强；空气外边界面设置为壁面，温度为环境温度，考虑环境吸收外壳的所有辐射，所以空气外边界面辐射系数设置为1；空气与输电管道交界面辐射系数根据输电管道外壳外表面的具体情况确定；
E3：绝缘气体由于只发生自然对流，所以选用层流模型；端面设置为对称面；内层绝缘气体不考虑辐射；外层绝缘气体边界的辐射系数根据输电管道外壳内表面及导体外表面的具体情况确定；
E4：设置求解器：包括迭代步数、时间步、收敛条件；设置温度监测点；然后进行求解，得到温度场的分布结果。
6.根据权利要求1所述的一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法，其特征在于，所述步骤F包括以下步骤：
F1：在workbench工作界面中建立有限元建模器Finite Element Modeler模块，并在该模块中输入由ICEM CFD生成的.uns格式的固体域网格文件；
F2：在workbench工作界面中建立稳态或瞬态热分析模块，将上述有限元建模器中的网格导入热分析模块中；
F3：在热分析模块中输入材料导热率；
F4：将CFX模块中的固体边界的温度计算结果导入到热分析模块中的模型的对应边界上，然后重新进行导体或外壳的温度场计算。
7.根据权利要求1所述的一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法，其特征在于，所述步骤G包括以下步骤：
G1：在workbench工作界面建立静态或瞬态结构分析模块，将有限元建模器中的网格导入结构分析模块中；
G2：在结构分析模块中输入材料的密度、热膨胀系数、零热应力的参考温度；
G3：施加的约束条件为在导体或外壳的一端截面施加位移约束，仅约束轴向位移，并将热分析模块中计算的温度导入结构分析模块；同时在分析选项中打开弱弹簧选项，然后进行热变形计算。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711149404.7A CN108052697B (zh) | 2017-11-18 | 2017-11-18 | 一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711149404.7A CN108052697B (zh) | 2017-11-18 | 2017-11-18 | 一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108052697A CN108052697A (zh) | 2018-05-18 |
| CN108052697B true CN108052697B (zh) | 2020-07-10 |
Family
ID=62119219
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711149404.7A Expired - Fee Related CN108052697B (zh) | 2017-11-18 | 2017-11-18 | 一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108052697B (zh) |
Families Citing this family (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109376487B (zh) * | 2018-12-04 | 2022-12-06 | 国网青海省电力公司电力科学研究院 | 一种高海拔地区gis长母线结构温差变形的计算方法 |
| CN110781589B (zh) * | 2019-10-18 | 2023-01-20 | 广东电网有限责任公司广州供电局 | 检测气体绝缘金属封闭输电线路搭接处过热故障的方法 |
| CN110781630A (zh) * | 2019-11-28 | 2020-02-11 | 青海电研科技有限责任公司 | 一种gis设备管线结构温差应力变形在线监测方法 |
| CN111460706B (zh) * | 2020-03-19 | 2023-06-06 | 国网江苏省电力有限公司电力科学研究院 | Gil管廊在线监测及温度状态判别方法、系统 |
| CN111400934B (zh) * | 2020-04-15 | 2023-04-07 | 天华化工机械及自动化研究设计院有限公司 | 基于workbench低温碳化炉氧含量分布模拟方法 |
| CN111754633B (zh) * | 2020-06-09 | 2024-09-27 | 中国石油天然气集团有限公司 | 一种含汞天然气泄漏造成的土壤污染面积的计算方法 |
| CN111680432B (zh) * | 2020-07-17 | 2023-04-07 | 天华化工机械及自动化研究设计院有限公司 | 一种基于workbench的低温碳化炉多耦合场应力分布的模拟方法 |
| CN112001101A (zh) * | 2020-07-20 | 2020-11-27 | 中国南方电网有限责任公司超高压输电公司广州局 | 基于有限元的gil电-磁-热-流-力多物理场耦合仿真建模方法 |
| CN112318888A (zh) * | 2020-11-05 | 2021-02-05 | 南京航空航天大学 | 一种电流辅助铆接的温度轴向分布预测装置及方法 |
| CN112528363B (zh) * | 2020-11-17 | 2023-04-25 | 栗怀广 | 建立位移响应预测模型的方法、装置及电子设备 |
| CN114235881A (zh) * | 2021-11-26 | 2022-03-25 | 南方电网科学研究院有限责任公司 | 一种环保绝缘气体导热与耐热性能的评价方法及装置 |
| CN114218833B (zh) * | 2021-12-16 | 2023-11-10 | 西北工业大学太仓长三角研究院 | 一种二级轻气炮内流场性能预测方法及系统 |
| CN115906555B (zh) * | 2022-10-12 | 2023-09-19 | 中建二局安装工程有限公司 | 基于comsol的焊接过程中残余应力及其漏磁场的预测方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101819072A (zh) * | 2009-05-12 | 2010-09-01 | 中国水电顾问集团西北勘测设计研究院 | 超高压气体绝缘金属封闭输电线路温度监测方法 |
| GB2537932A (en) * | 2015-05-01 | 2016-11-02 | Fujitsu Ltd | Method and apparatus for use in thermal coupled analysis |
| CN106124928A (zh) * | 2016-06-21 | 2016-11-16 | 中国电力科学研究院 | 特高压交流混合气体gil线段温度监测方法 |
| CN106202610B (zh) * | 2016-06-22 | 2019-05-14 | 华南理工大学 | 一种基于ansys cfx的架空线径向温度场仿真方法 |
-
2017
- 2017-11-18 CN CN201711149404.7A patent/CN108052697B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN108052697A (zh) | 2018-05-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108052697B (zh) | 一种基于workbench的廊管GIL三维温度场及胀缩形变计算方法 | |
| Biglarian et al. | A numerical model for transient simulation of borehole heat exchangers | |
| Yadav et al. | Modeling and Simulation of Turbulent Flows through a Solar Air Heater Having Square‐Sectioned Transverse Rib Roughness on the Absorber Plate | |
| CN105045966A (zh) | 一种变压器内部温度场的混合计算方法 | |
| Kim et al. | A numerical study of the effect of a hybrid cooling system on the cooling performance of a large power transformer | |
| Rejeb et al. | Statistical study of a solar nanofluid photovoltaic thermal collector performance using response surface methodology | |
| Sedaghat et al. | Enhanced thermal model of power cables installed in ducts for ampacity calculations | |
| Minaei et al. | Thermal resistance capacity model for short-term borehole heat exchanger simulation with non-stiff ordinary differential equations | |
| CN106021676A (zh) | 一种基于转移矩阵的多回电缆稳态温升获取方法 | |
| Fu et al. | Numerical study of convection and radiation heat transfer in pipe cable | |
| CN108038268B (zh) | 一种廊管内gil运行温度评估方法 | |
| Dutt et al. | Thermo‐hydraulic performance investigation of solar air heater duct having staggered D‐shaped ribs: Numerical approach | |
| Lee | A modified three-dimensional numerical model for predicting the short-time-step performance of borehole ground heat exchangers | |
| Mortazavinejad et al. | Numerical investigation of two-dimensional heat transfer of an absorbing plate of a flat-plate solar collector using dual-reciprocity method based on boundary element | |
| Olakoyejo et al. | Constructal heat transfer and fluid flow enhancement optimization for cylindrical microcooling channels with variable cross‐section | |
| Jose et al. | An Analysis of the Effects of Nanofluid‐Based Serpentine Tube Cooling Enhancement in Solar Photovoltaic Cells for Green Cities | |
| Li et al. | Comparison of heat transfer characteristics between natural ester oil and mineral oil in large oil-immersed transformer | |
| Trokhaniak et al. | Research of thermal and hydrodynamic flows of heat exchangers for different air-cooling systems in poultry houses. | |
| He et al. | Applications of a dynamic three-dimensional numerical model for borehole heat exchangers | |
| Saboohi et al. | The effect of sinusoidal fins’ amplitude on the thermo-hydraulic performance of a solar air heater | |
| Alaie et al. | Investigation of thermal interaction between shallow boreholes in a GSHE using the FLS-STRCM model | |
| Wang et al. | Analysis of the Three-Dimensional Temperature Distribution and Ampacity of Forced Ventilation Multi-loop Cable Tunnel by Finite Element Method | |
| Rezvani et al. | Numerical study of entropy generation for natural convection in cylindrical cavities | |
| Fu et al. | Heat transfer simulation and temperature rapid prediction for trench laying cable | |
| Gelis et al. | Development of novel flow fields for pem fuel cells: numerical solution and experimental validation |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20190202 Address after: No. 28 Xianning West Road, Xianning stele area, Xi'an, Shaanxi Applicant after: XI'AN JIAOTONG University Applicant after: STATE GRID CORPORATION OF CHINA Applicant after: STATE GRID JIANGSU ELECTRIC POWER Co.,Ltd. Applicant after: STATE GRID JIANGSU ELECTRIC POWER Co.,Ltd. MAINTENANCE BRANCH Address before: No. 28 Xianning West Road, Xianning stele area, Xi'an, Shaanxi Applicant before: Xi'an Jiaotong University |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200710 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |