Modelado de Diodos en SPICE: Guía Práctica para Simulación Precisa

SPICE es el estándar de la industria para simulación de circuitos, y los diodos son uno de sus componentes más modelados. El comportamiento de un diodo se captura en SPICE mediante dos declaraciones: el elemento del dispositivo y la definición del modelo. La declaración del elemento conecta el diodo al circuito, mientras la declaración del modelo especifica sus características eléctricas.

Definiendo un Diodo en SPICE

Todo elemento diodo debe iniciar con la letra d seguida por un identificador opcional (ej. d1, dtest, d101). Los siguientes dos números son los números de nodo de ánodo y cátodo. El token final es el nombre del modelo que se definirá después.

Ejemplo:

 d1 1 2 mod1 

La declaración de modelo correspondiente inicia con .model, seguido del nombre del modelo, la palabra clave d, y luego asignaciones de parámetros opcionales en forma Parametro=Valor. Los parámetros omitidos default a valores mostrados en las tablas abajo.

Ejemplo:

 .model mod1 d 

Ejemplo con parámetros:

 D2 1 2 Da1N4004 .model Da1N4004 D (IS=18.8n RS=0 BV=400 IBV=5.00u CJO=30 M=0.333 N=2) 

Parámetros Clave del Diodo

Obteniendo un Modelo Fiable

La forma más precisa de obtener un modelo SPICE es descargarlo directamente del sitio web del fabricante. Si no está disponible, puede construir un modelo desde el datasheet extrayendo los parámetros clave listados arriba. Como tercera opción, puede medir el dispositivo usted mismo y ajustar los parámetros SPICE a la curva I-V medida.

Ejemplo de Extracción de Parámetros – 1N4004

Para ilustrar el proceso, derivaremos un modelo para el rectificador común 1N4004. Comenzamos seleccionando un coeficiente de emisión N de 2,0, típico para rectificadores de silicio. Usando la ecuación del diodo

 I_D = I_S (e^{V_D/(N V_T)} - 1) 

y el punto del datasheet (V_D = 0,925 V a I_D = 1 A, V_T = 26 mV a 25 °C) resolvemos para I_S:

 I_S = 18,8 nA 

Otros parámetros se establecen a sus valores default a menos que el datasheet especifique lo contrario:

• RS = 0 Ω (inicialmente, luego ajustado para mejor ajuste)
• TT = 0 s (sin datos Q_RR disponibles)
• CJO = 30 pF (de gráfico V_R-C_J)
• M = 0,333 (unión linealmente graduada)
• VJ = 1 V, EG = 1,11 eV, XTI = 3,0
• BV = 400 V, IBV = 5 µA (especificaciones ruptura inversa)

Declaración de modelo resultante:

 .model Da1N4004 D (IS=18.8n RS=0 BV=400 IBV=5.00u CJO=30 +M=0.333 N=2.0 TT=0) 

Ajuste Fino del Modelo

Al comparar el modelo derivado con el modelo del fabricante y la curva I-V del datasheet, notamos que la corriente a 1,4 V excedía la especificación (12 A). Ajustando la resistencia en serie a 28,6 mΩ se alineó la curva a 1,4 V manteniendo el punto a 1 A. Mayor ajuste fino puede lograrse iterando N y RS hasta que la curva simulada coincida con la del datasheet en todo el rango.

Comparando Modelos

El siguiente circuito SPICE demuestra los tres modelos en paralelo para fácil comparación:

Extracto de netlist:

 *SPICE circuit <03468.eps> from XCircuit v3.20 D1 1 5 DI1N4004 V1 5 0 0 D2 1 3 Da1N4004 V2 3 0 0 D3 1 4 Default V3 4 0 0 V4 1 0 1 .DC V4 0 1400mV 0.2m .MODEL Da1N4004 D (IS=18.8n RS=28.6m BV=400 IBV=5.00u CJO=30 +M=0.333 N=2.0 TT=0) .MODEL DI1N4004 D (IS=76.9n RS=42.0m BV=400 IBV=5.00u CJO=39.8p +M=0.333 N=1.45 TT=4.32u) .MODEL Default D .END 

Diodos Zener, Tunel y Gunn

Para diodos Zener, existen dos estrategias de modelado:

• Establecer el parámetro BV al voltaje Zener (ej. .model D1N4469 D (BV=15 IBV=17m)). Este es el enfoque más simple.
• Usar un subcircuito con un diodo sujetador que fuerza el voltaje de ruptura (ver subcircuito abajo para modelo 1N4477A).

Ejemplo de subcircuito:

 .SUBCKT DI-1N4744A 1 2 * Terminales A K D1 1 2 DF DZ 3 1 DR VZ 2 3 13.7 .MODEL DF D (IS=27.5p RS=0.620 N=1.10 +CJO=78.3p VJ=1.00 M=0.330 TT=50.1n) .MODEL DR D (IS=5.49f RS=0.804 N=1.77) .END 

Para diodos tunel, un par de JFETs puede emular la región de resistencia diferencial negativa, como se muestra en referencia [KHM]. Similarmente, un diodo Gunn puede modelarse con dos JFETs para capturar su comportamiento de oscilación de microondas ([ISG]).

Conclusiones Clave

• Use el elemento diode para conectar el dispositivo; la declaración .model define su comportamiento eléctrico.
• Parámetros estáticos (N, IS, RS) dominan rendimiento a baja frecuencia; especificaciones de ruptura inversa (BV, IBV) importan para aplicaciones de alto voltaje.
• Comportamiento dinámico depende de TT y CJO; inclúyalos para simulaciones transitorias precisas.
• Modelos de fabricante son preferidos; si no disponibles, derive de datos de datasheet o mediciones empíricas.