Phenomenological invariants and their application to geometrically nonlinear formulation of triangular finite elements of shear deformable shells

A phenomenological definition of classical invariants of strain and stress tensors is considered. Based on this definition, the strain and stress invariants of a shell obeying the assumptions of the Reissner–Mindlin plate theory are determined using only three normal components of the corresponding tensors associated with three independent directions at the shell middle surface. The relations obtained for the invariants are employed to formulate a 15-dof curved triangular finite element for geometrically nonlinear analysis of thin and moderately thick elastic transversely isotropic shells undergoing arbitrarily large displacements and rotations. The question of improving nonlinear capabilities of the finite element without increasing the number of degrees of freedom is solved by assuming that the element sides are extensible planar nearly circular arcs. The shear locking is eliminated by approximating the curvature changes and transverse shear strains based on the solution of the Timoshenko beam equations. The performance of the finite element is studied using geometrically linear and nonlinear benchmark problems of plates and shells.     

The scattering of harmonic elastic anti-plane shear waves by a finite crack in infinitely long strip using the non-local theory

In this paper, the scattering of harmonic anti-plane shear waves by a finite crack in infinitely long strip is studied using the non-local theory. The Fourier transform is applied and a mixed boundary value problem is formulated. Then a set of dual integral equations is solved using the Schmidt method instead of the first or the second integral equation method. A one-dimensional non-local kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress occurring at the crack tips. Contraty to the classical elasticity solution, it is found that no stress singularity is present at the crack tip. The non-local dynamic elastic solutions yield a finite hoop stress at the crack tip, thus allowing for a fracture criterion based on the maximum dynamic stress hypothesis. The finite hoop stress at the crack tip depends on the crack length, the width of the strip and the lattice parameter. Supported by the Post Doctoral Science Foundation of Heilongjiang Province, the Natural Science Foundation of Heilongjiang Province and the National Foundation for Excellent Young Investigators.     

Closed-form solution for two collinear mode-III cracks in an orthotropic elastic strip of finite width

The problem of an orthotropic strip containing two collinear cracks normal to the strip boundaries is considered. The Fourier series method is used to reduce the associated boundary value problem to triple series equations, then to a singular integral equation, which can be solved analytically. Under remote uniform antiplane shear loading, the stress field and the crack sliding displacement are determined analytically and stress intensity factors are also given in a closed form.     

Stress intensity factors for an edge interface crack in a bonded semi-infinite plate for arbitrary material combination

Although a lot of interface crack problems were previously treated, few solutions are available under arbitrary crack lengths and material combinations. In this paper the stress intensity factors of an edge interface crack in a bonded strip are considered under tension with varying the crack length and material combinations systematically. Then, the limiting solutions are provided for an edge interface crack in a bonded semi-infinite plate under arbitrary material combinations. In order to calculate the stress intensity factors accurately, exact solutions in an infinite bonded plate are also considered to produce proportional singular stress fields in the analysis of FEM by superposing specific tensile and shear stresses at infinity. The details of this new numerical solution are described with clarifying the effect of the element size on the stress intensity factor. It is found that for the edge interface crack the normalized stress intensity factors are not always finite depending upon Dunders’ parameters. This behavior can be explained from the condition of the singular stress at the end of bonded strip. Convenient formulas are also given by fitting the computed results.     

双剪统一弹塑性有限差分方法研究

基于拉格朗日有限差分方法,建立了双剪统一弹塑性有限差分计算格式,并利用VC++语言编写动态链接库文件将双剪统一弹塑性模型导入拉格朗日有限差分程序FLAC(Fast Lagrangian Analysis of Continua)中进行计算分析。双剪统一弹塑性有限差分方法可以模拟复杂应力状态下结构的渐进破坏,无需形成刚度矩阵,对于材料非线性问题无需进行迭代计算,因此在理论和工程应用中都有积极的意义。本文利用双剪统一弹塑性有限差分方法对拉压强度不等材料的厚壁圆筒受内压、中心带孔板条受拉压、条形基础下的地基极限分析及边坡问题进行了数值分析并与滑移线场等解析方法计算结果进行对比,结果均吻合较好。     

