vCANN: Material Constitutive Modeling and Finite Element Simulation Automation Based on Physics-Constrained Neural Networks

This thread introduces vCANN (viscoelastic Constitutive Artificial Neural Networks), a thermodynamically consistent physics-informed neural network framework. It enables end-to-end automated modeling from experimental data to finite element simulation, addressing core challenges in material science and engineering. Key highlights include embedding physical constraints (thermodynamic consistency, objectivity) to ensure model validity, and direct export to FEM software like Abaqus.

Traditional material constitutive modeling relies on empirical mathematical formulas derived by physicists, then fitting parameters with experimental data. However, complex viscoelastic materials (e.g., polymers, biological tissues) exhibit high nonlinearity and time dependence, making traditional models inaccurate. Moreover, integrating models into FEM software requires extensive manual coding, creating a "gap" between experiments and simulation that hinders efficiency in new material development and engineering design.

vCANN uses physics-informed neural networks (PINNs) to learn constitutive relations from experimental data while embedding thermodynamic constraints:

1. Energy conservation (1st law of thermodynamics)
2. Entropy increase principle (2nd law, non-negative dissipation)
3. Objectivity (material frame indifference)

Network Structure: Inputs include deformation state (gradient tensor), deformation history (internal variables), and environment parameters (temperature). Outputs are Cauchy stress tensors, with architecture designed to satisfy constraints automatically.

End-to-End Flow:

1. Import experimental data (e.g., uniaxial/biaxial tension, shear tests; example uses VHB4905 material data).
2. Train with Keras/TensorFlow: loss function combines data fitting error and physical constraint penalties.
3. Export to FEM user material subroutines (UMAT/VUMAT) for direct use in Abaqus.

Using VHB4905 viscoelastic material data (different temperatures, loading-unloading cycles), vCANN demonstrates:

• Accurate stress prediction: High agreement with experimental stress-strain curves.
• Hysteresis capture: Correctly reproduces loading-unloading hysteresis loops of viscoelastic materials.
• Temperature dependence: Predicts material behavior across different temperatures via learned internal variables.
• Generalization: Reasonable predictions on untrained loading paths.
• Physical validity: Avoids non-physical oscillations in sparse data regions due to embedded constraints.

vCANN brings significant value to material engineering and simulation:

1. Accelerate new material development: Reduces constitutive modeling time from months/years to days (standard experiments → training → FEM-ready model).
2. Lower simulation barriers: Enables engineers without deep constitutive theory background to build high-quality models.
3. Multi-scale analysis: Extensible to integrate micro-scale data (e.g., DIC measurements) with macro FEM simulations.
4. Uncertainty quantification: Bayesian neural network extensions provide prediction uncertainty estimates for safety assessment.

Current Limitations:

• Data requirements: Needs sufficient experimental data covering operating condition ranges; extreme conditions (e.g., impact) data are hard to obtain.
• Computational cost: Higher than traditional analytical models in FEM simulations.
• Interpretability: Lack of intuitive explanation for internal network learning mechanisms.

Future Directions:

• Combine active learning to reduce data needs.
• Develop more efficient network architectures to lower online computation cost.
• Integrate with molecular dynamics simulations for seamless multi-scale modeling.
• Extend to complex material behaviors (damage, fracture, anisotropy).

Dependencies:

• TensorFlow/Keras (neural network training/inference)
• NumPy/SciPy (numerical computation/optimization)
• scikit-learn (data preprocessing/validation)
• matplotlib/plotly (visualization)
• kormos (constitutive model optimization library)

Code Structure:

• Data loading/preprocessing modules
• Neural network architecture definition
• Thermodynamic constraint implementation
• Training scripts & hyperparameter config
• FEM model export tools (Abaqus UMAT generator)
• Validation/visualization tools