Field Solver Examples
The following example shows you how to use the Star-Hspice field solver. All the examples shown in this section are run with the
HIGH
accuracy mode and
GRIDFACTOR
= 1.
Example 1: Cylindrical Conductor Above a Ground Plane
In the first example, consider a copper cylindrical conductor above an ideal (lossless) ground plane. Cylindrical Conductor Above a Ground Plane shows the geometry. Input File Listing for Example 1 lists the corresponding Star-Hpsice netlist.
In this case, you can derive the exact analytical formulas for all of the transmission line parameters:
Table 18-4: Input File Listing for Example 1
Header, options and sources
|
*Example 1: cylindrical conductor
.OPTION PROBE POST
VIMPULSE in1 gnd PULSE 4.82v 0v 5n 0.5n 0.5n 25n
|
W Element
|
W1 in1 gnd out1 gnd FSmodel=cir_trans N=1 l=0.5
|
Materials
|
.MATERIAL diel_1 DIELECTRIC ER=4,
+ LOSSTANGENT=1.2e-3
.MATERIAL copper METAL
+ CONDUCTIVITY=57.6meg
|
Shapes
|
.SHAPE circle_1 CIRCLE RADIUS=0.5mm
|
Defines a half-space
|
.LAYERSTACK halfSpace BACKGROUND=diel_1,
+ LAYER=(PEC,1mm)
|
Option settings
|
.FSOPTIONS opt1 PRINTDATA=YES,
+ COMPUTERS=yes, COMPUTEGD=yes
|
Model definition
|
.MODEL cir_trans W MODELTYPE=FieldSolver
+ LAYERSTACK=halfSpace, FSOPTIONS=opt1,
+ RLGCFILE=ex1.rlgc
+ CONDUCTOR=(SHAPE=circle_1,
+ ORIGIN=(0,4mm), MATERIAL=copper)
|
Analysis, outputs and end
|
.TRAN 0.5n 100n
.PROBE v(out1)
.END
|
Compare the computed results with the analytical solutions in Comparison Result for Example 1. The resistance and conductance are computed at the frequency of 200 MHz, and the DC resistance (R0) and conductance (G0) are not included in the computed values.
Table 18-5: Comparison Result for Example 1
Value
|
Exact
|
Computed
|
C (pF/m)
|
89.81
|
89.66
|
L (nH/m)
|
494.9
|
495.7
|
G (mS/m)
|
0.1354
|
0.1352
|
R (
/m)
|
1.194
|
1.178
|
Example 2: Stratified Dielectric Media
Three Traces Immersed in Stratified Dielectric Media shows an example of three traces immersed in stratified dielectric media. The input file listing is shown in Input File Listing for Example 2.
Comparison Result for Example 2 compares the computed capacitance matrix with results from two other numerical methods.
Table 18-6: Input File Listing for Example 2
Header, options and sources
|
*Example 2, three traces in dielectric
.OPTION PROBE POST
+ VIMPULSE in1 gnd PULSE 4.82v 0v 5n 0.5n 0.5n
+ 25n
|
W Element
|
W1 in1 in2 in3 gnd out1 out2 out3 gnd
+ FSmodel=cond3_sys N=3 l=0.5
|
Materials
|
.MATERIAL diel_1 DIELECTRIC ER=4.3
.MATERIAL diel_2 DIELECTRIC ER=3.2
|
Shapes
|
.SHAPE rect_1 RECTANGLE WIDTH=0.35mm,
+ HEIGHT=0.07mm
|
Uses the default AIR background
|
.LAYERSTACK stack_1
+ LAYER=(PEC,1um),LAYER=(diel_1,0.2mm),
+ LAYER=(diel_2,0.1mm)
|
Option settings
|
.FSOPTIONS opt1 PRINTDATA=YES
|
Three conductors share the same shape
|
.MODEL cond3_sys W MODELTYPE=FieldSolver,
+ LAYERSTACK=stack1, FSOPTIONS=opt1,
+ RLGCFILE=ex2.rlgc
+ CONDUCTOR=(SHAPE=rect_1,ORIGIN=
+ (0,0.201mm)),
+ CONDUCTOR=(SHAPE=rect_1,
+ ORIGIN=(0.5mm,0.301mm)),
+ CONDUCTOR=(SHAPE=rect_1,ORIGIN=
+ (1mm,0.301mm))
|
Analysis, outputs and end
|
.TRAN 0.5n 100n
.PROBE v(out1)
.END
|
Table 18-7: Comparison Result for Example 2
Computed
|
(pF/m)
|
Raphael
(Finite-Difference Solver)
|
(pF/m)
|
Reference
|
(pF/m)
|
|
|
|
|
Convergence of Accuracy Modes shows the results of convergence analysis performed based on the total capacitance of the first conductor with respect to the
GRIDFACTOR
parameter.
Example 3: Two Traces Between Two Ground Planes
The following example uses the coupled strip line case shown in Two Traces Between Two Ground Planes. Input File Listing for Example 3 lists the complete input netlist. Comparison Result for Example 3 shows the comparison between the computed result and the Finite Element (FEM) solver result.
Table 18-8: Input File Listing for Example 3
Header, options and sources
|
*Example 3: two traces between gnd planes
.OPTION PROBE POST
+ IMPULSE in1 gnd PULSE 4.82v 0v 5n 0.5n 0.5n
+ 25n
|
W Element
|
W1 in1 in2 gnd out1 out2 gnd FSmodel=cond2_sys
+N=2 l=0.5
|
Materials
|
.MATERIAL diel_1 DIELECTRIC ER=10.0
.MATERIAL diel_2 DIELECTRIC ER=2.5
|
Shapes
|
.SHAPE rect RECTANGLE WIDTH=1mm,
+ HEIGHT=0.2mm,
|
Top and bottom ground planes
|
.LAYERSTACK stack_1,
+ LAYER=(PEC,1mm), LAYER=(diel_1,2mm),
+ LAYER=(diel_2,3mm), LAYER=(PEC,1mm)
|
Option settings
|
.FSOPTIONS opt1 PRINTDATA=YES
|
Two conductors share the same shape
|
.MODEL cond2_sys W MODELTYPE=FieldSolver,
+ LAYERSTACK=stack1, FSOPTIONS=opt1
+ RLGCFILE=ex3.rlgc
+ CONDUCTOR=(SHAPE=rect, ORIGIN=
+ (0,3mm)),
+ CONDUCTOR=(SHAPE=rect,
+ ORIGIN=(1.2mm,3mm))
|
Analysis, outputs and end
|
.TRAN 0.5n 100n
.PROBE v(out1)
.END
|
Table 18-9: Comparison Result for Example 3
Computed
|
(pF/m)
|
FEM Solver
|
(pF/m)
|
Example 4: Using Field Solver with Monte Carlo Analysis
The following example shows how to perform transient analysis using Monte Carlo analysis to model variations in the manufacturing of the microstrip. Input File Listing for Example 4 shows the Star-Hspice input listing with the W Element. Monte Carlo Analysis with a Field Solver and W Element shows the transient output waveforms.
Table 18-10: Input File Listing for Example 4
Header, options and sources
|
*PETL Example 4: example 2 with Monte-Carlo
.OPTION PROBE POST
+ VIMPULSE in1 gnd AC=1v PULSE 4.82v 0v 5ns
+ 0.5ns 0.5ns 25ns
|
Parameter definitions
|
.PARAM x1=Gauss(0,0.02,1)
+ x2=Gauss(0.5mm,0.02,1) x3=Gauss(1mm,0.02,1)
.PARAM dRef=1u dY1=Gauss(2mm,0.02,1)
+ dY2=Gauss(1mm,0.02,1)
|
W Element
|
W1 in1 in2 in3 0 out1 out2 out3 0
+ FSMODEL=cond3_sys N=3 l=0.5
|
Materials
|
.MATERIAL diel_1 DIELECTRIC ER=4.3
.MATERIAL diel_2 DIELECTRIC ER=3.2
|
Shapes
|
.SHAPE r1 RECTANGLE WIDTH=0.35mm,
+ HEIGHT=0.070mm
|
Uses the default AIR background
|
.LAYERSTACK stack_1
+ LAYER= (PEC,dRef),LAYER=(diel_1,dY1),
+ LAYER= (diel_2,dY2)
|
Three conductors share the same shape
|
.MODEL cond3_sys W MODELTYPE=FieldSolver,
+ LAYERSTACK=stack1,
+ CONDUCTOR=(SHAPE=r1,ORIGIN=
+ (x1,`dRef+dY1')),
+ CONDUCTOR=(SHAPE=r1,ORIGIN=
+ (x2,`dRef+dY1+dY2')),
+ CONDUCTOR=(SHAPE=r1,ORIGIN=
+ (x3,`dRef+dY1+dY2'))
|
Analysis, outputs and end
|
.PROBE TRAN v(in1) v(out1) v(in3)
.PROBE AC v(out1) v(out3)
.PROBE DC v(in1) v(out1) v(out3)
.AC LIN 200 0Hz 0.3GHz
.DC v1 0v 5v 0.01v
.TRAN 0.5ns 100ns SWEEP MONTE=3
.END
|
Star-Hspice Manual - Release 2001.2 - June 2001