GURTIN-TYPE VARIATIONAL PRINCIPLES FOR DYNAMICS OF A NON-LOCAL THERMAL EQUILIBRIUM SATURATED POROUS MEDIUM

Based on the porous media theory and by taking into account the efects of the pore fuid viscidity, energy exchanges due to the additional thermal conduction and convection between solid and fuid phases, a mathematical model for the dynamic-thermo-hydro-mechanical coupling of a non-local thermal equilibrium fuid-saturated porous medium, in which the two constituents are assumed to be incompressible and immiscible, is established under the assumption of small de- formation of the solid phase, small velocity of the fuid phase and small temperature changes of the two constituents. The mathematical model of a local thermal equilibrium fuid-saturated porous medium can be obtained directly from the above one. Several Gurtin-type variational principles, especially Hu-Washizu type variational principles, for the initial boundary value problems of dy- namic and quasi-static responses are presented. It should be pointed out that these variational principles can be degenerated easily into the case of isothermal incompressible fuid-saturated elastic porous media, which have been discussed previously.

A SIMPLIFIED CALCULATING METHOD OF NONLINEAR FREQUENCY OF CABLE NET UNDER MEAN WIND LOAD

The cable net supported glass curtain wallas the most advanced technique in dot point supported glass curtain wall, is widely used in China. Because of its large deflection and high nonlinearity under wind load, the dynamic performance of the cable net is greatly different from that of the conventional linear structures. The continuous membrane theory is used to construct the nonlinear vibration differential equation of the cable net, and the harmonic balance method is used to solve the analytic formula of the nonlinear frequency. In order to verify the accuracy of the above analytic formula, the results of the formula and the nonlinear FEM time-history method are compared and found to be in good agreement. Furthermore, the nonlinear vibration differential equation and the nonlinear frequency obtained in this paper are the basis for the wind-induced response analysis of a cable net under fluctuating wind load.

BENDING OF FUNCTIONALLY GRADED PIEZOELECTRIC RECTANGULAR PLATES

By introducing two displacement functions as well as two stressfunctions, two independent state equations with variable coefficientsare derived from the three-dimensional theory equations of piezo-elasticity for transverse isotropy. A laminated approximation is usedto transform the state equations to those with constant coefficientsin each sub-layer. The bending problem of a functionally gradedrectangular plate is then analyzed based on the state equations.Numerical results are presented and the effect of material gradi- entindex is discussed.

RELIABILITY ANALYSIS OF FRP LAMINATED PLATES WITH CONSIDERATION OF BOTH INITIAL IMPERFECTION AND FAILURE SEQUENCE

A probabilistic progressive failure analyzing method is applied to estimating the reliability of a simply supported laminated composite plate with an initial imperfection under bi-axial compression load. The initial imperfection and the strength parameters are considered as random variables. Ply-level failure probability is evaluated by the first order reliability method (FORM) together with the Tsai-Wu strength criterion and Tan criterion. Current stresses in the laminated structure are calculated by the classical lamination theory with the stiffness modified based on the last step ply failure. Probabilistically dominant ply-level failure sequences leading to overall system failure are identified, based on which the system failure probability is estimated. A numerical example is presented to demonstrate the methodology proposed. Through parameter studies it is shown that the deviation of the initial imperfection and some of the strength parameters largely influence the system reliability.

ELASTIC MODULUS REDUCTION METHOD FOR LIMIT LOAD EVALUATION OF FRAME STRUCTURES

A new strategy for elastic modulus adjustment is proposed based on the element bearing ratio （EBR）,and the elastic modulus reduction method （EMRM） is proposed for limit load evaluation of frame structures. The EBR is defined employing the generalized yield criterion,and the reference EBR is determined by introducing the extrema and the degree of uniformity of EBR in the structure. The elastic modulus in the element with an EBR greater than the reference one is reduced based on the linear elastic finite element analysis and the equilibrium of strain energy. The lower-bound of limit-loads of frame structures are analyzed and the numerical example demonstrates the flexibility,accuracy and effciency of the proposed method.

ZONAL DISINTEGRATION MECHANISM OF DEEP CRACK-WEAKENED ROCK MASSES UNDER DYNAMIC UNLOADING

