Brevet FR2860086A1 — publié le 20050325

DÉTERMINATION DES TENSIONS D'UN CIRCUIT ÉLECTRIQUE La présente invention concerne la simulation d'un circuit électrique et plus particulièrement des procédés de détermination des variations temporelles des tensions et des courants d'un circuit électrique en fonction de stimuli externes. La présente invention s'applique tout particulièrement aux circuits électriques comprenant des composants analogiques tels que des transistors, des résistances, des condensateurs, des diodes et des inductances. La figure 1 est un exemple de circuit simple compre- nant une diode d et une résistance variable r en série entre une borne d'alimentation supérieure vdd et une borne d'alimentation inférieure gnd. La valeur R de la résistance r est commandée par un signal de couunande C. Dans cet exemple, on souhaite déterminer la tension au noeud N entre l'anode de la diode d et la résistance r ainsi que le courant traversant la résistance r et la diode d. La figure 2 est un graphique représentant les varia- tions du courant dans chaque composant en fonction de la tension VN au noeud N. La courbe Cd représente les variations du courant id traversant la diode d en fonction de la tension VN. La droite Cr, représente les variations du courant ir traversant la 2860086 2 résistance r en fonction de la tension VN pour une tension VC donnée. Au vu de ce graphique, il est possible de déterminer le courant traversant la résistance r et la diode d ainsi que la tension VN au noeud N. En effet, d'après la loi des mailles, le courant id traversant la diode d est égal au courant ir traversant la résistance r. Le point de croisement PO de la droite Cr et de la courbe Cd correspond à un courant i0 et à une tension VNO pour lesquels le courant dans la diode d est égal au courant dans la résistance r. L'exemple de circuit décrit ci-dessus est extrêmement simple et les modèles mathématiques du fonctionnement de la diode d et de la résistance r sont également très simples. De façon générale, dans le domaine de la simulation des circuits électriques, les circuits analysés peuvent comprendre plusieurs centaines de composants et les modèles mathématiques du fonctionnement de chaque composant peuvent inclure plusieurs dizaines d'équations. Contrairement à l'exemple précédent où la détermination du courant i0 et de la tension VNO consiste à résoudre une équation simple, la détermination des tensions et des courants d'un circuit plus complexe nécessite de résoudre un grand nombre d'équations différentielles généralement non linéaires. En pratique, la détermination des tensions et des courants d'un circuit électrique est réalisée selon un procédé itératif décrit ci-après à partir de l'exemple simple illustré en figure 1. La figure 3 est un graphique semblable au graphique de la figure 2. Les droites Cr0 et Crl représentent les variations du courant ir traversant la résistance r en fonction de la tension VN pour deux valeurs différentes de la tension VC, respectivement VCO et Vol, VC1 étant supérieure à VCO. La courbe Cd représente, canne précédemment, les variations du courant id dans la diode d en fonction de la tension VN. A un instant t0, la tension VC est égale à VCO. On considère que la tension VN et les courants id et ir sont connus et qu'ils sont égaux respectivement à VNO, i0 et i0, la tension VNO et le courant i0 correspondant au point de croisement PO entre la droite Cr0 et la courbe Cd. A un instant tl suivant t0, la tension VC change et devient égale à Vol. La détermination de la tension VN et des courants id et ir est alors effectuée selon le procédé suivant. On ne connaît pas a priori le sens de variation de la tension VN, et on fait l'hypothèse que la tension VN à l'instant t1 est inchangée et égale à VNO. On calcule alors le courant id à l'instant t1 à partir du modèle de la diode d. Le courant id est égal au courant i0 étant donné que les tensions aux bornes de la diode n'ont pas changé. En faisant toujours l'hypothèse que la tension VN est inchangée à l'instant t1, on calcule le courant ir à l'instant t1 à partir du modèle de la résistance r pour une tension de commande VC égale à Vol. Le courant ir est égal au courant ir,a correspondant au point Pria de la droite Cri pour une tension VN égale à VNO. On vérifie ensuite la loi des mailles au noeud N. En d'autres termes, on regarde si le courant id est égal au courant ir. Dans cet exemple, les courants ip et irl,a sont manifes- tement différents, comme cela apparaît en figure 3. On calcule ensuite la dérivée du courant id en fonction de la tension VN pour une tension VN égale à VNO. On peut ainsi définir l'équation de la tangente Td,a à la courbe Cd au point P0. On calcule de même la dérivée du courant ir en fonction de la tension VN pour une tension égale à VNO. Puis on définit l'équation de la tangente Tr,a à la courbe Crl au point Pr,a. La tangente Tr,a est en fait confondue avec la droite Cri. Les tangentes Td,a et Tr,a précédemment calculées indiquent le sens de variation des courants ir et id quand la tension VN varie autour de VNO. La tangente Td,a est une repré- sentation simplifiée de la courbe Cd qui permet d'effectuer des estimations sur la valeur du courant en fonction du sens de variation de la tension VN. On identifie alors le point de croisement PC,a des deux tangentes Td,a et Tr,a. Le point de croisement PC,a permet d'identifier une tension VN égale à VNC,a pour laquelle les courants ir et id sont susceptibles d'être identiques. On fait alors une nouvelle hypothèse sur la valeur de la tension VN à l'instant tl en la prenant égale à la tension VNC,a. On calcule les courants id et ir correspondant à la tension VNC,aÉ Le courant ir est égal au courant iC,a correspondant au point PC,a. Le courant id est égal au courant id,b correspondant au point Pd,b de la courbe Cd pour une tension VN égale à VNC,aÉ On vérifie ensuite de nouveau la loi des mailles au noeud N. Les courants iC,a et id,b ne sont pas strictement égaux mais l'écart est faible. Si la précision souhaitée est peu élevée, on considère que la tension VN à l'instant tl est égale à la tension VNC,a et que le courant traversant la diode d et la résistance r est égal à ic,a. Dans le cas inverse, on affine la recherche de la tension VN et du courant traversant la diode d et la résistance r en calculant les équations des tangentes Td,b et Tr,b respectivement aux points Pd,b et PC,a pour faire une nouvelle hypothèse quant à la tension VN. On calcule ensuite les courants pour cette nouvelle tension VN égale à VNC,b et on vérifie de nouveau la loi des mailles. On continue la recherche jusqu'à converger vers des courants suffisamment proches. De façon générale, une simulation d'un circuit tel que celui de la figure 1 a pour but de déterminer les tensions et les courants d'un circuit recevant des stimuli. On détermine tout d'abord l'état initial du circuit à partir des valeurs initiales des stimuli à un instant initial t0 de simulation. Pour ce faire, on donne des valeurs de tension par défaut à chaque noeud du circuit puis on procède à une première série d'itérations de "convergence" selon un procédé identique à celui décrit précédemment, jusqu'à obtenir des tensions initiales et des courants entrant ou sortant de chacun des composants du circuit qui satisfassent à la loi des mailles pour ces tensions initiales. On détermine ensuite l'état du circuit à un instant tl suivant l'instant t0, le temps écoulé entre l'instant t0 et l'instant t1 correspondant au "pas" de simulation. L'état du circuit à l'instant t1 est déterminé par une deuxième série d'itérations telle que celle décrite précédemment en relation à la figure 3, en faisant l'hypothèse que les tensions sont inchangées. On détermine ensuite l'état du circuit aux instants suivants t2 à tn, où n est un entier positif. La détermination des tensions à un instant donné tj, j étant compris entre 2 et n, consiste à définir des tensions "de départ" préalablement à la mise en oeuvre d'une série d'itérations de convergence. La définition des tensions "de départ" tient compte des tensions calculées aux instants précédents. On peut par exemple faire une extrapolation en tenant compte des tensions calculées pour les deux derniers instants. L'utilisation de procédés de simulation classiques tels que celui décrit précédemment amène à des durées de simulation longues. Ces durées deviennent prohibitives quand on veut simuler des circuits de plus en plus complexes en utilisant des modèles de composants de plus en plus précis et de plus en plus complexes. On est donc amené à réduire le nombre de pas de simulation. Un objet de la présente invention est de prévoir un procédé plus rapide de détermination des variations temporelles des tensions et des courants d'un circuit de façon à permettre des simulations complètes en des durées raisonnables. Pour atteindre cet objet, la présente invention prévoit un procédé de simulation d'un circuit électrique comprenant des composants élémentaires dont chacune des bornes est reliée à un noeud du circuit, certains noeuds recevant des stimuli externes consistant en une suite de valeurs de tensions appliquées à des instants successifs de simulation, dans lequel la détermination à un instant donné des tensions aux noeuds du circuit consiste à effectuer une ou plusieurs itérations de calcul, chaque itération consistant à définir une tension probable pour chaque noeud, à calculer les courants entrant ou sortant dans chacun des composants en utilisant pour chaque composant un modèle de calcul des courants en fonction des tensions aux bornes du composant, puis à vérifier la loi des mailles pour chaque noeud, la suite d'itérations continuant jusqu'à ce que la loi des mailles soit vérifiée, et dans lequel pour le premier instant de simulation suivant l'instant initial de simulation, lors d'une itération de calcul, un courant entrant ou sortant d'un composant donné est calculé à partir de l'un de trois modèles, simplifié, mixte ou complet, le modèle complet permettant de connaître très précisément la valeur du courant en fonction des tensions aux bornes du composant, le modèle simplifié correspondant à une linéarisation tangentielle du modèle complet autour de la valeur du courant correspondant aux tensions aux bornes du composant déterminées pour l'instant initial, le modèle mixte constituant un raccordement des modèles simplifié et complet en assurant la continuité des courants, le choix d'un modèle étant fonction de l'écart entre les tensions aux bornes du composant déterminées pour l'instant initial de la simulation et celles jugées probables pour le premier instant de simulation, de sorte que pour chaque composant: si l'écart est inférieur à un premier seuil défini pour ce composant, le modèle simplifié est utilisé ; si l'écart est supérieur à un second seuil défini pour ce composant, ce second seuil étant supérieur au premier seuil, le modèle complet est utilisé ; et si l'écart est compris entre les premier et second seuils pour ce composant, le modèle mixte est utilisé. Selon une variante de mise en oeuvre du procédé susmentionné, pour un instant de simulation donné, la première itération de calcul comprend, préalablement au calcul des courants, une étape supplémentaire consistant à choisir pour chaque composant un instant de simulation antérieur à l'instant donné pour lequel les courants entrant ou sortant du composant ont été calculés à partir du modèle complet du composant, le calcul d'un courant entrant ou sortant d'un composant donné, lors de la première itération et des suivantes, étant effectué à partir de l'un de trois modèles, simplifié, mixte ou complet, le modèle simplifié étant défini autour de la valeur du courant correspondant aux tensions aux bornes du composant déterminées pour ledit instant de simulation antérieur, le choix d'un modèle étant fonction de l'écartement entre les tensions aux bornes du composant déterminées pour l'instant antérieur choisi et celles jugées probables pour l'instant de simulation donné comme pour le premier instant. Selon une variante de mise en oeuvre du procédé susmentionné, pour un composant à n bornes, ledit écart est égal à la distance entre des premier et second points d'un espace à n dimensions, le premier point ayant pour coordonnées les tensions aux bornes du composant jugées probables pour l'instant donné, et le second point ayant pour coordonnées les tensions aux bornes du composant déterminées pour l'instant antérieur choisi. Selon une variante de mise en oeuvre du procédé susmentionné, l'ensemble des valeurs possibles d'un courant entrant ou sortant d'un composant donné correspondant au modèle complet du composant constitue une première hypersurface représentée dans un espace à n+l dimensions comprenant un axe des courants et ledit espace de tensions, les ensembles de points, dudit espace de tensions, éloignés du second point d'une valeur égale à la première tension de seuil ou égale à la seconde tension de seuil constituent des première et seconde hypersphères, l'ensemble des valeurs de courant représentatives du modèle simplifié du composant constitue une portion d'un hyperplan tangent à la première hypersurface pour un courant correspondant au second point, chacun des courants de la portion d'hyperplan correspondant à un point situé à l'intérieur de la première hypersphère, l'ensemble des valeurs de courant représentatives du modèle mixte du composant constitue une seconde hypersurface de raccordement entre la portion d'hyperplan et un sous-ensemble de valeurs de courant de la première hypersurface correspondant aux points de la seconde hypersphère. Selon un mode de mise en oeuvre du procédé susmen-35 tionné, le modèle mixte d'un composant est défini partiellement de la façon suivante lorsque ledit écart est compris entre les premier et second seuils: définir des troisième et quatrième points correspondant à l'intersection entre un axe ayant pour origine le second point et passant par le premier point et les première et seconde hypersphères; calculer des troisième et quatrième courants correspondant aux troisième et quatrième points à partir respectivement des modèles simplifié et complet du composant, et déterminer un troisième ensemble de valeurs de courant pour chacun des points dudit axe compris entre les troisième et quatrième points, le troisième ensemble de valeurs de courant reliant les troisième et quatrième valeurs de courant, la valeur du courant calculé pour le premier point étant la valeur de courant du troisième ensemble correspondant au premier point. Selon une variante de mise en oeuvre du procédé susmentionné, ledit instant de simulation antérieur est l'instant le plus proche de l'instant donné pour lequel les courants entrant ou sortant du composant ont été calculés à partir du modèle complet du composant. Selon une variante de mise en oeuvre du procédé susmentionné, ledit instant de simulation antérieur est tel que les tensions aux bornes du composant déterminées pour cet instant antérieur soient les plus proches des tensions aux bornes du composant jugées probables pour l'instant donné et tel que pour les tensions aux bornes du composant déterminées pour cet instant, les courants entrant ou sortant du composant correspondant aient été calculés à partir du modèle complet du composant. Cet objet, ces caractéristiques et avantages, ainsi que d'autres de la présente invention seront exposés en détail dans la description suivante de modes de réalisation particuliers faite à titre non-limitatif en relation avec les figures jointes parmi lesquelles: la figure 1 est un schéma, précédemment décrit, d'un 35 exemple de circuit électrique; la figure 2 est un graphique représentant le courant de chaque composant du circuit de la figure 1 en fonction des tensions à ses bornes; la figure 3 est un graphique semblable à celui de la 5 figure 2 sur lequel sont représentés divers éléments du procédé mis en oeuvre lors de la simulation; la figure 4 est un schéma d'un autre exemple de circuit; la figure 5 est une représentation en trois dimensions des valeurs d'un courant en fonction de tensions aux deux bornes d'un composant correspondant au modèle "complet" du composant; la figure 6 est une représentation en deux dimensions de l'ensemble des couples de tensions sur les bornes d'un composant; la figure 7 est une représentation en trois dimensions des valeurs d'un courant en fonction de tensions aux deux bornes d'un composant représentant les trois modèles utilisés dans le procédé de la présente invention; la figure 8 est une représentation en deux dimensions comprenant les éléments de la figure 6 plus d'autres points; et la figure 9 est une vue en coupe d'une partie de la représentation de la figure 7. La présente invention s'appuie sur la constatation faite par l'inventeur qu'une grande partie du temps de simu- lation est utilisée pour le calcul des courants à partir des modèles mathématiques souvent très complexes de chaque composant. En effet, le calcul des courants entrant ou sortant par un plot de connexion d'un composant en fonction des tensions estimées pour les noeuds auxquels sont reliés ses plots de connexion consiste à effectuer une série d'opérations souvent très longue. A titre d'exemple, le calcul d'un courant d'un transistor en fonction des tensions sur ses électrodes nécessite plusieurs centaines d'opérations de calcul. L'inventeur a en outre noté que, lors du calcul des 35 courants d'un composant, la suite d'opérations effectuée comprend toujours le calcul des dérivées de chaque courant en fonction des diverses tensions aux noeuds auxquels est relié le composant. Cette constatation permet d'envisager un procédé de simulation plus rapide comme cela apparaîtra ci- après. La figure 4 est un exemple de circuit électrique comprenant trois composants à deux bornes Z1, Z2 et Z3. Le composant Zl est placé entre un noeud N1 et un noeud N2, le composant Z2 est placé entre le noeud N2 et un noeud N3 et le composant Z3 est placé entre le noeud N2 et la masse. Une source de tension X1 est placée entre le noeud N3 et la masse. Une source de tension X2 est placée entre le noeud Nl et le noeud N3. Lors de la simulation du circuit de la figure 4, on cherche à déterminer la tension VN2 au noeud N2 en fonction des tensions VN1, VN3 appliquées sur les noeuds Nl et N3 par l'intermédiaire des sources de tension Xl et X2. De plus, on souhaite connaître les courants entrant ou sortant des composants Z1 à Z3. Par convention dans la suite de la description, on cherche à connaître la valeur algébrique des courants entrant dans un composant par un de ses plots de connexion. Les courants entrant dans les composants Z1, Z2 et Z3 en provenance du noeud N2 sont respectivement appelés il, i2 et i3. La figure 5 est une représentation en trois dimensions d'un des courants, dans cet exemple il, entrant dans un des composants, Z1, en fonction des tensions, VN1 et VN2, sur les bornes du composant. Le courant il varie selon un axe vertical. La tension VN1 varie selon un axe horizontal. La tension VN2 varie selon un axe perpendiculaire au plan de la feuille et représenté schématiquement avec un angle à 45 par rapport à chacun des axes, comme cela est courant dans une représentation à trois dimensions. A titre d'exemple, les sources de tension Xl et X2 sont telles que les tensions VN1 et VN2 peuvent prendre une valeur comprise entre 0 et 3 volts. L'ensemble des couples de tensions (VN1,VN2) est représenté par une portion de plan 30 ayant une forme de parallélogramme. L'ensemble des valeurs de courant il est représenté par une surface incluant une demi-sphère 31 posée sur la portion de plan 30. L'intersection de la demi-sphère 31 et de la portion de plan 30 est un cercle 32 ayant pour centre le couple de tensions (VN1,VN2) égal à 1,5/1, 5 volts, et pour rayon 1,5 volts de sorte que le cercle touche les points milieux de chacun des côtés de la portion de plan 30. Le courant il est maximal pour le couple de tensions (VN1,VN2) égal à 1,5/1,5 volts. Le courant il est nul pour l'ensemble des couples de tensions (VN1,VN2) placés à l'extérieur du cercle 32. Le modèle de composant illustré en figure 5 est donné à titre d'exemple. De façon générale, à chaque type de composant d'un circuit électrique est associé un modèle "complet" permettant de connaître de façon précise la valeur de chacun des courants entrant dans le composant en fonction des tensions à ses bornes. Dans le cas où le composant comprend k bornes, la représentation d'un courant en fonction des tensions sur chacune des k bornes est effectuée dans un espace à k+l dimensions. Le procédé de simulation de la présente invention consiste, comme précédemment décrit en relation aux figures 1 à 3, à déterminer la tension VN2 pour chaque instant ti en effectuant une série d'itérations de convergence. Lors de chaque itération, le simulateur détermine une tension VN2 probable puis il "interroge" les modèles de chacun des composants afin de connaître la valeur des courants entrant dans chacun des composants. La série d'itérations continue tant que la loi des mailles n'est pas vérifiée au noeud N2. La présente invention vise à perfectionner le procédé d'interrogation des modèles lors de chaque itération. Le procédé de la présente invention est décrit ci-après pour le calcul du courant il entrant dans le composant Z1. Le procédé de la présente invention vise à accélérer le calcul du courant il en utilisant un modèle "simplifié" du composant Z1 lorsqu'un couple de tensions (VN1,VN2) précédeuunent déterminé lors d'un instant antérieur de la simulation est proche d'un couple de tensions (VN1,VN2) estimé conne étant un couple probable pour l'instant de simulation tj. Bien que la tension VN1 à l'instant tj soit imposée et donc connue, le calcul des courants est effectué à partir des valeurs de tension de l'ensemble des noeuds, N1 et N2, auxquels est relié Z1. La figure 6 représente l'ensemble des couples VN1,VN2 appartenant à la portion de plan 30 représentée en figure 5. La tension VN1 varie selon un axe horizontal et la tension VN2 varie selon un axe vertical. La portion de plan 30 a en figure 6 la forme d'un carré. Le cercle 32 est placé à l'intérieur de la portion de plan 30. Dans le procédé de la présente invention, on estime lors de la première itération d'une série d'itérations permet-tant de déterminer l'état du circuit à l'instant tj, si le couple de tensions (VN1,VN2)j estimé probable pour l'instant de simulation tj est éloigné ou non d'un couple de tensions (VN1,VN2)j_x déterminé pour un instant de simulation antérieur tj_x, où x est compris entre 1 et j. Le choix de l'instant antérieur tj_x sera décrit ci-après. Le couple (VN1,VN2)j correspond à un point Pj de la portion de plan 30. Le couple de tensions (VN1lVN2)j-x correspond à un point Pj_x de coordonnées 1/1 dans cet exemple. L'estimation de l'éloignement entre les points Pj_x et Pj peut consister à déterminer la distance entre ces deux points. La distance entre les deux points est représentative d'un écart de tensions égal à : J(VN1, p - VN1, r)2 + (VN2, p - VN2, r)2 où VN1,p est la tension VN1 estimée pour l'instant tj, VN2,p est la tension VN2 estimée pour l'instant tj, VNl,r est la tension VN1 retenue pour l'instant tj_x, et VN2,r est la tension VN2 retenue pour l'instant tj_x. Dans le cas où la distance calculée est inférieure à une tension de seuil Vsl, la présente invention prévoit que le calcul du courant sera effectué à partir d'un modèle "simplifié" du composant défini ci-après. Dans cet exemple, la tension de seuil Vsl est choisie égale à 0,25 volt. L'ensemble des points du plan 30 dont la distance au point Pj_x est inférieur à la tension de seuil Vsl est une surface délimitée par un cercle 35 ayant pour centre le point Pi _x et pour rayon Vsl. Dans le cas où la distance calculée est supérieure à une seconde tension de seuil Vs2, la tension de seuil Vs2 étant supérieure à la tension de seuil Vsl, la présente invention prévoit que le calcul du courant sera effectué à partir du modèle complet du composant tel que représenté en figure 5. La tension de seuil Vs2 est dans cet exemple égale à 0,5 volt. L'ensemble des points de la portion du plan 30 dont la distance au point Pj_x est supérieure à la tension de seuil Vs2 sont placés à l'extérieur d'un cercle 36 de centre Pj _x et de rayon Vs2. Dans le cas où la distance calculée est supérieure à la tension de seuil Vsl et inférieure à la tension de seuil Vs2, la présente invention prévoit que le calcul du courant sera effectué à partir d'un modèle "mixte" défini ci-après. La figure 7 est une représentation du courant il en fonction des tensions VN1 et VN2 comprenant l'ensemble des éléments représentés en figure 5 ainsi que la représentation des modèles "simplifié" et "mixte" cités précédemment. La définition du modèle simplifié est la suivante. L'instant antérieur tjx a été choisi de sorte que le modèle complet du composant a été utilisé pour calculer le courant ij_x correspondant au couple de tensions (VN1, VN2)j-x déterminé pour l'instant antérieur tj_x. Le courant ij_x est indiqué en figure 7 sur la demi-sphère 31 au-dessus du point Pj_x. Lors du calcul du courant ij_x à partir du modèle complet, les dérivées du courant il en fonction de VN1 et VN2 sont aussi calculées. Les valeurs des dérivées du courant il permettent de définir un plan tangent à la demi-sphère 31 au point de la demi-sphère 31 correspondant au courant ij_x. Ce plan tangentiel indique le sens de variation du courant il en fonction des variations des tensions VN1 et VN2 autour du point Pj_x. La portion du plan tangentiel placée au-dessus de la surface délimitée par le cercle 35 est une surface plane elliptique 40 en contact avec la demi-sphère 31 uniquement au-dessus du point Pj_x. La surface 2860086 14 plane elliptique 40 est une représentation de l'ensemble des valeurs de courant il du modèle simplifié défini pour le courant ij_x. Lorsque le modèle simplifié est utilisé, le calcul du courant il est effectué à partir de l'équation du plan tan- gentiel. Ce calcul nécessite une unique opération de calcul ce qui permet d'accélérer la simulation. Le modèle mixte est quant à lui un modèle de "raccordement" entre le modèle "simplifié" et le modèle "complet". Les valeurs de courant il de la demi-sphère 31 correspondant à l'ensemble des couples de tensions (VN1, VN2) du cercle 36 forment une courbe fermée 41 sur la demi-sphère 31. Une surface 42 représentative de l'ensemble des valeurs de courant il du modèle mixte a une forme sensiblement conique "tronquée" ayant pour base la courbe fermée 41 et pour surface supérieure tronquée la surface plane elliptique 40. Les valeurs de courant de la surface 42 correspondent à des couples de tensions (VN1,VN2) placés dans une région 45 de la portion de plan 30 comprise entre les cercles 35 et 36. L'utilisation d'un modèle mixte pour le calcul du courant il permet d'assurer que la suite d'itérations effectuée pour déterminer le couple de tensions (VN1,VN2) converge. Le modèle mixte permet en effet d'assurer la continuité entre les valeurs de courant de la surface plane elliptique 40 et les valeurs de courant de la demi-sphère 31 située à l'extérieur de la courbe fermée 41. Lorsque le modèle mixte est utilisé pour le calcul du courant il, il est nécessaire de définir un jeu d'équations qui permette de calculer le courant il en fonction d'un couple de tensions (VN1,VN2). Afin de simplifier la détermination de ce jeu d'équations, la présente inventionprévoit de déterminer une unique équation permettant de calculer le courant il pour un sous-ensemble choisi de couples (VN1,VN2) appartenant à la région 45. La figure 8 est une représentation identique à celle de la figure 6. Les coordonnées du point Pi sont dans cet exemple 1,2/0,6. Lorsque le simulateur détecte que le point Pi appartient à la région 45 d'utilisation d'un modèle mixte, il détermine les coordonnées de deux points Ps et PC placés respectivement sur les cercles 35 et 36 et placés sur un axe 50 ayant pour origine le point Pi_x et passant par le point Pi. La figure 9 représente les valeurs de courant il correspondant à l'un des trois modèles pour un couple de tensions (VN1,VN2) appartenant à l'axe 50. Le courant il est représenté en ordonnées, l'axe 50 correspond à l'axe des abscisses. Le simulateur calcule la valeur is du courant il correspondant au point Ps à partir du modèle simplifié tel que défini ci- dessus. Le simulateur calcule ensuite la valeur iC du courant il correspondant au point PC à partir du modèle complet. La portion de droite de la surface plane elliptique 40 passant par les valeurs de courant et is est tangente à une portion circulaire de la demi-sphère 31 représentée en figure 9. Les valeurs de courant représentatives du modèle mixte pour les couples de tensions (VN1,VN2) de la région 45 situés sur l'axe 50 entre les points Ps et PC forment une portion de droite 51 passant par les valeurs de courant is et iC. Le calcul du courant il correspondant au point Pi est alors effectué à partir de l'équation de la portion de droite 51. On obtient un courant ij appartenant à la portion de droite 51. De façon générale, les valeurs de courant représentatives du modèle mixte pour les couples de tensions (VN1, VN2) appartenant à l'axe 50, peuvent former n'importe quel type de courbe entre les valeurs de courant is et iC. Afin que la série d'itérations de calcul du courant il converge plus facilement et plus rapidement, on pourra définir une courbe telle que les dérivées du courant il au point Pg soient identiques pour le modèle simplifié et mixte et telle que les dérivées du courant il au point PC soient identiques pour les modèles mixte et complet. A la suite ou parallèlement au calcul du courant il, on effectue de façon similaire le calcul des courants i2 et i3. 35 Après vérification de la loi des mailles au noeud N2, on effectue ou non une nouvelle itération de calcul pour déterminer les tensions/courants à l'instant de simulation tj. Dans l'exemple du composant Z1, lors d'une nouvelle itération de calcul pour l'instant tj, on détermine si le nouveau couple de tensions (VN1,VN2) estimé probable est éloigné ou non du couple de tensions (VN1/VN2)j-x déterminé pour l'instant de simulation tj_x choisi préalablement à la première itération de calcul. L'éloignement est apprécié selon un procédé identique à celui décrit précédemment en relation à la figure 6. On calcule ensuite le courant il comme précédemment décrit à partir d'un des trois modèles simplifié, mixte ou complet. A la suite ou parallèlement, on effectue de même le calcul des courants i2 et i3 pour cette nouvelle itération de calcul. Après vérification de la loi des mailles, on continue la série d'itérations selon un procédé identique jusqu'à ce que la loi des mailles soit vérifiée. De façon générale, le procédé de la présente invention s'applique à tout type de circuit intégré comprenant plusieurs composants élémentaires. Chaque composant a une ou plusieurs bornes, de façon générale n bornes, reliées chacune à un noeud du circuit. Pour chaque composant, un modèle complet, permet de calculer de façon très précise chacun des courants entrant ou sortant du composant par l'une de ses bornes en fonction d'un jeu de n tensions aux noeuds auxquels est relié le composant. Chaque courant entrant sur une des bornes du composant peut être représenté par une hypersurface dans un espace à n+l dimensions en fonction des n tensions aux bornes du composant. Préalablement à la première itération de calcul effectuée pour déterminer les tensions et courants du circuit à un instant donné tj, on définit pour chaque composant un instant antérieur de simulation tj_x pour lequel le modèle complet a été utilisé pour calculer les valeurs de courant correspondant au jeu de tensions déterminé pour cet instant antérieur. On définit ensuite pour chaque composant un espace de jeu de tensions possibles à n dimensions, puis on calcule un écart entre un point Pj_x correspondant au jeu de n tensions déterminé pour l'instant tj_x et un point Pj correspondant au jeu de n tensions estimé probable pour l'instant tj. On détermine ensuite pour chaque composant si l'écart calculé est inférieur à un premier seuil Vsl, est compris entre le premier seuil Vsl et un second seuil Vs2 supérieur à Vsl, ou est supérieur au second seuil Vs2. Selon le cas, on utilise respectivement un modèle simplifié, un modèle mixte ou un modèle complet pour calculer chacun des courants entrant ou sortant du composant considéré. Pour chaque courant entrant ou sortant d'un composant donné, les ensembles de jeux de tensions éloignés du point Pj_x d'une valeur de seuil Vsl et d'une valeur de seuil Vs2 constituent respectivement des première et seconde hypersphères "limites". L'ensemble des valeurs de courant représentatif du modèle simplifié correspond à une portion d'un hyperplan tangent à l'hypersurface représentative du modèle complet pour la valeur de courant correspondant au point Pj _x. Chacun des courants de la portion de l'hyperplan correspond à un jeu de tensions situé à l'intérieur de la première hypersphère limite. Dans le cas où le point Pj est situé à l'intérieur ou sur cette première hypersphère limite, le calcul du courant est effectué à partir de l'équation de cet hyperplan "tangent". L'ensemble des valeurs de courant représentatif du modèle mixte correspond à une hypersurface de "raccordement" entre la portion d'hyperplan susmentionnée et un ensemble de valeurs de courant de l'hypersurface représentative du modèle complet correspondant aux jeux de tensions de la seconde hypersphère limite. Chacun des points de l'hypersurface de "raccordement" correspond à un jeu de tensions situé dans une région délimitée par les première et seconde hypersphères. Le jeu d'équations permettant de connaître la valeur d'un courant de l'hypersurface de "raccordement" correspondant à un jeu de tensions donné de la région susmentionnée peut être défini partiellement pour une portion de la région selon un procédé similaire à celui décrit en relation aux figures 8 et 9. Une fois que tous les courants du circuit ont été calculés, on vérifie la loi des mailles pour chacun des noeuds du circuit. Dans le cas où la loi des mailles n'est pas vérifiée, on procède à une nouvelle itération de calcul pour l'instant tj. En conservant pour chaque composant le même point de "référence" Pj_x que celui défini pour la première itération, on détermine comme précédemnent pour chaque composant si le nouveau jeu de tensions estimé probable est éloigné ou non du point Pj_x. Selon l'éloignement on utilise un des trois modèles simplifié, mixte ou complet pour calculer les courants entrant dans chaque composant. Le choix pour chaque composant d'un instant antérieur de simulation tj_x pour lequel le modèle complet a été utilisé pour calculer les valeurs de courant correspondant au jeu de tensions déterminé pour cet instant antérieur peut être effectué de diverses façons. On pourra choisir l'instant antérieur le plus proche de l'instant donné tj pour lequel le modèle complet a été utilisé. On pourra aussi choisir, parmi les instants antérieurs pour lesquels le modèle complet a été utilisé, celui pour lequel le jeu de tensions déterminé est le plus proche du jeu de tensions estimé probable pour l'instant tj. Les tensions de seuil Vsl et Vs2 peuvent être définies pour chaque type de composant. Dans le cas d'une diode, on pourra définir des tensions de seuil Vsl et Vs2 respectivement égales à une ou deux fois la tension de seuil intrinsèque de la diode, soit par exemple quelques centaines de millivolts. Dans le cas d'une résistance on pourra éventuellement choisir des tensions de seuil "infinies" étant donné que les modèles simplifié et mixte sont confondus avec le modèle complet quelque soient les valeurs de tension. Bien entendu, la présente invention est susceptible de diverses variantes et modifications qui apparaîtront à l'homme de l'art. En particulier, l'homme de l'art pourra définir d'autres procédés d'estimation de l'écart des tensions aux noeuds auxquels est relié un composant entre deux instants de simulation. REVENDICATIONS 1. Procédé de simulation d'un circuit électrique comprenant des composants élémentaires (Z1, Z2, Z3) dont chacune des bornes est reliée à un noeud du circuit (Ni, N2, N3), certains noeuds recevant des stimuli externes consistant en une suite de valeurs de tensions (VN1, VN3) appliquées à des instants successifs de simulation, dans lequel la détermination à un instant donné des tensions (VN2) aux' noeuds du circuit consiste à effectuer une ou plusieurs itérations de calcul, chaque itération consistant à définir une tension probable pour chaque noeud, à calculer les courants entrant ou sortant dans chacun des composants en utilisant pour chaque composant un modèle de calcul des courants en fonction des tensions aux bornes du composant, puis à vérifier la loi des mailles pour chaque noeud, la suite d'itérations continuant jusqu'à ce que la loi des mailles soit vérifiée, et dans lequel pour le premier instant de simulation (tl) suivant l'instant initial de simulation (t0), lors d'une itération de calcul, un courant entrant ou sortant d'un composant donné est calculé à partir de l'un de trois modèles, simplifié, mixte ou complet, le modèle complet permettant de connaître très précisément la valeur du courant en fonction des tensions aux bornes du composant, le modèle simplifié correspondant à une linéarisation tangentielle du modèle complet autour de la valeur du courant correspondant aux tensions aux bornes du composant déterminées pour l'instant initial, le modèle mixte constituant un raccordement des modèles simplifié et complet en assurant la continuité des courants, le choix d'un modèle étant fonction de l'écart entre les tensions aux bornes du composant déterminées pour l'instant initial de la simulation et celles jugées pro- bables pour le premier instant de simulation, de sorte que pour chaque composant: -- si l'écart est inférieur à un premier seuil (Vs?) défini pour ce composant, le modèle simplifié est utilisé ; - - si l'écart est supérieur à un second seuil (Vs2) défini pour ce composant, ce second seuil étant supérieur au premier seuil, le modèle complet est utilisé ; et - - si l'écart est compris entre les premier et second seuils pour ce composant, le modèle mixte est utilisé. 2. Procédé selon la revendication 1, dans lequel pour un instant de simulation donné (tj), la première itération de calcul comprend, préalablement au calcul des courants, une étape supplémentaire consistant à choisir pour chaque composant un instant de simulation (t3.-x) antérieur à l'instant donné pour lequel les courants entrant ou sortant du composant ont été calculés à partir du modèle complet {31) du composant, le calcul d'un courant (il) entrant ou sortant d'un composant donné (Zi), lors de la première itération et des suivantes, étant effectué à partir de l'un de trois modèles, simplifié, mixte ou complet, le modèle simplifié étant défini autour de la valeur du courant (ij_x) correspondant aux tensions ((VN1,VN2)j-x) aux bornes du composant déterminées pour ledit instant de simulation antérieur, le choix d'un modèle étant fonction de l'écartement entre les tensions aux bornes du composant déterminées pour l'instant antérieur choisi et celles jugées probables ((VN1,VN2)j) pour l'instant de simulation donné comme pour le premier instant. 3. Procédé selon la revendication 2, dans lequel, pour un composant à n bornes, ledit écart est égal à la distance entre des premier et second points d'un espace à n dimensions, le premier point ayant pour coordonnées les tensions aux bornes du composant jugées probables ((VN1, VN2)j) pour l'instant donné, et le second point ayant pour coordonnées les tensions aux bornes du composant ((VN1,VN2)j-x) déterminées pour l'instant antérieur choisi (tj_x). 4. Procédé selon la revendication 3, dans lequel: l'ensemble des valeurs possibles d'un courant (i2) entrant ou sortant d'un composant donné (21) correspondant au 35 modèle complet du composant constitue une première hypersurface (31) représentée dans un espace à n+l dimensions comprenant un axe des courants et ledit espace de tensions, les ensembles de points, dudit espace de tensions, éloignés du second point (Pi_x) d'une valeur égale à la première tension de seuil (Vsl) ou égale à la seconde tension de seuil (Vs2) constituent hypersphères, l'ensemble modèle simplifié du hyperplan tangent à (i3. _x) correspondantdes première (35) et seconde (36) des valeurs de courant représentatives du composant constitue une portion (40) d'un la première hypersurface pour un courant au second point, chacun des courants de la portion d'hyperplan correspondant à un point situé à l'intérieur de la première hypersphère, l'ensemble des valeurs de courant représentatives du modèle mixte du composant constitue une seconde hypersurface (42) de raccordement entre la portion d'hyperplan et un sous- ensemble de valeurs de courant (41) de la première hypersurface correspondant aux points de la seconde hypersphère. 5. Procédé selon la revendication 4, dans lequel le 20 modèle mixte d'un composant est défini partiellement de la façon suivante lorsque ledit écart est compris entre les premier et second seuils (Vsl, Vs2) : définir des troisième et quatrième points (Ps, PC) correspondant à l'intersection entre un axe (50) ayant pour 25 origine le second point (Pi_x) et passant par le premier point (Pi) et les première et seconde hypersphères, - calculer des troisième et quatrième courants (is, i0) correspondant aux troisième et quatrième points à partir respectivement des modèles simplifié (40) et complet (31) du 30 composant, et déterminer un troisième ensemble de valeurs de courant pour chacun des points dudit axe compris entre les troisième et quatrième points, le troisième ensemble de valeurs de courant reliant les troisième et quatrième valeurs de courant, la valeur du courant calculé pour le premier point (Pi) étant la valeur de courant (is) du troisième ensemble correspondant au premier point. 6. Procédé selon la revendication 2, dans lequel ledit instant de simulation antérieur (tj _x) est l'instant le plus proche de l'instant donné (tj) pour lequel les courants entrant ou sortant du composant ont été calculés à partir du modèle complet du composant. 7. Procédé selon la revendication 2, dans lequel ledit instant de simulation antérieur (tj_x) est tel que les tensions aux bornes du composant déterminées pour cet instant antérieur soient les plus proches des tensions aux bornes du composant jugées probables pour l'instant donné et tel que pour les tensions aux bornes du composant déterminées pour cet instant, les courants entrant ou sortant du composant correspondant aient été calculés à partir du modèle complet du composant.