Crack tip field andJ-integral with strain gradient effect

The mode I plane strain crack tip field with strain gradient effects is presented in this paper based on a simplified strain gradient theory within the framework proposed by Acharya and Bassani. The theory retains the essential structure of the incremental version of the conventionalJ ₂ deformation theory. No higher-order stress is introduced and no extra boundary value conditions beyond the conventional ones are required. The strain gradient effects are considered in the constitutive relation only through the instantaneous tangent modulus. The strain gradient measures are included into the tangent modulus as internal parameters. Therefore the boundary value problem is the same as that in the conventional theory. Two typical crack problems are studied: (a) the crack tip field under the small scale yielding condition induced by a linear elastic mode-IK-field and (b) the complete field for a compact tension specimen. The calculated results clearly show that the stress level near the crack tip with strain gradient effects is considerable higher than that in the classical theory. The singularity of the strain field near the crack tip is nearly equal to the square-root singularity and the singularity of the stress field is slightly greater than it. Consequently, theJ-integral is no longer path independent and increases monotonically as the radius of the calculated circular contour decreases. The project supported by the National Natural Science Foundation of China (19704100 and 10202023) and the Natural Science Foundation of Chinese Academy of Sciences (KJ951-1-20)     

Two dimensional stress intensity factor computations by multi-domain boundary element method for crack closure with friction

A multi-domain boundary element method is used to compute the stress intensity factor of plane stress/plane strain crack problems with friction. The analysis is performed by using traction-singular quarter-point boundary elements on each side of the crack tips. The increment iteration is given. The technique is applied to some specific examples in order to show that the results will be with good accuracy.The Project 13 supported by National Natural Science Foundation of China.     

Crack problems of piezoelectric materials

This paper presents an analysis of crack problems in homogeneous piezoelectrics or on the interfaces between two dissimilar piezoelectric materials based on the continuity of normal electric displacement and electric potential across the crack faces. The explicit analytic solutions are obtained for a single crack in piezoelectrics or on the interfaces of piezoelectric bimaterials. A class of boundary problems involving many cracks is also solved. For homogeneous materials it is found that the normal electric displacementD ₂ induced by the crack is constant along the crack faces which depends only on the applied remote stress field. Within the crack slit, the electric fields induced by the crack are also constant and not affected by the applied electric field. For the bimaterials with realH, the normal electric displacementD ₂ is constant along the crack faces and electric fieldE ₂ has the singularity ahead of the crack tip and a jump across the interface. The project is supported by the National Natural Science Foundation of China(No. 19704100) and the Natural Science Foundation of Chinese Academy of Sciences(No. KJ951-1-201).     

Positive solutions to （n-1,1） m-point boundary value problems with dependence on the first order derivative

Using the extension of Krasnoselskii＇s fixed point theorem in a cone, we prove the existence of at least one positive solution to the nonlinear nth order m-point boundary value problem with dependence on the first order derivative. The associated Green＇s function for the nth order m-point boundary value problem is given, and growth conditions are imposed on the nonlinear term f which ensures the existence of at least one positive solution. A simple example is presented to illustrate applications of the obtained results.     

Weight functions for interface cracks in dissimilar anisotropic materials

Bueckner‘s work conjugate integral customarily adopted for linear elastic materials is established for an interface crack in dissimilar anisotropic materials. The difficulties in separating Stroh‘s six complex arguments involved in the integral for the dissimilar materials are overcome and then the explicit function representations of the integral are given and studied in detail. It is found that the pseudo-orthogonal properties of the eigenfunction expansion form (EEF) for a crack presented previously in isotropic elastic cases, in isotopic bimaterial cases, and in orthotropic cases are also valid in the present dissimilar arbitrary anisotropic cases. The relation between Bueckner‘s work conjugate integral and the J-integral in these cases is obtained by introducing a complementary stressdisplacement state. Finally, some useful path-independent integrals and weight functions are proposed for calculating the crack tip parameters such as the stress intensity factors.     

Further investigation of interaction between interface macrocrack and parallel microcracks in bimaterial anisotropic solids

