Source Code

Auto3mmRadius2.py

@@ -0,0 +1,83 @@
+import math
+import os
+
+def read_data(file_name):
+    with open(file_name, 'r') as file:
+        lines = file.readlines()
+        data = []
+        for line in lines:
+            node_data = [float(value) for value in line.split()]
+            data.append(node_data)
+    return data
+
+def euclidean_distance(coord1, coord2):
+    return math.sqrt(sum([(a - b)**2 for a, b in zip(coord1, coord2)]))
+
+def find_nodes_within_distance(data, reference_node, distance):
+    nodes_within_distance = []
+    ref_node_coord = reference_node[1:4]
+    for node in data:
+        node_coord = node[1:4]
+        if euclidean_distance(ref_node_coord, node_coord) <= distance:
+            nodes_within_distance.append(node)
+    return nodes_within_distance
+
+def process_data(data):
+    distance = 3
+    average_min_principal_list = []
+    nodes_within_distance_list = []
+
+    for specified_node in data:
+        nodes_within_distance = find_nodes_within_distance(data, specified_node, distance)
+        nodes_within_distance_list.append(nodes_within_distance)
+
+        total_min_principal = sum(node[4] for node in nodes_within_distance)
+        average_min_principal = total_min_principal / len(nodes_within_distance)
+        average_min_principal_list.append(average_min_principal)
+
+    return average_min_principal_list, nodes_within_distance_list
+
+def write_results_to_file(file_name, data, average_min_principal_list):
+    MFS_coefficient = 500 * 0.0104
+    with open(file_name, 'w') as file:
+        for i, avg_min_principal in enumerate(average_min_principal_list):
+            MFS = MFS_coefficient / avg_min_principal
+            file.write(f"Node, {int(data[i][0])}, Average E. Min. Principal for 3mm Radius, {avg_min_principal:.6f},     MFS, {MFS:.6f}\n")
+
+def write_nodes_within_distance_to_file(file_name, data, nodes_within_distance_list):
+    with open(file_name, 'w') as file:
+        MFS_coefficient = 500 * 0.0104
+        for i, nodes_within_distance in enumerate(nodes_within_distance_list):
+            center_node = data[i]
+            file.write(f"Center Node: {int(center_node[0])}\n")
+            total_min_principal = 0
+
+            for node in nodes_within_distance:
+                line = "{:12d}{:16.6f}{:16.6f}{:16.6f}{:16.6f}\n".format(int(node[0]), node[1], node[2], node[3], node[4])
+                file.write(line)
+                total_min_principal += node[4]
+
+            if len(nodes_within_distance) > 0:
+                average_min_principal = total_min_principal / len(nodes_within_distance)
+                MFS = MFS_coefficient / average_min_principal
+                file.write(f"Average E. Min. Principal: {average_min_principal:16.6f},\nMFS: {MFS:.6f}\n\n")
+
+import os
+import math
+
+def process_files_auto_3mm_radius(input_folder, output_folder_1, output_folder_2):
+    # Iterate over the 28 input files
+    for file_number in range(1, 29):
+        input_file = os.path.join(input_folder, f"{file_number}_averaged_combined.txt")
+        data = read_data(input_file)
+
+        # Process the data and find the average E. Min. Principal for each node and the nodes within 3mm distance
+        average_min_principal_list, nodes_within_distance_list = process_data(data)
+
+        # Write the average E. Min. Principal results to a text file
+        output_file_1 = os.path.join(output_folder_1, f"{file_number}resultsimple.txt")
+        write_results_to_file(output_file_1, data, average_min_principal_list)
+
+        # Write the nodes within 3mm distance for each node to a different text file
+        output_file_2 = os.path.join(output_folder_2, f"{file_number}resultadvanced.txt")
+        write_nodes_within_distance_to_file(output_file_2, data, nodes_within_distance_list)

AverageEMinforeachnode.py

