EFFECT OF Z-PINS’ DIAMETER,SPACING AND OVERLAP LENGTH ON CONNECTING PERFORMANCE OF CMC SINGLE LAP JOINT

The effect of through-thickness reinforcement by composite pins (Z-pins) on the static tensile strength and failure mechanisms of the joints made from ceramic matrix composite (CMC) is investigated. Overlap length of the single lap joint is 15 mm, 20 mm, 23 mm, 37 mm, and 60 mm, respectively. The experimental results indicate that the final failure modes of the joints can be divided into two groups, (a) the bond-line stops debonding until crack encounters Z-pins; and then the adherends break at the location of Z-pins, when overlap length is more than 20 mm; (b) the bond-line detaches entirely and Z-pins are drawn from adherends, when overlap length is equal to 15 mm. A simple effcient computational approach is presented for analyzing the benefit of through-thickness pins for restricting failure in the single lap joints. Here, the mechanics problem is simplified by representing the effect of the pins by tractions acting on the fracture surfaces of the cracked bond-line. The tractions are prescribed as functions of the crack displacement, which are available in simple forms that summarize the complex deformations to a reasonable accuracy. The resulting model can be used to track the evolution of complete failure mechanisms, for example, bond-line initial delamination and ultimate failure associated with Z-pin pullout, ultimate failure of the adherends. The paper simulates connecting performance of the single lap joints with different Z-pins’ diameter, spacing and overlap length; the numerical results agree with the experimental results; the numerical results indicate enlarging diameter and decreasing spacing of Z-pins are in favor of improving the connecting performance of the joints. By numerical analysis method, the critical overlap length that lies between two final failure modes is between 18 mm and 19 mm, when Z-pins’ diameter and spacing are 0.4 mm, 5 mm, respectively.     

THE CONSTITUTIVE RELATION OF CRACK-WEAKENED ROCK MASSES UNDER AXIAL-DIMENSIONAL UNLOADING

An accurate and efficient numerical method for solving the crack-crack interaction problem is presented. The method is mainly by means of the dislocation model, stress superposition principle and Chebyshev polynomial expansion of the pseudo-traction. This method can be applied to compute the stress intensity factors of multiple kinked cracks and multiple rows of periodic cracks as well as the overall strains of rock masses containing multiple kinked cracks under complex loads. Many complex computational examples are given. The dependence of the crack-crack interaction on the crack configuration, the geometrical and physical parameters, and loads pattern, is investigated. By comparison with numerical results under confining pressure unloading, it is shown that the crack-crack interaction under axial-dimensional unloading is weaker than those under confining pressure unloading. Numerical results for single faults and crossed faults show that the single faults are more unstable than the crossed faults. It is found from numerical results for different crack lengths and different crack spacing that the interaction among kinked cracks decreases with an increase in length of the kinked cracks and the crack spacing under axial-dimensional unloading.     

Convective heat transfer from two rotating circular cylinders in tandem arrangement by using lattice Boltzmann method

The numerical investigation of the two-dimensional laminar flow past two ro- tating circular cylinders in the tandem arrangement is conducted by the lattice Boltzmann method. The numerical strategy is used for dealing with curved and moving boundaries of the second-order accuracy for velocity and temperature fields. The effects of various rotational speed ratios and gap spacing are studied with the Reynolds number of 100 and the Prandtl number of 0.71. A varied range of rotational speed ratios are investigated for four different gap spacing, i.e., 3.0, 1.5, 0.7, and 0.2. The results show that, for the first cylinder, the lift and drag coefficients for large gap spacing are similar to those for a single cylinder; for the second cylinder, the lift coefficient descends with the increase in the angular velocity for all gap spacing, while the drag coefficient ascends except for the gap spacing of 3.0. The results of the averaged periodic Nusselt number on the surface of the cylinders show that, for small distances between the cylinders and low angular velocities, conduction is a dominant mechanism of heat transfer, but for large distances and high angular velocities, convection is the main mechanism of heat transfer.     

