struct_post 0.1.2__tar.gz

struct_post-0.1.2/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2025 CHUNHAO LYU
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

struct_post-0.1.2/PKG-INFO

@@ -0,0 +1,17 @@
+Metadata-Version: 2.4
+Name: struct_post
+Version: 0.1.2
+Summary: A module designed to analyse common structural test results.
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: matplotlib>=3.10.6
+Requires-Dist: openpyxl>=3.1.5
+Requires-Dist: pandas>=2.3.2
+Dynamic: license-file
+
+# struct_post
+A module designed to analyze common structural test results. The common stuctural tests include:
+    - Material coupon test
+    - Four-pont bending test
+    - Three-point bending test

struct_post-0.1.2/README.md

@@ -0,0 +1,5 @@
+# struct_post
+A module designed to analyze common structural test results. The common stuctural tests include:
+    - Material coupon test
+    - Four-pont bending test
+    - Three-point bending test

struct_post-0.1.2/postanalysis.py

@@ -0,0 +1,280 @@
+def coupon_test_analysis (testdata_file_name: str,
+                 thickness: float = 2.5, 
+                 width: float = 10, 
+                 showfig: bool = True,
+                 savefig: bool = False,
+                 lower_bound: float = 0.1, 
+                 upper_bound: float = 0.3,
+                 ):
+    """
+    Post-process a tensile coupon test and plot stress-strain curve.
+     
+    To active with Jupyter Lab, '%matplotlib widget' is required
+    
+    Args:
+        Thickness (float): Specimen thickness in mm.
+        Width (float): Specimen width in mm.
+        file_name (str): CSV file containing test data.
+        low_bound (float): Lower bound of elastic region as fraction of UTS.
+        upper_bound (float): Upper bound of elastic region as fraction of UTS.
+
+    Returns:
+        fig (matplotlib.figure.Figure): Figure object containing the plot.
+    """
+    import pandas as pd
+    import numpy as np
+    import matplotlib.pyplot as plt
+    #%matplotlib widget
+    # Constants
+    thickness # mm
+    width # mm
+    area = thickness * width  # Calculate the area of the specimen
+    
+    # Load tensile test data
+    df = pd.read_csv(testdata_file_name, header=[0])
+    
+    #df.columns = [f"{col[0]} {col[1]}" for col in df.columns]
+    
+    # Extract relevant columns
+    #time = df["Time (sec)"]
+    #displacement = df["Crosshead separation (mm)"]
+    force = df[df.columns[1]]
+    elongation = df[df.columns[2]]
+    strain = elongation # Strain in mm/mm
+    
+    # Calculate stress and strain
+    force = force * 1000 # Convert kN to N
+    stress = (force / area)  # N/m^2 or Pa
+    uts = stress.max()
+    
+    #find the data before uts
+    idx_peak = np.argmax(stress)
+    strain_up = strain[:idx_peak+1]
+    stress_up = stress[:idx_peak+1]
+    
+    #Boundary for 20% - 50% of UTS
+    lower_bound = 0.1 * uts
+    upper_bound = 0.3 * uts
+    
+    elastic_reg = (lower_bound <= stress_up) & (stress_up <= upper_bound)
+    
+    stress_ela = stress_up[elastic_reg]
+    strain_ela = strain_up[elastic_reg] 
+    
+    E, intercept = np.polyfit(strain_ela, stress_ela, 1)
+    #print(f"Young's Modulus is: {E} MPa",)
+    E_GPa = E / 1000  # Convert MPa to GPa
+    #print(f"Intercept: {intercept} MPa")
+   
+    # Select over 30% of UTS, as yield stress will over 30% uts
+    strain_new = elongation
+    stress_new = force / area
+    mask = (lower_bound <= stress)
+    strain_mask = strain_new[mask]
+    stress_mask = stress_new[mask]
+    
+    offset_decimal = 0.002  # 0.2% in decimal
+    offset_line = E * (strain_new - offset_decimal) + intercept
+
+    #Find the Yield strength
+    diff = stress_mask - offset_line
+    cross_index = np.where(diff <= 0)[0][0] 
+    x1 = strain_mask[cross_index-1]
+    x2 = strain_mask[cross_index]
+    y1 = diff[cross_index-1]
+    y2 = diff[cross_index]
+    yield_strain = x1 - y1 * (x2 - x1) / (y2 - y1)
+    yield_strength = np.interp(yield_strain, strain_new, stress_new)
+
+    #Plot
+    fig, ax = plt.subplots(figsize=(10,6))
+    
+    ax.plot(strain, stress, label='Stress-Strain Curve', color='blue')
+    ax.plot(strain_new, offset_line, label='0.2% Offset Strain Line', color='yellow',linestyle = '-.')
+    
+    ax.axhline(y=yield_strength, label=f'Yield Strength = {yield_strength:.2f} MPa', color='green', linestyle = 'dotted')
+    ax.axhline(y=uts, color='red', linestyle = '--', label=f'UTS = {uts:.2f} MPa')
+    ax.plot(yield_strain, yield_strength, 'ro', label='Yield Point')
+    
+    ax.set_xlabel('Strain (mm/mm)')
+    ax.set_ylabel('Stress (MPa)')
+    ax.set_title(testdata_file_name[:-4]+' Stress-Strain Curve with Mechanical Properties')
+    ax.legend()
+    ax.grid(True)
+    ax.set_ylim(0, 1.2 *uts)
+
+    #Show fig or not
+    if showfig == True:
+        plt.show()
+    else:
+        plt.close(fig)
+
+    # Save fig or not
+    if savefig == True:
+        fig.savefig(testdata_file_name[:-4], dpi=300, bbox_inches='tight')
+    
+    # Print results
+    print(f"Sample: {testdata_file_name[:-4]}")
+    print(f"Young's Modulus (E): {E:.2f} MPa")
+    print(f"Ultimate Tensile Strength (UTS): {uts:.2f} MPa")
+    print(f"Yield Strength: {yield_strength:.2f} MPa")
+    print('-' * 40)
+    
+    # Prepare results dictionary
+    results = {
+        "E_MPa": E,
+        "UTS_MPa": uts,
+        "Yield_Strength_MPa": yield_strength
+    }
+    return testdata_file_name[:-4], results
+
+def coupon_sample_geodata_read (Excelfile_name: str):
+    '''
+    Reads sample geometric properties from an Excel file.
+
+    Parameters
+    ----------
+    Excelfile_name : str
+        The name of the Excel file to read, located in the current working directory.
+
+    Returns
+    -------
+    list
+        A list containing:
+        - sample_file_name (sample name from cell A2 with ".csv" extension)
+        - thickness (value from cell C2)
+        - width (value from cell B2)
+    '''
+    
+    from openpyxl import load_workbook
+    from pathlib import Path
+    sample_data_file = Path.cwd()/Excelfile_name
+    wb = load_workbook(sample_data_file)
+    ws = wb.worksheets[0]
+    samples = []
+    for row in ws.iter_rows(min_row=2, values_only=True):  
+        if all(cell is None for cell in row):  # skip empty
+            continue
+        # 取前三列数据
+        sample_file_name = str (row[0] + ".csv")
+        sample_name = coupon_SampleDetails(row[0], row[1], row[2], sample_file_name)
+        samples.append(sample_name)
+    return samples
+
+def coupon_batch_analysis(Coupon_geodata: str,
+                   showfig: bool = True,
+                   savefig: bool = False):
+    """
+    Perform batch analysis on a list of samples and return the results.
+
+    Parameters
+    ----------
+    Coupon_geodata : list
+        A list of SampleDetails objects, each containing:
+        - sample_file_name : Name of the sample file
+        - thickness        : Sample thickness
+        - width            : Sample width
+    savefig : bool, optional
+        Whether to save the figures generated during analysis. Default is False.
+
+    Returns
+    -------
+    list
+        A list of SampleAnalysisResults objects, each containing:
+        - modulus_of_elasticity       : Elastic modulus (E_GPa)
+        - ultimate_tensile_strength   : Ultimate tensile strength (UTS_MPa)
+        - yield_Strength              : Yield strength (Yield_Strength_MPa)
+    """
+    SARS = []
+    for Coupon_detail in Coupon_geodata:
+        result = coupon_test_analysis(Coupon_detail.sample_file_name, Coupon_detail.thickness, Coupon_detail.width, showfig, savefig)
+        SAR = coupon_SampleAnalysisResults(result[0], result[1]['E_MPa'], result[1]['UTS_MPa'], result[1]['Yield_Strength_MPa'])
+        SARS.append(SAR)
+    return SARS
+
+def coupon_results_save(Excelfile_name: str, analysis_results: list):
+    """
+    Save a list of sample analysis results into an Excel file, matching by sample name.
+
+    Parameters
+    ----------
+    Excelfile_name: str
+        Path to the Excel file to save the results.
+    analysis_results : list
+        A list of SampleAnalysisResults objects, each containing:
+        - sample_name
+        - modulus_of_elasticity
+        - ultimate_tensile_strength
+        - yield_Strength
+
+    Notes
+    -----
+    - Assumes the Excel file has sample names in column A, starting from row 2.
+    - Data will be written starting from column D (fourth column).
+    - Overwrites the original Excel file.
+    """
+    from openpyxl import load_workbook
+
+    wb = load_workbook(Excelfile_name)
+    ws = wb.worksheets[0]
+
+    # Build a dictionary for faster lookup
+    result_dict = {res.sample_name: res for res in analysis_results}
+
+    for row_idx, row in enumerate(ws.iter_rows(min_row=2, max_col=ws.max_column), start=2):
+        excel_sample_name = row[0].value
+        if excel_sample_name in result_dict:
+            analysis_result = result_dict[excel_sample_name]
+            values_to_write = [
+                analysis_result.modulus_of_elasticity,
+                analysis_result.ultimate_tensile_strength,
+                analysis_result.yield_Strength
+            ]
+            for col_idx, value in enumerate(values_to_write, start=4):  # D列开始
+                ws.cell(row=row_idx, column=col_idx, value=float(value))
+
+    wb.save(Excelfile_name)
+    print('The coupon test data analysis is complete.')
+
+from dataclasses import dataclass
+@dataclass # this thing is called a "decorator"
+class coupon_SampleDetails:
+    """
+    Holds basic information for a sample.
+
+    Attributes
+    ----------
+    sample_name : str
+        Name or ID of the sample.
+    width : float
+        Width of the sample (geometric property, in consistent units).
+    thickness : float
+        Thickness of the sample (geometric property, in consistent units).
+    sample_file_name : str
+        Name of the associated data file for the sample, e.g., CSV file.
+    """
+    sample_name: str
+    width: float
+    thickness: float
+    sample_file_name: str
+
+@dataclass # this thing is called a "decorator"
+class coupon_SampleAnalysisResults:
+    """
+    Stores the analysis results for a sample after mechanical testing.
+
+    Attributes
+    ----------
+    sample_name : str
+        Name or ID of the sample corresponding to the analysis.
+    modulus_of_elasticity : float
+        Elastic modulus of the sample (E, in GPa).
+    ultimate_tensile_strength : float
+        Maximum tensile strength of the sample (UTS, in MPa).
+    yield_Strength : float
+        Yield strength of the sample (in MPa).
+    """
+    sample_name: str
+    modulus_of_elasticity: float
+    ultimate_tensile_strength: float
+    yield_Strength: float

