Simple and efficient full-wave modeling of electromagnetic coupling in realistic RF multilayer PCB layouts

A simple method to model a printed circuit board (PCB) that takes advantage of the unique features found in PCBs is proposed. This method is capable of analyzing coupling between any nets in the entire multilayer PCB. Using the equivalence principle, the PCB is modeled as a cascade of parallel-plate waveguides with half-space regions residing above and below the PCB. Instead of formulating the problem in terms of electric currents in the horizontal metal layers, it is formulated using equivalent magnetic currents in the nonmetallic regions of layer interfaces. The equivalent magnetic currents at the dielectric interfaces are expressed in terms of the Rao-Wilton-Glisson (RWG) basis functions. The electric currents flowing on the vias inside dielectric layers are assumed constant in the vertical direction. These vertical electric currents radiate TEM modes in the parallel-plate environment. Integral equations based on simple parallel-plate and free-space Green's functions enforcing the boundary conditions are set up and solved using the method of moments (MoM). The equivalent magnetic currents in each layer interact with currents in the adjacent layers only, thereby resulting in a "chained-block-banded" matrix. Excitation is provided through ports defined at each pair of pads, or between a pad and nearby ground. These ports are located only on the top and the bottom layers of the PCB where the circuit components and integrated-circuit pins are mounted. Two different localized excitation schemes, one with a current loop injection and the other with a strip current excitation, are proposed. This formulation requires the computation of the MoM matrix only once per frequency for any number of ports. Further, the solution for only those unknown equivalent magnetic currents around the port regions is required to obtain the N-port impedance parameter characterization of the PCB. Consequently, a memory-efficient block matrix solution process can be used to solve problems of large size for a given memory. Simple and realistic examples are given to illustrate the applicability of this approach.

[1]  R. F. Milsom,et al.  Equivalent circuit model for PCB ground planes , 1992 .

[2]  David C. Chang,et al.  Electromagnetic modeling of passive circuit elements in MMIC , 1992 .

[3]  O. Fordham,et al.  Modeling planar arbitrarily shaped microstrip elements in multilayered media , 1997 .

[4]  Roger F. Harrington,et al.  A generalized network formulation for aperture problems , 1976 .

[5]  Ahmed A. Kishk,et al.  Green's functions analysis of planar circuits in a two-layer grounded medium , 1992 .

[6]  Linda P. B. Katehi,et al.  Application of system-level EM modeling to high-speed digital IC packages and PCBs , 1997 .

[7]  Yongshik Lee,et al.  Improved complex image method for a horizontal magnetic dipole in a parallel-plate waveguide , 1997 .

[8]  Roger F. Harrington,et al.  Field computation by moment methods , 1968 .

[9]  Jian-Ming Jin,et al.  An efficient algorithm for analyzing large-scale microstrip structures using adaptive integral method combined with discrete complex image method , 1998, ICMMT'98. 1998 International Conference on Microwave and Millimeter Wave Technology. Proceedings (Cat. No.98EX106).

[10]  T. Vaupel,et al.  Electrodynamic analysis of combined microstrip and coplanar/slotine structures with 3-D components based on a surface/volume integral-equation approach , 1999 .

[11]  Nicolaos G. Alexopoulos,et al.  An asymptotic extraction technique for evaluating Sommerfeld-type integrals , 1986 .

[12]  C. Balanis Advanced Engineering Electromagnetics , 1989 .

[13]  J. Mosig Arbitrarily shaped microstrip structures and their analysis with a mixed potential integral equation , 1988 .

[14]  Michael J. Lancaster,et al.  Aperture-coupled microstrip open-loop resonators and their applications to the design of novel microstrip bandpass filters , 1999 .

[15]  Albert E. Ruehli,et al.  Retarded models for PC board interconnects-or how the speed of light affects your SPICE circuit simulation , 1991, 1991 IEEE International Conference on Computer-Aided Design Digest of Technical Papers.

[16]  Ronghong Jin,et al.  Accelerant method for Green's function of horizontal dipole in parallel plate waveguide , 1998 .

[17]  Kapur,et al.  IES/sup 3/: a fast integral equation solver for efficient 3-dimensional extraction , 1997, 1997 Proceedings of IEEE International Conference on Computer Aided Design (ICCAD).

[18]  G. E. Howard,et al.  A closed-form spatial Green's function for the thick microstrip substrate , 1991 .

[19]  S. Kim,et al.  Fast capacitance extraction of general three-dimensional structures , 1991, [1991 Proceedings] IEEE International Conference on Computer Design: VLSI in Computers and Processors.

[20]  E. R. Pillai,et al.  Coax via-A technique to reduce crosstalk and enhance impedance match at vias in high-frequency multilayer packages verified by FDTD and MoM modeling , 1997 .

[21]  Y. Chow,et al.  Simulated image method for Green's function of multilayer media , 1999 .

[22]  Mattan Kamon,et al.  FastHenry: A Multipole-Accelerated 3-D Inductance Extraction Program , 1993, 30th ACM/IEEE Design Automation Conference.

[23]  D. Wilton,et al.  Potential integrals for uniform and linear source distributions on polygonal and polyhedral domains , 1984 .

[24]  H.-I. Liaw,et al.  Simulation and modeling of mode conversion at vias in multilayer interconnections , 1995, 1995 Proceedings. 45th Electronic Components and Technology Conference.

[25]  Jian-Ming Jin,et al.  Multilevel fast multipole algorithm for analysis of large-scale microstrip structures , 1999 .

[26]  R. Collin Field theory of guided waves , 1960 .

[27]  N. Kinayman,et al.  Efficient use of closed-form Green's functions for the analysis of planar geometries with vertical connections , 1997 .

[28]  A. Ruehli,et al.  Efficient Capacitance Calculations for Three-Dimensional Multiconductor Systems , 1973 .

[29]  W. Weeks Calculation of Coefficients of Capacitance of Multiconductor Transmission Lines in the Presence of a Dielectric Interface , 1970 .

[30]  Jonathan Richard Shewchuk,et al.  Triangle: Engineering a 2D Quality Mesh Generator and Delaunay Triangulator , 1996, WACG.

[31]  R. Jackson Notes on microwave circuits , 1986, IEEE Antennas and Propagation Society Newsletter.

[32]  A. E. Ruehii Inductance Calculations in a Complex Integrated Circuit Environment , 2002 .

[33]  M. I. Aksun A robust approach for the derivation of closed-form Green's functions , 1996 .

[34]  D. Wilton,et al.  Electromagnetic scattering by surfaces of arbitrary shape , 1980 .

[35]  K. Michalski,et al.  Multilayered media Green's functions in integral equation formulations , 1997 .

[36]  G. Y. Delisle,et al.  Discrete image theory for horizontal electric dipoles in a multilayered medium , 1988 .

[37]  Jen-Tsai Kuo,et al.  Fast full-wave analysis of planar microstrip circuit elements in stratified media , 1998 .

[38]  J. Mosig,et al.  Modeling the coplanar transmission line excitation of planar antennas in the method of moments , 1997 .

[39]  Albert E. Ruehli,et al.  Challenges and advances in electrical interconnect analysis , 1992, [1992] Proceedings 29th ACM/IEEE Design Automation Conference.