GEKKO使用B样条时出现访问冲突错误求助

算力阿峰

2026-7-9

问题描述

搭建了一个带有流体状态初始值的模型，希望用CoolProp计算状态变量（此为之前相关问题的后续）。尝试计算压力和温度时，出现Exception: Access Violation错误，无法定位问题，请求协助解决。

代码示例

m = GEKKO(remote=False)
m.time = np.arange(1,3,0.5)
p_scale = 1e5
p_range = np.arange(0.9,2.5, 0.1) #bar
D_range = np.arange(PropsSI('D','P',p_range[-1]*p_scale,'Q',0,"Nitrogen"),PropsSI('D','P',p_range[-1]*p_scale,"T",75,"Nitrogen"), 0.5)
T_range = np.arange(75,PropsSI('T','P',p_range[-1]*p_scale,"Q",0,"Nitrogen"), 0.1) 
u_range = np.arange(PropsSI('U','P',p_range[-1]*p_scale,'T',75,"Nitrogen")/1e3,
                    (PropsSI('D','P',p_range[-1]*p_scale,"Q",0,"Nitrogen")+1)/1e3, 1) #kJ/kg   
p_data,T_data = [], []
for n in D_range:
    row_p,row_T = [],[]
    for i in T_range:
        row_p.append(PropsSI('P','D',n,'T',i,"Nitrogen")/p_scale)
    for i in u_range: row_T.append(PropsSI('T','D',n,'U',i*1e3,"Nitrogen"))
    p_data.append(row_p)
    T_data.append(row_T)
T = m.Var(value=80)
p = m.Var()
D=m.Param(value=793.97)
u=m.Param(value=-116.57)
m.bspline(D,u,T,D_range,u_range,T_data) #T
m.bspline(D,T,p,D_range,T_range,p_data) #p
m.options.SOLVER = 3 # solver (IPOPT) 
m.options.IMODE = 4
m.solve(disp=False)