KR101616036B1 - 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델 - Google Patents

인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델 Download PDF

Info

Publication number
KR101616036B1
KR101616036B1 KR1020157032574A KR20157032574A KR101616036B1 KR 101616036 B1 KR101616036 B1 KR 101616036B1 KR 1020157032574 A KR1020157032574 A KR 1020157032574A KR 20157032574 A KR20157032574 A KR 20157032574A KR 101616036 B1 KR101616036 B1 KR 101616036B1
Authority
KR
South Korea
Prior art keywords
inductor
idc
current
voltage
direct current
Prior art date
Application number
KR1020157032574A
Other languages
English (en)
Other versions
KR20150138396A (ko
Inventor
세이지 히다카
코 야마나가
Original Assignee
가부시키가이샤 무라타 세이사쿠쇼
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 가부시키가이샤 무라타 세이사쿠쇼 filed Critical 가부시키가이샤 무라타 세이사쿠쇼
Publication of KR20150138396A publication Critical patent/KR20150138396A/ko
Application granted granted Critical
Publication of KR101616036B1 publication Critical patent/KR101616036B1/ko

Links

Images

Classifications

    • G06F17/5036
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/5045

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Filters And Equalizers (AREA)
  • Networks Using Active Elements (AREA)
  • Measurement Of Resistance Or Impedance (AREA)
  • Semiconductor Integrated Circuits (AREA)

Abstract

직류 중첩 시의 비선형 특성을 정밀도 높게 동적으로 시뮬레이션 가능한 인덕터의 시뮬레이션 방법 및 비선형 등가 회로 모델을 단순한 구성으로 간단하게 제공한다. 그렇기 때문에 본 발명은 인덕터의 등가 회로를 수동 회로 소자(L1, R1)의 직렬 회로를 사용하여 나타낸다. 또한, 직류 전류(Idc)의 중첩 시에서의 수동 회로 소자(L1, R1)의 특성 변화율(kL1(Idc), kR1(Idc))을, 실측값을 기초로, 근사 함수 exp(f(x))로서 나타낸다. 수동 회로 소자(L1, R1)에 직렬로 접속되는 제어 전압원(B1)에 의해 전압원 모델(V0, V1)에 따라 측정되는 참조 전류(Iref)가 참조되고, 참조 전류(Iref)에 대응하여 특성 변화율(kL1(Idc), kR1(Idc))이 산출된다. 그리고 특성 변화율(kL1(Idc), kR1(Idc)) 및 무 중첩시 전압(VL1, VR1)에 기초하여, 차분 전압(ΔVL1(Idc), ΔVR1(Idc))이 발생되고, 무 중첩시 전압(VL1, VR1)에 중첩됨으로써 비선형 특성이 시뮬레이트된다.

Description

인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델{INDUCTOR SIMULATION METHOD AND NONLINEAR EQUIVALENT CIRCUIT MODEL FOR INDUCTOR}
본 발명은 인덕터의 직류 전류 중첩 시의 비선형 특성을 시뮬레이트하는 인덕터의 시뮬레이션 방법 및 그 시뮬레이트를 할 때에 이용되는 인덕터의 비선형 등가 회로 모델에 관한 것이다.
종래, 이 종류의 전자부품의 시뮬레이션 방법 및 등가 회로 모델은 전자 회로 설계에서의 회로 시뮬레이트 시에 이용되고 있다. 회로 시뮬레이트에는, SPICE(Simulation Program with Integrated Circuit Empahsis) 등의 회로 시뮬레이터가 이용되고, 회로 시뮬레이터 중에는, 전자부품 제조 메이커의 홈페이지 사이트에서 이용할 수 있는 것도 있다. 이용자는 퍼스널 컴퓨터(personal computer) 등의 단말로부터 인터넷망을 통해서 전자부품 제조 메이커의 홈페이지 사이트에 액세스하여, 회로 시뮬레이터를 이용한다.
종래, 이 종류의 시뮬레이션 방법 및 등가 회로 모델로서는, 예를 들면, 특허문헌 1에 개시된 콘덴서에 대한 것이 있다.
이 시뮬레이션에서는, 동 문헌의 도 1에 표시되는 바와 같이, 제1 스텝에 있어서, 콘덴서가 주어진 주파수 특성을 입력하고, 제2 스텝에 있어서, 주파수에 의존하지 않는 저항(R), 커패시턴스(C), 인덕턴스(L)를 이용한, 시간 영역에서의 시뮬레이션이 가능한 회로로서, RC회로 및 RL회로와, RCL회로 중 어느 하나를 등가 회로 모델로서 형성한다. 그리고 제3 스텝에 있어서, 제2 스텝에서 형성한 등가 회로 모델의 정밀도를 판정하기 위한 평가 함수를 합성하고, 제4 스텝에 있어서, 제3 스텝에서 합성한 평가 함수를 최소화함으로써 회로 정수를 결정한다.
특허문헌 1에서는, 상기의 구성에 의해, 주파수 영역에서 임피던스를 나타낸 콘덴서의, 시간 영역에서의 시뮬레이션이 가능한 등가 회로 모델을 도출하고, 콘덴서의 주파수 영역 또는 시간 영역에서의 전기적 특성을, 회로 시뮬레이션에 의해 예측한다.
또한, 종래, 인덕터의 시뮬레이션 방법 및 등가 회로 모델로서, 예를 들면, 특허문헌 2에 개시된 것이 있다.
이 시뮬레이션에서는, 동 문헌의 도 1(C)에 나타내는 바와 같이, 내부 도체의 표피 효과를 고려하는 인덕턴스(L1)와 레지스탕스(resistance)(R1)의 직렬 회로에, 직류에 대한 인덕턴스(L0)와 인덕턴스(L1) 사이의 상호 인덕턴스(Lm)가 병렬 접속되고, 이것에 직류에 대한 인덕턴스(L0)와 내부 도체의 직류 저항(Rdc1)이 직렬로 접속된 등가 회로 모델이 사용된다. 이 등가 회로 모델에서는, 더욱 외부전극의 인덕턴스 및 레지스탕스가 동시에 고려되어, 인덕턴스(L0)에 직렬로 외부전극의 인덕턴스(Ls)가 직렬로 접속됨과 함께, 내부 도체의 직류 저항(Rdc1)에 외부전극의 직류 저항(Rdc2)이 직렬로 접속된다. 또한, 적층 칩 인덕터의 칩을 구성하는 유전체의 기생 커패시턴스(Cp)와, 유전체의 손실을 나타내는 저항(Rp)을 직렬로 접속한 직렬 회로가 외부전극의 등가소자(Ls, Rdc2)의 내측에 병렬로 접속된다.
특허문헌 2에서는, 상기의 등가 회로 모델을 사용한 회로 시뮬레이션에 의해, 회로 설계와 실제의 회로 성능 사이에 발생하는 오차를 억제하고 있다.
또한, 종래, 콘덴서의 시뮬레이션 방법 및 등가 회로 모델로서, 예를 들면, 특허문헌 3에 개시된 이상(理想)(C) 회로 모델 및 광대역 고정밀도 등가 회로 모델을 이용한 것이 있다.
이상(C) 회로 모델은, 동 문헌의 도 1(A)에 나타내는 바와 같이, 하나의 커패시턴스 엘리먼트(capacitive element)(C)를 회로 소자로 하는 등가 회로로 표시된다. 커패시턴스 엘리먼트(C)의 양단(兩端)에 인가되는 전압(v)은, 그 양단에 인가된 시간 변화하는 신호 전압과 노이즈 전압을 vac, 양단에 인가된 DC 바이어스 전압을 Vdc이라고 한 경우, 다음의 (1)식으로 표시된다.
Figure 112015110688576-pct00001
커패시턴스 엘리먼트(C)가 DC 바이어스 전압(Vdc)에 의해 변화되는 특성은 다음 (2)식으로 표시되는 다항식으로 표현된다.
Figure 112015110688576-pct00002
또한, 커패시턴스 엘리먼트(C)에 흐르는 전류(i)는 다음 (3)식으로 표시된다.
Figure 112015110688576-pct00003
이 (3)식을 계산하기 위해, 동 문헌의 도 1(B)에 나타내는 바와 같이, 연산회로가 구축된다. 이 연산회로에서는 커패시턴스 엘리먼트(C)가 DC 바이어스 전압(Vdc)에 의해 컨트롤되는 논 리니어(nonlinear) 전압 제어 전압원(UA3)으로 변환된다. 또한, 콘덴서의 양단에 인가된 총 전압(v)이 리니어(linear) 전압 제어 전압원(E1)을 경유하여, 매우 낮은 컷오프(cutoff) 주파수를 가지는 로우 패스 필터(low pass filters)(L1 및 R1)를 통과함으로써 DC 바이어스 전압(VdC)이 얻어져서, 논 리니어 전압 제어 전압원(UA3)에 주어진다. 또한, 총 전압(v)이 리니어 전압 제어 전압원(E2)을 통해서 미분 디바이스(UA1)의 입력 단자에 공급됨으로써 미분(dv/dt)이 이루어진다. 또한, 미분 디바이스(UA1)의 출력 전압(v1)이, 커패시턴스 엘리먼트(C)를 대체하는 논 리니어 전압 제어 전압원(UA3)의 출력 전압(C(Vdc))과 함께, 3단자 승산 디바이스(multiplier)(UA2)에 입력됨으로써 승산(C(Vdc)·dv/dt)이 이루어진다. 이에 따라, 승산 디바이스(UA2)의 출력 단자에 승산 결과가 출력된다. 승산 디바이스(UA2)의 출력 전압(v2)은, 콘덴서에 흐르고 있는 전류(i)와 단위 저항의 곱과 동일하므로, 출력 전압(v2)에 의해 제어되는 리니어 전압 제어 전류원(G)을 이용하여 콘덴서와 치환된다.
이러한 이상(C) 회로 모델은 실제 부품의 임피던스 특성과의 차이, 특히 고주파 대역에서의 차이가 지나치게 크기 때문에 회로 시뮬레이션에는 적절하지 않지만, 회로 설계의 초기 단계 혹은 회로 특성의 예측에는 편리하다.
또한, 특허문헌 3에 개시된 광대역 고정밀도 등가 회로 모델은 MLCC(적층 세라믹 콘덴서)의 시뮬레이션에 적용된다. 이 시뮬레이션에서는 동 문헌의 도 5(A)에 나타내는 회로 구성의 등가 회로 모델이 사용된다. 동 문헌의 도 5(B)에 나타내는 바와 같이, 적층 칩 콘덴서(10)에서는 복수의 내부전극(20)이 적층되어 있고, 번갈아 전극의 인출이 이루어지고 있다. 동 문헌의 도 5(A)에 나타내는 등가 회로는 적층 칩 콘덴서(10)의 복수의 내부전극(20)의 두께를 고려한 것이며, 복수의 내부전극(20)의 각각 상부면(22) 및 하부면(24)에서의 전자(電磁) 효과 외에 복수의 내부전극(20)의 한쪽 측면(26) 및 다른 쪽의 측면(28)과, 오픈 단면(30)의 전자 효과도 고려한 것이다.
이 등가 회로에서의 각종 회로 소자의 값은 모두 DC 바이어스 전압에 의해 변화한다. DC 바이어스 전압에 의한 각 회로 소자의 특성 변화는 다항식으로 표현되고, 이 특성 변화를 예상한 경우의 MLCC의 등가 회로 모델은 동 문헌의 도 12에 나타낸다. 이 모델에서는 미분 디바이스, 승산 디바이스, 3단자나 4단자의 가산 디바이스 등 외에, 제산 디바이스와 5단자의 가산 디바이스도 사용된다. DC 바이어스 전압에 의한 특성 변화를 예상한 이러한 광대역 고정밀도 모델은 넓은 주파수 대역에 있어서 양호한 시뮬레이션 정밀도를 얻을 수 있다.
공개특허공보 2002-259482호 공개특허공보 2010-204869호 공개특허공보 2012-150579호
콘덴서나 인덕터와 같은 전자부품은 그 정전용량이나 인덕턴스와 같은 특성값이, 중첩 인가되는 DC 바이어스 전압이나 직류 전류에 의해 변화되고, 그 변화는 회로 시뮬레이션에 무시할 수 없는 것이 된다. 그러나 특허문헌 1이나 특허문헌 2에 개시된 상기 종래의 전자부품의 시뮬레이션 방법에 이용되는 등가 회로 모델은, 모두 저항 소자(R), 유도 소자(L) 및 용량 소자(C)의 수동 회로 소자만의 조합에 의해 구성되어 있고, DC 바이어스 전압 또는 직류 전류가 중첩 인가됐을 때에 전자부품에 생기는 특성 변화가 반영되지 않는다.
또한, 특허문헌 3에 개시된 상기 종래의 전자부품의 시뮬레이션 방법에서는, 전압원 모델이나 전류원 모델이 등가 회로 모델에 이용됨으로써 DC 바이어스 전압이 중첩 인가되었을 때에 전자부품에 생기는 특성 변화가 반영되어 있다. 그러나 이 등가 회로 모델은 미분, 승산, 가산 등의 연산 회로를 포함한 복잡한 구성을 하고 있기 때문에, 어떤 규칙에 준하여 규칙적으로 등가 회로 모델을 도출할 수 없어서 등가 회로 모델의 도출 작업은 매우 어렵다. 또한, 이러한 복잡한 구성을 하고 있기 때문에 회로의 동작을 읽는 것도 어렵고, 회로를 예측하기 어렵다. 게다가, (2)식에 나타낸 바와 같이, 직류 바이어스의 의존성을 반영하는 근사식에 홀수차의 거듭제곱항이 포함되기 때문에, 직류 바이어스의 부호가 반대로 된 경우에 대응할 수 없어서 모델에 극성의 문제가 생긴다. 또한, 직류 바이어스의 값이 돌발적으로 변화하면 발산한 값으로 변환되는 문제도 생긴다.
본 발명은 이러한 과제를 해결하기 위해 이루어진 것으로,
인덕터의 등가 회로를 수동 회로 소자를 사용하여 나타내고,
직류 전류 중첩 시에서의 수동 회로 소자의 특성 변화율을, 실측값을 기초로 전류를 변수로 하는 근사 함수로서 나타내고,
인덕터에 흐르는 전류를 참조하여, 참조한 전류에 대응해서 근사 함수에 의해 산출되는 특성 변화율, 및 직류 전류 무중첩 시에 수동 회로 소자에 생기는 무중첩 시 전압에 기초하고, 직류 전류의 중첩에 의해 특성이 변화하는 수동 회로 소자에 직렬로 접속되는 제어 전압원에 의해, 직류 전류 중첩 시에 수동 회로 소자에 생기는 중첩 시 전압과 무중첩 시 전압과의 차분 전압을 발생시켜, 무중첩 시 전압에 차분 전압을 중첩시킴으로써,
인덕터의 직류 전류 중첩 시의 비선형 특성을 시뮬레이트하는 인덕터의 시뮬레이션 방법을 구성했다.
또한, 인덕터의 등가 회로를 나타내는 수동 회로 소자와,
인덕터에 흐르는 전류를 참조하는 전류 참조 수단과,
직류 전류 중첩 시의 수동 회로 소자의 특성 변화율을, 실측값을 기초로, 전류를 변수로서 나타내는 근사 함수에 의해, 전류 참조 수단에 의해 참조되는 전류에 대응하여 산출되는 특성 변화율, 및 직류 전류 무중첩 시에 수동 회로 소자에 생기는 무중첩 시 전압에 기초하여, 직류 전류 중첩 시에 수동 회로 소자에 생기는 중첩 시 전압과 무중첩 시 전압과의 차분 전압을 발생시키는, 직류 전류의 중첩에 의해 특성이 변화되는 수동 회로 소자에 직렬로 접속되는 제어 전압원을 포함하여, 인덕터의 비선형 등가 회로 모델을 구성했다.
본 구성의 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델에 있어서는, 직류 전류 중첩 시에서의 수동 회로 소자의 특성 변화율이, 실측값을 기초로, 참조하는 전류를 변수로 한 근사 함수에 의해 표시된다. 따라서, 수동 회로 소자의 특성 변화율은 참조하는 전류에 따라, 이 근사 함수에 의해 산출된다. 또한, 직류 전류 중첩 시에 수동 회로 소자에 생기는 중첩 시 전압은 직류 전류 무중첩 시에 수동 회로 소자에 생기는 무중첩 시 전압에, 중첩 시 전압과 무중첩 시 전압과의 차분 전압을 중첩시킴으로써 구할 수 있다. 이렇기 때문에, 상기의 특성 변화율 및 무중첩 시 전압에 기초하여, 제어 전압원에 의해 중첩 시 전압과 무중첩 시 전압과의 차분 전압을 발생시켜서, 제어 전압원에 수동 회로 소자를 직렬로 접속하여 이 차분 전압을 무중첩 시 전압에 중첩시킴으로써 수동 회로 소자의 중첩 시 전압을 시뮬레이트 할 수 있다.
즉, 인덕터에 흐르는 전류를 참조하여 근사 함수에 의해 수동 회로 소자의 특성 변화율을 산출하고, 특성 변화율 및 무중첩 시 전압에 기초하여 제어 전압원에 의해 차분 전압을 발생시킴으로써 임의의 직류 중첩 전류에 대하여 동적으로 추종 가능한 시뮬레이션을 할 수 있다. 이 결과, 인덕터의 직류 전류 중첩 시의 비선형 특성을 정밀도 높게 동적으로 시뮬레이트 할 수 있는 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델을 단순한 구성으로 간단히 제공할 수 있다. 또한, 인덕터의 비선형 등가 회로 모델은, 상기와 같이 무중첩 시 전압을 기준으로 하여 무중첩 시 전압에 차분 전압을 제어 전압원에 의해 단지 중첩시킴으로써 얻어지므로, 반대로 본 등가 회로 모델로부터 제어 전압원을 제외함으로써, 무중첩 시 전압에 대응한, 즉 직류 전류가 중첩하지 않을 경우의 인덕터의 등가 회로 모델을 간단히 얻을 수 있다.
또한, 본 발명은 차분 전압이, 차분 전압을 ΔV, 무중첩 시 전압을 V0, 근사 함수를 참조하는 전류 x를 변수로 하는 함수 exp(f(x))로 한 경우에, 다음 식
ΔV=V0×[exp(f(x))-1]
의 함수형으로 주어지는 것을 특징으로 한다.
본 구성에 의하면, 참조하는 전류 x의 값이 제로일 때에는 지수 함수인 함수 exp(f(x))의 값이 1이 되고, 무중첩 시 전압(V0)에 승산되는 계수 [exp(f(x))-1]의 값이 제로가 되므로, 차분 전압(ΔV)의 값도 제로가 된다. 또한, 참조하는 전류 x의 값이 제로가 아닐 때에는 함수 exp(f(x))의 값이 1보다 커지고, 무중첩 시 전압(V0)에 승산되는 계수 [exp(f(x))-1]의 값이 제로보다 커진다. 이렇기 때문에, 실용 시에 입각하여 차분 전압(ΔV)의 값이 산출되고, 참조하는 전류(x)의 값이 제로가 아닐 때에는 차분 전압(ΔV)이 전류 x에 대응하여 반드시 산출되어서 인덕터의 특성이 시뮬레이트되어, 인덕터의 정성(定性) 이해에 제공된다.
또한, 본 발명은 근사 함수가 홀수차의 거듭제곱항을 포함하지 않는 다항식 형식의 우함수로 주어지는 것을 특징으로 한다.
본 구성에 의하면, 근사 함수가 홀수차의 거듭제곱항을 포함하지 않는 다항식 형식으로 표현되기 때문에, 종래의 인덕터의 시뮬레이션과는 달리 직류 바이어스의 부호가 반대로 된 경우나 직류 바이어스의 값이 돌발적으로 변화한 경우 등에도 수동 회로 소자의 특성 변화율은 근사 함수에 의해 적절하게 근사된다.
또한, 본 발명은 인덕터에 흐르는 전류가 등가 회로의 입력단 또는 출력단으로 참조되어, 무중첩 시 전압이 수동 회로 소자의 양단에서 참조되는 것을 특징으로 한다.
본 구성에 의하면, 참조하는 전류 또는 무중첩 시 전압을 등가 회로 모델과 별도로 설정하여 계산하는 방식과는 달리, 등가 회로 모델 내의 입력단 혹은 출력 단에 생기는 순시 전류 또는 등가 회로 모델 내의 수동 회로 소자 양단에 생기는 순시 전압이 참조되어서 차분 전압의 연산이 이루어진다. 이렇기 때문에, 차분 전압의 연산에 이용되는 전류 및 무중첩 시 전압은 시간이 늦어지지 않게 참조되어, 인덕터의 비선형 특성에 대한 과도 응답 해석을 고속이면서 고정밀도로 실시할 수 있게 된다.
또한, 본 발명은 제어 전압원에 직렬로 접속되는 수동 회로 소자가, 유도 소자 단체 또는 유도 소자와 저항 소자의 직렬 회로, 또는 유도 소자와 저항 소자와 용량 소자의 직렬 회로인 것을 특징으로 한다.
본 구성에 의하면, 직류 전류가 중첩되지 않을 경우의 수동 회로 소자의 특성은 유도 소자 단체 또는 유도 소자와 저항 소자의 직렬 회로 또는 유도 소자와 저항 소자와 용량 소자의 직렬 회로에 의해 표시된다. 그리고 이들 회로에 직렬로 제어 전압원이 접속됨으로써 수동 회로 소자의 직류 전류가 중첩될 경우의 특성이 시뮬레이트된다.
또한, 본 발명은 제어 전압원과 수동 회로 소자와의 직렬 회로가 복수 병렬로 접속되는 것을 특징으로 한다.
본 구성에 의하면, 제어 전압원과 수동 회로 소자와의 간단한 직렬 회로가 단지 병렬로 접속되고, 그 병렬수가 늘어남으로써 등가 회로 모델의 시뮬레이트 정밀도가 높아진다. 이렇기 때문에, 시뮬레이트 정밀도가 높은 등가 회로 모델을 규칙적이면서 예측할 수 있게 구성할 수 있다. 또한, 제어 전압원과 수동 회로 소자의 직렬 회로를 복수 병렬로 접속하는 것일 뿐이므로, 계통적인 연산 수속에 의해 수동 회로 소자의 직류 전류 중첩 시의 특성을 시뮬레이트 할 수 있다.
또한, 본 발명은 등가 회로가 인덕터에 대한 직류 전류의 중첩에 의해 특성이 변화되지 않는 수동 회로 소자를 포함하여 구성되는 것을 특징으로 한다.
본 구성에 의하면, 직류 전류의 중첩에 의해 특성이 변화하는 수동 회로 소자와, 직류 전류의 중첩에 의해 특성이 변화하지 않는 수동 회로 소자가 조합되어서 등가 회로 모델이 구성되므로, 인덕터의 비선형 특성의 시뮬레이션을 더욱 고정밀도화함과 함께, 시뮬레이션의 주파수 대역을 광대역화 할 수 있게 된다.
또한, 본 발명은,
인덕터의 종류를 입력하는 제1 스텝과,
인덕터에 인가하는 전압 또는 인덕터에 흘려보내는 전류를 입력하는 제2 스텝과,
제2 스텝에서 입력된 전압 또는 전류에 의해 인덕터에 흐르는 전류를 참조하고, 제1 스텝에서 입력된 종류의 인덕터에 대해서 미리 준비된 근사 함수에 의해 참조 전류에 대응하여 산출되는 특성 변화율 및 무중첩 시 전압에 기초하여, 제어 전압원에 의해 차분 전압을 발생시켜, 무중첩 시 전압에 차분 전압을 중첩시킴으로써 인덕터의 직류 전류 중첩 시의 비선형 특성을 시뮬레이트하는 제3 스텝을 포함하고,
상기 중 어느 하나의 인덕터의 시뮬레이션 방법을 실시하는, 또는 상기 중 어느 하나의 인덕터의 비선형 등가 회로 모델을 기능시키는 컴퓨터 프로그램을 구성했다.
본 구성에 의하면, 시뮬레이트하는 인덕터의 종류, 및 인덕터에 인가하는 전압 또는 인덕터로 흘려보내는 전류의 값이 컴퓨터 프로그램에 입력됨으로써, 입력된 종류의 인덕터의 비선형 특성은, 컴퓨터 프로그램에 의해 수동 회로 소자의 무중첩 시 전압에 차분 전압이 중첩되어 자동적으로 시뮬레이션된다. 이렇기 때문에, 본 시뮬레이션 방법 또는 본 비선형 등가 회로 모델의 이용자는 시뮬레이트하는 인덕터의 종류, 및 인덕터에 인가하는 전압 또는 인덕터로 흘려보내는 전류의 값을 컴퓨터 프로그램에 입력하는 것만으로 정확한 회로 시뮬레이션을 고정밀도이면서 간단히 실시할 수 있다. 이 결과, 회로 시뮬레이트에 관한 전문지식을 가지지 않는 일반의 이용자여도, 인덕터를 이용한 전자회로의 정확한 회로 시뮬레이션을 고정밀도이면서 간단히 실시할 수 있다.
또한, 본 발명은 상기 컴퓨터 프로그램을 포함하는 서버에 인터넷망을 통해서 액세스하고, 인터넷망에 접속된 단말로부터 상기 컴퓨터 프로그램을 사용하는 컴퓨터 프로그램의 사용 방법을 구성했다.
본 구성에 의하면, 이용자는 인터넷망에 접속된 단말로부터 상기 컴퓨터 프로그램을 포함하는 서버에 액세스함으로써 상기 컴퓨터 프로그램을 용이하게 사용할 수 있다. 이렇기 때문에, 본 발명에 의한 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델을 다수의 이용자에게 제공할 수 있게 된다.
본 발명에 의하면, 상기와 같이 인덕터의 직류 전류 중첩 시의 비선형 특성을 정밀도 높게 동적으로 시뮬레이트 할 수 있는 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델을 단순한 구성으로 간단히 제공할 수 있다.
도 1(a)는 본 발명의 제1 실시형태에서의 인덕터의 직류 전류 무중첩 시의 수동 등가 회로 모델을 나타내는 회로도이고, 도 1(b)는 제1 실시형태에서의 직류 전류 중첩 시의 인덕터의 비선형 등가 회로 모델을 나타내는 회로도이다.
도 2(a)는 가변 저항 소자(RX1(Idc)) 및 가변 유도 소자(LX1(Idc))를 사용하여 표시되는 직류 전류 중첩 시의 인덕터의 비선형 등가 회로 모델을 나타내는 회로도이고, 도 2(b)는 제1 실시형태에 있어서, 가변 저항 소자(RX1(Idc)) 및 가변 유도 소자(LX1(Idc))가 등가 모델화되어서 표시되는 직류 전류 중첩 시의 인덕터의 비선형 등가 회로 모델을 나타내는 회로도이다.
도 3(a)는 제1 실시형태에 있어서, 제어 전압원(BR)과 제어 전압원(BL)을 이용하여 직류 전류 중첩 시의 비선형 특성이 표시된, 저항 소자(R1)와 유도 소자(L1)의 직렬 회로이고, 도 3(b)는 1개의 제어 전압원(BR+BL)을 이용하여 직류 전류 중첩 시의 비선형 특성이 표시된, 저항 소자(R1)와 유도 소자(L1)의 직렬 회로이다.
도 4(a)는 도 1(b)에 나타내는 비선형 등가 회로 모델로부터 산출되는 인덕터의 임피던스(Z)의 크기(MagZ)에 대한 주파수 특성과, 도 1(a)에 나타내는 수동 등가 회로 모델로부터 산출되는 특성과 비교하여 나타내는 그래프이고, 도 4(b)는 도 1(b)에 나타내는 비선형 등가 회로 모델로부터 산출되는 인덕터의 등가 직렬 저항(ESR)의 주파수 특성과, 도 1(a)에 나타내는 수동 등가 회로 모델로부터 산출되는 특성을 비교하여 나타내는 그래프이다.
도 5(a)는 본 발명의 제2 실시형태에서의 인덕터의 직류 전류 무중첩 시의 수동 등가 회로 모델을 나타내는 회로도이고, 도 5(b)는 제2 실시형태에서의 직류 전류 중첩 시의 인덕터의 비선형 등가 회로 모델을 나타내는 회로도이다.
도 6(a)는 도 5(b)에 나타내는 비선형 등가 회로 모델로부터 산출되는 인덕터의 임피던스(Z)의 크기(MagZ)에 대한 주파수 특성과, 도 5(a)에 나타내는 수동 등가 회로 모델로부터 산출되는 특성과 비교하여 나타내는 그래프이고, 도 6(b)는 도 5(b)에 나타내는 비선형 등가 회로 모델로부터 산출되는 인덕터의 등가 직렬 저항(ESR)의 주파수 특성과, 도 5(a)에 나타내는 수동 등가 회로 모델로부터 산출되는 특성을 비교하여 나타내는 그래프이다.
도 7의 (a), (b), (c)는 본 발명의 각 실시형태에서의 수동 등가 회로 모델에서 사용되는, 직류 전류(Idc)가 중첩되지 않을 경우의 특성을 나타내는 수동 회로 소자의 구성을 나타내는 도면이고, (d), (e), (f)는 각 실시형태에서의 비선형 등가 회로 모델로 사용되는, 직류 전류(Idc)가 중첩될 경우의 비선형 특성을 나타내는 수동 회로 소자의 구성을 나타내는 도면이다.
도 8(a)는 인덕터에 대한 직류 전류(Idc)의 중첩에 의해 특성이 변화하지 않는 수동 회로 소자(r, c)를 포함하여 구성되는 어드미턴스 전개형 수동 등가 회로 모델을 나타내는 회로도이고, 도 8(b)는 이들 수동 회로 소자(r, c)를 포함하여 구성되는 어드미턴스 전개형 비선형 등가 회로 모델을 나타내는 회로도이다.
도 9는 도 8 (b)에 나타내는 비선형 등가 회로 모델을 일반형으로 하여 나타낸, 본 발명의 제3 실시형태에서의 비선형 등가 회로 모델의 회로도이다.
도 10은 도 9에 나타내는 비선형 등가 회로 모델의 제1 변형예를 나타내는 회로도이다.
도 11은 도 9에 나타내는 비선형 등가 회로 모델의 제2 변형예를 나타내는 회로도이다.
도 12는 직류 전류(Idc)의 중첩에 의해 특성이 변화하지 않는 수동 회로 소자가 조합되어 구성된, 본 발명의 제4 실시형태에서의 인덕터의 어드미턴스 전개형 등가 회로 모델의 구체적인 일례를 나타내는 회로도이다.
도 13(a)는 등가 회로 모델이 도 12에 나타내는 인덕터의 등가 직렬 저항(ESR)의 주파수 특성을 각 직류 중첩 전류(Idc)에 대해서 측정한 결과를 나타내는 그래프이고, 도 13(b)는 등가 회로 모델이 도 12에 나타내는 인덕터의 등가직렬 인덕턴스(ESL)의 주파수 특성을 각 직류 중첩 전류(Idc)에 대해서 측정한 결과를 나타내는 그래프이다.
도 14(a)는 도 12에 나타내는 등가 회로 모델을 이용하여 계산되는 직류 전류 무중첩 시의 인덕터의 등가 직렬 저항(ESR)의 주파수 특성과, 도 13(a)에 나타내는 측정된 주파수 특성과 비교하여 나타내는 그래프이고, 도 14(b)는 도 12에 나타내는 등가 회로 모델을 이용하여 계산되는 직류 전류 무중첩 시의 인덕터의 등가직렬 인덕턴스(ESL)의 주파수 특성과, 도 13(b)에 나타내는 측정된 주파수 특성을 비교하여 나타내는 그래프이다.
도 15(a), (b)는 도 12에 나타내는 수동 회로 소자(LX1, RX1)의 특성 변화율(kL1(Idc)), (kR1(Idc))을 나타내는 그래프이다.
도 16(a), (b)는 도 12에 나타내는 수동 회로 소자(LX2, RX2)의 특성 변화율(kL2(Idc), (kR2(Idc))을 나타내는 그래프이다.
도 17(a)는 도 12에 나타내는 등가 회로 모델에 기초한, 제4 실시형태에서의 인덕터의 직류 전류 무중첩 시의 수동 등가 회로 모델을 나타내는 회로도이고, 도 17(b)는 동일하게 도 12에 나타내는 등가 회로 모델에 기초한, 제4 실시형태에서의 직류 전류 중첩 시의 인덕터의 비선형 등가 회로 모델을 나타내는 회로도이다.
도 18(a)는 도 17에 나타내는 등가 회로 모델을 이용하여 계산되는, 각 직류 전류(Idc)의 중첩 시 및 무중첩 시의 인덕터의 등가 직렬 저항(ESR)의 주파수 특성과, 도 13(a)에 나타내는 측정된 주파수 특성과 비교하여 나타내는 그래프이고, 도 18(b)는 도 17에 나타내는 등가 회로 모델을 이용하여 계산되는, 각 직류 전류(Idc)의 중첩 시 및 무중첩 시의 인덕터의 등가직렬 인덕턴스(ESL)의 주파수 특성과, 도 13(b)에 나타내는 측정된 주파수 특성과 비교하여 나타내는 그래프이다.
도 19는 등가 회로 모델이 도 17에 표시되는 인덕터에 대한, 인가 전압에 대한 전류의 과도 응답특성을 나타내는 그래프이다.
다음으로, 본 발명에 의한 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델을, 리니어 테크놀로지사가 제공하는 LTspice에 적용한 실시형태에 대해서 설명한다.
도 1(a)는 제1 실시형태에서의 인덕터의 직류 전류 무중첩 시의 수동 등가 회로 모델을 나타내는 회로도이고, 도 1(b)는 제1 실시형태에서의 직류 전류 중첩 시의 인덕터의 비선형 등가 회로 모델을 나타내는 회로도이다.
각 등가 회로 모델의 유도 소자(L1)와 저항 소자(R1)와의 직렬 회로는, 시뮬레이션의 대상으로 하는 인덕터의 등가 회로를 나타내는 수동 회로 소자를 구성한다. 도 1(a)에 나타내는 수동 등가 회로 모델에서는 직류 전류(Idc)가 중첩하지 않는 교류 전류(Iac)가, LTspice의 전류원 모델(I0)에 의해 수동 회로 소자(L1, R1)에 전류(I)로서 흐른다. 또한, 도 1(b)에 나타내는 비선형 등가 회로 모델에서는 직류 전류(Idc)가 중첩한 교류 전류(Iac)가, LTspice의 전류원 모델(I1)에 의해 수동 회로 소자(L1, R1)에 전류(I)로서 흐른다. 여기서, 각 등가 회로에서의 유도 소자(L1)의 회로 정수는 5[μH], 저항 소자(R1)의 회로 정수는 0.2[Ω], 직류 중첩 전류(Idc)는 1[A]로 설정했다.
도 1(a)에 나타내는 수동 등가 회로 모델에서의 전압원 모델(V0) 및 도 1(b)에 나타내는 비선형 등가 회로 모델의 전압원 모델(V1)은 각 등가 회로 모델에 흐르는 전류(I)를 계측한다. 이들 전압원 모델(V0, V1)은 각각의 회로에서의 전류를 계측하기 위해서 설정한, LTspice에서의 편의상 컴포넌트이며, 설정 전압(V)이 0[V]로 설정되어서 전류계로서 대용된다. 본 실시형태에서는, 계측된 전류(I) 중 직류 전류(Idc)가 참조 전류(Iref)이다. 전압원 모델(V0, V1)은 시뮬레이션의 대상으로 하는 인덕터에 흐르는 전류(Iref)를 참조하는 전류 참조 수단을 구성한다. 또한, 본 실시형태에서는, 전류(I) 중 직류 전류(Idc)만을 참조하고 있지만, 직류 전류(Idc)와 교류 전류(Iac)의 쌍방을 참조하여 참조 전류(Iref)로 하도록 해도 된다.
도 1(b)에 나타내는 제어 전압원(B1)은 수동 회로 소자(L1, R1)의 회로 정수의 변화율, 즉 특성 변화율(kL1(Idc), kR1(Idc)) 및 무중첩 시 전압(VL1, VR1)에 기초하여, 후술하는 바와 같이 중첩 시 전압(VL1(Idc), VR1(Idc))과 무중첩 시 전압(VL1, VR1)의 차분 전압(ΔVL1, ΔVR1)을 발생시킨다. 여기서, 특성 변화율(kL1(Idc), kR1(Idc))은 직류 전류(Idc)의 무중첩 시에서의 수동 회로 소자(L1, R1)의 회로 정수에 대한, 직류 전류(Idc)의 중첩 시에서의 회로 정수의 비이다. 또한, 무중첩 시 전압(VL1, VR1)은 직류 전류(Idc)의 무중첩 시에 수동 회로 소자(L1, R1)에 생기는 전압이며, 중첩 시 전압(VL1(Idc), VR1(Idc))은 직류 전류(Idc)의 중첩 시에 수동 회로 소자(L1, R1)에 생기는 전압이다. 이 제어 전압원(B1)은 도시하는 바와 같이, 직류 전류(Idc)의 중첩에 의해 특성이 변화하는 수동 회로 소자(L1, R1)에 직렬로 접속된다.
각 등가 회로 모델에서의 수동 회로 소자(L1, R1)는 그 회로 정수가 직류 중첩 전류(Idc)에 의존하지 않는 불변인 것을 나타내고 있고, 그 특성 변화에 의한 전압 변동분은 제어 전압원(B1)에 의해 표시된다. 제어 전압원(B1)은 LTspice에 있어서 비헤이비어 전압원 모델(behavior voltage source model)로서 취급되는 LTspice의 컴포넌트이며, 참조하는 전류(Iref) 및 무중첩 시 전압(VL1, VR1)에 의존하여 그 자신의 값이 종속적으로 결정된다.
직류 전류(Idc)의 중첩에 의해 저항 소자(R1)와 유도 소자(L1)의 회로 정수가 변화될 경우, 본래라면 도 2(a)의 회로도에 나타내는 바와 같이, 저항 소자(R1)는 직류 중첩 전류(Idc)의 값에 의존하여 저항값이 변화하는 가변 저항 소자(RX1(Idc)), 유도 소자(L1)는 직류 중첩 전류(Idc)의 값에 의존하여 리액턴스(reactance)가 가변하는 가변 유도 소자(LX1(Idc))로 표시된다. 여기서, 저항 소자(R1)와 유도 소자(L1)의 직렬 회로에는 저항 소자(R2)와 용량 소자(C2)의 직렬 회로가 병렬로 접속되어 있는 것으로 한다. 이들 저항 소자(R2) 및 용량 소자(C2)는 직류 전류(Idc)의 중첩에 의해 회로 정수가 변화되지 않고 불변인 것으로 한다. 이렇기 때문에, 회로도에서의 저항 소자(R2) 및 용량 소자(C2)의 소자 기호에는 가변 유도 소자(LX1(Idc)) 및 가변 저항 소자(RX1(Idc))와 달리 가변을 나타내는 화살표가 없다.
이 회로의 입력 단자에 도시하는 바와 같이, 전류(I)가 흐르고, 가변 저항 소자(RX1(Idc))와 가변 유도 소자(LX1(Idc))와의 직렬 회로에 전류(I1), 저항 소자(R2)와 용량 소자(C2)의 직렬 회로에 전류(I2)가 흐르는 것이라고 가정한다. 이 경우, 회로에 흐르는 전류(I) 및 회로의 입출력 단자간에 생기는 전압(V)은 각각 다음의 (4)식 및 (5)식으로 표시된다.
Figure 112015110688576-pct00004
본 실시형태에서는 직류 전류(Idc)의 중첩에 의해 회로 정수가 변화하는 가변 저항 소자(RX1(Idc))를, 도 3(a)에 나타내는 바와 같이 LTspice에서 비헤이비어 전압원 모델로서 취급되는 제어 전압원(BR)과 저항 소자(R1)의 직렬 회로(1)로서 나타낸다. 또한, 직류 전류(Idc)의 중첩에 의해 회로 정수가 변화하는 가변 유도 소자(LX1(Idc))를, 도 3(a)에 나타내는 바와 같이, LTspice에서 비헤이비어 전압원 모델로서 취급되는 제어 전압원(BL)과 유도 소자(L1)의 직렬 회로(2)로서 나타낸다. 이 저항 소자(R1) 및 유도 소자(L1)는 도 1에 나타낸 것과 동일하며, 그 회로 정수가 직류 중첩 전류(Idc)에 의존하지 않는 불변인 것을 나타내고 있다.
제어 전압원(BR)은 직류 전류(Idc)의 중첩에 의해 저항 소자(R1)에 생기는 전압의 변동분을, 직류 전류(Idc)의 값에 따라서 차분 전압(ΔVR1(Idc))으로서 발생시킨다. 제어 전압원(BL)은 직류 전류(Idc)의 중첩에 의해 유도 소자(L1)에 생기는 전압의 변동분을, 직류 전류(Idc)의 값에 따라 차분 전압(ΔVL1(Idc))으로서 발생시킨다. 직렬로 접속되는 이들 제어 전압원(BR 및 BL)은 도 3(b)에 나타내는 바와 같이, 1개의 제어 전압원(BR+BL)으로서 나타낼 수 있다. 단, 제어 전압원(BR+BL)의 발생하는 전압의 값은 각 제어 전압원(BR 및 BL)이 발생하는 전압의 값의 합이 된다.
가변 저항 소자(RX1(Idc)) 및 가변 유도 소자(LX1(Idc))를 상기와 같이 등가 모델화함으로써, 본래 도 2(a)에 표시되는 회로는, 도 2(b)에 나타내는 본 실시형태의 비선형 등가 회로 모델로 치환된다. 즉, 도 2(a)에 표시되는 가변 저항 소자(RX1(Idc))와 가변 유도 소자(LX1(Idc))의 직렬 회로는, 저항 소자(R1)와 유도 소자(L1)의 직렬 회로에 제어 전압원(B1)을 직렬로 접속한 도 2(b)에 나타내는 회로로 치환된다. 이 제어 전압원(B1)은 도 3(b)에서의 제어 전압원(BR+BL)에 상당하고, 직류 전류(Idc)의 중첩에 의해 저항 소자(R1)와 유도 소자(L1)의 양쪽에 생기는 전압의 변동분을, 직류 중첩 전류(Idc)의 값에 따라 차분 전압(ΔVX1(Idc))으로서 발생시킨다. 또한, 저항 소자(R1)의 양단에는 무중첩 시 전압(VR1), 유도 소자(L1)의 양단에는 무중첩 시 전압(VL1)이 생긴다.
이 경우, 비선형 등가 회로 모델의 입출력 단자간의 전압(V)은 다음 (6)식으로 표시된다.
Figure 112015110688576-pct00005
또한, 차분 전압(ΔVX1(Idc))은 저항 소자(R1)에 생기는 전압의 변동분인 차분 전압(ΔVR1(Idc))과 유도 소자(L1)에 생기는 전압의 변동분인 차분 전압(ΔVLI(Idc))의 합이기 때문에, 다음 (7)식으로 표시된다.
Figure 112015110688576-pct00006
또한, 차분 전압(ΔVR1(Idc))은 저항 소자(R1)의 중첩 시 전압(VR1(Idc))과 무중첩 시 전압(VR1)의 차이이기 때문에, 다음 (8)식으로 표시된다.
Figure 112015110688576-pct00007
여기서, 특성 변화율(kR1(Idc))은 직류 전류(Idc)의 무중첩 시에서의 수동 회로 소자(R1)의 회로 정수(R1)에 대한, 직류 전류(Idc)의 중첩 시에서의 회로 정수(RX1(Idc))의 비이며, 다음 (9)식으로 표시된다.
Figure 112015110688576-pct00008
제어 전압원(B1)은 (8)식에 나타내는 바와 같이 특성 변화율(kR1(Idc))로부터 1을 감산한 값에 무중첩 시 전압(VR1)을 승산함으로써, 즉 특성 변화율(kR1(Idc)) 및 무중첩 시 전압(VR1)에 기초하여 중첩 시 전압(VR1(Idc))과 무중첩 시 전압(VR1)의 차분 전압(ΔVR1(Idc))을 발생시킨다.
또한, 차분 전압(ΔVL1(Idc))은 유도 소자(L1)의 중첩 시 전압(VL1(Idc))과 무중첩 시 전압(VL1)의 차이이기 때문에, 다음의 (10)식에 의해 표시된다.
Figure 112015110688576-pct00009
여기서, 특성 변화율(kL1(Idc))은 직류 전류(Idc)의 무중첩 시에서의 수동 회로 소자(L1)의 회로 정수(L1)에 대한, 직류 전류(Idc)의 중첩 시에서의 회로 정수(LX1(Idc))의 비이며, 다음 (11)식으로 표시된다.
Figure 112015110688576-pct00010
제어 전압원(B1)은 (10)식에 나타내는 바와 같이 특성 변화율(kL1(Idc))로부터 1을 감산한 값에 무중첩 시 전압(VL1)을 승산함으로써, 즉 특성 변화율(kL1(Idc)) 및 무중첩 시 전압(VL1)에 기초하여 중첩 시 전압(VL1(Idc))과 무중첩 시 전압(VL1)의 차분 전압(ΔVL1(Idc))을 발생시킨다.
또한 본 실시형태에서는, 도 1에 나타내는 바와 같이, 전류 참조 수단을 구성하는 전압원 모델(V0, V1)이 등가 회로 모델의 입력단에 마련되어 있다. 제어 전압원(B1)은 차분 전압(ΔVL1(Idc), ΔVR1(Idc))을 발생시킬 때, 등가 회로 모델의 입력단에서 전류(Iref)를 참조한다. 한편, 전압원 모델(V0, V1)을 등가 회로 모델의 출력단에 마련하고, 제어 전압원(B1)이 등가 회로 모델의 출력단에서 전류(Iref)를 참조하도록 해도 된다. 또한, 제어 전압원(B1)은, 차분 전압(ΔVL1, ΔVR1)을 발생시킬 때, 수동 회로 소자(L1, R1)의 양단에서 무중첩 시 전압(VL1, VR1)을 참조한다.
특성 변화율(kR1(Idc), kL1(Idc))은 근사 함수 exp(f(x))에 의해 전압원 모델(V0, V1)에 따라 계측되는 참조 전류(Iref(=Idc))에 대응하여 산출된다. 근사 함수 exp(f(x))는 실측값을 기초로, 인덕터에 흐르는 참조 전류(Iref)를 변수 x로 ㅎ하여 후술하는 바와 같이(도 15, 도 16 참조) 표시된다. 본 실시형태에서 근사 함수 exp(f(x))는 홀수차의 거듭제곱항을 포함하지 않는 다항식 형식의 우함수로 주어진다. 또한, (8)식 및 (10)식으로 표시되는 차분 전압(ΔVR1(Idc), ΔVL1(Idc))은 이 근사 함수 exp(f(x))를 이용하여 다음 (12)식 및 (13)식의 함수형으로 주어진다.
Figure 112015110688576-pct00011
본 실시형태의 인덕터의 시뮬레이션 방법에서는, 우선 인덕터의 등가 회로를 수동 회로 소자(L1, R1)의 직렬 회로를 사용하여 표시하고, 도 1(b)에 나타내는 비선형 등가 회로 모델을 구축한다. 그리고, 직류 전류(Idc)의 중첩 시에서의 수동 회로 소자(L1, R1)의 특성 변화율(kL1(Idc), kR1(Idc))을, 실측값을 기초로, 참조 전류(Iref)를 변수 x로 하는 근사 함수 exp(f(x))로서 나타낸다. 다음으로, 수동 회로 소자(L1, R1)에 직렬로 접속되는 제어 전압원(B1)에 의해 전압원 모델(V0, V1)에 따라 계측되는 참조 전류(Iref)가 참조되고, 근사 함수 exp(f(x))가 이용되어, 참조 전류(Iref)에 대응하여 특성 변화율(kL1(Idc), kR1(Idc))이 산출된다. 또한, 제어 전압원(B1)에 의해 수동 회로 소자(L1, R1)의 양단에 생기는 무중첩 시 전압(VL1, VR1)이 참조된다. 그리고, 특성 변화율(kL1(Idc), kR1(Idc)) 및 무중첩 시 전압(VL1, VR1)에 기초하여, 제어 전압원(B1)에 의해 중첩 시 전압(VL1(Idc), VR1(Idc))과 무중첩 시 전압(VL1, VR1)의 차분 전압(ΔVL1(Idc), ΔVR1(Idc))이 발생되어, 무중첩 시 전압(VL1, VR1)에 차분 전압(ΔVL1(Idc), ΔVR1(Idc))이 중첩됨으로써 인덕터의 직류 전류(Idc)의 중첩 시의 비선형 특성이 시뮬레이트된다.
도 4(a)는 상기와 같이 인덕터의 비선형 특성이 시뮬레이트되어서, 도 1(b)에 나타내는 비선형 등가 회로 모델로부터 산출되는 인덕터의 임피던스(Z)의 크기(MagZ)에 대한 주파수 특성과, 도 1(a)에 나타내는 수동 등가 회로 모델로부터 산출되는 동일한 크기(MagZ)에 대한 주파수 특성을 비교하여 나타내는 그래프이다. 이 그래프의 가로축은 주파수[Hz], 세로축은 크기(MagZ)의 값[Ω]을 나타낸다. 또한, 실선으로 나타내는 주파수 특성(A1)은 비선형 등가 회로 모델로부터 산출되는 직류 중첩 전류(Idc)=1[A]일 때의 특성이고, 파선으로 나타내는 주파수 특성(A0)은 수동 등가 회로 모델로부터 산출되는 직류 중첩 전류(Idc)=0일 때의 특성이다.
도 4(b)는 상기와 같이 인덕터의 비선형 특성이 시뮬레이트되어서, 도 1(b)에 나타내는 비선형 등가 회로 모델로부터 산출되는 인덕터의 등가 직렬 저항(ESR)의 주파수 특성과, 도 1(a)에 나타내는 수동 등가 회로 모델로부터 산출되는 동일 등가 직렬 저항(ESR)에 관한 주파수 특성을 비교하여 나타내는 그래프이다. 이 그래프의 가로축은 주파수[Hz], 세로축은 등가 직렬 저항(ESR)의 값[Ω]을 나타낸다. 또한, 실선으로 나타내는 주파수 특성(B1)은 비선형 등가 회로 모델로부터 산출되는 직류 중첩 전류(Idc)=1[A]일 때의 특성이고, 파선으로 나타내는 주파수 특성(B0)은 수동 등가 회로 모델로부터 산출되는 직류 중첩 전류(Idc)=0일 때의 특성이다.
도 4(a)의 그래프에 나타내는 바와 같이 임피던스의 크기(MagZ)에 대한 주파수 특성(A1)은 직류 전류(Idc)의 중첩에 의해 MagZ의 값이 직류 전류 무중첩 시의 주파수 특성(A0)보다 작고, 직류 전류(Idc)의 중첩에 의해 임피던스가 변동되고 있다. 또한, 저주파 영역에서는, 도 4(b)에 나타내는 등가 직렬 저항(ESR)의 특성에 점점 가까워진다. 또한, 도4(b)의 그래프에 나타내는 바와 같이 등가 직렬 저항(ESR)에 대한 주파수 특성(B1)도, 직류 전류(Idc)의 중첩에 의해 ESR의 값이 직류 전류 무중첩 시의 주파수 특성(B0)보다 작고, 직류 전류(Idc)의 중첩에 의해 변동되고 있다. 그러나 주파수 특성 B1 및 B0 모두, 등가 직렬 저항(ESR)의 값은 주파수와 관계없이 일정 값으로 되어 있다.
이러한 제1 실시형태에 의한 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델에 있어서는, 상기와 같이, 직류 전류(Idc)의 중첩 시의 수동 회로 소자(L1, R1)의 특성 변화율(kL1(Idc), kR1(Idc))이, 실측값을 기초로, 인덕터에 흐르는 참조 전류(Iref)를 변수 x로 한 근사 함수 exp(f(x))에 의해 표시된다. 따라서, (9)식, (11)식으로 표시되는 특성 변화율(kL1(Idc), kR1(Idc))은 참조하는 전류(Iref)에 따라 이 근사 함수 exp(f(x))에 의해 산출된다. 또한, 중첩 시 전압(VL1(Idc), VR1(Idc))은 무중첩 시 전압(VL1, VR1)에, (8)식, (10)식으로 표시되는 차분 전압(ΔVL1(Idc), ΔVR1(Idc))을 중첩시킴으로써 구할 수 있다. 이렇기 때문에, 상기의 특성 변화율(kL1(Idc), kR1(Idc)) 및 무중첩 시 전압(VL1, VR1)에 기초하여 제어 전압원(B1)에 의해 차분 전압(ΔVL1(Idc), ΔVR1(Idc))을 발생시켜, 제어 전압원(B1)에 수동 회로 소자(L1, R1)를 직렬로 접속하고, 이 차분 전압(ΔVL1(Idc), ΔVR1(Idc))을 무중첩 시 전압(VL1, VR1)에 중첩시킴으로써 수동 회로 소자(L1, R1)의 중첩 시 전압(VL1(Idc), VR1(Idc))을 시뮬레이트 할 수 있다.
즉, 인덕터에 흐르는 전류(Iref)를 참조하여 근사 함수 exp(f(x))에 의해 수동 회로 소자(L1, R1)의 특성 변화율(kL1(Idc), kR1(Idc))을 산출하고, 특성 변화율(kL1(Idc), kR1(Idc)) 및 무중첩 시 전압(VL1, VR1)에 기초하여 제어 전압원(B1)에 의해 (8)식, (10)식으로 표시되는 차분 전압(ΔVL1(Idc), ΔVR1(Idc))을 발생시킴으로써, 임의의 직류 중첩 전류(Idc)에 대하여 동적으로 추종 가능한 시뮬레이션을 할 수 있다. 이 결과, 인덕터의 직류 전류 중첩 시의 비선형 특성을 정밀도 높게 동적으로 시뮬레이트 할 수 있는 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델을 단순한 구성으로 간단히 제공할 수 있다. 또한, 인덕터의 비선형 등가 회로 모델은, 상기와 같이, 무중첩 시 전압(VL1, VR1)을 기준으로 하여, 무중첩 시 전압(VL1, VR1)에 차분 전압(ΔVL1(Idc), ΔVR1(Idc))을 제어 전압원(B1)에 의해 단지 중첩시킴으로써 얻어지므로, 반대로 도 1(b)에 표시되는 비선형 등가 회로 모델로부터 제어 전압원(B1)을 제외함으로써, 무중첩 시 전압(VL1, VR1)에 대응한, 즉 직류 전류(Idc)가 중첩되지 않을 경우의 도 1(a)에 표시되는 인덕터의 등가 회로 모델을 간단히 얻을 수 있다.
또한, 제1 실시형태에 의한 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델에 의하면, 참조 전류(Iref)의 값이 제로(x=0)일 때에는 지수 함수인 근사 함수 exp(f(x))의 값이 1이 되고, 무중첩 시 전압(VL1, VR1)에 승산되는, (12)식 및 (13)식으로 표시되는 계수 [exp(f(x))-1]의 값이 제로가 되므로, 차분 전압(ΔVL1(Idc), ΔVR1(Idc))의 값도 제로가 된다. 또한, 참조 전류(Iref)의 값이 제로가 아닐 때에는 근사 함수 exp(f(x))의 값이 1보다 커지고, 무중첩 시 전압(V0)에 승산되는 계수 [exp(f(x))-1]의 값이 제로보다 커지기 때문에, 차분 전압(ΔVL1(Idc), ΔVR1(Idc))이 참조 전류(Iref)에 대응하여 반드시 산출되어서 인덕터의 특성이 시뮬레이트되어, 인덕터의 정성 이해에 제공된다.
또한, 제1 실시형태에 의한 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델에 의하면, 근사 함수 exp(f(x))가 홀수차의 거듭제곱항을 포함하지 않는 다항식 형식의 우함수로 표현되기 때문에 종래의 인덕터의 시뮬레이션과는 달라서, 직류 바이어스의 부호가 반대로 된 경우나 직류 바이어스의 값이 돌발적으로 변화된 경우 등에도, 특성 변화율(kL1(Idc), kR1(Idc))은 근사 함수 exp(f(x))에 의해 적절하게 근사된다.
또한, 제1 실시형태에 의한 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델에 의하면, 참조 전류(Iref) 또는 무중첩 시 전압(VL1, VR1)을 등가 회로 모델과 별도로 설정하여 계산하는 방식과는 달리, 등가 회로 모델 내의 입력단 혹은 출력단에 생기는 순시 전류, 또는 등가 회로 모델 내의 수동 회로 소자(L1, R1) 양단에 생기는 순시 전압이 참조되어서, 차분 전압(ΔVL1(Idc), ΔVR1(Idc))의 연산이 이루어진다. 이렇기 때문에, 차분 전압(ΔVL1(Idc), ΔVR1(Idc))의 연산에 이용되는 참조 전류(Iref) 및 무중첩 시 전압(VL1, VR1)은 시간이 늦어지지 않게 참조되어, 인덕터의 비선형 특성에 관한 과도 응답 해석을 고속이면서 고정밀도로 실시할 수 있게 된다.
도 5(a)는 본 발명의 제2 실시형태에서의 인덕터의 직류 전류 무중첩 시의 수동 등가 회로 모델, 도 5(b)는 제2 실시형태에서의 직류 전류 중첩 시의 인덕터의 비선형 등가 회로 모델을 나타내는 회로도이다. 한편, 도 5에 있어서 도 1과 동일 또는 상응하는 부분에는 동일부호를 붙이고 그 설명은 생략한다.
제2 실시형태에서의 수동 등가 회로 모델 및 비선형 등가 회로 모델에서는, 유도 소자(L1)와 저항 소자(R1)의 직렬 회로에 병렬로, 유도 소자(L2)와 저항 소자(R2)의 직렬 회로가 접속되어서, 시뮬레이션 대상으로 하는 인덕터의 등가 회로를 나타내는 수동 회로 소자가 구성되어 있다. 그리고 도 5(b)에 나타내는 비선형 등가 회로 모델에서는, 유도 소자(L1)와 저항 소자(R1)의 직렬 회로에 직렬로 제어 전압원(B1)이 접속되고, 유도 소자(L2)와 저항 소자(R2)의 직렬 회로에 직렬로 제어 전압원(B1)과 같은 제어 전압원(B2)이 접속되어 있다. 즉, 이 제2 실시형태에서의 각 등가 회로 모델에서는, 제어 전압원(B1)과 수동 회로 소자(L1, R1)와의 직렬 회로와, 제어 전압원(B2)과 수동 회로 소자(L2, R2)와의 직렬 회로가 복수 병렬로 접속되어 있다.
이 제2 실시형태에 있어서도, 도 5(a)에 나타내는 수동 등가 회로 모델에서는 직류 전류(Idc)가 중첩하지 않고 있는 교류 전류(Iac)가 전류원 모델(I0)에 의해, 수동 회로 소자(L1, R1)와 수동 회로 소자(L2, R2)의 병렬 회로에 전류(I)로서 흐른다. 또한, 도 5(b)에 나타내는 비선형 등가 회로 모델에서는 직류 전류(Idc)가 중첩한 교류 전류(Iac)가 전류원 모델(I1)에 의해, 수동 회로 소자(L1, R1)와 수동 회로 소자(L2, R2)의 병렬 회로에 전류(I)로서 흐른다. 여기서, 각 등가 회로에서의 유도 소자(L1)의 회로 정수는 4[μH], 저항 소자(R1)의 회로 정수는 0.2[Ω], 유도 소자(L2)의 회로 정수는 100[μH], 저항 소자(R2)의 회로 정수는 1000[Ω], 직류 중첩 전류(Idc)는 1[A]으로 설정했다.
도 5(b)에 나타내는 제어 전압원(B1)은 수동 회로 소자(L1, R1)의 특성 변화율(kL1(Idc), kR1(Idc)) 및 무중첩 시 전압(VL1, VR1)에 기초하여 제1 실시형태의 시뮬레이션 방법과 마찬가지로, 중첩 시 전압(VL1(Idc), VR1(Idc))과 무중첩 시 전압(VL1, VR1)과의 차분 전압(ΔVL1(Idc), ΔVR1(Idc))을 발생시킨다. 또한, 제어 전압원(B2)은 수동 회로 소자(L2, R2)의 특성 변화율(kL2(Idc), kR2(Idc)) 및 무중첩 시 전압(VL2, VR2)에 기초하여 제어 전압원(B1)과 마찬가지로, 중첩 시 전압(VL2(Idc), VR2(Idc))과 무중첩 시 전압(VL2, VR2)과의 차분 전압(ΔVL2(Idc), ΔVR2(Idc))을 발생시킨다.
도 6(a)는 도 5(b)에 나타내는 비선형 등가 회로 모델로부터 산출되는 인덕터의 임피던스(Z)의 크기(MagZ)에 관한 주파수 특성과, 도 5(a)에 나타내는 수동 등가 회로 모델로부터 산출되는 동일 크기(MagZ)에 대한 주파수 특성을, 비교해서 나타내는 그래프이다. 이 그래프의 가로축은 주파수[Hz], 세로축은 크기(MagZ)의 값[Ω]을 나타낸다. 또한, 실선으로 나타내는 주파수 특성(C1)은 비선형 등가 회로 모델로부터 산출되는 직류 중첩 전류(Idc)=1[A]일 때의 특성이고, 파선으로 나타내는 주파수 특성(C0)은 수동 등가 회로 모델로부터 산출되는 직류 중첩 전류(Idc)=0일 때의 특성이다.
도 6(b)는 도 5(b)에 나타내는 비선형 등가 회로 모델로부터 산출되는 인덕터의 등가 직렬 저항(ESR)의 주파수 특성과, 도 5(a)에 나타내는 수동 등가 회로 모델로부터 산출되는 동일 등가 직렬 저항(ESR)에 대한 주파수 특성을 비교해서 나타내는 그래프이다. 이 그래프의 가로축은 주파수[Hz], 세로축은 등가 직렬 저항(ESR)의 값[Ω]을 나타낸다. 또한, 실선으로 나타내는 주파수 특성(D1)은 비선형 등가 회로 모델로부터 산출되는 직류 중첩 전류(Idc)=1[A]일 때의 특성이고, 파선으로 나타내는 주파수 특성(D0)은 수동 등가 회로 모델로부터 산출되는 직류 중첩 전류(Idc)=0일 때의 특성이다.
도 6(a)의 그래프에 나타내는 바와 같이 임피던스의 크기(MagZ)에 대한 주파수 특성(C1)도, 도 4(a)의 그래프에 표시되는 주파수 특성(A1)과 마찬가지로 직류 전류(Idc)의 중첩에 의해 MagZ의 값이 직류 전류 무중첩 시의 주파수 특성(C0)보다 작고, 직류 전류(Idc)의 중첩에 의해 임피던스가 변동하고 있다. 또한, 저주파 영역에서는, 도 6(b)에 나타내는 등가 직렬 저항(ESR)의 특성에 점점 가까워진다. 또한, 도 6(b)의 그래프에 나타내는 바와 같이 등가 직렬 저항(ESR)에 대한 주파수 특성(D1)도, 도 4(b)의 그래프에 나타내는 주파수 특성(B1)과 마찬가지로 직류 전류(Idc)의 중첩에 의해 ESR의 값이 직류 전류 무중첩 시의 주파수 특성(D0)보다 작고, 직류 전류(Idc)의 중첩에 의해 변동하고 있다. 그러나 주파수 특성 D1 및 D0 모두, 도 4(b)의 그래프에 나타내는 주파수 특성(B1 및 B0)과 달리 등가 직렬 저항(ESR)의 값이 주파수에 의해 변화되고, 주파수 특성을 가지는 값으로 되어 있다.
이러한 제2 실시형태에 의한 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델에 의하면, 제어 전압원(B1)과 수동 회로 소자(L1, R1)의 간단한 직렬 회로와, 제어 전압원(B2)과 수동 회로 소자(L2, R2)의 간단한 직렬 회로가 단지 병렬로 접속되고, 그 병렬수가 늘어남으로써 등가 회로 모델의 시뮬레이트 정밀도가 높아진다. 이렇기 때문에, 시뮬레이트 정밀도가 높은 등가 회로 모델을 규칙적이며 예측할 수 있게 구성할 수 있다. 또한, 제어 전압원(B1, B2)과 수동 회로 소자(L1, R1, L2, R2)에 의한 각 직렬 회로를 복수 병렬로 접속하는 것일 뿐이므로, 계통적인 연산 수속에 의해 수동 회로 소자(L1, R1, L2, R2)의 직류 전류 중첩 시의 특성을 시뮬레이트 할 수 있다.
한편, 상기의 제1 실시형태 및 제2 실시형태에 의한 수동 등가 회로 모델에서는, 도 7(b)에 나타내는 바와 같이 유도 소자(L)와 저항 소자(R)의 직렬 회로에 의해 수동 회로 소자를 구성함으로써 직류 전류(Idc)가 중첩되지 않을 경우의 수동 회로 소자의 특성을 나타냈다. 그리고 비선형 등가 회로 모델에서는, 도 7(e)에 나타내는 바와 같이 이 직렬 회로에 직렬로 제어 전압원(B)을 접속함으로써 수동 회로 소자의 직류 전류(Idc)가 중첩되는 경우의 특성을 시뮬레이트 했다.
그러나 도 7(a)에 나타내는 수동 등가 회로 모델과 같이, 유도 소자(L)의 단체에 직렬로 수동 회로 소자를 구성하도록 해도 된다. 이 경우, 비선형 등가 회로 모델은 도 7(d)에 나타내는 바와 같이 유도 소자(L)의 단체에 직렬로 제어 전압원(B)이 접속된다. 또한, 도 7(c)에 나타내는 수동 등가 회로 모델과 같이, 유도 소자(L)와 저항 소자(R)와 용량 소자(C)의 직렬 회로에서 수동 회로 소자를 구성하도록 해도 된다. 이 경우, 비선형 등가 회로 모델은 도 7(f)에 나타내는 바와 같이 이 직렬 회로에 직렬로 제어 전압원(B)이 접속된다. 이렇게 수동 회로 소자를 구성해도, 직류 전류(Idc)가 중첩되지 않을 경우의 수동 회로 소자의 특성은 유도 소자(L)의 단체 또는 유도 소자(L)와 저항 소자(R)와 용량 소자(C)의 직렬 회로에 의해 표시된다. 그리고, 이 회로에 직렬로 제어 전압원(B)이 접속됨으로써 수동 회로 소자의 직류 전류(Idc)가 중첩되는 경우의 특성이 시뮬레이트된다.
또한, 상기의 제1 실시형태 또는 제2 실시형태에서는 모든 수동 회로 소자(L1, R1, L2, R2)가 직류 전류(Idc)의 중첩에 의해 그 특성이 변화될 경우에 대해서 설명했다. 그러나, 도 8(a)에 나타내는 어드미턴스 전개형 수동 등가 회로 모델과 같이, 인덕터에 대한 직류 전류(Idc)의 중첩에 의해 특성이 변화되지 않는 수동 회로 소자(r, c)를 포함하여 등가 회로가 구성되도록 해도 된다. 이 경우, 도 8(b)에 나타내는 어드미턴스 전개형 비선형 등가 회로 모델에서는, 직류 전류(Idc)의 중첩에 의해 특성이 변화되는 수동 회로 소자(Rx, Lx, Cx(x=1, 2, 3,…))에 직렬로 제어 전압원(Bx(x=1, 2, 3,…))이 접속되어, 수동 회로 소자(r, c)에 제어 전압원(Bx)이 접속되지 않음으로써 수동 회로 소자(Rx, Lx, Cx)의 직류 전류(Idc)가 중첩되는 경우의 특성이 시뮬레이트된다. 한편, 도 8에 있어서 도 7과 동일 또는 상응하는 부분에는 동일부호를 붙이고 그 설명은 생략한다. 또한, 여기서는, 등가 회로 모델을 어드미턴스 전개형에 구성한 경우에 대해서 설명하고 있지만, 수동 회로 소자를 직렬로 접속하여 임피던스 전개형에 각 등가 회로 모델을 구성하는 것도 가능하다.
도 9는 도 8(b)에 나타내는 비선형 등가 회로 모델을 일반형으로 하여 나타낸, 본 발명의 제3 실시형태에서의 비선형 등가 회로 모델의 회로도이다. 한편, 도 9에 있어서 도 8(b)과 동일 또는 상응하는 부분에는 동일부호를 붙이고 그 설명은 생략한다.
도 9에 나타내는 일반형의 비선형 등가 회로 모델을 사용한 인덕터의 시뮬레이션 방법에서는, 우선 인덕터의 등가 회로를 수동 회로 소자(R, L, C 및 r, l, c)를 사용하여 나타내어, 도 9에 나타내는 바와 같은 비선형 등가 회로 모델을 구축한다. 유도 소자(l)는 저항 소자(r) 및 용량 소자(c)와 마찬가지로 인덕터에 대한 직류 전류(Idc)의 중첩에 의해 특성이 변화되지 않는 수동 회로 소자이다. 그리고 도 1과 동일한 전류원 모델(I1)에 의해 교류 전류(Iac)에 직류 전류(Idc)가 중첩한 전류I(=Idc+Iac)를 비선형 등가 회로 모델로 흘려보낸다..
전류(I)는 STEP 1에 있어서, 도 1에서의 것과 마찬가지로 전류계로서 대용되는 전압원 모델(V1)에 의해 직류 전류(Idc)가 참조 전류(Iref)로서 참조된다(Iref=Idc). 한편, 여기서는 전류(I) 중 직류 전류(Idc)만을 참조하고 있지만, 직류 전류(Idc)와 교류 전류(Iac)의 양쪽을 참조하도록 해도 된다. 다음으로, STEP 2에 있어서, 직류 전류(Idc)의 중첩 시에서의 수동 회로 소자(Rx, Cx, Lx(x=1,2,3,…))의 특성 변화율(kRX(Iref), kCX(Iref), kLX(Iref)(x=1,2,3,…))을, 측정값을 기초로, 참조 전류(Iref)를 변수 x로 하는 근사 함수 exp(f(x))로서 나타낸다. 그리고, 수동 회로 소자(Rx, Cx, Lx)에 직렬로 접속되는 제어 전압원(BX(x=1, 2, 3,…))에 의해, 전압원 모델(V1)에 의해 계측되는 참조 전류(Iref)가 참조되어, 근사 함수 exp(f(x))가 이용되고, 참조 전류(Iref)에 대응하여 특성 변화율(kRX(Iref), kCX(Iref), kLX(Iref))이 산출된다.
다음으로, STEP 3에 있어서, 제어 전압원(BX)에 의해, 수동 회로 소자(Rx, Cx, Lx)의 양단에 생기는 무중첩 시 전압(VRX, VCX, VLX(x=1, 2, 3,…))이 참조된다. 그리고, 특성 변화율(kRX(Iref), kCX(Iref), kLX(Iref)) 및 무중첩 시 전압(VRX, VCX, VLX)에 기초하여, 제어 전압원(BX)에 의해, 중첩 시 전압(VRX(Iref), VCX(Iref), VLX(Iref))과 무중첩 시 전압(VRX, VCX, VLX)의 차분 전압(ΔVRX(Iref), ΔVCX(Iref), ΔVLX(Iref) (x=1, 2, 3,…))이 보정 전압으로서 발생되어, 무중첩 시 전압(VRX, VCX, VLX)에 차분 전압(ΔVRX(Iref), ΔVCX(Iref), ΔVLX(Iref))이 중첩됨으로써 인덕터의 직류 전류(Idc)의 중첩 시의 비선형 특성이 시뮬레이트된다.
이러한 제3 실시형태에 의한 인덕터의 시뮬레이션에 의하면, 직류 전류(Idc)의 중첩에 의해 특성이 변화하는 수동 회로 소자(Rx, Cx, Lx)와, 직류 전류(Idc)의 중첩에 의해 특성이 변화하지 않는 수동 회로 소자(r, c, l)가 조합되어 등가 회로 모델이 구성되므로, 인덕터의 비선형 특성의 시뮬레이션을 더욱 고정밀도화함과 함께, 시뮬레이션의 주파수 대역을 광대역화 할 수 있게 된다.
한편, 도 9에 나타내는 일반형의 비선형 등가 회로 모델에서는, 전류계로서 대용되는 전압원 모델(V1)을 등가 회로의 입력단에 마련하고, 등가 회로의 입력단측에 있어서 참조 전류(Iref)를 계측한 경우에 대해서 설명했다. 그러나 도 10에 나타내는 바와 같이, 등가 회로의 출력단에 평활 콘덴서(CF)를 통해서 전압원 모델(V1)을 마련하고, 등가 회로의 출력단측에 있어서, 평활 콘덴서(CF)에서 평활한 후의 참조 전류(Iref)를 계측하도록 해도 된다. 한편, 도 10에 있어서, 도 9와 동일 또는 상응하는 부분에는 동일부호를 붙이고 그 설명은 생략한다.
또한, 도 11에 나타내는 바와 같이, 수동 회로 소자 단체의 회로 및 수동 회로 소자의 직렬 회로의 각 입력단에 전류계가 되는 전압원 모델(VX(x=1, 2, 3,…))을 각각 마련하고, 등가 회로의 입력단측에 있어서 참조 전류(Iref)를 계측하도록 해도 된다. 한편, 도 11에 있어서, 도 9와 동일 또는 상응하는 부분에는 동일부호를 붙이고 그 설명은 생략한다.
도 12는 직류 전류(Idc)의 중첩에 의해 특성이 변화하지 않는 수동 회로 소자(r1, c1, r2, c2)가 조합되어 구성된, 본 발명의 제4 실시형태에서의 인덕터의 어드미턴스 전개형 등가 회로 모델의 구체적인 일례를 게시하는 회로도이다. 특성을 시뮬레이트하는 인덕터는, 인덕턴스가 4.7[μH], 정격 전류가 1100[㎃]인 파워 인덕터 제품이다. 이하, 이 파워 인덕터 제품을 공시(供試) 인덕터라고 기재한다. 이 등가 회로 모델에서는, 직류 전류(Idc)의 중첩에 의해 특성이 변화하는 수동 회로 소자가 RX1, RX2, LX1, LX2로서 표시되어 있다. 직류 전류(Idc)의 중첩에 의해 특성이 변화할지 하지 않을지는 실제의 측정값의 분석에 의해 결정된다.
이 어드미턴스 전개형 등가 회로 모델은, 유도 소자(LX1)와 저항 소자(RX1)의 직렬 회로 №.1, 용량 소자(c1)의 단체 회로 №.2, 저항 소자(r1)의 단체 회로 №.3, 유도 소자(LX2)와 저항 소자(RX2)의 직렬 회로 №.4 및 용량 소자(c2)와 저항 소자(r2)의 직렬 회로 №.5가 병렬로 접속되어 구성되어 있다. 직렬 회로 №.1, 회로 №.2 및 회로 №.3의 병렬 회로는 시뮬레이트하는 공시 인덕터의 주공진주파수 부근의 주파수 특성을 실제의 특성에 피팅한다. 직렬 회로 №.4는 시뮬레이트하는 공시 인덕터의 유도성 대역의 주파수 특성을 실제 특성에 피팅한다. 직렬 회로 №.5는 시뮬레이트하는 공시 인덕터의 용량성 대역의 주파수 특성을 실제의 특성에 피팅한다.
상기의 피팅에 의해 도출되는 각 회로 소자의 회로 정수는 직렬 회로 №.1을 구성하는 유도 소자(LX1)가 3.68×10-6[H], 저항 소자(RX1)가 1.93×10-1[Ω]이다. 또한, 회로 №.2의 용량 소자(c1)가 3.65×10-12[F], 회로 №.3의 저항 소자(r1)가 3.89×104[Ω], 직렬 회로 №.4의 유도 소자(LX2)가 1.31×10-4[H], 저항 소자(RX2)가 1.37×103[Ω], 직렬 회로 №.5의 용량 소자(c2)가 2.66×10-12[F], 저항 소자(r2)가 1.75×103[Ω]이다.
도 13(a)는 상기 공시 인덕터의 등가 직렬 저항(ESR)의 주파수 특성을 100[kHz]~100[MHz]의 주파수 대역에서 각 직류 중첩 전류(Idc)에 대해서 측정한 결과를 나타내는 그래프이다. 이 그래프의 가로축은 주파수[Hz], 세로축은 등가 직렬 저항(ESR)의 값[Ω]을 나타낸다. 여기서, 눈금에는 10의 거듭제곱값으로 나타내고 있으며, 예를 들면 1E+5는 1×105(=100000)을 나타내고, "E"는 밑수인 10, +5는 지수를 나타낸다. 또한, 1E-1은 마찬가지로 1×10-1(=0.1)을 나타낸다. 또한, 파선으로 나타내는 주파수 특성(f0)은 직류 중첩 전류(Idc)=0(zero)[㎃]일 때의 주파수 특성이고, 실선으로 나타내는 주파수 특성(f1)은 직류 중첩 전류(Idc)=400[㎃], 일점 쇄선으로 나타내는 주파수 특성(f2)은 직류 중첩 전류(Idc)=800[㎃], 점으로 나타내는 주파수 특성(f3)은 직류 중첩 전류(Idc)=1100[㎃]일 때의 주파수 특성이다.
도 13(b)는 상기 공시 인덕터의 등가 직렬 인덕턴스(ESL)의 주파수 특성을 100[kHz]~100[MHz]과 동일한 주파수 대역에서 각 직류 중첩 전류(Idc)에 대해서 측정한 결과를 나타내는 그래프이다. 이 그래프의 가로축은 주파수[Hz], 세로축은 등가 직렬 인덕턴스(ESL)의 값[H]을 나타낸다. 또한, 파선으로 나타내는 주파수 특성(g0)은 직류 중첩 전류(Idc)=0(zero)[㎃]일 때의 주파수 특성이고, 실선으로 나타내는 주파수 특성(g1)은 직류 중첩 전류(Idc)=400[㎃], 일점 쇄선으로 나타내는 주파수 특성(g2)은 직류 중첩 전류(Idc)=800[㎃], 점으로 나타내는 주파수 특성(g3)은 직류 중첩 전류(Idc)=1100[㎃]일 때의 주파수 특성이다.
각 그래프에 나타내는 바와 같이, 직류 중첩 전류(Idc)의 증대에 따라, 등가 직렬 저항(ESR) 및 등가직렬 인덕턴스(ESL) 모두 값이 작아져서 변동한다.
도 14(a)는 도 12에 나타내는 등가 회로 모델을 이용하여 계산되는 직류 전류(Idc)의 무중첩 시의 공시 인덕터의 등가 직렬 저항(ESR)의 주파수 특성을, 도 13(a)에 나타내는 측정된 주파수 특성과 비교하여 나타내는 그래프이다. 이 그래프의 가로축은 주파수[Hz], 세로축은 등가 직렬 저항(ESR)의 값[Ω]을 나타낸다. 또한, 파선으로 나타내는 주파수 특성(f0)은 직류 중첩 전류(Idc)=0(zero)[㎃]일 때에 측정된 상술한 특성이고, 실선으로 나타내는 주파수 특성(F0)은 도 12에 나타내는 등가 회로 모델을 이용하여 계산된 직류 중첩 전류(Idc)=0(zero)[㎃]일 때의 특성이다.
도 14(b)는 도 12에 나타내는 등가 회로 모델을 이용하여 계산되는 직류 전류(Idc)의 무중첩 시의 공시 인덕터의 등가 직렬 인덕턴스(ESL)의 주파수 특성을, 도 13(b)에 나타내는 측정된 주파수 특성과 비교하여 나타내는 그래프이다. 이 그래프의 가로축은 주파수[Hz], 세로축은 등가 직렬 인덕턴스(ESL)의 값[H]을 나타낸다. 또한, 파선으로 나타내는 주파수 특성(g0)은 직류 중첩 전류(Idc)=0(zero)[㎃]일 때에 측정된 상술한 특성이고, 실선으로 나타내는 주파수 특성(G0)은 도 12에 나타내는 등가 회로 모델을 이용하여 계산된 직류 중첩 전류(Idc)=0(zero)[㎃]일 때의 특성이다.
각 그래프에 나타내는 바와 같이, 등가 직렬 저항(ESR) 및 등가 직렬 인덕턴스(ESL) 모두, 도 12에 나타내는 등가 회로 모델을 이용한 계산값은 실제의 측정값과 정밀도가 좋게 일치하고 있다.
도 15(a), (b), 도 16(a), (b)는, 도 12에 나타내는 수동 회로 소자(LX1, RX1, LX2, RX2)의 특성 변화율(kL1(Idc), kR1(Idc), kL2(Idc), kR2(Idc))을 나타내는 그래프이다. 각 그래프의 가로축은 공시 인덕터에 중첩시키는 직류 전류(Idc)(DCbias)[㎃]이고, 세로축은 특성 변화율이다. 상술한 바와 같이, 특성 변화율은 유도 소자(LX1, LX2)에 대해서는 (11)식에 나타내는 바와 같이, 직류 전류 무중첩 시의 수동 회로 소자(LX1, LX2)의 회로 정수(L1, L2)에 대한, 직류 전류 중첩 시의 회로 정수(LX1(Idc), LX2(Idc))의 비이다. 저항 소자(RX1, RX2)에 대해서는 (9)식에 나타내는 바와 같이, 직류 전류 무중첩 시의 수동 회로 소자(RX1, RX2)의 회로 정수(R1, R2)에 대한, 직류 전류 중첩 시의 회로 정수(RX1(Idc), RX2(Idc))의 비이다.
또한, 각 그래프에서의 사각형의 플롯(p)은 특성 변화율의 측정값이며, 각 플롯(p)을 잇는 특성선(H1, H2, H3, H4)은 이 측정값을 기초로 도출되는, 상술한 근사 함수 exp(f(x))로 표시된다.
도 17(a)는 도 12에 나타내는 등가 회로 모델에 기초하는, 제4 실시형태에서의 인덕터의 직류 전류 무중첩 시의 수동 등가 회로 모델, 도 17(b)는 마찬가지로 도 12에 나타내는 등가 회로 모델에 기초하는 제4 실시형태에서의 직류 전류 중첩 시의 인덕터의 비선형 등가 회로 모델을 나타내는 회로도이다. 한편, 도 17에 있어서 도 5 및 도 12와 동일 또는 상응하는 부분에는 동일부호를 붙이고 그 설명은 생략한다.
제4 실시형태에서의 도 17(a)에 나타내는 수동 등가 회로 모델에서는, 도 12에서는 회로 소자(RX1, RX2, LX1, LX2)로서 표시되어 있던 저항 소자(R), 유도 소자(L)가 직류 중첩 전류(Idc)에 따르지 않고, 정수라고 간주할 수 있는 수동 회로 소자(R1, R2, L1, L2)로서 표시되어 있다. 또한, 도 17(b)에 비선형 등가 회로 모델에서는, 상기의 유도 소자(L1)와 저항 소자(R1)의 직렬 회로에 직렬로 제어 전압원(B1)이 접속되고, 위 모델에서는 상기의 유도 소자(L1)와 저항 소자(R1)의 직렬 회로에 직렬로 제어 전압원(B1)이 접속되고, 상기의 유도 소자(L2)와 저항 소자(R2)와의 직렬 회로에 직렬로 제어 전압원(B1)과 동일한 제어 전압원(B2)이 접속되고, 도 12의 회로 소자(RX1, RX2, LX1, LX2)가 표시되어 있다. 이 제4 실시형태에서의 각 등가 회로 모델이라도, 도 5에 표시되는 제2 실시형태에서의 각 등가 회로 모델과 마찬가지로 제어 전압원(B1)과 수동 회로 소자(L1, R1)의 직렬 회로와 제어 전압원(B2)과 수동 회로 소자(L2, R2)의 직렬 회로가 복수 병렬로 접속되어 있다. 또한, 도 9에 표시되는 제3 실시형태에서의 등가 회로 모델과 마찬가지로 인덕터에 대한 직류 전류(Idc)의 중첩에 의해 특성이 변화하지 않는 수동 회로 소자(r1, c1, r2, c2)를 포함하여 구성되어 있다.
이 제4 실시형태에 있어서도, 도 17(a)에 나타내는 수동 등가 회로 모델에서는 직류 전류(Idc)가 중첩하고 있지 않은 교류 전류(IaC)가, 전류원 모델(I0)에 의해 회로에 전류(I)로서 흐른다. 또한, 도 17(b)에 나타내는 비선형 등가 회로 모델에서는 직류 전류(Idc)가 중첩한 교류 전류(Iac)가, 전류원 모델(I1)에 의해 회로에 전류(I)로서 흐른다.
도 17에 나타내는 등가 회로 모델을 이용한 인덕터의 시뮬레이션에 있어서도, 우선 수동 회로 소자(L1, R1)의 특성 변화율(kL1(Idc), kR1(Idc)) 및 무중첩 시 전압(VL1, VR1)에 기초하여, 제1 실시형태의 시뮬레이션 방법과 동일하게 해서, 중첩 시 전압(VL1(Idc), VR1(Idc))과 무중첩 시 전압(VL1, VR1)의 차분 전압(ΔVL1(Idc), ΔVR1(Idc))을 발생시킨다. 또한, 제어 전압원(B2)은, 산출한 수동 회로 소자(L2, R2)의 특성 변화율(kL2(Idc), kR2(Idc)) 및 무중첩 시 전압(VL2, VR2)에 기초하여, 제어 전압원(B1)과 동일하게 해서, 중첩 시 전압(VL2(Idc), VR2(Idc))과 무중첩 시 전압(VL2, VR2)의 차분 전압(ΔVL2(Idc), ΔVR2(Idc))을 발생시킨다.
도 18(a)는 도 17에 나타내는 등가 회로 모델을 이용하여 계산되는, 각 직류 전류(Idc)의 중첩 시 및 무중첩 시의 공시 인덕터의 등가 직렬 저항(ESR)의 주파수 특성을, 도 13(a)에 나타내는 측정된 주파수 특성과 비교하여 나타내는 그래프이다. 이 그래프의 가로축은 주파수[Hz], 세로축은 등가 직렬 저항(ESR)의 값[Ω]을 나타낸다.
또한, 파선으로 나타내는 주파수 특성(f0, f1, f2, f3)은 직류 중첩 전류 Idc=0(zero), 400, 800, 1100[㎃]일 때에 측정된, 상술한 도 13(a)에 나타내는 특성이다. 또한, 실선으로 나타내는 주파수 특성(F0)은 도 17(b)에 나타내는 비선형 등가 회로 모델을 이용하여, 직류 중첩 전류 Idc=0(zero)[㎃]일 때에 계산된 특성이고, 실선으로 나타내는 주파수 특성(F1, F2, F3)은 직류 중첩 전류 Idc=400, 800, 1100[㎃]일 때에 계산된 특성이다. 도 13(a)에 나타내는 주파수 특성(f0~f3)의 측정 주파수 대역은 상술한 바와 같이, 100[kHz]~100[MHz]이지만, 이 시뮬레이션에서의 계산 주파수 대역은 100[kHz]~1000[MHz]이다.
도 18(b)는 도 17에 나타내는 등가 회로 모델을 이용하여 계산되는, 각 직류 전류(Idc)의 중첩 시 및 무중첩 시의 공시 인덕터의 등가 직렬 인덕턴스(ESL)의 주파수 특성을 도 13(b)에 나타내는 측정된 주파수 특성과 비교하여 나타내는 그래프이다. 이 그래프의 가로축은 주파수[Hz], 세로축은 등가 직렬 인덕턴스(ESL)의 값[H]을 나타낸다.
또한, 파선으로 나타내는 주파수 특성(g0, g1, g2, g3)은, 직류 중첩 전류Idc=0(zero), 400, 800, 1100[㎃]일 때에 측정된, 전술한 도 13 (b)에 나타내는 특성이다. 또한, 실선으로 나타내는 주파수 특성(G0)은, 도 17 (b)에 나타내는 비선형 등가 회로 모델을 이용하고, 직류 중첩 전류 Idc=0(zero)[㎃]일 때에 계산된 특성, 실선으로 나타내는 주파수 특성(G1, G2, G3)은, 직류 중첩 전류 Idc=400, 800, 1100[㎃]일 때에 계산된 특성이다. 도 13 (b)에 나타내는 주파수 특성(g0~g3)의 측정 주파수 대역은 전술한 바와 같이, 100[kHz]~100[MHz]이지만, 이 시뮬레이션에 서의 계산 주파수 대역도 100[kHz]~1000[MHz]이다.
각 그래프에 나타내는 바와 같이, 등가 직렬 저항(ESR) 및 등가직렬 인덕턴스(ESL) 모두, 도 17에 나타내는 등가 회로 모델을 이용한 계산값은 실제 측정값과 정밀도가 좋게 일치하고 있다.
도 19는 등가 회로 모델이 도 17에 표시되는 공시 인덕터에 대한, 인가 전압에 대한 전류의 과도 응답 특성을 나타내는 그래프이다.
도 19(a)에 나타내는 그래프는 공시 인덕터에 인가한, 스위칭 레귤레이터IC(고 집적화 회로)의 출력 전압 파형을 나타낸다. 이 그래프의 가로축은 시간[㎲], 세로축은 전압[V]을 나타낸다. 또한, 실선으로 나타내는 전압 파형(Va)은 측정값으로부터 얻어진 파형이며, 파선으로 나타내는 전압 파형(Vb)은 사용한 스위칭 레귤레이터IC의 제조 회사가 제공하는 디바이스 모델을 사용하여 계산한 계산값으로부터 얻어진 파형이다.
도 19(b)에 나타내는 그래프는 도 19(a)에 나타내는 전압이 인가됨으로써 공시 인덕터에 흐르는 전류의 과도 응답 특성을 나타내는 전류 파형이다. 이 그래프의 가로축은 시간[㎲], 세로축은 전류[A]를 나타낸다. 또한, 실선으로 나타내는 전류 파형(Ia)은 측정값으로부터 얻어진 파형이며, 파선으로 나타내는 전류 파형(Ib)은 도 17에 나타내는 공시 인덕터의 등가 회로 모델을 사용하여 계산한 계산값으로부터 얻어진 파형이다.
도 19(b)의 그래프에 나타내는 바와 같이, 도 17에 나타내는 공시 인덕터의 등가 회로 모델을 사용한 과도 응답 특성의 계산값은, 실제 측정값과 정밀도가 좋게 일치하고 있다.
즉, 제4 실시형태의 도 17에 나타내는 공시 인덕터의 등가 회로 모델을 이용한 시뮬레이션에 의해서도, 직류 중첩 전류(Idc)에 따라 변동하는 공시 인덕터의 특성을 광대역이면서 고정밀도로 재현해서 시뮬레이트 할 수 있고, 제1부터 제3의 각 실시형태와 동일한 작용 효과가 발휘된다.
또한, 상기의 각 실시형태에서는, 각 등가 회로 모델을 리니어 테크놀로지사가 제공하는 LTspice에 적용한 경우에 대해서 설명했지만, 적용하는 회로 시뮬레이터는 이에 한정되지는 않는다. 예를 들면, Applied Wave Reserch사(AWR사)가 제공하는 Microwave Office나, Agilent Technologies Inc.사(Agilent사)가 제공하는 ADS와 같은 다른 회로 시뮬레이터에도 동일하게 적용할 수 있다.
상술한 각 실시형태의 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델은, 다음 컴퓨터 프로그램을 이용함으로써 간단히 이용할 수 있다. 이 컴퓨터 프로그램은 제1 스텝과 제2 스텝과 제3 스텝을 포함한다. 제1 스텝에서는, 전자회로 설계에 이용하는 인덕터의 종류를 입력한다. 제2 스텝에서는, 인덕터에 인가하는 전압(V) 또는 인덕터에 흐르는 전류(I)를 입력한다. 제3 스텝에서는, 제2 스텝에서 입력된 전압(V) 또는 전류(I)에 의해 인덕터에 흐르는 전류(I)를 계측하고, 참조 전류(Iref)를 참조한다. 그리고 제1 스텝에서 입력된 종류의 인덕터에 대해서 미리 준비된 근사 함수 exp(f(x))에 의해 참조 전류(Iref)에 대응해서 산출되는 특성 변화율(kRX(Iref), kCX(Iref), kLX(Iref)) 및 무중첩 시 전압(VRX, VCX, VLX)에 기초하여, 제어 전압원(BX)에 의해서 참조 전류(Iref)에 따라 차분 전압(ΔVRX(Iref), ΔVCX(Iref), ΔVLX(Iref))을 발생시켜, 무중첩 시 전압(VRX, VCX, VLX)에 중첩시킴으로써, 인덕터의 직류 전류 중첩 시의 비선형 특성을 시뮬레이트 한다. 컴퓨터 프로그램은 이들 각 스텝을 실행하는 연산 처리에 의해, 상기의 각 실시형태의 인덕터 시뮬레이션 방법을 실시하는 또는 상기의 각 실시형태의 인덕터 비선형 등가 회로 모델을 기능시킨다.
본 구성에 의하면, 시뮬레이트하는 인덕터의 종류 및 인덕터에 인가하는 전압(V) 또는 인덕터에 유실시키는 전류(I)의 값이 컴퓨터 프로그램에 입력됨으로써, 입력된 종류의 인덕터의 비선형 특성은 컴퓨터 프로그램에 의해, 수동 회로 소자(Rx, Cx, Lx)의 무중첩 시 전압(VRX, VCX, VLX)에 차분 전압(ΔVRX(Iref), ΔVCX(Iref), ΔVLX(Iref))이 중첩되어 자동적으로 시뮬레이션된다. 이렇기 때문에, 본 시뮬레이션 방법 또는 본 비선형 등가 회로 모델의 이용자는 시뮬레이트하는 인덕터의 종류 및 인덕터에 인가하는 전압(V) 또는 인덕터에 흐르는 전류(I)의 값을 컴퓨터 프로그램에 입력하는 것만으로, 정확한 회로 시뮬레이션을 고정밀도이면서 간단히 실시할 수 있다. 이 결과, 회로 시뮬레이트에 관한 전문지식을 가지지 않는 일반 이용자여도 인덕터를 이용한 전자회로의 정확한 회로 시뮬레이션을 고정밀도이면서 간단히 실시할 수 있다.
또한, 상기의 컴퓨터 프로그램은, 상기 컴퓨터 프로그램을 포함하는 전자부품 제조 메이커 등의 서버에 인터넷망을 통해서 액세스함으로써 인터넷망에 접속된 퍼스널 컴퓨터 등의 단말로부터 사용할 수 있다. 본 구성에 의하면, 이용자는, 인터넷망에 접속된 단말로부터 상기 컴퓨터 프로그램을 포함하는 서버에 액세스함으로써 상기 컴퓨터 프로그램을 용이하게 사용할 수 있다. 이렇기 때문에, 상기의 각 실시형태에 의한 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델을 다수의 이용자에게 제공할 수 있게 된다.
Rx, Lx, Cx, r, l, c: 수동 회로 소자
BX: 제어 전압원(비헤이비어 전압원 모델)
VX: 전압원 모델(전류계)
IX: 전류원 모델
RX1(Idc): 가변 저항 소자
LX1(Idc): 가변 유도 소자

Claims (11)

  1. 수동 회로 소자를 사용하여 나타낸 인덕터의 등가 회로를 이용하여 인덕터의 특성을 컴퓨터에 의해 시뮬레이트하는 인덕터의 시뮬레이션 방법에 있어서,
    컴퓨터는,
    인덕터에 흐르는 전류를 참조하여,
    참조한 전류에 대응하여, 실측값을 기초로 전류를 변수로서 나타낸 근사 함수에 의해 산출되는 직류 전류 중첩 시에서의 상기 수동 회로 소자의 특성 변화율, 및 직류 전류 무중첩 시의 상기 수동 회로 소자의 특성에 따라 직류 전류 중첩 시에 상기 수동 회로 소자에 생기는 무중첩 시 전압에 기초하여, 직류 전류의 중첩에 의해 특성이 변화하는 상기 수동 회로 소자에 직렬로 접속되는 제어 전압원에 의해, 직류 전류 중첩 시의 상기 수동 회로 소자의 특성에 따라 직류 전류 중첩 시에 상기 수동 회로 소자에 생기는 중첩 시 전압과 상기 무중첩 시 전압의 차분 전압을 발생시켜, 상기 무중첩 시 전압에 상기 차분 전압을 중첩시킴으로써,
    인덕터의 직류 전류 중첩 시의 비선형 특성을 시뮬레이트하는 것을 특징으로 하는 인덕터의 시뮬레이션 방법.
  2. 인덕터의 등가 회로를 나타내는 수동 회로 소자와,
    상기 인덕터에 흐르는 전류를 참조하는 전류 참조 수단과,
    직류 전류 중첩 시에서의 상기 수동 회로 소자의 특성 변화율을, 실측값을 기초로, 전류를 변수로서 나타내는 근사 함수에 의해, 상기 전류 참조 수단에 의해 참조되는 전류에 대응하여 산출되는 상기 특성 변화율, 및 직류 전류 무중첩 시의 상기 수동 회로 소자의 특성에 따라 직류 전류 중첩 시에 상기 수동 회로 소자에 생기는 무중첩 시 전압에 기초하여, 직류 전류 중첩 시의 상기 수동 회로 소자의 특성에 따라 직류 전류 중첩 시에 상기 수동 회로 소자에 생기는 중첩 시 전압과 상기 무중첩 시 전압의 차분 전압을 발생시키는, 직류 전류의 중첩에 의해 특성이 변화하는 상기 수동 회로 소자에 직렬로 접속되는 제어 전압원을 포함하여 구성되는 인덕터의 비선형 등가 회로 모델을 이용하여 인덕터의 직류 전류 중첩 시의 비선형 특성을 시뮬레이트하는 것을 특징으로 하는 인덕터의 시뮬레이션 방법.
  3. 제1항 또는 제2항에 있어서,
    상기 차분 전압은, 상기 차분 전압을 ΔV, 상기 무중첩 시 전압을 V0, 상기 근사 함수를 참조하는 전류 x를 변수로 하는 함수 exp(f(x))로 한 경우에, 다음 식
    ΔV=V0×[exp(f(x))-1]
    의 함수형으로 주어지는 것을 특징으로 하는 인덕터의 시뮬레이션 방법.
  4. 제1항 또는 제2항에 있어서,
    상기 근사 함수는 홀수차의 거듭제곱항을 포함하지 않는 다항식 형식의 우함수로 주어지는 것을 특징으로 하는 인덕터의 시뮬레이션 방법.
  5. 제1항 또는 제2항에 있어서,
    상기 인덕터에 흐르는 전류는 상기 등가 회로의 입력단 또는 출력단에서 참조되고, 상기 무중첩 시 전압은 상기 수동 회로 소자의 양단(兩端)에서 참조되는 것을 특징으로 하는 인덕터의 시뮬레이션 방법.
  6. 제1항 또는 제2항에 있어서,
    상기 제어 전압원에 직렬로 접속되는 상기 수동 회로 소자는, 유도 소자 단체, 또는 유도 소자와 저항 소자의 직렬 회로, 또는 유도 소자와 저항 소자와 용량 소자의 직렬 회로인 것을 특징으로 하는 인덕터의 시뮬레이션 방법.
  7. 제1항 또는 제2항에 있어서,
    상기 제어 전압원과 상기 수동 회로 소자의 직렬 회로가 복수 병렬로 접속되는 것을 특징으로 하는 인덕터의 시뮬레이션 방법.
  8. 제1항 또는 제2항에 있어서,
    상기 등가 회로는 직류 전류의 중첩에 의해 특성이 변화하지 않는 상기 수동 회로 소자를 포함하여 구성되는 것을 특징으로 하는 인덕터의 시뮬레이션 방법.
  9. 인덕터의 등가 회로를 나타내는 수동 회로 소자와,
    상기 인덕터에 흐르는 전류를 참조하는 전류 참조 수단과,
    직류 전류 중첩 시에서의 상기 수동 회로 소자의 특성 변화율을, 실측값을 기초로, 전류를 변수로 하여 나타내는 근사 함수에 의해, 상기 전류 참조 수단에 의해 참조되는 전류에 대응하여 산출되는 상기 특성 변화율, 및 직류 전류 무중첩 시의 상기 수동 회로 소자의 특성에 따라 직류 전류 중첩 시에 상기 수동 회로 소자에 생기는 무중첩 시 전압에 기초하여, 직류 전류 중첩 시의 상기 수동 회로 소자의 특성에 따라 직류 전류 중첩 시에 상기 수동 회로 소자에 생기는 중첩 시 전압과 상기 무중첩 시 전압의 차분 전압을 발생시키는, 직류 전류의 중첩에 의해 특성이 변화하는 상기 수동 회로 소자에 직렬로 접속되는 제어 전압원을 포함하여 구성되는 인덕터의 비선형 등가 회로 모델을 격납하는 격납 수단과,
    상기 인덕터의 종류를 입력하는 제1 스텝,
    상기 인덕터에 인가하는 전압 또는 상기 인덕터에 흘려보내는 전류를 입력하는 제2 스텝,
    상기 제2 스텝에서 입력된 전압 또는 전류에 의해 상기 인덕터에 흐르는 전류를 참조하고, 상기 제1 스텝에서 입력된 종류의 상기 인덕터에 대해서 미리 준비된 상기 근사 함수에 의해 참조 전류에 대응하여 산출되는 상기 특성 변화율 및 상기 무중첩 시 전압에 기초하여, 상기 제어 전압원에 의해 상기 차분 전압을 발생시켜, 상기 무중첩 시 전압에 상기 차분 전압을 중첩시킴으로써, 상기 인덕터의 직류 전류 중첩 시의 비선형 특성을 시뮬레이트하는 제3 스텝을 포함하는 컴퓨터 프로그램을 컴퓨터에 실행시켜서 제1항 또는 제2항에 기재된 인덕터의 시뮬레이션 방법을 실행하는 실행 수단으로부터 구성되는 것을 특징으로 하는 인덕터의 시뮬레이션 장치.
  10. 상기 격납 수단 및 상기 컴퓨터 프로그램을 포함하는 서버에 인터넷망을 통해서 액세스하고, 인터넷망에 접속된 단말로부터 상기 실행 수단에 컴퓨터 프로그램을 실행시키는 것을 특징으로 하는 제9항에 기재된 인덕터의 시뮬레이션 장치의 사용 방법.
  11. 삭제
KR1020157032574A 2013-05-14 2014-05-02 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델 KR101616036B1 (ko)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JPJP-P-2013-102542 2013-05-14
JP2013102542 2013-05-14
PCT/JP2014/062157 WO2014185294A1 (ja) 2013-05-14 2014-05-02 インダクタのシミュレーション方法およびインダクタの非線形等価回路モデル

Publications (2)

Publication Number Publication Date
KR20150138396A KR20150138396A (ko) 2015-12-09
KR101616036B1 true KR101616036B1 (ko) 2016-04-27

Family

ID=51898273

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020157032574A KR101616036B1 (ko) 2013-05-14 2014-05-02 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델

Country Status (5)

Country Link
US (1) US10521533B2 (ko)
JP (1) JP5773102B2 (ko)
KR (1) KR101616036B1 (ko)
CN (1) CN105229645B (ko)
WO (1) WO2014185294A1 (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20190067016A (ko) 2017-12-06 2019-06-14 삼성전기주식회사 인덕터 등가 회로를 저장하는 저장 장치 및 인덕터 등가 회로를 제공하는 서버

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6245187B2 (ja) * 2015-02-02 2017-12-13 株式会社村田製作所 パワーインダクタの評価装置、及び、パワーインダクタの評価プログラム
JP6657822B2 (ja) * 2015-11-13 2020-03-04 株式会社村田製作所 インダクタのシミュレーションモデル
WO2018061016A1 (en) * 2016-10-02 2018-04-05 Ariel Scientific Innovations Ltd. Active electronic emulation of a passive circuit component
JP7043178B2 (ja) * 2017-03-23 2022-03-29 太陽誘電株式会社 受動素子の等価回路のシミュレーション方法及びその装置
CN107404241B (zh) * 2017-07-04 2019-06-18 华中科技大学 一种非线性电感分析方法和三相pwm变流器控制方法
US10755015B2 (en) * 2017-08-21 2020-08-25 Semiconductor Components Industries, Llc Agnostic model of semiconductor devices and related methods
CN108490344B (zh) * 2018-03-28 2020-05-22 电子科技大学 基于训练样本优化的模拟电路故障诊断方法
US11501049B1 (en) * 2018-09-28 2022-11-15 Cadence Design Systems, Inc. Systems and methods for modeling interactions of power and signals in a multi-layered electronic structure
JP7107339B2 (ja) * 2019-09-30 2022-07-27 株式会社村田製作所 非線形特性計算の方法、非線形特性計算プログラムとその使用方法、並びに記録媒体
CN115980493B (zh) * 2023-01-03 2023-11-07 广州市德珑电子器件有限公司 多电感的光伏逆变器测试方法、装置、设备以及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006163606A (ja) 2004-12-03 2006-06-22 Nec Toppan Circuit Solutions Inc 回路設計支援方法、回路設計支援システム、装置及びプログラム
JP2012003409A (ja) 2010-06-15 2012-01-05 Taiyo Yuden Co Ltd 積層チップコンデンサの等価回路モデル,その回路定数解析方法・プログラム・装置,及び回路シミュレータ
WO2012090602A1 (ja) 2010-12-29 2012-07-05 株式会社村田製作所 等価回路モデル,プログラム及び記録媒体

Family Cites Families (32)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4485356A (en) * 1982-06-04 1984-11-27 Lorenzo Fassino Variable low frequency sinusoidal oscillator using simulated inductance
JPH0714916Y2 (ja) * 1988-03-04 1995-04-10 横河・ヒユーレット・パツカード株式会社 回路定数測定装置
JPH02128573U (ko) 1989-03-29 1990-10-23
US5343177A (en) * 1992-12-01 1994-08-30 Hughes Aircraft Company Simulated coupled inductors using generalized impedance converters
US6064808A (en) * 1997-08-01 2000-05-16 Lucent Technologies Inc. Method and apparatus for designing interconnections and passive components in integrated circuits and equivalent structures by efficient parameter extraction
JP3250503B2 (ja) * 1997-11-11 2002-01-28 株式会社村田製作所 可変インダクタ素子
US6324493B1 (en) * 1999-04-01 2001-11-27 Agere Systems Guardian Corp. Method and apparatus for modeling electromagnetic interactions in electrical circuit metalizations to simulate their electrical characteristics
US6633154B1 (en) * 2000-01-04 2003-10-14 William B. Duff, Jr. Method and circuit for using polarized device in AC applications
JP3435636B2 (ja) * 2000-04-06 2003-08-11 株式会社村田製作所 可変インダクタンス素子
JP4507421B2 (ja) 2001-02-27 2010-07-21 パナソニック株式会社 受動素子の等価回路モデル導出方法、シミュレータ、及び記憶媒体
US8712544B2 (en) * 2001-04-13 2014-04-29 Greatbatch Ltd. Electromagnetic shield for a passive electronic component in an active medical device implantable lead
JP2003122810A (ja) 2001-10-12 2003-04-25 Matsushita Electric Ind Co Ltd シミュレーション方法
JP4391263B2 (ja) * 2004-02-20 2009-12-24 Okiセミコンダクタ株式会社 半導体素子、その製造方法及びその半導体素子を用いた高周波集積回路
US7562000B2 (en) * 2004-05-04 2009-07-14 Lorentz Solution, Inc. Electromagnetic simulator systems and methods
US20060217948A1 (en) * 2005-03-23 2006-09-28 Tdk Corporation Component for a simulation tool
US7546230B2 (en) * 2005-04-18 2009-06-09 International Business Machines Corporation Electromagnetic response model with improved high frequency stability
JP2007207168A (ja) 2006-02-06 2007-08-16 Nec Electronics Corp Emiシミュレーションモデル、emiシミュレーションシステムと方法
TWI319581B (en) * 2006-08-08 2010-01-11 Murata Manufacturing Co Laminated coil component and method for manufacturing the same
KR100957446B1 (ko) * 2007-12-24 2010-05-11 순천향대학교 산학협력단 3차원 구조의 직렬 l-c 공진기 및 이를 이용한광대역통과여파기
WO2010038799A1 (ja) 2008-09-30 2010-04-08 日立金属株式会社 インダクタンス素子の直流重畳特性の解析方法及び電磁界シミュレータ
US8610311B1 (en) * 2009-01-27 2013-12-17 Energy Magnification Corporation Passive power generation system
JP5121757B2 (ja) * 2009-03-02 2013-01-16 太陽誘電株式会社 積層チップインダクタの等価回路モデルの回路定数解析方法及び回路シミュレーション方法
JP5084801B2 (ja) * 2009-08-31 2012-11-28 株式会社村田製作所 インダクタおよびdc−dcコンバータ
US8260435B2 (en) * 2010-03-17 2012-09-04 Greatbatch Ltd. Implantable lead for an active medical device having an inductor design minimizing eddy current losses
CN102236728B (zh) * 2010-04-30 2013-08-07 国际商业机器公司 一种集成电路设计方法和设计仿真系统
JP2012150579A (ja) 2011-01-17 2012-08-09 Taiyo Yuden Co Ltd コンデンサの回路シミュレーションモデル及びその構築方法,回路シミュレーション方法及び回路シミュレータ
WO2012131863A1 (ja) * 2011-03-28 2012-10-04 富士通株式会社 設計支援方法、設計支援プログラムおよび設計支援装置
JP2012231595A (ja) * 2011-04-26 2012-11-22 Canon Inc 振動装置の駆動回路、塵埃除去装置及び振動型アクチュエータにおける振動装置の駆動回路
WO2012173147A1 (ja) * 2011-06-15 2012-12-20 株式会社 村田製作所 積層コイル部品、及び該積層コイル部品の製造方法
JP6012960B2 (ja) * 2011-12-15 2016-10-25 太陽誘電株式会社 コイル型電子部品
US9595893B2 (en) * 2012-07-23 2017-03-14 Schlumberger Technology Corporation Non-stationary multi-frequency vibration energy harvesting with tunable electrical impedance
JP6148630B2 (ja) 2014-02-13 2017-06-14 ブリヂストンフローテック株式会社 加締装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006163606A (ja) 2004-12-03 2006-06-22 Nec Toppan Circuit Solutions Inc 回路設計支援方法、回路設計支援システム、装置及びプログラム
JP2012003409A (ja) 2010-06-15 2012-01-05 Taiyo Yuden Co Ltd 積層チップコンデンサの等価回路モデル,その回路定数解析方法・プログラム・装置,及び回路シミュレータ
WO2012090602A1 (ja) 2010-12-29 2012-07-05 株式会社村田製作所 等価回路モデル,プログラム及び記録媒体

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Takeshi Horigome, "An actual waveform is reproduced by the electronic circuit simulator LTspice the 4nd time transistor technology", CQ publication CO., LTD. 48th volume No.10(pp.175-1822), 2011.10.01

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20190067016A (ko) 2017-12-06 2019-06-14 삼성전기주식회사 인덕터 등가 회로를 저장하는 저장 장치 및 인덕터 등가 회로를 제공하는 서버

Also Published As

Publication number Publication date
CN105229645A (zh) 2016-01-06
JP5773102B2 (ja) 2015-09-02
JPWO2014185294A1 (ja) 2017-02-23
US20160070837A1 (en) 2016-03-10
CN105229645B (zh) 2018-05-15
US10521533B2 (en) 2019-12-31
WO2014185294A1 (ja) 2014-11-20
KR20150138396A (ko) 2015-12-09

Similar Documents

Publication Publication Date Title
KR101616036B1 (ko) 인덕터의 시뮬레이션 방법 및 인덕터의 비선형 등가 회로 모델
KR101616037B1 (ko) 콘덴서의 시뮬레이션 방법 및 콘덴서의 비선형 등가 회로 모델
US9563728B2 (en) Equivalent circuit model, program, and recording medium
US8718987B2 (en) Circuit simulation model of capacitor, constructing method of simulation model, method of circuit simulation, circuit simulator
JP2010204869A (ja) 積層チップインダクタの等価回路モデルの回路定数解析方法及び回路シミュレーション方法
JP5861774B2 (ja) コンデンサの等価回路モデルの導出方法
WO2003023661A2 (en) Method and apparatus for analysing and modeling of analog systems
US20180096088A1 (en) Method of constructing and method of simulating equivalent circuit for capacitor, and simulation device therefor
Elg et al. Optimization of the design of a wideband 1000 kV resistive reference divider
CN105453088B (zh) 考虑到施加交流电压的电容器的静电电容值决定方法及程序
Verma Practical simulation and modelling of lightning impulse voltage generator using marx circuit
Onikienko et al. Conductive EMI of class D audio amplifiers prediction system
JP2017091346A (ja) インダクタのシミュレーションモデル
JP2018160132A (ja) 等価回路の構築方法,シミュレーション方法及びその装置
Ekman et al. On characterizing artifacts observed in PEEC based modeling
Maffezzoni et al. Steady-state analysis of strongly nonlinear oscillators by means of runge-kutta methods
TWI383162B (zh) Fault location method
CN115438461A (zh) 一种动态相量表征的直流电力系统随机暂态分析方法
Abdulahovic et al. Improvement of admittance measurements for modeling of an ungrounded reactor/transformer winding

Legal Events

Date Code Title Description
A201 Request for examination
A302 Request for accelerated examination
E701 Decision to grant or registration of patent right
GRNT Written decision to grant