CN112683965A - 一种复合材料介电常数计算方法及系统 - Google Patents

一种复合材料介电常数计算方法及系统 Download PDF

Info

Publication number
CN112683965A
CN112683965A CN202011405208.3A CN202011405208A CN112683965A CN 112683965 A CN112683965 A CN 112683965A CN 202011405208 A CN202011405208 A CN 202011405208A CN 112683965 A CN112683965 A CN 112683965A
Authority
CN
China
Prior art keywords
equivalent
dielectric constant
unit
circuit
filler
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
CN202011405208.3A
Other languages
English (en)
Other versions
CN112683965B (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.)
North China Electric Power University
Original Assignee
North China Electric Power University
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 North China Electric Power University filed Critical North China Electric Power University
Priority to CN202011405208.3A priority Critical patent/CN112683965B/zh
Publication of CN112683965A publication Critical patent/CN112683965A/zh
Application granted granted Critical
Publication of CN112683965B publication Critical patent/CN112683965B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Measurement Of Resistance Or Impedance (AREA)

Abstract

本发明公开了一种复合材料介电常数计算方法及系统。将复合材料划分为独立的立方体单元结构,对立方体结构进行了区域划分,对界面极化不规则区域进行了电路等效,考虑了不同材料间的界面极化,不但能计算高频段的介电谱,也能计算复合材料的中、低频段的介电谱;并且利用第二关系式,能够根据基体和填充物的介电谱线计算复合材料的介电常数频谱,亦可根据复合材料及其中一种材料的介电谱线反推另外一种材料的介电常数频谱,极大降低了工作量,节约了测试时间。

Description

