Deformation Gradient Research Articles

Restraining geometrically-necessary dislocations to the active slip systems in a crystal plasticity-based finite element framework

Strain gradients have been cast in the form of geometrically-necessary dislocations (GND) to relate the length-scale dependence of strength and to determine potential sites for failure initiation. The literature contains various different incompatibility measures, the main ones being: the total form (∇×Fp), the rate form for large displacements (∇×γ̇anaFp), and the slip gradient form (∇γ̇a). Here, these different approaches are compared rigorously for the first time. Obtaining GND densities when using the total form is a rank-deficit linear problem, solved by singular value decomposition (SVD) known as the Least Squares Minimization (L2 method). Alternative methods for finding GND densities such as Karush–Kuhn–Tucker (KKT) optimization are also investigated. Both L2 and KKT methods predict unrealistic GND densities on inactive slip systems leading to excessive strain hardening; even for a single crystal single slip case. Therefore, the restriction of GNDs to the active slip systems by using a threshold based on the total slip is proposed. This restriction reveals relatively consistent results for various single crystal single slip cases including: simple shear, uniaxial tension, and four-point bending. In addition, the small numerical differences in the slip leads to large discrepancies in the flow stress due to error accumulation, even for strain-gradient-free uniaxial tension, hence a threshold for the GND density increment (2×102 m−2) is used in all models to avoid formation of erroneous GND densities. Finally, the proposed method is applied to the evolution of the GND density for a grain inside a polycrystal aggregate that posses a complex stress state. The lowest incompatibility error is obtained by both of the total forms that use the curl of the plastic deformation gradient with the active slip system restriction suggesting them to be the most reliable GND measures.

Finite strain modelling of Shape Memory Alloys in the logarithmic strain space: A comparative study with other finite strain approaches

Finite strain models for Shape Memory Alloys (SMAs) offer significant advantages over small strain models in that they excel in their ability to characterize SMA behaviour under large strain and finite rotations while eliminating the inherent assumptions associated with the small strain framework. The paper at hand presents a numerical model for SMAs at finite strain, employing the logarithmic strain measure to obtain an additive split of the total strain into its elastic and inelastic (transformation) parts. Notably, the model maintains the thermodynamic consistency validated using the Clausius–Duhem inequality. Additionally, we compare this model with two other finite strain SMA models employing logarithmic strain measures: a hypoelastic model, which splits the rate of deformation tensor additively, and a hyperelastic model, which decomposes the deformation gradient multiplicatively. The numerical formulation for the three finite strain models is specifically tailored to address superelasticity in SMAs. It is derived from the established small strain model and seamlessly extended into the domain of finite strain framework. Several numerical examples are carried out to showcase the logarithmic strain space approach applied to SMAs and to compare its response with that of the other aforementioned finite strain models in SMAs. The meticulous examination reveals that the logarithmic strain space approach aligns well with the existing finite strain models, demonstrating its effectiveness and compatibility.

An SPH formulation for general plate and shell structures with finite deformation and large rotation

In this paper, we propose a reduced-dimensional smoothed particle hydrodynamics (SPH) formulation for quasi-static and dynamic analyses of plate and shell structures undergoing finite deformation and large rotation. By exploiting Uflyand–Mindlin plate theory, the present surface-particle formulation is able to resolve the thin structures by using only one layer of particles at the mid-surface. To resolve the geometric non-linearity and capture finite deformation and large rotation, two reduced-dimensional linear-reproducing correction matrices are introduced, and weighted non-singularity conversions between the rotation angle and pseudo normal are formulated. A new non-isotropic Kelvin-Voigt damping is proposed especially for the both thin and moderately thick plate and shell structures to increase the numerical stability. In addition, a shear-scaled momentum-conserving hourglass control algorithm with an adaptive limiter is introduced to suppress the mismatches between the particle position and pseudo normal and those estimated with the deformation gradient. A comprehensive set of test problems, for which the analytical or numerical results from literature or those of the volume-particle SPH model are available for quantitative and qualitative comparison, are examined to demonstrate the accuracy and stability of the present method.

