Numpy example from Flocode Newsletter #023 - The Best Python Libraries for Civil/Structural Engineering

023_numpy_example.py

import numpy as np

# Constants
M = 140  # Moment at the cross-section (kN.m)
I = 102_000_000  # Moment of inertia in mm^4
y = np.array([-150, -100, -50, 0, 50, 100, 150])  # Distance from the neutral axis in mm

# Bending stress calculation
sigma = ((M * 1e6) * y) / I  # Convert M from kN.m to N.mm

# Print bending stress values
print("Bending Stress (N/mm^2) at different y positions (mm):") # Print header
for y_val, stress in zip(y, sigma): # Loop through the y values and corresponding stress values
    print(f"At y = {y_val} mm: {stress:.1f} N/mm^2") # Print stress values to 1 decimal places