文档章节

Python制作NTF传递函数工况文件和后处理文件

o
 osc_w9s1w4o0
发布于 2019/03/31 23:39
字数 1885
阅读 3
收藏 0

精选30+云产品,助力企业轻松上云!>>>

     摘要:在平时工作中,TB车身的传递函数分析,涉及到大量重复行的工作,费时费力。在学习python基础后,希望通过代码解决这部分重复工作。基础入门级操作,但是能够解决很大一部分工作内容。日后,待python学习达到一定水平后,在寻找更优的解决方案。

     NVH频响,接附点一般在30个左右,分x,y,z三个方向,激励工况在90个上下,每个工况对应DAREA,RLOAD,DLOAD;在实际项目中,习惯运用工况文本文件进行提交运算。NTF计算,90个工况,对应90个响应;VTF计算,90个工况对应270个响应;IPI计算90个工况对应270个响应。前处理工况加载事件费事费力的工作,结果处理更是重复繁杂,而且容易出错。基于这样的情况,以下运用python基础技巧,进行解决。

 

一.文件管理

1.项目阶段TG1,TG2,TG3;   

  1.1每一阶段(如TG1)分别用三个文件夹来存储计算模型和结果

    1.TG1_base:文件夹存放该论分析的基础模型;备注好数据日期,数据输入缺失项,以及商定缺失数据处理方案;

    2.TG1_Engineering_scheme:文件夹存放工程技术方案(由你向设计部分提供概念优化方案,并且认定为可实施的工程方案,设计部分修改后的数据,逐一替换到该模型之中,每一次替换保存一个hm),在该文件夹中创建以日期命名的文件夹,用于存储设计部分反馈数据并且替换好的模型,并且以Excel备注好反馈的简明信息,并且每版数据替换后保存好hm文件,以免connector丢失;如有必要,在对应的日期的工程反馈数据模型中,进行问题诊断。

    3.TG1_conceptual_scheme:文件夹存放概念优化方案(case1代表第一个概念方案,case2代表第二个概念方案········;一旦方案被认定为可工程化方案,移交设计部分进行设计画数据),并且以Excel备注好,方案基于那一版工程数据,方案的截图;路径及文件夹命名最好不要用汉字,用字母、数字、下划线;】 

  1.2运用python处理的文件为工况文件fem及结果处理文件mvw 

    

 

 

 

二.工况文件fem及结果处理文件mvw

  2.1.NTF计算工况文件fem

         30个这样的组合(甚至更多,后续默认30个接附点)

                           

 

 

  2.NTF后处理文件mvw

    90个这样的图片,一个工况subcase对应一个page包含3条曲线,红色代表基础模型计算结果,绿色代表可工程化的优化方案,蓝色虚线代表概念优化方案。

 

三.运用python处理重复工作(工况文件fem,结果处理文件mvw)

  1.计算工况文件及结果处理文件mvw的制作

    1.1subcase以及DAREA,RLOAD,DLOAD的编写

      

    1.2所用到的python的命令

      函数:def    函数定义好后,不调用就不会执行,当需要制作计算工况文件时,调用函数 header_sub();当需要制作后处理MVW文件时,调用NTF_result()

 

 

 