ANALYTICAL SOLUTIONS FOR RESPONSE OF COLLISION OF PARTICLE WITH CONICAL ROD CAUSED BY LONGITUDINAL VIBRATION

The objective is to present exact analytical solutions of longitudinal impact analysis for slender conical rods struck by a particle and a new method is proposed for conical rod-particle impact analysis, in which the superposition method is used and the response of the rod is presented. These analytical results are exact and can be used to validate the numerical methods or other analytical results. The numerical example shows that one of the advantages of the present method is that the analytical form is very simple. The result is that mass ratio and some variables describing the geometrical shape of rods such as taper, length and radius play an important role in impact dynamic system.     

STUDY ON TENSILE FAILURE MECHANISM AND HYBRID EFFECT OF INTRAPLY HYBRID COMPOSITES

The random critical-core model is adapted to investigate the tensile failure mech-anism and hybrid effect of unidirectionally arrayed hybrid composites with alternating low andhigh elongation fibers.By utilizing the model in conjunction with the results of the stress con-centration analysis in which the interfacial damage between fiber and matrix is considered,amicroscopic statistical analysis of both the first failure and ultimate failure of hybrids is per-formed.The variations of the first failure strain,the ultimate failure strain and the hybrid effectas the interfacial shear strength are obtained quantitatively.The concept of the hybrid effect forstrains has been clarified.The present results are compared with available experiment data anda reasonable agreement is found between the analytical predictions and the experimental results.     

NUMERICAL INVESTIGATION OF THE LIMIT LOADS FOR PRESSURE VESSELS WITH PART-THROUGH SLOTS

A numerical investigation of the limit loads is carried out for pressurevessels with part-through slots using a general computational method for the limitanalysis of 3-D structures.The limit pressures are given for a comprehensive range ofgeometric parameters.Some of the calculated results are compared with the results of3-D elastic-plastic finite element analysis and existing numerical solutions.The effectsof various shapes and sizes of part-through slots on the load carrying capacity ofcylindrical shells are investigated and evaluated.Two kinds of typical failure modescorresponding to different dimensions of slots are studied.Based on the numericalresults,a geometric parameter G which combines the slot dimensions and the cylindergeometry is presented.It reasonably reflects the overall effect of slots on the limit loadsof cylinders.An empirical formula for estimating the limit pressures of cylindricalshells with part-through slots is obtained.     

ANALYSIS OF ADHESIVE LAP JOINT

This paper is to solve the interlaminar stresses of adhesive lap joint by the energymethod without considering the adhesive layer. The joint is made of two identical narrowplates.Two cases are discussed:one is for the isotropic material and the other is fororthotropic material.Because of the different materials forming the joint,the length ofdistribution and the magnitude of the interlaminar stresses for the two cases will be verydifferent.     

THE APPLICATION OF DISCRETE ELEMENT METHOD IN SOLVING THREE-DIMENTIONAL IMPACT DYNAMICS PROBLEMS

A three-dimensional discrete element model of the connective type is presented. Moreover, a three- dimensional numerical analysis code, which can carry out the transitional process from connective model (for continuum) to contact model (for non-continuum), is developed for simulating the mechanical process from continuum to non-continuum. The wave propagation process in a concrete block (as continuum) made of cement grout under impact loading is numerically simulated with this code. By comparing its numerical results with those by LS-DYNA, the calculation accuracy of the model and algorithm is proved. Furthermore, the failure process of the concrete block under quasi-static loading is demonstrated, showing the basic dynamic transitional process from continuum to non-continuum. The results of calculation can be displayed by animation. The damage modes are similar to the experimental results. The two numerical examples above prove that our model and its code are powerful and efficient in simulating the dynamic failure problems accompanying the transition from continuum to non-continuum. It also shows that the discrete element method (DEM) will have broad prospects for development and application.     

Pendulum with linear damping and variable length