一种复合材料介电常数计算方法及系统
技术领域
本发明涉及材料性能评价技术领域,特别是涉及一种复合材料介电常数计算方法及系统。
背景技术
随着现代科学技术的发展,在电力系统中出现了众多新型复合电介质材料,包括纳米粒子填充材料、发泡材料、纤维薄膜等,在其进行设计的过程中,组成两相复合介质的基体与填充物的介电常数具有很重要的参考价值,但由于基体材料孔隙的影响无法将基体与填充物分离分别进行介电响应测试来获得介电常数,所以需要通过整体的介电频谱反推组成部分的介电频谱。同时,介电响应测试时间很长,为了迅速对合成后的材料性能进行初步了解,需要通过复合材料的组成部分的介电频谱反推整体的介电频谱以提高研发效率,节约测试时间。
在电力系统中为保证电气设备的安全稳定可靠运行,则需要对这些材料进行绝缘性能监测,其中材料的介电性能是一项关键指标。在材料的绝缘状态诊断中发现,低频下的介电参数更能反应材料的缺陷,但是现有的预测多相复合材料介电特性模型主要表征材料的高频介电性能,并没有考虑不同材料的界面极化,仅能与实测谱在高频段内相吻合,而在中、低频段内则存在较大差异,无法分析材料在中、低频下的介电性能。并且不同结构的电介质材料需要分别进行介电响应测试,不能进行理论分析,操作繁琐、耗费时间。
发明内容
本发明的目的是提供一种复合材料介电常数计算方法及系统,考虑了不同材料间的界面极化,能够高效计算高、中、低频段的介电谱,极大降低了工作量,节约了测试时间。
为实现上述目的,本发明提供了如下方案:
一种复合材料介电常数计算方法,包括:
对复合材料进行单元划分,得到多个相同的立方体单元;每一个所述立方体单元均包括一个位于立方体单元中心位置处的填充物和设置在所述填充物外表面的基体;
将所述立方体单元划分为第一区域和第二区域,根据所述第一区域和所述第二区域生成等效电路;其中,所述第一区域包括第一基体,所述第二区域包括上基体、填充物和下基体;所述等效电路包括并联的第一等效电路单元和第二等效电路单元;所述第一等效电路单元根据所述第一基体确定,所述第二等效电路单元根据所述上基体、所述填充物和所述下基体确定;
根据所述等效电路的参数,建立基体介电常数和填充物介电常数与所述立方体单元的等效阻抗之间的关系式,得到第一关系式;
根据所述第一关系式,建立基体介电常数、填充物介电常数和复合材料介电常数之间的关系式,得到第二关系式;
获取基体、填充物、复合材料三种材料中的任两种材料的介电常数频谱,根据所述第二关系式采用最小二乘法确定另一种材料的介电常数频谱。
可选的,所述根据所述第一区域和所述第二区域生成等效电路,具体包括:
根据所述第一区域生成第一等效电路单元;所述第一等效电路单元包括多个相互并联的第一RC并联电路,所述第一RC并联电路的参数根据所述第一基体确定;
在垂直于电场方向,将所述第二区域划分为多个第二子区域;
将所述第二子区域等效为第二等效电路子单元;所述第二等效电路子单元为相互串联的第二RC并联电路、第三RC并联电路和第四RC并联电路,所述第二RC并联电路的参数根据所述上基体确定,所述第三RC并联电路的参数根据所述填充物确定,所述第四RC并联电路的参数根据所述下基体确定;
将多个所述第二等效电路子单元并联,生成第二等效电路单元;
将所述第一等效电路单元和所述第二等效电路单元并联,生成等效电路。
可选的,所述根据所述等效电路的参数,建立基体介电常数和填充物介电常数与所述立方体单元的等效阻抗之间的关系式,得到第一关系式,具体包括:
根据所述第一RC并联电路的参数,采用如下公式确定第一区域的第一RC并联电路的等效阻抗:
Figure BDA0002813866540000031
其中,
Figure BDA0002813866540000032
Figure BDA0002813866540000033
式中,ΔZ1为第一区域的第一RC并联电路的等效阻抗,ΔR1为第一区域的第一RC并联电路的等效电阻,ΔC1为第一区域的第一RC并联电路的等效电容,ω为角频率,ε'm为基体的复介电常数的实部,ε”m为基体的复介电常数的虚部,a为立方体单元的边长,Δs为为第一区域离散化为多个微元的微元底面积;
根据所述第二RC并联电路的参数、所述第三RC并联电路的参数和第四RC并联电路的参数,采用如下公式确定第二子区域的等效阻抗:
ΔZ2=Z2a+Z2m+Z2b
其中,
Figure BDA0002813866540000034
Figure BDA0002813866540000035
Figure BDA0002813866540000036
Figure BDA0002813866540000037
Figure BDA0002813866540000038
Figure BDA0002813866540000039
式中,ΔZ2为第二子区域的等效阻抗,Z2a为上基体等效阻抗,Z2b为下基体等效阻抗,ΔR2a为上基体的等效电阻,ΔR2b为下基体的等效电阻,ΔC2a为上基体的等效电容,ΔC2b为下基体的等效电容,R为填充物的半径,r为第二区域横截面上任一点到该横截面中心的距离,r<R,Z2m为填充物的等效阻抗,ΔR2m为填充物的等效电阻,ΔC2m为填充物的等效电容,ε'f为填充物的介电常数实部,ε”f为填充物的介电常数虚部;
根据所述第一区域的第一RC并联电路的等效阻抗和所述第二子区域的等效阻抗,采用如下公式确定所述立方体单元的等效阻抗,得到第一关系式:
Figure BDA0002813866540000041
式中,Zeq为立方体单元的等效阻抗。
可选的,所述根据所述第一关系式,建立基体介电常数、填充物介电常数和复合材料介电常数之间的关系式,得到第二关系式,具体包括:
对所述立方体单元的等效阻抗采用如下公式进行分解,得到所述立方体单元的等效电阻和所述立方体单元的等效电容:
Figure BDA0002813866540000042
Figure BDA0002813866540000043
其中,
Zeq=Z'eq-jZ”eq
式中,Ceq为立方体单元的等效电容,Req为立方体单元的等效电阻,Z'eq为立方体单元的等效阻抗实部,Z”eq为立方体单元的等效阻抗虚部;
根据所述立方体单元的等效电阻和所述立方体单元的等效电容,采用如下公式计算复合材料的介电常数,得到第二关系式:
Figure BDA0002813866540000044
Figure BDA0002813866540000045
式中,ε'eff为复合材料的介电常数实部,ε”eff为复合材料的介电常数虚部。
本发明还提供一种复合材料介电常数计算系统,包括:
复合材料划分模块,用于对复合材料进行单元划分,得到多个相同的立方体单元;每一个所述立方体单元均包括一个位于立方体单元中心位置处的填充物和设置在所述填充物外表面的基体;
等效电路生成模块,用于将所述立方体单元划分为第一区域和第二区域,根据所述第一区域和所述第二区域生成等效电路;其中,所述第一区域包括第一基体,所述第二区域包括上基体、填充物和下基体;所述等效电路包括并联的第一等效电路单元和第二等效电路单元;所述第一等效电路单元根据所述第一基体确定,所述第二等效电路单元根据所述上基体、所述填充物和所述下基体确定;
第一关系式确定模块,用于根据所述等效电路的参数,建立基体介电常数和填充物介电常数与所述立方体单元的等效阻抗之间的关系式,得到第一关系式;
第二关系式确定模块,用于根据所述第一关系式,建立基体介电常数、填充物介电常数和复合材料介电常数之间的关系式,得到第二关系式;
介电常数计算模块,用于获取基体、填充物、复合材料三种材料中的任两种材料的介电常数频谱,根据所述第二关系式采用最小二乘法确定另一种材料的介电常数频谱。
可选的,所述等效电路生成模块,具体包括:
第一等效电路单元生成单元,用于根据所述第一区域生成第一等效电路单元;所述第一等效电路单元包括多个相互并联的第一RC并联电路,所述第一RC并联电路的参数根据所述第一基体确定;
第二区域划分单元,用于在垂直于电场方向,将所述第二区域划分为多个第二子区域;
第二等效单元,用于将所述第二子区域等效为第二等效电路子单元;所述第二等效电路子单元为相互串联的第二RC并联电路、第三RC并联电路和第四RC并联电路,所述第二RC并联电路的参数根据所述上基体确定,所述第三RC并联电路的参数根据所述填充物确定,所述第四RC并联电路的参数根据所述下基体确定;
第二等效电路单元生成单元,用于将多个所述第二等效电路子单元并联,生成第二等效电路单元;
等效电路生成单元,用于将所述第一等效电路单元和所述第二等效电路单元并联,生成等效电路。
可选的,所述第一关系式确定模块,具体包括:
第一区域的第一RC并联电路的等效阻抗确定单元,用于根据所述第一RC并联电路的参数,采用如下公式确定第一区域的第一RC并联电路的等效阻抗:
Figure BDA0002813866540000061
其中,
Figure BDA0002813866540000062
Figure BDA0002813866540000063
式中,ΔZ1为第一区域的第一RC并联电路的等效阻抗,ΔR1为第一区域的第一RC并联电路的等效电阻,ΔC1为第一区域的第一RC并联电路的等效电容,ω为角频率,ε'm为基体的复介电常数的实部,ε”m为基体的复介电常数的虚部,a为立方体单元的边长,Δs为第一区域离散化为多个微元的微元底面积;
第二子区域的等效阻抗确定单元,用于根据所述第二RC并联电路的参数、所述第三RC并联电路的参数和第四RC并联电路的参数,采用如下公式确定第二子区域的等效阻抗:
ΔZ2=Z2a+Z2m+Z2b
其中,
Figure BDA0002813866540000064
Figure BDA0002813866540000065
Figure BDA0002813866540000066
Figure BDA0002813866540000067
Figure BDA0002813866540000068
Figure BDA0002813866540000071
式中,ΔZ2为第二子区域的等效阻抗,Z2a为上基体等效阻抗,Z2b为下基体等效阻抗,ΔR2a为上基体的等效电阻,ΔR2b为下基体的等效电阻,ΔC2a为上基体的等效电容,ΔC2b为下基体的等效电容,R为填充物的半径,r为第二区域横截面上任一点到该横截面中心的距离,r<R,Z2m为填充物的等效阻抗,ΔR2m为填充物的等效电阻,ΔC2m为填充物的等效电容,ε'f为填充物的介电常数实部,ε”f为填充物的介电常数虚部;
第一关系式确定单元,用于根据所述第一区域的第一RC并联电路的等效阻抗和所述第二子区域的等效阻抗,采用如下公式确定所述立方体单元的等效阻抗,得到第一关系式:
Figure BDA0002813866540000072
式中,Zeq为立方体单元的等效阻抗。
可选的,所述第二关系式确定模块,具体包括:
等效阻抗分解单元,用于对所述立方体单元的等效阻抗采用如下公式进行分解,得到所述立方体单元的等效电阻和所述立方体单元的等效电容:
Figure BDA0002813866540000073
Figure BDA0002813866540000074
其中,
Zeq=Z'eq-jZ”eq
式中,Ceq为立方体单元的等效电容,Req为立方体单元的等效电阻,Z'eq为立方体单元的等效阻抗实部,Z”eq为立方体单元的等效阻抗虚部;
第二关系式确定单元,用于根据所述立方体单元的等效电阻和所述立方体单元的等效电容,采用如下公式计算复合材料的介电常数,得到第二关系式:
Figure BDA0002813866540000075
Figure BDA0002813866540000081
式中,ε'eff为复合材料的介电常数实部,ε”eff为复合材料的介电常数虚部。
与现有技术相比,本发明的有益效果是:
本发明提出了一种复合材料介电常数计算方法及系统,将复合材料划分为立方体单元结构,对立方体结构进行了区域划分,对界面极化不规则区域进行了电路等效,考虑了不同材料间的界面极化,不但能计算高频段的介电谱,也能计算复合材料的中、低频段的介电谱;并且利用第二关系式,能够根据基体和填充物的介电谱线计算复合材料的介电常数频谱,亦可根据复合材料及其中一种材料的介电谱线反推另外一种材料的介电常数频谱,极大降低了工作量,节约了测试时间。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例中复合材料介电常数计算方法流程图;
图2为本发明实施例中立方体单元示意图;
图3为本发明实施例中等效电路示意图;
图4为本发明实施例中复合材料介电常数计算系统结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种复合材料介电常数计算方法及系统,考虑了不同材料间的界面极化,能够高效计算高、中、低频段的介电谱,极大降低了工作量,节约了测试时间。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例
图1为本发明实施例中复合材料介电常数计算方法流程图,如图1所示,一种复合材料介电常数计算方法,包括:
步骤101:对复合材料进行单元划分,得到多个相同且独立的立方体单元;每一个立方体单元均包括一个位于立方体单元中心位置处的填充物和在填充物外表面设置的基体。
如图2所示,把复合材料分为m×n×h个相同的独立立方单元,大小为a×b×c。每个立方体都包含一个位于中心的任意形状介质。本发明实施例为边长为a的立方体,其中含有半径为R的球形介质,立方体中间的球形介质是充填物,球形介质体积在立方体中的占比等于充填物在复合材料的体积占比,其他部分的介质是基体,体积占比可采用浸渍法等方法进行求取。把与球体相切的圆柱体确定为区域2(Area2),立方体其余部分为区域1(Area1)。
步骤102:将立方体单元划分为第一区域(区域1)和第二区域(区域2),根据第一区域和第二区域生成等效电路;其中,第一区域包括第一基体,第二区域包括上基体、填充物和下基体;等效电路包括并联的第一等效电路单元和第二等效电路单元;第一等效电路单元根据第一基体确定,第二等效电路单元根据上基体、填充物和下基体确定。
步骤102,具体包括:
根据第一区域生成第一等效电路单元;第一等效电路单元包括多个相互并联的第一RC并联电路,第一RC并联电路的参数根据第一基体确定;
在垂直于电场方向,将第二区域划分为多个第二子区域;
将第二子区域等效为第二等效电路子单元;第二等效电路子单元为相互串联的第二RC并联电路、第三RC并联电路和第四RC并联电路,第二RC并联电路的参数根据上基体确定,第三RC并联电路的参数根据填充物确定,第四RC并联电路的参数根据下基体确定;
将多个第二等效电路子单元并联,生成第二等效电路单元;
将第一等效电路单元和第二等效电路单元并联,生成等效电路。
假设立方体单元中的介质均为线性各向同性、均匀介质。沿z轴方向施加交变电场。将每个立方体单元离散成两个区间,如图2所示。立方单元中的各种介质可以用电容及电阻来进行等效。将第二区域的上下表面离散成若干个边长为Δa(Δa→0)的正方形。如图3所示,第二子区域个数为n个,ΔR1n为第一区域的等效电阻,ΔC1n为第一区域的等效电容,ΔRan为第n个第二子区域的上基体的等效电阻,ΔCan为第n个第二子区域的上基体的等效电容,ΔRmn为第n个第二子区域的填充物的等效电阻,ΔCmn为第n个第二子区域的填充物的等效电容,ΔRbn为第n个第二子区域的下基体的等效电阻,ΔCbn为第n个第二子区域的下基体的等效电容。
步骤103:根据等效电路的参数,建立基体介电常数和填充物介电常数与立方体单元的等效阻抗之间的关系式,得到第一关系式;基体包括第一基体、上基体和下基体。
步骤103,具体包括:
根据第一RC并联电路的参数,采用如下公式确定第一区域的第一RC并联电路的等效阻抗:
Figure BDA0002813866540000101
其中,
Figure BDA0002813866540000102
Figure BDA0002813866540000103
式中,ΔZ1为第一区域的第一RC并联电路的等效阻抗,ΔR1为第一区域的第一RC并联电路的等效电阻,ΔC1为第一区域的第一RC并联电路的等效电容,ω为角频率,ε'm为在角频率ω对应的频率下的基体的复介电常数的实部,ε”m为在角频率ω对应的频率下的基体的复介电常数的虚部,a为立方体单元的边长,Δs为第一区域离散化为多个微元的微元底面积,即将区域的上下表面离散成若干个边长为Δa(Δa→0)的正方形的面积;Δs的总和为(a2-πR2),即
Figure BDA0002813866540000104
根据第二RC并联电路的参数、第三RC并联电路的参数和第四RC并联电路的参数,采用如下公式确定第二子区域的等效阻抗:
ΔZ2=Z2a+Z2m+Z2b
其中,
Figure BDA0002813866540000111
Figure BDA0002813866540000112
Figure BDA0002813866540000113
Figure BDA0002813866540000114
Figure BDA0002813866540000115
Figure BDA0002813866540000116
式中,ΔZ2为第二子区域的等效阻抗,Z2a为上基体等效阻抗,Z2b为下基体等效阻抗,ΔR2a为上基体的等效电阻,ΔR2b为下基体的等效电阻,ΔC2a为上基体的等效电容,ΔC2b为下基体的等效电容,R为填充物的半径,r为第二区域横截面上任一点到该横截面中心的距离,r<R,Z2m为填充物的等效阻抗,ΔR2m为填充物的等效电阻,ΔC2m为填充物的等效电容,ε'f为填充物的介电常数实部,ε”f为填充物的介电常数虚部;
根据第一区域的第一RC并联电路的等效阻抗和第二子区域的等效阻抗,采用如下公式确定立方体单元的等效阻抗,得到第一关系式:
Figure BDA0002813866540000117
式中,Zeq为立方体单元的等效阻抗。
步骤104:根据第一关系式,建立基体介电常数、填充物介电常数和复合材料介电常数之间的关系式,得到第二关系式。
步骤104,具体包括:
对立方体单元的等效阻抗采用如下公式进行分解,得到立方体单元的等效电阻和立方体单元的等效电容:
Figure BDA0002813866540000121
Figure BDA0002813866540000122
其中,
Zeq=Z'eq-jZ”eq
式中,Ceq为立方体单元的等效电容,Req为立方体单元的等效电阻,Z'eq为立方体单元的等效阻抗实部,Z”eq为立方体单元的等效阻抗虚部;
根据立方体单元的等效电阻和立方体单元的等效电容,采用如下公式计算复合材料的介电常数,得到第二关系式:
Figure BDA0002813866540000123
Figure BDA0002813866540000124
式中,ε'eff为复合材料的介电常数实部,ε”eff为复合材料的介电常数虚部。
其中,
由m×n×h个相同的立方体单元组成的复合材料的总等效电容和总等效电阻为:
Figure BDA0002813866540000125
Figure BDA0002813866540000126
由于样品的介电常数等于单个立方体的介电常数,因此可以使用单个立方体单元进行计算:
Figure BDA0002813866540000127
Figure BDA0002813866540000128
步骤105:获取基体、填充物、复合材料三种材料中的任两种材料的介电常数频谱,根据第二关系式采用最小二乘法确定另一种材料的介电常数频谱。
通过python、matlab等能够进行拟合计算的数学软件,根据上述建立的基体、充填物及复合材料的介电常数之间的关系,输入任意两者的介电常数及电导率,通过最小二乘法计算另外一种材料的介电常数并绘制介电频谱。
图4为本发明实施例中复合材料介电常数计算系统结构图。如图4所示,一种复合材料介电常数计算系统,包括:
复合材料划分模块201,用于对复合材料进行单元划分,得到多个相同的立方体单元;每一个立方体单元均包括一个位于立方体单元中心位置处的填充物和在填充物外表面设置的基体。
等效电路生成模块202,用于将立方体单元划分为第一区域和第二区域,根据第一区域和第二区域生成等效电路;其中,第一区域包括第一基体,第二区域包括上基体、填充物和下基体;等效电路包括并联的第一等效电路单元和第二等效电路单元;第一等效电路单元根据第一基体确定,第二等效电路单元根据上基体、填充物和下基体确定。
等效电路生成模块202,具体包括:
第一等效电路单元生成单元,用于根据第一区域生成第一等效电路单元;第一等效电路单元包括多个相互并联的第一RC并联电路,第一RC并联电路的参数根据第一基体确定;
第二区域划分单元,用于在垂直于电场方向,将第二区域划分为多个第二子区域;
第二等效单元,用于将第二子区域等效为第二等效电路子单元;第二等效电路子单元为相互串联的第二RC并联电路、第三RC并联电路和第四RC并联电路,第二RC并联电路的参数根据上基体确定,第三RC并联电路的参数根据填充物确定,第四RC并联电路的参数根据下基体确定;
第二等效电路单元生成单元,用于将多个第二等效电路子单元并联,生成第二等效电路单元;
等效电路生成单元,用于将第一等效电路单元和第二等效电路单元并联,生成等效电路。
第一关系式确定模块203,用于根据等效电路的参数,建立基体介电常数和填充物介电常数与立方体单元的等效阻抗之间的关系式,得到第一关系式;基体包括第一基体、上基体和下基体。
第一关系式确定模块203,具体包括:
第一区域的第一RC并联电路的等效阻抗确定单元,用于根据第一RC并联电路的参数,采用如下公式确定第一区域的第一RC并联电路的等效阻抗:
Figure BDA0002813866540000141
其中,
Figure BDA0002813866540000142
Figure BDA0002813866540000143
式中,ΔZ1为第一区域的第一RC并联电路的等效阻抗,ΔR1为第一区域的第一RC并联电路的等效电阻,ΔC1为第一区域的第一RC并联电路的等效电容,ω为角频率,ε'm为在角频率ω对应的频率下的基体的复介电常数的实部,ε”m为在角频率ω对应的频率下的基体的复介电常数的虚部,a为立方体单元的边长,Δs为第一区域离散化为多个微元的微元底面积;
第二子区域的等效阻抗确定单元,用于根据第二RC并联电路的参数、第三RC并联电路的参数和第四RC并联电路的参数,采用如下公式确定第二子区域的等效阻抗:
ΔZ2=Z2a+Z2m+Z2b
其中,
Figure BDA0002813866540000144
Figure BDA0002813866540000145
Figure BDA0002813866540000146
Figure BDA0002813866540000147
Figure BDA0002813866540000148
Figure BDA0002813866540000149
式中,ΔZ2为第二子区域的等效阻抗,Z2a为上基体等效阻抗,Z2b为下基体等效阻抗,ΔR2a为上基体的等效电阻,ΔR2b为下基体的等效电阻,ΔC2a为上基体的等效电容,ΔC2b为下基体的等效电容,R为填充物的半径,r为第二区域横截面上任一点到该横截面中心的距离,r<R,Z2m为填充物的等效阻抗,ΔR2m为填充物的等效电阻,ΔC2m为填充物的等效电容,ε'f为填充物的介电常数实部,ε”f为填充物的介电常数虚部;
第一关系式确定单元,用于根据第一区域的第一RC并联电路的等效阻抗和第二子区域的等效阻抗,采用如下公式确定立方体单元的等效阻抗,得到第一关系式:
Figure BDA0002813866540000151
式中,Zeq为立方体单元的等效阻抗。
第二关系式确定模块204,用于根据第一关系式,建立基体介电常数、填充物介电常数和复合材料介电常数之间的关系式,得到第二关系式。
第二关系式确定模块204,具体包括:
等效阻抗分解单元,用于对立方体单元的等效阻抗采用如下公式进行分解,得到立方体单元的等效电阻和立方体单元的等效电容:
Figure BDA0002813866540000152
Figure BDA0002813866540000153
其中,
Zeq=Z'eq-jZ”eq
式中,Ceq为立方体单元的等效电容,Req为立方体单元的等效电阻,Z'eq为立方体单元的等效阻抗实部,Z”eq为立方体单元的等效阻抗虚部;
第二关系式确定单元,用于根据立方体单元的等效电阻和立方体单元的等效电容,采用如下公式计算复合材料的介电常数,得到第二关系式:
Figure BDA0002813866540000154
Figure BDA0002813866540000161
式中,ε'eff为复合材料的介电常数实部,ε”eff为复合材料的介电常数虚部。
介电常数计算模块205,用于获取基体、填充物、复合材料三种材料中的任两种材料的介电常数频谱,根据第二关系式采用最小二乘法确定另一种材料的介电常数频谱。
对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上,本说明书内容不应理解为对本发明的限制。

