Skip to content

Cure-dependent inclusion fractions and chemical shrinkage as an eigenstrain - #137

Merged
petlenz merged 44 commits into
mainfrom
feature/curing-eigenstrain
Oct 27, 2026
Merged

Cure-dependent inclusion fractions and chemical shrinkage as an eigenstrain#137
petlenz merged 44 commits into
mainfrom
feature/curing-eigenstrain

Conversation

@petlenz

@petlenz petlenz commented Oct 16, 2026 •
edited
Loading

Copy link
Copy Markdown
Member

Closes #115. Part of the uvwmaterial port (#94); stacked on #136.

uvwmaterial's eigenstrain_curing_material_base and degree_of_cure_dependent_material become two wirings rather than classes:

  • Chemical shrinkage ε_ch = β ξ I is thermal_eigenstrain with the degree of cure as its scalar source ("temperature_source": "cure", "alpha": beta, reference 0); the integrator's state is a history like a temperature source's. Nothing new to port; the test pins it.
  • Cure-dependent phase fractions: a phase reads its fraction from any producer (fraction_source / fraction_property, e.g. "cure", "state" of a scalar_stepper or a weighted_sum output); its volume_fraction output follows that property every update, scalar_complement_weight produces the matrix fraction 1 − c, and mean_field_tangent reads both. Cure-dependent moduli already go through weighted_sum.

Tests (test_curing_eigenstrain, 3): the phase fraction follows its source and the Mori–Tanaka bulk modulus follows the analytic formula at each value; constrained shrinkage gives −3Kβξ; fraction_source from JSON.

