王 嵐 惠小強(qiáng)

(1西安郵電大學(xué)理學(xué)院； 2西安郵電大學(xué)物聯(lián)網(wǎng)與兩化融合研究院，陜西 西安 710061)

薄圓盤直徑端點(diǎn)之間電阻的研究

王 嵐1惠小強(qiáng)2

(1西安郵電大學(xué)理學(xué)院；2西安郵電大學(xué)物聯(lián)網(wǎng)與兩化融合研究院，陜西 西安 710061)

針對(duì)薄圓盤直徑端點(diǎn)之間電阻難以利用傳統(tǒng)直接積分法或數(shù)值解法測(cè)量的問(wèn)題，本文利用有限元分析法研究了薄圓盤直徑端點(diǎn)之間電阻的測(cè)量方法。利用Comsol multiphysics設(shè)計(jì)薄圓盤電阻測(cè)量平臺(tái)，通過(guò)在薄圓盤直徑端點(diǎn)設(shè)計(jì)電極，研究了薄圓盤直徑端點(diǎn)之間電阻隨薄圓盤半徑、電極位置、薄圓盤上小孔位置和大小的變化規(guī)律；利用最小二乘法擬合薄圓盤直徑端點(diǎn)之間電阻表達(dá)式，實(shí)現(xiàn)薄圓盤直徑端點(diǎn)之間電阻間接測(cè)量；在典型條件下，利用直接積分測(cè)量法對(duì)上述測(cè)量方法進(jìn)行驗(yàn)證，結(jié)果表明本文提出的方法誤差較小。本文研究工作對(duì)薄圓盤電阻的設(shè)計(jì)具有指導(dǎo)意義，且該研究方法還可以推廣用于其他形狀電阻的設(shè)計(jì)。

薄圓盤；電阻；Comsol multiphysics

穩(wěn)恒條件下導(dǎo)體內(nèi)電勢(shì)/電流的分布、不規(guī)則或不均勻?qū)w電阻的測(cè)量等都是電磁學(xué)應(yīng)用研究中的基本問(wèn)題[1，2]。目前，針對(duì)高對(duì)稱性物體主要通過(guò)直接積分測(cè)量法計(jì)算電阻值，而針對(duì)低對(duì)稱性物體一般采用數(shù)值解法[3，4]。

縱觀國(guó)內(nèi)外公開文獻(xiàn)[5-8]，尚未見(jiàn)到有關(guān)薄圓盤直徑端點(diǎn)之間電阻研究的相關(guān)報(bào)道或文獻(xiàn)，其主要原因是薄圓盤中電流方向和等勢(shì)線的關(guān)系復(fù)雜，其電阻難以直接利用直接積分測(cè)量法或數(shù)值解法計(jì)算。

針對(duì)上述問(wèn)題，本文研究了薄圓盤直徑端點(diǎn)之間電阻的測(cè)量方法。首先，利用Comsol multiphysics[9]軟件設(shè)計(jì)了薄圓盤電阻測(cè)量平臺(tái)，研究了薄圓盤直徑端點(diǎn)之間電阻隨薄圓盤半徑、電極位置、薄圓盤上小孔位置和大小等因素的變化規(guī)律；其次，利用最小二乘法擬合薄圓盤直徑端點(diǎn)之間電阻表達(dá)式，實(shí)現(xiàn)薄圓盤直徑端點(diǎn)之間電阻間接測(cè)量；最后，在典型條件下，利用直接積分測(cè)量法對(duì)上述測(cè)量方法進(jìn)行驗(yàn)證。

1 薄圓盤直徑端點(diǎn)之間電阻測(cè)量平臺(tái)設(shè)計(jì)

為實(shí)現(xiàn)薄圓盤直徑端點(diǎn)之間電阻測(cè)量平臺(tái)設(shè)計(jì)，本文利用有限元分析法建立電導(dǎo)率固定、半徑為a的薄圓盤；在薄圓盤上挖一個(gè)半徑為r的小孔，并補(bǔ)上不同電導(dǎo)率的介質(zhì)；在薄圓盤邊緣設(shè)計(jì)電極，如圖1所示。

