Python学习与实践（1）ansys载荷文件转换

之前的文章里讲了设计波和aqwa的载荷传递，接下来介绍夜歌Python在载荷文件转换中的应用。

 众所周知，aqwa经过载荷传递后得到的载荷文件是以dat结尾的文件，而且里面只有波浪载荷文件、加速地、加速度作用的坐标和角加速度。  

 虽然，经典ANSYS可以直接读取，但就算计算强度工况也要上百个，计算时域疲劳上百万个都有可能，都手动读取是不现实的。  

 所以，就要把这些.dat文件转换为ANSYS计算载荷步的.sxxx文件。  

 下面的代码就是Python用来做载荷文件转换的一个代码的例子。  

# coding: utf8

import os

i=1

while i<=32:

    n=i+1000

    f = open('OC5_aqld'+str(10000+i)+'.dat', 'r')

    content = f.read()  # 读取文件内容

    modified_content = content.replace('SOLVE', '')# 替换 "SOLVE" 为空字符串

    f_new = open('b.dat', 'w')

    f_new.write(modified_content)    # 写入原文件内容

    f_new.write('\r\nF,35859,FX,2648400\r\nF,35859,FY,2743100\r\nF,35859,FZ,-20139800\r\nF,35859,MX,-339363900\r\nF,35859,MY,245666200\r\nF,35859,MZ,8300600') 

    f_new.write('\r\nd,111066, , , , , ,UX,UY,UZ, , ,'

                '\r\nd,55533  , , , , , ,UX,UY,UZ, , ,  '

                '\r\nd,166599, , , , , ,UX,UY,UZ, , ,  '

                '\r\nd,108188, , , , , ,,UX,UY,UZ, , ,'

                '\r\nd,52655, , , , , ,,UX,UY,UZ, , ,'

                '\r\nd,163721, , , , , ,,UX,UY,UZ, , , '

                '\r\nd,101630, , , , , ,,,UX,UY,UZ, , ,'

                '\r\nd,46097, , , , , ,,,UX,UY,UZ, , , '

                '\r\nd,157163, , , , , ,,,UX,UY,UZ, , ,')  # 写入边界条件

    f.close()

    f_new.close()

    os.remove('SPAR_aqld'+str(10000+i)+'.dat')  # 移除老文件

    os.rename('b.dat', 'SPAR.s'+str(int(n))) # 新文件重新命名

i=i+1

一共循环的32次也就处理了32个载荷文件。

为了防止过早的求解，删除了SOLVE命令，加入了一组载荷，并加入了边界条件。