CN109408927A - 一种基于黑盒传输线模型的二维静磁场并行有限元加速方法 - Google Patents
一种基于黑盒传输线模型的二维静磁场并行有限元加速方法 Download PDFInfo
- Publication number
- CN109408927A CN109408927A CN201811192762.0A CN201811192762A CN109408927A CN 109408927 A CN109408927 A CN 109408927A CN 201811192762 A CN201811192762 A CN 201811192762A CN 109408927 A CN109408927 A CN 109408927A
- Authority
- CN
- China
- Prior art keywords
- linear
- finite element
- transmission line
- circuit
- iteration
- 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
Links
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]
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)
- Complex Calculations (AREA)
Abstract
本发明公开了一种基于黑盒传输线模型的二维非线性静磁场模型的有限元求解方法,所述方法包括如下步骤:一、确定待求解的变量以及求解域;二、建立rz直角坐标系;三、列出二维轴对称非线性静磁场中的控制方程和边界条件式并组成一微分方程组;四、对求解域进行分网;五、计算有限元的系数矩阵;六、建立等效的黑盒电路模型;七、在黑盒电路模型与线性电路之间***一段传输线线段;八、进行传输线法迭代;九、重复进行迭代,直到迭代结果收敛到固定误差,结束求解。本发明能够对所求问题并行地进行求解计算,并且能够用于复杂的有限元分网模型求解当中,从而解决牛顿迭代法求解有限元非线性问题时带来的求解时间长、效率低的问题。
Description
技术领域
本发明属于电磁场数值计算技术领域,涉及一种二维非线性静磁 场模型的有限元求解方法,具体涉及一种将二维非线性静磁场的有限 元问题转化为等价的黑盒电路问题,并借助传输线迭代法实现对二维 非线性静磁场的有限元并行加速求解的方法。
背景技术
有限元法是工业领域中应用最广泛的数值计算方法,被诸多商用 仿真软件采用。然而,随着求解模型的日益复杂化以及分网单元数目 的不断增多,以传统的牛顿迭代法为核心的非线性有限元求解方法面 临着求解耗时严重的问题,这直接关系到产品的仿真研发的速度和效 率。
有限元问题的求解的核心在于求解线性方程组,而对于非线性问 题来说,传统的牛顿迭代法每一步都要利用新的迭代结果重新生成有 限元模型的全局矩阵,随着模型分网的不断增大,全局矩阵的维度不 断变大,每一步矩阵的LU分解等消耗的时间会相应的增大,总体的 求解时间可能随着分网的变密而成几何式增大。使用传输线迭代法能 够避免重复的LU分解过程,求解时间。然而,传统的传输线迭代法 只能用于一阶三角形单元的有限元模型当中,无法用于复杂的分网单 元模型。因此,需要研究一种新的迭代方法,来提高传输线法的适用 范围和应对复杂模型的能力,从而解决牛顿迭代法求解有限元非线性 问题时带来的求解时间长、效率低的问题。
发明内容
本发明的目的是提供一种基于黑盒传输线模型的二维非线性静 磁场模型的有限元求解方法,该方法能够对所求问题并行地进行求解 计算,并且能够用于复杂的有限元分网模型求解当中,包括三角形单 元和四边形单元等,从而解决牛顿迭代法求解有限元非线性问题时带 来的求解时间长、效率低的问题。
本发明的目的是通过以下技术方案实现的:
一种基于黑盒传输线模型的二维非线性静磁场模型的有限元求 解方法,包括如下步骤:
一、确定待求解的变量以及求解域,其中:所述待求解的变量为 一二维轴对称非线性静磁场的磁势A,所述求解域为二维轴对称非线 性静磁场所在的区域;
二、建立一个rz直角坐标系;
三、列出二维轴对称非线性静磁场中的控制方程和边界条件式并 组成一微分方程组,其中:
所述控制方程为:
式中,J为电流密度变量,μ′=rμ,μ为三角单元的磁导率,A′=rA, A为磁势,r为横坐标,z为纵坐标;
所述边界条件式为:
Γ1:A′=0;
式中,m为求解边界上的法向量;
四、采用分网程序对求解域进行分网,将求解域离散为三角形单 元或者四边形单元,在每一个有限元单元中,变量式中, n为有限元单元的节点数目,Nj为单元当中的形函数,为对应节点 的变量值的大小;
五、按照下式计算有限元的系数矩阵:
其中,每一项的计算表达式如下:
六、将步骤五当中的矩阵Ke看作电路的导纳矩阵,be看作电流源 向量,建立等效的黑盒电路模型;
七、在黑盒电路模型与线性电路之间的连接导线任意位置***一 段传输线线段;
八、进行传输线法迭代:
将电压信号入射到线性网络中,求解该电路,即方程组: (Ylinear+YTL)A′=b+2ViYTL,在整个迭代过程当中,矩阵(Ylinear+YTL)保持不 变,在迭代的第1步执行一次矩阵的LU分解操作,后续步骤无需再 次执行,求解完成后,计算反射电压Vr=A′-Vi;
将电压信号入射到非线性有限元单元中,求解该电路,即非线性 方程组:采用牛顿迭代法进行求解,第k步的迭 代公式如下:其中:
每个非线性单元内的方程组求解都单独放到独立的计算核心中 求解,实现并行计算,求解完成后,再次计算反射电压
式中,Ylinear为入射到线性网络内的线性电路的导纳矩阵;YTL为入 射到线性网络内的传输线电路的导纳矩阵;A′为待求节点电压;Vr为 反射到非线性电路的反射电压;Vi为入射到线性电路的入射电压;Ke为非线性单元的系数矩阵;为反射回非线性单元内的传输线的导纳 矩阵;为反射回非线性单元内的反射电压;Je为雅克比矩阵;为 第k次迭代时的雅克比矩阵;f为计算当中的中间变量;Ae′为单元内 每个节点的磁势值。
九、重复进行步骤八的迭代,直到迭代结果收敛到固定误差,结 束求解。
相比于现有技术,本发明具有如下优点:
1、采用有限元方法对静磁场进行离散建模,将有限元中的每个 单元的单元系数矩阵等效为一个密封的黑盒电路模型,从而实现将有 限元方程组等效为一个非线性的电路网络问题。
2、在非线性单元的黑盒模型与线性单元的网络之间***传输线 线段,利用传输线的迭代规律,实现传输线迭代,从而实现静磁场问 题的求解。
3、在传输线迭代过程中,每个非线性黑盒模型都是相互独立的, 它的求解计算可以采用单独的计算资源(CPU或者GPU核心)进行 计算,实现并行加速。
4、采用传输线法和黑盒模型能够进行有限元的迭代求解,能够 处理多种复杂的有限元模型,例如三角形单元和四边形单元。
5、在迭代求解过程中,全局矩阵Y能够保持不变,在矩阵求解 过程当中,采用LU分解法,只需要在计算的第一步进行LU分解, 由于LU分解一般占用矩阵求解的95%左右的时间,使用这种方法将 大幅减小每一步的迭代时间。
6、使用传输线法能够将非线性单元与线性求解区域进行隔离, 实现并行加速效果。
附图说明
图1为黑盒电路模型;
图2为***传输线后的黑盒电路网络;
图3为电压入射到线性网络后的等效电路图;
图4为电压反射回非线性单元后的等效电路图;
图5为迭代求解流程图;
图6为某型号接触器的结构图;
图7为某型号接触器的剖面图;
图8为接触器简化模型;
图9为接触器计算模型;
图10为分网结果,(a)三角形单元,(b)四边形单元;
图11为采用传输线法求解得到不同状态下的磁场分布情况;
图12为单步计算时间加速效果。
具体实施方式
下面结合附图对本发明的技术方案作进一步的说明,但并不局限 于此,凡是对本发明技术方案进行修改或者等同替换,而不脱离本发 明技术方案的范围,均应涵盖在本发明的保护范围中。
具体实施方式一:本实施方式提供了一种基于黑盒传输线模型的 二维非线性静磁场模型的有限元求解方法,所述方法包括如下步骤:
一、确定待求解的变量以及求解域,待求解的变量为一二维轴对 称非线性静磁场的磁势A,二维轴对称非线性静磁场由通电线圈中的 电流产生,通电线圈周围的各元件均为铁磁材料,求解域为二维轴对 称非线性静磁场所在的区域。
二、建立一个rz直角坐标系,r为横轴,z为纵轴。
三、列出二维轴对称非线性静磁场中的控制方程和边界条件式并 组成一微分方程组,其控制方程为:其中, J为电流密度变量,μ′=rμ,μ为三角单元的磁导率,A′=rA,A为磁 势,r为横坐标,z为纵坐标,边界条件为:Γ1:A′=0; m为求解边界上的法向量。
四、采用分网程序对求解区域进行分网,将求解域离散为三角形 单元或者四边形单元。在每一个有限元单元中,变量其 中,n为有限元单元的节点数目,Nj为单元当中的形函数,为对应 节点的变量值的大小。
五、计算有限元的系数矩阵。根据伽辽金法,对控制方程进行加 权积分,得到方程:
写成矩阵形式KeAe′=be,即为:
其中,每一项的计算表达式如下:
六、将步骤五当中的矩阵Ke看作电路的导纳矩阵,be看作电流源 向量,建立等效的黑盒电路模型,如图1所示。
七、黑盒模型与线性电路之间是通过导线相连接的,在黑盒模型 与线性电路之间的这些导线的任意位置***一段传输线线段,如图2 所示。
八、进行传输线法迭代。
电压信号入射到线性网络中,如图3所示,求解该电路,即方程 组:(Ylinear+YTL)A′=b+2ViYTL,在整个迭代过程当中,矩阵(Ylinear+YTL)保 持不变,在迭代的第1步执行一次矩阵的LU分解操作,后续步骤无 需再次执行,求解完成后,计算反射电压Vr=A′-Vi。
电压信号入射到非线性有限元单元中,如图4所示,求解该电路, 即非线性方程组:该矩阵采用牛顿迭代法进行求 解,第k步的迭代公式如下:其中:
由于每个非线性单元相互之间都是隔离的,每个单元内的方程组求解 都单独放到独立的计算核心中求解,实现并行计算。求解完成后,再 次计算反射电压
式中,Ylinear为图3当中线性电路的导纳矩阵;YTL为图3虚线框内 传输线电路的导纳矩阵;A′为待求节点电压;Vr为反射到非线性电路 的反射电压;Vi为入射到线性电路的入射电压;Ke为非线性单元的系 数矩阵;为图4中传输线的导纳矩阵;为反射回单元内的反射电 压;Je为雅克比矩阵;为第k次迭代时的雅克比矩阵。
不断地进行迭代,直到迭代结果收敛到固定误差,结束求解。求 解流程图如图5。
具体实施方式二:本实施方式提供了一种基于黑盒传输线模型的 二维非线性静磁场模型的有限元求解方法,所述方法具体实施步骤如 下:
一、确定待求解的变量以及求解域,待求解的变量为一二维轴对 称非线性静磁场的磁势A,二维轴对称非线性静磁场由通电线圈中的 电流产生,通电线圈周围的各元件均为铁磁材料,求解域为二维轴对 称非线性静磁场所在的区域。本实施方式针对某型号轴对称结构接触 器进行计算,该模型的外观图如图6所示,剖面图如图7所示。对其 进行简化,并建模,简化模型如图8所示,计算模型如图9所示。
二、建立一个rz直角坐标系,r为横轴,z为纵轴。
三、列出二维轴对称非线性静磁场中的控制方程和边界条件式并 组成一微分方程组,其控制方程为:其中, J为电流密度变量,μ′=rμ,μ为三角单元的磁导率,A′=rA,A为磁 势,r为横坐标,z为纵坐标,边界条件为:Γ1:A′=0; m为求解边界上的法向量。
四、采用分网程序对求解区域进行分网,将求解域离散为三角形 单元或者四边形单元,分网结果如图10所示。在每一个有限元单元 中,变量其中,n为有限元单元的节点数目,Nj为单元 当中的形函数,为对应节点的变量值的大小。
对于三角形单元,Δ为三角形单元 的面积,
对于四边形单元,
五、计算有限元的系数矩阵。根据伽辽金法,对控制方程进行加 权积分,得到方程:
写成矩阵形式KeAe′=be,即为:
其中,每一项的计算表达式如下:
六、将步骤五当中的矩阵Ke看作电路的导纳矩阵,be看作电流源 向量,建立等效的黑盒电路模型,如图1所示。
七、在黑盒模型与线性电路之间***一段传输线线段,如图2所 示。
八、进行传输线法迭代。
电压信号入射到线性网络中,如图3所示,求解该电路,即方程 组:(Ylinear+YTL)A′=b+2ViYTL,然后,计算反射电压Vr=A′-Vi。
电压信号入射到非线性有限元单元中,如图4所示,求解该电路, 即非线性方程组:该矩阵采用牛顿迭代法进行求 解,第k步的迭代公式如下:其中:
求解完成后,再次计算反射电压
不断地进行迭代,直到迭代结果收敛到固定误差,结束求解。求 解流程图如图5。
采用传输线法求解得到不同状态下的磁场分布情况,结果如图11 所示。
表1不同的求解模型
case 1-5为三角形单元分网;case 6-10为四边形单元分网。
表2 N-R迭代法和本发明BB-TLM迭代法的误差对比
No.1-5为表1case 5模型当中任取五个点的误差对比;No.6-10为表1 case 10模型任取五个点的误差对比。
对表1当中case10模型分别使用1、4、8、16、20个CPU核心 进行加速,对比牛顿迭代法与传输线法的单步计算时间,由图12可 见,传输线法的单步计算时间比牛顿迭代法快几倍。
对表1当中case10模型分别使用1、4、8、16、20个CPU核心 进行加速,对比牛顿迭代法与传输线法的总计算时间,当CPU核心 增大时,传输线法的加速效果比牛顿迭代法更加明显。
表3针对表1中的模型分别采用牛顿迭代法和传输线法求解(CPU 核心为20)
由表3可知,当分网模型增大时,传输线法的计算时间比牛顿迭 代法快。
Claims (1)
1.一种基于黑盒传输线模型的二维非线性静磁场模型的有限元求解方法,其特征在于所述方法包括如下步骤:
一、确定待求解的变量以及求解域,其中:所述待求解的变量为一二维轴对称非线性静磁场的磁势A,所述求解域为二维轴对称非线性静磁场所在的区域;
二、建立一个rz直角坐标系;
三、列出二维轴对称非线性静磁场中的控制方程和边界条件式并组成一微分方程组,其中:
所述控制方程为:
式中,J为电流密度变量,μ′=rμ,μ为三角单元的磁导率,A′=rA,A为磁势,r为横坐标,z为纵坐标;
所述边界条件式为:
Γ1:A′=0;
式中,m为求解边界上的法向量;
四、采用分网程序对求解域进行分网,将求解域离散为三角形单元或者四边形单元,在每一个有限元单元中,变量式中,n为有限元单元的节点数目,Nj为单元当中的形函数,为对应节点的变量值的大小;
五、按照下式计算有限元的系数矩阵:
其中,每一项的计算表达式如下:
六、将步骤五当中的矩阵Ke看作电路的导纳矩阵,be看作电流源向量,建立等效的黑盒电路模型;
七、在黑盒电路模型与线性电路之间的连接导线任意位置***一段传输线线段;
八、进行传输线法迭代:
将电压信号入射到线性网络中,求解该电路,即方程组:(Ylinear+YTL)A′=b+2ViYTL,在整个迭代过程当中,矩阵(Ylinear+YTL)保持不变,在迭代的第1步执行一次矩阵的LU分解操作,后续步骤无需再次执行,求解完成后,计算反射电压Vr=A′-Vi;
将电压信号入射到非线性有限元单元中,求解该电路,即非线性方程组:采用牛顿迭代法进行求解,第k步的迭代公式如下:每个非线性单元内的方程组求解都单独放到独立的计算核心中求解,实现并行计算,求解完成后,再次计算反射电压
式中,Ylinear为入射到线性网络内的线性电路的导纳矩阵;YTL为入射到线性网络内的传输线电路的导纳矩阵;A′为待求节点电压;Vr为反射到非线性电路的反射电压;Vi为入射到线性电路的入射电压;Ke为非线性单元的系数矩阵;为反射回非线性单元内的传输线的导纳矩阵;为反射回非线性单元内的反射电压;为第k次迭代时的雅克比矩阵;
九、重复进行步骤八的迭代,直到迭代结果收敛到固定误差,结束求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811192762.0A CN109408927B (zh) | 2018-10-13 | 2018-10-13 | 一种基于黑盒传输线模型的二维静磁场并行有限元加速方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811192762.0A CN109408927B (zh) | 2018-10-13 | 2018-10-13 | 一种基于黑盒传输线模型的二维静磁场并行有限元加速方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109408927A true CN109408927A (zh) | 2019-03-01 |
CN109408927B CN109408927B (zh) | 2022-01-18 |
Family
ID=65467151
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811192762.0A Active CN109408927B (zh) | 2018-10-13 | 2018-10-13 | 一种基于黑盒传输线模型的二维静磁场并行有限元加速方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109408927B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112016232A (zh) * | 2020-08-31 | 2020-12-01 | 中国原子能科学研究院 | 一种撕裂有限元过程处理方法及*** |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050228601A1 (en) * | 2003-06-05 | 2005-10-13 | Fujitsu Limited | Device and method for calculating electro-magnetic field intensity by cooperation of circuit analyses and electro-magnetic wave analyses |
CN105631105A (zh) * | 2015-12-24 | 2016-06-01 | 北京航空航天大学 | 一种含外加场耦合的线缆等效电路建模方法 |
CN106649939A (zh) * | 2016-09-28 | 2017-05-10 | 哈尔滨工业大学 | 基于传输线迭代的2d轴对称非线性静磁场模型的求解方法 |
CN106991263A (zh) * | 2017-05-25 | 2017-07-28 | 清华大学 | 一种中低压配电变压器绕组雷电宽频模型的建模方法 |
CN107609274A (zh) * | 2017-09-14 | 2018-01-19 | 哈尔滨工业大学 | 基于传输线与级别调度法的二维静磁场并行有限元方法 |
US20180246154A1 (en) * | 2015-09-30 | 2018-08-30 | Osaka University | Method and program for calculating potential, current, and peripheral electromagnetic field in electric circuit |
-
2018
- 2018-10-13 CN CN201811192762.0A patent/CN109408927B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050228601A1 (en) * | 2003-06-05 | 2005-10-13 | Fujitsu Limited | Device and method for calculating electro-magnetic field intensity by cooperation of circuit analyses and electro-magnetic wave analyses |
US20180246154A1 (en) * | 2015-09-30 | 2018-08-30 | Osaka University | Method and program for calculating potential, current, and peripheral electromagnetic field in electric circuit |
CN105631105A (zh) * | 2015-12-24 | 2016-06-01 | 北京航空航天大学 | 一种含外加场耦合的线缆等效电路建模方法 |
CN106649939A (zh) * | 2016-09-28 | 2017-05-10 | 哈尔滨工业大学 | 基于传输线迭代的2d轴对称非线性静磁场模型的求解方法 |
CN106991263A (zh) * | 2017-05-25 | 2017-07-28 | 清华大学 | 一种中低压配电变压器绕组雷电宽频模型的建模方法 |
CN107609274A (zh) * | 2017-09-14 | 2018-01-19 | 哈尔滨工业大学 | 基于传输线与级别调度法的二维静磁场并行有限元方法 |
Non-Patent Citations (4)
Title |
---|
LI BING等: ""Influence of electromagnetic pulse on electromagnetic relay by using equivalent circuit model"", 《26TH INTERNATIONAL CONFERENCE ON ELECTRICAL CONTACTS (ICEC 2012)》 * |
周建美: ""各向异性地层中可控源电磁法一维全参数反演及三维有限体积正演算法研究"", 《中国博士学位论文全文数据库 基础科学辑》 * |
张旭锋: ""传输线理论及电磁兼容计算的半解析方法研究"", 《中国博士学位论文全文数据库 信息科技辑》 * |
张雷: ""基于巨磁阻抗效应的弱磁传感器技术研究"", 《中国博士学位论文全文数据库 信息科技辑》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112016232A (zh) * | 2020-08-31 | 2020-12-01 | 中国原子能科学研究院 | 一种撕裂有限元过程处理方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN109408927B (zh) | 2022-01-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103902785B (zh) | 一种基于多元不确定性的结构有限元模型修正方法 | |
Shivanian | Meshless local Petrov–Galerkin (MLPG) method for three-dimensional nonlinear wave equations via moving least squares approximation | |
Wei et al. | Dynamical analysis of the generalized Sprott C system with only two stable equilibria | |
CN106777827A (zh) | 一种机电‑电磁混合仿真方法及*** | |
CN110275733A (zh) | 基于有限体积法求解声子玻尔兹曼方程的gpu并行加速方法 | |
CN105005675A (zh) | 基于多目标遗传算法的复合绝缘子电场优化方法 | |
Liu et al. | Bifurcation of a predator–prey model with disease in the prey | |
CN108921288A (zh) | 神经网络激活处理装置和基于该装置的神经网络处理器 | |
CN110502785A (zh) | 一种三维时域计算波导s参数的电磁数值方法 | |
CN107038302A (zh) | 一种含柱状缺陷的涡流无损检测磁场的半解析计算方法 | |
CN109408927A (zh) | 一种基于黑盒传输线模型的二维静磁场并行有限元加速方法 | |
CN106777729A (zh) | 一种基于fpga的算法库仿真验证平台实现方法 | |
Chandar et al. | A GPU-based incompressible Navier–Stokes solver on moving overset grids | |
CN106649939B (zh) | 基于传输线迭代的2d轴对称非线性静磁场模型的求解方法 | |
Nicolas et al. | Improved adaptive mesh refinement for conformal hexahedral meshes | |
De Su et al. | Numerical solution for a kind of nonlinear telegraph equations using radial basis functions | |
Yin et al. | Multi-step Prediction Algorithm of Traffic Flow Chaotic Time Series Based on Volterra Neural Network. | |
CN103942376A (zh) | 一种应用于实时头发运动仿真的碰撞修正方法 | |
CN106651979A (zh) | 基于投影的空间曲面四边形网格剖分方法 | |
Moirangthem et al. | GIS integrated automation of a near real-time power-flow service for electrical grids | |
CN106096176A (zh) | 屏蔽线缆的软辫线的天线模型的建模方法及装置 | |
Monsegny et al. | Fast marching method in seismic ray tracing on parallel GPU devices | |
CN107609274B (zh) | 基于传输线与级别调度法的二维静磁场并行有限元方法 | |
CN112417794B (zh) | 一种散射参数计算方法 | |
CN105589678A (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 |