圖1 薄圓盤直徑端點(diǎn)之間電阻測(cè)量平臺(tái)示意圖

在Comsol multiphysics軟件AC/DC模塊中通過(guò)設(shè)置薄圓盤屬性(半徑、厚度、電導(dǎo)率)、電極位置和半徑、小孔屬性(半徑、位置及電導(dǎo)率)、指定電流輸入端和輸出端建立實(shí)驗(yàn)?zāi)Ｐ?。在Comsol multiphysics軟件后處理模式下研究薄圓盤中電勢(shì)、電流分布的模擬效果，并在此基礎(chǔ)上分析不同實(shí)驗(yàn)條件下薄圓盤的電勢(shì)、電流變化，實(shí)現(xiàn)薄圓盤電阻變化規(guī)律研究。

2 薄圓盤直徑端點(diǎn)之間電阻測(cè)量方法研究

針對(duì)特定材質(zhì)和厚度的薄圓盤，影響薄圓盤直徑端點(diǎn)之間電阻的主要因素有： (1)薄圓盤半徑；(2)小孔半徑和位置；(3)電極位置。因此，為實(shí)現(xiàn)薄圓盤直徑端點(diǎn)之間電阻測(cè)量，本文分別研究了薄圓盤直徑端點(diǎn)之間電阻隨薄圓盤半徑、電極位置、薄圓盤上小孔位置和大小的變化規(guī)律，并在此基礎(chǔ)上利用最小二乘法擬合薄圓盤直徑端點(diǎn)之間電阻表達(dá)式。

2.1 薄圓盤半徑對(duì)薄圓盤直徑端點(diǎn)之間電阻的影響分析

為研究薄圓盤半徑對(duì)薄圓盤直徑端點(diǎn)之間電阻的影響，本文設(shè)置薄圓盤厚度(0.1cm)和電導(dǎo)率(15S/m)、端點(diǎn)電極的位置(如圖1所示)和半徑(0.25cm)、電極輸入電流(1A)不變，圓盤半徑a變化范圍為0≤a≤10cm。根據(jù)電極輸入電流及薄圓盤直徑端點(diǎn)之間電勢(shì)模擬結(jié)果，計(jì)算得到薄圓盤直徑端點(diǎn)之間電阻變化規(guī)律如圖2所示。

圖2 薄圓盤直徑端點(diǎn)之間電阻隨薄圓盤半徑的變化規(guī)律

從計(jì)算結(jié)果可以看出：薄圓盤直徑端點(diǎn)之間電阻隨薄圓盤半徑的增大而增大。根據(jù)計(jì)算結(jié)果，利用最小二乘法可以得到薄圓盤直徑端點(diǎn)之間電阻R與薄圓盤半徑a之間關(guān)系為

R=-1.66×10-3a6+7.63×10-2a5-1.49a4+16.03a3-106.10a2+480.70a+520.10

(1)

擬合結(jié)果擬合度為1.0(擬合度越接近1.0表示擬合結(jié)果越準(zhǔn)確)。

2.2 薄圓盤上小孔半徑及位置對(duì)薄圓盤直徑端點(diǎn)之間電阻的影響分析

1) 小孔半徑對(duì)薄圓盤直徑端點(diǎn)之間電阻的影響分析

為研究小孔半徑對(duì)薄圓盤直徑端點(diǎn)之間電阻的影響，本文設(shè)置薄圓盤尺寸(厚度為0.1cm、半徑為10cm)和電導(dǎo)率(15S/m)、端點(diǎn)電極的位置(如圖1所示)和半徑(0.25cm)、電極輸入電流(1A)、小孔位置(位于薄圓盤中心)不變，小孔半徑變化范圍為0≤r≤9cm。根據(jù)電極輸入電流及薄圓盤直徑端點(diǎn)之間電勢(shì)模擬結(jié)果，計(jì)算得到薄圓盤直徑端點(diǎn)之間電阻變化規(guī)律如圖3所示。

圖3 薄圓盤直徑端點(diǎn)之間電阻隨小孔半徑的變化規(guī)律

