A three-dimensional hybrid smoothed finite element method (H-SFEM) for nonlinear solid mechanics problems

This paper presents a novel three-dimensional hybrid smoothed finite element method (H-SFEM) for solid mechanics problems. In 3D H-SFEM, the strain field is assumed to be the weighted average between compatible strains from the finite element method (FEM) and smoothed strains from the node-based smoothed FEM with a parameter α equipped into H-SFEM. By adjusting α, the upper and lower bound solutions in the strain energy norm and eigenfrequencies can always be obtained. The optimized α value in 3D H-SFEM using a tetrahedron mesh possesses a close-to-exact stiffness of the continuous system, and produces ultra-accurate solutions in terms of displacement, strain energy and eigenfrequencies in the linear and nonlinear problems. The novel domain-based selective scheme is proposed leading to a combined selective H-SFEM model that is immune from volumetric locking and hence works well for nearly incompressible materials. The proposed 3D H-SFEM is an innovative and unique numerical method with its distinct features, which has great potential in the successful application for solid mechanics problems.

A multiscale study on the structural and mechanical properties of the luffa sponge from Luffa cylindrica plant

Cellular materials that are often observed in biological systems exhibit excellent mechanical properties at remarkably low densities. Luffa sponge is one of such materials with a complex interconnecting porous structure. In this paper, we studied the relationship between its structural and mechanical properties at different levels of its hierarchical organization from a single fiber to a segment of whole sponge. The tensile mechanical behaviors of three single fibers were examined by an Instron testing machine and the ultrastructure of a fractured single fiber was observed in a scanning electronic microscope. Moreover, the compressive mechanical behaviors of the foam-like blocks from different locations of the sponge were examined. The difference of the compressive stress-strain responses of four sets of segmental samples were also compared. The result shows that the single fiber is a porous composite material mainly consisting of cellulose fibrils and lignin/hemicellulose matrix, and its Young's modulus and strength are comparable to wood. The mechanical behavior of the block samples from the hoop wall is superior to that from the core part. Furthermore, it shows that the influence of the inner surface on the mechanical property of the segmental sample is stronger than that of the core part; in particular, the former's Young's modulus, strength and strain energy absorbed are about 1.6 times higher. The present work can improve our understanding of the structure-function relationship of the natural material, which may inspire fabrication of new biomimetic foams with desirable mechanical efficiency for further applications in anti-crushing devices and super-light sandwich panels.

Nonlinear mechanical property of tracheal cartilage: A theoretical and experimental study

Background: Despite being the stiffest airway of the bronchial tree, the trachea undergoes significant deformation due to intrathoracic pressure during breathing. The mechanical properties of the trachea affect the flow in the airway and may contribute to the biological function of the lung. Method: A Fung-type strain energy density function was used to investigate the nonlinear mechanical behavior of tracheal cartilage. A bending test on pig tracheal cartilage was performed and a mathematical model for analyzing the deformation of tracheal cartilage was developed. The constants included in the strain energy density function were determined by fitting the experimental data. Result: The experimental data show that tracheal cartilage is a nonlinear material displaying higher strength in compression than in tension. When the compression forces varied from -0.02 to -0.03 N and from -0.03 to -0.04 N, the deformation ratios were 11.03±2.18% and 7.27±1.59%, respectively. Both were much smaller than the deformation ratios (20.01±4.49%) under tension forces of 0.02 to 0.01 N. The Fung-type strain energy density function can capture this nonlinear behavior very well, whilst the linear stress-strain relation cannot. It underestimates the stability of trachea by exaggerating the displacement in compression. This study may improve our understanding of the nonlinear behavior of tracheal cartilage and it may be useful for the future study on tracheal collapse behavior under physiological and pathological conditions.

Simulation of the interaction between blood flow and atherosclerotic plaque

It has been well accepted that over 50% of cerebral ischemic events are the result of rupture of vulnerable carotid atheroma and subsequent thrombosis. Such strokes are potentially preventable by carotid interventions. Selection of patients for intervention is currently based on the severity of carotid luminal stenosis. It has been, however, widely accepted that luminal stenosis alone may not be an adequate predictor of risk. To evaluate the effects of degree of luminal stenosis and plaque morphology on plaque stability, we used a coupled nonlinear time-dependent model with flow-plaque interaction simulation to perform flow and stress/strain analysis for stenotic artery with a plaque. The Navier-Stokes equations in the Arbitrary Lagrangian-Eulerian (ALE) formulation were used as the governing equations for the fluid. The Ogden strain energy function was used for both the fibrous cap and the lipid pool. The plaque Principal stresses and flow conditions were calculated for every case when varying the fibrous cap thickness from 0.1 to 2mm and the degree of luminal stenosis from 10% to 90%. Severe stenosis led to high flow velocities and high shear stresses, but a low or even negative pressure at the throat of the stenosis. Higher degree of stenosis and thinner fibrous cap led to larger plaque stresses, and a 50% decrease of fibrous cap thickness resulted in a 200% increase of maximum stress. This model suggests that fibrous cap thickness is critically related to plaque vulnerability and that, even within presence of moderate stenosis, may play an important role in the future risk stratification of those patients when identified in vivo using high resolution MR imaging.