Claims (8)

1.一种复合材料介电常数计算方法,其特征在于,包括:
对复合材料进行单元划分,得到多个相同的立方体单元;每一个所述立方体单元均包括一个位于立方体单元中心位置处的填充物和设置在所述填充物外表面的基体;
将所述立方体单元划分为第一区域和第二区域,根据所述第一区域和所述第二区域生成等效电路;其中,所述第一区域包括第一基体,所述第二区域包括上基体、填充物和下基体;所述等效电路包括并联的第一等效电路单元和第二等效电路单元;所述第一等效电路单元根据所述第一基体确定,所述第二等效电路单元根据所述上基体、所述填充物和所述下基体确定;
根据所述等效电路的参数,建立基体介电常数和填充物介电常数与所述立方体单元的等效阻抗之间的关系式,得到第一关系式;
根据所述第一关系式,建立基体介电常数、填充物介电常数和复合材料介电常数之间的关系式,得到第二关系式;
获取基体、填充物、复合材料三种材料中的任两种材料的介电常数频谱,根据所述第二关系式采用最小二乘法确定另一种材料的介电常数频谱。
2.根据权利要求1所述的复合材料介电常数计算方法,其特征在于,所述根据所述第一区域和所述第二区域生成等效电路,具体包括:
根据所述第一区域生成第一等效电路单元;所述第一等效电路单元包括多个相互并联的第一RC并联电路,所述第一RC并联电路的参数根据所述第一基体确定;
在垂直于电场方向,将所述第二区域划分为多个第二子区域;
将所述第二子区域等效为第二等效电路子单元;所述第二等效电路子单元为相互串联的第二RC并联电路、第三RC并联电路和第四RC并联电路,所述第二RC并联电路的参数根据所述上基体确定,所述第三RC并联电路的参数根据所述填充物确定,所述第四RC并联电路的参数根据所述下基体确定;
将多个所述第二等效电路子单元并联,生成第二等效电路单元;
将所述第一等效电路单元和所述第二等效电路单元并联,生成等效电路。
3.根据权利要求2所述的复合材料介电常数计算方法,其特征在于,所述根据所述等效电路的参数,建立基体介电常数和填充物介电常数与所述立方体单元的等效阻抗之间的关系式,得到第一关系式,具体包括:
根据所述第一RC并联电路的参数,采用如下公式确定第一区域的第一RC并联电路的等效阻抗:
Figure FDA0002813866530000021
其中,
Figure FDA0002813866530000022
Figure FDA0002813866530000023
式中,ΔZ1为第一区域的第一RC并联电路的等效阻抗,ΔR1为第一区域的第一RC并联电路的等效电阻,ΔC1为第一区域的第一RC并联电路的等效电容,ω为角频率,ε'm为基体的复介电常数的实部,ε″m为基体的复介电常数的虚部,a为立方体单元的边长,Δs为第一区域离散化为多个微元的微元底面积;
根据所述第二RC并联电路的参数、所述第三RC并联电路的参数和第四RC并联电路的参数,采用如下公式确定第二子区域的等效阻抗:
ΔZ2=Z2a+Z2m+Z2b
其中,
Figure FDA0002813866530000024
Figure FDA0002813866530000025
Figure FDA0002813866530000026
Figure FDA0002813866530000027
Figure FDA0002813866530000028
Figure FDA0002813866530000029
式中,ΔZ2为第二子区域的等效阻抗,Z2a为上基体等效阻抗,Z2b为下基体等效阻抗,ΔR2a为上基体的等效电阻,ΔR2b为下基体的等效电阻,ΔC2a为上基体的等效电容,ΔC2b为下基体的等效电容,R为填充物的半径,r为第二区域横截面上任一点到该横截面中心的距离,r<R,Z2m为填充物的等效阻抗,ΔR2m为填充物的等效电阻,ΔC2m为填充物的等效电容,ε'f为填充物的介电常数实部,ε″f为填充物的介电常数虚部;
根据所述第一区域的第一RC并联电路的等效阻抗和所述第二子区域的等效阻抗,采用如下公式确定所述立方体单元的等效阻抗,得到第一关系式:
Figure FDA0002813866530000031
式中,Zeq为立方体单元的等效阻抗。
4.根据权利要求3所述的复合材料介电常数计算方法,其特征在于,所述根据所述第一关系式,建立基体介电常数、填充物介电常数和复合材料介电常数之间的关系式,得到第二关系式,具体包括:
对所述立方体单元的等效阻抗采用如下公式进行分解,得到所述立方体单元的等效电阻和所述立方体单元的等效电容:
Figure FDA0002813866530000032
Figure FDA0002813866530000033
其中,
Zeq=Z′eq-jZ″eq
式中,Ceq为立方体单元的等效电容,Req为立方体单元的等效电阻,Z'eq为立方体单元的等效阻抗实部,Z″eq为立方体单元的等效阻抗虚部;
根据所述立方体单元的等效电阻和所述立方体单元的等效电容,采用如下公式计算复合材料的介电常数,得到第二关系式:
Figure FDA0002813866530000034
Figure FDA0002813866530000035
式中,ε'eff为复合材料的介电常数实部,ε″eff为复合材料的介电常数虚部。
5.一种复合材料介电常数计算系统,其特征在于,包括:
复合材料划分模块,用于对复合材料进行单元划分,得到多个相同的立方体单元;每一个所述立方体单元均包括一个位于立方体单元中心位置处的填充物和设置在所述填充物外表面的基体;
等效电路生成模块,用于将所述立方体单元划分为第一区域和第二区域,根据所述第一区域和所述第二区域生成等效电路;其中,所述第一区域包括第一基体,所述第二区域包括上基体、填充物和下基体;所述等效电路包括并联的第一等效电路单元和第二等效电路单元;所述第一等效电路单元根据所述第一基体确定,所述第二等效电路单元根据所述上基体、所述填充物和所述下基体确定;
第一关系式确定模块,用于根据所述等效电路的参数,建立基体介电常数和填充物介电常数与所述立方体单元的等效阻抗之间的关系式,得到第一关系式;
第二关系式确定模块,用于根据所述第一关系式,建立基体介电常数、填充物介电常数和复合材料介电常数之间的关系式,得到第二关系式;
介电常数计算模块,用于获取基体、填充物、复合材料三种材料中的任两种材料的介电常数频谱,根据所述第二关系式采用最小二乘法确定另一种材料的介电常数频谱。
6.根据权利要求5所述的复合材料介电常数计算系统,其特征在于,所述等效电路生成模块,具体包括:
第一等效电路单元生成单元,用于根据所述第一区域生成第一等效电路单元;所述第一等效电路单元包括多个相互并联的第一RC并联电路,所述第一RC并联电路的参数根据所述第一基体确定;
第二区域划分单元,用于在垂直于电场方向,将所述第二区域划分为多个第二子区域;
第二等效单元,用于将所述第二子区域等效为第二等效电路子单元;所述第二等效电路子单元为相互串联的第二RC并联电路、第三RC并联电路和第四RC并联电路,所述第二RC并联电路的参数根据所述上基体确定,所述第三RC并联电路的参数根据所述填充物确定,所述第四RC并联电路的参数根据所述下基体确定;
第二等效电路单元生成单元,用于将多个所述第二等效电路子单元并联,生成第二等效电路单元;
等效电路生成单元,用于将所述第一等效电路单元和所述第二等效电路单元并联,生成等效电路。
7.根据权利要求6所述的复合材料介电常数计算系统,其特征在于,所述第一关系式确定模块,具体包括:
第一区域的第一RC并联电路的等效阻抗确定单元,用于根据所述第一RC并联电路的参数,采用如下公式确定第一区域的第一RC并联电路的等效阻抗:
Figure FDA0002813866530000051
其中,
Figure FDA0002813866530000052
Figure FDA0002813866530000053
式中,ΔZ1为第一区域的第一RC并联电路的等效阻抗,ΔR1为第一区域的第一RC并联电路的等效电阻,ΔC1为第一区域的第一RC并联电路的等效电容,ω为角频率,ε'm为基体的复介电常数的实部,ε″m为基体的复介电常数的虚部,a为立方体单元的边长,Δs为第一区域离散化为多个微元的微元底面积;
第二子区域的等效阻抗确定单元,用于根据所述第二RC并联电路的参数、所述第三RC并联电路的参数和第四RC并联电路的参数,采用如下公式确定第二子区域的等效阻抗:
ΔZ2=Z2a+Z2m+Z2b
其中,
Figure FDA0002813866530000054
Figure FDA0002813866530000055
Figure FDA0002813866530000061
Figure FDA0002813866530000062
Figure FDA0002813866530000063
Figure FDA0002813866530000064
式中,ΔZ2为第二子区域的等效阻抗,Z2a为上基体等效阻抗,Z2b为下基体等效阻抗,ΔR2a为上基体的等效电阻,ΔR2b为下基体的等效电阻,ΔC2a为上基体的等效电容,ΔC2b为下基体的等效电容,R为填充物的半径,r为第二区域横截面上任一点到该横截面中心的距离,r<R,Z2m为填充物的等效阻抗,ΔR2m为填充物的等效电阻,ΔC2m为填充物的等效电容,ε'f为填充物的介电常数实部,ε″f为填充物的介电常数虚部;
第一关系式确定单元,用于根据所述第一区域的第一RC并联电路的等效阻抗和所述第二子区域的等效阻抗,采用如下公式确定所述立方体单元的等效阻抗,得到第一关系式:
Figure FDA0002813866530000065
式中,Zeq为立方体单元的等效阻抗。
8.根据权利要求7所述的复合材料介电常数计算系统,其特征在于,所述第二关系式确定模块,具体包括:
等效阻抗分解单元,用于对所述立方体单元的等效阻抗采用如下公式进行分解,得到所述立方体单元的等效电阻和所述立方体单元的等效电容:
Figure FDA0002813866530000066
Figure FDA0002813866530000067
其中,
Zeq=Z′eq-jZ″eq
式中,Ceq为立方体单元的等效电容,Req为立方体单元的等效电阻,Z'eq为立方体单元的等效阻抗实部,Z″eq为立方体单元的等效阻抗虚部;
第二关系式确定单元,用于根据所述立方体单元的等效电阻和所述立方体单元的等效电容,采用如下公式计算复合材料的介电常数,得到第二关系式:
Figure FDA0002813866530000071
Figure FDA0002813866530000072
式中,ε'eff为复合材料的介电常数实部,ε″eff为复合材料的介电常数虚部。
CN202011405208.3A 2020-12-03 2020-12-03 一种复合材料介电常数计算方法及系统 Active CN112683965B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011405208.3A CN112683965B (zh) 2020-12-03 2020-12-03 一种复合材料介电常数计算方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011405208.3A CN112683965B (zh) 2020-12-03 2020-12-03 一种复合材料介电常数计算方法及系统

