CN114462264A - 一种基于多物理场耦合的电缆接头状态评估方法及装置 - Google Patents
一种基于多物理场耦合的电缆接头状态评估方法及装置 Download PDFInfo
- Publication number
- CN114462264A CN114462264A CN202111527853.7A CN202111527853A CN114462264A CN 114462264 A CN114462264 A CN 114462264A CN 202111527853 A CN202111527853 A CN 202111527853A CN 114462264 A CN114462264 A CN 114462264A
- Authority
- CN
- China
- Prior art keywords
- field grid
- temperature
- grid node
- electromagnetic
- field
- 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.)
- Pending
Links
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/04—Power grid distribution networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
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)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本申请涉及一种基于多物理场耦合的电缆接头状态评估方法、装置、计算机设备、存储介质和计算机程序产品。该方法基于待测电缆接头的工作参数和电磁场数学模型计算得到电磁场网格节点的单位体积电磁损耗,并利用径向基函数在电磁场网格节点、温度场网格节点和应力场网格节点之间进行参量映射,并进行参量的计算更新,耦合计算直至达到迭代收敛条件时,得到待测电缆接头的内部物理量分布结果并得到待测电缆接头的状态评估结果。该方法利用径向基函数进行待测电缆接头内部多物理场网格节点数据的快速映射,实现待测电缆接头内部电磁场‑温度场‑应力场快速耦合计算,进而为电缆接头状态评估以及电缆接头缺陷故障机理提供有效途径。
Description
技术领域
本申请涉及电力技术领域,特别是涉及一种基于多物理场耦合的电缆接头 状态评估方法、装置、计算机设备、存储介质和计算机程序产品。
背景技术
近年来,随着城市配电网规模的不断扩大和架空线入地工程的不断加速, 配电电缆的数量迅速上涨,电缆的健康稳定已成为城市电网安全可靠运行的重 要保证。
在电力系统中,电缆之间常常用于电力传输和电能分配,受电缆长度和运 行维护影响,两段电缆常常使用电缆中间接头连接。电缆接头由现场人工制作 安装,其质量受现场施工环境和施工人员技术素质影响较大。因此,电缆接头 在运行过程中难免存在一些潜伏性的缺陷。
电缆接头在缺陷下的运行往往是一个涉及电场集中、热场集中、应力场集 中的复杂过程,很难通过单一物理场进行评估,该过程会使得电缆接头绝缘状 态加速恶化,并最终导致电缆接头事故的发生,影响电网安全可靠性。
发明内容
基于此,有必要针对上述技术问题,提供一种能够基于多物理场耦合得到 对电缆接头较为准确的状态评估结果的基于多物理场耦合的电缆接头状态评估 方法、装置、计算机设备、计算机可读存储介质和计算机程序产品。
第一方面,本申请提供了一种基于多物理场耦合的电缆接头状态评估方法。 方法包括:
基于待测电缆接头的工作参数和电磁场数学模型计算得到电磁场网格节点 的单位体积电磁损耗;
由计算得到的电磁场网格节点的单位体积电磁损耗开始,利用径向基函数 在电磁场网格节点与温度场网格节点之间进行参量映射、以及利用径向基函数 在温度场网格节点和应力场网格节点之间进行参量映射,利用电磁场数学模型 基于电磁场网格节点接收到的参量计算更新得到电磁场网格节点的单位体积电 磁损耗,利用温度场数学模型基于温度场网格节点接收到的参量计算更新得到 温度场网格节点的温度,利用应力场数学模型基于应力场网格节点接收到的参 量计算更新得到应力场网格节点的应力张量;
在电磁场、温度场和应力场之间进行迭代的参量映射以及耦合计算直至达 到迭代收敛条件时,得到待测电缆接头的内部物理量分布结果并得到待测电缆 接头的状态评估结果,内部物理量分布结果包括内部电磁损耗分布结果、内部 温度分布结果和内部应力分布结果中的至少一种。
在其中一个实施例中,利用径向基函数在电磁场网格节点与温度场网格节 点之间进行参量映射,包括:利用径向基函数将电磁场网格节点的单位体积电 磁损耗映射得到温度场网格节点的单位体积电磁损耗、并利用温度场数学模型 基于温度场网格节点的单位体积电磁损耗计算更新得到温度场网格节点的温 度;以及,利用径向基函数将温度场网格节点的温度映射得到电磁场网格节点 的温度、并利用电磁场数学模型基于电磁场网格节点的温度计算更新得到电磁 场网格节点的单位体积电磁损耗。
在其中一个实施例中,利用径向基函数在温度场网格节点和应力场网格节 点之间进行参量映射,包括:利用径向基函数将温度场网格节点的温度映射得 到应力场网格节点的温度、并利用应力场数学模型基于应力场网格节点的温度 计算更新得到应力场网格节点的应力张量;利用径向基函数将应力场网格节点 的应力张量映射得到温度场网格节点的应力张量、并利用温度场数学模型基于 温度场网格节点的应力张量计算更新得到温度场网格节点的温度。
在其中一个实施例中,利用径向基函数将第一物理场网格节点中的参量映 射得到第二物理场网格节点中的参量表示为:
AS=HSfAf;
其中,表示n个第二物理场网格节点中的参量,对于任意的 1≤i≤n,表示任意第i个第二物理场网格节点si中的参量;表示 n个第一物理场网格节点中的参量,对于任意的1≤j≤n,表示任意第j个第一 物理场网格节点fj中的参量;
表示任意第j个第一物理场网格节点fj的坐标,表示任意第i个第二物 理场网格节点si中的参量;对于任意的1≤k≤n,表示任意的第j个第一物理 场网格节点fj和第k个第一物理场网格节点fk之间的径向基函数关系;表示 任意第j个第一物理场网格节点fj与第i个第二物理场网格节点si之间的径向基 函数关系;
其中,第一物理场网格节点和第二物理场网格节点中一个是电磁场网格节 点、另一个是温度场网格节点,或者,第一物理场网格节点和第二物理场网格 节点中一个是温度场网格节点、另一个是应力场网格节点。
在其中一个实施例中,电磁场数学模型包括有外加激励源存在的电流区域 的第一控制方程和非电流区域的第二控制方程,且第一控制方程为 第二控制方程为其中,A表示矢量 磁位,μ表示材料磁导率,σ表示材料电导率,Js表示外加电流密度,ω为角 频率;基于待测电缆接头的工作参数和电磁场数学模型计算得到电磁场网格节 点的单位体积电磁损耗,包括:建立待测电缆接头的几何模型并离散为若干个 电磁场网格节点;基于待测电缆接头的工作参数计算得到各个电磁场网格节点 内的外加电流密度Js和材料电导率σ并代入对应的电磁场数学模型中,计算得 到矢量磁位A;按照计算得到总电流密度J,并按照计 算得到各个电磁场网格节点的单位体积电磁损耗Qv。
在其中一个实施例中,待测电缆接头的工作参数包括施加电流数据以及温 度传感数据,基于待测电缆接头的工作参数计算得到各个电磁场网格节点内的 外加电流密度Js和材料电导率σ,包括:将施加电流数据作为外加电流密度Js, 将温度传感数据按照计算得到当前的材料电导率σ,σ20是20℃ 下的材料电导率,α是温度系数。
第二方面,本申请还提供了一种基于多物理场耦合的电缆接头状态评估装 置。装置包括:
初始计算模块,用于基于待测电缆接头的工作参数和电磁场数学模型计算 得到电磁场网格节点的单位体积电磁损耗;
迭代计算模块,用于由计算得到的电磁场网格节点的单位体积电磁损耗开 始,利用径向基函数在电磁场网格节点与温度场网格节点之间进行参量映射、 以及利用径向基函数在温度场网格节点和应力场网格节点之间进行参量映射, 利用电磁场数学模型基于电磁场网格节点接收到的参量计算更新得到电磁场网 格节点的单位体积电磁损耗,利用温度场数学模型基于温度场网格节点接收到 的参量计算更新得到温度场网格节点的温度,利用应力场数学模型基于应力场 网格节点接收到的参量计算更新得到应力场网格节点的应力张量;
结果获取模块,用于在电磁场、温度场和应力场之间进行迭代的参量映射 以及耦合计算直至达到迭代收敛条件时,得到待测电缆接头的内部物理量分布 结果并得到待测电缆接头的状态评估结果,内部物理量分布结果包括内部电磁 损耗分布结果、内部温度分布结果和内部应力分布结果中的至少一种。
第三方面,本申请还提供了一种计算机设备。计算机设备包括存储器和处 理器,存储器存储有计算机程序,处理器执行计算机程序时实现第一方面提供 的基于多物理场耦合的电缆接头状态评估方法的步骤。
第四方面,本申请还提供了一种计算机可读存储介质。计算机可读存储介 质,其上存储有计算机程序,计算机程序被处理器执行时实现第一方面提供的 基于多物理场耦合的电缆接头状态评估方法的步骤。
第五方面,本申请还提供了一种计算机程序产品。计算机程序产品,包括 计算机程序,该计算机程序被处理器执行时实现第一方面提供的基于多物理场 耦合的电缆接头状态评估方法的步骤。
上述基于多物理场耦合的电缆接头状态评估方法、装置、计算机设备、存 储介质和计算机程序产品,通过确定待测电缆接头内部电磁场-温度场-应力场耦 合作用形式,利用径向基函数进行待测电缆接头内部多物理场网格节点数据的 快速映射,实现待测电缆接头内部电磁场-温度场-应力场快速耦合计算,进而为 电缆接头状态评估以及电缆接头缺陷故障机理提供有效途径。
附图说明
图1为一个实施例中基于多物理场耦合的电缆接头状态评估方法的流程示 意图。
图2为另一个实施例中基于多物理场耦合的电缆接头状态评估方法的流程 示意图。
图3为一个实施例中基于多物理场耦合的电缆接头状态评估装置的结构框 图。
图4为另一个实施例中基于多物理场耦合的电缆接头状态评估装置的结构 框图。
图5为一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实 施例,对本申请进行进一步详细说明。应当理解,此处描述的具体实施例仅仅 用以解释本申请,并不用于限定本申请。
在一个实施例中,如图1所示,提供了一种基于多物理场耦合的电缆接头 状态评估方法,以该方法应用于终端进行举例说明,可以理解的是,该方法也 可以应用于服务器,还可以应用于包括终端和服务器的系统,并通过终端和服 务器的交互实现,包括以下步骤:
步骤102,基于待测电缆接头的工作参数和电磁场数学模型计算得到电磁场 网格节点的单位体积电磁损耗。待测电缆接头的工作参数用于反映待测电缆接 头实际的工作情况,主要包含初始参数和运行参数,初始参数反映待测电缆接 头的初始基础特性,运行参数反映待测电缆接头在运行过程中的各项物理参数。 在一个实施例中,初始参数包括待测电缆接头电压等级、安装时间和耐压情况 中的至少一种。在一个实施例中,运行参数包括施加电流数据以及温度传感数 据,施加电流数据是待测电缆接头在实际工作时内部的电流,温度传感数据是 待测电缆接头在实际工作时的现场温度。由于待测电缆接头在实际工作时内部 的电流不是恒定值,因此一般取一段时间内的均值作为施加电流数据。
电磁场数学模型是单一电磁场的数学模型、根据麦克斯韦方程组变换得到。 为了计算得到电磁场网格节点的单位体积电磁损耗,首先还需要建立待测电缆 接头的几何模型并离散为若干个电磁场网格节点,具体离散得到的电磁场网格 节点可以是三角形、四边形等形状的有限元网格。而在建立待测电缆接头的几 何模型时,需要结合待测电缆接头实际的工作情况,也即基于待测电缆接头的 工作参数来建立待测电缆接头的几何模型并离散得到电磁场网格节点,具体的 建模离散网格的方法,可以结合现有的有限元软件实现,本申请不展开详述。
步骤104,由计算得到的电磁场网格节点的单位体积电磁损耗开始,利用径 向基函数在电磁场网格节点与温度场网格节点之间进行参量映射、以及利用径 向基函数在温度场网格节点和应力场网格节点之间进行参量映射。
任意两个物理场网格节点之间的参量映射包括双向的映射,也即对于电磁 场网格节点和温度场网格节点之间的参量映射,包括电磁场网格节点向温度场 网格节点的参量映射,以及温度场网格节点向电磁场网格节点的参量映射。同 样的,对于温度场网格节点和应力场网格节点之间的参量映射,包括温度场网 格节点向应力场网格节点的参量映射,以及应力场网格节点向温度场网格节点 的参量映射。而上述双向映射过程最初由步骤102计算得到的电磁场网格节点 的单位体积电磁损耗开始,按照电磁场-温度场-应力场的耦合方式进行耦合映 射,也即首先将步骤102计算得到的电磁场网格节点的单位体积电磁损耗映射 向温度场网格节点进行参量映射,再进行后续的不同物理场之间的迭代映射。
在参量映射过程中,接收到映射的参量的物理场还基于自身物理场的数学 模型进行参量的更新计算,包括:利用电磁场数学模型基于电磁场网格节点接 收到的参量计算更新得到电磁场网格节点的单位体积电磁损耗,利用温度场数 学模型基于温度场网格节点接收到的参量计算更新得到温度场网格节点的温 度,利用应力场数学模型基于应力场网格节点接收到的参量计算更新得到应力 场网格节点的应力张量。
温度场数学模型是单一温度场的数学模型,根据傅里叶传热定律和能量守 恒定律得到。应力场数学模型是单一应力场的数学模型,根据电缆接头内部应 力场分析方程组推导得到。
步骤106,在电磁场、温度场和应力场之间进行迭代的参量映射以及耦合计 算直至达到迭代收敛条件时,得到待测电缆接头的内部物理量分布结果并得到 待测电缆接头的状态评估结果。
其中迭代收敛条件为预先设定的条件,可以是迭代次数达到预设阈值,或 者是同样的参量在迭代过程中的取值误差在误差范围内,无论具体如何设定迭 代收敛条件,当达到迭代收敛条件时,可以获得较为准确的电磁场、温度场和 应力场的内部物理量分布结果。
其中,内部物理量分布结果包括内部电磁损耗分布结果、内部温度分布结 果和内部应力分布结果中的至少一种,可以利用三个物理场的分布结果中的任 意一种或多种的结合得到最终的状态评估结果,状态评估结果用于指示待测电 缆接头的运行状态是否存在故障风险。具体得到状态评估结果的方法有多种, 比如以仅根据内部温度分布结果得到状态评估结果为例,比如可以是在内部温 度分布结果指示存在温度达到温度阈值时得到指示有故障风险的状态评估结 果,在所有温度都未达到温度阈值时得到指示没有故障风险的状态评估结果。
上述基于多物理场耦合的电缆接头状态评估方法中,通过确定待测电缆接 头内部电磁场-温度场-应力场耦合作用形式,利用径向基函数进行待测电缆接头 内部多物理场网格节点数据的快速映射,实现待测电缆接头内部电磁场-温度场- 应力场快速耦合计算,进而为电缆接头状态评估以及电缆接头缺陷故障机理提 供有效途径。
在上述实施例中,主要使用径向基函数实现待测电缆接头内部多物理场网 格节点的数据快速映射,该方法的核心是根据一系列采样点上的输入-输出关系, 利用一定的径向基函数构造输入-输出系统逼近真实输入输出关系的近似解析 解。可选的,径向基函数是距离的函数,设x=(x1,x2,…,xn)∈Rn,则其欧几里得 长度可以表示为则定义径向基函数为φ(r),即距离为式中||·||为n维空间的欧式范数,c为径向基函数中心,函数在 某点的函数值取决于该点到径向基函数中心的距离。
利用径向基函数进行物理场网格数据快速映射的方法是一种局部插值的方 法,通过在已知物理场网格所有单元内搜索,寻找未知物理场网格每一个待求 物理场网格节点在已知物理场网格上对应的主单元。径向基函数表示如下:
上式中,x1,x2,…,xn是在d维欧几里德空间内给定的一组位置不同的点; f(x)是关于未知物理场网格节点x的连续函数;φ代表径向基函数;λi代表每个 中心点的权重。p(x)一般采用线性多项式的形式为p(x)=a0+a1x+a2y。
基于径向基函数的原理和插值公式,可以得到在利用径向基函数将第一物 理场网格节点中的参量映射得到第二物理场网格节点中的参量时的表达式:
AS=HSfAf;
其中,表示n个第二物理场网格节点中的参量,对于任意的1≤i≤n,表示任意第i个第二物理场网格节点si中的参量。表示 n个第一物理场网格节点中的参量,对于任意的1≤j≤n,表示任意第j个第一 物理场网格节点fj中的参量。
其中,表示任意第j个第一物理场网格节点fj的坐标,表示任意第i个 第二物理场网格节点si中的参量。对于任意的1≤k≤n,表示任意的第j个第 一物理场网格节点fj和第k个第一物理场网格节点fk之间的径向基函数关系。 表示任意第j个第一物理场网格节点fj与第i个第二物理场网格节点si之间 的径向基函数关系。
而在上述实施例中,第一物理场网格节点和第二物理场网格节点中一个是 电磁场网格节点、另一个是温度场网格节点,实现电磁场网格节点向温度场网 格节点的映射,或者实现温度场网格节点向电磁场网格节点的映射。或者,第 一物理场网格节点和第二物理场网格节点中一个是温度场网格节点、另一个是 应力场网格节点,实现温度场网格节点向应力场网格节点的映射,或者实现应 力场网格节点向温度场网格节点的映射。
在一个实施例中,建立的三个物理场的单一的数学模型分别如下:
(1)电磁场数学模型。
电磁场数学模型包括有外加激励源存在的电流区域的第一控制方程和非电 流区域的第二控制方程,引入矢量磁位A,两个控制方程分别对矢量磁位A的 控制方程。
有外加激励源存在的电流区域的第一控制方程为:
非电流区域的第二控制方程为:
其中,A表示矢量磁位。μ表示材料磁导率、是材料特性参数,可以认为 是已知量。ω为角频率、可以认为是已知量。σ表示材料电导率、会随着发生 变化。Js表示外加电流密度,与电流数据相关,当没有外加电流时,Js=0。
(2)温度场数学模型。
温度场数学模型包括第一温度方程和第二温度方程,第一温度方程反映受 电磁场影响所导致的温度变化情况,第二温度方程反映受应力场影响所导致的 温度变化情况。
第一温度方程在直角坐标系中可写成:
其中,ρ为材料密度,c为材料比热容,可以认为是已知的材料特性参数。 t为时间。λx、λy和λz分别为材料沿x、y和z方向的导热率,对于同性材料有 λx=λy=λz。T为待求的温度。Qv为物体内部单位体积产生的热量。
在该方程中,物体温度升高所需要的热量应与外界传入物体内部的热量及 物体内部热源产生的热量之和相等。上述方程中第1项为单位时间内温度升高 所需要的热量;第2、3和4项分别为外界由x、y和z方向传入物体内部的热 量;第5项为物体内部单位体积产生的热量。
第二温度方程可以写成为:
εTh=βΔT=β(T-Tref);
εTh表示应力场的应力张量方程组,β表示线膨胀系数,T为温度,ΔT是 温度T与参考温度Tref相比的温度变化量。
(3)应力场数学模型
应力场数学模型的表达式为:
其中,ε为应力张量,u为位移,D为应变系数,β(T)为温度T下的线 膨胀系数,ΔT是温度T与参考温度Tref相比的温度变化量,δ是狄拉克函数, m、n、p、q表示方向,且分别为直角坐标系中的x、y、z中的任意一个。
基于上述三个物理场的单一的数学模型,在一个实施例中,该方法包括如 下步骤:
步骤202,建立待测电缆接头的几何模型并离散为若干个电磁场网格节点。 该步骤的具体操作可以参考上述步骤102部分的介绍,该实施例不再详述。
步骤204,基于待测电缆接头的工作参数计算得到各个电磁场网格节点内的 外加电流密度Js和材料电导率σ,并代入对应的电磁场数学模型中,计算得到 矢量磁位A。
具体的,将获取到的工作参数中的施加电流数据作为外加电流密度Js。将 获取到的工作参数中的温度传感数据按照计算得到当前温度下 的材料电导率σ,σ20是20℃下的材料电导率、为一个作为参考值的已知量。α 是温度系数,表示材料电导率随着温度的变化关系,为已知量。
步骤210,利用径向基函数将电磁场网格节点的单位体积电磁损耗映射Qv得 到温度场网格节点的单位体积电磁损耗Qv。并将其作为热源,代入温度场数学 模型中的第一温度方程中,将其作为物体内部单位体积产生的热量Qv这一项, 可以利用温度场数学模型计算得到各个温度场网格节点的温度T,也即利用温 度场数学模型基于温度场网格节点的单位体积电磁损耗可以计算更新得到温度 场网格节点的温度。
在得到各个温度场网格节点的温度T后,可以利用径向基函数将温度场网 格节点的温度映射得到电磁场网格节点的温度,也可以利用径向基函数将温度 场网格节点的温度映射得到应力场网格节点的温度,从而使得迭代映射的顺序 可以有多种,该实施例以电磁场-温度场-应力场-温度场-电磁场的依次映射顺序 为例进行说明。
步骤212,利用径向基函数将温度场网格节点的温度映射得到应力场网格节 点的温度,利用应力场数学模型基于应力场网格节点的温度计算更新得到应力 场网格节点的应力张量。也即将映射得到的应力场网格节点的温度代入应力场 数学模型中的温度T,可以计算得到应力张量ε。
步骤214,利用径向基函数将应力场网格节点的应力张量映射得到温度场网 格节点的应力张量、并利用温度场数学模型基于温度场网格节点的应力张量计 算更新得到温度场网格节点的温度。
具体的,在映射得到的温度场网格节点的应力张量后,代入温度场数学模 型中的第二温度方程中作为应力张量方程组εTh这一项,从而计算更新得到温度 场网格节点的温度。
步骤216,利用径向基函数将温度场网格节点的温度映射得到电磁场网格节 点的温度、并利用电磁场数学模型基于电磁场网格节点的温度计算更新得到电 磁场网格节点的单位体积电磁损耗。
在映射得到的电磁场网格节点的温度后,重新按照计算得到 当前温度下的材料电导率σ,并代入电磁场数学模型中计算得到电磁场网格节 点的单位体积电磁损耗。然后可以重复执行上述步骤210-216进行循环迭代,直 至达到迭代收敛条件。
应该理解的是,虽然如上的各实施例所涉及的流程图中的各个步骤按照箭 头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。 除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤 可以以其它的顺序执行。而且,如上的各实施例所涉及的流程图中的至少一部 分步骤可以包括多个步骤或者多个阶段,这些步骤或者阶段并不必然是在同一 时刻执行完成,而是可以在不同的时刻执行,这些步骤或者阶段的执行顺序也 不必然是依次进行,而是可以与其它步骤或者其它步骤中的步骤或者阶段的至 少一部分轮流或者交替地执行。
基于同样的发明构思,本申请实施例还提供了一种用于实现上述所涉及的 基于多物理场耦合的电缆接头状态评估方法的基于多物理场耦合的电缆接头状 态评估装置。该装置所提供的解决问题的实现方案与上述方法中所记载的实现 方案相似,故下面所提供的一个或多个基于多物理场耦合的电缆接头状态评估 装置实施例中的具体限定可以参见上文中对于基于多物理场耦合的电缆接头状 态评估方法的限定,在此不再赘述。
在一个实施例中,如图3所示,提供了一种基于多物理场耦合的电缆接头 状态评估装置,包括:初始计算模块、迭代计算模块和结果获取模块,其中:
初始计算模块310,用于基于待测电缆接头的工作参数和电磁场数学模型计 算得到电磁场网格节点的单位体积电磁损耗;
迭代计算模块320,用于由计算得到的电磁场网格节点的单位体积电磁损耗 开始,利用径向基函数在电磁场网格节点与温度场网格节点之间进行参量映射、 以及利用径向基函数在温度场网格节点和应力场网格节点之间进行参量映射, 利用电磁场数学模型基于电磁场网格节点接收到的参量计算更新得到电磁场网 格节点的单位体积电磁损耗,利用温度场数学模型基于温度场网格节点接收到 的参量计算更新得到温度场网格节点的温度,利用应力场数学模型基于应力场 网格节点接收到的参量计算更新得到应力场网格节点的应力张量;
结果获取模块330,用于在电磁场、温度场和应力场之间进行迭代的参量映 射以及耦合计算直至达到迭代收敛条件时,得到待测电缆接头的内部物理量分 布结果并得到待测电缆接头的状态评估结果,内部物理量分布结果包括内部电 磁损耗分布结果、内部温度分布结果和内部应力分布结果中的至少一种。
在另一个实施例中,迭代计算模块320还用于利用径向基函数将电磁场网 格节点的单位体积电磁损耗映射得到温度场网格节点的单位体积电磁损耗、并 利用温度场数学模型基于温度场网格节点的单位体积电磁损耗计算更新得到温 度场网格节点的温度;以及,利用径向基函数将温度场网格节点的温度映射得 到电磁场网格节点的温度、并利用电磁场数学模型基于电磁场网格节点的温度 计算更新得到电磁场网格节点的单位体积电磁损耗。
在另一个实施例中,迭代计算模块320还用于利用径向基函数将温度场网 格节点的温度映射得到应力场网格节点的温度、并利用应力场数学模型基于应 力场网格节点的温度计算更新得到应力场网格节点的应力张量;利用径向基函 数将应力场网格节点的应力张量映射得到温度场网格节点的应力张量、并利用 温度场数学模型基于温度场网格节点的应力张量计算更新得到温度场网格节点 的温度。
在另一个实施例中,请参考图4,初始计算模块310还内包括模型建立单元 311、矢量磁位计算单元312和电磁损耗计算单元313,其中:
模型建立单元311,用于建立待测电缆接头的几何模型并离散为若干个电磁 场网格节点。
矢量磁位计算单元312,用于基于待测电缆接头的工作参数计算得到各个电 磁场网格节点内的外加电流密度Js和材料电导率σ并代入对应的电磁场数学模 型中,计算得到矢量磁位A。
上述基于多物理场耦合的电缆接头状态评估装置中的各个模块可全部或部 分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于 计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中, 以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中,提供了一种计算机设备,该计算机设备可以是服务器, 其内部结构图可以如图5所示。该计算机设备包括通过系统总线连接的处理器、 存储器和网络接口。其中,该计算机设备的处理器用于提供计算和控制能力。 该计算机设备的存储器包括非易失性存储介质和内存储器。该非易失性存储介 质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中 的操作系统和计算机程序的运行提供环境。该计算机设备的数据库用于存储待 测电缆接头的工作参数、电磁场数学模型、温度场数学模型、应力场数学模型 等数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计 算机程序被处理器执行时以实现一种基于多物理场耦合的电缆接头状态评估方 法。
本领域技术人员可以理解,图5中示出的结构,仅仅是与本申请方案相关 的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定, 具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件, 或者具有不同的部件布置。
在一个实施例中,提供了一种计算机设备,包括存储器和处理器,存储器 中存储有计算机程序,该处理器执行计算机程序时实现以下步骤:
基于待测电缆接头的工作参数和电磁场数学模型计算得到电磁场网格节点 的单位体积电磁损耗;
由计算得到的电磁场网格节点的单位体积电磁损耗开始,利用径向基函数 在电磁场网格节点与温度场网格节点之间进行参量映射、以及利用径向基函数 在温度场网格节点和应力场网格节点之间进行参量映射,利用电磁场数学模型 基于电磁场网格节点接收到的参量计算更新得到电磁场网格节点的单位体积电 磁损耗,利用温度场数学模型基于温度场网格节点接收到的参量计算更新得到 温度场网格节点的温度,利用应力场数学模型基于应力场网格节点接收到的参 量计算更新得到应力场网格节点的应力张量;
在电磁场、温度场和应力场之间进行迭代的参量映射以及耦合计算直至达 到迭代收敛条件时,得到待测电缆接头的内部物理量分布结果并得到待测电缆 接头的状态评估结果,内部物理量分布结果包括内部电磁损耗分布结果、内部 温度分布结果和内部应力分布结果中的至少一种。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:利用径向基 函数将电磁场网格节点的单位体积电磁损耗映射得到温度场网格节点的单位体 积电磁损耗、并利用温度场数学模型基于温度场网格节点的单位体积电磁损耗 计算更新得到温度场网格节点的温度;以及,利用径向基函数将温度场网格节 点的温度映射得到电磁场网格节点的温度、并利用电磁场数学模型基于电磁场 网格节点的温度计算更新得到电磁场网格节点的单位体积电磁损耗。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:利用径向基 函数将温度场网格节点的温度映射得到应力场网格节点的温度、并利用应力场 数学模型基于应力场网格节点的温度计算更新得到应力场网格节点的应力张 量;利用径向基函数将应力场网格节点的应力张量映射得到温度场网格节点的 应力张量、并利用温度场数学模型基于温度场网格节点的应力张量计算更新得 到温度场网格节点的温度。
在一个实施例中,处理器执行计算机程序时还实现以下步骤:建立待测电 缆接头的几何模型并离散为若干个电磁场网格节点;基于待测电缆接头的工作 参数计算得到各个电磁场网格节点内的外加电流密度Js和材料电导率σ并代入 对应的电磁场数学模型中,计算得到矢量磁位A;按照计算得到 总电流密度J,并按照计算得到各个电磁场网格节点的单位体积电磁 损耗Qv。
在一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程 序,计算机程序被处理器执行时实现以下步骤:
基于待测电缆接头的工作参数和电磁场数学模型计算得到电磁场网格节点 的单位体积电磁损耗;
由计算得到的电磁场网格节点的单位体积电磁损耗开始,利用径向基函数 在电磁场网格节点与温度场网格节点之间进行参量映射、以及利用径向基函数 在温度场网格节点和应力场网格节点之间进行参量映射,利用电磁场数学模型 基于电磁场网格节点接收到的参量计算更新得到电磁场网格节点的单位体积电 磁损耗,利用温度场数学模型基于温度场网格节点接收到的参量计算更新得到 温度场网格节点的温度,利用应力场数学模型基于应力场网格节点接收到的参 量计算更新得到应力场网格节点的应力张量;
在电磁场、温度场和应力场之间进行迭代的参量映射以及耦合计算直至达 到迭代收敛条件时,得到待测电缆接头的内部物理量分布结果并得到待测电缆 接头的状态评估结果,内部物理量分布结果包括内部电磁损耗分布结果、内部 温度分布结果和内部应力分布结果中的至少一种。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:利用径向 基函数将电磁场网格节点的单位体积电磁损耗映射得到温度场网格节点的单位 体积电磁损耗、并利用温度场数学模型基于温度场网格节点的单位体积电磁损 耗计算更新得到温度场网格节点的温度;以及,利用径向基函数将温度场网格 节点的温度映射得到电磁场网格节点的温度、并利用电磁场数学模型基于电磁 场网格节点的温度计算更新得到电磁场网格节点的单位体积电磁损耗。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:利用径向 基函数将温度场网格节点的温度映射得到应力场网格节点的温度、并利用应力 场数学模型基于应力场网格节点的温度计算更新得到应力场网格节点的应力张 量;利用径向基函数将应力场网格节点的应力张量映射得到温度场网格节点的 应力张量、并利用温度场数学模型基于温度场网格节点的应力张量计算更新得 到温度场网格节点的温度。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:建立待测 电缆接头的几何模型并离散为若干个电磁场网格节点;基于待测电缆接头的工 作参数计算得到各个电磁场网格节点内的外加电流密度Js和材料电导率σ并代 入对应的电磁场数学模型中,计算得到矢量磁位A;按照计算得 到总电流密度J,并按照计算得到各个电磁场网格节点的单位体积电 磁损耗Qv。
在一个实施例中,提供了一种计算机程序产品,包括计算机程序,该计算 机程序被处理器执行时实现以下步骤:
基于待测电缆接头的工作参数和电磁场数学模型计算得到电磁场网格节点 的单位体积电磁损耗;
由计算得到的电磁场网格节点的单位体积电磁损耗开始,利用径向基函数 在电磁场网格节点与温度场网格节点之间进行参量映射、以及利用径向基函数 在温度场网格节点和应力场网格节点之间进行参量映射,利用电磁场数学模型 基于电磁场网格节点接收到的参量计算更新得到电磁场网格节点的单位体积电 磁损耗,利用温度场数学模型基于温度场网格节点接收到的参量计算更新得到 温度场网格节点的温度,利用应力场数学模型基于应力场网格节点接收到的参 量计算更新得到应力场网格节点的应力张量;
在电磁场、温度场和应力场之间进行迭代的参量映射以及耦合计算直至达 到迭代收敛条件时,得到待测电缆接头的内部物理量分布结果并得到待测电缆 接头的状态评估结果,内部物理量分布结果包括内部电磁损耗分布结果、内部 温度分布结果和内部应力分布结果中的至少一种。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:利用径向 基函数将电磁场网格节点的单位体积电磁损耗映射得到温度场网格节点的单位 体积电磁损耗、并利用温度场数学模型基于温度场网格节点的单位体积电磁损 耗计算更新得到温度场网格节点的温度;以及,利用径向基函数将温度场网格 节点的温度映射得到电磁场网格节点的温度、并利用电磁场数学模型基于电磁 场网格节点的温度计算更新得到电磁场网格节点的单位体积电磁损耗。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:利用径向 基函数将温度场网格节点的温度映射得到应力场网格节点的温度、并利用应力 场数学模型基于应力场网格节点的温度计算更新得到应力场网格节点的应力张 量;利用径向基函数将应力场网格节点的应力张量映射得到温度场网格节点的 应力张量、并利用温度场数学模型基于温度场网格节点的应力张量计算更新得 到温度场网格节点的温度。
在一个实施例中,计算机程序被处理器执行时还实现以下步骤:建立待测 电缆接头的几何模型并离散为若干个电磁场网格节点;基于待测电缆接头的工 作参数计算得到各个电磁场网格节点内的外加电流密度Js和材料电导率σ并代 入对应的电磁场数学模型中,计算得到矢量磁位A;按照计算得 到总电流密度J,并按照计算得到各个电磁场网格节点的单位体积电 磁损耗Qv。
需要说明的是,本申请所涉及的用户信息(包括但不限于用户设备信息、 用户个人信息等)和数据(包括但不限于用于分析的数据、存储的数据、展示 的数据等),均为经用户授权或者经过各方充分授权的信息和数据。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程, 是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于 一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述 各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、 数据库或其它介质的任何引用,均可包括非易失性和易失性存储器中的至少一 种。非易失性存储器可包括只读存储器(Read-Only Memory,ROM)、磁带、 软盘、闪存、光存储器、高密度嵌入式非易失性存储器、阻变存储器(ReRAM)、 磁变存储器(Magnetoresistive Random Access Memory,MRAM)、铁电存储器 (Ferroelectric Random Access Memory,FRAM)、相变存储器(Phase ChangeMemory,PCM)、石墨烯存储器等。易失性存储器可包括随机存取存储器(Random AccessMemory,RAM)或外部高速缓冲存储器等。作为说明而非局 限,RAM可以是多种形式,比如静态随机存取存储器(Static Random Access Memory,SRAM)或动态随机存取存储器(Dynamic Random Access Memory, DRAM)等。本申请所提供的各实施例中所涉及的数据库可包括关系型数据库 和非关系型数据库中至少一种。非关系型数据库可包括基于区块链的分布式数 据库等,不限于此。本申请所提供的各实施例中所涉及的处理器可为通用处理器、中央处理器、图形处理器、数字信号处理器、可编程逻辑器、基于量子计 算的数据处理逻辑器等,不限于此。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述 实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特 征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细, 但并不能因此而理解为对本申请专利范围的限制。应当指出的是,对于本领域 的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和 改进,这些都属于本申请的保护范围。因此,本申请的保护范围应以所附权利 要求为准。
Claims (10)
1.一种基于多物理场耦合的电缆接头状态评估方法,其特征在于,所述方法包括:
基于待测电缆接头的工作参数和电磁场数学模型计算得到电磁场网格节点的单位体积电磁损耗;
由计算得到的所述电磁场网格节点的单位体积电磁损耗开始,利用径向基函数在电磁场网格节点与温度场网格节点之间进行参量映射、以及利用径向基函数在温度场网格节点和应力场网格节点之间进行参量映射,利用电磁场数学模型基于电磁场网格节点接收到的参量计算更新得到电磁场网格节点的单位体积电磁损耗,利用温度场数学模型基于温度场网格节点接收到的参量计算更新得到温度场网格节点的温度,利用应力场数学模型基于应力场网格节点接收到的参量计算更新得到应力场网格节点的应力张量;
在电磁场、温度场和应力场之间进行迭代的参量映射以及耦合计算直至达到迭代收敛条件时,得到所述待测电缆接头的内部物理量分布结果并得到所述待测电缆接头的状态评估结果,所述内部物理量分布结果包括内部电磁损耗分布结果、内部温度分布结果和内部应力分布结果中的至少一种。
2.根据权利要求1所述的方法,其特征在于,所述利用径向基函数在电磁场网格节点与温度场网格节点之间进行参量映射,包括:
利用径向基函数将电磁场网格节点的单位体积电磁损耗映射得到温度场网格节点的单位体积电磁损耗、并利用温度场数学模型基于温度场网格节点的单位体积电磁损耗计算更新得到温度场网格节点的温度;
以及,利用径向基函数将温度场网格节点的温度映射得到电磁场网格节点的温度、并利用电磁场数学模型基于电磁场网格节点的温度计算更新得到电磁场网格节点的单位体积电磁损耗。
3.根据权利要求1所述的方法,其特征在于,所述利用径向基函数在温度场网格节点和应力场网格节点之间进行参量映射,包括:
利用径向基函数将温度场网格节点的温度映射得到应力场网格节点的温度、并利用应力场数学模型基于应力场网格节点的温度计算更新得到应力场网格节点的应力张量;
利用径向基函数将应力场网格节点的应力张量映射得到温度场网格节点的应力张量、并利用温度场数学模型基于温度场网格节点的应力张量计算更新得到温度场网格节点的温度。
4.根据权利要求1所述的方法,其特征在于,利用径向基函数将第一物理场网格节点中的参量映射得到第二物理场网格节点中的参量表示为:
AS=HSfAf;
其中,表示n个第二物理场网格节点中的参量,对于任意的1≤i≤n,表示任意第i个第二物理场网格节点si中的参量;表示n个第一物理场网格节点中的参量,对于任意的1≤j≤n,表示任意第j个第一物理场网格节点fj中的参量;
表示任意第j个第一物理场网格节点fj的坐标,表示任意第i个第二物理场网格节点si中的参量;对于任意的1≤k≤n,表示任意的第j个第一物理场网格节点fj和第k个第一物理场网格节点fk之间的径向基函数关系;表示任意第j个第一物理场网格节点fj与第i个第二物理场网格节点si之间的径向基函数关系;
其中,所述第一物理场网格节点和所述第二物理场网格节点中一个是电磁场网格节点、另一个是温度场网格节点,或者,所述第一物理场网格节点和所述第二物理场网格节点中一个是温度场网格节点、另一个是应力场网格节点。
5.根据权利要求1-4任一所述的方法,其特征在于,所述电磁场数学模型包括有外加激励源存在的电流区域的第一控制方程和非电流区域的第二控制方程,且第一控制方程为第二控制方程为其中,A表示矢量磁位,μ表示材料磁导率,σ表示材料电导率,Js表示外加电流密度,ω为角频率;
所述基于待测电缆接头的工作参数和电磁场数学模型计算得到电磁场网格节点的单位体积电磁损耗,包括:
建立所述待测电缆接头的几何模型并离散为若干个电磁场网格节点;
基于所述待测电缆接头的工作参数计算得到各个电磁场网格节点内的外加电流密度Js和材料电导率σ并代入对应的电磁场数学模型中,计算得到矢量磁位A;
7.一种基于多物理场耦合的电缆接头状态评估装置,其特征在于,所述装置包括:
初始计算模块,用于基于待测电缆接头的工作参数和电磁场数学模型计算得到电磁场网格节点的单位体积电磁损耗;
迭代计算模块,用于由计算得到的所述电磁场网格节点的单位体积电磁损耗开始,利用径向基函数在电磁场网格节点与温度场网格节点之间进行参量映射、以及利用径向基函数在温度场网格节点和应力场网格节点之间进行参量映射,利用电磁场数学模型基于电磁场网格节点接收到的参量计算更新得到电磁场网格节点的单位体积电磁损耗,利用温度场数学模型基于温度场网格节点接收到的参量计算更新得到温度场网格节点的温度,利用应力场数学模型基于应力场网格节点接收到的参量计算更新得到应力场网格节点的应力张量;
结果获取模块,用于在电磁场、温度场和应力场之间进行迭代的参量映射以及耦合计算直至达到迭代收敛条件时,得到所述待测电缆接头的内部物理量分布结果并得到所述待测电缆接头的状态评估结果,所述内部物理量分布结果包括内部电磁损耗分布结果、内部温度分布结果和内部应力分布结果中的至少一种。
8.一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至6中任一项所述的方法的步骤。
9.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至6中任一项所述的方法的步骤。
10.一种计算机程序产品,包括计算机程序,其特征在于,该计算机程序被处理器执行时实现权利要求1至6中任一项所述的方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111527853.7A CN114462264A (zh) | 2021-12-14 | 2021-12-14 | 一种基于多物理场耦合的电缆接头状态评估方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111527853.7A CN114462264A (zh) | 2021-12-14 | 2021-12-14 | 一种基于多物理场耦合的电缆接头状态评估方法及装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114462264A true CN114462264A (zh) | 2022-05-10 |
Family
ID=81405591
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111527853.7A Pending CN114462264A (zh) | 2021-12-14 | 2021-12-14 | 一种基于多物理场耦合的电缆接头状态评估方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114462264A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116052820A (zh) * | 2023-03-22 | 2023-05-02 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种材料热性能评估方法、装置、设备及可读存储介质 |
-
2021
- 2021-12-14 CN CN202111527853.7A patent/CN114462264A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116052820A (zh) * | 2023-03-22 | 2023-05-02 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种材料热性能评估方法、装置、设备及可读存储介质 |
CN116052820B (zh) * | 2023-03-22 | 2023-06-02 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种材料热性能评估方法、装置、设备及可读存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112327166B (zh) | 基于分数阶平方根无迹卡尔曼滤波的锂电池soc估计方法 | |
CN110023966B (zh) | 使用量子计算仿真材料 | |
CN109755683A (zh) | 一种基于压缩感知理论的电池包内部温度实时监测方法 | |
Bargiela | On-line monitoring of water distribution networks | |
CN114462264A (zh) | 一种基于多物理场耦合的电缆接头状态评估方法及装置 | |
Zhang et al. | Battery internal temperature estimation via a semilinear thermal PDE model | |
Rahman et al. | Investigating effects of weather parameter uncertainty on transmission line power handling capabilities using affine arithmetic | |
Gao et al. | Physics embedded graph convolution neural network for power flow calculation considering uncertain injections and topology | |
CN111507591B (zh) | 电力系统状态确定方法、装置、计算机介质及存储介质 | |
Pawlowski et al. | Flow-based density of states for complex actions | |
CN114168574A (zh) | 一种面向工业负荷的数据缺失处理方法及装置 | |
CN116400247B (zh) | 电池软短路故障的确定方法及装置 | |
Al-Othman | A fuzzy state estimator based on uncertain measurements | |
Stemler et al. | A guide to using shadowing filters for forecasting and state estimation | |
Xie et al. | Grounding grid corrosion detection based on mini-batch gradient descent and greedy method | |
CN117114161A (zh) | 一种基于元学习的输电线路风偏闪络风险的预测方法 | |
CN112731150B (zh) | 电压暂降状态估计方法、装置、计算机设备和存储介质 | |
CN115859783A (zh) | 用于核反应堆仿真测试的多保真度网络构建方法和装置 | |
CN114036806A (zh) | 基于热导率各向异性介质的三维地温场数值模拟方法 | |
CN114036805A (zh) | 三维稳态热传导地温场的正演方法、装置、设备及介质 | |
CN113688465A (zh) | 一种基于载荷与状态结合的飞行器结构强度数字孪生方法 | |
CN114487552A (zh) | 基于微型电场传感芯片的电压测量方法、装置和设备 | |
Cai et al. | Low rank matrix completion for recovering missing load data in power system | |
Girault et al. | Parametric low-order models in transient heat diffusion by MIM. Estimation of thermal conductivity in a 2D slab | |
CN112182739A (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 |