structdyn 0.5.1__py3-none-any.whl

docs/source/conf.py

@@ -0,0 +1,37 @@
+# Configuration file for the Sphinx documentation builder.
+#
+# For the full list of built-in configuration values, see the documentation:
+# https://www.sphinx-doc.org/en/master/usage/configuration.html
+
+# -- Project information -----------------------------------------------------
+# https://www.sphinx-doc.org/en/master/usage/configuration.html#project-information
+import os
+import sys
+
+sys.path.insert(0, os.path.abspath("../../"))
+
+project = "structdyn"
+copyright = "2026, Abinash Mandal"
+author = "Abinash Mandal"
+release = "0.1.0"
+
+# -- General configuration ---------------------------------------------------
+# https://www.sphinx-doc.org/en/master/usage/configuration.html#general-configuration
+
+extensions = [
+    "sphinx.ext.autodoc",
+    "sphinx.ext.napoleon",  # for NumPy-style docstrings
+    "sphinx.ext.viewcode",
+]
+
+
+templates_path = ["_templates"]
+exclude_patterns = []
+
+
+# -- Options for HTML output -------------------------------------------------
+# https://www.sphinx-doc.org/en/master/usage/configuration.html#options-for-html-output
+
+html_theme = "sphinx_rtd_theme"
+
+html_static_path = ["_static"]

examples/mdf/eg1_mdf.py

@@ -0,0 +1,21 @@
+# Example 10.4; Chopra A. K., Dynamics of structure, 5th edn
+import numpy as np
+from structdyn.mdf.mdf import MDF
+
+# Define mdf system: 2-story shear building
+masses = [2, 1]
+stiffness = [2, 1]
+mdf = MDF.from_shear_building(masses, stiffness)
+# print(mdf.M, mdf.K)
+
+# Modal analysis
+omega, phi = mdf.modal.modal_analysis(
+    dof_normalize=-1
+)  # -1 means: mode shape normalized based on roof dof
+
+print("Natural Frequencies (rad/s):", omega)  # omegas = [0.70710678 1.41421356]
+print("Mode shapes:", phi)
+
+# u = np.array([1, 1])
+# qn = mdf.modal.modal_coordinates(u)
+# print("Modal coordinates:", qn)

examples/mdf/eg2_mdf.py

@@ -0,0 +1,26 @@
+# Example 10.9; Chopra A. K., Dynamics of structure, 5th edn
+import numpy as np
+from structdyn.mdf.mdf import MDF
+import matplotlib.pyplot as plt
+
+# Define mdf mdf: 2-story shear building
+masses = [2, 1]
+stiffness = [2, 1]
+mdf = MDF.from_shear_building(masses, stiffness)
+
+# Modal analysis
+omega, phi = mdf.modal.modal_analysis(
+    dof_normalize=-1
+)  # -1 means: mode shape normalized based on roof dof
+
+# Initial conditions
+u0 = [0.5, 1]
+v0 = [0, 0]
+time = np.arange(0, 20, 0.01)
+
+# Run analysis
+response = mdf.modal.free_vibration_response(u0, v0, time)
+print(response["u1"].abs().max())  # result = 0.49999999999999994
+
+response.plot(x="time", y="u1")
+plt.show()

examples/mdf/eg_cdm.py

@@ -0,0 +1,19 @@
+import numpy as np
+from structdyn.mdf.mdf import MDF
+import matplotlib.pyplot as plt
+
+# Define MDF system
+masses = [2 * 45594, 45594]
+stiffness = [2 * 18 * 10**5, 18 * 10**5]
+mdf = MDF.from_shear_building(masses, stiffness)
+
+# Define external load
+dt = 0.01
+time = np.arange(0, 1.01, dt)
+load = 50 * np.sin(np.pi * time / 0.6) * 1000
+load[time >= 0.6] = 0
+load = np.column_stack((load, np.zeros_like(load)))
+
+res = mdf.find_response(time, load)
+print(res)
+print(res.shape)

examples/mdf/eg_damping_mat.py

@@ -0,0 +1,14 @@
+# Example 11.4; Chopra A. K., Dynamics of structure, 5th edn
+from structdyn.mdf.mdf import MDF
+
+# Define MDF system
+masses = [2e5, 2e5, 1e5]
+stiffness = [1e8, 1e8, 1e8]
+mdf = MDF.from_shear_building(masses, stiffness)
+
+# Create damping matrix based on modal damping ratios; zeta
+mdf.set_modal_damping(zeta=[0.05, 0.05, 0.05])
+
+# Damping matrix
+print(mdf.C)
+print(mdf.C.flatten()[-1])  # result = 287983.44117400126

examples/mdf/eg_newmark.py

