A new method for evaluation of stress intensities for interface cracks

A new method is presented for calculating the values of K-I and K-II in the elasticity solution at the tip of an interface crack. The method is based on an evaluation of the J-integral by the virtual crack extension method. Expressions for calculating K-I and K-II by using the displacements and the stiffness derivative of the finite element solution and asymptotic crack tip displacements are derived. The method is shown to produce very accurate solutions even with coarse element mesh.

Variable power singular interface elements for a crack normal to the bimaterial interface

Zero thickness crack tip interface elements for a crack normal to the interface between two materials are presented. The elements are shown to have the desired r(lambda-1) (0 < lambda < 1) singularity in the stress field at the crack tip and are compatible with other singular elements. The stiffness matrices of the quadratic and cubic interface element are derived. Numerical examples are given to demonstrate the applicability of the proposed interface elements for a crack perpendicular to the bimaterial interface.

Micromechanics of elastic solids weakened by a large number of cracks

This paper presents a micromechanics analysis of the elastic solids weakened by a large number of microcracks in a plane problem. A new cell model is proposed. Each cell is an ellipse subregion and contains a microcrack. The effective moduli and the stress intensity factors for an ellipse cell are obtained. The analytic closed formulas of concentration factor tensor for an isotropic matrix containing an anisotropic inclusion are derived. Based on a self-consistent method, the effective elastic moduli of the solids weakened by randomly oriented microcracks are obtained.

Local stresses and strains in axially cracked cylindrical shells

A perturbation solution is obtained for the local stress-strain fields in an axially cracked cylindrical shell. The tenth-order differential equations are used that take into account the transverse shear deformation. The perturbation of a curvature parameter, λ, is employed, where . The stress intensity factors for finite size cylindrical shells subjected to bending and internal pressure are evaluated. Sufficient accuracy can be obtained without using fine mesh sizes in regions near the crack tip. Also analyzed are the influence of cylinder diameter and shearing stiffness on bulging.

Structures on sliding base subject to horizontal and vertical motions

In this paper particular investigation is directed towards the combined effects of horizontal and vertical motions of real earthquakes to structures resting on sliding base. A simplified method is presented to treat the nonlinear effects of time dependent frictional force of the sliding base as a function of the vertical reaction produced by the foundation. As an example, the El Centro 1940 earthquake record is used on a structural model to show the structural responses due to a sliding base with different frictional and stiffness characteristics. The study shows that vertical ground motion does affect both the superstructure response and the base sliding displacement. Nevertheless, the sliding base isolator is shown to be effective for the reduction of seismic response of a superstructure.

Study on cracked spherical shells

The local-global anatysis method is systematically extended to the fracture analysis of spherical shells. On the basis of the shallow shell theory, which takes into account transverse shear deformations, governing equations for cracked spherical shells expressed in displacement and stress functions f, F and φ are proposed, and then a general solution including Modes, Ⅰ, Ⅱ, Ⅲ for stress-strain fields at crack tip in a spherical shell is obtained, which plays the same role as Williams's expansion in plane elasticity. The numerical results for finite-size spherical shells under different boundary conditions have been obtained. Furthermore, the bulging factors are analyzed with regard to shearing stiffness and an approximate formula is given.

Dynamic-response of a laminated plate with friction damping

考虑刚度及边界条件的索力精确求解

目前索力的测试大多不考虑结构刚度的影响，这导致测试的结果偏大。本文中不仅考虑了索的刚度，而且将索的边界条件假定为简支以及简支和固支耦合两种情况，得到不同假设下的理论公式，分析了不同方法下索力的区别，并利用实测的数据进行分析，给出了末考虑刚度情况下和仅考虑为简支时的误差，证明在索力测试中，特别是短索的测试中，特别是短索的测试中刚度以及边界条件考虑的必要性。

Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment

Study of the damage-induced anisotropy of quasi-brittle materials using the component assembling model

Damage-induced anisotropy of quasi-brittle materials is investigated using component assembling model in this study. Damage-induced anisotropy is one significant character of quasi-brittle materials coupled with nonlinearity and strain softening. Formulation of such complicated phenomena is a difficult problem till now. The present model is based on the component assembling concept, where constitutive equations of materials are formed by means of assembling two kinds of components' response functions. These two kinds of components, orientational and volumetric ones, are abstracted based on pair-functional potentials and the Cauchy - Born rule. Moreover, macroscopic damage of quasi-brittle materials can be reflected by stiffness changing of orientational components, which represent grouped atomic bonds along discrete directions. Simultaneously, anisotropic characters are captured by the naturally directional property of the orientational component. Initial damage surface in the axial-shear stress space is calculated and analyzed. Furthermore, the anisotropic quasi-brittle damage behaviors of concrete under uniaxial, proportional, and nonproportional combined loading are analyzed to elucidate the utility and limitations of the present damage model. The numerical results show good agreement with the experimental data and predicted results of the classical anisotropic damage models.