A Multiphysics Thermoelastoviscoplastic Damage Internal State Variable Constitutive Model including Magnetism.

We present a macroscale constitutive model that couples magnetism with thermal, elastic, plastic, and damage effects in an Internal State Variable (ISV) theory. Previous constitutive models did not include an interdependence between the internal magnetic (magnetostriction and magnetic flux) and mechanical fields. Although constitutive models explaining the mechanisms behind mechanical deformations caused by magnetization changes have been presented in the literature, they mainly focus on nanoscale structure-property relations. A fully coupled multiphysics macroscale ISV model presented herein admits lower length scale information from the nanoscale and microscale descriptions of the multiphysics behavior, thus capturing the effects of magnetic field forces with isotropic and anisotropic magnetization terms and moments under thermomechanical deformations. For the first time, this ISV modeling framework internally coheres to the kinematic, thermodynamic, and kinetic relationships of deformation using the evolving ISV histories. For the kinematics, a multiplicative decomposition of deformation gradient is employed including a magnetization term; hence, the Jacobian represents the conservation of mass and conservation of momentum including magnetism. The first and second laws of thermodynamics are used to constrain the appropriate constitutive relations through the Clausius-Duhem inequality. The kinetic framework employs a stress-strain relationship with a flow rule that couples the thermal, mechanical, and magnetic terms. Experimental data from the literature for three different materials (iron, nickel, and cobalt) are used to compare with the model's results showing good correlations.

A peridynamic-informed deep learning model for brittle damage prediction

Physics-informed Neural Network (PINN) has been introduced recently to predict and understand complex physical phenomena by directly incorporating feedback from governing equations. While PINNs offer remarkable capabilities, challenges arise when addressing discontinuities and heterogeneity, potentially compromising accuracy and reliability in prediction outcomes. In this study, a novel approach that combines the principles of peridynamic (PD) theory with PINN is presented to predict quasi-static damage and crack propagation in brittle materials. To achieve high prediction accuracy and convergence rate, the linearized PD governing equation is enforced in the PINN’s residual-based loss function. The proposed PD-INN is able to learn and capture complicated displacement patterns associated with different geometrical parameters, such as pre-crack position and length. The computational method presented in this study can be implemented to any PD constitutive model like Bond-Based PD or State-Based PD as well as any horizon factor. Several enhancements like cyclical annealing schedule and deformation gradient aware optimization technique are proposed to ensure the model would not get stuck in its trivial solution. The model’s performance assessment is conducted by monitoring the behavior of loss function throughout the training process. The results obtained through PD-INN demonstrate a satisfactory level of agreement in 2D and 3D when compared with the ground truth deformation from high-fidelity techniques like PD direct numerical method and Extended finite element method (X-FEM). Our results show the ability of the nonlocal PD-INN to predict damage and crack propagation accurately and efficiently.

Membrane Tilt Drives Phase Separation of Adhesion Receptors.

Cell adhesion receptors are transmembrane proteins that bind cells to their environment. These proteins typically cluster into disk-shaped or linear structures. Here, we show that such clustering patterns spontaneously emerge when the receptor senses the membrane deformation gradient, for example, by reaching a lower-energy conformation when the membrane is tilted relative to the underlying binding substrate. Increasing the strength of the membrane gradient-sensing mechanism first yields isolated disk-shaped clusters and then long linear structures. Our theory is coherent with experimental estimates of the parameters, suggesting that a tilt-induced clustering mechanism is relevant in the context of cell adhesion.

Folding of a growing hyperelastic sheet in a viscous fluid

