Stiffness and Shear Stress Distribution of Glulam Beams in Elastic-Plastic Stage:Theory,Experiments and Numerical Modelling

Traditional methods focus on the ultimate bending moment of glulam beams and the fracture failure of materials with defects,which usually depends on empirical parameters.There is no systematic theoretical method to predict the stiffness and shear distribution of glulam beams in elastic-plastic stage,and consequently,the failure of such glulam beams cannot be predicted effectively.To address these issues,an analytical method considering material nonlinearity was proposed for glulam beams,and the calculating equations of deflection and shear stress distribution for different failure modes were established.The proposed method was verified by experiments and numerical models under the corresponding conditions.Results showed that the theoretical calculations were in good agreement with experimental and numerical results,indicating that the equations proposed in this paper were reliable and accurate for such glulam beams with wood material in the elastic-plastic stage ignoring the influence of mechanic properties in radial and tangential directions of wood.Furthermore,the experimental results reported by the previous studies indicated that the method was applicable and could be used as a theoretical reference for predicting the failure of glulam beams.

GENERAL FORM OF MATCHING EQUATION OF ELASTIC-PLASTIC FIELD NEAR CRACK LINE FOR MODE I 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 I crack in elastic-perfectly plastic solids under plane stress condition, and in turn given out the solving process and result for a specific problem.

Interaction of underwater explosion bubble with complex elastic-plastic structure

Considering the elastic-plasticity of the structure, the combination of boundary element method and finite element method （FEM） is employed to present the calculation method for solving the complex coupled dynamic problem of bubble, elastic-plastic structure and the free surface, and the complete three-dimensional calculation program is developed as well. The error between the calculated result and the experimental result is within 10%. Taking a surface ship for example, the three-dimensional calculation program is extended to engineering filed. By employing the program, the response of the ship under the bubble loading is analyzed. From the stress-time history curves of typical elements of the structure, it can be seen that the pressure reaches its maximum when the bubble collapses and this validates that the pressure generated by the bubble collapse and the jet can cause serious damage on the ship structure. From the dynamic process of the interaction between the three-dimensional bubble and the ship, the low order vertical mode of the ship is provoked and the ship presents whip-shaped motion. And the ship does elevation and subsidence movement with the expansion and shrinkage of the bubble. Some rules and conclusions which can be applied to the engineering problems are obtained from the analysis in this paper.

ELASTIC-PLASTIC FINITE ELEMENT ANALYSIS OF THREE-DIMENSIONAL CONTACT-IMPACT AT RAIL JOINT

Using the finite element code ANSYS/LS-DYNA, a dynamic finite element modelwith an elastic-linear-kinematic-hardening plastic material is established to analyzeelastic-plastic stresses in the railhead in the impact process of wheel and rail occurring at thegap of rail joint. The model is based on the discrete elastic support condition of the rails, whichis suitable for the actual situation of wheel/track rolling contact. In the analysis the influencesof axle load, yield stress and tangent modulus of rail material on the stresses and strains areinvestigated in detail. The distribution of stresses and strains in the jointed railhead are given.It is found that the axle load, yield stress and tangent modulus of rail material greatly affect thestresses and strains in the railhead during impacting. The study provides a reliable method anduseful datum for the further research on fatigue and wear of railhead and improving the rail jointmode.

A NOVEL METHOD FOR MEASURING AND CHARACTERIZING DYNAMIC FRACTURE-INITIATION TOUGHNESS OF ELASTIC-PLASTIC MATERIALS

The characterization and testing methods of the dynamic fractureinitiation toughness of elas- tic-plastic materials under tensileimpact are studied. By using the self-designed bar-bar tensile impactappa- ratus, a novel test method for studying dynamicfracture-initiation ahs been proposed based on the one-di- mensionaltest principle. The curve of average load v. s. displacement (P-δ)is smooth until unstable crack propagation, and the kinetic energywhich does not contribute to the crack growth can be removed fromtotal work done by external-force to the specimen.