DSC studies on the phase transformation of hydrazonium sulfate

The high temperature phase transformation of hydrazonium sulfate, N2H6SO4 has been studied using DSC. The enthalpy of phase transition is found to be 3.63 ± 0.1 kJ mole−1. The phase transition temperature is found to decrease with the increase of particle size. It appears that the strain energy and not surface energy, is responsible for the phase transformation. The molar volume of the salt increases during the transformation as found by the dilatometric experiment involving percentage of linear thermal expansion. On cooling, the transformation from the high temperature modification to orthorhombic form is incomplete and extends over a wide range of temperature.

Macromechanical behaviour of composite laminates

An improved higher order transverse shear deformation theory is employed to arrive at modified constitutive relations which can be used in the flexural, buckling and vibration analysis of laminated plates and shells. The strain energy for such systems is then expressed in terms of the displacements and the rotations for ready reference and use. Numerical values of vibration frequencies are obtained using this formulation employing Ritz's method of analysis. The results are compared with those available in the literature to validate the analysis presented.

Fracture Behavior of Multidirectional DCB Specimen: Higher-Order Beam Theories

Mathematical models, for the stress analysis of symmetric multidirectional double cantilever beam (DCB) specimen using classical beam theory, first and higher-order shear deformation beam theories, have been developed to determine the Mode I strain energy release rate (SERR) for symmetric multidirectional composites. The SERR has been calculated using the compliance approach. In the present study, both variationally and nonvariationally derived matching conditions have been applied at the crack tip of DCB specimen. For the unidirectional and cross-ply composite DCB specimens, beam models under both plane stress and plane strain conditions in the width direction are applicable with good performance where as for the multidirectional composite DCB specimen, only the beam model under plane strain condition in the width direction appears to be applicable with moderate performance. Among the shear deformation beam theories considered, the performance of higher-order shear deformation beam theory, having quadratic variation for transverse displacement over the thickness, is superior in determining the SERR for multidirectional DCB specimen.

Fracture analysis of stiffened panels under combined tensile, bending, and shear loads

The objective is to present the formulation of numerically integrated modified virtual crack closure integral technique for concentrically and eccentrically stiffened panels for computation of strain-energy release rate and stress intensity factor based on linear elastic fracture mechanics principles. Fracture analysis of cracked stiffened panels under combined tensile, bending, and shear loads has been conducted by employing the stiffened plate/shell finite element model, MQL9S2. This model can be used to analyze plates with arbitrarily located concentric/eccentric stiffeners, without increasing the total number of degrees of freedom, of the plate element. Parametric studies on fracture analysis of stiffened plates under combined tensile and moment loads have been conducted. Based on the results of parametric,studies, polynomial curve fitting has been carried out to get best-fit equations corresponding to each of the stiffener positions. These equations can be used for computation of stress intensity factor for cracked stiffened plates subjected to tensile and moment loads for a given plate size, stiffener configuration, and stiffener position without conducting finite element analysis.

Elastic behaviour of matter under very high pressures - Uniform compression

In order to study the elastic behaviour of matter when subjected to very large pressures, such as occur for example in the interior of the earth, and to provide an explanation for phenomena like earthquakes, it is essential to be able to calculate the values of the elastic constants of a substance under a state of large initial stress in terms of the elastic constants of a natural or stress-free state. An attempt has been made in this paper to derive expressions for these quantities for a substance of cubic symmetry on the basis of non-linear theory of elasticity and including up to cubic powers of the strain components in the strain energy function. A simple method of deriving them directly from the energy function itself has been indicated for any general case and the same has been applied to the case of hydrostatic compression. The notion of an effective elastic energy-the energy require to effect an infinitesimal deformation over a state of finite strain-has been introduced, the coefficients in this expression being the effective elastic constants. A separation of this effective energy function into normal co-ordinates has been given for the particular case of cubic symmetry and it has been pointed out, that when any of such coefficients in this normal form becomes negative, elastic instability will set in, with associated release of energy.

Effect of interaction of macrocracks on the stress intensity factor in a beam

Beams with a central edge crack, as well as other noncentral vertical and inclined edge cracks distributed symmetrically, subjected to three-point as well as four-point bending, are analysed using the finite element technique. Values of stress intensity factor K1 at the central crack tip for a crack-to-beam depth ratio Image equal to 0.5, are calculated for various cracked-beam configurations, using the compliance calibration technique as well as method of strain energy release rate. These are compared with the value of K1 for the case of a beam with central edge crack alone. Results of the present parametric study are used to specify the range of values pertaining to basic parameters such as crack-to-beam depth ratios, geometry and position with respect to central edge crack, of other macrocracks for which interaction exists. Accordingly, the macrocracks are classified as either interacting or noninteracting types. Hence for noninteracting types of cracks, analytical expressions available for the determination of K1 in the case of beam with a central edge crack alone, are applicable.