有限宽板中心裂纹断裂过程区的长度和位移公式

对岩石、混凝土和陶瓷等准脆性材料进行断裂分析,有必要研究裂纹前端的断裂过程区所起的作用以及各种因素对它的影响。本文提出分析有限宽中心裂纹板剪切断裂过程区的方法,此方法基于D－B模型的叠加原理,考虑了压剪断裂的摩擦阻力,并将有限板宽影响简化为载荷的修正,以非常简便的方法推导出计算断裂过程区长度和位移的公式,这些公式补充了现有文献只有“无限大”板的解的不足。利用这些公式,分析各种参数对断裂过程区长度和位移的影响变得非常方便。     

ANTI-PLANE FRACTURE ANALYSIS OF FUNCTIONALLY GRADIENT MATERIAL INFINITE STRIP WITH FINITE WIDTH

The special case of a crack under mode III conditions was treated, lying parallel to the edges of an infinite strip with finite width and with the shear modulus varying exponentially perpendicular to the edges. By using Fourier transforms the problem was formulated in terms of a singular integral equation. It was numerically solved by representing the unknown dislocation density by a truncated series of Chebyshev polynomials leading to a linear system of equations. The stress intensity factor (SIF) results were discussed with respect to the influences of different geometric parameters and the strength of the non-homogeneity. It was indicated that the SIF increases with the increase of the crack length and decreases with the increase of the rigidity of the material in the vicinity of crack. The SIF of narrow strip is very sensitive to the change of the non-homogeneity parameter and its variation is complicated. With the increase of the non-homogeneity parameter, the stress intensity factor may increase, decrease or keep constant, which is mainly determined by the strip width and the relative crack location. If the crack is located at the midline of the strip or if the strip is wide, the stress intensity factor is not sensitive to the material non-homogeneity parameter.     

功能梯度条共线Griffith裂纹反平面剪切冲击

研究正交各向异性功能梯度条中多个共线Griffith裂纹的反平面剪切冲击问题．材料两个方向的剪切模量假定按比例同时以特定的梯度变化．采用Laplace和Fourier变换及引进位错密度函数将问题化为求解Cauchy奇异积方程，进而化为代数方程数值求解．考查材料非均匀性、正交性和功能梯度条高度对裂尖动态断裂特性的影响．动应力强度因子的数值结果显示：增加剪切模量的梯度和(或)增加垂直于裂纹面方向的剪切模量，可以抑制动应力强度因子的幅度；若功能梯度条较薄，增大条形域的高度也可抑制裂纹扩展．     

A MOVING CRACK IN A NONHOMOGENEOUS MATERIAL STRIP

This paper considers an anti-plane moving crack in a nonhomogeneous material strip of finite thickness. The shear modulus and the mass density of the strip are considered for a class of functional forms for which the equilibrium equation has analytical solutions. The problem is solved by means of the singular integral equation technique. The stress field near the crack tip is obtained. The results are plotted to show the effect of the material non-homogeneity and crack moving velocity on the crack tip field. Crack bifurcation behaviour is also discussed. The paper points out that use of an appropriate fracture criterion is essential for studying the stability of a moving crack in nonhomogeneous materials. The prediction whether the unstable crack growth will be enhanced or retarded is strongly dependent on the type of the fracture criterion used. Based on the analysis, it seems that the maximum 'anti-plane shear' stress around the crack tip is a suitable failure criterion for moving cracks in nonhomogeneous materials.     

A finite deformation theory of higher-order gradient crystal plasticity