Elastic-plastic contact force history and response characteristics of circular plate subjected to impact by a projectile

A new elastic-plastic impact-contact model is proposed in this paper. By adopting the principle of minimum acceleration for elastic-plastic continue at finite deformation, and with the aid of finite difference method, the proposed model is applied in the problem of dynamic response of a clamped thin circular plate subjected to a projectile impact centrally. The impact force history and response characteristics of the target plate is studied in detail. The theoretical predictions of the impact force and plate deflection are in good agreements with those of LDA experimental data. Linear expressions of the maximum impact force/transverse deflection versus impact velocity are given on the basis of the theoretical results.

CRACK INITIATION AND STEADY PROPAGATION OF ELASTIC-PLASTIC MATERIALS IN PLANE STRESS CONDITION

The crack tip stress-strain fields of the elastic-plastic cracked specimens have been analyzed using finite element calculations.The crack initiation and steady propagation behaviours have also been investigated by means of slip line pattern etching technique and mechanical tests. The results show that there are HRR near field and distant field in the crack tip region,and the later depends on the specimen configuration.The crack initiation behaviour is controlled by a single parameter J.In contrast,the steady crack propagation is affected by the distant strain field and can not be described by single parameter only.

Harten-Lax-van Leer-contact(HLLC)approximation Riemann solver with elastic waves for one-dimensional elastic-plastic problems

A Harten-Lax-van Leer-contact （HLLC） approximate Riemann solver is built with elastic waves （HLLCE） for one-dimensional elastic-plastic flows with a hypo- elastic constitutive model and the von Mises＇ yielding criterion. Based on the HLLCE, a third-order cell-centered Lagrangian scheme is built for one-dimensional elastic-plastic problems. A number of numerical experiments are carried out. The numerical results show that the proposed third-order scheme achieves the desired order of accuracy. The third-order scheme is used to the numerical solution of the problems with elastic shock waves and elastic rarefaction waves. The numerical results are compared with a reference solution and the results obtained by other authors. The comparison shows that the pre- sented high-order scheme is convergent, stable, and essentially non-oscillatory. Moreover, the HLLCE is more efficient than the two-rarefaction Riemann solver with elastic waves （TRRSE）

An elastic-plastic analysis of short-leg shear wall structures during earthquakes

Short-leg shear wall structures are a new form of building structure that combine the merits of both frame and shear wall structures. Its architectural features, structure bearing and engineering cost are reasonable. To analyze the elastic-plastic response of a short-leg shear wall structure during an earthquake, this study modified the multiple-vertical-rod element model of the shear wall, considered the shear lag effect and proposed a multiple-vertical-rod element coupling beam model with a new local stiffness domain. Based on the principle of minimum potential energy and the variational principle, the stiffness matrixes of a short-leg shear wall and a coupling beam are derived in this study. Furthermore, the bending shear correlation for the analysis of different parameters to describe the structure, such as the beam height to span ratio, short-leg shear wall height to thickness ratio, and steel ratio are introduced. The results show that the height to span ratio directly affects the structural integrity; and the short-leg shear wall height to thickness ratio should be limited to a range of approximately 6.0 to 7.0. The design of short-leg shear walls should be in accordance with the ＂strong wall and weak beam＂ principle.

Numerical Simulation and Experimental Studies on Elastic-Plastic Fatigue Crack Growth

A elastic-plastic fatigue crack growth(FCG)finite element model was developed for predicting crack growth rate under cyclic load.The propagation criterion for this model was established based on plastically dissipated energy.The crack growth simulation under cyclic computation was implemented through the ABAQUS scripting interface.The predictions of this model are in good agreement with the results of crack propagation experiment of compact tension specimen made of 304 stainless steel.Based on the proposed model,the single peak overload retardation effect of elastic-plastic fatigue crack was analyzed.The results shows that the single peak overload will reduce the accumulation rate of plastic energy dissipation of elements at crack tip plastic zone,so that crack growth will be arrested.The crack growth rate will not recover until the crack tip exceed the affected region.Meanwhile,the crack growth rate is mainly determined by the amplitude rather than the mean load under the condition of small scale yielding.The proposed model would be helpful for predicting the growth rate of mode I elastic-plastic fatigue crack.

