Add Modular UMAT integration

.github/workflows/build.yml

@@ -54,8 +54,15 @@ jobs:
       - name: Editable install
         run: |
           uv venv
-          uv pip install scikit-build-core pybind11 numpy
-          uv pip install -e .[dev]
+          # scikit-build-core pinned <1: the 1.x editable install changed where
+          # the compiled extension lands, which breaks the Windows DLL-vendoring
+          # step below (green on 0.12.2, DLL load failure on 1.0.2). Lift the pin
+          # together with a rework of that step. --no-build-isolation so the
+          # pinned backend in the venv is the one that actually builds the
+          # editable wheel (isolation would re-resolve pyproject's
+          # scikit-build-core>=0.11 to 1.x).
+          uv pip install "scikit-build-core<1" pybind11 numpy
+          uv pip install -e .[dev] --no-build-isolation
 
       - name: Vendor conda BLAS/LAPACK runtime next to the extension (Windows)
         if: runner.os == 'Windows'

conda.recipe/meta.yaml

@@ -60,4 +60,4 @@ test:
   imports:
     - simcoon
     - simcoon._core
-
+    - simcoon.modular

examples/data/MODUL_path.txt

@@ -0,0 +1,68 @@
+#Initial_temperature
+293
+#Number_of_blocks
+2
+
+#Block
+1
+#Loading_type
+1
+#Control_type(NLGEOM)
+1
+#Repeat
+1
+#Steps
+1
+
+#Mode
+1
+#Dn_init 1.
+#Dn_mini 1.
+#Dn_inc 0.005
+#time
+1
+#prescribed_mechanical_state
+E 0.02
+S 0 S 0
+S 0 S 0 S 0
+#prescribed_temperature_state
+T 293
+
+#Block
+2
+#Loading_type
+1
+#Control_type(NLGEOM)
+1
+#Repeat
+2
+#Steps
+2
+
+#Mode
+1
+#Dn_init 1.
+#Dn_mini 1.
+#Dn_inc 0.01
+#time
+1
+#prescribed_mechanical_state
+E -0.02
+S 0 S 0
+S 0 S 0 S 0
+#prescribed_temperature_state
+T 293
+
+#Mode
+1
+#Dn_init 1.
+#Dn_mini 1.
+#Dn_inc 0.01
+#time
+1
+#prescribed_mechanical_state
+E 0.02
+S 0 S 0
+S 0 S 0 S 0
+#prescribed_temperature_state
+T 293

examples/data/SHEAR_path.txt

@@ -0,0 +1,29 @@
+#Initial_temperature
+293.5
+#Number_of_blocks
+1
+
+#Block
+1
+#Loading_type
+1
+#Control_type(NLGEOM)
+1
+#Repeat
+1
+#Steps
+1
+
+#Mode
+1
+#Dn_init 1.
+#Dn_mini 1.
+#Dn_inc 0.005
+#time
+1
+#prescribed_mechanical_state
+S 0
+E 0.03 S 0
+S 0 S 0 S 0
+#prescribed_temperature_state
+T 293.5

examples/mechanical/MODUL.py