Understanding the growth-induced deformation of soft materials in viscous environments is important for a variety of problems in nature and engineering. Here, we focus on the fluid–structure interaction of a hyperelastic sheet growing in an incompressible Newtonian fluid in the Stokes flow regime. We develop a computational framework for simulating this problem, where the isogeometric boundary integral method is used with the Kirchhoff–Love shell formulation and elastic–plastic decomposition of the deformation gradient tensor. We quantify the relative effects of the growth rate, the sheet bending rigidity, and the fluid viscosity on the fold formation and development of the growing sheet. Our results suggest that the viscous resistance to in-plane deformation promotes fold formation, whereas the viscous resistance to out-of-plane deformation suppresses fold development. We also investigate the effects of the thickness and aspect ratio of the rectangular sheet. Finally, we compare the growth- and prestrain-induced deformations to find a common behavior of sheets under viscous environments.

Micro-mechanisms of anisotropic deformation in the presence of notch in commercially pure Titanium: An in-situ study with CPFEM simulations

This study investigates the impact of notch severity and initial texture on the micro-mechanisms during low strain deformation in commercially pure Titanium, exploiting in-situ electron back scatter diffraction (EBSD) experiments and crystal plasticity finite element model (CPFEM). In-situ tests were performed at different strain steps of un-notched and notched samples of transverse direction (TD) and rolling direction (RD) orientations. CPFEM, based on initial EBSD microstructures, predicted profuse prismatic slip traces and early activation of prismatic slip in notched sample, with RD orientation exhibiting higher activity. Further, CPFEM results revealed early activation of high CRSS slip systems as well as evidence of early twin activity at notch tip due to severely localized plastic deformation and steep strain gradient, as observed by GND maps causing higher stress at grain boundaries. At the notch tip, digital images correlation (DIC) at microscale indicated strain localization at 45∘ and 90∘ to the tensile axis for TD and RD orientations, respectively. Furthermore, 2.5D and 3D CPFEM confirmed distinct strain patterns at notch tip: TD orientation exhibited combined basal and prismatic slip influences, while RD orientation displayed dominant prismatic slip systems at localized strain. The model also successfully predicted anisotropic surface roughness, contributing to early necking in RD orientation.

Connecting weakly nonlinear elasticity theories of isotropic hyperelastic materials

Soft materials exhibit significant nonlinear geometric deformations and stress–strain relationships under external forces. This paper explores weakly nonlinear elasticity theories, including Landau’s and Murnaghan’s formulations, advancing understanding beyond linear elasticity. We establish connections between these methods and extend strain-energy functions to the third and fourth orders in power of ε , where ε = H · H and 0 < ε ≤ 1 , and H is the perturbation to the deformation gradient tensor F = I + H . Furthermore, we address simplified strain-energy functions applicable to incompressible materials. Through this work, we contribute to a comprehensive understanding of nonlinear elasticity and its relationship to weakly nonlinear elasticity, facilitating the study of moderate deformations in soft material behavior and its practical applications.

A first‐order hyperbolic arbitrary Lagrangian Eulerian conservation formulation for non‐linear solid dynamics

SummaryThe paper introduces a computational framework using a novel Arbitrary Lagrangian Eulerian (ALE) formalism in the form of a system of first‐order conservation laws. In addition to the usual material and spatial configurations, an additional referential (intrinsic) configuration is introduced in order to disassociate material particles from mesh positions. Using isothermal hyperelasticity as a starting point, mass, linear momentum and total energy conservation equations are written and solved with respect to the reference configuration. In addition, with the purpose of guaranteeing equal order of convergence of strains/stresses and velocities/displacements, the computation of the standard deformation gradient tensor (measured from material to spatial configuration) is obtained via its multiplicative decomposition into two auxiliary deformation gradient tensors, both computed via additional first‐order conservation laws. Crucially, the new ALE conservative formulation will be shown to degenerate elegantly into alternative mixed systems of conservation laws such as Total Lagrangian, Eulerian and Updated Reference Lagrangian. Hyperbolicity of the system of conservation laws will be shown and the accurate wave speed bounds will be presented, the latter critical to ensure stability of explicit time integrators. For spatial discretisation, a vertex‐based Finite Volume method is employed and suitably adapted. To guarantee stability from both the continuum and the semi‐discretisation standpoints, an appropriate numerical interface flux (by means of the Rankine–Hugoniot jump conditions) is carefully designed and presented. Stability is demonstrated via the use of the time variation of the Hamiltonian of the system, seeking to ensure the positive production of numerical entropy. A range of three dimensional benchmark problems will be presented in order to demonstrate the robustness and reliability of the framework. Examples will be restricted to the case of isothermal reversible elasticity to demonstrate the potential of the new formulation.

