Search results

In this paper, a new rate type endochronic constitutive model is introduced to describe deformations in the finite strain range. A new material dependent objective rate of Cauchy stress is suggested based on the general form of spin tensors, defining objective stress rates. The endochronic constitutive equations are extended using the concept of corotational stress rates and additive decomposition of deformation rate. The constitutive relations are specialized for thin‐walled tubes under torsion and a procedure for solving the ordinary differential equations for cases of simple and pure torsion is developed. The axial effects for various materials, subjected to simple and pure torsion, are simulated and compared with experimental data. The results clearly indicate that the new combined rate endochronic model can be effectively used to describe the behavior of material in the finite strain range.

In recent years, the partially averaged Navier–Stokes (PANS) methodology has earned acceptability as a viable scale-resolving bridging method of turbulence. To further enhance its capabilities, especially for simulating separated flows past bluff bodies, this paper aims to combine PANS with a non-linear eddy viscosity model (NLEVM).

The authors first extract a PANS closure model using the Shih’s quadratic eddy viscosity closure model [originally proposed for Reynolds-averaged Navier–Stokes (RANS) paradigm (Shih et al., 1993)]. Subsequently, they perform an extensive evaluation of the combination (PANS + NLEVM).

The NLEVM + PANS combination shows promising result in terms of reduction of the anisotropy tensor when the filter parameter (f_k) is reduced. Further, the influence of PANS filter parameter f on the magnitude and orientation of the non-linear part of the stress tensor is closely scrutinized. Evaluation of the NLEVM + PANS combination is subsequently performed for flow past a square cylinder at Reynolds number of 22,000. The results show that for the same level of reduction in f_k, the PANS + NLEVM methodology releases significantly more scales of motion and unsteadiness as compared to the traditional linear eddy viscosity model (LEVM) of Boussinesq (PANS + LEVM). The authors further demonstrate that with this enhanced ability the NLEVM + PANS combination shows much-improved predictions of almost all the mean quantities compared to those observed in simulations using LEVM + PANS.

Based on these results, the authors propose the NLEVM + PANS combination as a more potent methodology for reliable prediction of highly separated flow fields.

Combination of a quadratic eddy viscosity closure model with PANS framework for simulating flow past bluff bodies.

The purpose of this paper is to examine quadratic convergence of finite element analysis for hyperelastic material at finite strains via Abaqus‐UMAT as well as classification of the rates of convergence for iterative solutions in regular cases.

Different formulations for stiffness – Hessian form of the free energy functionals – are systematically given for getting the rate‐independent analytical tangent and the numerical tangent as well as rate‐dependent tangents using the objective Jaumann rate of Kirchoff stress tensor as used in Abaqus. The convergence rates for available element types in Abaqus are computed and compared for simple but significant nonlinear elastic problems, such as using the 8‐node linear brick (B‐bar) element – also with hybrid pressure formulation and with incompatible modes – further the 20‐node quadratic brick element with corresponding modifications as well as the 6‐node linear triangular prism element and 4‐node linear tetrahedral element with modifications.

By using the Jaumann rate of Kirchoff stress tensor for both, rate dependent and rate independent problems, quadratic or nearly quadratic convergence is achieved for most of the used elements using Abaqus‐UMAT interface. But in case of using rate independent analytical tangent for rate independent problems, even convergence at all is not assured for all elements and the considered problems.

First time the convergence properties of 3D finite elements available in Abaqus sre systematically treated for elastic material at finite strain via Abaqus‐UMAT.

The model of elastoviscoplastic medium is based on the concept of slip proposed by Batdorf and Budiansky. The conditions of slip on a slip plane take into account the dependence of tangential stresses on slip velocity. But when the viscosity is low, they are almost similar to the dry friction conditions. Under some assumptions we succeeded in integrating the plastic shear rates over all possible slip planes in case of arbitrary three‐dimensional stress state and obtained an expression for the plastic strain rate tensor.

Derives a formulation for spatial stress tensors and spatial material tensors of hyperelastic material. Looks at a class of materials with the strain energy function decomposed into a volumetric and a deviatoric part. Separate terms formulate the strain energy with respect to the invariants of the left Cauchy‐Green tensor. Stress and material tensors, which play a crucial role in the solution process of the finite element formulation, are derived solely in the current configuration. Applies the described framework to several different constitutive models based on phenomenologically and physically motivated material descriptions. Proposes a formulation for the finite element implementation of van der Waals material. Compares numerical results with experimental investigations given in the literature. For three‐dimensional finite element computations standard elements and mixed elements, based on a three‐field variational principle where displacements, the hydrostatic pressure and the dilatations are independent variables, are used.