No report yet (time-driven point, #126).

924/924 tests pass (GCC 14, reports on).

This was referenced Oct 16, 2026
Open
Merged
@petlenz
petlenz force-pushed the feature/kamal-sourour-rate branch from f56825b to 515796b Compare October 16, 2026 23:43
@petlenz
petlenz force-pushed the feature/curing-eigenstrain branch from c660cef to 4a060f0 Compare October 16, 2026 23:43
@petlenz

petlenz commented Oct 17, 2026

Copy link
Copy Markdown
Member Author

Review findings (design, tmech expressions instead of component loops, modern C++, comments, tests, correctness). Severity: must = wrong result, hang, build break or unreachable feature; should = design, duplication, accuracy; nit = style.

curing fractions / chemical shrinkage — ready with one should

  • [should] volume_fraction stays required with fraction_source (dummy in JSON, can disagree until first update) → default 0 / exactly one of the two.
  • [nit] sourced fraction never range-checked; pre-existing direction loop.

This was referenced Oct 17, 2026
Merged
Merged
@petlenz
petlenz force-pushed the feature/kamal-sourour-rate branch 2 times, most recently from e4f80c1 to 316263a Compare October 17, 2026 14:31
@petlenz
petlenz force-pushed the feature/curing-eigenstrain branch from 4a060f0 to bd0004d Compare October 17, 2026 14:40
@petlenz

petlenz commented Oct 17, 2026

Copy link
Copy Markdown
Member Author

Rebased on the Eshelby-node refactor (#128/#131): the curing fixture and the JSON test wire an eshelby_sphere node to the matrix and the inclusion names it; fraction_source is merged into the new inclusion. The volume_fraction-beside-fraction_source finding is still open.

@petlenz
petlenz force-pushed the feature/kamal-sourour-rate branch from 316263a to 9928a69 Compare October 17, 2026 15:57
@petlenz
petlenz force-pushed the feature/curing-eigenstrain branch from bd0004d to 6eff511 Compare October 17, 2026 16:05
@petlenz
petlenz force-pushed the feature/kamal-sourour-rate branch from 9928a69 to 971cf5a Compare October 17, 2026 17:04
@petlenz
petlenz force-pushed the feature/curing-eigenstrain branch from 6eff511 to 30fbcba Compare October 17, 2026 17:12
@petlenz

petlenz commented Oct 17, 2026

Copy link
Copy Markdown
Member Author

Producer/consumer pass: the curing inclusion gives its fraction only through fraction_source: "cure::state" (no dummy volume_fraction any more), the shrinkage coefficient comes through alpha_source.

@petlenz
petlenz force-pushed the feature/kamal-sourour-rate branch from 971cf5a to 6966212 Compare October 17, 2026 21:36
@petlenz
petlenz force-pushed the feature/curing-eigenstrain branch from 30fbcba to bc568df Compare October 17, 2026 21:36
@petlenz
petlenz force-pushed the feature/kamal-sourour-rate branch from 6966212 to 5463017 Compare October 17, 2026 22:46
@petlenz
petlenz force-pushed the feature/curing-eigenstrain branch from bc568df to f964442 Compare October 17, 2026 22:46
@petlenz
petlenz force-pushed the feature/kamal-sourour-rate branch from 5463017 to 8558f90 Compare October 18, 2026 12:32
@petlenz
petlenz force-pushed the feature/curing-eigenstrain branch from f964442 to a4bd728 Compare October 18, 2026 12:51
@petlenz
petlenz force-pushed the feature/kamal-sourour-rate branch from 8558f90 to e2a4434 Compare October 18, 2026 13:38
@petlenz
petlenz force-pushed the feature/curing-eigenstrain branch from a4bd728 to 4c1aac1 Compare October 18, 2026 13:38
@petlenz
petlenz force-pushed the feature/kamal-sourour-rate branch from e2a4434 to 39adb32 Compare October 18, 2026 18:55
@petlenz
petlenz force-pushed the feature/curing-eigenstrain branch from 4c1aac1 to 93a893f Compare October 18, 2026 18:56
@petlenz
petlenz force-pushed the feature/kamal-sourour-rate branch from 39adb32 to 2afbfb9 Compare October 18, 2026 20:03
@petlenz
petlenz force-pushed the feature/curing-eigenstrain branch from 93a893f to 8a07c8c Compare October 18, 2026 20:03
@petlenz
petlenz force-pushed the feature/kamal-sourour-rate branch from 2afbfb9 to e7d2ffe Compare October 18, 2026 20:22
@petlenz
petlenz force-pushed the feature/curing-eigenstrain branch from 8a07c8c to e1819c4 Compare October 18, 2026 20:22
petlenz added 11 commits October 19, 2026 05:49
@petlenz
Isotropisation, Lebedev quadrature and depolarisation tensors
207536a 
core/isotropize.h: K = Pvol::C / Dim, G = Pdev::C / (2 (Dim(Dim+1)/2 - 1)),
and the isotropic part built from them. core/lebedev.h with generated
tables (tools/extract_lebedev.py from uvwmaterial's rule, precisions 3 to
23, up to 194 points): spherical means and integrals of tensor-valued
functions of the direction. eshelby::depolarisation for ellipsoids via
Carlson's R_J, the elliptic cylinder in closed form, and the cylinder along
a direction.
@petlenz
rule_of_mixture: Voigt, Reuss and Hill effective stiffness
6520437 
Phases are materials providing a tangent, with volume fractions; outputs
the effective tangent, the stress and the strain concentration tensors
(sum c_i A_i = Isym). Report with the closed forms of the three kernels
and their ordering. Report generator: derived values resolve in
dependency order, independent of the key order.
@petlenz
mean_field_tangent: dilute, Mori-Tanaka and self-consistent as a tang…
655efa4 
…ent producer over phases

The composite consumes the phases' produced tangents, fractions and Eshelby
tensors and produces the effective tangent and the concentrations; the stress
is a linear_stress on the tangent. It reaches into no other node: no
material_ref, no eshelby(nu), no isotropisation of the matrix.

Self-consistent is a graph cycle closed by one local edge: isotropic_moduli
on the composite tangent produces nu of the effective medium, the Eshelby
nodes read it with nu_edge = local (from nu_initial on the first pass), and
one update() is one relaxed Picard step with "residual" as its error. The
matrix is a phase with its own Eshelby node, so its shape is no longer a
hard-coded sphere. material_point_evaluator closes the fixed point per host
call (residual_property, residual_tolerance, max_passes; the report point
names them). isotropic_moduli with a non-positive tolerance projects without
the isotropy check, which the effective medium needs.
@petlenz
short_fibre_composite: Mori-Tanaka over the fibre orientation distrib…
618831b 
…ution

Spheroidal fibres in an isotropic matrix with the second-order orientation
tensor: Mori-Tanaka with the orientation-averaged dilute concentration, or
the averaged unidirectional stiffness. Eshelby tensors per direction are
cached for the matrix Poisson ratio. Report with aligned, injection-moulded
and random states against the unidirectional composite and isotropy.
@petlenz
short_fibre_composite report: the aligned check as a ratio against zero
369385f
@petlenz
Orientation averaging with the angular central Gaussian
d855f53 
core/orientation.h: fit_acg solves a_i = R_D(b_j, b_k, b_i) / 3 with
det b = 1 by Newton in the principal frame of the orientation tensor (any
A, not only diagonal), the aligned case as a single direction, planar
cases by an eigenvalue floor; nodes maps the Lebedev rule through B^-1/2;
average integrates a function of the direction.
@petlenz
short_fibre_composite: matrix and fibre as phases, a tangent producer…
02524f8 
… in the fibre frame

Both phases come in as produced properties (tangent, fraction, the fibre's
Eshelby tensor along x1); the composite consumes them, forms the dilute
concentration once in the fibre frame, rotates onto the quadrature directions
with otimesu(R, R) fixed per direction, and produces the tangent and the two
concentrations. The stress is a linear_stress on the tangent. No material_ref.
@petlenz
incremental_linear_elasticity: hypoelastic update with the stress as …
108389e 
…history

sigma_{n+1} = sigma_n + C : (eps_{n+1} - eps_n) with the tangent read from
another material, for path-dependent tangents. Report: equals the total
form on every path, a closed path returns to zero, exact for any step.
Drivers: the convergence test also scales with |C| |d eps| (|A| |dF|), so
a hypoelastic stress crossing zero under stress control converges.
@petlenz
thermal_eigenstrain and an eigenstrain_source for linear_elasticity
25909f3 
eps_th = alpha (T - T_ref) I from a scalar temperature source; linear
elasticity subtracts an optional eigenstrain: sigma = C : (eps - eps_eig).
@petlenz
kamal_sourour_rate: nth-order plus autocatalytic cure kinetics
9037f62 
dz/dt = (k1 + k2 z^m) (1 - z)^n with Arrhenius rate constants, as a rate
law for rk_integrator next to curing_rate (the autocatalytic term alone;
nth order is curing_rate with m = 0).
@petlenz
Cure-dependent inclusion fractions; chemical shrinkage as an eigenstrain
9af23f2 
inclusion takes a fraction_source ("material::property") so the volume
fraction can follow the degree of cure; the chemical shrinkage
beta xi I is thermal_eigenstrain with the cure state as its source.
@petlenz
petlenz force-pushed the feature/kamal-sourour-rate branch from e7d2ffe to 9037f62 Compare October 19, 2026 04:07
@petlenz
petlenz force-pushed the feature/curing-eigenstrain branch from e1819c4 to 9af23f2 Compare October 19, 2026 04:07
@petlenz
Merge main into feature/material-reports
ea56e0d
@petlenz
Merge main into feature/material-reports
ea3a120 
Main has the reviewed split of this branch (#199-#203), so main's version is
taken everywhere; the branch no longer differs from main.
@petlenz
Merge branch 'feature/material-reports' into feature/elliptic-integra…
fe57dcf 
…ls-and-constants

# Conflicts:
#	tests/CMakeLists.txt
@petlenz
petlenz changed the base branch from feature/kamal-sourour-rate to main October 24, 2026 21:16
petlenz added 11 commits October 24, 2026 23:22
@petlenz
Merge branch 'feature/elliptic-integrals-and-constants' into feature/…
b3414d0 
…eshelby-tensors

# Conflicts:
#	tests/CMakeLists.txt
@petlenz
Merge branch 'feature/eshelby-tensors' into feature/isotropize-lebede…
ed3f23e 
…v-depolarisation

# Conflicts:
#	tests/CMakeLists.txt
@petlenz
Merge branch 'feature/isotropize-lebedev-depolarisation' into feature…
64b5d15 
…/rule-of-mixture

# Conflicts:
#	tests/CMakeLists.txt
@petlenz
Merge branch 'feature/rule-of-mixture' into feature/mean-field-composite
9a5c937 
# Conflicts:
#	include/numsim-materials/report/report_json.h
#	tests/CMakeLists.txt
@petlenz
Drivers pass the residual settings through to the material point (#210)
b2ed115 
residual_source, residual_tolerance and max_passes are driver parameters
(material_point's defaults) and keys of the report's "point"; the mean-field
tangent report uses them for the self-consistent fixed point.
@petlenz
Merge branch 'feature/mean-field-composite' into feature/orientation-…
3d07c94 
…averaging

# Conflicts:
#	tests/CMakeLists.txt
@petlenz
Merge branch 'feature/orientation-averaging' into feature/short-fibre…
47eb27b 
…-composite

# Conflicts:
#	tests/CMakeLists.txt
@petlenz
Merge branch 'feature/short-fibre-composite' into feature/incremental…
5531189 
…-elasticity

# Conflicts:
#	include/numsim-materials/umat/deformation_gradient_driver.h
#	include/numsim-materials/umat/material_point_driver.h
#	tests/CMakeLists.txt
@petlenz
Merge branch 'feature/incremental-elasticity' into feature/thermal-ei…
427897e 
…genstrain

# Conflicts:
#	tests/CMakeLists.txt
@petlenz
Merge branch 'feature/thermal-eigenstrain' into feature/kamal-sourour…
e6f9623 
…-rate

# Conflicts:
#	tests/CMakeLists.txt
@petlenz
Merge branch 'feature/kamal-sourour-rate' into feature/curing-eigenst…
b24212a 
…rain

# Conflicts:
#	tests/CMakeLists.txt
@petlenz
petlenz marked this pull request as ready for review October 27, 2026 17:50
@petlenz
petlenz merged commit 6090947 into main Oct 27, 2026
2 checks passed
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Reviewers

No reviews

Assignees

No one assigned

Labels

None yet

Projects

None yet

Milestone

No milestone

Development

Successfully merging this pull request may close these issues.

curing_eigenstrain and degree-of-cure dependent properties

1 participant

@petlenz