ON A CLASS OF METHOD FOR SOLVING PROBLEMS WITH RANDOM BOUNDARY NOTCHES AND/OR CRACKS—(Ⅲ) COMPUTATIONS FOR BOUNDARY CRACKS

This paper continues the discussions to a class of method for solving problems withrandom boundary notches and for cracks in refs.[1] and [2].Using the method developed in[1].[2]with important modifications about inclusion of singularities in the formulation. wearrive at a very effective computational process for problems with random boundarycracks. Actual computations for boundary cracks with or without apptied tractions in theirsurfaces. Show that the present method is quite workable for the problems consideredwithin proper range of characteristic parameters. The results obtained here extend thecontents of “Handbook of Stress Intensity Factors” given by G. C. Sih.     

ON A CLASS OF METHOD FOR SOLVING PROBLEMS WITH RANDOM BOUNDARY NOTCHES AND/OR CRACKS (Ⅱ) COMPUTATIONS FOR BOUNDARY NOTCHES

This work is a continuation of the discussion of [1], "On a class of method for solving problems with random boundary notches and/or cracks, (Ⅰ)" by C. Ouyang (Appl. Math. & Mech.,Vol. 1, No. 2, 1980). Here computations for boundary notches are made by using the theory and formulas presented in [1]. In the computation modification is also made for the boundary conditions in parametric plane in [1]. Numerical results for examples show that within ranges of parameter considered in the paper, for example L, the present method in quite workable in practical computations.     

Thermoelastic analysis of multiple defects with the extended finite element method

In this paper, the extended finite element method (XFEM) is adopted to analyze the interaction between a sin-gle macroscopic inclusion and a single macroscopic crack as well as that between multiple macroscopic or micro-scopic defects under thermal/mechanical load. The effects of different shapes of multiple inclusions on the material thermomechanical response are investigated, and the level set method is coupled with XFEM to analyze the interaction of multiple defects. Further, the discretized extended finite element approximations in relation to thermoelastic prob-lems of multiple defects under displacement or temperature field are given. Also, the interfaces of cracks or materials are represented by level set functions, which allow the mesh assignment not to conform to crack or material interfaces. Moreover, stress intensity factors of cracks are obtained by the interaction integral method or the M-integral method, and the stress/strain/stiffness fields are simulated in the case of multiple cracks or multiple inclusions. Finally, some numer-ical examples are provided to demonstrate the accuracy of our proposed method.     

THE CONVERGENT CONDITION AND UNITED FORMULA OF STEP REDUCTION METHOD

The step reduction method was first suggested by Prof. Yeh Kai-yuan. This method has more advantages than other numerical methods. By this method, the analytic expression of solution can he obtained for solving nonuniform elastic mechanics. At the same time. its ealculuting time is very short and convergent speed very fast. In this paper, the convergent condition and nited formula of step reduction method are given by mathemutical method. It is proved that the solution of displacement and stress resultants obtained by this method can eonverge to exact solution uniformly, when the convergent condition is sutisfied. By united formula, the analytic solution solution can be expressed as matrix form, and therefore the former complicated expression can be avoMed. Two numerical examples are given at the end of this paper which indicate that. by the theory in this paper, a right model can be obtained for step reduction method.     

The Solution for Axisymmetrical Shells with Abrupt Curvature Change(Corrugated Shells)by the Finite Element Method

It has been noted in the present paper that the finite element method using linear elements for solving axisymmetrical shells cannot be applied to the analysis of axisymmetrical shells with abrupt curvature change,owing to the fact that the influence of the curvature upon the angular displace-ments has been neglected.The present paper provides a finite element method using linear elements in which the influence of curvature is considered and the angular displacements are treated as continuous parameters.This method has been applied to the calculation of corrugated shells of the type C,and compared with the experimental results obtained by Turner-Ford as well as with the analytical solution given by Prof.Chien Wei-zang.The compari-sons have been proved that this theory is correct.     

A SEMI－ANALYTICAL AND SEMI－NUMERICAL METHOD FOR SOLVING 2－D SOUND－STRUCTURE INTERACTION PROBLEMS

Based on the transfer matrix method and the virtual source simulation technique,this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound-structure interaction problems under a harmonic excitation. Within any integration segment,as long as its length is small enough, along the circumferential curvilinear coordinate, the nonhomogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper. For the exterior fluid domain, the multi-circular virtual source simulation technique is adopted. The source density distributed on each virtual circular curve may be ex-panded as the Fourier‘‘s series. Combining with the inverse fast Fourier transformation, a higher accuracy and efficiency method for solving 2-D exterior Helmholtz‘s problems is presented in this paper. In the aspect of solution to the coupling equations, the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper. According to the superposition principle and compatibility conditions at the interface between the elastic ring and fluid, the algebraic equation of system can be directly constructed by using the least square approximation. Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented. Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.     

