Ⅰ. 서 론
다중도체 전송선로(multiconductor transmission lines, MTLs)는 n+1개의 평행 도체로 구성된 구조를 의미하며, 전기 및 전자 시스템에서 신호를 전송하거나 전력을 공급하는 데 사용된다. 이러한 다중도체 전송선로는 고주파 및 고속회로 설계에 있어 중요한 역할을 하고 있으며, 이를 분석하기 위해 많은 연구가 이루어졌다[1]~[5].
다중도체 전송선로의 특성은 전파상수(propagation constants)와 특성 임피던스(characteristic impedances)로 결정되며, 이러한 매개변수는 impedance 행렬로부터 얻을 수 있는 resistance 및 inductance 행렬로 계산된다[5]. 또한, 신호지연(signal delay)이나 crosstalk와 같은 현상을 예측하고 제어하기 위해서는 capacitance 행렬을 추출하는 것은 필수적이다[6]~[8].
다중도체 선로와 같은 전자기 구조의 해석을 위해서는 Maxwell 방정식에 대한 해를 구해야 한다. 하지만, 이 방정식의 참 해(exact solution)를 구하는 것은 매우 복잡하다. 따라서, 컴퓨터를 통해 Maxwell 방정식에 대한 수치 해(numerical solution)를 구하는 전자기 수치해석(computational electromagnetics, CEM)의 도입이 필수적으로 필요하다[9]. 이러한 전자기 수치해석 기법의 하나인 유한요소법(finite element method, FEM)은 곡면과 같은 기하학적으로 복잡한 구조의 해석이 가능하며, 비균질성 매질의 해석에 강점을 가진다. 게다가, 유한요소 과정을 통해 얻게 되는 행렬은 0을 많이 포함하고 있는 sparse matrix가 되기 때문에 계산 부하가 적으며, 반복법 및 분할 알고리듬을 적용한 병렬 연산을 적용하기에 적합하다는 장점이 있다[10],[11].
유한요소법을 통한 다중도체 전송선로 해석에는 full wave 해석과 준정적 해석(quasistatic analysis 또는 quasi-TEM analysis)이 존재한다. Full wave 해석의 경우 정확한 결과를 얻을 수 있지만, eigenvalue problem을 주파수마다 풀어야 하므로 상당한 계산을 요구한다[3],[4]. 반면에 준정적 해석의 경우 full wave 해석보다 적은 미지수로 계산이 가능하며, eigenvalue problem이 아닌 matrix equation을 풀기 때문에 계산이 비교적 간단하며, full wave 해석에서 발생하는 저주파 불안정성 문제를 피할 수 있다는 장점이 있다[3],[4]. 하지만 준정적 해석의 경우 비균질 매질에서 주파수가 증가함에 따라 정확도가 감소한다는 단점이 있어, 두 방법을 결합하여 해를 구하기도 한다[3],[4].
본 논문에서는 100 % 직접 작성한 유한요소 코드를 사용하여 2-D shielded stripline 및 3-D interconnect 구조의 capacitance 행렬을 추출하는 시뮬레이션을 수행하였다. 이와 함께, 2-D 및 3-D 정전기(electrostatic) 문제 해결을 위한 유한요소법 적용 과정과 후처리 과정을 통해 capacitance 행렬을 추출하는 과정에 대해 자세히 다루었다. 아울러, 유한요소법을 이용한 다중도체 전송선로 구조 해석에 있어 full wave 해석 대신 준정적 해석을 통해 다중도체 및 microstrip line 구조에 대한 per-unit-length resistance와 inductance 행렬을 계산하는 방법에 대해 논의하고, 이를 이용하여 직접 시뮬레이션을 수행하여 자체 유한요소 코드의 정확성을 확인하였다.
본 논문에서 구현된 기법은 국내 전파해석 소프트웨어 개발에 있어 해석 시나리오 및 포트폴리오를 추가하는 데 상당한 기여를 할 것으로 사료된다. 추가로, 본 논문은 학술대회에서 발표된 참고문헌[12],[13]을 수정 보완한 확장판임을 밝힌다.
Ⅱ. 2-D 유한요소해석
본 절에서는 2-D 유한요소법을 적용하여 그림 1과 같은 shielded stripline 구조에 대한 전압분포 해석과 후처리 과정을 통해 capacitance 행렬을 추출하는 시뮬레이션을 수행하였다.
![jkiees-35-1-7-g1](/journal/jkiees/jkiees-35-1/gif/jkiees-35-1-7-g1.gif)
2-D 유한요소해석에서 경곗값 문제(boundary-value problem)는 식 (1)로 정의되며, 그림 1과 같은 정전기 문제에서 전압(V)은 식 (2)와 같은 Poisson 방정식의 지배를 받는다[5],[12],[15].
유한요소 적용을 위해 식 (2)의 형태를 식 (1) 형태의 일반적인 2계 미분 방정식과 맞춰줘야 하며, 그 결과 식 (3)과 같은 관계를 도출할 수 있다[5],[12],[15].
유한요소 과정을 통해 전압을 구하면 식 (4)를 통해 저장된 에너지(W)를 구할 수 있으며, 이를 통해 coefficients of capacitance() 및 induction(,i≠j)을 계산할 수 있다[8],[12],[15]. 아울러, 식 (4)는 유한요소 적용을 통해 식 (5)와 같이 계산된다[12],[15].
Capacitance 행렬의 모든 요소를 구하기 위해서는 해석 구조의 도체의 수를 Nc개라고 했을 때, 도체 간의 상호작용을 고려하여 총 Nc+Nc(Nc−1)/2번의 계산이 필요하며, 계산과정에서의 경계조건은 식 (6)과 같다[8],[12],[15].
Capacitance 행렬의 원소가 되는 및 는 식 (7)을 통해 얻을 수 있으며, per-unit-length capacitance는 식 (8)을 통해 계산된다[1],[8],[15].
본 절에서는 2-D 유한요소법과 준정적 해석을 적용하여 그림 2 및 그림 3과 같은 다중도체 전송선로 구조의 per-unit-length resistance 및 inductance 행렬을 추출하는 문제에 대해 다룬다.
손실 도체들로 구성된 전송선로의 경우 변위전류(displacement current)는 무시하고 전도전류(conduction current)만 고려하는 준정적 해석을 적용할 수 있다[3]~[5],[16]. 준정적 해석을 적용하면 도체 내에서의 transverse electric field는 무시할 수 있다고 가정하며, 이에 따라 도체에 z 방향의 전기장을 인가할 경우 Maxwell 방정식으로부터 식 (9)를 얻을 수 있다[3]~[5],[16].
식 (9)에서 σ는 손실 도체의 도전율을, Jimp 는 인가된 전류를 의미한다. 이때, Jimp는 인가된 전기장(Eimp)에 의해 z 방향을 따라 유도된다. 따라서 식 (9)를 식 (10)과 같이 표현할 수 있으며, 유한요소 적용을 위해 식 (1)과 식 (10)의 형태를 맞추면 식 (11)과 같은 관계 식을 얻을 수 있다[3]~[5].
그림 2 및 그림 3과 같은 개방 전송선로의 경우 potential Az가 전송선로에 집중되기 때문에 해당 구조를 매우 큰 상자로 감싼 뒤, homogeneous Dirichlet 또는 Neumann 경계조건을 적용하여 해석할 수 있다[3]~[5].
Per-unit-length resistance 및 inductance 행렬은 per-unit -length impedance 행렬 [Z]=[R]+jω[L]을 통해 얻을 수 있다[3]~[5]. 이러한 행렬을 계산하기 위해서는 i번째 도체에 Eimp=1을 인가하고 유한요소 과정을 통해 Az를 계산한 뒤, 식 (12)를 이용하여 Ik를 계산해야 한다[3]~[5].
도체의 개수만큼 Eimp=1을 인가하여 유한요소 과정을 반복함으로써 식 [Z]=[V][I]−1를 통해 per-unit-length impedance 행렬을 얻을 수 있다[3]~[5]. 이때, [V]는 Eimp=1로 이루어진 대각행렬이며, [I]는 i번째 열이 식 (12)를 통해 계산된 전류로 이루어진 정방행렬이다.
Ⅲ. 3-D 유한요소 해석
3-D 유한요소 해석에서는 그림 4와 같은 PEC box 안에 있는 3차원 interconnect 구조에서의 전압분포를 해석하고, 이를 이용하여 capacitance 행렬을 추출하였다.
그림 4와 같은 3차원 구조에서 전압은 식 (13)과 같은 일반화된 Poisson 방정식(generalized Poisson’s equation)의 지배를 받는다[5],[13],[15].
![jkiees-35-1-7-g4](/journal/jkiees/jkiees-35-1/gif/jkiees-35-1-7-g4.gif)
유한요소 해석을 위해서는 식 (13)을 3차원 경곗값 문제의 일반적인 형태인 식 (14)와 형태를 맞춰줘야 하며, 그 결과 식 (15)와 같은 관계 식을 얻을 수 있다[5],[13],[15].
아울러, capacitance 행렬 계산은 Net에 서로 다른 전압 인가를 고려하여 2차원과 같게 총 Nc+Nc(Nc-1)/2번의 계산이 필요하며, 경계조건 또한 식 (6)과 같다[8],[12],[13],[15].
이때, 3차원 구조의 경우 저장된 에너지는 식 (16)을 통해 계산할 수 있다[8],[13],[15]. 여기서, ve는 요소의 부피를 의미한다.
이후에, 저장된 에너지를 통해 계산되는 및 , 그리고 per-unit-length capacitance의 계산과정은 2차원(2-1-1절)과 같다[8],[13],[15].
Ⅳ. 시뮬레이션 결과
본 장에서는 in-house FEM code와 참고문헌[5] 및 Q3D[17]와의 결과 비교를 통해 해석의 정확도와 유효성을 검증하였다. 상용 소프트웨어의 경우 어떤 구체적 알고리듬을 사용하는지 공개되어 있지 않아 in-house FEM code와 엄정한 비교는 현실적으로 쉽지 않다. 따라서 본 논문에서는 상용 소프트웨어와 최대한 같은 상황을 모사하고 이에 대한 정확도와 실행시간 비교 중 정확도에 초점을 맞춰 시뮬레이션을 수행하였다. 또한, 본 장에서 다루는 시뮬레이션의 솔버는 MATLAB에서 built-in으로 제공하는 direct 솔버인 backslash 연산을 선택하여 사용하였다. 채택의 주된 이유는 최종 시스템 행렬의 크기가 아주 크지 않으며, 손실 도체 및 비균일 매질로 인한 condition number의 증가로 iterative 솔버의 수렴속도가 감소하는 현상을 회피하기 위함이다. 추가적으로, 본 장의 2-D 시뮬레이션에서는 node-based triangular basis, 3-D 시뮬레이션에서는 node-based brick basis를 사용하여 해석을 진행하였다.
2-D 노드 기반 유한요소법을 사용한 shielded stripline 구조(그림 1) 해석은 직접 작성한 in-house code(MATLAB R2022b 버전)를 이용하여 진행되었으며, 전자기 상용 소프트웨어 Q3D(2D extractor)[17]와의 결과 비교를 통해 시뮬레이션의 정확성을 검증하였다. 아울러, 삼각형 모양의 요소를 사용하여 해석 영역을 나누었으며, 총 15,600개의 요소와 7,991개의 노드가 사용되었다.
그림 5는 Net 2와 Net 3에 전압을 인가하였을 때의 전압분포 시뮬레이션 결과이며, 후처리 과정을 통해 계산한 capacitance 행렬과 Q3D[17]와의 결과 비교를 표 1에 정리하였다.
In-house | 167.524 | 173.637 | 167.524 |
Q3D[17] (2D extractor) |
165.514 | 171.374 | 165.628 |
Rel. error (%) | 1.214 | 1.321 | 1.145 |
In-house | −29.19 | −0.435 | −29.19 |
Q3D[17] (2D extractor) |
−28.38 | −0.431 | −28.219 |
Rel. error (%) | 2.854 | 0.928 | 3.441 |
다중도체 구조(그림 2) 및 coupled microstrip line(그림 3)의 시뮬레이션은 MATLAB R2022b 버전을 기반으로 실행되었으며, 참고문헌[3]~[5]과의 결과 비교를 통해 시뮬레이션의 유효성을 검증하였다.
그림 2 구조의 경우 총 107,800개의 삼각형 요소와 54,426개의 노드를 사용하였으며, 도체의 도전율은 aluminum-oxide와 같은 σ=3.6×107 S/m로 설정하였다. 또한, 그림 2 구조를 x 방향으로 33 um, y 방향으로 12 um의 크기를 가진 큰 상자로 감싼 뒤 상자 경계에 PEC 경계조건을 적용하여 계산영역을 한정시켰다.
도체 수가 4개인 그림 2 구조는 식 (17)과 같은 telegrapher equation을 만족한다[4]. 따라서 본 시뮬레이션에서는 4×4 크기의 per-unit-length impedance 행렬을 구한 뒤, 4번 도체를 reference conductor로 설정하고, 식 I1+I2+I3+I4=0 을 이용하여 식 (18)과 같은 3×3 크기의 Ztrue행렬을 통해 per-unit-length resistance 및 inductance 행렬을 계산하였다[4].
그림 6 및 그림 7은 그림 2 구조에 대한 주파수에 따른 per-unit-length resistance 및 inductance 행렬 시뮬레이션 결과를 나타내며, 참고문헌 [5]의 Figure 4.12 (b) 및 Figure 4.12 (c)의 결과와 잘 일치함이 확인된다.
![jkiees-35-1-7-g6](/journal/jkiees/jkiees-35-1/gif/jkiees-35-1-7-g6.gif)
![jkiees-35-1-7-g7](/journal/jkiees/jkiees-35-1/gif/jkiees-35-1-7-g7.gif)
그림 3 구조는 2차원 노드 기반 삼각형 요소를 사용하여 총 126,360개의 요소와 63,685개의 노드를 사용하였으며, 도체의 도전율은 각각 σcu=5.8×107 S/m와 σcr=7.75× 106 S/m로 설정하였다. 아울러, 그림 3 구조를 감싼 PEC 상자의 크기는 x 방향으로 135 um, y 방향으로 117 um이다. 그림 8은 그림 3에 대한 per-unit-length resistance 및 inductance 행렬 시뮬레이션 결과를 나타낸다. In-house code의 결과와 참고문헌 [5]의 Figure 4.13 (b) 결과가 비교적 잘 일치하지만, 0 근처의 log-scale 그래프의 특성으로 인해 다소간의 오차가 관찰되는 것을 확인할 수 있다.
![jkiees-35-1-7-g8](/journal/jkiees/jkiees-35-1/gif/jkiees-35-1-7-g8.gif)
3-D 노드 기반 유한요소법을 사용한 interconnect 구조(그림4) 해석은 MATLAB R2022b 환경에서 실행되었으며, Q3D[17]와의 결과 비교를 통해 시뮬레이션의 정확성을 검증하였다. interconnect 구조는 크기가 6×5×30 um인 PEC 상자 내부에 들어가 있으며, 해석을 위해 총 72,000개의 brick 요소와 77,469개의 노드가 사용되었다.
그림 9는 Net 1과 Net 2에 전압을 인가하였을 때, xy 평면(z=15 um)에서의 전압분포를 나타낸다. 또한, 총 6번의 유한요소 과정을 통해 capacitance 행렬을 계산하였으며, Q3D[17]와의 결과 비교를 표 2에 정리하였다.
In-house | 1.92 | 0.736 | 0.736 |
Q3D[17] | 1.863 | 0.704 | 0.704 |
Error (%) | 3.06 | 4.546 | 4.546 |
In-house | −0.199 | −0.199 | −0.282 |
Q3D[17] | −0.193 | −0.193 | −0.266 |
Error (%) | 3.109 | 3.109 | 6.015 |
Ⅴ. 결 론
본 논문에서는 2-D 및 3-D 정전기 문제에 있어, 전압분포 해석과 capacitance 행렬 추출이 가능한 유한요소 자체 코드를 확보하였다. 이 과정에서, 대표적인 전자기해석 시뮬레이션 소프트웨어 중 하나인 Q3D를 사용한 시뮬레이션 결과와 직접 작성한 in-house 코드의 결과를 비교하여 자체 코드의 유효성을 검증하였다. 또한, 손실 도체들로 구성된 전송선로 구조를 간단하게 해석하기 위해서 준정적 해석을 적용하였으며, per-unit-length impedance 행렬을 통해 per-unit-length resistance 및 inductance 행렬을 구하고, 이를 참고문헌과 비교하여 시뮬레이션 결과의 정확성을 확인하였다.
Capacitance 행렬을 구하기 위해선 총 Nc+Nc(Nc-1)/2번, per-unit-length resistance 및 inductance 행렬을 구하기 위해선 총 Nc번의 유한요소 계산이 필요하게 된다. 따라서 Nc의 수가 증가할수록 계산량이 증가하며, 이에 따라 해석시간 또한 증가하게 될 것이다. 그뿐만 아니라, 개방 전송선로 해석을 위해 큰 PEC 상자를 설정하는 과정에서, 계산영역의 증가로 인해 해석시간이 더욱 증가할 수 있다. 이러한 해석시간의 증가는 추후에 domain decomposition, CUDA를 이용한 GPU 병렬처리 등의 병렬화 기법을 통해 더 나은 속도 향상을 기대해 볼 수 있으며, 개방 전송선로의 경우 PML(perfectly matched layer) 등의 적용으로 계산영역을 한정시켜 속도 향상을 도모할 수 있을 것으로 예상한다.