Publications (2)

Publication Number Publication Date
CN112683965A true CN112683965A (zh) 2021-04-20
CN112683965B CN112683965B (zh) 2021-11-09

Family

ID=75446033

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011405208.3A Active CN112683965B (zh) 2020-12-03 2020-12-03 一种复合材料介电常数计算方法及系统

Country Status (1)

Country Link
CN (1) CN112683965B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113504571A (zh) * 2021-07-23 2021-10-15 吉林大学 一种多相导电媒质的极化等效环装置和设计方法
CN116151064A (zh) * 2022-12-30 2023-05-23 西南大学 一种油纸绝缘套管芯体的介电参数多分区反演方法

Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04198847A (ja) * 1990-11-29 1992-07-20 Oyama Kogyo Koutou Senmon Gatsukouchiyou 含水率の測定方法
US20100217576A1 (en) * 2009-02-23 2010-08-26 Georgia Tech Research Corporation Multi-Layer Finite Element Method for Modeling of Package Power and Ground Planes
US20110102075A1 (en) * 2006-06-06 2011-05-05 Teraspeed Consulting Group Llc Power distribution system for integrated circuits
CN102281708A (zh) * 2011-07-22 2011-12-14 西安电子科技大学 基于交指电容电磁带隙结构的电路板
US20120116683A1 (en) * 2010-11-09 2012-05-10 General Electric Company Highly selective chemical and biological sensors
CN102510658A (zh) * 2011-09-26 2012-06-20 北京邮电大学 一种面向多频天线衬底的h型槽分形uc-ebg结构的实现方法
EP2809130A1 (en) * 2012-01-25 2014-12-03 Konica Minolta, Inc. Evaluation method, evaluation device, evaluation program, recording medium, and manufacturing method for organic electroluminescent element
CN105419277A (zh) * 2015-12-11 2016-03-23 南京华格电汽塑业有限公司 一种低介电性能复合材料
CN105825037A (zh) * 2016-06-13 2016-08-03 华北电力大学 变电站气体绝缘开关设备管道支架的等效电路建模方法
CN106751238A (zh) * 2016-11-17 2017-05-31 成都市创斯德机电设备有限公司 一种介电复合材料及其制备方法
US20170287792A1 (en) * 2016-03-30 2017-10-05 Tsinghua University General four-port on-wafer high frequency de-embedding method
CN108170950A (zh) * 2017-12-27 2018-06-15 电子科技大学 基于神经网络的多层频率选择表面吸波材料建模优化方法
JP2018132458A (ja) * 2017-02-17 2018-08-23 桓達科技股▲ふん▼有限公司FINETEK Co.,Ltd. 検知装置及び物質検知方法
CN109142865A (zh) * 2018-07-27 2019-01-04 福州大学 考虑油纸绝缘界面极化等效电路参数的频域谱辨识方法
US20190187072A1 (en) * 2017-12-18 2019-06-20 United Arab Emirates University Plant water sensor
CN110456208A (zh) * 2019-09-12 2019-11-15 河海大学 一种压电振子电学特性分布的估算方法及可读存储介质
CN111398117A (zh) * 2020-03-15 2020-07-10 长江大学 一种岩心电频散分析方法及装置

