Search results

This paper presents a comparison of various models for strain‐softening due to damage such as cracking or void growth, as proposed recently in the literature. Continuum‐based models expressed in terms of softening stress—strain relations, and fracture‐type models expressed in terms of softening stress—displacement relations are distinguished. From one‐dimensional wave propagation calculations, it is shown that strain‐localization into regions of finite size cannot be achieved. The previously well‐documented spurious convergence is obtained with continuum models, while stress—displacement relations cannot model well smeared‐crack situations. Continuum models may, however, be used in general if a localization limiter is implemented. Gradient‐type localization limiters appear to be rather complicated; they require solving higher‐order differential equations of equilibrium with additional bourdary conditions. Non‐local localization limiters, especially the non‐local continuum with local strain, in which only the energy dissipating variables are non‐local, is found to be very effective, and also seems to be physically realistic. This formulation can correctly model the transition between homogeneous damage states and situations in which damage localizes into small regions that can be viewed as cracks. The size effect observed in the experimental and numerical response of specimens in tension or compression is shown to be a consequence of this progressive transition from continuum‐type to fracture‐type formulations.

The behaviour of cracked finite elements is investigated. It is shown that spurious kinematic modes may emerge when softening type constitutive laws are employed. These modes are not always suppressed by surrounding elements. This is exemplified for a double‐notched concrete beam and for a Crack‐Line‐Wedge‐Loaded Double‐Cantilever‐Beam (CLWL—DCB). The latter example has been analysed for a large variety of finite elements and integration schemes. To investigate the phenomenon in greater depth an eigenvalue analysis has been carried out for some commonly used finite elements.

Plane strain constitutive behaviour of von Mises and isotropic Hoffman materials is examined using single element tests. Two kinds of tests are conducted – (a) prescribed displacement tests; and (b) tests with a mixture of displacements and boundary tractions prescribed. While (a) are used to understand the manner of stress traversal on the yield surface in principal stress space, (b) are employed to study the load displacement response and the possibility of ensuing localization. Associated plasticity is assumed throughout. The tests are conducted using perfect and strain softening plasticity. It is found that for the von Mises criterion limited exact solutions can be evolved even under softening (or hardening) conditions. For isotropic Hoffman materials the nature of the stress traversal, load deflection response and the satisfaction of the localization conditions are strongly influenced by the ratio and difference of uniaxial yield strengths, in tension and compression, as well as by the softening parameters.

The localized strain softening behaviour of concrete has been modelled by two approaches: (i) the stiffness degrading model based on the total stress‐strain constitutive relationship, and (ii) the tangent softening model based on the incremental stress‐strain relationship. The models are implemented using a new softening initiation criterion proposed for application in multi‐dimensional finite element analysis. Parametric analyses on plain concrete beams, tested experimentally by other researchers, have been carried out to investigate the required numerical efforts, the mesh objectivity, and the energy dissipation characteristics of the structures. The stiffness degrading model is very stable even when applied with relatively coarse finite element meshes. However, the computational demand of this model is relatively high. The combination of a total stress‐strain constitutive relationship to compute the element responses, and an incremental relationship to formulate the stiffness matrix, appears to be computationally efficient and stable, provided that adequately refined finite element mesh is used to model the structure.

This paper aims to develop a finite element analysis strategy, which is suitable for the analysis of progressive failure that occurs in pressure-dependent materials in practical engineering problems.

The numerical difficulties stemming from the strain-softening behaviour of the frictional material, which is represented by a non-associated Drucker–Prager material model, is tackled using the Cosserat continuum theory, while the mixed finite element formulation based on Hu–Washizu variational principle is adopted to allow the utilization of low-order finite elements.

The effectiveness and robustness of the low-order finite element are verified, and the simulation for a real-world landslide which occurred at the upstream side of Carsington embankment in Derbyshire reconfirms the advantages of the developed elastoplastic Cosserat continuum scheme in capturing the entire progressive failure process when the strain-softening and the non-associated plastic law are involved.

The permit of using low-order finite elements is of great importance to enhance computational efficiency for analysing large-scale engineering problems. The case study reconfirms the advantages of the developed elastoplastic Cosserat continuum scheme in capturing the entire progressive failure process when the strain-softening and the non-associated plastic law are involved.