New a posteriori error estimator and adaptive mesh refinement strategy for 2-D crack problems

A posteriori error estimation and adaptive refinement technique for fracture analysis of 2-D/3-D crack problems is the state-of-the-art. The objective of the present paper is to propose a new a posteriori error estimator based on strain energy release rate (SERR) or stress intensity factor (SIF) at the crack tip region and to use this along with the stress based error estimator available in the literature for the region away from the crack tip. The proposed a posteriori error estimator is called the K-S error estimator. Further, an adaptive mesh refinement (h-) strategy which can be used with K-S error estimator has been proposed for fracture analysis of 2-D crack problems. The performance of the proposed a posteriori error estimator and the h-adaptive refinement strategy have been demonstrated by employing the 4-noded, 8-noded and 9-noded plane stress finite elements. The proposed error estimator together with the h-adaptive refinement strategy will facilitate automation of fracture analysis process to provide reliable solutions.

A Bilinear Constitutive Model for Isotropic Bimodulus Materials

Proper formulation of stress-strain relations, particularly in tension-compression situations for isotropic biomodulus materials, is an unresolved problem. Ambartsumyan's model [8] and Jones' weighted compliance matrix model [9] do not satisfy the principle of coordinate invariance. Shapiro's first stress invariant model [10] is too simple a model to describe the behavior of real materials. In fact, Rigbi [13] has raised a question about the compatibility of bimodularity with isotropy in a solid. Medri [2] has opined that linear principal strain-principal stress relations are fictitious, and warned that the bilinear approximation of uniaxial stress-strain behavior leads to ill-working bimodulus material model under combined loading. In the present work, a general bilinear constitutive model has been presented and described in biaxial principal stress plane with zonewise linear principal strain-principal stress relations. Elastic coefficients in the model are characterized based on the signs of (i) principal stresses, (ii) principal strains, and (iii) on the value of strain energy component ratio ER greater than or less than unity. The last criterion is used in tension-compression and compression-tension situations to account for different shear moduli in pure shear stress and pure shear strain states as well as unequal cross compliances.

Fracture Properties of Concrete-Concrete Interfaces Using Digital Image Correlation

The mode I and mode II fracture toughness and the critical strain energy release rate for different concrete-concrete jointed interfaces are experimentally determined using the Digital Image Correlation technique. Concrete beams having different compressive strength materials on either side of a centrally placed vertical interface are prepared and tested under three-point bending in a closed loop servo-controlled testing machine under crack mouth opening displacement control. Digital images are captured before loading (undeformed state) and at different instances of loading. These images are analyzed using correlation techniques to compute the surface displacements, strain components, crack opening and sliding displacements, load-point displacement, crack length and crack tip location. It is seen that the CMOD and vertical load-point displacement computed using DIC analysis matches well with those measured experimentally.

A general procedure for modified crack closure integral in 3D problems with cracks

In linear elastic fracture mechanics (LEFM), Irwin's crack closure integral (CCI) is one of the signficant concepts for the estimation of strain energy release rates (SERR) G, in individual as well as mixed-mode configurations. For effective utilization of this concept in conjunction with the finite element method (FEM), Rybicki and Kanninen [Engng Fracture Mech. 9, 931 938 (1977)] have proposed simple and direct estimations of the CCI in terms of nodal forces and displacements in the elements forming the crack tip from a single finite element analysis instead of the conventional two configuration analyses. These modified CCI (MCCI) expressions are basically element dependent. A systematic derivation of these expressions using element stress and displacement distributions is required. In the present work, a general procedure is given for the derivation of MCCI expressions in 3D problems with cracks. Further, a concept of sub-area integration is proposed which facilitates evaluation of SERR at a large number of points along the crack front without refining the finite element mesh. Numerical data are presented for two standard problems, a thick centre-cracked tension specimen and a semi-elliptical surface crack in a thick slab. Estimates for the stress intensity factor based on MCCI expressions corresponding to eight-noded brick elements are obtained and compared with available results in the literature.

Cracks emanating from pin-loaded lugs

Pin-loaded lugs were analysed in the presence of cracks emanating from circular holes. The analysis presents a unified treatment of interference, push or clearance fit pins. Both metallic (isotropic) and composite (orthotropic) plates were dealt with. The finite element model used special singular six-noded quadrilateral elements at the crack tip. The non-linear load contact behaviour at the pin-hole interface was dealt with by an inverse technique. A modified crack closure integral (MCCI) technique was used to evaluate the strain energy release rates (SERRs) and stress intensity factors (SIFs) at the crack tips. Numerical results are presented showing the non-linear variation of SIF with applied stress, and the influence of the amount of interference or clearance and the interfacial friction on SIF.