Patent Citations (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04198847A (ja) * 1990-11-29 1992-07-20 Oyama Kogyo Koutou Senmon Gatsukouchiyou 含水率の測定方法
US20110102075A1 (en) * 2006-06-06 2011-05-05 Teraspeed Consulting Group Llc Power distribution system for integrated circuits
US20100217576A1 (en) * 2009-02-23 2010-08-26 Georgia Tech Research Corporation Multi-Layer Finite Element Method for Modeling of Package Power and Ground Planes
US20120116683A1 (en) * 2010-11-09 2012-05-10 General Electric Company Highly selective chemical and biological sensors
CN102281708A (zh) * 2011-07-22 2011-12-14 西安电子科技大学 基于交指电容电磁带隙结构的电路板
CN102510658A (zh) * 2011-09-26 2012-06-20 北京邮电大学 一种面向多频天线衬底的h型槽分形uc-ebg结构的实现方法
EP2809130A1 (en) * 2012-01-25 2014-12-03 Konica Minolta, Inc. Evaluation method, evaluation device, evaluation program, recording medium, and manufacturing method for organic electroluminescent element
CN105419277A (zh) * 2015-12-11 2016-03-23 南京华格电汽塑业有限公司 一种低介电性能复合材料
US20170287792A1 (en) * 2016-03-30 2017-10-05 Tsinghua University General four-port on-wafer high frequency de-embedding method
CN105825037A (zh) * 2016-06-13 2016-08-03 华北电力大学 变电站气体绝缘开关设备管道支架的等效电路建模方法
CN106751238A (zh) * 2016-11-17 2017-05-31 成都市创斯德机电设备有限公司 一种介电复合材料及其制备方法
JP2018132458A (ja) * 2017-02-17 2018-08-23 桓達科技股▲ふん▼有限公司FINETEK Co.,Ltd. 検知装置及び物質検知方法
US20190187072A1 (en) * 2017-12-18 2019-06-20 United Arab Emirates University Plant water sensor
CN108170950A (zh) * 2017-12-27 2018-06-15 电子科技大学 基于神经网络的多层频率选择表面吸波材料建模优化方法
CN109142865A (zh) * 2018-07-27 2019-01-04 福州大学 考虑油纸绝缘界面极化等效电路参数的频域谱辨识方法
CN110456208A (zh) * 2019-09-12 2019-11-15 河海大学 一种压电振子电学特性分布的估算方法及可读存储介质
CN111398117A (zh) * 2020-03-15 2020-07-10 长江大学 一种岩心电频散分析方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
MUFEI XIAO: ""Equivalent-circuit interpretation of the dielectric dispersion in ferroelectric superlattice capacitor"", 《MATERIALS LETTERS》 *
董明 等: ""油纸绝缘频域介电谱解释方法研究"", 《中国电机工程学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113504571A (zh) * 2021-07-23 2021-10-15 吉林大学 一种多相导电媒质的极化等效环装置和设计方法
CN113504571B (zh) * 2021-07-23 2022-05-03 吉林大学 一种多相导电媒质的极化等效环装置和设计方法
CN116151064A (zh) * 2022-12-30 2023-05-23 西南大学 一种油纸绝缘套管芯体的介电参数多分区反演方法
CN116151064B (zh) * 2022-12-30 2023-08-01 西南大学 一种油纸绝缘套管芯体的介电参数多分区反演方法

Also Published As

Publication number Publication date
CN112683965B (zh) 2021-11-09

Similar Documents

Publication Publication Date Title
CN112683965B (zh) 一种复合材料介电常数计算方法及系统
Brosseau et al. Dielectric and microstructure properties of polymer carbon black composites
Mauritz et al. State of understanding of Nafion
Zhu et al. A theoretically rigorous full-wave finite-element-based solution of Maxwell's equations from DC to high frequencies
Koledintseva et al. PCB conductor surface roughness as a layer with effective material parameters
CN104142117A (zh) 用于测量印刷电路板中的层的厚度的方法及设备
CN108519261B (zh) 一种基于三明治结构的半导电材料介电性能测试方法
Chen et al. Analysis of a capacitive sensor for the evaluation of circular cylinders with a conductive core
Penaranda-Foix et al. Full-wave analysis of dielectric-loaded cylindrical waveguides and cavities using a new four-port ring network
CN110472356B (zh) 一种电磁波多入射角度下的复合材料电磁参数计算方法
Flandin et al. A 3‐D numerical simulation of AC electrical properties of short fiber composites
Vadlamudi et al. Multiphysics based simulation of damage progression in composites
CN112364524A (zh) 一种多铺层碳纤维复合材料的宽频段电磁参数获取方法
Afrasiab et al. Electrical and mechanical characterization of the gas diffusion layer during compression in PEM fuel cells
Casieri et al. Evidence of temperature-induced subdiffusion of water on the micrometer scale in a Nafion membrane
Lachini et al. The simulation of electric field distribution on cable under the presence of moisture and air voids
Gil’mutdinov et al. Physical implementation of elements with fractal impedance: state of the art and prospects
Ma Modeling of conductor surface roughness effect
Mounkaila et al. Modelling the electrical behaviour of carbon/epoxy composites and monitoring changes in their microstructure during oven and autoclave curing using electrical impedancemetry
Pan et al. An experimental study on the complex resistivity of fractured rock under different saturation conditions
CN109917254A (zh) 一种油浸式套管内绝缘受潮的频域介电谱建模方法
Cerato et al. Dielectric measurement of soil-electrolyte mixtures in a modified oedometer cell using 400 kHz to 20 MHz electromagnetic waves
Han et al. Dielectric behavior of suspensions of polystyrene− zinc oxide composite microspheres
Hamar et al. Permittivity determination of composite materials based on a 3D electric field model
Wu et al. Thin Absorber with Dispersive Medium Simulated by Conformal Hybrid Explicit Implicit FDTD Method

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