Limit analysis of orthotropic plates

The limit analysis problem for plates in bending is considered. The failure criterion for the material is assumed as orthotropic, with possible non-symmetric strength properties. According to Kirchhoff’s hypothesis, the plate is conceived as a superposition of layers, individually in plane stress situation, and continuity is enforced by means of a kinematic assumption. By exploiting previous results, recently established by the authors, the expression of the dissipation power per unit plate area is defined on this basis and the kinematic (upper bound) theorem of limit analysis is cast in a form suitable for numerical solutions. To this purpose, efficient algorithms successfully employed in the isotropic case can be used with minor modifications. The effectiveness of the procedure is demonstrated by solving some homogeneous plate examples. Results permit the assessment of the influence of different aspects, such as the ratio between strengths along the orthotropy directions, the tensile to compressive strength differential and the inclination of the orthotropy axes with respect to the sides. The effects of in-plane edge constraints are also discussed and it appears that they are emphasized considerably by anisotropy. Even if referred to specific cases, some conclusions can be regarded as fairly general.     

Development of a Hybrid Method of Reflection Photoelasticity for Crack Problems in Anisotropic Plates

Reflection type photoelastic experiment can be used more effectively than a transmission type photoelastic experiment especially in industrial fields. Moreover, composite materials have been widely used in engineering applications and structures because of their outstanding advantages which individual isotropic components do not have. The development of these materials requires a promising technique such as reflection photoelasticity to analyze their behaviors in service. Unfortunately, there are few experimental studies based on this technique. Therefore, a hybrid method based on this technique was developed in this research to analyze the fracture behavior of opaque anisotropic materials. The application of this method will help to understand the fracture behaviors of anisotropic materials used in engineering components and structures. The validity of this method was verified by comparison of the results obtained from this method with ones obtained from the hybrid methods for isotropic material on the same isotropic specimen. The reflection type photoelastic experiment for orthotropic materials was then applied to orthotropic plates with a central crack of various inclination angles. Using this hybrid method for anisotropic materials, stress intensity factors and separated stress components were obtained at the vicinity of the crack-tip in orthotropic plates from only the isochromatic fringe patterns of the isotropic coating material.     

Determining material true stress–strain curve from tensile specimens with rectangular cross-section

The uniaxial true stress logarithmic strain curve for a thick section can be determined from the load–diameter reduction record of a round tensile specimen. The correction of the true stress for necking can be performed by using the well-known Bridgman equation. For thin sections, it is more practical to use specimens with rectangular cross-section. However, there is no established method to determine the complete true stress–logarithmic strain relation from a rectangular specimen. In this paper, an extensive three-dimensional numerical study has been carried out on the diffuse necking behaviour of tensile specimens made of isotropic materials with rectangular cross-section, and an approximate relation is established between the area reduction of the minimum cross-section and the measured thickness reduction. It is found that the area reduction can be normalized by the uniaxial strain at maximum load which represents the material hardening and also the section aspect ratio. Furthermore, for the same material, specimens with different aspect ratio give exactly the same true average stress–logarithmic strain curve. This finding implies that Bridgmans correction can still be used for necking correction of the true average stress obtained from rectangular specimens. Based on this finding, a method for determining the true stress–logarithmic strain relation from the load–thickness reduction curve of specimens with rectangular cross-section is proposed.     

正交各向异性韧性材料应力-应变关系

采用大变形弹塑性有限元方法分析了各向同性和正交各向异性韧性材料光滑圆棒拉伸试件的颈缩问题．首先给出了采用计算机模拟确定各向同性韧性材料真实应力-应变曲线的具体方法；对正交各向异性韧性材料的分析表明，颈缩截面呈椭圆形，其长短轴方向的等效塑性应变基本上均匀分布，与Bridgman假设一致；轴向拉伸载荷-位移曲线与其它两方向的各向异性参数关系不大．在此基础上，建议了一种确定正交各向异性韧性材料真实应力-应变曲线的方法．     

Singular behavior at the tip of a crack in an elastic-plastic material with plastic orthotropy

The two-dimensional stress field at the tip of a crack in a plastically orthotropic material is analyzed by the total deformation theory of plasticity in conjunction with the J-integral. A model of a plastically orthotropic material is constructed by the use of the theory proposed by R. Hill (1950) and the uniaxial stress-strain relation suggested by W. Ramberg and W.R. Osgood (1943). It is found that the stresses in the vicinity of the crack tip have a singularity of the same order as that in the case of isotropic materials, but their amplitudes are greatly influenced by the plastic orthotropy. Numerical work is carried out for two typical metals, and the effect of the plastic orthotropy is examined for the stress field, the crack opening displacement, the strain energy density, and the shape of the elastic-plastic boundary.     