Effects of fundamental structure parameters on dynamic responses of submerged floating tunnel under hydrodynamic loads

This paper investigates the effects of structure parameters on dynamic responses of submerged floating tunnel (SFT) under hydrodynamic loads. The structure parameters includes buoyancy-weight ratio (BWR), stiffness coefficients of the cable systems, tunnel net buoyancy and tunnel length. First, the importance of structural damp in relation to the dynamic responses of SFT is demonstrated and the mechanism of structural damp effect is discussed. Thereafter, the fundamental structure parameters are investigated through the analysis of SFT dynamic responses under hydrodynamic loads. The results indicate that the BWR of SFT is a key structure parameter. When BWR is 1.2, there is a remarkable trend change in the vertical dynamic response of SFT under hydrodynamic loads. The results also indicate that the ratio of the tunnel net buoyancy to the cable stiffness coefficient is not a characteristic factor affecting the dynamic responses of SFT under hydrodynamic loads.

Prediction of Effective Moduli of Carbon Nanotube-Reinforced Composites with Waviness and Debonding

Carbon nanotubes (CNTs) have been regarded as ideal reinforcements of high-performance composites with enormous applications. However, the waviness of the CNTs and the interfacial bonding condition between them and the matrix are two key factors that influence the reinforcing efficiency. In this paper, the effects of the waviness of the CNTs and the interfacial debonding between them and the matrix on the effective moduli of CNT-reinforced composites are studied. A simple analytical model is presented to investigate the influence of the waviness on the effective moduli. Then, two methods are proposed to examine the influence of the debonding. It is shown that both the waviness and debonding can significantly reduce the stiffening effect of the CNTs. The effective moduli are very sensitive to the waviness when the latter is small, and this sensitivity decreases with the increase of the waviness. (C) 2008 Elsevier Ltd. All rights reserved.

纤维复合材料刚度设计的力学原理及其应用

针对航空用纤维复合材料,基于纤维与基体的协调变形规律和材料应变能密度的计算,导出了纤维复合材料的刚度表达式,研究了刚度矩阵的基本力学特性,给出了复合材料刚度矩阵元素与基体材料的刚度、纤维束的刚度及铺设方向的关系,提供了平面铺层和三维取向纤维复合材料的基本刚度设计公式和设计方法.最后给出工程应用设计方法的几个具体例子,并与其他理论和实验结果进行了对比.研究结果表明,在考虑材料制备工艺缺陷造成的刚度折减后,本文提出的设计公式及方法可应用于纤维复合材料的优化设计.

Numerical investigation of crack growth in concrete subjected to compression by the generalized beam lattice model

The beam lattice-type models, such as the Euler-Bernoulli (or Timoshenko) beam lattice and the generalized beam (GB) lattice, have been proved very effective in simulating failure processes in concrete and rock due to its simplicity and easy implementation. However, these existing lattice models only take into account tensile failures, so it may be not applicable to simulation of failure behaviors under compressive states. The main aim in this paper is to incorporate Mohr-Coulomb failure criterion, which is widely used in many kinds of materials, into the GB lattice procedure. The improved GB lattice procedure has the capability of modeling both element failures and contact/separation of cracked elements. The numerical examples show its effectiveness in simulating compressive failures. Furthermore, the influences of lateral confinement, friction angle, stiffness of loading platen, inclusion of aggregates on failure processes are respectively analyzed in detail.

The elasto-damage theory of the components assembling model

The potential energy in materials is well approximated by pair functional which is composed of pair potentials and embedding energy. During calculating material potential energy, the orientational component and the volumetric component are derived respectively from pair potentials and embedding energy. The sum of energy of all these two kinds of components is the material potential. No matter how microstructures change, damage or fracture, at the most level, they are all the changing and breaking atomic bonds. As an abstract of atomic bonds, these components change their stiffness during damaging. Material constitutive equations have been formulated by means of assembling all components' response functions. This material model is called the component assembling model. Theoretical analysis and numerical computing indicate that the proposed model has the capacity of reproducing some results satisfactorily, with the advantages of great conceptual simplicity, physical explicitness, and intrinsic induced anisotropy, etc.