struct_post 0.1.4.3__py3-none-any.whl → 0.1.7__py3-none-any.whl

struct_post/beam.py

@@ -1,20 +1,133 @@
+from pandas import DataFrame
 import pandas as pd
-def four_point_bending (data: pd.DataFrame, 
+
+def beam_lvmdata_file(file_name, force_index=2, displ_index=3,
+             mid_lvdt_ind=None, left_lvdt_ind=None, right_lvdt_ind=None):
+    """
+    Read a .lvm file and extract force, displacement, and LVDT data in groups.
+
+    Parameters
+    ----------
+    file_name : str
+        Path to the .lvm file to read.
+    force_index : int, default=2
+        1-based column index for the force data.
+    displ_index : int, default=3
+        1-based column index for the displacement (stroke) data.
+    mid_lvdt_ind : int or list of int, optional
+        Column index/indices for the middle LVDTs.
+    left_lvdt_ind : int or list of int, optional
+        Column index/indices for the left LVDTs.
+    right_lvdt_ind : int or list of int, optional
+        Column index/indices for the right LVDTs.
+
+    Returns
+    -------
+    df_new : pandas.DataFrame
+        Extracted data as a DataFrame with columns renamed:
+        ['Force (kN)', 'Stroke (mm)', 'Mid_1 (mm)', 'Mid_2 (mm)',
+         'Left_1 (mm)', 'Left_2 (mm)', 'Right_1 (mm)', 'Right_2 (mm)'].
+        Only existing columns are included; missing columns are ignored.
+    file_base_name : str
+        Base name of the file without the '.lvm' extension.
+
+    Notes
+    -----
+    - Automatically detects the start of the data by locating the last
+      occurrence of '***End_of_Header***' in the file.
+    - Converts all columns to numeric values; non-numeric entries
+      are coerced to NaN.
+    - The last column is dropped assuming it is empty.
+    - Single integer indices are automatically converted to lists.
+    - LVDT columns are combined in the order: mid -> left -> right.
+    """
+    
+    mid_lvdt_ind = mid_lvdt_ind or []
+    left_lvdt_ind = left_lvdt_ind or []
+    right_lvdt_ind = right_lvdt_ind or []
+    
+    if isinstance(mid_lvdt_ind, int):
+        mid_lvdt_ind = [mid_lvdt_ind]
+    if isinstance(left_lvdt_ind, int):
+        left_lvdt_ind = [left_lvdt_ind]
+    if isinstance(right_lvdt_ind, int):
+        right_lvdt_ind = [right_lvdt_ind]
+    
+    lvdt_indices = mid_lvdt_ind + left_lvdt_ind + right_lvdt_ind
+    
+    with open(file_name, 'r') as f:
+        lines = f.readlines()
+    start_idx = max(i for i, line in enumerate(lines) if "***End_of_Header***" in line) + 1
+    df = pd.read_csv(file_name, sep="\t", skiprows=start_idx)
+    df = df.apply(pd.to_numeric, errors='coerce')
+    df = df.drop(df.columns[-1], axis=1)
+    
+    user_indices = [force_index, displ_index] + lvdt_indices
+    existing_cols = [df.columns[i-1] for i in user_indices if i-1 < len(df.columns)]
+    
+    df_new = df[existing_cols].copy()
+    
+    col_names = ['Force (kN)','Stroke (mm)','Mid_1 (mm)','Mid_2 (mm)','Left_1 (mm)','Left_2 (mm)','Right_1 (mm)','Right_2 (mm)']
+    df_new.columns = col_names[:len(existing_cols)]
+    
+    return df_new, file_name[:-4]
+
+def four_point_bending (data: DataFrame, 
                         width: float, 
                         depth: float, 
                         beam_span: float):
+    """
+    Perform four-point bending analysis to calculate apparent and true modulus of elasticity.
+
+    Parameters
+    ----------
+    data : tuple
+        A tuple containing:
+        - data[0] : pandas.DataFrame
+            Experimental dataset with columns:
+            - 'Moog Force_kN' : Applied load (kN).
+            - 'LVDT 1_mm' ... 'LVDT 6_mm' : Deflections from six LVDTs (mm).
+        - data[1] : str
+            Sample name identifier.
+    width : float
+        Beam specimen width (mm).
+    depth : float
+        Beam specimen depth (mm).
+    beam_span : float
+        Beam span length (mm).
+
+    Returns
+    -------
+    tuple
+        (sample_name, results) where:
+        - sample_name : str
+            Name of the processed sample.
+        - results : dict
+            Dictionary containing:
+            - "E_app" : float
+                Apparent modulus of elasticity (MPa).
+            - "E_true" : float
+                True modulus of elasticity (MPa).
+
+    Notes
+    -----
+    - Apparent modulus (E_app) is calculated from mid-span deflection (LVDT 3 & 4).
+    - True modulus (E_true) is calculated from relative deflection (mid-span vs. supports).
+    - Load range for regression is limited to 10–40% of ultimate load.
+    """
+    
     import pandas as pd
     import numpy as np
 
-    #geo
+    #Experimental test data post-process
     sample_name = data[1]
-    force = data[0]['Moog Force_kN'] * 1000
-    delta_1 = abs(data[0]['LVDT 1_mm'])
-    delta_2 = abs(data[0]['LVDT 2_mm'])
-    delta_3 = abs(data[0]['LVDT 3_mm'])
-    delta_4 = abs(data[0]['LVDT 4_mm'])
-    delta_5 = abs(data[0]['LVDT 5_mm'])
-    delta_6 = abs(data[0]['LVDT 6_mm'])
+    force = data[0]['Force (kN)'] * 1000
+    delta_1 = abs(data[0]['Left_1 (mm)'])
+    delta_2 = abs(data[0]['Left_2 (mm)'])
+    delta_3 = abs(data[0]['Mid_1 (mm)'])
+    delta_4 = abs(data[0]['Mid_2 (mm)'])
+    delta_5 = abs(data[0]['Right_1 (mm)'])
+    delta_6 = abs(data[0]['Right_2 (mm)'])
 
     F_ult = force.max()
     f_b = (F_ult * beam_span) / (width * depth **2) #MPa
@@ -31,10 +144,10 @@ def four_point_bending (data: pd.DataFrame,
     
     F_ms = force[calcs_reg]
     delta_ms_calcs = delta_ms[calcs_reg]
-    delat_rel_calcs = delta_rel[calcs_reg]
+    delta_rel_calcs = delta_rel[calcs_reg]
     
     Delta_ms, intercept_ms = np.polyfit(delta_ms_calcs,F_ms,1)
-    Delta_rel, intercept_rel = np.polyfit(delat_rel_calcs,F_ms,1)
+    Delta_rel, intercept_rel = np.polyfit(delta_rel_calcs,F_ms,1)
     
     E_app = (23/108) * (beam_span/depth)**3 * Delta_ms * (1/width)
     E_true = (1/36) *  (beam_span/depth)**3 * Delta_rel * (1/width)
@@ -47,3 +160,4 @@ def four_point_bending (data: pd.DataFrame,
     print(f"Sample Name: {sample_name}")
     print('-' * 40)
     return sample_name, results
+

struct_post/coupon.py

@@ -1,10 +1,44 @@
-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,
+def coupon_csvdata_read(testdata_file_name: str,
+                         force_index: int,
+                         strain_index: int):
+    """
+    Read tensile coupon test data from a CSV file and extract force and strain columns.
+
+    Args:
+        testdata_file_name (str): Path to the CSV file containing the test data.
+        force_index (int): 1-based column index for the force data.
+        strain_index (int): 1-based column index for the strain data.
+
+    Returns:
+        tuple:
+            df_new (pd.DataFrame): DataFrame containing only 'Force' and 'Strain' columns.
+            sample_name (str): Name of the sample extracted from the file name (without extension).
+
+    Notes:
+        - The CSV file is expected to have a header in the first row.
+        - The second row is skipped during reading (`skiprows=[1]`).
+        - The first column is treated as the index column (`index_col=0`).
+        - Force and strain columns are extracted based on the provided 1-based indices.
+    """
+    
+    import pandas as pd
+    df = pd.read_csv(testdata_file_name, skiprows=[1])
+    Force = df[df.columns[force_index - 1]]
+    Strain = df[df.columns[strain_index - 1]]
+    df_new = pd.concat([Force, Strain], axis=1)
+    df_new.columns = ['Force', 'Strain']
+    sample_name = testdata_file_name[:-4]
+    return sample_name, df_new
+
+def coupon_test_analysis (sample_name: str,
+                          Force: float,
+                          Strain: float,
+                          thickness: float = 2.5,
+                          width: float = 10,
+                          showfig: bool = True,
+                          savefig: bool = False,
+                          low_bound: float = 0.1,
+                          up_bound: float = 0.3,
                  ):
     """
     Post-process a tensile coupon test and plot stress-strain curve.
@@ -25,36 +59,27 @@ def coupon_test_analysis (testdata_file_name: str,
     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
+
+    strain = Strain # Strain in mm/mm
     
     # Calculate stress and strain
-    force = force * 1000 # Convert kN to N
-    stress = (force / area)  # N/m^2 or Pa
+    force = Force * 1000 # Convert kN to N
+    stress = (force / area)  # N/mm^2 or MPa
     uts = stress.max()
     
     #find the data before uts
     idx_peak = np.argmax(stress)
-    strain_up = strain[:idx_peak+1]
+    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
+    #Boundary for [low bound] - [up bound] of uts
+    lower_bound = low_bound * uts
+    upper_bound = up_bound * uts
     
     elastic_reg = (lower_bound <= stress_up) & (stress_up <= upper_bound)
     
@@ -66,8 +91,8 @@ def coupon_test_analysis (testdata_file_name: str,
     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
+    # Select over [lower bound] of UTS, as yield stress will over [lower bound] uts
+    strain_new = strain
     stress_new = force / area
     mask = (lower_bound <= stress)
     strain_mask = strain_new[mask]
@@ -78,7 +103,7 @@ def coupon_test_analysis (testdata_file_name: str,
 
     #Find the Yield strength
     diff = stress_mask - offset_line
-    cross_index = np.where(diff <= 0)[0][0] 
+    cross_index = np.where(diff <= 0)[0][0] +1
     x1 = strain_mask[cross_index-1]
     x2 = strain_mask[cross_index]
     y1 = diff[cross_index-1]
@@ -98,7 +123,7 @@ def coupon_test_analysis (testdata_file_name: str,
     
     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.set_title(sample_name + ' Stress-Strain Curve with Mechanical Properties')
     ax.legend()
     ax.grid(True)
     ax.set_ylim(0, 1.2 *uts)
@@ -111,10 +136,10 @@ def coupon_test_analysis (testdata_file_name: str,
 
     # Save fig or not
     if savefig == True:
-        fig.savefig(testdata_file_name[:-4], dpi=300, bbox_inches='tight')
+        fig.savefig(sample_name, dpi=300, bbox_inches='tight')
     
     # Print results
-    print(f"Sample: {testdata_file_name[:-4]}")
+    print(f"Sample: {sample_name}")
     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")
@@ -126,7 +151,7 @@ def coupon_test_analysis (testdata_file_name: str,
         "UTS_MPa": uts,
         "Yield_Strength_MPa": yield_strength
     }
-    return testdata_file_name[:-4], results
+    return sample_name, results
 
 def coupon_sample_geodata_read (Excelfile_name: str):
     '''
@@ -155,39 +180,61 @@ def coupon_sample_geodata_read (Excelfile_name: str):
     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):
+                          force_index: float,
+                          strain_index:float,
+                          showfig: bool = True,
+                          savefig: bool = False):
     """
-    Perform batch analysis on a list of samples and return the results.
+    Perform batch analysis on a list of tensile coupon samples and return their results.
 
     Parameters
     ----------
-    Coupon_geodata : list
+    Coupon_geodata : list of SampleDetails
         A list of SampleDetails objects, each containing:
-        - sample_file_name : Name of the sample file
-        - thickness        : Sample thickness
-        - width            : Sample width
+        - sample_file_name : str
+            Path to the sample CSV file.
+        - thickness : float
+            Sample thickness (mm).
+        - width : float
+            Sample width (mm).
+    force_index : int
+        Column index of the Force data in the CSV file (1-based).
+    strain_index : int
+        Column index of the Strain data in the CSV file (1-based).
+    showfig : bool, optional
+        Whether to display stress-strain plots during analysis. Default is True.
     savefig : bool, optional
-        Whether to save the figures generated during analysis. Default is False.
+        Whether to save the stress-strain plots to files. Default is False.
 
     Returns
     -------
-    list
+    list of SampleAnalysisResults
         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)
+        - sample_name : str
+            Name of the processed sample.
+        - modulus_of_elasticity : float
+            Elastic modulus (E_GPa).
+        - ultimate_tensile_strength : float
+            Ultimate tensile strength (UTS_MPa).
+        - yield_strength : float
+            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)
+        csvdata = coupon_csvdata_read(Coupon_detail.sample_file_name,force_index,strain_index)
+        result = coupon_test_analysis(Coupon_detail.sample_file_name[:-4], 
+                                      csvdata[1]['Force'],
+                                      csvdata[1]['Strain'],
+                                      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
@@ -237,7 +284,7 @@ def coupon_results_save(Excelfile_name: str, analysis_results: list):
     print('The coupon test data analysis is complete.')
 
 from dataclasses import dataclass
-@dataclass # this thing is called a "decorator"
+@dataclass
 class coupon_SampleDetails:
     """
     Holds basic information for a sample.
@@ -258,7 +305,7 @@ class coupon_SampleDetails:
     thickness: float
     sample_file_name: str
 
-@dataclass # this thing is called a "decorator"
+@dataclass
 class coupon_SampleAnalysisResults:
     """
     Stores the analysis results for a sample after mechanical testing.

{struct_post-0.1.4.3.dist-info → struct_post-0.1.7.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: struct_post
-Version: 0.1.4.3
+Version: 0.1.7
 Summary: A module designed to analyse common structural test results.
 Requires-Python: >=3.10
 Description-Content-Type: text/markdown
@@ -8,6 +8,7 @@ License-File: LICENSE
 Requires-Dist: matplotlib>=3.10.6
 Requires-Dist: openpyxl>=3.1.5
 Requires-Dist: pandas>=2.3.2
+Requires-Dist: requests>=2.32.5
 Dynamic: license-file
 
 # struct_post

struct_post-0.1.7.dist-info/RECORD

@@ -0,0 +1,9 @@
+struct_post/__init__.py,sha256=aNaAaoGGLodY5FTsVGCnUhzuXW9F50Mi_7xm-nQ-xFg,182
+struct_post/beam.py,sha256=ftoh-J_7SCJUdptUnLgCXrvvHSKyACr0D8ZFJT0358s,5748
+struct_post/coupon.py,sha256=vLHHSc37DaZFYQmx8zDHpa0S0CEgPlfbtedx9gGZY6o,11886
+struct_post/read_lvm_file.py,sha256=3x-vbkgEQ-G3f939euIMqwjzVEt9zHtvNneYd1-ViDw,2539
+struct_post-0.1.7.dist-info/licenses/LICENSE,sha256=0KCwPGoSiY1zA6U97FZJob9MJbaJOo6sRZF-xv6A1YI,1089
+struct_post-0.1.7.dist-info/METADATA,sha256=i34uSxS4iIRHSFZHE6qXm2b2OxY9nNjaJ0Wc8oZj6v0,689
+struct_post-0.1.7.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+struct_post-0.1.7.dist-info/top_level.txt,sha256=jgiqPN0SMTMUomzOey9HfIrScTFmOzyuLpwIJ-Sz1-0,12
+struct_post-0.1.7.dist-info/RECORD,,

struct_post-0.1.4.3.dist-info/RECORD

@@ -1,9 +0,0 @@
-struct_post/__init__.py,sha256=aNaAaoGGLodY5FTsVGCnUhzuXW9F50Mi_7xm-nQ-xFg,182
-struct_post/beam.py,sha256=y-dV-taU1JNcyZYVTeJDe0i2aOkDoCUp4Ar5I8QC9iw,1495
-struct_post/coupon.py,sha256=pi-aGs-EnKosq7QwPR1unskTVoNZ8kZ3o3gQjlFSf34,9835
-struct_post/read_lvm_file.py,sha256=3x-vbkgEQ-G3f939euIMqwjzVEt9zHtvNneYd1-ViDw,2539
-struct_post-0.1.4.3.dist-info/licenses/LICENSE,sha256=0KCwPGoSiY1zA6U97FZJob9MJbaJOo6sRZF-xv6A1YI,1089
-struct_post-0.1.4.3.dist-info/METADATA,sha256=R4qNtqJsWhPf2MLAK6z063xNdw3gbGOsUTUg52nozfE,658
-struct_post-0.1.4.3.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
-struct_post-0.1.4.3.dist-info/top_level.txt,sha256=jgiqPN0SMTMUomzOey9HfIrScTFmOzyuLpwIJ-Sz1-0,12
-struct_post-0.1.4.3.dist-info/RECORD,,