The invention relates to the field of Electronic Design Automation (EDA) software for designing inductors, transformers and related circuits particularly for integrated circuits (ICs).
In electrical and electronic circuit engineering, a transformer is a versatile and ubiquitous component that accomplishes impedance level transformation between its primary and secondary windings through magnetic coupling, while providing direct current (DC) isolation between the primary and secondary windings. These properties are useful today in the field of IC design, and especially for radio frequency (RF) IC design, because inductive components such as spiral coils functioning as primary and secondary windings can be fabricated with good performance characteristics using known semiconductor fabrication techniques.
What is needed then are methods and systems for improving and/or simplifying the design of components such as integrated spiral transformers for use in ICs fabricated using modern semiconductor fabrication techniques that take into account the different functionality of the designed-for transformers and the interaction of such transformers with other components of an IC.
Embodiments of the invention also provide a method of optimizing the design of a transformer or circuit for an integrated circuit. In an embodiment, a method for designing a transformer in an integrated circuit comprises receiving one or more desired characteristics of the transformer from user input and iteratively determining a design solution for the transformer through one or more simulations and modifications using a rule-set. The method combines the one or more desired characteristics with other preset characteristics of the transformer or the integrated circuit. A first model of the transformer is defined with typical load impedances and simulated having the combined characteristics to determine performance. Results of the simulation are processed to calculate performance with the load impedances specified by the user. The results are further processed to obtain a mathematical model that includes tuning capacitors. The first and subsequent models are modified by drawing on a rule-set of expert knowledge relating to general dependency of at least one design criterion, such as a physical, geometrical or performance characteristic, with another design criterion.
In an embodiment, the method further comprises receiving one or more optimization targets from user input. A constrained optimization problem is solved to determine optimal values of the tuning capacitors based on the one or more desired characteristics and the one or more optimization targets. In an embodiment, the one or more optimization targets include lowest possible insertion loss.
In an embodiment, the method is usable to design a transformer including a plurality of taps connecting the primary coil to either a fixed potential, a signal, or a load. In an embodiment, the transformer includes a plurality of taps connecting the secondary coil to either a fixed potential, a signal, or a load.
In an embodiment, an apparatus for designing a transformer in an integrated circuit, the apparatus comprises a processor and a user interface adapted to allow a user to input one or more desired characteristics of the transformer. The one or more desired characteristics can include load impedances at primary and secondary coils of the transformer, and the processor is configured to determine sufficient optimal characteristics of the transformer to define the transformer for fabrication in the integrated circuit by executing embodiment of methods for designing a transformer in an integrated circuit comprises receiving one or more desired characteristics of the transformer from user input and iteratively determining a design solution for the transformer through one or more simulations and modifications using a rule-set.
Embodiments of the invention will now be described, by way of example only, with reference to the accompanying drawings in which:
As disclosed in U.S. application Ser. No. 12/575,410 entitled “EXPERT SYSTEM BASED INTEGRATED INDUCTOR SYNTHESIS AND OPTIMIZATION” by Bantas et al., incorporated herein by reference, the physical design of an integrated electrical component such as a spiral inductor can be effectively determined using a computer implemented method. The disclosed method iteratively uses electrical circuit simulation to determine the performance of a candidate physical component and draws on a rule-set of expert knowledge to enable modification of its physical properties, until a set of user-provided objectives and constraints are met. The model of a candidate physical design is preferably provided by an efficient electromagnetics (EM) model extraction engine. The disclosed method brings several advantages such as rapid execution and a high degree of computer-assisted automation; it has been incorporated in the VeloceRF™ EDA product offered by Helic, Inc. and is being used by several semiconductor design companies in production.
Having efficiently resolved the problem of designing integrated components such as spiral inductors, in the invention disclosed herein we concentrate on the related problem of designing integrated transformers.
Referring to
A. Differential to differential.
B. Differential with center tap, to differential.
C. Differential with center tap, to differential with center tap.
D. Differential to single-ended, sometimes referred-to as a ‘balun’.
E. Differential with center tap, to single-ended.
The problem of optimally designing an integrated transformer presents some difficulties, which are known to those skilled in the art, such as:
The method and apparatus disclosed in U.S. application Ser. No. 12/575,410 helps to overcome some of the limitations described above. However, these need to be adapted and extended for tackling the problem of integrated transformers. The additional steps necessary are disclosed herein.
As disclosed in U.S. application Ser. No. 12/575,410, an expert system-based inductor synthesis system can comprise the following modules:
Specific functions need to be introduced at the Model Extraction and Simulation steps to enable optimal transformer design according to this invention.
We first define a set of formulae that are used in the invention:
Scattering parameters (S-parameters)
Let T be a given transformer instance. By performing S-parameter analysis for a set of M frequencies:
F={f(1), . . . ,f(M)}
we obtain a set s of S-parameter matrices:
for each i in {1, . . . , M}.
Impedance Matrix
For any given S-parameter matrix S(i), we can obtain the corresponding impedance matrix,
from the known formulae (as taught in “S-Parameter Techniques”, Application Note 95-1, by Agilent Technologies, incorporated herein by reference):
Z11=Z0·((1+S11)·(1−S22)+S12·S21)/det(I−S)
Z12=Z0·2·S12/det(I−S)
Z21=Z0·2·S21/det(I−S)
Z22=Z0·((1−S11)·(1+S22)+S12·S21)/det(I−S)
where z0 the characteristic impedance (e.g. 50 Ohm).
Inductance
If z(i) represents the impedance matrix of a transformer, the primary and secondary inductances can be calculated by:
L1(i)=Im(Z11(i))/(2·π·f(i))
L2(i)=Im(Z22(i))/(2·π·f(i))
k(i)=Im(Z12(i))/sqrt(Im(Z11(i))·Im(Z22(i)))
where index 1 is assumed to correspond to the primary coil port of the transformer, and index 2 to the secondary coil port.
The resonance frequency Fres of the transformer, such that,
f(k)<=Fres<f(k+1)
L1(k)>0 and L2(k)>0
L1(k+1)<0 or L2(k+1)<0
can be calculated by linear interpolation, after finding the first value of i for which L1(i)L1(i+1)<0, or L2(i)L2(i+1)<0.
Quality Factor
If z(i) represents the impedance matrix of a transformer, the primary and secondary quality factor can be calculated by:
QI(i)=Im(Z11(i))/Re(Z11(i))
Q2(i)=Im(Z22(i))/Re(Z22(i))
The maximum quality factor QImax, Q2max and the corresponding frequencies F1max, F2max can also be calculated by evaluating the discrete functions Q1(i), Q2(i) for each sampled frequency f(1)<=f(i)<Fres.
Bandwidth
We define the bandwidth of a transformer in terms of its insertion loss S12. Let S12_ref be the insertion loss at a given frequency f_ref. Given a tolerance value S12_tol (e.g. −1 dB or ˜3 dB relative to s12_ref) the bandwidth can then be described as a frequency interval:
BW=[f_ref−f_BW1,f_ref+f_BW2]
such that:
S12(f_ref−f_BW1)=S12_ref−S12_tol
S12(f_ref+f_BW2)=S12_ref−S12_tol
and:
S12(f)>S12_ref−S12_tol
for all f in BW.
The bandwidth of a transformer at given frequency f(i) can be numerically approximated as follows:
Arbitrary Loads
The S-parameter matrix S(i) is measured with both ports connected to reference loads of Z0=50 ohm. If the circuit is connected to arbitrary complex loads Zp, Zs, its electrical behavior is described by a transformation of S(i) as described in the book “Microwave Amplifiers and Oscillators” by Christian Gentili, incorporated herein by reference:
where:
S11′=((1−Gs·S22)·(S11−Conj(Gp))+S12·S21·Gs·A11)/DS
S12′=(S12·(1.0−|Gp|^2)·A21)/DS
S21′=(S21·(1.0−|Gs|^2)·A12)/DS
S22′=((1−Gp·S11)·(S22−Conj(Gs))+S12·S21·Gp·A22)/DS
and where the reflection coefficients are defined as:
Gp=(Zp−Z0)/(Zp+Z0)
Gs=(Zs−Z0)/(Zs+Z0)
and the remaining matrices are defined as:
Ap=((1−Conj(Gp))·Sqrt(1.0−|Gp|^2)))|1−Gp|
As=((1−Conj(Gs))·Sqrt(1.0−|Gs|^2)))|1−Gs|
DS=(1−Gp·S11)·(1−Gs·S22)−(S12·S21·Gp·Gs)
Cascaded Circuits
For any given S-parameter matrix s, we can obtain the corresponding scattering transfer parameter (T-parameter) matrix (as taught in “S-Parameter Design”, Application Note 154, by Agilent Technologies):
using the formulae:
T11=−det(S)/S12
T12=SII/S12
T21=−S22/S12
T22=1/S12
Consider a pair of 2-port circuits described by their S-parameter matrices S1 and S2. By definition, the T-parameter matrix T of the cascaded circuit is given by:
T=T1·T2
where T1, T2 the T-parameter matrices corresponding to S1 and S2 respectively.
The S-parameter matrix of the cascaded circuit can then be calculated by the reverse transformation:
SII=T12/T22
S12=det(T)/T22
S21=1/T22
S22=−T21/T22
Tuning Capacitors
Consider a 2-port network consisting of a shunt capacitor C. The S-parameter matrix of this circuit would be:
where y(C,f) is the normalized admittance of C and f the frequency, whereby:
y(C,f)=j·2·π·f·C
The corresponding T-parameter matrix is given by:
Now, consider a transformer described by a set of S-parameter matrices S(i) and a pair of tuning capacitors c1, connected in parallel to port 1 (primary), and c2, connected in parallel to port 2 (secondary).
Then, the T-parameter matrix of the whole (tuned) transformer circuit, for any sampled frequency f(i), is given by:
Ttun(cl,c2,i)=Tc(cl,f(i))·T(i)·Tc(c2,f(i))
The S-parameter matrix stun (c1, c2, i) of the tuned transformer can easily be obtained by a simple T-to-S transformation.
Transformer Inductance and Quality Factor Calculation
Given a transformer described by:
Transformer Optimization to Arbitrary Loads
We next disclose a method according to the invention, for optimizing an instance of a transformer design coupled to arbitrary loads.
Given a transformer described by:
The above equations define a bounded constrained optimization problem, which is solvable by known methods, such as Powell's COBYLA (Constrained Optimization BY Linear Approximations) algorithm for derivative-free optimization with nonlinear inequality constraints, which is described in the paper by M. J. D. Powell, “A direct search optimization method that models the objective and constraint functions by linear interpolation”, Advances in Optimization and Numerical Analysis, 1994, p. 51-67, incorporated herein by reference.
The COBYLA method constructs linear approximations to the objective and the constraint functions by interpolating between the vertices of a n-dimensional simplex, where n the number of variables (in this case, the variables being c1 and c2, we have n=2).
We then define, according to the invention, a callback function callback(c1, c2) which performs the following operations:
The COBYLA optimization engine calls callback(cl, c2) repeatedly for various values of c1, c2, attempting to
The invention can be adapted by those skilled in the art, to incorporate alternative constrained optimization methods, equivalent to COBYLA and/or alternative constraints to meet the desired objectives.
Although the invention is presented here with embodiments involving transformers with primary and secondary coils, and with optional taps at either primary and/or secondary coils connected at fixed potential, those skilled in the art will realize that it can be extended and applied to the synthesis and optimization of ‘multi-tap’ transformers 400 as well (illustrated in
The above constrained optimization problem may be combined with a simulation based expert system method to accomplish physical design synthesis and optimization of transformers, including the design of tuning capacitors, in an advantageous manner. Those skilled in the art will appreciate upon reflecting on the teaching included herein that an EDA method thus implemented will automate the process of implementing integrated transformers and will moreover provide a valuable tool to circuit designers for optimizing their circuits employing transformers.
Referring to
To further enhance the usefulness of the invention, we propose a user interface whereby circuit designers may enter their desired parameters and accomplish automated transformer synthesis and optimization. Such parameters may include electrical specifications for the transformation, geometrical constraints for the transformer's physical layout, and other parameters, objectives and constraints that are relevant to the circuit being designed.
The present invention may be conveniently implemented using one or more conventional general purpose or specialized digital computer, computing device, machine, or microprocessor, including one or more processors, memory and/or computer readable storage media programmed according to the teachings of the present disclosure. Appropriate software coding can readily be prepared by skilled programmers based on the teachings of the present disclosure, as will be apparent to those skilled in the software art.
In some embodiments, the present invention includes a computer program product which is a storage medium or computer readable medium (media) having instructions stored thereon/in which can be used to program a computer to perform any of the processes of the present invention. The storage medium can include, but is not limited to, any type of disk including floppy disks, optical discs, DVD, CD-ROMs, microdrive, and magneto-optical disks, ROMs, RAMs, EPROMs, EEPROMs, DRAMs, VRAMs, flash memory devices, magnetic or optical cards, nanosystems (including molecular memory ICs), or any type of media or device suitable for storing instructions and/or data.
The previous description of the preferred embodiments is provided to enable any person skilled in the art to make or use the embodiments of the present invention. While the invention has been particularly shown and described with reference to preferred embodiments thereof, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention.
This application is a continuation application of U.S. Pat. No. 9,032,355 titled “Integrated Transformer Synthesis and Optimization, issued May 12, 2015, the entire disclosure of which is incorporated herein by reference in its entirety.
Number | Name | Date | Kind |
---|---|---|---|
5838723 | Mack et al. | Nov 1998 | A |
5867397 | Koza | Feb 1999 | A |
6026286 | Long | Feb 2000 | A |
6629044 | Papallo, Jr. | Sep 2003 | B1 |
6816012 | Aoki et al. | Nov 2004 | B2 |
6832182 | Wilson, Jr. | Dec 2004 | B1 |
6956944 | Koren | Oct 2005 | B1 |
7020452 | Kobayashi | Mar 2006 | B1 |
7243321 | Wang et al. | Jul 2007 | B2 |
7425869 | Aoki et al. | Sep 2008 | B2 |
7496871 | Suaya et al. | Feb 2009 | B2 |
7539965 | Bromberger et al. | May 2009 | B2 |
8032338 | Akselrod | Oct 2011 | B2 |
8108163 | Kashiwakura | Jan 2012 | B2 |
8181140 | Kireev et al. | May 2012 | B2 |
8250506 | Bantas et al. | Aug 2012 | B2 |
8667452 | Verghese et al. | Mar 2014 | B2 |
8769447 | Bantas et al. | Jul 2014 | B2 |
9171118 | Bantas | Oct 2015 | B2 |
20030182639 | Lehner et al. | Sep 2003 | A1 |
20050240381 | Seiler | Oct 2005 | A1 |
20060251221 | Rosenberg | Nov 2006 | A1 |
20090064060 | Yoshitomi | Mar 2009 | A1 |
20090193370 | Bantas et al. | Jul 2009 | A1 |
20090235208 | Nakayama et al. | Sep 2009 | A1 |
20100023896 | Rastogi et al. | Jan 2010 | A1 |
20100037189 | Bickel | Feb 2010 | A1 |
20100088657 | Bantas et al. | Apr 2010 | A1 |
20100175035 | He et al. | Jul 2010 | A1 |
20100259110 | Kurs et al. | Oct 2010 | A1 |
20110025435 | Gorbachov | Feb 2011 | A1 |
20110113401 | Kireev et al. | May 2011 | A1 |
20120022713 | Deaver, Sr. et al. | Jan 2012 | A1 |
20120112531 | Kesler et al. | May 2012 | A1 |
20120239117 | Kesler et al. | Sep 2012 | A1 |
20120242159 | Lou et al. | Sep 2012 | A1 |
20120248981 | Karalis et al. | Oct 2012 | A1 |
20120280765 | Kurs | Nov 2012 | A1 |
20130159956 | Verghese et al. | Jun 2013 | A1 |
20140181782 | Verghese et al. | Jun 2014 | A1 |
20140250417 | Bantas et al. | Sep 2014 | A1 |
Entry |
---|
Yang et. all; “Technique Improving and Its Simulation Analysis on Load Test of Transformer”; Year: 2012 ; 2012 Asia-Pacific Power and Energy Engineering Conference; pp. 1-4. |
Chiang et al., “Optimal capacitor placement, replacement and control in large-scale unbalanced distribution systems: modeling and a new formulation”, Power Systems, IEEE Transactions on: vol. 10, Issue 1: Publication Year: 1995, pp. 356-362. |
Number | Date | Country | |
---|---|---|---|
20150261900 A1 | Sep 2015 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 14045305 | Oct 2013 | US |
Child | 14710320 | US |