Localization patterns emerging in CPTu tests in a saturated natural clay soil

In this work, the mechanics of cone penetration during CPTu tests in Osaka clay, a natural structured silty clay, was investigated with the Particle Finite Element Method (PFEM). To describe the behavior of the soil, the FD_MILAN model, recently proposed by the Authors, was adopted. The model, based on the multiplicative decomposition of the deformation gradient, incorporates a scalar internal variable, the bond strength Pt, to quantify the effects of structure. The macroscopic description of mechanical destructuration effects is provided by the hardening law of Pt, according to which the bond strength decreases monotonically with accumulated plastic deformations. The brittle behavior resulting from the softening process associated to destructuration make the soil quite susceptible to the spontaneous development of strain localization in the form of shear bands. In order to deal with strain localization phenomena, the model is equipped with a non-local version of the hardening laws, incorporating a material constant – the characteristic length ℓc – which provides the material with an internal length scale. The objective of this work was to get some insight on the effects of the characteristic length and of the initial degree of structure on: the kinematics of the soil deformation around the advancing cone tip; the evolution of soil structure with accumulated plastic deformations; the excess pore pressure field generated in the soil as a result of the hydro-mechanical coupling between the solid skeleton and the pore water, and the net cone resistance measured at the base of the cone. The results of the PFEM simulations show that the deformations around the piezocone are strongly affected by the characteristic length. For heavily structured soils, when ℓc is relatively small with respect to cone radius, the accumulated plastic deviatoric deformation field is characterized by clearly visible shear bands while, as ℓc increases, the detection of localized deformation regions becomes more difficult, if not impossible. This depends on the fact that, for large values of ℓc, the shear band width may be of comparable size to the cone radius. In all cases, the soil around the advancing piezocone is subjected to a very strong destructuration process, which leads to the complete loss of bond strength in a large region around the cone tip and shaft. As a consequence, the use of conventional Nc values from the literature in the interpretation of the CPTu test performed in heavily structured soils could provide undrained strength values closer to the ultimate undrained shear strength. At the same time, the use of empirical correlation to estimate the yield stress in oedometric conditions would lead to a significant underestimation of the overconsolidation ratio.

Phase field modeling of hyperelastic material interfaces –Theory, implementation and application to phase transformations

Interface mechanics can significantly govern the evolution of multiple phases on smaller scales, e.g., determining the properties of TWIP- and TRIP-steels, geopolymers or Li-ion batteries. The present contribution is specifically centered around the influence of interface elasticity on mechanically induced phase transformations. A geometrically exact finite element framework is developed for this purpose based on incremental energy minimization. It is set up in three-dimensional space and comprises an Allen-Cahn type phase field formulation based on a Modica-Mortola double-well functional, as well as a Ginzburg-Landau type viscosity formulation. The phase field approximation of the interface deformation gradient is derived by means of homogenization theory. A monolithic solver is employed for the resulting set of non-linear coupled equations. The example of coalescing spheres shows that larger interface energies can amplify and accelerate coalescence. Lower interface energies, in contrast, allow for distinct spreading of the energetically favorable bulk phase. Deformation-dependent interface energies particularly add stress peaks to this process that are localized at the phase transition zone. The relaxation of a stepped interface shape showed that deformation-dependent interface energies favor a plane partition of the bulk phases. Yet, they also delay breakup events compared to simple constant interface energies. This peculiarity can be beneficial for calibration purposes. Strain dependence moreover causes an additional elastic stiffness in the direction of the interface tangent plane and constitutes thus another origin for anisotropy of the overall structure.

