107 lines
4.7 KiB
Python
107 lines
4.7 KiB
Python
__author__ = 'dmy'
|
||
|
||
import math
|
||
import shutil
|
||
class Instantce:#仿真的实例
|
||
def __init__(self,instance,ipsoData):
|
||
self._ipsoData=ipsoData
|
||
self._genInstance(instance)
|
||
pass
|
||
def _genInstance(self,instance):
|
||
#序号 投入点 功率因数
|
||
#DGnode=DGnode.split(',')
|
||
#iPsodf=iPsoData(r'D:\Project\佛山项目\大学路西Ⅰ线_3447519_2013-10-29\大学路西Ⅰ线_3447519_2013-10-29_iPso_newFile.txt')
|
||
iPsodf=self._ipsoData
|
||
#22500
|
||
nodePow=iPsodf.GetNodePower()
|
||
#先统计一下总的负荷是多少
|
||
sumS=0
|
||
for n in nodePow:
|
||
P=float(n[6])
|
||
Q=float(n[7])
|
||
S=math.sqrt(P*P+Q*Q)
|
||
sumS+=S
|
||
#按系数重新分配负荷
|
||
for n in nodePow:
|
||
P=float(n[6])
|
||
Q=float(n[7])
|
||
if math.fabs(P)<1e-5 or math.fabs(Q)<1e-5:
|
||
continue
|
||
S=math.sqrt(P*P+Q*Q)
|
||
loadFactor=1#负荷比例
|
||
newS=S/sumS
|
||
newP=math.cos(math.atan(Q/P))*newS
|
||
newQ=math.sin(math.atan(Q/P))*newS
|
||
n[6]=newP
|
||
n[7]=newQ
|
||
self._ipsoData=iPsodf
|
||
return iPsodf
|
||
pass
|
||
def _copyGAMS(self,dirPath):
|
||
shutil.copyfile(r"D:\Project\佛山项目\Code\GAMS\xiangmu\try8.gms", dirPath+'/try8.gms')
|
||
shutil.copyfile(r"D:\Project\佛山项目\Code\GAMS\xiangmu\opf.gms", dirPath+'/opf.gms')
|
||
shutil.copyfile(r"D:\Project\佛山项目\Code\GAMS\xiangmu\main.m", dirPath+'/main.m')
|
||
shutil.copyfile(r"D:\Project\佛山项目\Code\GAMS\xiangmu\readOPF.m", dirPath+'/readOPF.m')
|
||
shutil.copyfile(r"D:\Project\佛山项目\Code\GAMS\xiangmu\FormY.m", dirPath+'/FormY.m')
|
||
shutil.copyfile(r"D:\Project\佛山项目\Code\GAMS\xiangmu\pfxiugai.m", dirPath+'/pfxiugai.m')
|
||
shutil.copyfile(r"D:\Project\佛山项目\Code\GAMS\xiangmu\FormSet.m", dirPath+'/FormSet.m')
|
||
shutil.copyfile(r"D:\Project\佛山项目\Code\GAMS\xiangmu\FormParameter.m", dirPath+'/FormParameter.m')
|
||
shutil.copyfile(r"D:\Project\佛山项目\Code\GAMS\xiangmu\doMain.m", dirPath+'/doMain.m')
|
||
shutil.copyfile(r"D:\Project\佛山项目\Code\GAMS\xiangmu\generateGDX.m", dirPath+'/generateGDX.m')
|
||
|
||
def ToDataFile(self,fileName):
|
||
newFilePath=fileName
|
||
iPsodf=self._ipsoData
|
||
with open(newFilePath,'w') as f:
|
||
for b in iPsodf.GetBasicInfo():
|
||
#f.write('%d,%d,%d,%.10f,%d\n'%(b[0],b[1],1,b[3],b[4]))
|
||
f.write('%d,%d,%d,%.10f,%d\n'%(b[0],b[1],1,b[3],1))
|
||
f.write('0\n')
|
||
for b in iPsodf.GetBalance():
|
||
f.write('%d,%d,%.10f\n'%(b[0],b[1],b[2]))
|
||
f.write('0\n')
|
||
for b in iPsodf.GetLine():
|
||
f.write('%d,%d,%d,%.10f,%.10f,%.10f,%.10f\n'%(b[0],b[1],b[2],b[3],b[4],b[5],b[6]))
|
||
f.write('0\n')
|
||
for b in iPsodf.GetTrans():
|
||
f.write('%d,%d,%d,%d,%.10f,%.10f,%.10f,%.10f,%.10f,%d\n'%(
|
||
b[0],b[1],b[2],b[3],b[4],b[5],b[6],b[7],b[8],b[9]))
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
for b in iPsodf.GetGround():
|
||
f.write('%d,%d,%.10f,%.10f\n'%(b[0],b[1],b[2],b[3]))
|
||
f.write('0\n')
|
||
for b in iPsodf.GetCap():
|
||
f.write('%d,%d,%d,%.10f,%.10f,%d,%d,%d\n'%(b[0],b[1],b[2],b[3],b[4],b[5],b[6],b[7]))
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
for b in iPsodf.GetNodePower():
|
||
f.write('%d,%d,%d,%d,%.10f,%.10f,%.10f,%.10f,%f,%.10f\n'%(
|
||
b[0],b[1],b[2],b[3],b[4],b[5],b[6],b[7],b[8],b[9]))
|
||
f.write('0\n')
|
||
for b in iPsodf.GetGen():
|
||
f.write('%d,%d,%d,%d,%.10f,%.10f,%.10f,%.10f,%.10f,%.10f\n'%(
|
||
b[0],b[1],b[2],b[3],b[4],b[5],b[6],b[7],b[8],b[9]))
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
for b in iPsodf.GetReactiveGen():
|
||
f.write('%d,%d,%d,%d,%.10f,%.10f,%.10f\n'%(b[0],b[1],b[2],b[3],b[4],b[5],b[6]))
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
#for b in iPsodf.GetVoltLimit():
|
||
#f.write('%d,%d,%.10f,%.10f\n'%(b[0],b[1],b[2],b[3]))
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
f.write('0\n')
|
||
copyGAMSPath=os.path.dirname(newFilePath)
|
||
self._copyGAMS(copyGAMSPath) |