@@ -0,0 +1,18 @@
+import numpy as np
+from structdyn.mdf.mdf import MDF
+import matplotlib.pyplot as plt
+
+# Define MDF system
+masses = [2 * 45594, 45594]
+stiffness = [2 * 18 * 10**5, 18 * 10**5]
+mdf = MDF.from_shear_building(masses, stiffness)
+
+# Define external load
+dt = 0.01
+time = np.arange(0, 1.01, dt)
+load = 50 * np.sin(np.pi * time / 0.6) * 1000
+load[time >= 0.6] = 0
+load = np.column_stack((load, np.zeros_like(load)))
+
+res = mdf.find_response(time, load, method="newmark_beta")
+print(res)

examples/mdf/eg_newmark_gm.py

@@ -0,0 +1,32 @@
+# Example 16.1; Chopra A. K., Dynamics of structure, 5th edn
+import numpy as np
+from structdyn.mdf.mdf import MDF
+import matplotlib.pyplot as plt
+from structdyn.ground_motions.ground_motion import GroundMotion
+
+# Define MDF system
+masses = np.ones(5) * 0.2591 * 1e5
+stiffness = np.ones(5) * 100 * 1e5
+mdf = MDF.from_shear_building(masses, stiffness)
+mdf.set_modal_damping(zeta=np.ones(5) * 0.05)  # set damping matrix
+
+# Define ground motion as given in the same example
+time = np.arange(0, 2.01, 0.1)
+acc = (1.93 / 9.81) * np.sin(2 * np.pi * time)
+acc[time > 1.0] = 0
+gm = GroundMotion.from_arrays(acc, 0.1)
+
+# Influence vector
+inf_vec = np.ones(5)
+
+# Run analysis
+response = mdf.find_response_ground_motion(
+    gm, inf_vec, method="newmark_beta", use_modal=True, n_modes=2
+)
+# Note we can get more accurate results by considering more modes. Also if we remove use_modal, then direct integration of coupled equations would give more accurate results.
+print(response)
+print(response["u5"][20])  # result = 0.05888401249961952
+
+# Plot response
+# res.plot(x="time", y="u5")
+# plt.show()

examples/sdf/eg_analytical.py

@@ -0,0 +1,18 @@
+from structdyn import SDF
+import matplotlib.pyplot as plt
+from structdyn.sdf.analytical_methods.analytical_response import AnalyticalResponse
+
+# Create SDF object
+sdf = SDF(m=1.0, k=100.0, ji=0.05)
+analytical = AnalyticalResponse(sdf)  # analytical solver
+
+# Free vibration
+df_free = analytical.free_vibration(u0=0.01, v0=0.0)
+# print(df_free)
+df_free.plot(x="time", y="displacement")
+plt.show()
+
+# Harmonic sine forcing
+df_harm = analytical.harmonic_response(p0=10.0, w=5.0, excitation="sine")
+df_harm.plot(x="time", y="displacement")
+plt.show()

examples/sdf/eg_cdm.py

@@ -0,0 +1,23 @@
+# Example 5.2; Chopra A. K., Dynamics of structure, 5th edn
+from structdyn import SDF
+import numpy as np
+import matplotlib.pyplot as plt
+
+# Define external load
+dt = 0.1
+time_steps = np.arange(0, 1.01, dt)
+load_values = 50 * np.sin(np.pi * time_steps / 0.6) * 1000
+load_values[time_steps >= 0.6] = 0
+
+# Create SDF object
+sdf = SDF(45594, 18 * 10**5, 0.05)
+# sdf = SDF(45594, 18 * 10**5, 0.05, fd="elastoplastic", uy=0.02, fy=36000)     #for non-linear analysis
+
+# Analysis
+responses = sdf.find_response(
+    time_steps,
+    load_values,
+    method="central_difference",
+)
+print(responses)
+print(responses["displacement"][10])  # result is -0.03812718853606678

examples/sdf/eg_el_centro.py

@@ -0,0 +1,21 @@
+# Section 6.4 (& Figure 6.4.1a); Chopra A. K., Dynamics of structure, 5th edn
+from structdyn import SDF
+import numpy as np
+from structdyn.ground_motions.ground_motion import GroundMotion
+from structdyn.utils.helpers import elcentro_chopra
+
+# Define el centro ground motion from Chopra's book- Appendix 6
+elc = elcentro_chopra()
+gm = GroundMotion.from_arrays(elc["acc (g)"], 0.02)
+
+# Create SDF object
+m = 1
+ji = 0.02
+t_n = 0.5
+w_n = 2 * np.pi / t_n
+k = w_n**2 * m
+sdf = SDF(m, k, ji)
+
+# Analysis
+results = sdf.find_response_ground_motion(gm, method="interpolation")
+print(results["displacement"].abs().max())  # result is 0.0679400697200714

examples/sdf/eg_ground_motion.py