The methods of multiple scales and approximate potential are used to study pendulums with linear damping and variable length. According to the order of the coefficient of friction compared with that of the slowly varying parameter of length, three different cases are discussed in details. Asymptotic analytical expressions of amplitude, frequency and solution are obtained. The method of approximate potential makes the results effective for large oscillations. A modified multiple scales method is used to get more accurate leading order approximations when the coefficient friction is not small. Comparisons are also made with numerical results to show the efficiency of the present method.     

MESO-MECHANICAL ANALYSIS OF SHAPE MEMORY ALLOY REINFORCED SMART STRUCTURE WITH DAMAGE

The mechanical behaviors of shape memory alloy (SMA) wires reinforced smart structure with damage were analyzed through the variational principle, a governing equation for the structure was derived, mathematical expressions for the meso-displacement field, stress-strain field of typical element with damage were presented, and a failure criterion for interface failure between SMA wires and matrix was established under two kinds of actuation which are dead-load and temperature, where the temperature is included in effective free restoring strain. In addition, there are some other composing factors in the failure criterion such as the interface properties, dynamical properties of SMA, initial debonding length L-l etc. The results are significant to understand structural strength self-adapted control and failure mechanism of SMA wires reinforced smart structure with damage.     

REFLECTION CHARACTERISTICS OF LONGITUDINAL WAVES IN A SEMI-INFINITE CYLINDRICAL ROD CONNECTED TO AN INFINITE VISCOELASTIC STRATUM

Reflection characteristics of longitudinal strain waves in a semi-infinite elastic rod con-nected to a viscoelastic stratum are investigated analytically.The three-dimensional viscoelasticity the-ory is applied to the stratum,and the Laplace transform with respect to time and the numerical inverseLaplace transform by means of Laguerre function are used.The time histories for the longitudinalstrain of an arbitrary point of the rod are presented.Two typical viscoelastic models are considered,one is the usual Maxwell-Voigt model,the other is whose relaxation function is given by a power law.The numerical results for the two models are presented and compared each other and also with previ-ously published results for the elastic stratum.     

AN ESTIMATION METHOD OF BEARING STRENGTH OF BOLTED JOINTS IN FIBRE REINFORCED COMPOSITE

Some researchers have estimated the strength of bolted joints in fibre reinforcedcomposite,using simple and efficient engineering procedures.However,for theseprocedures the effect of clamping due to the strength of bolted joints is not considered.Inthis paper,a method is presented for predicating critical bearing strength of single-holebolted joints in composite on the basis of observing and analysing the results of experiments.The clamping effect of bolts is considered.The calculated results correspond to the test dataon Glaphic/Epoxy laminates.     

DAMPING OF VERTICALLY EXCITED SURFACE WAVE IN WEAKLY VISCOUS FLUID

In a vertically oscillating circular cylindrical container, singular perturbation theory of two-time scale expansions is developed in weakly viscous fluids to investigate the motion of single free surface standing wave by linearizing the Navier-Stokes equation. The fluid field is divided into an outer potential flow region and an inner boundary layer region. The solutions of both two regions are obtained and a linear amplitude equation incorporating damping term and external excitation is derived. The condition to appear stable surface wave is obtained and the critical curve is determined. In addition, an analytical expression of damping coefficient is determined. Finally, the dispersion relation, which has been derived from the inviscid fluid approximation, is modified by adding linear damping. It is found that the modified results are reasonably closer to experimental results than former theory. Result shows that when forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when forcing frequency is high, the surface tension of the fluid is prominent.     

Validation of a CIP-based tank for numerical simulation of free surface flows

A constrained interpolation profile CIP-based numerical tank is developed to simulate violent free surface flows.The numerical simulation is performed by the CIP-based Cartesian grid method,which is described in the present paper.The tangent of hyperbola for interface capturing(THINC) scheme is applied for capturing complex free surfaces.The new model is capable of simulating a flow with violently varied free surface.A series of computations are conducted to assess the developed algorithm and its versatility.These tests include the collapse of water column with and without an obstacle,sloshing in a fixed tank,the generation of regular waves in a tank,the generation of extreme waves in a tank.Excellent agreements are obtained when numerical results are compared with available analytical,experimental,and other numerical results.     

