2

I have the following circuit:

schematic

simulate this circuit – Schematic created using CircuitLab

Where \$\text{V}_\text{i}\left(t\right)=\hat{\text{u}}\sin\left(\omega t+\varphi\right)\$ and the relation between the current through and voltage across the diode is given by the Shockley diode equation:

$$\text{I}_\text{D}=\text{I}_\text{S}\left(\exp\left(\frac{\text{q}\text{V}_\text{D}}{\eta\text{k}\text{T}}\right)-1\right)\tag2$$

Question: What is the mathematical function that describes the voltage \$\text{V}_1\$ using Ohm's law and the Shockley diode equation?

JRE
  • 67,678
  • 8
  • 104
  • 179
Jan Eerland
  • 7,203
  • 12
  • 21
  • 1
    This is a few lines that asks for mountains of effort in reply. The site is not for free personal tutoring. Please edit your question and detail you already know, show all that you have discovered for yourself on the subject. – TonyM Feb 02 '22 at 14:02
  • 1
    \$I_D = \dfrac{V_1}{R_1}\$ sounds like a good place to start. – Andy aka Feb 02 '22 at 14:07
  • Given the history, Jan might soon post an excellent answer. – Justme Feb 02 '22 at 14:13
  • 1
    This is a nonlinear circuit so you can use iteration (numerical solution) or you could try to use the Lambert W function. Is this what you want? – G36 Feb 02 '22 at 14:57
  • @G36 Yes, that is what I want. – Jan Eerland Feb 02 '22 at 14:58
  • 1
    Did you see this https://en.wikipedia.org/wiki/Diode_modelling#Large-signal_modelling and this one https://paklaunchsite.jimdofree.com/app/download/7813805654/PAK-Course-101-D%2BR-W-analysis-ALL.pdf?t=1618348827 – G36 Feb 02 '22 at 15:08
  • 1
    @JanEerland Don't know if [this helps](https://electronics.stackexchange.com/a/592785/38098), or not. Might be. From the diode current you can use that to multiply by the external resistance and subtract from the supply voltage to find the diode voltage, for example. – jonk Feb 02 '22 at 17:31

1 Answers1

1

It is given that

\$V_i=u\sin(\omega t +\varphi)\$ and \$I_D=I_S \left[ \exp \left( \frac{qV_D}{\eta k T} \right)-1 \right]\$

From the circuit, we can say: \$V_1 = V_i-V_D\$ and \$I_D=\frac{V_1}{R_1}\$

So, we get:

\$\frac{V_1}{R_1}=I_S \left[ \exp \left( \frac{qV_D}{\eta k T} \right)-1 \right]\$

\$\Rightarrow \frac{V_1}{R_1I_S}= \exp \left( \frac{qV_D}{\eta k T} \right)-1\$

Defining \$C_1 := \frac{1}{R_1I_S}\$ and \$C_2:=\frac{q}{\eta k T}\$, we get

\$C_1V_1= \exp \left( C_2V_D \right)-1\$

\$\Rightarrow \frac{V_1}{C_1}+1=\exp(C_2V_D)\$ \$\Rightarrow \ln \left( \frac{V_1}{C_1}+1 \right)=C_2V_D\$ \$\Rightarrow V_D = \frac{1}{C_2} \ln \left( \frac{V_1}{C_1}+1 \right)\$

Hence voltage relation can be re-written as:

\$V_1 = V_i - \frac{1}{C_2} \ln \left( \frac{V_1}{C_1}+1 \right)\$

\$\Rightarrow C_2V_1 = C_2V_i - \ln \left( \frac{V_1}{C_1}+1 \right)\$

\$\Rightarrow \ln \left( \frac{V_1}{C_1}+1 \right) = C_2V_i - C_2V_1 \$

\$\Rightarrow \frac{V_1}{C_1}+1 = \exp \left( C_2V_i - C_2V_1 \right) \$

\$\Rightarrow \frac{V_1}{C_1}+1 = \exp \left( C_2V_i \right) \exp \left(- C_2V_1 \right) \$

\$\Rightarrow \frac{V_1}{C_1} = \exp \left( C_2V_i \right) \exp \left(- C_2V_1 \right) -1 \$

\$\Rightarrow V_1 = C_1\exp \left( C_2V_i \right) \exp \left(- C_2V_1 \right) - C_1 \$

The general solution using Lambert W function is given by (see Wikipedia for quick reference and a book for compelte reference): \$x = a+b\exp \left( cx \right) \Rightarrow x=a-\frac{1}{c}W \left( -bc\exp(ac) \right)\$

Using the above solution considering \$V_1=x\$, \$a= -C_1\$, \$b=C_1\exp(C_2V_i)\$, and \$c=-C_2\$, we obtain

\$V_1=-C_1+\frac{1}{C_2}W \left( C_1C_2\exp(C_2V_i) \exp(C_1C_2) \right)\$

And hence:

\$V_1=-C_1+\frac{1}{C_2}W \left( C_1C_2\exp(C_2u\sin(\omega t +\varphi)) \exp(C_1C_2) \right)\$

SKGadi
  • 862
  • 3
  • 13
  • Plotting this gives the same voltage as the input which can't be right. – Jan Eerland Feb 02 '22 at 17:17
  • What values are you using for the numerical simulation? Remember when the diode is in the forward bias, the circuit can be approximated by replacing the voltage drop by a diode is same as its knee voltage. This is how the LED resistance is calculated. – SKGadi Feb 02 '22 at 17:21