Seismic elastic-plastic time history analysis and reliability study of quayside container crane

Quayside container crane is a kind of huge dimension steel structure,which is the major equipment used for handling container at modern ports.With the aim to validate the safety and reliability of the crane under seismic loads,besides conventional analysis,elastic-plastic time history analysis under rare seismic intensity is carried out.An ideal finite element（FEM） elastic-plastic mechanical model of the quayside container crane is presented by using ANSYS codes.Furthermore,according to elastic-plastic time history analysis theory,deformation,stress and damage pattern of the structure under rare seismic intensity are investigated.Based on the above analysis,the established reliability model according to the reliability theory,together with seismic reliability analysis based on Monte-Carlo simulation is applied to practical analysis.The results show that the overall structure of the quayside container crane is generally unstable under rare seismic intensity,and the structure needs to be reinforced.

A NEW APPROACH OF PREDICTING BIFURCATION POINTS OF ELASTIC-PLASTIC BUCKLING OF PLATES AND SHELLS

The paper proposes a new approach of predicting the bifurcation points of elastic-plastic buckling of plates and shells,which is obtained from the natural combination of the Lyaponov's dy- namic criterion on stability and the modified adaptive Dynamic Relaxation(maDR)method developed recently by the authors.This new method can overcome the difficulties in the applications of the dy- namic criterion.Numerical results show that the theoretically predicted bifurcation points are in very good agreement with the corresponding experimental ones.The paper also provides a new means for further research on the plastic buckling paradox of plates and shells.

Elastic-plastic behavior of a semicircular frame being pressed against a rigid plane

As a simplified structural model, a semicircular frame is used to study the crashworthiness behavior of an aircraft fuselage. The quasi-static large elastic-plastic deformation of a semicircular frame in the process of its being pressed against a rigid ground is analyzed. First, based on the linear elastic assumption, the quasi-static large deformation contact process of the frame can be divided into three phases, i.e., point contact, line contact and post-buckling. By means of a shooting method, the relations between the displacement and contact force as well as the distribution of bending moment in the three phases are obtained. Then, by assuming an elastic, perfectly-plastic moment-curvature relationship for the semi-circular frame, the contact process is analyzed in detail to reveal the plastic collapse mechanism, the traveling of plastic hinge and the force-displacement relationship. In order to verify the analysis, a preliminary experiment was conducted, in which two types of half rings with clamped ends were pressed by a rigid plate. In addition, a numerical simulation is also conducted by employing ABAQUS to analyze both rectangular cross-sectional beam and I-beam. Finally, the theoretical predictions are compared with the experimental results and numerical solutions, showing that the elastic-plastic analysis can predict the contact process very well.

A new subregion mesh method for the investigation of the elastic-plastic impact in flexible multibody systems

Impact processes between flexible bodies often lead to local stress concentration and wave propagation of high frequency. Therefore, the modeling of flexible multibody systems involving impact should consider the local plastic deformation and the strict requirements of the spatial discretization. Owing to the nonlinearity of the stiffness matrix, the reduction of the element number is extremely important. For the contact-impact problem, since different regions have different requirements regarding the element size, a new subregion mesh method is proposed to reduce the number of the unnecessary elements. A dynamic model for flexible multibody systems with elastic-plastic contact impact is established based on a floating frame of reference formulation and complete Lagrange incremental nonlinear finite-element method to investigate the effect of the elastic-plastic deformation as well as spatial discretization. Experiments on the impact between two bodies are carried out to validate the correctness of the elastic-plastic model. The proposed formulation is applied to a slider-crank system with elastic-plastic impact.

The acoustoelastic effect on Rayleigh waves in elastic-plastic deformed layered rocks