從計(jì)算結(jié)果可以看出：薄圓盤直徑端點(diǎn)之間電阻隨小孔半徑增大呈遞增的趨勢(shì)；小孔半徑越大，薄圓盤直徑端點(diǎn)之間電阻增大的速度越快。根據(jù)計(jì)算結(jié)果，利用最小二乘法可以得到薄圓盤直徑端點(diǎn)之間電阻R與小孔半徑r之間關(guān)系為

R=0.43r6-9.84r5+85.52r4-346.10r3+

666.20r2-472.20r+1929

(2)

擬合結(jié)果擬合度為0.9993。

2) 小孔位置對(duì)薄圓盤直徑端點(diǎn)之間電阻的影響分析

為研究小孔位置對(duì)薄圓盤直徑端點(diǎn)之間電阻的影響，本文設(shè)置薄圓盤尺寸(厚度為0.1cm、半徑為10cm)和電導(dǎo)率(15S/m)、端點(diǎn)電極的位置(如圖1所示)和半徑(0.25cm)、電極輸入電流(1A)不變，分別針對(duì)小孔位于薄圓盤x軸或y軸上、小孔位于薄圓盤上任意位置兩種情況，研究小孔位置變化對(duì)薄圓盤直徑端點(diǎn)之間電阻的影響。

(1) 小孔位于薄圓盤x軸或y軸上

設(shè)置小孔圓心分別位于x軸、y軸上，小孔半徑r滿足0.5≤r≤2.5cm。根據(jù)電極輸入電流及薄圓盤直徑端點(diǎn)之間電勢(shì)模擬結(jié)果，計(jì)算得到薄圓盤直徑端點(diǎn)之間電阻變化規(guī)律如圖4所示。

圖4 薄圓盤直徑端點(diǎn)電阻隨小孔位置的變化規(guī)律

圖5 薄圓盤等分示意圖

從計(jì)算結(jié)果可以看出： ①當(dāng)小孔位于x軸上時(shí)，小孔越靠近薄圓盤邊緣位置，薄圓盤直徑端點(diǎn)之間電阻越大；針對(duì)同一小孔位置，薄圓盤直徑端點(diǎn)之間電阻隨小孔半徑的增大而增大；小孔越靠近薄圓盤邊緣位置，電阻隨小孔半徑遞增的趨勢(shì)越明顯；②當(dāng)小孔位于y軸上時(shí)，小孔越靠近薄圓盤邊緣位置，薄圓盤直徑端點(diǎn)之間電阻越大?。会槍?duì)同一小孔位置，薄圓盤直徑端點(diǎn)之間電阻隨小孔半徑的增大而增大；小孔越靠近薄圓盤中心位置，電阻隨小孔半徑遞增的趨勢(shì)越明顯。

(2) 小孔位于薄圓盤上任意位置

為研究任意位置小孔對(duì)薄圓盤直徑端點(diǎn)之間電阻的影響，本文將薄圓盤分為4個(gè)象限，每個(gè)象限等分為8份，如圖5所示。設(shè)置小孔圓心與x軸的夾角θ(如圖1所示)變化范圍為π/32≤θ≤61π/32、小孔圓心到大圓圓心的距離ρ變化范圍為3.4cm≤ρ≤5.8cm。根據(jù)電極輸入電流及薄圓盤直徑端點(diǎn)之間電勢(shì)模擬結(jié)果，計(jì)算得到薄圓盤直徑端點(diǎn)之間電阻變化規(guī)律如圖6所示。

圖6 薄圓盤直徑端點(diǎn)電阻隨小孔位置的變化規(guī)律

從計(jì)算結(jié)果可以看出： ①小孔圓心與x坐標(biāo)軸的夾角越大，薄圓盤直徑端點(diǎn)之間電阻在Ⅰ/Ⅲ象限呈遞減的趨勢(shì)，而在Ⅱ、Ⅳ象限呈遞增趨勢(shì)。即小孔越靠近y軸，薄圓盤直徑端點(diǎn)之間電阻越?。虎诋?dāng)小孔圓心與x坐標(biāo)軸的夾角相同時(shí)，小孔距離薄圓盤圓心的距離越小，薄圓盤直徑端點(diǎn)之間電阻越小。

