Search results

The purpose of this paper is to describe a general theoretical and finite element implementation framework for the constitutive modelling of biological soft tissues.

The model is based on continuum fibers reinforced composites in finite strains. As an extension of the isotropic hyperelasticity, it is assumed that the strain energy function is decomposed into a fully isotropic component and an anisotropic component. Closed form expressions of the stress tensor and elasticity tensor are first established in the general case of fully incompressible plane stress which orthotropic and transversely isotropic hyperelasticity. The incompressibility is satisfied exactly.

Numerical examples are presented to illustrate the model's performance.

The paper presents a constitutive model for incompressible plane stress transversely isotropic and orthotropic hyperelastic materials.

A general first-invariant constitutive model has been derived in literature for incompressible, isotropic hyperelastic materials, known as Marlow model, which reproduces test data exactly without the need of curve-fitting procedures. This paper aims to describe how to extend Marlow’s constitutive model to the more general case of compressible hyperelastic materials.

The isotropic constitutive model is based on a strain energy function, whose isochoric part is solely dependent on the first modified strain invariant. Based on Marlow’s idea, a principle of energetically equivalent deformation states is derived for the compressible case, which is used to determine the underlying strain energy function directly from measured test data. No particular functional of the strain energy function is assumed. It is shown how to calibrate the volumetric and isochoric strain energy functions uniquely with uniaxial or biaxial test data only. The constitutive model is implemented into a finite element program to demonstrate its applicability.

The model is well suited for use in finite element analysis. Only one set of test data is required for calibration without any need for curve-fitting procedures. These test data are reproduced exactly, and the model prediction is reasonable for other deformation modes.

Marlow’s basic concept is extended to the compressible case and applied to both the volumetric and isochoric part of the compressible strain energy function. Moreover, a novel approach is described on how both compressive and tensile test data can be used simultaneously to calibrate the model.

This paper aims to describe a shape optimization for hyperelastic axisymmetric structure with an exact sensitivity method.

The whole shape optimization process is carried out by integrating a closed geometric shape in the real space R2 with boundaries defined by B-splines curves. An exact sensitivity analysis and a mathematical programming method (SQP: Sequential Quadratic Programming) are implemented. The design variables are the control points' coordinates which minimize the Von-Mises criteria, with a constraint that the total material volume of the structure remains constant. The feasibility of the proposed methods is carried out by two numerical examples. Results show that the exact Jacobian has an important computing time reduction.

Numerical examples are presented to illustrate its performance.

In this work, the sensitivity performance is computed using two numerical methods: the efficient finite difference scheme and the exact Jacobian.

The purpose of this paper is to analyze, computationally, the kinematic response (including large‐scale rotation and deformation, buckling, plastic yielding, failure initiation, fracture and fragmentation) of a pick‐up truck to the detonation of a landmine (shallow‐buried in one of six different soils, i.e. either sand, clay‐laden sand or sandy gravel, each in either dry or water‐saturated conditions, and detonated underneath the vehicle) using ANSYS/Autodyn, a general‐purpose transient non‐linear dynamics analysis software.

The computational analysis, using ANSYS/Autodyn, a general‐purpose transient non‐linear dynamics analysis software, included the interactions of the gaseous detonation products and the sand ejecta with the vehicle and the transient non‐linear dynamics response of the vehicle.

The results obtained clearly show the differences in the blast loads resulting from the landmine detonation in dry and saturated sand, as well as the associated kinematic response of the vehicle. It was also found that the low frequency content of the blast loads which can match the whole‐vehicle eigen modes is quite small so that resonance plays a minor role in the kinematic/ballistic response of the vehicle. Furthermore, it was demonstrated that mine blast analytical loading functions which are often used in transient non‐linear dynamic analyses have limited value when used in the analyses of a complete vehicle.

This is the first time that the kinematic response of a pick‐up truck to the detonation of a shallow‐buried landmine (using a full‐scale/complete model) has been analyzed computationally.

Several a posteriori error indicators and error estimators for frictionless contact problems are compared. In detail, residual based error estimators, error indicators relying on superconvergence properties and error estimators based on duality principles are investigated. Applications are to 2D solids under the hypothesis of nonlinear elastic material behaviour associated with finite deformations. A penalization technique is applied to enforce multilateral boundary conditions due to contact. The approximate solution of the problem is obtained by using the finite element method. Several numerical results are reported to show the applicability of the adaptive algorithm to the considered problems.

A direct and unified approach is proposed toward simultaneously simulating large strain elastic behaviors of gellan gels with different gellan polymer concentrations. The purpose of this paper is to construct an elastic potential with certain parameters of direct physical meanings, based on well-designed invariants of Hencky’s logarithmic strain.

For each given value of the concentration, the values of the parameters incorporated may be determined in the sense of achieving accurate agreement with large strain uniaxial extension and compression data. By means of a new interpolating technique, each parameter as a function of the concentration is then obtained from a given set of parameter values for certain concentration values.

