Ngspice to PSpice conversion

captainkirk86

Hi All,

I have a piece of Ngspice code that I need to convert to PSpice. I feel it is a simple matter of syntax but I have so far failed to make it work.

The code models a Memristor with thresholding ability (the fourth fundamental component along with resistors, capacitors, and inductors)

The code is as follows:

.subckt memristor plus minus PARAMS: Ron=1K Roff=10K Rinit=5K alpha=0 beta=1E13 Vt=4.6
Bx 0 x I='(f1(V(plus)-V(minus))>0) && (V(x)<Roff) ? {f1(V(plus)-V(minus))}: (f1(V(plus)-V(minus))<0) && (V(x)>Ron) ? {f1(V(plus)-V(minus))}: {0}'
Cx x 0 1 IC={Rinit}
Rmem plus minus r={V(x)}
.func f1(y)={beta*y+0.5*(alpha-beta)*(abs(y+Vt)-abs(y-Vt))}
.ends


The code explanation is as follows:
- The memristive system is realized as a sub-circuit combining a behavioral resistor R, a current source ↑, and a capacitor C.
- The second line of code (Bx ...) defines the current source with the current specified through ternary functions. (A ternary function is defined in the code as a ? b : c , which means ”IF a, THEN b, ELSE c” [17].) The purpose of these functions is to limit RM between Ron and Roff.
- The third line of the code specifies the capacitor C (Cx ...) with an initial condition.
- The fourth line (Rmem ...) defines the behavioral resistor whose resistance takes the same numerical value as the voltage across the capacitor.
- The next line (.func ...) provides the function f according to Equation (not important here)


If you can help me in any way that would be great!!

Many thanks
Captainkirk

Kholdstare

I can help you. But can you be a little more clear with simple circuit diagram and simple expressions?

captainkirk86

Hi there, thanks for replying.

I've attached the circuit diagram as image.

PLEASE NOTE - I dont need help with the Model, I just need help with syntax. Thats the code I gave in my first post just needs its syntax tweaked so it is compatible with Pspice.

The research paper that I'm referring to states the following for the Spice model:

"SPICE model: The memristive devices with threshold that we consider are those described by
I = X^-1V_M
X ̇ = (βV_M +0.5(α−β)[|V_M +Vt|−|V_M −Vt|])

where Vt is the threshold voltage, R_on and R_off are limiting values of the memristance R_M ≡ X, and the θ-functions

Note - These equations are represented in ngspice code I gave in my first post.

- The important model parameters are the coefficients α and β that characterize the rate of memristance change at |V_M|<Vt and |V_M|>Vt , respectively.
- These two coefficients define the slopes of the f (V_M) curve below and above the threshold. When α = 0, the device state changes only if |V_M| > Vt"

Please let me know if this makes it clearer.

I appreciate any help you can provide here!

Many thanks again
CaptainKirk

captainkirk86

attached file

Attached Thumbnails

 

Kholdstare

Captainkirk, I really did not try to understand the model. Actually I could not make out the ngspice syntax. I am working on it. You'll get your spice model after sometime.

h_vogt

Kholdstare,

see the ngspice-users forum (I am not allowed to give the link here).

captainkirk is engaging at least three forums ...

But is was nice just to get it up and running.

Regards

Holger

Kholdstare

Here's your PSpice macromodel as requested. To use it in your circuit just copy the "* memres ... begins" to "* memres ... ends" portion to your circuit and call the model using "X(ref) (node1) (node2) R_MEMR PARAMS: VCUT=(param) ALPHA=(param) BETA=(param)" device card.

The test parameters I used are Vcut = 4.6V, alpha = 0, beta = 1e-3
The test result is in the attachment.

I simplified your expression as written below,
I = a*V for |V| < Vt
= a*V+(b-a)(V-Vt) for V>Vt
= a*V+(b-a)(V+Vt) for V<Vt

I could not model capacitance as you did not give its expression clearly.

Attached Thumbnails

 

Kholdstare

Sorry for the late. I had actually forgotten about it after your last post. Shows kinda lazy person I am. haha

captainkirk86

Kholdstare

Thank you were much with your reply!!

I can hardly call you lazy, after all you did come up with something very interesting.

I'm very impressed with your model. Its is different from what I came up with but I still like the output graph that you attached!

Thank you once again

Kind regards,
Kirk