2.3 電極半徑對(duì)薄圓盤直徑端點(diǎn)之間電阻的影響分析

為研究電極半徑對(duì)薄圓盤直徑端點(diǎn)之間電阻的影響，本文設(shè)置薄圓盤尺寸(厚度為0.1cm、半徑為10cm)和電導(dǎo)率(15S/m)、端點(diǎn)電極的位置(如圖1所示)、電極輸入電流(1A)、小孔位置(位于薄圓盤中心)和半徑(1.5cm)不變，電極半徑變化范圍為0.1≤re≤0.25cm。根據(jù)電極輸入電流及薄圓盤直徑端點(diǎn)之間電勢(shì)模擬結(jié)果，計(jì)算得到薄圓盤直徑端點(diǎn)之間電阻變化規(guī)律如圖7所示。

圖7 薄圓盤直徑端點(diǎn)電阻隨電極半徑的變化規(guī)律

從計(jì)算結(jié)果可以看出：薄圓盤直徑端點(diǎn)之間電阻隨電極半徑的減小而增大。根據(jù)計(jì)算結(jié)果，利用最小二乘法可以得到薄圓盤直徑端點(diǎn)之間電阻R與小孔半徑re之間關(guān)系為

(3)

擬合結(jié)果擬合度為1.0。

3 驗(yàn)證與分析

為驗(yàn)證上文所述薄圓盤直徑端點(diǎn)之間電阻的測(cè)量方法，本文在特定條件下對(duì)上述測(cè)量方法進(jìn)行驗(yàn)證：設(shè)置薄圓盤厚度為0.1cm、半徑為10cm、電導(dǎo)率為15S/m，電極輸入電流(1A)，小孔位于圓心、半徑為4.5cm，端點(diǎn)電極半徑為0.25cm，電極內(nèi)切于薄圓盤、外切于小孔(如圖8(a)所示)。此時(shí)薄圓盤直徑端點(diǎn)之間的電阻R相當(dāng)于兩個(gè)并聯(lián)的薄長(zhǎng)方體的電阻(如圖8(b)所示)，利用上述方法可計(jì)算得到直徑端點(diǎn)之間的阻值R為19898.61847Ω。

圖8 特定條件下電極位置示意圖(a) 電極位置示意圖; (b) 薄圓盤等效并聯(lián)電阻示意圖

(4)

則上文所述薄圓盤直徑端點(diǎn)之間電阻測(cè)量方法相對(duì)誤差為

(5)

由計(jì)算結(jié)果可以看出，本文提出的薄圓盤直徑端點(diǎn)之間電阻測(cè)量方法的誤差較小。

4 結(jié)語(yǔ)

本文借助Comsolmultiphysics軟件，利用有限元分析法實(shí)現(xiàn)了薄圓盤直徑端點(diǎn)之間電阻的測(cè)量，并在此基礎(chǔ)上研究了不同因素對(duì)薄圓盤直徑端點(diǎn)之間電阻的影響。研究結(jié)果表明：①本文提出的薄圓盤直徑端點(diǎn)之間電阻測(cè)量方法的誤差較??；②當(dāng)電極參數(shù)固定時(shí)，薄圓盤直徑端點(diǎn)之間電阻隨薄圓盤半徑的增大而增大；③當(dāng)薄圓盤參數(shù)固定時(shí)，薄圓盤直徑端點(diǎn)之間電阻隨電極半徑增大而減小，隨電極與薄圓盤圓心之間距離的增大而減??；④薄圓盤上小孔越靠近y軸，薄圓盤直徑端點(diǎn)之間電阻越小。本文研究工作對(duì)薄圓盤電阻的設(shè)計(jì)具有指導(dǎo)意義，且該研究方法還可以推廣用于其他形狀電阻的設(shè)計(jì)。

[1] 李靜.有限元方法在電磁場(chǎng)分析中的研究[D]. 重慶: 重慶師范大學(xué), 2011. Li Jing. The Finite element method be used in electromagnetic field analysis[D]. Chongqing: Chongqing Normal University, 2011. (in Chinese)