An anisotropic elastoplastic constitutive formulation generalised for orthotropic materials

This paper presents a finite strain constitutive model to predict a complex elastoplastic deformation behaviour that involves very high pressures and shockwaves in orthotropic materials using an anisotropic Hill’s yield criterion by means of the evolving structural tensors. The yield surface of this hyperelastic–plastic constitutive model is aligned uniquely within the principal stress space due to the combination of Mandel stress tensor and a new generalised orthotropic pressure. The formulation is developed in the isoclinic configuration and allows for a unique treatment for elastic and plastic orthotropy. An isotropic hardening is adopted to define the evolution of plastic orthotropy. The important feature of the proposed hyperelastic–plastic constitutive model is the introduction of anisotropic effect in the Mie–Gruneisen equation of state (EOS). The formulation is further combined with Grady spall failure model to predict spall failure in the materials. The proposed constitutive model is implemented as a new material model in the Lawrence Livermore National Laboratory (LLNL)-DYNA3D code of UTHM’s version, named Material Type 92 (Mat92). The combination of the proposed stress tensor decomposition and the Mie–Gruneisen EOS requires some modifications in the code to reflect the formulation of the generalised orthotropic pressure. The validation approach is also presented in this paper for guidance purpose. The \({\varvec{\psi }}\) tensor used to define the alignment of the adopted yield surface is first validated. This is continued with an internal validation related to elastic isotropic, elastic orthotropic and elastic–plastic orthotropic of the proposed formulation before a comparison against range of plate impact test data at 234, 450 and \({\mathrm {895\,ms}}^{\mathrm {-1}}\) impact velocities is performed. A good agreement is obtained in each test.     

Effect of transverse shear and material orthotropy in a cracked spherical cap

The elastostatic problem for a relatively thin-walled spherical cap containing a through crack is considered. The problem is formulated for a specially orthotropic material within the confines of a linearized, shallow shell theory. The theory used is equivalent to Reissner's theory of flat plates and hence permits the consideration of all five physical conditions on the shell boundaries separately. The solution of the problem is reduced to that of a pair of singular integral equations and the asymptotic stress state around the crack tips is investigated. The numerical solution of the problem is given for an isotropic shell and for two specially orthotropic shells. The results indicate that the material orthotropy as well as the shell curvature and thickness may have a considerable effect on the stress intensity factors at the crack tips.     

Membrane Curvatures and Stress-strain Full Fields of Axisymmetric Bulge Tests from 3D-DIC Measurements. Theory and Validation on Virtual and Experimental results

The bulge test is mostly used to analyze equibiaxial tensile stress state at the pole of inflated isotropic membranes. Three-dimensional digital image correlation (3D-DIC) technique allows the determination of three-dimensional surface displacements and strain fields. In this paper, a method is proposed to determine also the membrane stress tensor fields for in-plane isotropic materials, independently of any constitutive equation. Stress-strain state is then known at any surface point which enriches greatly experimental data deduced from the axisymmetric bulge tests. Our method consists, first in calculating from the 3D-DIC experimental data the membrane curvature tensor at each surface point of the bulge specimen. Then, curvature tensor fields are used to investigate axisymmetry of the test. Finally in the axisymmetric case, membrane stress tensor fields are determined from meridional and circumferential curvatures combined with the measurement of the inflating pressure. Our method is first validated for virtual 3D-DIC data, obtained by numerical simulation of a bulge test using a hyperelastic material model. Afterward, the method is applied to an experimental bulge test performed using as material a silicone elastomer. The stress-strain fields which are obtained using the proposed method are compared with results of the finite element simulation of this overall bulge test using a neo-Hookean model fitted on uniaxial and equibiaxial tensile tests.     

Orthotropic semi-infinite medium with a crack under thermal shock

A cracked orthotropic semi-infinite plate under thermal shock is investigated. The thermal stresses are generated due to sudden cooling of the boundary by ramp function temperature change. The superposition technique is used to solve the problem. The crack problem is formulated by applying the thermal stresses obtained from the uncracked plate with opposite sign to be the only external loads on the crack surfaces as the crack surface tractions. The Fourier transform technique is used to solve the problem leading to a singular equation of the Cauchy type. The singular integral equation is solved numerically using the expansion method. The influence of the material orthotropy on the stress intensity factors is shown by comparing the results obtained for different orthotropic materials and isotropic materials in the case of plane stress. The numerical results of the stress intensity factors are demonstrated as a function of time, crack length, location of the crack and the duration of the cooling rate.     

