世纪电源网社区logo
社区
Datasheet
标题
返回顶部
未解决

自己建立的saber模型出现问题a**x with "a" less than or equal to zero in

[复制链接]
查看: 1270 |回复: 3
1
wanerlover
  • 积分:230
  • |
  • 主题:10
  • |
  • 帖子:28
积分:230
LV3
助理工程师
  • 2014-7-8 09:09:15
10问答币
*** ERROR "MAST_POWER_OP1" *** a**x with "a" less than or equal to zero in这个问题怎么解决?求助大神
andy_zhou_nuaa
  • 积分:2029
  • |
  • 主题:3
  • |
  • 帖子:120
积分:2029
LV8
副总工程师
  • 2014-7-8 09:41:45
  • 倒数3
 
附上你的MAST源码,并简要介绍一你的模型;这样才好分析。
wanerlover
  • 积分:230
  • |
  • 主题:10
  • |
  • 帖子:28
积分:230
LV3
助理工程师
  • 2014-7-9 16:08:13
  • 倒数2
 
template sct2080ke m_1 m_2 m_3
electrical m_1,m_2,m_3
{
foreign sstmath
foreign sstput,sstget
electrical m_11,m_12,m_21,m_22,m_31,m_32,m_103,m_104,m_41,m_42,m_43,m_111,m_112,
m_51,m_52
external number temp
number t0=25,sstputrv

function (y) = r1(i)
number y,i
{
number temp,t0
temp = sstget(instance(),"temp")
t0 = sstget(instance(),"t0")
y=(((0.05956*i)*(1+(0.5563*(limexp(((temp-t0)/128.7))-1))))+(((4.3e-006*i)*(abs(i)**(1.898/(1+(0.5341*(limexp(((temp-t0)/(-17.66)))-1))))))*(1+(0.9978*(limexp(((temp-t0)/(-6.998)))-1)))))
} # end function r1
function (y) = v1(v,w)
number y,v,w
{
number temp,t0
temp = sstget(instance(),"temp")
t0 = sstget(instance(),"t0")
y=((v-((0.08701*w)/(1+(0.2299*(limexp(((temp-t0)/(-71.48)))-1)))))-(0.3955*asinh(((w/0.3073)/(1+(0.1881*(limexp(((temp-t0)/(-47.01)))-1)))))))
} # end function v1
function (y) = v2(v)
number y,v
{
number f_0_0
number temp,t0
temp = sstget(instance(),"temp")
t0 = sstget(instance(),"t0")
f_0_0=sstmath("psp_if",(v>0),
((2.053e-006*(v**7.933))/(1+(0.9745*(limexp(((temp-t0)/(-32.48)))-1)))),
0)
y=f_0_0
} # end function v2
function (y) = i1(v,w)
number y,v,w
{
y=(((((v*(abs(w)+1000))/1010)*1.2)*w)/(abs(w)+2))
} # end function i1
function (y) = c1(v)
number y,v
{
number f_0_0
f_0_0=sstmath("psp_if",(v>(-1.891)),(1526+(193*v)),
(1905*((1-(v/1.119))**(-0.5006))))
y=f_0_0
} # end function c1
function (y) = c2(v)
number y,v
{
number f_0_0
f_0_0=c1(v)
y=((f_0_0*((0.4187*tanh(((v+1.426)*1.264)))+0.5813))-(0.005147*v))
} # end function c2
function (y) = r101(i)
number y,i
{
number temp,t0
temp = sstget(instance(),"temp")
t0 = sstget(instance(),"t0")
y=(((1.395*asinh((i/6.151)))*limexp(((temp-t0)/1035)))+((0.0004045*i)*abs(i)))
} # end function r101
function (y) = i101(v)
number y,v
{
number f_0_0
number temp,t0
temp = sstget(instance(),"temp")
t0 = sstget(instance(),"t0")
f_0_0=sstmath("psp_if",(v>0),
((0.0003094*(v**(10.21*limexp(((temp-t0)/2494)))))*limexp(((temp-t0)/43.22))),
0)
y=f_0_0
} # end function i101
function (y) = i102(v)
number y,v
{
number f_0_0
number temp,t0
temp = sstget(instance(),"temp")
t0 = sstget(instance(),"t0")
f_0_0=sstmath("psp_if",(v<0),
((((1.528e-009*(limexp((v/10))-1))*limexp(((-v)/969.2)))*limexp(((temp-t0)/267.4)))*(limexp((((-v)-(1720*limexp(((temp-t0)/1250))))/10))+1)),
0)
y=f_0_0
} # end function i102
function (y) = c101(v)
number y,v
{
number f_0_0
f_0_0=sstmath("psp_if",(v>0.5735),(896.6+(988.5*v)),
(1119*((1-(v/1.147))**(-0.3874))))
y=((f_0_0-50.2)-(0.02923*v))
} # end function c101
val nu l_0_0,l_0_1,l_1_0,l_2_0,l_2_1,l_3_0,l_3_1,l_4_0,l_5_0,l_6_0,l_6_1,l_7_0,
l_7_1,l_7_2,l_7_3,l_7_4,l_7_5
var i ri_0,ri_1,ri_2,ri_3,ri_4,ri_5,ri_7,ri_8,ri_10,ri_11,ri_12
val i li_6,li_9,li_13

parameters {
sstputrv = sstput(instance(),"temp",temp)
sstputrv = sstput(instance(),"t0",t0)
sstputrv = sstput(instance(),"temp",temp)
sstputrv = sstput(instance(),"t0",t0)
sstputrv = sstput(instance(),"temp",temp)
sstputrv = sstput(instance(),"t0",t0)
sstputrv = sstput(instance(),"temp",temp)
sstputrv = sstput(instance(),"t0",t0)
sstputrv = sstput(instance(),"temp",temp)
sstputrv = sstput(instance(),"t0",t0)
sstputrv = sstput(instance(),"temp",temp)
sstputrv = sstput(instance(),"t0",t0)
}
pspr.1 m_12 m_32 = r=1e+015
v.1 m_1 m_11 = dc=0
v.2 m_2 m_21 = dc=0
pspc.1 m_12 m_32 = c=1e-012
pspc.2 53 0 = c=1e-012
l_0_0=sstmath("psp_limit",i(v.1),(-1000000),1000000)
l_0_1=r1(l_0_0)
i(m_11->m_12) += ri_0 ; ri_0 : v(m_11,m_12) = l_0_1 # e1
i(m_21->m_22) += ri_1 ; ri_1 : v(m_21,m_22) = (i(v.2)*11.5) # e2
v.3 m_3 m_31 = dc=0
i(m_31->m_32) += ri_2 ; ri_2 : v(m_31,m_32) = (i(v.3)*0) # e3
l_1_0=sstmath("psp_limit",(v(m_22)-v(m_32)),0,22)
i(m_41) += ri_3 ; ri_3 : v(m_41) = l_1_0 # e4
l_2_0=sstmath("psp_limit",v(m_43),0,200)
l_2_1=v1(v(m_41),l_2_0)
i(m_42) += ri_4 ; ri_4 : v(m_42) = l_2_1 # e5
l_3_0=sstmath("psp_limit",v(m_42),0,20)
l_3_1=v2(l_3_0)
i(m_43) += ri_5 ; ri_5 : v(m_43) = l_3_1 # e6
l_4_0=i1(v(m_43),(v(m_12)-v(m_32)))
li_6=l_4_0
i(m_12->m_32) += li_6 # g1
i(m_51) += ri_7 ; ri_7 : v(m_51) = (v(m_22)-v(m_1)) # e7
i(m_52) += ri_8 ; ri_8 : v(m_52) = (v(m_22)-v(m_1)) # e8
v.4 m_52 53
l_5_0=c2(v(m_51))
li_9=(i(v.4)*l_5_0)
i(m_22->m_1) += li_9 # g2
pspc.3 m_22 m_32 = c=2.064e-009
pspr.2 m_22 m_32 = r=1000000000
pspc.101 113 0 = c=1e-012
#******** ********* ********* ********* ********* ********* ********* *********
# ********* *********
v.101 m_3 m_103 = dc=0
l_6_0=sstmath("psp_limit",i(v.101),(-1000000),1000000)
l_6_1=r101(l_6_0)
i(m_103->m_104) += ri_10 ; ri_10 : v(m_103,m_104) = l_6_1 # e101
i(m_111) += ri_11 ; ri_11 : v(m_111) = (v(m_104)-v(m_1)) # e102
i(m_112) += ri_12 ; ri_12 : v(m_112) = (v(m_104)-v(m_1)) # e103
v.102 m_112 113 = dc=0
l_7_0=sstmath("psp_limit",v(m_111),0,20)
l_7_1=i101(l_7_0)
l_7_2=sstmath("psp_limit",v(m_111),(-3000),0)
l_7_3=i102(l_7_2)
l_7_4=sstmath("psp_limit",v(m_111),(-3000),20)
l_7_5=c101(l_7_4)
li_13=((l_7_1+l_7_3)+(i(v.102)*l_7_5))
i(m_104->m_1) += li_13 # g101
} # end sct2080ke

这是sct2080ke的代码,是一个mosfet。转换为saber模型后就出现了这个问题
andy_zhou_nuaa
  • 积分:2029
  • |
  • 主题:3
  • |
  • 帖子:120
积分:2029
LV8
副总工程师
最新回复
  • 2014-7-11 13:43:47
  • 倒数1
 
不好意思,已超出我能力范围;爱莫能助。
热门技术、经典电源设计资源推荐

世纪电源网总部

地 址:天津市南开区黄河道大通大厦8层

电 话:400-022-5587

传 真:(022)27690960

邮 编:300110

E-mail:21dy#21dianyuan.com(#换成@)

世纪电源网分部

广 东:(0755)82437996 /(138 2356 2357)

北 京:(010)69525295 /(15901552591)

上 海:(021)24200688 /(13585599008)

香 港:HK(852)92121212

China(86)15220029145

网站简介 | 网站帮助 | 意见反馈 | 联系我们 | 广告服务 | 法律声明 | 友情链接 | 清除Cookie | 小黑屋 | 不良信息举报 | 网站举报

Copyright 2008-2024 21dianyuan.com All Rights Reserved    备案许可证号为:津ICP备10002348号-2   津公网安备 12010402000296号