data_analyze_platform/LIB/MIDDLE/CellStateEstimation/BatSafetyWarning/main.py

import pandas as pd
import pymysql
from LIB.BACKEND import DBManager, Log
from apscheduler.schedulers.blocking import BlockingScheduler
import time, datetime
from LIB.MIDDLE.CellStateEstimation.Common.V1_0_1 import DBDownload
from LIB.MIDDLE.CellStateEstimation.Common.V1_0_1 import log
import CBMSBatInterShort
import CBMSSafetyWarning

#电池热安全预警核心算法函数......................................................................................................................
def saftywarning_cal():
    global SNnums
    global df_warning_ram
    global df_warning_ram1
    global df_warning_ram2
    global df_warning_ram3
    global df_lfp_ram
    global now_time
    
    start=time.time()
    # now_time=datetime.datetime.now()
    # start_time=now_time-datetime.timedelta(hours=12)
    # start_time1=now_time-datetime.timedelta(days=7)
    # start_time2=now_time-datetime.timedelta(days=90)
    # start_time=start_time.strftime('%Y-%m-%d %H:%M:%S')
    # start_time1=start_time1.strftime('%Y-%m-%d %H:%M:%S')
    # end_time=now_time.strftime('%Y-%m-%d %H:%M:%S')

    now_time=now_time+datetime.timedelta(hours=12)
    start_time=now_time-datetime.timedelta(hours=12)
    start_time1=now_time-datetime.timedelta(days=7)
    start_time2=now_time-datetime.timedelta(days=90)
    start_time=start_time.strftime('%Y-%m-%d %H:%M:%S')
    start_time1=start_time1.strftime('%Y-%m-%d %H:%M:%S')
    end_time=now_time.strftime('%Y-%m-%d %H:%M:%S')

    #数据库配置
    host='rm-bp10j10qy42bzy0q77o.mysql.rds.aliyuncs.com'
    port=3306
    user='qx_read'
    password='Qx@123456'

    #读取故障结果库中code==110且end_time='0000-00-00 00:00:00'...............................
    db='safety_platform'
    mysql = pymysql.connect (host=host, user=user, password=password, port=port, database=db)
    cursor = mysql.cursor()
    param='start_time,end_time,product_id,code,level,info,advice'
    tablename='all_fault_info'
    sql =  "select %s from %s where code=110 and end_time='0000-00-00 00:00:00'" %(param,tablename)
    cursor.execute(sql)
    res = cursor.fetchall()
    df_fault_ram= pd.DataFrame(res,columns=param.split(','))
    cursor.close()
    mysql.close()

    for sn in SNnums:
        if 'PK500' in sn:
            celltype=1 #6040三元电芯
        elif 'PK502' in sn:
            celltype=2 #4840三元电芯
        elif 'K504B' in sn:
            celltype=99    #60ah林磷酸铁锂电芯
        elif 'MGMLXN750' in sn:
            celltype=3 #力信50ah三元电芯
        elif 'MGMCLN750' or 'UD' in sn: 
            celltype=4 #CATL 50ah三元电芯
        else:
            print('SN：{}，未找到对应电池类型！！！'.format(sn))
            continue
            # sys.exit()

        #读取原始数据库数据........................................................................................................................................................
        dbManager = DBManager.DBManager()
        df_data = dbManager.get_data(sn=sn, start_time=start_time, end_time=end_time, data_groups=['bms'])
        df_bms = df_data['bms']
        # tim=start_time.replace(':','_')
        # df_bms.to_csv(r'D:\Platform\platform_python\data_analyze_platform\USER\spf\01qixiang\98Download\\'+sn+'_BMS_'+tim+'.csv',encoding='GB18030')

        #读取结果数据库数据........................................................................................................................................................
        db='qx_cas'
        mode=1
        tablename='cellstateestimation_soh'
        DBRead=DBDownload.DBDownload(host, port, db, user, password,mode)
        with DBRead as DBRead:
            df_soh=DBRead.getdata('time_st,sn,soh,cellsoh', tablename=tablename, sn=sn, timename='time_sp', st=start_time, sp=end_time)

        if not df_bms.empty:
            #电池内短路计算................................................................................................................................................................
            df_warning_ram_sn=df_warning_ram[df_warning_ram['sn']==sn]
            df_warning_ram_sn1=df_warning_ram1[df_warning_ram1['sn']==sn]
            df_warning_ram_sn2=df_warning_ram2[df_warning_ram2['sn']==sn]
            df_warning_ram_sn3=df_warning_ram3[df_warning_ram3['sn']==sn]
            df_warning_ram_sn.reset_index(inplace=True,drop=True)     #重置索引
            df_warning_ram_sn1.reset_index(inplace=True,drop=True)     #重置索引
            df_warning_ram_sn2.reset_index(inplace=True,drop=True)     #重置索引
            df_warning_ram_sn3.reset_index(inplace=True,drop=True)     #重置索引
            if celltype>50:
                df_lfp_ram=df_lfp_ram.reindex(columns=df_bms.columns.tolist()+['sn'])
                df_lfp_ram_sn=df_lfp_ram[df_lfp_ram['sn']==sn]
                df_lfp_ram_sn.reset_index(inplace=True,drop=True)     #重置索引
            else:
                df_lfp_ram_sn=pd.DataFrame()

            BatShort=CBMSBatInterShort.BatInterShort(sn,celltype,df_bms,df_soh,df_warning_ram_sn,df_warning_ram_sn1,df_warning_ram_sn2,df_warning_ram_sn3,df_lfp_ram_sn)
            df_short_res, df_ram_res, df_ram_res1, df_ram_res2, df_ram_res3, df_ram_res4=BatShort.intershort() 

            if not df_short_res.empty:
                with open(r'D:\Platform\platform_python\data_analyze_platform\USER\spf\01qixiang\07BatSafetyWarning\内短路.txt','a') as file:
                    file.write(str(df_short_res)+'\n')

            #内短路ram处理.........................................................
            df_warning_ram=df_warning_ram.drop(df_warning_ram[df_warning_ram.sn==sn].index)
            df_warning_ram1=df_warning_ram1.drop(df_warning_ram1[df_warning_ram1.sn==sn].index)
            df_warning_ram2=df_warning_ram2.drop(df_warning_ram2[df_warning_ram2.sn==sn].index)
            df_warning_ram3=df_warning_ram3.drop(df_warning_ram3[df_warning_ram3.sn==sn].index)

            df_warning_ram=pd.concat([df_warning_ram,df_ram_res],ignore_index=True)
            df_warning_ram1=pd.concat([df_warning_ram1,df_ram_res1],ignore_index=True)
            df_warning_ram2=pd.concat([df_warning_ram2,df_ram_res2],ignore_index=True)
            df_warning_ram3=pd.concat([df_warning_ram3,df_ram_res3],ignore_index=True)
            
            if celltype>50:
                df_lfp_ram=df_lfp_ram.drop(df_lfp_ram[df_lfp_ram.sn==sn].index)
                df_lfp_ram=pd.concat([df_lfp_ram,df_ram_res4],ignore_index=True)
            

        #电池热安全预警..............................................................................................................................................................
        #读取内短路、析锂和一致性结果数据库数据
        db='qx_cas'
        mode=2
        tablename1='cellstateestimation_intershort'
        tablename2=''   #析锂表单
        tablename3='cellstateestimation_uniform_socvoltdiff'
        DBRead=DBDownload.DBDownload(host, port, db, user, password,mode)
        with DBRead as DBRead:
            df_short=DBRead.getdata('time_sp,sn,short_current', tablename=tablename1, sn=sn, timename='time_sp', st=start_time1, sp=end_time)
            # df_liplated=DBRead.getdata('time,sn,liplated,lipltd_amount', tablename=tablename2, sn=sn, timename='time_sp', st=start_time2, sp=end_time)
            df_uniform=DBRead.getdata('time,sn,cellsoc_diff,cellmin_num', tablename=tablename3, sn=sn, timename='time', st=start_time1, sp=end_time)
            df_liplated=pd.DataFrame()
        
        #获取sn的故障RAM
        df_fault_ram_sn=df_fault_ram[df_fault_ram['product_id']==sn]
        
        #热安全预警
        BatWarning=CBMSSafetyWarning.SafetyWarning(sn,celltype,df_short,df_liplated,df_uniform,df_fault_ram_sn)
        df_warning_res=BatWarning.diag()

        df_warning_end=df_warning_res[df_warning_res['end_time'] != '0000-00-00 00:00:00']
        if not df_warning_end.empty:    #df_warning_end历史故障筛选并更改数据库故障结束时间
            pass   
        elif df_fault_ram_sn.empty and (not df_warning_res.empty):
            with open(r'D:\Platform\platform_python\data_analyze_platform\USER\spf\01qixiang\07BatSafetyWarning\热失控预警.txt','a') as file:
                file.write(str(tuple(df_warning_res.iloc[-1]))+'\n')
        else:
            pass

        
        print(df_warning_ram)
        print(df_warning_ram1)
        print(df_warning_ram2)
        print(df_warning_ram3)
        print(df_lfp_ram)

        end=time.time()
        print(end-start)
        # print(df_soh)