How does the uncut chip thickness affect the deformation states within the primary shear zone during metal cutting?

The deformation states within the primary shear zone (PSZ) significantly affect material removal during machining. Uncut chip thickness (UCT) is an important factor that influences the material deformation states. However, the specific mechanism by which UCT influences the deformation states within PSZ remains unknown. This study aims to investigate the relationship between the deformation states in PSZ and UCTs via in-situ measurement and microscopic characterization techniques. Using the digital image correlation (DIC) technique, strain and strain rate distributions were derived to reveal the discrepant deformation in PSZ with increasing UCT. Furthermore, velocity vector fields and Electron Back-Scattered Diffraction (EBSD) characterizations were employed to examine the heterogeneity of deformation modes. To determine the specific deformation information, a deformation extraction framework based on the deformation gradient tensor theory was developed. Thus, strong and weak shear modes within PSZ were revealed based on the full-field deformation information of compression and extension. As the UCT increased, the transition of deformation states from a strong shear state to a hybrid shear state was determined. This work presents a new understanding of the deformation mechanism within PSZ in a ductile material of pure iron. A critical UCT was proposed to guide the cutting process to avoid inefficient weak shear mode.

Evaluation method of the synergistic service performance between ultrahigh-performance concrete and normal strength concrete based on the parameter of environmental equivalent force ratio

This study investigated the synergistic service performances of ultrahigh-performance concrete (UHPC) and normal-strength concrete (NSC). Restrained specimens used in this study were designed to assess the non-equivalent deformation of a substrate and UHPC under specific service conditions and shrinkage. Various UHPC repair materials characterized by different shrinkage deformation gradients were developed using expansive agents and lightweight aggregates. The study explored the effects of shrinkage on the synergistic service properties of UHPC-NSC. In addition, the bond fatigue performance of UHPC-NSC was investigated, and a bond fatigue-life equation was formulated. Under static loads and a stress ratio of 0.9, specimens exhibited both substrate and bond failures. By contrast, specimens under stress ratios of 0.8 and 0.7 displayed substrate failure. The synergistic effect of a 10 % dosage of expansive agent and a 15 % dosage of lightweight aggregate resulted in a 67.2 % reduction in UHPC shrinkage after 28 d. Internal stress derived from the non-equivalent deformation of the substrate and UHPC due to shrinkage was quantified. A reduction in UHPC shrinkage led to a 33.0–67.3 % decrease in tensile stress for repair materials. Internal stresses in the UHPC-NSC under alternating environments, including wet–dry, hot–cool, and freeze–thaw cycles, were calculated, where the ratio of internal stress to interfacial strength was defined as the interfacial stress-intensity ratio. This study presents an evaluation method for the synergistic service performance of UHPC-NSC based on an environmental equivalent-force ratio, thereby providing guidance for practical concrete-repair service-life design.

The Mechanics of Elastomeric Sheet Reinforced With Bidirectional Fiber Mesh Subjected to Flexure on Boundaries