@@ -0,0 +1,24 @@
+import numpy as np
+from structdyn.ground_motions.ground_motion import GroundMotion
+from structdyn.utils.helpers import elcentro_chopra
+from pathlib import Path
+
+# Method 1: directly from Chopra's book- Appendix 6
+elc = elcentro_chopra()
+gm = GroundMotion.from_arrays(elc["acc (g)"], 0.02)
+print(gm.to_dataframe())
+
+# Method 2: from numpy arrays
+dt = 0.01
+time = np.arange(0, 2.01, dt)
+acc = (1.93 / 9.81) * np.sin(2 * np.pi * time)
+gm = GroundMotion.from_arrays(acc, dt)
+print(gm.to_dataframe())
+
+# Method 3: from event (in-built PEER database)
+gm = GroundMotion.from_event("northridge_sylmar_1994", component="hor1")
+print(gm.to_dataframe())
+
+# Method 4: from .at2 file (PEER database format)
+# file_path = "path to your file"
+# gm = GroundMotion.from_at2(file_path)

examples/sdf/eg_interpolation.py

@@ -0,0 +1,21 @@
+# Example 5.1; Chopra A. K., Dynamics of structure, 5th edn
+from structdyn import SDF
+import numpy as np
+
+# Define external load
+dt = 0.1
+time_steps = np.arange(0, 1.01, dt)
+load_values = 50 * np.sin(np.pi * time_steps / 0.6) * 1000
+load_values[time_steps >= 0.6] = 0
+
+# Create SDF object
+sdf = SDF(45594, 18 * 10**5, 0.05)
+
+# Run analysis
+responses = sdf.find_response(
+    time_steps,
+    load_values,
+    method="interpolation",
+)
+print(responses)
+print(responses["displacement"][10])  # result is -0.034534260954800985

examples/sdf/eg_newmark.py

@@ -0,0 +1,19 @@
+# Example 5.4; Chopra A. K., Dynamics of structure, 5th edn
+from structdyn import SDF
+import numpy as np
+
+# Define external load
+dt = 0.1
+time_steps = np.arange(0, 1.01, dt)
+load_values = 50 * np.sin(np.pi * time_steps / 0.6) * 1000
+load_values[time_steps >= 0.6] = 0
+
+# Create SDF object
+sdf = SDF(45594, 18 * 10**5, 0.05)
+
+# Analysis
+responses = sdf.find_response(
+    time_steps, load_values, method="newmark_beta", acc_type="linear"
+)
+print(responses)
+print(responses["displacement"][10])  # result is -0.03391195470077432

examples/sdf/eg_newmark_non_linear.py

@@ -0,0 +1,19 @@
+# Example 5.5; Chopra A. K., Dynamics of structure, 5th edn
+from structdyn import SDF
+import numpy as np
+
+# Define external load
+dt = 0.1
+time_steps = np.arange(0, 1.01, dt)
+load_values = 50 * np.sin(np.pi * time_steps / 0.6) * 1000
+load_values[time_steps >= 0.6] = 0
+
+# Create SDF object
+sdf = SDF(45594, 18 * 10**5, 0.05, fd="elastoplastic", uy=0.02, fy=36000)
+
+# Analysis
+responses = sdf.find_response(
+    time_steps, load_values, method="newmark_beta", acc_type="average"
+)
+print(responses)
+print(responses["displacement"][10])  # result is 0.03606328101158249

examples/sdf/eg_res_spec_el_centro.py

@@ -0,0 +1,19 @@
+# Section 6.4 (& Figure 6.6.2); Chopra A. K., Dynamics of structure, 5th edn
+import numpy as np
+from structdyn.ground_motions.ground_motion import GroundMotion
+from structdyn.sdf.response_spectrum import ResponseSpectrum
+from structdyn.utils.helpers import elcentro_chopra
+
+# Define el centro ground motion from Chopra's book- Appendix 6
+elc = elcentro_chopra()
+gm = GroundMotion.from_arrays(elc["acc (g)"], 0.02)
+
+# Define the period range of interest
+periods = np.arange(0, 5.01, 0.1)
+
+# Create response spectrum object
+rs = ResponseSpectrum(periods, 0.02, gm)
+
+# Analysis
+spectra = rs.compute()
+print(spectra["Sd"][20])  # result is 0.1896749378231744

structdyn/__init__.py

@@ -0,0 +1,23 @@
+"""
+structdyn: A Python library for solving structural dynamics problems numerically.
+
+Author: Abinash Mandal
+
+"""
+
+from .sdf.sdf import SDF
+from .sdf.numerical_methods.central_difference import CentralDifference
+from .sdf.numerical_methods.interpolation import Interpolation
+from .sdf.numerical_methods.newmark_beta import NewmarkBeta
+from .utils.helpers import plot_displacement, elcentro_chopra
+from .sdf.analytical_methods.analytical_response import AnalyticalResponse
+
+__all__ = [
+    "SDF",
+    "CentralDifference",
+    "Interpolation",
+    "NewmarkBeta",
+    "AnalyticalResponse",
+    "plot_displacement",
+    "elcentro_chopra",
+]

structdyn/ground_motions/__init__.py

@@ -0,0 +1,3 @@
+from .ground_motion import GroundMotion
+
+__all__ = ["GroundMotion"]