#coding=utf-8
#Version:python3.7.0
#tools:Pycharm 2017.3.2
##先通过hypergraph按自己要求制作出一张page,保存mvw格式文件作为需要重复的输入文件
mvw_onepage_input = "D:/NVH/NVH/NVH_scripts/TB_NTF/1.header_result/1.input_mvw/lne3_target3.mvw"
##python通过接收输入的文件,逐行读取,修改关键参数,在写入到后处理文件中
mvw_NTF_output = "D:/NVH/NVH/NVH_scripts/TB_NTF/1.header_result/2.output_mvw/Base_Engineering_conceptual.mvw"
##计算工况文件,内含在python代码中,直接输出即可。其他参数比较单一,直接复制粘贴即可
header_NTF_output = "D:/NVH/NVH/NVH_scripts/TB_NTF/1.header_result/3.output_header/NTF_header.fem"
##需要加载工况的接附点,制作成列表,方便代码调用
input_point = [101,102,103,201,202,203,204,205,206,301,302,303,304,305,306,307,308,309,310,311,312,401,402,403,404,405,406,501, 502,503,504,505,506,601,602,603,604]
##在工况文件中一个接附点(如101),对应3个方向的激励,定义方向列表
direc = ['X','Y','Z']
##定义后处理mvw文件制作,替换onepage中的所有SUBCASE 1 = 101:+X,p1w1,Page 1,page1,WindowIDs(2),然后逐行左对齐写入Base_Engineering_conceptual.mvw文件中
def NTF_result():
with open(mvw_NTF_output, 'a') as NTF_object: #以追加形式打开文件mvw_NTF_output,如果文件不存在,创建文件
NTF_object.write("{ safe_quotes_on }\n") #写入文本
NTF_object.write('{ PLOT_FILE_1 = "TG1_base.op2"}\n')
NTF_object.write('{ PLOT_FILE_2 = "TG1_Engineering_scheme.op2" }\n')
NTF_object.write('{ PLOT_FILE_3 = "TG1_conceptual_scheme.op2" }\n')
NTF_object.write('*Id("HyperWorks", "12.*")\n')
NTF_object.write('*BeginPalette()\n')
NTF_object.write('*EndPalette()\n')
i = 0
while i < len(input_point): #开始处理重复page
j = 0
while j < 3:
with open(mvw_onepage_input) as file_object: #打开事先调试好的一页page文件mvw_onepage_input
for line in file_object:
sub_id = 3 * i + j + 1
if sub_id < 10 : #此处是因为前任制作的工况文件的sub是1-9 19-120,所以此处用if语句跳过10-18
s1 = line.replace('SUBCASE 1 = 101:+X',
'SUBCASE {n1} = {n2}:+{n3}'.format(n1=3 * i + j + 1, n2=input_point[i],
n3=direc[j]))
s2 = s1.replace('p1w1', 'p{n1}w1'.format(n1=3 * i + j + 1))
s3 = s2.replace('Page 1', 'Page {n1}'.format(n1=3 * i + j + 1))
s4 = s3.replace('WindowIDs(2)', 'WindowIDs({n1})'.format(n1=3 * i + j + 2))
s5 = s4.replace('page1', 'page{n1}'.format(n1=3 * i + j + 1))
with open(mvw_NTF_output, 'a') as NTF_object:
NTF_object.write(s5.rstrip())
NTF_object.write('\n')
else :
sub_id = 3 * i + j + 10
s1 = line.replace('SUBCASE 1 = 101:+X',
'SUBCASE {n1} = {n2}:+{n3}'.format(n1=sub_id, n2=input_point[i],
n3=direc[j]))
s2 = s1.replace('p1', 'p{n1}'.format(n1=3 * i + j + 1))
s3 = s2.replace('Page 1', 'Page {n1}'.format(n1=3 * i + j + 1))
s4 = s3.replace('WindowIDs(2)', 'WindowIDs({n1})'.format(n1=3 * i + j + 2))
s5 = s4.replace('page1', 'page{n1}'.format(n1=3 * i + j + 1))
with open(mvw_NTF_output, 'a') as NTF_object:
NTF_object.write(s5.rstrip())
NTF_object.write('\n')
j += 1
i += 1

def NTF_sub(): #定义工况文件,注意int数据类型转变为str
i = 0
while i < len(input_point): #写入sucase
j = 0
while j < 3:
with open(header_NTF_output, 'a') as header_object:
header_object.write('SUBCASE {n1}\n'.format(n1=3 * i + j + 1))
header_object.write(' LABEL {n1}:+{n2}\n'.format(n1=input_point[i],n2=direc[j]))
header_object.write(' DLOAD = {n1}\n'.format(n1=input_point[i]*10+j+1))
j += 1
i += 1
a = 0
while a < len(input_point): #写入DAREA,RLOAD,DLOAD
b = 0
while b < 3:
with open(header_NTF_output, 'a') as header_object:
dload_id = input_point[a]*10 + b+1
rLoad_id = 20000+input_point[a]*10 + b+1
darea_id = 10000+input_point[a]*10 + b+1
point = '#'+str(input_point[a])
xyz = b+1
header_object.write('\n')
header_object.write(point.ljust(8))
header_object.write('\n')
header_object.write('DLOAD'.ljust(8))
header_object.write(str(dload_id).rjust(8))
header_object.write('1.0'.rjust(8))
header_object.write('1.0'.rjust(8))
header_object.write(str(rLoad_id).rjust(8))
header_object.write('\n')
header_object.write('RLOAD1'.ljust(8))
header_object.write(str(rLoad_id).rjust(8))
header_object.write(str(darea_id).rjust(8))
header_object.write('3'.rjust(24))
header_object.write('0'.rjust(8))
header_object.write('LOAD'.rjust(8))
header_object.write('\n')
header_object.write('DAREA'.ljust(8))
header_object.write(str(darea_id).rjust(8))
header_object.write(str(input_point[a]).rjust(8))
header_object.write(str(xyz).rjust(8))
header_object.write('\n')
b += 1
a += 1
NTF_sub() #工况文件值写入了sub DAREA RLOAD DLOAD,其余相关设置复制粘贴即可
NTF_result() #调用定义的函数,不调用函数不执行


 四.运用

  1.NTF_header.fem

    不重复的部分,python代码并没有写入;相关参数直接用复制以往计算文件,用NTF_header中的sub替换即可

  2.Base_Engineering_conceptual.mvw

    用UE打开,对应结果文件的路径及名称替换后就看可以用hypergraph打开,就出现三条曲线对比的图

    { PLOT_FILE_1 = "`````/TG1/1.TG1_NTF_base/TG1_base.op2"}
    { PLOT_FILE_2 = "`````/TG1/2.TG1_NTF_engineering_scheme/20190325/TG1_Engineering_scheme.op2" }
    { PLOT_FILE_3 = "`````/TG1/3.TG1_NTF_conceptual_scheme/case1/TG1_conceptual_scheme.op2" }