Abstract We investigate the concurrent three-dimensional deformations of fiber-reinforced composite sheets subjected to out-of-plane bending moments via a continuum model, where we invoke the neo-Hookean strain energy model for the matrix material of fiber-reinforced composite, and assimilate the strain energy of fiber reinforcements into the matrix material model by accounting for stretching, bending, and twisting kinematics of the fibers through the computations of the first-order and second-order gradient of deformation. Emphasis is placed on deriving the Euler equation and boundary conditions of bending moment within the framework of the variational principle and configuring composite surfaces using differential geometry. Significant attention has been given to illustrating the concurrent three-dimensional deformation of fiber composite, meshwork deformation, and fiber kinematics. The simulation results reveal that for a square fiber composite subjected to the out-of-plane bending moment, the maximum in-plane deformation of matrix material occurs along the diagonal direction of the domain while the center of the domain experiences weak in-plane deformation. Notably, the matrix material performs isotropic/anisotropic properties depending on the domain size/shape. In addition, the simulated unit fiber deformations reasonably validate the overall deformation of the network, underscoring that the deformations of the embedded fiber units govern the overall mechanical performance of the fiber meshwork. More importantly, the continuum model qualitatively provides reasonable predictions on the damage patterns of construction materials by demonstrating the kinematics of matrix material and meshwork deformation.

Modelling ductile damage in metals and alloys through Weyl condition exploiting local gauge symmetries

Local translational and scaling symmetries in space–time is exploited for modelling ductile damage in metals and alloys over wide ranges of strain rate and temperature. The invariant energy density corresponding to the ductile deformation is constructed through the gauge invariant curvature tensor by imposing the Weyl like condition. In contrast, the energetics of the plastic deformation is brought in through the gauge compensating field emerged due to local translation and attempted to explore the geometric interpretation of certain internal variables often used in classical viscoplasticity models. Invariance of the energy density under the local action of translation and scaling is preserved through minimally replaced space–time gauge covariant operators. Minimal replacement introduces two non-trivial gauge compensating fields pertaining to local translation and scaling. These are used to describe ductile damage, including plastic flow and micro-crack evolution in the material. A space–time pseudo-Riemannian metric is used to lay out the kinematics in a finite-deformation setting. Recognizing the available insights in classical theories of viscoplasticity, we also establish a correspondence of the gauge compensating field due to spatial translation with Kröner’s multiplicative decomposition of the deformation gradient. Thermodynamically consistent coupling between viscoplasticity and ductile damage is ensured through an appropriate degradation function. Non-ordinary state-based (NOSB) peridynamics (PD) discretization of the model is used for numerical implementation. We conduct simulations of uniaxial deformation to validate the model against available experimental evidence and to assess its predictive features. The model’s viability is tested in reproducing a few experimentally known facts, viz., strain rate locking in the stress–strain response, whose origin is traced to a nonlinear microscopic inertia term arising out of the space–time translation symmetry. Finally, we solved 2D and axisymmetric deformation problems for qualitatively validating the model’s viability. NOSB peridynamics axisymmetric formulation in finite deformation setup is also presented.

Extension of the a posteriori finite element method (APFEM) to geometrical alterations and application to stochastic homogenisation

AbstractWe recently proposed an efficient method facilitating the parametric study of a finite element mechanical simulation as a postprocessing step, that is, without the need to run multiple simulations: the a posteriori finite element method (APFEM). APFEM only requires the knowledge of the vertices of the parameter space and is able to predict accurately how the degrees of freedom of a simulation, i.e., nodal displacements, and other outputs of interests, for example, element stress tensors, evolve when simulation parameters vary within their predefined ranges. In our previous work, these parameters were restricted to material properties and loading conditions. Here, we extend the APFEM to additionally account for changes in the original geometry. This is achieved by defining an intermediary reference frame whose mapping is defined stochastically in the weak form. Subsequent deformation is then reached by correcting for this stochastic variation in the reference frame through multiplicative decomposition of the deformation gradient tensor. The resulting framework is shown here to provide accurate mechanical predictions for relevant applications of increasing complexity: (i) quantifying the stress concentration factor of a plate under uniaxial loading with one and two elliptical holes of varying eccentricities, and (ii) performing the stochastic homogenisation of a composite plate with uncertain mechanical properties and geometry inclusion. This extension of APFEM completes our original approach to account parametrically for geometrical alterations, in addition to boundary conditions and material properties. The advantages of this approach in our original work in terms of stochastic prediction, uncertainty quantification, structural and material optimisation and Bayesian inferences are all naturally conserved.