In this paper, with the aid of superimposing technique and the Pseudo Traction Method (PTM), the interaction problem between an interface macrocrack and parallel microcracks in the process zone in bimaterial anisotropic solids is reduced to a system of integral equations. After the integral equations are solved numerically, a conservation law among three kinds ofJ-integrals is obtained which are induced from the interface macrocrack tip, the microcrack and the remote field, respectively. This conservation law reveals that the microcrack shielding effect in such materials could be considered as the redistribution of the remoteJ-integral. The project supported by the National Natural Science Foundation of China, and the Doctorate Foundation of Xi'an Jiaotong University     

Strong convergence theorem for a generalized equilibrium problem and a <Emphasis Type="Italic">k</Emphasis>-strict pseudocontraction in Hilbert spaces

The main purpose of this paper is to study a new iterative algorithm for finding a common element of the set of solutions for a generalized equilibrium problem and the set of fixed points for a k-strict pseudocontractive mapping in the Hilbert space. The presented results extend and improve the corresponding results reported in the lit-erature.     

Conservation Law and Application of J-Integral in Multi-Materials

The conservation law of J-integral in two-media with a crack paralleling to the interface of the two media was firstly proved by analytical and numerical finite element method. Then a schedule model was established that an interface crack is inserted in four media. According to the J-integral conservation law on multi-media, the energy release ratio of Ⅰ-type crack was considered to be conservation when the middle medium layers are very thin. And the conservation law was also convinced by numerical method. By means of the dimension analysis on the model, the asymptotic results and formula calculating the energy release ratio and complex stress intensity factor are presented.     

Lower bound limit analysis of three-dimensional elastoplastic structures by boundary element method

IntroductionThelimitanalysisofstructuresisoneofthemostpracticalandusefulbranchesinplasticity .Ithasimportantapplicationbackgroundforproblemssuchasthedeterminationofloadcarryingcapacityandplasticformingofmetal.Thepurposeofthelimitanalysisofstructuresistoprovidereliabletheoreticalbasesforengineeringdesignandsafetyassessment.Asasimplifiedmethodforelastoplasticproblems,limitanalysisneednotrequirethehistoryofloadandcancomputethelimitloadsdirectlyinsteadofelastoplasticincrementalcomputationwhichisus…     

A viscoplastic model based on no-yield-surface concept

The elastic/viscoplastic constitutive equation which describes deformation law of metal materials was suggested based on no-yield-surface concept and thermal activation theory of dislocation. The equation which takes account of effects of strain-rate, strain history, strain-rate history, hardening and temperature has stronger physical basis. Comparison of the theoretical prediction with experimental results of mechanical behaviours of Ti under conditions of uniaxial stress and room temperature shows good consistency.The Projects Supported by National Natural Science Foundation of China     

用于设计灵敏度分析的弹粘塑性边界元法

基于一致切线算子概念的弹粘塑性隐式边界元方法,进行了弹粘塑性设计灵敏度分析.采用了Perzyna经典粘塑性本构模型,针对包含各向同性硬化和运动硬化的混合硬化模型,导出了弹粘塑性灵敏度分析的径向返回算法和一致切线算子.利用直接微分的方法,建立了设计灵敏度分析的弹粘塑性边界元增量方程,导出了弹粘塑性径向返回的灵敏度公式.给出了在不同粘塑性流动参数下的三个典型算例的结果,与ANSYS结果相吻合,证明了方法是正确的.     

Implicit BEM formulations for usual and sensitivity problems in elasto-plasticity using the consistent tangent operator concept

This paper presents boundary element method (BEM) formulations for usual and sensitivity problems in (small strain) elasto-plasticity using the concept of the local consistent tangent operator (CTO). “Usual” problems here refer to analysis of nonlinear problems in structural and solid continua, for which Simo and Taylor first proposed the use of the CTO within the context of the finite element method (FEM). A new implicit BEM scheme for such problems, using the CTO, is presented first. A formulation for sensitivity analysis follows. It is shown that the sensitivity of the strain increment, associated with an infinitesimal variation of some design parameter, solves a linear problem which is governed by the (converged value of the) same global CTO as the one that appears in the usual problem. Numerical results for both usual and sensitivity problems are shown for a one-dimensional example. They demonstrate the effectiveness of the present approach. In particular, accurate sensitivities with respect to material parameters (e.g., exponent of the power-type hardening law) are obtained even with few integration cells and for large load increments.     