Application of subregion function method for solving beam-board structure

In this paper,the solution to the structure consisting of a bead and a board is given as aresult of the application of the subregion function method which was suggested in ref.[1].The same problem is also computed with finite element method.The comparison betweenthe two results shows that the application of the subregion function in the method of weightedresiduals is practical and effective,especially for solving compound structures.     

NEW BOUNDARY ELEMENT METHOD FOR TORSION PROBLEMS OF CYLINDER WITH CURVILINEAR CRACKS

The Saint-Venant torsion problems of a cylinder with curvilinear cracks were considered and reduced to solving the boundary integral equations only on cracks. Using the interpolation models for both singular crack tip elements and other crack linear elements, the boundary element formulas of the torsion rigidity and stress intensity factors were given. Some typical torsion problems of a cylinder involving a straight, kinked or curvilinear crack were calculated. The obtained results for the case of straight crack agree well with those given by using the Gauss-Chebyshev integration formulas, which demonstrates the validity and applicability of the present boundary element method.     

The solution of rectangular plates with large deflection by spline functions

In this paper, Von Karman ’s set of nonlinear equations for rectangular plates with large deflection is divided into several sets of linear equations by perturbation method, the dimensionless center deflection being taken as a perturbation parameter. These sets of linear equations are solved by the spline finite-point (SFP) method and by the spline finite element (SFE) method. The solutions for rectangular plates having any length-to-width ratios under a uniformly distributed load and with various boundary conditions are presented, and the analytical formulas for displacements and deflections are given in the paper. The computer programs are worked out by ourselves. Comparison of the results with those in other papers indicates that the results of this paper are satisfactorily better.     

THE PERTURBATION FINITE ELEMENT METHOD FOR SOLVING THE PLANE PROBLEM IN CONSIDERATION OF LINEAR CREEP

In this paper, amethod (PFMC) for solving plane problem of linear creep is presented by using perturbation finite element. It can be used in plane problem in consideration of creep, such as reinforced concrete beam, prestressed concrete beam, reinforced concrete cylinder and reinforced concrete tunnel in elastic or visco-elastic medium, as well as underground building and so on. In the presented method, the assumption made in the general increment method that variables remain constant in a divtded time interval is not taken. The accuracy is improved and the length of time step becomes larger. The computer storage can be reduced and the calculating efficiency can be increased. Perturbation finite element formulae for four-node quadrilateral isoparametric element including reinforcement are established and five numerical examples are given. As contrasted with the analytical solution, the accuracy is satisfactory.     

GENERAL FORM OF MATCHING EQUATION OF ELASTIC-PLASTIC FIELD NEAR CRACK LINE FOR MODEⅠCRACK UNDER PLANE STRESS CONDITION

Crack line field analysis method has become an independent method for crack elastic-plastic analysis, which greatly simplifies the complexity of crack elastic-plastic problems and overcomes the corresponding mathematical difficulty. With this method, the precise elastic-plastic solutions near crack lines for variety of crack problems can be obtained. But up to now all solutions obtained by this method were for different concrete problems, no general steps and no general form of matching equations near crack line are given out. With crack line analysis method, this paper proposes the general steps of elastic-plastic analysis near crack line for mode Ⅰ crack in elastic-perfectly plastic solids under plane stress condition, and in turn given out the solving process and result for a specific problem.     

APPLICATION OF THE METHOD OF SPLIT RIGIDITIES TO ANlSOTROPIC LAMINATED SHALLOW SHELLS

In this paper,according to the method stated by Hu Hai-chang in [3],on the basis of [1],the method of split rigi-dities is generalized for the purpose of solving the problemsof lateral deflection,stability and lateral vibration foranisotropic laminated shallow shells.and a simple and prac-tical approximate method is obtained.in which the errors andcomputing work are comparatively small.     

Precise integration method for a class of singular two-point boundary value problems

In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations(ODEs) are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.     

SOLVE A CONTACT PROBLEM BY OPTIMIZATION METHOD—TO CALCULATE THE STRESSES FOR THE SOFTWHEEL IN A HARMONIC GEAR