For higher-order gradient crystal plasticity, a finite deformation formulation is presented. The theory does not deviate much from the conventional crystal plasticity theory. Only a back stress effect and additional differential equations for evolution of the geometrically necessary dislocation (GND) densities supplement the conventional theory within a non-work-conjugate framework in which there is no need to introduce higher-order microscopic stresses that would be work-conjugate to slip rate gradients. We discuss its connection to a work-conjugate type of finite deformation gradient crystal plasticity that is based on an assumption of the existence of higher-order stresses. Furthermore, a boundary-value problem for simple shear of a constrained thin strip is studied numerically, and some characteristic features of finite deformation are demonstrated through a comparison to a solution for the small deformation theory. As in a previous formulation for small deformation, the present formulation applies to the context of multiple and three-dimensional slip deformations.     

Crack problem for superconducting strip with finite thickness

The inclined crack problems are considered for a thin strip and a strip with finite thickness in a perpendicular magnetic field. The critical current density is assumed to be a constant. The crack orientation is varied and the effect of crack on the magnetic field distribution is neglected. Based on the analytical results and variational inequality, the field and current distributions are computed for both thin strip and strip with finite thickness cases, respectively. Then, the stress intensity factors at the crack tip are determined using the finite element method for magnetic field loads. The numerical results are presented for different inclined crack angles, magnetization processes and geometry parameters of the strip. The results show that the fracture behavior of the strip with finite thickness is more complicated than that of the thin strip. With the numerical results, we can predict the largest possibility of cracking as the strip is in an external field.     

Conservation properties for plane deformations of anisotropic linearly elastic curved strips

Plane deformations of a curved strip, composed of an homogeneous cylindrically anisotropic linearly elastic material, are considered. The strip is in equilibrium under the action of end loads, with the lateral sides traction-free. Two conservation properties for certain cross-sectional stress measures are established, generalizing previously known results for the case of a rectangular strip. Such conservation properties are useful in assessing the influence of material anisotropy on Saint-Venant's principle, as well as in establishing convexity properties for cross-sectional stress measures. In particular, it is anticipated that the results should be useful in determining the extent of edge effects in the testing of anisotropic and composite curved strips.     

断裂过程区D—B模型的叠加原理和截荷修正法

提出分析中心裂纹板剪切断裂的断裂过程区D－B模型的叠加原理，考虑了压剪断裂的摩擦阻力和体现有限权宽影响的载荷修正，以非常简便的方法推导出计算断裂过程区长度和位移的公式。     

An anti-plane propagating crack in a functionally graded piezoelectric strip bonded to a homogeneous piezoelectric strip

Summary A finite crack propagating at constant speed in a functionally graded piezoelectric strip (FGPS) bonded to a homogeneous piezoelectric strip is considered. It is assumed that the electroelastic material properties of the FGPS vary exponentially across the thickness of the strip, and that the bimaterial strip is under combined anti-plane mechanical shear and in-plane electrical loads. The analysis is conducted for the electrically unified crack boundary condition, which includes both the traditional permeable and the impermeable ones. By using the Fourier transform, the problem is reduced to the solution of Fredholm integral equations of the second kind. Numerical results for the stress intensity factor and the crack sliding displacement are presented to show the influences of the crack propagation speed, electric loads, FGPS gradation, crack length, electromechanical coupling coefficient, properties of the bonded homogeneous piezoelectric strip and crack location.     

Analysis of simple shear endochronic equations for finite plastic deformation

IntroductionThestress_strainbehaviorofmaterialswithfiniteplasticdeformationisaninterestingissue ,onwhichsignificantprogresshasbeenmadethroughboththephenomenologicalandphysicalapproaches.Thephenomenologicalapproachisbasedoncontinuummechanicsofplasticity .Ithasitsadvantageinsolvingcomplicatedproblemsbecauseofitssimplicity .Mostofphenomenologicaltheoriesareinvolvedintheconceptofcorotationalrates.Thematerialderivativeofstresswasnotobjectiveunderfinitedeformation .TheJaumannratewasusuallyusedbefo…     