<Emphasis Type="Italic">J</Emphasis><Subscript>2</Subscript> invariant relative orbits via differential correction algorithm

This paper describes a practical method for finding the invariant orbits in J ₂ relative dynamics. Working with the Hamiltonian model of the relative motion including the J ₂ perturbation, the effective differential correction algorithm for finding periodic orbits in three-body problem is extended to formation flying of Earth’s orbiters. Rather than using orbital elements, the analysis is done directly in physical space, which makes a direct connection with physical requirements. The asymptotic behavior of the invariant orbit is indicated by its stable and unstable manifolds. The period of the relative orbits is proved numerically to be slightly different from the ascending node period of the leader satellite, and a preliminary explanation for this phenomenon is presented. Then the compatibility between J ₂ invariant orbit and desired relative geometry is considered, and the design procedure for the initial values of the compatible configuration is proposed. The influences of measure errors on the invariant orbit are also investigated by the Monte–Carlo simulation. The project supported by the Innovation Foundation of Beihang University for Ph.D. Graduates, and the National Natural Science Foundation of China (60535010).     

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.     

Boundary element method to solve torsion problem of cracked cylinder

Separating the discontinuous solution by use of the single crack solution, together with the regular solution of harmonic function, the torsion problem of a cracked cylinder is reduced to solving a set of mixed-type integral equations and its numerical technique is then proposed by combining the numerical method of singular integral equation with the boundary element method. Several numerical examples are calculated which will be useful to engineering practice. The method proposed is characterized by its fine accuracy and convenience for using, which can be extended to the cases of multiple crack.The project supported by National Natural Science Foundation of China.     

Computer simulated analyses on deformation and fracture of non-cracked and cracked specimens

A unified damage and fracture model, the combinatory work density model, which is suitable for either non-cracked body or cracked body has been suggested^[t−7]. In the present paper, the deformation and fracture of the two kinds of tensile spceimen and TPB specimen made of 40Cr steel have been simulated by using the new model together with the large elastic-plastic deformation finite element method. The results give a good picture of the whole deformation and fracture processes of the specimens in experiments; especially, the results on the TPB specimen can be used to obtain the relationship between load and displacement at the loading pointP-Δ, and between crack extension and displacement at the loading point Δa-Δ, the resistance curveJ _R -Δa and the fracture toughnessJ _1C . All the results are in remarkable agreement with those obtained by experiments. Therefore the model suggested here can be used to simulate crack initiation and propagation in non-cracked body and fracture initiation and crack stable propagation in cracked body. The project supported by National Natural Science Foundation of China     

Stress analysis for an infinite strip weakned by periodic cracks

Stress analysis for an infinite stripcracks were assumed in a horizontal position,weakened by periodic cracks is studied. The and the strip was applied by tension “p“ in y-direction. The boundary value problem can be reduced into a complex mixed one. It is found that the EEVM ( eigenfunction expansion variational method) is efficient to solve the problem. The stress intensity factor at the crack tip and the T-stress were evaluated. From the deformation response under tension the cracked strip can be equivalent to an orthotropic strip without cracks. The elastic properties in the equivalent orthotropic strip were also investigated. Finally, numerical examples and results were given.     

Fixed point indexes and its applications to nonlinear integral equations modelling infectious diseases

In this paper the fixed point index problem for a class of positive operators with boundary control conditions is discussed, and some sufficient conditions for the fixed point index to be equal to 1 or 0 are given. Moreover, a general fixed point theorem of expansions and compressions for cone is obtained, which generalizes and improves the corresponding results of [3,8,9]. As an application, we utilize the results presented above to study the existence conditions of positive solutions of nonlinear integral equations modelling infectious diseases.Project supported by National Natural Science Foundation of China.     

Boundary element method for buckling eigenvalue problem and its convergence analysis

The conditions for determining solution of buckling eigenvalue problem are discussed. The corresponding system of integral equations with constraint conditions and boundary variational equations with Lagrange multiplier are established. The theorems on the existence and uniqueness of the solution for these problems are given. The corresponding boundary element method is constructed and the error estimation for the approximation solution is obtained. Finally the numerical example is given. Foundation item: the National Natural Science Foundation Pre-research Project (T4107015) Biography: Ding Rui (1969-)