On the asymptotic crack-tip stress fields in nonlocal orthotropic elasticity

The crack-tip stress fields in orthotropic bodies are derived within the framework of Eringen’s nonlocal elasticity via the Green’s function method. The modified Bessel function of second kind and order zero is considered as the nonlocal kernel. We demonstrate that if the localisation residuals are neglected, as originally proposed by Eringen, the asymptotic stress tensor and its normal derivative are continuous across the crack. We prove that the stresses attained at the crack tip are finite in nonlocal orthotropic continua for all the three fracture modes (I, II and III). The relative magnitudes of the stress components depend on the material orthotropy. Moreover, non-zero self-balanced tractions exist on the crack edges for both isotropic and orthotropic continua. The special case of a mode I Griffith crack in a nonlocal and orthotropic material is studied, with the inclusion of the T-stress term.     

Material symmetry and the evolution of anisotropies in a simple material—I. Change of reference configuration

Noll's rule is used to determine the structure of a material symmetry group written with respect to one reference configuration when the representation of the symmetry with respect to another configuration is the traditional material symmetry group associated with isotropy, transverse isotropy or orthotropy, and for an arbitrary deformation gradient relating the two configurations. It is shown that the former symmetry group can contain an orthogonal subgroup. It is determined whether this subgroup is that for isotropic, transversely isotropic, orthotropic, monoclinic, or triclinic response, and the preferred directions of the symmetry are determined.     

Dependence of stress intensity factors on elastic constants for cracks in an orthotropic bimaterial with a thin film

An interface crack and a subinterface crack in an orthotropic bimaterial structure consisting of a thin film and a half plane substrate are analyzed. The orthotropic bimaterial structure is subjected to compressive load and bending moment per unit thickness. Complete expressions of stress intensity factors for the two cracks are obtained based on the path independence of the J integral, apart from one dimensionless parameter undetermined each. The dependence of the dimensionless parameters on material constants is examined. A reduction of the number of necessary material parameters for the parameters is made based upon the modified Stroh formalism. The explicit dependence of the dimensionless parameters on one orthotropic parameter for the film is determined by using the orthotropy rescaling technique. Variations of the dimensionless parameters with the other material parameters are also obtained through numerical computations.     

Optimization of the strain energy density in linear anisotropic elasticity

The problem considered here is that of extremizing the strain energy density of a linear anisotropic material by varying the relative orientation between a fixed stress state and a fixed material symmetry. It is shown that the principal axes of stress must coincide with the principal axes of strain in order to minimize or maximize the strain energy density in this situation. Specific conditions for maxima and minima are obtained. These conditions involve the stress state and the elastic constants. It is shown that the symmetry coordinate system of cubic symmetry is the only situation in linear anisotropic elasticity for which a strain energy density extremum can exist for all stress states. The conditions for the extrema of the strain energy density for transversely isotropic and orthotropic materials with respect to uniaxial normal stress states are obtained and illustrated with data on the elastic constants of some composite materials. Not surprisingly, the results show that a uniaxial normal stress in the grain direction in wood minimizes the strain energy in the set of all uniaxial stress states. These extrema are of interest in structural and material optimization.     

Specification of crack opening and surface sliding in orthotropic viscoelastic material

Modelling of crack opening and surface sliding in an orthotropic viscoelastic material is made by introducing two coefficients: one for the surface displacement and surface friction. The material possesses orthotropy in two dimensions and viscoelastic property consisting of a Kelvin element in series with a spring. The method of Laplace transform is applied to obtain a closed form solution to the problem. Explicit expressions of Mode I and II stress intensity factors are computed together with crack surface opening. Trade-off between the Mode I and II stress intensity factors depends on the nature of material orthotropy.     

Full field methods and residual stress analysis in orthotropic material. I Linear approach