@@ -0,0 +1,133 @@
+"""
+Modular UMAT Example — Composable Elasto-Plasticity
+====================================================
+
+Demonstrates the ``simcoon.modular`` high-level Python interface that composes a
+constitutive model declaratively and runs it through the standard ``sim.solver``
+load-path driver. The C++ ``ModularUMAT`` infrastructure (ElasticityModule,
+YieldCriterion, hardening, ...) is internal — the user only builds a
+``ModularMaterial`` and hands its ``.props`` / ``.nstatev`` to the ``"MODUL"``
+UMAT code registered in simcoon's UMAT table.
+
+The path file applies a monotonic tensile ramp to 2% strain, then two
+strain-controlled cycles between -2% and +2%, exposing isotropic-hardening
+growth and the initial yield plateau.
+"""
+
+import os
+import numpy as np
+import matplotlib.pyplot as plt
+import simcoon as sim
+from simcoon.modular import (
+    ModularMaterial,
+    IsotropicElasticity,
+    Plasticity,
+    VonMisesYield,
+    VoceHardening,
+)
+
+plt.rcParams["figure.figsize"] = (14, 6)
+
+###################################################################################
+# 1. Compose the constitutive model
+# ----------------------------------
+# Isotropic elasticity + von Mises yield + Voce isotropic hardening.
+# Parameters: E=210 GPa, nu=0.3, sigma_Y=300 MPa, Q=200 MPa, b=10.
+
+mat = ModularMaterial(
+    elasticity=IsotropicElasticity(C1=210000.0, C2=0.3, alpha=1.2e-5),
+    mechanisms=[
+        Plasticity(
+            sigma_Y=300.0,
+            yield_criterion=VonMisesYield(),
+            isotropic_hardening=VoceHardening(Q=200.0, b=10.0),
+        ),
+    ],
+)
+
+print(mat.summary())
+
+###################################################################################
+# 2. Run the solver
+# ------------------
+# ``mat.umat_name`` is ``"MODUL"``, the UMAT code registered at
+# ``umat_smart.cpp:316`` (id 200). ``mat.props`` serializes the composition into
+# the flat array that the C++ ``umat_modular`` deserializes.
+
+umat_name = mat.umat_name
+props = mat.props
+nstatev = mat.nstatev
+
+psi_rve = 0.0
+theta_rve = 0.0
+phi_rve = 0.0
+solver_type = 0
+corate_type = 1
+
+path_data = "../data"
+path_results = "results"
+pathfile = "MODUL_path.txt"
+outputfile = "results_MODUL.txt"
+
+os.makedirs(path_results, exist_ok=True)
+
+sim.solver(
+    umat_name,
+    props,
+    nstatev,
+    psi_rve,
+    theta_rve,
+    phi_rve,
+    solver_type,
+    corate_type,
+    path_data,
+    path_results,
+    pathfile,
+    outputfile,
+)
+
+###################################################################################
+# 3. Plot the stress-strain curve
+# --------------------------------
+
+outputfile_macro = os.path.join(path_results, "results_MODUL_global-0.txt")
+
+e11, e22, e33, e12, e13, e23, s11, s22, s33, s12, s13, s23 = np.loadtxt(
+    outputfile_macro,
+    usecols=(8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19),
+    unpack=True,
+)
+time, T, Q_out, r = np.loadtxt(outputfile_macro, usecols=(4, 5, 6, 7), unpack=True)
+Wm, Wm_r, Wm_ir, Wm_d = np.loadtxt(
+    outputfile_macro, usecols=(20, 21, 22, 23), unpack=True
+)
+
+fig = plt.figure()
+
+# Stress-strain curve
+ax1 = fig.add_subplot(1, 2, 1)
+plt.grid(True)
+plt.tick_params(axis="both", which="major", labelsize=13)
+plt.xlabel(r"Strain $\varepsilon_{11}$", size=14)
+plt.ylabel(r"Stress $\sigma_{11}$ (MPa)", size=14)
+plt.plot(e11, s11, c="royalblue", lw=1.5, label="MODUL: iso-elastic + VM + Voce")
+plt.axhline(y=300.0, color="0.6", linestyle="--", lw=0.8, label=r"initial $\sigma_Y$")
+plt.axhline(y=-300.0, color="0.6", linestyle="--", lw=0.8)
+plt.legend(loc="best")
+plt.title("Stress-strain response")
+
+# Work terms vs time
+ax2 = fig.add_subplot(1, 2, 2)
+plt.grid(True)
+plt.tick_params(axis="both", which="major", labelsize=13)
+plt.xlabel("time (s)", size=14)
+plt.ylabel("Work (MPa)", size=14)
+plt.plot(time, Wm, c="black", label=r"$W_m$ (total)")
+plt.plot(time, Wm_r, c="green", label=r"$W_m^r$ (recoverable)")
+plt.plot(time, Wm_ir, c="blue", label=r"$W_m^{ir}$ (irreversible)")
+plt.legend(loc="best")
+plt.title("Energy decomposition")
+
+plt.tight_layout()
+plt.savefig("MODUL_stress_strain.png", dpi=120)
+plt.show()