Size and quantity of fractured zone and non-fractured zone are controlled by cracks contained in deep rock masses. Zonal disintegration mechanism is strongly dependent on the interaction among cracks. The strong interaction among cracks is investigated using stress superposition principle and the Chebyshev polynomials expansion of the pseudo-traction. It is found from numerical results that crack nucleation, growth and coalescence lead to failure of deep crack- weakened rock masses. The stress redistribution around the surrounding rock mass induced by unloading excavation is studied. The effect of the excavation time on nucleation, growth, interaction and coalescence of cracks was analyzed. Moreover, the influence of the excavation time on the size and quantity of fractured zone and non-fractured zone was given. When the excavation time is short, zonal disintegration phenomenon may occur in deep rock masses. It is shown from numerical results that the size and quantity of fractured zone increase with decreasing excavation time, and the size and quantity of fractured zone increase with the increasing value of in-situ geostress.

WAVELET METHOD FOR CALCULATING THE DEFECT STATES OF TWO-DIMENSIONAL PHONONIC CRYSTALS

Based on the variational theory, a wavelet-based numerical method is developed to calculate the defect states of acoustic waves in two-dimensional phononic crystals with point and line defects. The supercell technique is applied. By expanding the displacement field and the material constants （mass density and elastic stiffness） in periodic wavelets, the explicit formulations of an eigenvalue problem for the plane harmonic bulk waves in such a phononic structure are derived. The point and line defect states in solid-liquid and solid-solid systems are calculated. Comparisons of the present results with those measured experimentally or those from the plane wave expansion method show that the present method can yield accurate results with faster convergence and less computing time.

COHESIVE ZONE FINITE ELEMENT-BASED MODELING OF HYDRAULIC FRACTURES

Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.

AN ANALYTICAL SOLUTION TO LONGITUDINAL VIBRATION OF A PILE OF ARBITRARY SEGMENTS WITH VARIABLE MODULUS

The vibration problem of a pile of arbitrary segments withvariable modulus modules under ex- citing force is established, inwhich the influence of the soil under pile toe and the surroundingsis tak- en into account. With Laplace transforms, the transmitfunctions for velocity and displacement of pile are derived.Furthermore, in terms of the convolution theorem and inversed Laplacetransform, an analytical solution for the time domain response of apile subjected to semi-sine impulse is developed, which is thetheoretical basis of the sonic method in pile integrity testing. Basedon the solution, the vibration properties of pile with sharp orcontinuous modulus are studied. The validity of this approach isverified through field dynamic tests on some engineering piles. Itshows that the theoretical predic- tion ad the response of the pileare in good agreement.

NON-LINEAR FORCED VIBRATION OF AXIALLY MOVING VISCOELASTIC BEAMS

The non-linear forced vibration of axially moving viscoelastic beams excited by the vibration of the supporting foundation is investigated. A non-linear partial-differential equation governing the transverse motion is derived from the dynamical, constitutive equations and geometrical relations. By referring to the quasi-static stretch assumption, the partial-differential non-linearity is reduced to an integro-partial-differential one. The method of multiple scales is directly applied to the governing equations with the two types of non-linearity, respectively. The amplitude of near- and exact-resonant steady state is analyzed by use of the solvability condition of eliminating secular terms. Numerical results are presented to show the contributions of foundation vibration amplitude, viscoelastic damping, and nonlinearity to the response amplitude for the first and the second mode.

DIGITAL SPECKLE CORRELATION METHOD IMPROVED BY GENETIC ALGORITHM

The digital speckle correlation method is an important optical metrology for sur- face displacement and strain measurement.With this technique,the whole field deformation in- formation can be obtained by tracking the geometric points on the speckle images based on a correlation-matching search technique.However,general search techniques suffer from great com- putational complexity in the processing of speckle images with large deformation and the large random errors in the processing of images of bad quality.In this paper,an advanced approach based on genetic algorithms (GA) for correlation-matching search is developed.Benefiting from the abilities of global optimum and parallelism searching of GA,this new approach can complete the correlation-matching search with less computational consumption and at high accuracy.Two experimental results from the simulated speckle images have proved the efficiency of the new approach.