Investigation of radial basis function dynamic mesh method with rotation correction based on adaptive background mesh

Strategies to improve the quality, robustness, and efficiency of dynamic mesh deformation are investigated. Radial basis function (RBF) interpolation is used for the background mesh, and volume-weighted interpolation is utilized for the computational mesh. An adaptive background mesh with an octree structure is automatically generated. It is refined according to the resolution near the boundaries and the wall's interior is also covered by mesh which is of sufficient resolution and prevents a negative volume due to large deformations of complex configurations. Secondly, the RBF method based on adaptive background mesh, the RBF interpolation angle method based on adaptive background mesh, and the damping function-RBF interpolation angle method based on adaptive background mesh have been developed. This method is applicable to complex shapes with multiple gaps, ensuring a high-quality mesh and high computational efficiency for large deformations. The numerical validation of the interpolation rotation angles for different RBFs indicates that the CTPS C0 function is suitable for large gradients in the boundary deformation. The RBF interpolation angle method ensures a smoother deformation for rigid and large deformations. The damping function-RBF interpolation angle method ensures mesh orthogonality near the boundary during elastic deformation, improving the boundary layer's quality and deformation robustness. Several examples show that the proposed method improves the mesh quality by about 30% and increases the time consumption by about 10%. In addition, the methods are applicable to rigid/elastic deformation of objects with multiple boundaries.

Visualization of Deformation and Stress Waves in Wooden Solid and Glued Elements of Building Structures.

The basis of the research hypothesis is the assumption that in wooden structures, deformations and stresses propagate in waves. The numerical experiment demonstrated a correct qualitative visual picture of the wave propagation of deformations, with wave manifestations and characteristic effects on the surface of the sample, at axial and corner points. Visually, the numerical model showed Rayleigh waves on the surface layer of the sample, depending on the ratio of the external geometric dimensions of the sample model, with pronounced wave interference on the outer shell. The visual manifestation of deformation on the outer sides (faces) and the reflection of deformation waves from the outer boundaries of the elastic medium of the sample in the form of Rayleigh waves confirm the correctness of the general hypothesis and the implemented model. Visualization of the process of emergence, propagation and attenuation of deformation waves on the surface of the sample shows that in the quantitative description of the deformation gradient, areas dangerous for the material can be identified.

Numerical simulation of three dimensional concrete printing based on a unified fluid and solid mechanics formulation

Deformation control constitutes one of the main technological challenges in three dimensional (3D) concrete printing, and it presents a challenge that must be addressed to achieve a precise and reliable construction process. Model-based information of the expected deformations and stresses is required to optimize the construction process in association with the specific properties of the concrete mix. In this work, a novel thermodynamically consistent finite strain constitutive model for fresh and early-age 3D-printable concrete is proposed. The model is then used to simulate the 3D concrete printing process to assess layer shapes, deformations, forces acting on substrate layers and prognoses of possible structural collapse during the layer-by-layer buildup. The constitutive formulation is based on a multiplicative split of the deformation gradient into elastic, aging and viscoplastic parts, in combination with a hyperelastic potential and considering evolving material properties to account for structural buildup or aging. One advantage of this model is the stress-update-scheme, which is similar to that of small strain plasticity and therefore enables an efficient integration with existing material routines. The constitutive model uses the particle finite element method, which serves as the simulation framework, allowing for modeling of the evolving free surfaces during the extrusion process. Computational analyses of three printed layers are used to create deformation plots, which can then be used to control the deformations during 3D concrete printing. This study offers further investigations, on the structural level, focusing on the potential structural collapse of a 3D printed concrete wall. The capability of the proposed model to simulate 3D concrete printing processes across the scales—from a few printed layers to the scale of the whole printed structure—in a unified fashion with one constitutive formulation, is demonstrated.