On the basis of the acoustoelastic theory for elastic-plastic materials, the influence of statically deformed states including both the elastic and plastic deformations induced by applied uniaxial stresses on the Rayleigh wave in layered rocks is investigated by using a transfer matrix method. The acoustoelastic effects of elastic plastic strains in rocks caused by static deformations, are discussed in detail. The Rayleigh-type and Sezawa modes exhibit similar trends in acoustoelastic effect： the acoustoelastic effect increasing rapidly with the frequency-thickness product and the phase velocity change approaching a constant value for thick layer and high frequency limit. Elastic-plastic deformations in the Castlegate layered rock obviously modify the phase velocity of the Rayleigh wave and the cutoff points for the Sezawa modes. The investigation may be useful for seismic exploration, geotechnical engineering and ultrasonic detection.

ELASTIC-PLASTIC CONSTITUTIVE RELATION OF PARTICLE REINFORCED COMPOSITES

In this paper, a systematic approach is proposed to obtain the macroscopic elastic-plastic constitutive relation of particle reinforced composites (PRC). The strain energy density of PRC is analyzed based on the cell model, and the analytical formula for the macro-constitutive relation of PRC is obtained. The strength effects of volume fraction of the particle and the strain hardening exponent of matrix material on the macro-constitutive relation are investigated, the relation curve of strain versus stress of PRC is calculated in detail. The present results are consistent with the results given in the existing references.

On the yield functions in Lagrangian stress space of elastic-plastic materials with internal constraints

In the present paper it is shown that the elastic range in the second Piola-Kirchhoff stress space can be chosen in a hyperplane which is through the origin of Lagrangian stress space and perpendicular to the normal of the constraint manifold at the plastic configuration, if the determinate stress response of the elastic-plastic material with simple internal constraints with some condition is correctly chosen, otherwise, it is in general in a hypersurface and the normal flow rule by Ⅱ yushin＇s postulate will have an indeterminate part. The choice of deterrninate stress response is probable because of its indeterminacy. Therefore the yield function should be a function of the second Piola-Kirchhoff stress lying in the hyperplane so that it is more simple and the back stress as the geometric center of the elastic range in general is inside the elastic range. Finally some examples are concerned.

PATH-INDEPENDENT J-INTEGRAL AND ITS DUAL FORM IN ELASTIC-PLASTIC SOLIDS WITH FINITE DISPLACEMENTS

In this paper, based on energy variational principles of elastic-plastic solids, the path-independent J-integral and its dual form in elastic-plastic solids with finite displacements are presented. Whose testification is given there after.

Analysis of the elastic-plastic indentation properties of materials with varying ratio of hardness to Young’s modulus

The elastic-plastic indentation properties of materials with varying ratio of hardness to Young’s modulus(H/E) were analyzed with the finite element method. And the indentation stress fields of materials with varying ratio H/E on the surface were studied by the experiment. The results show that the penetration depth, contact radius, plastic pile-up and the degree of elastic recovery depend strongly on the ratio H/E. Moreover, graphs were established to describe the relationship between the elastic-plastic indentation parameters and H/E. The established graphs can be used to predict the H/E of materials when compared with experimental data.

MIXED ENERGY METHOD FOR SOLUTION OF QUADRATIC PROGRAMMING PROBLEMS AND ELASTIC-PLASTIC ANALYSIS OF TRUSS STRUCTURES

A new algorithm for the solution of quadratic programming problemsis put forward in terms of the mixed energy theory and is furtherused for the incremental solution of elastic-plastic trussstructures. The method proposed is different from the traditionalone, for which the unknown variables are selected just in one classsuch as displacements or stresses. The present method selects thevariables in the mixed form with both displacement and stress. As themethod is established in the hybrid space, the information found inthe previous incremental step can be used for the solution of thepresent step, making the algorithm highly effi- cient in thenumerical solution process of quadratic programming problems. Theresults obtained in the exm- ples of the elastic-plastic solution ofthe truss structures verify what has been predicted in thetheoretical anal- ysis.