The dispersive behaviour of waves in softening problems is analysed. Attention is focused on the influence of the numerical scheme on the dispersion characteristics in the process of localization of deformation. Distinction has been made between softening models defined in a standard plasticity framework and in a gradient‐dependent plasticity theory. Waves in a standard softening plasticity continuum do not disperse but due to spatial discretization dispersion is introduced which results in a mesh size dependent length scale effect. On the other hand, wave propagation in a gradient‐dependent softening plasticity continuum is dispersive. By carrying out the dispersion analysis on the discretized system the influence of numerical dispersion on material dispersion can be quantified which enables us to determine the accuracy for the solution of the localization zone. For a modelling with and without the inclusion of strain gradients accuracy considerations with respect to mass discretization, finite element size, time integration scheme and time step have been carried out.

Classical continuum models, i.e. continuum models that do not incorporate an internal length scale, suffer from pathological mesh‐dependence when strain‐softening models are employed in failure analyses. In this contribution the governing field equations are regularized by adding rotational degrees‐of‐freedom to the conventional translational degrees‐of‐freedom. This so‐called elasto‐plastic Cosserat continuum model, for which an efficient and accurate integration algorithm and a consistent tangent operator are also derived in this contribution, warrants convergence of the load—deflection curve to a unique solution upon mesh refinement and a finite width of the localization zone. This is demonstrated for an infinitely long shear layer and a biaxial test of a strain‐softening elasto‐plastic von Mises material.

Classical continuum models, i.e. continuum models that do not incorporate an internal length scale, suffer from excessive mesh dependence when strain‐softening models are used in numerical analyses and cannot reproduce the size effect commonly observed in quasi‐brittle failure. In this contribution three different approaches will be scrutinized which may be used to remedy these two intimately related deficiencies of the classical theory, namely (i) the addition of higher‐order deformation gradients, (ii) the use of micropolar continuum models, and (iii) the addition of rate dependence. By means of a number of numerical simulations it will be investigated under which conditions these enriched continuum theories permit localization of deformation without losing ellipticity for static problems and hyperbolicity for dynamic problems. For the latter class of problems the crucial role of dispersion in wave propagation in strain‐softening media will also be highlighted.

This study aims to propose a novel viscoelastic–viscoplastic combined constitutive model for glassy amorphous polymers within the framework of thermodynamics at finite strain that is capable of capturing their rate-dependent inelastic mechanical behavior in wide ranges of deformation rate and amount.

The rheology model whose viscoelastic and viscoplastic elements are connected in series is set in accordance with the multi-mechanism theory. Then, the constitutive functions are formulated on the basis of the multiplicative decomposition of the deformation gradient implicated by the rheology model within the framework of thermodynamics. Dynamic mechanical analysis (DMA) and loading/unloading/no-load tests for polycarbonate (PC) are conducted to identify the material parameters and demonstrate the capability of the proposed model.

The performance was validated in comparison with the series of the test results with different rates and amounts of deformation before unloading together. It has been confirmed that the proposed model can accommodate various material behaviors empirically observed, such as rate-dependent elasticity, elastic hysteresis, strain softening, orientation hardening and strain recovery.

This paper presents a novel rheological constitutive model in which the viscoelastic element connected in series with the viscoplastic one exclusively represents the elastic behavior, and each material response is formulated according to the multiplicatively decomposed deformation gradients. In particular, the yield strength followed by the isotropic hardening reflects the relaxation characteristics in the viscoelastic constitutive functions so that the glass transition temperature could be variant within the wide range of deformation rate. Consequently, the model enables us to properly represent the loading process up to large deformation regime followed by unloading and no-load processes.

This study aims to explore the compressive behavior of hollow triply periodic minimal surface (HTPMS) cellular structures by selective laser melting (SLM).

This study presents a design method for gyroid hollow triply periodic minimal surfaces (G-HTPMS) and primitive hollow triply periodic minimal surfaces (P-HTPMS) cellular structures, and SLM technology was applied to manufacture these cellular structures. Compressive behaviors and energy absorption behaviors of hollow cellular structures were researched in this study.

Compared with normal gyroid triply periodic minimal surfaces (G-TPMS) and normal primitive triply periodic minimal surfaces (P-TPMS), the G-HTPMS and P-HTPMS have higher elastic modulus, plateau stress and effective energy absorption under uniaxial compression. The hollow design in HTPMS can enhance the mechanical properties and energy absorption of the cellular structure. Finite element analysis also demonstrates that the hollow design can reduce stress concentration, which improved the compressive curves from a severely fluctuating state to a relatively flat state and reduces fracture. According to compressive behaviors, G-TPMS and G-HTPMS are the bending-dominated cellular structures with strain hardening characteristics, and P-TPMS and P-HTPMS are the stretching-dominated cellular structures with strain softening characteristics.

This research provided a design method for HTPMS, and it was proved that the mechanical properties increased by hollow design inspired by bamboo.