@@ -0,0 +1,56 @@
+def read_data(file_name):
+    with open(file_name, 'r') as file:
+        lines = file.readlines()
+        data = []
+        for line in lines:
+            split_line = line.split()
+            if len(split_line) < 2:
+                print(f"Skipping line due to insufficient data: '{line.strip()}'")
+                continue
+            try:
+                node_data = [int(split_line[0]), float(split_line[1])]
+                data.append(node_data)
+            except ValueError:
+                print(f"Skipping line due to a conversion error: '{line.strip()}'")
+                continue
+    return data
+
+def calculate_average_strains(data):
+    strains = {}
+    for node_data in data:
+        node, strain = node_data
+        if node in strains:
+            strains[node].append(strain)
+        else:
+            strains[node] = [strain]
+
+    average_strains = {}
+    for node, strain_values in strains.items():
+        average_strains[node] = sum(strain_values) / len(strain_values)
+
+    return average_strains
+
+def write_data(file_name, data):
+    with open(file_name, 'w') as file:
+        for node, average_strain in data.items():
+            line = "{:d}\t{:15.10E}\n".format(node, average_strain)
+            file.write(line)
+
+import os
+
+def process_files_average_emin(input_dir, output_dir):
+    if not os.path.exists(output_dir):
+        os.makedirs(output_dir)
+
+    files = [f for f in os.listdir(input_dir) if os.path.isfile(os.path.join(input_dir, f))]
+
+    for file in files:
+        input_file_path = os.path.join(input_dir, file)
+        output_file_name = f"{os.path.splitext(file)[0]}_averaged.txt"
+        output_file_path = os.path.join(output_dir, output_file_name)
+
+        data = read_data(input_file_path)
+        average_strains = calculate_average_strains(data)
+        write_data(output_file_path, average_strains)
+
+    print(f"Processed {len(files)} files, and saved them to '{output_dir}'.")

ResultsPostProcessor1.1.0.py

