原文件切片,选出主接触界面的坐标和力
path = "C:/Users/Administrator/Desktop/"
file1 = "ncforc.txt"
with open(path+file1, "r") as fp:
data_origin = fp.readlines()
data_master = []
'''
for i in range(11):
data_master.append(data_origin[i])
'''
j = 20938
data_master.append(str(0)+"\n""forces\n")
k=1
while(j<len(data_origin)):
for i in range(9540):
if (j+i+4)<len(data_origin):
data_master.append(data_origin[j+i+4])
data_master.append(str(k)+"\n"+"forces\n")
j = j + 30474
k=k+1
file2 = "ncforce_maaster.txt"
with open(path+file2,"w") as fp2:
for i in range(len(data_master)):
fp2.write(data_master[i])
fp.close()
主界面的节点导出的是整个表面的,这段程序用来选出来y方向的力和节点坐标,下一个程序我们通过判断条件选出来同脚接触的那个面的节点的坐标和力
import numpy as