[2] 李為峰. 接地網(wǎng)電阻檢測(cè)方法的研究[D]. 西安: 西安理工大學(xué), 2009. Li Weifeng. The Research on resistance measuring method for grounding grids[D]. Xi’an: Xi’an University of Technology, 2009. (in Chinese)

[3] Miroslav ?imandl,Jaroslav ?vácha. Separation approach for numerical solution of the fokker-planck equation in estimation problem[J]. IFAC Proceedings Volumes, 2005, 38(1).

[4] 倪興. 常微分方程數(shù)值解法及其應(yīng)用[D]. 北京: 中國(guó)科學(xué)技術(shù)大學(xué), 2010. Ni Xing. Numerical solution of ordinary differential equations and its application[D]. Beijing: University of Science and Technology of China, 2010. (in Chinese)

[5] 侯雪.基于comsol的肺部電阻抗斷層成像仿真研究[D]. 天津:天津大學(xué), 2011. Hou Xue. Lung electrical impedance tomography simulation based on COMSOL[D]. Tianjin: Tianjin University, 2011. (in Chinese)

[6] 李忠良.基于壓電阻抗健康監(jiān)測(cè)的實(shí)驗(yàn)及有限元仿真分析[D]. 浙江:浙江大學(xué),2013. Li Zhongliang. Experimental and numerical simulation of health monitoring based on EMI[D]. Zhejiang: Zhejiang University, 2013. (in Chinese)

[7] 毛亞琴.電阻法研究碳纖維樹脂基單向拉擠復(fù)合材料的破壞行為[D]. 北京:北京化工大學(xué),2004. Mao Yaqin. Study on the failure behavior of cfrp unidirectional pultruded composites by electrical resistance method[D].Beijing: Beijing University of Chemical Technology, 2004. (in Chinese)

[8] 陳小燕.基于COMSOL軟件的大地電磁正演模型[D]. 西安:長(zhǎng)安大學(xué),2014. Chen Xiaoyan. Magnetotelluric forward modeling based on software of COMSOL[D]. Xi’an: Chang’an University, 2014. (in Chinese)

[9] Feng Zhongkui, Xu Ying, Li Xian, et al. Optimum electromagnetism design based on co-simulation of Matlab and COMSOL[J]. Electromagnetic Analysis and Applications, 2012, 01(01): 13-17.

[10] 許燦榮.基于COMSOL的近距離下保護(hù)層開采瓦斯流動(dòng)模型[D]. 河南:河南理工大學(xué)，2012. Xu Canrong. Methane flow model of the near distance lower protective coal seam based on COMSOL[D]. Henan: Henan Polytechnic University, 2012. (in Chinese)

■

STUDY ON THE RESISTANCE BETWEEN ENDPOINTS IN DIAMETER OF A THIN DISC CONDUCTOR

Wang Lan1Xi Xiaoqiang2

(1School of Science；2Institute of Internet of Things and IT-based Industrialization, Xi’an University of Posts and Telecommunications, Xi’an Shaanxi 710061)

It was difficult to apply immediate integral to measure the resistance of thin disc conductor. In order to solve this problem, this paper studied the method for measuring the resistance of thin disc by finite element method. Firstly, the measurement platform of thin disc resistance was designed by Comsol multiphysics software. Based on this design of electrode, we studied the relationship between the thin disc resistance and the disc radius, the electrode position, the position and size of the small hole on the thin disc. Then, resistance of thin disc could be measured indirectly by the disc resistance expression, which was fitted by the least square method. Finally, the above measurement method for the resistance of thin disc was verified with immediate integral under the typical condition, and the results show that the error of the method was quite small. This paper was of great significance for the design of resistance of thin disc. In addition, the resistance measurement method could also be promoted the design of resistance of other shapes.

thin disc conductor; resistance; Comsol multiphysics

2016-04-27；

2017-03-10

國(guó)家自然科學(xué)基金(11275099;11475135)。

王嵐，女，碩士研究生，研究方向?yàn)槲锢矸椒ㄔ谛畔⑻幚碇械膽?yīng)用，563859805@qq.com。

王嵐，惠小強(qiáng). 薄圓盤直徑端點(diǎn)之間電阻的研究[J]. 物理與工程，2017，27(4)：37-41.