This paper shows a generalization of the classic isotropic plasticity theory to be applied to orthotropic or anisotropic materials. This approach assumes the existence of a real anisotropic space, and other fictitious isotropic space where a mapped fictitious problem is solved. Both spaces are related by means of a linear transformation using a fourth order transformation tensor that contains all the information concerning the real anisotropic material. The paper describes the basis of the spaces transformation proposed and the expressions of the resulting secant and tangent constitutive equations. Also details of the numerical integration of the constitutive equation are provided. Examples of application showing the good performance of the model for analysis of orthotropic materials and fibre‐reinforced composites are given.

The present paper is directed towards elasto‐plastic large deformation analysis of thin shells based on the concept of degenerated solids. The main aspect of the paper is the derivation of an efficient computational strategy placing emphasis on consistent elasto‐plastic tangent moduli and stress integration with the radial return method under the restriction of ‘zero normal stress condition’ in thickness direction. The advantageous performance of the standard Newton iteration using a consistent tangent stiffness matrix is compared to the classical scheme with an iteration matrix based on the infinitesimal elasto‐plastic constitutive tensor. Several numerical examples also demonstrate the effectiveness of the standard Newton iteration with respect to modified and quasi‐Newton methods like BFGS and others.

This paper seeks to explore the fiber orientation distribution and rheological properties of turbulent fiber suspensions flowing through a contraction.

The Reynolds averaged Navier‐Stokes equation was solved with the Reynolds stress model to get the mean fluid velocity and the turbulent kinetic energy in the turbulent flow of a contraction with rectangular cross‐section. The turbulent velocity fluctuations were represented as a Fourier series with random coefficients. Then the slender‐body theory was used to predict the fiber orientation distribution, orientation tensor, additional shear stress and first normal stress difference of suspensions in the flow.

It is found that the longer fibers tend to align the streamline easily. Increased contraction ratio results in higher fiber alignment in the direction of flow. The fibers are weakly and strongly aligned in the direction of flow in the region near the inlet and the exit, respectively. Fibers are significantly more aligned in the plane of the contraction than in the x‐z plane. Contraction ratio and fiber length were shown to strongly and weakly affect the distributions of additional shear stress and first normal stress difference.

It is the first time that the fiber orientation distribution and rheological properties of turbulent fiber suspensions flowing through a contraction have been computed numerically. The computational approach and results are valuable to the design and operation of contraction used in the industrial processes.

Regeneratively cooled thrust chamber is a key component of reusable liquid rocket engines. Subjected to cyclic thermal-mechanical loadings, its failure can seriously affect the service life of engines. QCr0.8 copper alloy is widely used in thrust chamber walls due to its excellent thermal conductivity, and its mechanical and fatigue properties are essential for the evaluation of thrust chamber life. This paper contributes to the understanding of the damage mechanism and material selection of regeneratively cooled thrust chambers for reusable liquid rocket engines.

In this paper, tensile and low-cycle fatigue (LCF) tests were conducted for QCr0.8 alloy, and a Chaboche combined hardening model was established to describe the elastic-plastic behavior of QCr0.8 at different temperatures and strain levels. In addition, an LCF life prediction model was established based on the Manson–Coffin formula. The reliability and accuracy of models were then verified by simulations in ABAQUS. Finally, the service life was evaluated for a regenerative cooling thrust chamber, under the condition of cyclic startup and shutdown.

In this paper, a Chaboche combined hardening model was established to describe the elastoplastic behavior of QCr0.8 alloy at different temperatures and strain levels through LCF experiments. The parameters of the fitted Chaboche model were simulated in ABAQUS, and the simulation results were compared with the experimental results. The results show that the model has high reliability and accuracy in characterizing the viscoplastic behavior of QCr0.8 alloy.

(1)The parameters of a Chaboche combined hardening constitutive model and LCF life equation were optimized by tensile and strain-controlled fatigue tests of QCr0.8 copper alloy. (2) Based on the Manson–Coffin formula, the reliability and accuracy of constitutive model were then verified by simulations in ABAQUS. (3)Thermal-mechanical analysis was carried out for regeneratively cooled thrust chamber wall of a reusable liquid rocket engine, and the service life considering LCF, creep and ratcheting damage was analyzed.

Microstructure void volume fraction is taken into account in finite element models developed for large strain elastoplastic problems. Void nucleation rate is related to matrix effective strain rate, void growth to material strain rate and associated elastoplastic potential available for porous material, void coalescence to matrix effective strain rate. The related radial return algorithm is described. Three types of computations are proposed: first, axisymmetric Q4 element traction are given as validation example; second, collar cylinder compression are computed as reference example; third, bulk forming are analysed as large strain specific example. Void volume fraction and hydrostatic stress are mainly discussed according to microvoids nucleation, growth and coalescence. Finally, the main interests of those computations are enhanced.