[SOLVED] Pspice memristorproblem

Pspice memristor problem



Hi all,
I am facing problem :
**** INCLUDING SCHEMATIC1.net ****
* source SPICE_ABDALLA
X_AB_MEM1 N00244 0 modelmemristor PARAMS: PHIO=0.95 LM=0.0998 W1=0.1261
+ FOFF=3.5E-6 IOFF=115E-6 AOFF=1.2 FON=40E-6 ION=8.9E-6 AON=1.8 B=500E-6
+ WC=107E-3
V_V1 N00089 0
+SIN 0.2 0.1 1 0 0 0
R_R1 N00089 N00244 2.4k
**** RESUMING ab_mem.cir ****
.END
INTERNAL ERROR -- Invalid Number in device X_AB_MEM1.C1, Divide

ABORTING SIMULATION
**** 09/08/13 13:22:28 ******* PSpice 16.0.0 (July 2006) ****** ID# 0 *****
** Profile: "SCHEMATIC1-ab_mem" [ C:\OrCAD\OrCAD_16.0\tools\capture\library\pspice\memristors\spice_model_abdalla\SPICE_ABDALLA-PSp
**** JOB STATISTICS SUMMARY
******************************************************************************

Node counts:
Top level (NUNODS) = 3
External (NCNODS) = 17
Total (NUMNOD) = 17

Total device count (NUMEL) = 20
Capacitors (C) = 1
VCVS (E) = 12
VCCS (G) = 3
Resistors (R) = 3
Voltage Sources (V) = 1
Number of subcircuits (X) = 1
Matrix statistics:
Matrix size (NSTOP) = 30
Initial no. elements (NTTAR) = 74
No. elements w/ fillin (NTTBR) = 74
No. fillins (IFILL) = 0
No. overflows (NTTOV) = 0
No. LU operations (IOPS) = 0
Percent sparsity (PERSPA) = 91.778

Analysis statistics:
No. total time points (NUMTTP) = 0
No. rejected time points (NUMRTP) = 0
No. iterations (NUMNIT) = 0
ERROR -- Disk write error. The disk may be full.

ABORTING SIMULATION
**** 09/08/13 13:22:28 ******* PSpice 16.0.0 (July 2006) ****** ID# 0 ********
** Profile: "SCHEMATIC1-ab_mem" [ C:\OrCAD\OrCAD_16.0\tools\capture\library\pspice\memristors\spice_model_abdalla\SPICE_ABDALLA-PSp
**** JOB STATISTICS SUMMARY
CODE (subcircuit) of the device memristor is from pepar :"SPICE Modeling of Memristors ::by Hisham Abdalla "
*Abdalla model
.SUBCKT modelmemristor plus minus PARAMS:
+phio=0.95 Lm=0.0998 w1=0.1261 foff=3.5e-6
+ioff=115e-6 aoff=1.2 fon=40e-6 ion=8.9e-6
+aon=1.8 b=500e-6 wc=107e-3
G1 plus internal value={sgn(V(x))*(1/V(dw))^2*0.0617*(V(phiI)*exp(-V(B)*V(sr))-(V(phiI)+abs(V(x)))*exp(-V(B)*V(sr2)))}
Esr sr 0 value={sqrt(V(phiI))}
Esr2 sr2 0 value={sqrt(V(phiI)+abs(V(x)))}
Rs internal minus 215
Eg x 0 value={V(plus)-V(internal)}
Elamda Lmda 0 value={Lm/V(w)}
Ew2 w2 0 value={w1+V(w)-(0.9183/(2.85+4*V(Lmda)-2*abs(V(x))))}
EDw dw 0 value={V(w2)-w1}
EB B 0 value={10.246*V(dw)}
ER R 0 value={(V(w2)/w1)*(V(w)-w1)/(V(w)-V(w2))}
EphiI phiI 0 value={phio-abs(V(x))*((w1+V(w2))/(2*V(w))) - 1.15*V(Lmda)*V(w)*((log(V(R)))/V(dw))}
C1 w 0 1e-9 IC=1.2
R w 0 1e8MEG
Ec c 0 value={abs(V(internal)-V(minus))/215}
Emon1 mon1 0 value={((V(w)-aoff)/wc)-(V(c)/b)}
Emon2 mon2 0 value={(aon-V(w))/wc-(V(c)/b)}
Goff 0 w value={foff*sinh(stp(V(x))*V(c)/ioff)*exp(-exp(V(mon1))-V(w)/wc)}
Gon w 0 value={fon*sinh(stp(-V(x))*V(c)/ion)*exp(-exp(V(mon2))-V(w)/wc)}
.ENDS modelmemristor
Thanks for any help in advance .

My guess is that it doesn't like the IC=1.2 on the line for C1. Check the Pspice syntax for initial conditions and edit it. There are also quite a few old discussions on EDAboard about memristors - try the search function.

Keith

Thanks Keith for your response

I have changed code line

C1 w 0 1e-9
.ic V(w)=1.2
Which means initial voltage as node w will be 1.2 volt .
please read this paper carefully otherwise it would be difficult to simulate .
result is same as paper us piece wise linear voltage source .