This work analyzes the problem of residual stress determination in an orthotropic material using the hole drilling technique combined with non-contact, full field optical methods. Due to the complex behavior of the material, first a solution algorithm for the isotropic case is analyzed, then the procedure is extended to solving the more complex problem. In the first part of the work, the simplified Smith–Schajer solution to the through-hole problem for an orthotropic material is analyzed, showing that the same linear least square approach used in the isotropic case applies to a large set of orthotropic materials; based on this analysis a simple residual stress measurement algorithm is developed using either analytical or numerically estimated calibration coefficients.In the second part of the work, the general solution is discussed: since in this case the simplified Smith–Schajer solution cannot be used, the Lekhnitskii’s analysis of the through-hole plate in tension is introduced and extended to handle residual stresses. On this basis a solution algorithm using the nonlinear fit of the theoretical displacement field capable of treating all the orthotropic materials at the cost of a more complex numerical procedure is proposed. The performances of both algorithms are tested against numerically generated noisy fields and experimental ones and show a good reliability and accuracy.     

Acoustoelastic determination of stress in slightly orthotropic materials

A method is proposed to determine stresses in acoustoelasticity by making use of orthotropic stress-acoustic relations and the equations of equilibrium. It is derived theoretically that shear stress is determined ny ultrasonic data ofB and ?, which denote a magnitude of acoustic birefringence and its principal direction, respectively. Other stress components are obtained by numerical integration of the equilibrium equation with the shear stress thus determined. Experiments were carried out to show the validity and usefulment of the method. This method was applied to the measurement of stress field on a plate with a circular hole subjected to axial tension. Ultrasonic measurements were made by a Y-cut quartz transducer with 5-MHz fundamental frequency. The specimen was cut out from 1100 aluminum plate of 4-mm thickness, which shows a slight orthotropy due to roll working. The values ofB and ? were measured in both stressed and unstressed state. Then, stress distributions were determined by the method proposed here, and are compared with the known theoretical distributions.     

Interface crack problem in layered orthotropic materials under thermo-mechanical loading

In the present paper, the behavior of an interface crack for a homogeneous orthotropic strip sandwiched between two different functionally graded orthotropic materials subjected to thermal and mechanical loading is considered. It is assumed that interface crack is partly insulated, and the temperature drop across the crack surfaces is the result of the thermal resistance due to the heat conduction through the crack region. The elastic properties of the material are assumed to vary continuously along the thickness direction. The principal directions of orthotropy are parallel and perpendicular to the crack orientation. The complicated mixed boundary problems of equations of heat conduction and elasticity are converted analytically into singular integral equations, which are solved numerically. The main objective of the paper is to study the effects of material nonhomogeneity parameters and the dimensionless thermal resistance on the thermal stress intensity factors for the purpose of gaining better understanding of the thermal behavior of graded layer.     

Orthotropic elastic–plastic material model for paper materials

Paper and paperboard generally exhibit anisotropic and non-linear mechanical material behaviour. In this work, the development of an orthotropic elastic–plastic constitutive model, suitable for modelling of the material behaviour of paper is presented. The anisotropic material behaviour is introduced into the model by orthotropic elasticity and an isotropic plasticity equivalent transformation tensor. A parabolic stress–strain relation is adopted to describe the hardening of the material. The experimental and numerical procedures for evaluation of the required material parameters for the model are described. Uniaxial tensile testing in three different inplane material directions provides the calibration of the material parameters under plane stress conditions. The numerical implementation of the material model is presented and the model is shown to perform well in agreement with experimentally observed mechanical behaviour of paper.     

A large strain plasticity model for anisotropic materials — composite material application

In this work a generalized anisotropic model in large strains based on the classical isotropic plasticity theory is presented. The anisotropic theory is based on the concept of mapped tensors from the anisotropic real space to the isotropic fictitious one. In classical orthotropy theories it is necessary to use a special constitutive law for each material. The proposed theory is a generalization of classical theories and allows the use of models and algorithms developed for isotropic materials. It is based on establishing a one-to-one relationship between the behavior of an anisotropic real material and that of an isotropic fictitious one. Therefore, the problem is solved in the isotropic fictious space and the results are transported to the real field. This theory is applied to simulate the behavior of each material in the composite. The whole behavior of the composite is modeled by incorporating the anisotropic model within a model based on a modified mixing theory.     

Interaction of two off-axis cracks exhibiting orthotropic and general anisotropic behavior under arbitrary extension

Obtained in this work are the stress intensity factors for two interacting cracks, one of which coincides with the axis of material orthotropy and the other is oriented in a general direction. The direction of the applied stress can be varied with reference to the cracks and/or the axes of material anisotropy. Mode I and II stress intensity factors are displayed to show their variations with the geometric, material and load parameters.