总结:该过程只能解决整个分析过程的一些片段,不能系统的解决工程问题。本人也在不断学习改进中,希望通过解决实际工作中的问题,来学习python。没错,是为了学习python,把工作问题拿来练习。

如果各位大神有更好的解决方案,欢迎邮件联系:liuhuacai@outlook.com

 

o
粉丝 0
博文 500
码字总数 0
作品 0
私信 提问
加载中
请先登录后再评论。
Python异常处理(基础详解)

本文是python异常的基础知识,欢迎阅读,一起进步 :人生苦短-我学python 文章目录 一.异常简介 异常是指在语法正确的前提下,程序运行时报错就是异常。 当Python脚本发生异常时我们需要捕获处...

辰兮要努力
04/06
0
0
Django框架(四):视图、模板

视图 后台管理页面做好了,接下来就要做公共访问的页面了。当我们刚刚在浏览器中输入http://127.0.0.1:8000/admin/之后,浏览器显示出了后台管理的登录页面,那么服务器是如何找到这个页面并...

osc_a6lww194
01/17
1
0
使用python开发android应用程序:第四节 使用WebView制作程序界面

介绍 SL4A r1 版本引入了内置的HTML编译器和控制WebView的api。 UiFacade模块也提供了一系列直接操作对话框和菜单的方法。想了解更多信息,请参考Ui Help部分的介绍。 使用 WebView API 中的...

无鸯
2011/09/30
2.4K
0
跟上时代的步伐,学一波Python(二)

相关阅读:跟上时代的步伐,学一波Python(一) 6. 字典 6.1 一个简单的字典 字典alien_0 存储了外星人的颜色和点数。使用两条print 语句来访问并打印这些信息,如下所示: 6.2 使用字典 在Pyt...

朱敏_ITer
2018/01/17
0
0
2018年9月Python计算机二级复习攻略

考试方式:   上机考试,考试时长 120 分钟,满分 100 分。   1. 题型及分值   单项选择题 40 分( 含公共基础知识部分 10 分)。   操作题 60 分( 包括基本编程题和综合编程题)。   ...

安全牛课堂
2018/09/12
22
0

没有更多内容

加载失败,请刷新页面

加载更多

自从尝了 Rust,Java 突然不香了

Rust 是软件行业中相对而言比较新的一门编程语言,如果从语法上来比较,该语言与 C++ 其实非常类似,但从另一方面而言,Rust 能更高效地提供许多功能来保证性能和安全。而且,Rust 还能在无需...

osc_k3vwonkw
2分钟前
0
0
Java 高级 面试题 及 参考答案

一、面试题基础总结 1、 JVM结构原理、GC工作机制详解 答:具体参照:JVM结构、GC工作机制详解 ,说到GC,记住两点:1、GC是负责回收所有无任何引用对象的内存空间。 注意:垃圾回收回收的是无...

FH-Admin
3分钟前
0
0
机器学习中的AUC-ROC曲线

作者|ANIRUDDHA BHANDARI 编译|VK 来源|Analytics Vidhya AUC-ROC曲线 你已经建立了你的机器学习模型-那么接下来呢?你需要对它进行评估,并验证它有多好(或有多坏),这样你就可以决定是否...

osc_bg8v9gvf
4分钟前
0
0
音视频(消息)应用场景 :连麦交友例子

实现一个小例子: 效果类似唱吧APP里的 连麦交友功能,音视频,IM 及音视频 SDK参考融云服务商。 没有印象的可以搜索 ’连麦’ 关键字在 应用商店下载一款 连麦的软件 体验下 业务方面的需求...

T型人才追梦者
6分钟前
4
0
逛淘宝天猫想到SSO单点登录

我的原文地址:https://mp.weixin.qq.com/s/77xukPDlgkKnYpwu4LrqaA

osc_yy65eb2q
6分钟前
11
0

没有更多内容

加载失败,请刷新页面

加载更多

返回顶部
顶部