Then, the effects of gellan polymer concentrations on large strain elastic behaviors of gellan gels are studied in demonstrating how each parameter relies on the concentration. Plane-strain (simple shear) responses are also presented for gellan gels with different polymer concentrations.

A direct, unified approach was proposed toward achieving a simultaneous simulation of large elastic strain behaviors of gellan gels for different gellan polymer concentrations. Each parameter incorporated in the proposed elastic potential will be derived as a function of the polymer concentration in an explicit form, in the very sense of simultaneously simulating large strain data for different concentrations.

Mechanical properties testing of the hyperelastic thermoplastic polyurethane (TPU) produced by the continuous digital light processing (CDLP) method of additive manufacturing. Primarily, this paper aims to verify that 3D printed TPU still satisfies commonly assumed volumetric incompressibility and material isotropy in elastic range. The secondary aim is to investigate the accuracy and reproducibility of the CDLP method.

Cylindrical samples were printed and subjected to a volumetric compression test to reveal their bulk modulus K and maximal theoretical porosity (MTP). Dog bone specimens were oriented along different axes and printed. Their dimensions were measured, and they were subjected to cyclic uniaxial tests up to 100% strain to reveal the level of stress softening and possible anisotropy. The hyperelastic Yeoh model was fitted to the mean response.

The authors measured the bulk modulus of K = 1851 ± 184 MPa. The mean MTP was 0.9 ± 0.5%. The mean response was identical in both directions and the data could be fitted by the isotropic third order Yeoh function with R^2 = 0.996. The dimensions measurement revealed the largest error (above 5%) in the direction perpendicular to the direction of the digital light projection while the dimensions in other two dimensions were much more accurate (0.75 and 1%, respectively).

The TPU printed by CDLP can be considered and modelled as isotropic and practically volumetrically incompressible. The parts in the printing chamber should be positioned in a way that the important dimensions are not parallel to the direction of the digital light projection.

The authors experimentally confirmed the volumetric incompressibility and mechanical isotropy of the TPU printed using the CDLP method.

This paper aims to present the findings of a numerical investigation into the performance of the steel-concrete composite floor involved in Broadgate Phase 8 fire.

The investigation is conducted by carrying out a 3-D thermomechanical analysis of a composite floor similar to the one involved in the fire using ANSYS. Four fire scenarios are investigated, with each producing a unique stress – strain pattern. The results obtained are compared with the observations made after the fire and inferences drawn.

The results obtained are found to be correlated with the observations made after the fire. The performance of the composite floor is found to be dominated by development of large strains, leading to large deflections. Furthermore, colder parts of the structure, through redistribution of forces, are found to have a profound impact on the ability of a composite floor to resist heating effects. From the findings, it is concluded that connections’ design, occurrence of membrane action and thermal restraints were the key reasons the floor did not fail.

The study takes a more forensic approach. This is a departure from majority of published literature, where comparison is usually between experimental and numerical results. By comparing the findings from a real fire with those of a numerical investigation, the study provides an insight into the accuracy of applying numerical models for the prediction of effects of fire on structural behaviour.

This paper aims to develop a non-linear finite element model predicting the response of externally strengthened beams under a three-point flexure test.

The ANSYS software is used for modeling. SOILD65, LINK180, SHELL181 and SOLID185 elements are used, respectively, to model concrete, steel reinforcement, polymer and steel plate support. A parametric study was carried out. The effects of compressive strength, Young’s modulus, layers number and carbon fiber-reinforced polymer thickness on beam behavior are analyzed. A comparative study between the non-linear finite element and analytical models, including the ACI 440.2 R-08 model, and experimental data is also carried out.

A comparative study of the non-linear finite element results with analytical models, including the ACI 440.2 R-08 model and experimental data for different parameters, shows that the strengthened beams possessed better resistance to cracks. In general, the finite element model’s results are in good agreement with the experimental test data.

This model will predict the strengthened beams behavior and can describe the beams physical conditions, yielding the results that can be interpreted in the structural study context without using a laboratory testing.

On the basis of the results, a good match is found between the model results and experimental data at all stages of loading the tested samples. Crack models obtained in the non-linear finite element model in the beams are also presented. The submitted finite element model can be used to predict the behavior of the reinforced concrete beam. Also, the comparative study between an analytical model proposed by of current code of ACI 440.2 R-08 and finite element analysis is investigated.

The application of numerical techniques to the solution of practical problems which exist in rubber technology is described. Structures and components in the form of reinforced rubber shells are widely used in industry and prediction of their performance is complicated by both the anisotropic nature of composite construction and the incompressible behaviour of the basic material. A layered shell element is developed for the solution of such problems with general anisotropic behaviour independently permitted in each layer. The approach adopted permits the easy location of reinforcement patterns. Numerical solution is based on a single field formulation by eliminating at integrating point level the Lagrange multiplier imposing the incompressible constraint. Large deformation, including large rotation, behaviour is accommodated and a total Lagrangian solution process is adopted. The code developed permits the simulation of non‐conservative loading and its versatility is demonstrated by the solution of some practical examples.