python仿真报告自动化——excite TD齿轮角加速度级计算

python仿真报告自动化------excite TD齿轮角加速度级计算

  • [1 、问题-燃油泵相位优化](#1 、问题-燃油泵相位优化)
  • [2、难点-excite TD结果文件的提取](#2、难点-excite TD结果文件的提取)
  • 3、代码

1 、问题-燃油泵相位优化

用excite TD对齿轮系进行仿真,模拟不同燃油泵相位对齿轮传动振动的影响,用齿轮角加速度级作为评价指标,需要计算大量的数据,用IMPRESS Chart手动处理,太麻烦。

2、难点-excite TD结果文件的提取

excite TD结果文件.GID在用python读取时,有大量的空格,需合并连续空格,再用空格分割提取齿轮角角速度FFT的数据。
用正则化表达式处理:

line=re.sub(r'\s+',' ',line)

3、代码

计算各齿轮的角加速度级,绘图,并将数据写入excel中。

python 复制代码
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import re

plt.rcParams['font.sans-serif']=['SimHei'] #用来正常显示中文标签
plt.rcParams['axes.unicode_minus']=False #用来正常显示负号

file_excel_fullname=r'E:\NVH_sim_\ppt\ai_gear.xlsx'
cirv_name=['G-Crankshaft','G-SteppedidleBig','G-SteppedidleSmall','G-Camidle2','G-Camidle1','G-Intake','G-Exhaust',\
           'G-FIP','G-AirPumpidle','G-AirPump','G-OilPumpidle','G-OilPump']
file_name=['SHGE_164_F.GID','SHGE_191_F.GID','SHGE_201_F.GID','SHGE_212_F.GID','SHGE_215_F.GID','SHGE_219_F.GID','SHGE_218_F.GID',\
           'SHGE_193_F.GID','SHGE_169_F.GID','SHGE_171_F.GID','SHGE_180_F.GID','SHGE_183_F.GID']
NUM=12
fi=np.zeros((226,NUM))
ai=np.zeros((226,NUM))
a_db=np.zeros(NUM)
for j in range(NUM):
    file_path = r'E:\NVH_sim_\excite_td\TD-1600rpm\results'+'\\'+file_name[j]
    with open(file_path) as f:
        content=f.read()
        #拆分内容
        lines=content.split('\n')
        for i,line in enumerate(lines[149:375]):
            #利用正则表达式将连续的空格合并为一个空格
            line=re.sub(r'\s+',' ',line)
            item=line.split(' ')
            fi[i,j]=float(item[1])
            ai[i,j]=float(item[32])
        a_db[j]=10*np.log10(np.sum(ai[:,j]*ai[:,j]))
for j in range(NUM):
    plt.plot(fi[:,j],10*np.log10(ai[:,j]*ai[:,j]),label='$%s:%5.2f$'%(cirv_name[j],a_db[j]))
plt.legend()
plt.xlabel('fi/Hz')
plt.ylabel('ai/Db')
plt.show()
print(a_db)

# c创建一个DataFrame数据类型,并写入表头
dataFrame = pd.DataFrame(columns=("fi", 'G-Crankshaft','G-SteppedidleBig','G-SteppedidleSmall','G-Camidle2','G-Camidle1','G-Intake','G-Exhaust',\
                                  'G-FIP','G-AirPumpidle','G-AirPump','G-OilPumpidle','G-OilPump'))
#添加两行数据
for i in range(226):
    dataFrame.loc[i+1] = [fi[i,0], ai[i,0], ai[i,1],ai[i,2], ai[i,3],ai[i,4], ai[i,5],ai[i,6], \
                          ai[i,7],ai[i,8],ai[i,9], ai[i,10],ai[i,11]]
dataFrame.loc[i+2]=['overall',a_db[0], a_db[1],a_db[2], a_db[3], a_db[4],a_db[5], a_db[6],\
                    a_db[7], a_db[8],a_db[9], a_db[10], a_db[11]]
# 写入excel文件
#dataFrame.values[:,0]=fi[:,0]
dataFrame.to_excel(file_excel_fullname)
复制代码
相关推荐
小安运维日记几秒前
RHCA - DO374 | Day01:使用红帽Ansible自动化平台开发剧本
运维·服务器·云原生·自动化·云计算·ansible
秋空樱雨2 分钟前
C++入门
开发语言·c++
咬_咬22 分钟前
C++仿mudo库高并发服务器项目:Buffer模块
服务器·开发语言·c++·缓冲区·buffer·muduo库
江公望29 分钟前
Qt qmlplugindump浅谈
开发语言·qt·qml
曦樂~30 分钟前
【Qt】文件操作/事件--mainwindow做编辑器
开发语言·qt
zzzyulin35 分钟前
huggingface transformers调试问题--加载本地路径模型时pdb断点消失
python·transformer
敲代码的瓦龙36 分钟前
西邮移动应用开发实验室2025年二面题解
开发语言·c++·算法
laocooon52385788638 分钟前
一个适合新手的训练C题
c语言·开发语言
教练、我想打篮球38 分钟前
12 pyflink 的一个基础使用, 以及环境相关
python·flink·pyflink
飞翔的佩奇1 小时前
【完整源码+数据集+部署教程】【天线&运输】直升机战机类型识别目标检测系统源码&数据集全套:改进yolo11-CSP-EDLAN
前端·python·yolo·计算机视觉·数据集·yolo11·直升机战机类型识别目标检测系统