struct_post-0.1.2/pyproject.toml

@@ -0,0 +1,11 @@
+[project]
+name = "struct_post"
+version = "0.1.2"
+description = "A module designed to analyse common structural test results."
+readme = "README.md"
+requires-python = ">=3.10"
+dependencies = [
+    "matplotlib>=3.10.6",
+    "openpyxl>=3.1.5",
+    "pandas>=2.3.2",
+]

struct_post-0.1.2/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

struct_post-0.1.2/struct_post.egg-info/PKG-INFO

@@ -0,0 +1,17 @@
+Metadata-Version: 2.4
+Name: struct_post
+Version: 0.1.2
+Summary: A module designed to analyse common structural test results.
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+License-File: LICENSE
+Requires-Dist: matplotlib>=3.10.6
+Requires-Dist: openpyxl>=3.1.5
+Requires-Dist: pandas>=2.3.2
+Dynamic: license-file
+
+# struct_post
+A module designed to analyze common structural test results. The common stuctural tests include:
+    - Material coupon test
+    - Four-pont bending test
+    - Three-point bending test

struct_post-0.1.2/struct_post.egg-info/SOURCES.txt

@@ -0,0 +1,9 @@
+LICENSE
+README.md
+postanalysis.py
+pyproject.toml
+struct_post.egg-info/PKG-INFO
+struct_post.egg-info/SOURCES.txt
+struct_post.egg-info/dependency_links.txt
+struct_post.egg-info/requires.txt
+struct_post.egg-info/top_level.txt

struct_post-0.1.2/struct_post.egg-info/dependency_links.txt

@@ -0,0 +1 @@
+

struct_post-0.1.2/struct_post.egg-info/requires.txt

@@ -0,0 +1,3 @@
+matplotlib>=3.10.6
+openpyxl>=3.1.5
+pandas>=2.3.2

struct_post-0.1.2/struct_post.egg-info/top_level.txt

@@ -0,0 +1 @@
+postanalysis