#...............................................主函数起定时作用.......................................................................................................................
if __name__ == "__main__":
    
    excelpath=r'D:\Platform\platform_python\data_analyze_platform\USER\spf\01qixiang\sn-20210903.xlsx'
    SNdata_6060 = pd.read_excel(excelpath, sheet_name='科易6060')
    SNdata_6040 = pd.read_excel(excelpath, sheet_name='科易6040')
    SNdata_4840 = pd.read_excel(excelpath, sheet_name='科易4840')
    SNdata_L7255 = pd.read_excel(excelpath, sheet_name='格林美-力信7255')
    SNdata_C7255 = pd.read_excel(excelpath, sheet_name='格林美-CATL7255')
    SNdata_U7255 = pd.read_excel(excelpath, sheet_name='优旦7255')
    SNnums_6060=SNdata_6060['SN号'].tolist()
    SNnums_6040=SNdata_6040['SN号'].tolist()
    SNnums_4840=SNdata_4840['SN号'].tolist()
    SNnums_L7255=SNdata_L7255['SN号'].tolist()
    SNnums_C7255=SNdata_C7255['SN号'].tolist()
    SNnums_U7255=SNdata_U7255['SN号'].tolist()
    SNnums=SNnums_L7255 + SNnums_C7255 + SNnums_U7255, SNnums_6040 + SNnums_4840 + SNnums_6060
    SNnums=['MGMCLN750N215N099','PK504B10100004425']
    
    mylog=log.Mylog('log_warning.txt','error')
    mylog.logcfg()

    #............................模块运行前，先读取数据库中所有结束时间为0的数据，需要从数据库中读取...................................
    df_warning_ram=pd.DataFrame(columns=['sn','time','deltsoc','cellsoc'])
    df_warning_ram1=pd.DataFrame(columns=['sn','time1','deltsoc1'])
    df_warning_ram2=pd.DataFrame(columns=['sn','time2','deltAs2'])
    df_warning_ram3=pd.DataFrame(columns=['sn','standingtime','standingtime1'])
    df_lfp_ram=pd.DataFrame()

    now_time='2021-10-16 00:00:00'
    now_time=datetime.datetime.strptime(now_time,'%Y-%m-%d %H:%M:%S')

    #定时任务.......................................................................................................................................................................
    scheduler = BlockingScheduler()
    scheduler.add_job(saftywarning_cal, 'interval', seconds=5)

    try:  
        scheduler.start()
    except Exception as e:
        scheduler.shutdown()
        print(repr(e))
        mylog.logopt(e)