@@ -0,0 +1,187 @@
+import os
+import webbrowser
+import tkinter as tk
+from tkinter import ttk
+from tkinter import filedialog
+from tkinter import PhotoImage
+from AverageEMinforeachnode import process_files_average_emin
+from Auto3mmRadius2 import process_files_auto_3mm_radius
+
+def open_linkedin():
+    webbrowser.open_new(r"https://www.linkedin.com/in/mafazsyed")
+
+def open_github():
+    webbrowser.open_new(r"https://www.github.com/mafazsyed")
+
+def browse_folder():
+    folder = filedialog.askdirectory()
+    return folder
+
+def run_average_emin():
+    def browse_input_folder():
+        input_folder = filedialog.askdirectory()
+        input_entry.delete(0, tk.END)
+        input_entry.insert(0, input_folder)
+
+    def browse_output_folder():
+        output_folder = filedialog.askdirectory()
+        output_entry.delete(0, tk.END)
+        output_entry.insert(0, output_folder)
+
+    def start_processing():
+        input_dir = input_entry.get()
+        output_dir = output_entry.get()
+        process_files_average_emin(input_dir, output_dir)
+        emin_win.destroy()
+
+    emin_win = tk.Toplevel(root)
+    emin_win.title("Calculate the Average E. Min. Principal for Each Node")
+
+    emin_frame = ttk.Frame(emin_win, padding="20 20")
+    emin_frame.pack()
+
+    emin_label = ttk.Label(emin_frame, text="Average Minimum Principal Strain for Each Node", font=("Calibri", 12, "bold"), justify=tk.CENTER)
+    emin_label.grid(row=0, column=0, columnspan=3, pady=(0, 10))
+
+    description_label = ttk.Label(emin_frame, text="Calculate the Average E. Min. Principal for Each Node from E. Min. Principal values for Each Node Obtained from Adjacent Elements. Select an Input Folder Containing E. Min. Principal Values for Each Node & Output Folder for Averaged E. Min. Principal value for Each Node.", wraplength=450, justify=tk.CENTER)
+    description_label.grid(row=1, column=0, columnspan=2, pady=(0, 10))
+
+    input_label = ttk.Label(emin_frame, text="Input Folder:")
+    input_label.grid(row=2, column=0, padx=(0, 10), sticky=tk.W)
+    input_entry = ttk.Entry(emin_frame, width=50)
+    input_entry.grid(row=2, column=1)
+    input_browse_button = ttk.Button(emin_frame, text="Browse", command=browse_input_folder)
+    input_browse_button.grid(row=2, column=2, padx=(10, 0))
+
+    output_label = ttk.Label(emin_frame, text="Output Folder:")
+    output_label.grid(row=3, column=0, padx=(0, 10), sticky=tk.W)
+    output_entry = ttk.Entry(emin_frame, width=50)
+    output_entry.grid(row=3, column=1)
+    output_browse_button = ttk.Button(emin_frame, text="Browse", command=browse_output_folder)
+    output_browse_button.grid(row=3, column=2, padx=(10, 0))
+
+    process_button = ttk.Button(emin_frame, text="Start Processing", command=start_processing)
+    process_button.grid(row=4, column=1, pady=(20, 0))
+
+def run_auto_3mm_radius():
+    def browse_input_folder():
+        input_folder = filedialog.askdirectory()
+        input_entry.delete(0, tk.END)
+        input_entry.insert(0, input_folder)
+
+    def browse_output_folder_1():
+        output_folder = filedialog.askdirectory()
+        output_entry_1.delete(0, tk.END)
+        output_entry_1.insert(0, output_folder)
+
+    def browse_output_folder_2():
+        output_folder = filedialog.askdirectory()
+        output_entry_2.delete(0, tk.END)
+        output_entry_2.insert(0, output_folder)
+
+    def start_processing():
+        input_dir = input_entry.get()
+        output_dir_1 = output_entry_1.get()
+        output_dir_2 = output_entry_2.get()
+        process_files_auto_3mm_radius(input_dir, output_dir_1, output_dir_2)
+        auto3mm_win.destroy()
+
+    auto3mm_win = tk.Toplevel(root)
+    auto3mm_win.title("Determine Average E. Min. Principal & MFS When Considering Each Node as the Center Node")
+
+    auto3mm_frame = ttk.Frame(auto3mm_win, padding="20 20")
+    auto3mm_frame.pack()
+
+    auto3mm_label = ttk.Label(auto3mm_frame, text="Node-Centered E. Min. Principal & MFS Averaging", font=("Calibri", 12, "bold"), justify=tk.CENTER)
+    auto3mm_label.grid(row=0, column=0, columnspan=3, pady=(0, 10))
+
+    description_label = ttk.Label(auto3mm_frame, text="Determine the Average E. Min. Principal & MFS When Considering Each Node as the Center Node (of 3mm Radius). Select Input Folder Containing Files with Node ID, Node Coordinate, & E. Min. Principal value for Each Node & Output Folders for Simple & Detailed Results.", wraplength=450, justify=tk.CENTER)
+    description_label.grid(row=1, column=0, columnspan=2, pady=(0, 10))
+
+    input_label = ttk.Label(auto3mm_frame, text="Input Folder:")
+    input_label.grid(row=2, column=0, padx=(0, 10), sticky=tk.W)
+    input_entry = ttk.Entry(auto3mm_frame, width=50)
+    input_entry.grid(row=2, column=1)
+    input_browse_button = ttk.Button(auto3mm_frame, text="Browse", command=browse_input_folder)
+    input_browse_button.grid(row=2, column=2, padx=(10, 0))
+
+    output_label_1 = ttk.Label(auto3mm_frame, text="Output Folder 1, for Simple Results:")
+    output_label_1.grid(row=3, column=0, padx=(0, 10), sticky=tk.W)
+    output_entry_1 = ttk.Entry(auto3mm_frame, width=50)
+    output_entry_1.grid(row=3, column=1)
+    output_browse_button_1 = ttk.Button(auto3mm_frame, text="Browse", command=browse_output_folder_1)
+    output_browse_button_1.grid(row=3, column=2, padx=(10, 0))
+
+    output_label_2 = ttk.Label(auto3mm_frame, text="Output Folder 2, for Detailed Results:")
+    output_label_2.grid(row=4, column=0, padx=(0, 10), sticky=tk.W)
+    output_entry_2 = ttk.Entry(auto3mm_frame, width=50)
+    output_entry_2.grid(row=4, column=1)
+    output_browse_button_2 = ttk.Button(auto3mm_frame, text="Browse", command=browse_output_folder_2)
+    output_browse_button_2.grid(row=4, column=2, padx=(10, 0))
+
+    process_button = ttk.Button(auto3mm_frame, text="Start Processing", command=start_processing)
+    process_button.grid(row=5, column=1, pady=(20, 0))
+
+def about_window():
+    about_win = tk.Toplevel(root)
+    about_win.title("About")
+
+    about_frame = ttk.Frame(about_win, padding="20 20 20 20")
+    about_frame.grid(column=0, row=0, sticky=(tk.W, tk.E, tk.N, tk.S))
+
+    about_label = ttk.Label(about_frame, text="Results Post Processor", font=("Calibri", 12, "bold"))
+    about_label.grid(row=0, column=0, columnspan=2, pady=(0, 10))
+
+    description_label = ttk.Label(about_frame, text="This application:\n" "Calculates the Average of the E. Min. Principal values produced from Adjacent Elements for Each Node;\nDetermines the Average E. Min. Principal and MFS When Considering Each Node as the Center Node (of a 3mm Radius);\nFor More than One Simulation Result", wraplength=450, justify=tk.CENTER)
+    description_label.grid(row=1, column=0, columnspan=2, pady=(0, 10))
+
+    version_label = ttk.Label(about_frame, text="Version: 1.1.0", font=('Calibri', 9, 'italic'))
+    version_label.grid(row=2, column=0, columnspan=2, pady=(0, 10))
+
+    created_by_label = ttk.Label(about_frame, text="Created by Mafaz Syed", font=('Calibri', 9, 'italic'))
+    created_by_label.grid(row=3, column=0, columnspan=2, pady=(0, 10))
+
+    linkedin_logo = PhotoImage(file=r"C:\Users\mafaz\OneDrive\Desktop\ResultsPostProcessor - Copy\linkedin_logo.png.png")
+    linkedin_logo = linkedin_logo.subsample(3)  # Make the image 2 times smaller
+    linkedin_button = tk.Button(about_frame, image=linkedin_logo, command=open_linkedin, bd=0) # borderwidth set to 0
+    linkedin_button.image = linkedin_logo  # keep a reference
+    linkedin_button.grid(row=4, column=0, sticky='e')
+
+    github_logo = PhotoImage(file=r"C:\Users\mafaz\OneDrive\Desktop\ResultsPostProcessor - Copy\github_logo.png.png")
+    github_logo = github_logo.subsample(3)  # Make the image 2 times smaller
+    github_button = tk.Button(about_frame, image=github_logo, command=open_github, bd=0) # borderwidth set to 0
+    github_button.image = github_logo  # keep a reference
+    github_button.grid(row=4, column=1, sticky='w')
+
+    close_button = ttk.Button(about_frame, text="Close", command=about_win.destroy)
+    close_button.grid(row=5, column=0, columnspan=2, pady=(10, 0))
+
+root = tk.Tk()
+root.title("Results Post Processor")
+
+# Use ttk style
+style = ttk.Style()
+style.theme_use("xpnative") # choose a modern theme
+style.configure("TButton", padding=4, relief="flat",)
+
+frame = ttk.Frame(root)
+frame.pack(padx=10, pady=10)
+
+# Add a heading to the application
+heading_label = ttk.Label(frame, text="Results Post Processor", font=("Calibri", 12, "bold"))
+heading_label.pack()
+
+description_text = ("This application calculates the Average E. Min. Principal values from Adjacent Elements for Each Node, & Determines the Average E. Min. Principal & MFS When Each Node is Considered as the Center Node of a 3mm Radius, for Several Files at the Same Time. Specify Input & Output Folder Paths for Each Task.")
+description_label = ttk.Label(frame, text=description_text, wraplength=500, justify=tk.CENTER)
+description_label.pack()
+
+btn_average_emin = ttk.Button(frame, text="Calculate the Average E. Min. Principal for Each Node", command=run_average_emin)
+btn_average_emin.pack(pady=5)
+
+btn_auto_3mm_radius = ttk.Button(frame, text="Determine the Average E. Min. Principal & MFS When Considering Each Node as the Center Node (of 3mm Radius)", command=run_auto_3mm_radius)
+btn_auto_3mm_radius.pack(pady=5)
+
+about_button = ttk.Button(frame, text="About", command=about_window)
+about_button.pack(pady=5)
+
+root.mainloop()