Investigation of the Dynamic Behavior of a Griffith Crack in a Piezoelectric Material Strip Subjected to the Harmonic Elastic Anti-plane Shear Waves by Use of the Non-local Theory

In this paper, the dynamic behavior of a Griffith crack in a piezoelectric material strip subjected to the harmonic anti-plane shear waves is investigated by use of the non-local theory for impermeable crack surface conditions. To overcome the mathematical difficulties, a one-dimensional non-local kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress and the electric displacement near at the crack tip. By means of the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations. These equations are solved using the Schmidt method. Contrary to the classical solution, it is found that no stress and electric displacement singularity is present near the crack tip. The non-local dynamic elastic solutions yield a finite hoop stress near the crack tip, thus allowing for a fracture criterion based on the maximum dynamic stress hypothesis. The finite hoop stress at the crack tip depends on the crack length, the thickness of the strip, the circular frequency of incident wave and the lattice parameter.     

Rheological predictions of a transversely isotropic fluid model with extensible microstructure

A continuum extensible director theory was formulated to describe the isothermal, incompressible flow of uniaxial rodlike semiflexible liquid crystalline polymers. The model is strictly restricted to material that flow-align in shear, and that, in the absence of flow, are sufficiently far from the nematic-isotropic phase transition. The microstructure of the continuum is described by a variable length director, but the extensibility is finite. The model is an extension of the TIF (Transversely Isotropic Fluid) model of Ericksen (1960). The thermodynamic restrictions on the model parameters are found using the non-negative definiteness of the entropy production. The rheological material functions predicted by the model are calculated for steady simple shear and steady uniaxial extensional flows. In the rigid rod limit the model predictions agree with those of the TIF model, and for the finite extensibility case the model predictions are in agreement with those associated with flexible isotropic polymers: strong non-Newtonian shear viscosity, positive first normal stress differences, recoverable shear of order one, negative second normal stress differences, and a maximum in the steady uniaxial extensional viscosity.     

Finite anti-plane shear of an infinite slab with a traction-free elliptical cavity: Bounds on the stress concentration factor

This paper provides an analytical approach for obtaining bounds on elastic stress concentration factors in the theory of finite anti-plane shear of homogeneous, isotropic, incompressible materials. The problem of an infinite slab with traction-free elliptical cavity subject to a state of finite simple shear deformation is considered. Explicit estimates are obtained for the maximum shear stress in terms of the cavity geometry, applied stress at infinity and constitutive parameters. The analysis is based on application of maximum principles for second-order quasilinear uniformly elliptic equations.     

Anti-plane analysis of a functionally graded strip with multiple cracks

The stress fields are obtained for a functionally graded strip containing a Volterra screw dislocation. The elastic shear modulus of the medium is considered to vary exponentially. The stress components exhibit Cauchy as well as logarithmic singularities at the dislocation location. The dislocation solution is utilized to formulate integral equations for the strip weakened by multiple smooth cracks under anti-plane deformation. Several examples are solved and stress intensity factors are obtained.     

Interfacial stress transfer behavior in a specially-shaped fiber/matrix pullout test

Specially designed fibers are widely used in engineering practice because the specially-designed shape can help to improve the bonding strength of the fiber and the interface. Studied in this paper is the interfacial shear stress transfer behavior on both sides of the specially designed fiber when it is being pulled out; in which automatic analysis of three-dimensional photoelasticity is employed and the finite element method is adopted. The results show that the stress transfer occurs mainly in the region near the fiber＇s embedded end where the stress reaches its critical point, leading to debonding of the interface. Before debonding, as the pullout loading increases, the peak value of shear stress transfers along the fiber from the embedded end to the interior of the matrix, and then stops at the hooked part of the fiber because of its impediment. When the interface begins to debond as the load increases, the shear stress can be transferred to the hooked part.