It is difficult to solve the contact problem by usual finite element program. In this paper, we express the contact problem as an optimization problem. In this form we do not need to know all boundary condition in advance. We only need to know the constraint conditions. This method is especially good for solving contact problem. Using this method, we calculate the stresses of the softwheel in the harmonic gear given by Shanghai Jiaotong University, and the results are in good agreement with the experimental results.     

THE ASYMPTOTIC ANALYEICAL SOLUTION OF THE STRAINHARDENING ELASTO-PLASTIC PLATE CONTAINING A CIRCULAR HOLE UNDER SIMPLE TENSION

In this paper the general asymptotic analytical solution ofplane problem of elasto-plasticity with strain-hardening[2]isused in solving the problem of an infinitely large plate con-taining a circular hole under simple tension,and the analy-tical expressions of stress components of the first two approxi-mations,are given,These results are compared with the numeric-al and the experimental results given by other authors[4,5],and a good agreement is obtained.At the end of this paper theauthors inspect the correctness of Neuber’s formula[9]for thisproblem.     

The method of variation on parameters for integration of a generalized Birkhoffian system

This paper focuses on studying the integration method of a generalized Birkhoffian system.The method of variation on parameters for the dynamical equations of a generalized Birkhoffian system is presented.The procedure for solving the problem can be divided into two steps:the first step,a system of auxiliary equations is constructed and its general solution is given;the second step,the parameters are varied,and the solution of the problem is obtained by using the properties of generalized canonical transformation.The method of variation on parameters for the generalized Birkhoffian system is of universal significance,and we take a nonholonomic system and a nonconservative system as examples to illustrate the application of the results of this paper.     

The Schrödinger equation in theory of plates and shells with orthorhombic anisotropy

This work is the continuation of the discussion of Refs. [1-5]. In this paper:[A] The Love-Kirchhoff equations of vibration problem with small deflection for orthorhombic misotropic thin shells or orthorhombic anisotropic thin plates on Winkler’s base are classified as several of the same solutions of Schrodmger equation, and we can obtain the general solutions for the two above-mentioned problems by the method in Refs. [1] and [3-5].[B] The. von Karman-Vlasov equations of large deflection problem for shallow shells with orthorhombic anisotropy (their special cases are the von Harmon equations of large deflection problem for thin plates with orthorhombic anisotropy) are classified as the solutions of AKNS equation or Dirac equation, and we can obtain the exact solutions for the two abovementioned problems by the inverse scattering method in Refs. [4-5].The general solution of small deflection problem or the exact solution of large deflection problem for the corrugated or rib-reinforced plates and shells as special cases is included in this paper.     

THE SOLUTION OF RECTANGULAR PLATES WITH LARGE DEFLECTION BY SPLINE FUNCTIONS

In this paper,Von Kármán’s set of nonlinear equations for rectangular plates withlarge deflection is divided into several sets of linear equations by perturbation method,thedimensionless center deflection being taken as a perturbation parameter.These sets of linearequations are solved by the spline finite-point(SFP)method and by the spline finiteelement(SFE)method.The solutions for rectangular plates having any length-to-widthratios under a uniformly distributed load and with various boundary conditions arepresented,and the analytical formulas for displacements and deflections are given in thepaper.The computer programs are worked out by ourselves.Comparison of the results withthose in other papers indicates that the results of this paper are satisfactorily better.     

A <Emphasis Type="Italic">VU</Emphasis>-decomposition method for a second-order cone programming problem

A vu-decomposition method for solving a second-order cone problem is presented in this paper. It is first transformed into a nonlinear programming problem. Then, the structure of the Clarke subdifferential corresponding to the penalty function and some results of itsvu-decomposition are given. Under a certain condition, a twice continuously differentiable trajectory is computed to produce a second-order expansion of the objective function. A conceptual algorithm for solving this problem with a superlinear convergence rate is given.     

AN IMPROVED HYBRID BOUNDARY NODE METHOD IN TWO-DIMENSIONAL SOLIDS

The hybrid boundary node method （HBNM） is a promising method for solving boundary value problems with the hybrid displacement variational formulation and shape functions from the moving least squares（MLS） approximation. The main idea is to reduce the dimensionality of the former and keep the meshless advantage of the latter. Following its application in solving potential problems, it is further developed and numerically implemented for 2D solids in this paper. The rigid movement method is employed to solve the hyper-singular integrations. Numerical examples for some 2D solids have been given to show the characteristics. The computation results obtained by the present method are in excellent agreement with the analytical solution. The parameters that influence the performance of this method are studied through numerical examples.