Effects of size-dependent elasticity on stability of nanotweezers

It is well-recognized that the electromechanical response of a nanostructure is affected by its element size. In the present article, the size dependent stability behavior and nanotweezers fabricated from nanowires are investigated by modified couple stress elasticity （MCSE）. The governing equation of the nanotweezers is obtained by taking into account the presence of Coulomb and intermolecular attractions. To solve the equation, four techniques, i.e., the modified variational iteration method （MVIM）, the monotonic iteration method （MIM）, the MAPLE numerical solver, and a lumped model, are used. The variations of the arm displacement of the tweezers versus direct current （DC） voltage are obtained. The instability parameters, i.e., pull-in voltage and deflection of the system, are computed. The results show that size-dependency will affect the stability of the nanotweezers significantly if the diameter of the nanowire is of the order of the length scale. The impact of intermolecular attraction on the size-dependent stability of the system is discussed.     

ON THE NUMERICAL SOLUTION OF QUASILINEAR WAVE EQUATION WITH STRONG DISSIPATIVE TERM

The numerical solution for a type of quasilinear wave equation is studied. The three-level difference scheme for quasi-linear waver equation with strong dissipative term is constructed and the convergence is proved. The error of the difference solution is estimated. The theoretical results are controlled on a numerical example.     

Dynamic stress concentrations in thick plates with two holes based on refined theory

Based on complex variables and conformal mapping, the elastic wave scat- tering and dynamic stress concentrations in the plates with two holes are studied by the refined dynamic equation of plate bending. The problem to be solved is changed to a set of infinite algebraic equations by an orthogonM function expansion method. As examples, under free boundary conditions, the numerical results of the dynamic moment concen- tration factors in the plates with two circular holes are computed. The results indicate that the parameters such as the incident wave number, the thickness of plates, and the spacing between holes have great effects on the dynamic stress distributions. The results are accurate because the refined equation is derived without any engineering hypothese.     

Theoretical statistical solution and numerical simulation of heterogeneous brittle materials

The analytical stress-strain relation with heterogeneous parameters is derived for theheterogeneous brittle materials under a uniaxial extensional load, in which the distributions of theelastic modulus and the failure strength are assumed to be statistically independent. This theoreticalsolution gives an approximate estimate of the equivalent stress-strain relations for 3-D heterogeneousmaterials. In one-dimensional cases it may provide comparatively accurate results. The theoreticalsolution can help us to explain how the heterogeneity influences the mechanical behaviors, Further, anumerical approach is developed to model the non-linear behavior of three-dimensional heterogeneousbrittle materials. The lattice approach and statistical techniques are applied to simulate the initialheterogeneity of heterogeneous materials. The load increment in each loading stage is adaptivelydetermined so that the better approximation of the failure process can be realized. When the maximumtensile principal strain exceeds the failure strain, the elements are considered to be broken, which canbe carried out by replacing its Young‘s modulus with a very small value. A 3-D heterogeneous brittlematerial specimen is simulated during a full failure process. The numerical results are in good agreementwith the analytical solutions and experimental data.     

UNSTEADY／STEADY NUMERICAL SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON ARTIFICIAL COMPRESSIBILITY

A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.     

EQUIVALENT MODEL OF EXPANSION OF CEMENT MORTAR UNDER SULPHATE EROSION

The expansion property of cement mortar under the attack of sulfate ions is studied by experimental and theoretical methods. First, cement mortars are fabricated with the ratio of water to cement of 0.4, 0.6, and 0.8. Secondly, the expansion of specimen immerged in sulphate solution is measured at different times. Thirdly, a theoretical model of expansion of cement mortar under sulphate erosion is suggested by virtue of represent volume element method. In this model, the damage evolution due to the interaction between delayed ettringite and cement mortar is taken into account. Finally, the numerical calculation is results indicate that the model perfectly describes performed. The numerical and experimental the expansion of the cement mortar.