倾斜内锁型三维机织陶瓷基复合材料的细观力学模型

利用平均化方法提出了倾斜内锁型三维机织陶瓷基复合材料弹性性能分析的三维细观力学模型,对材料的弹性性能进行了预测。这个力学模型考虑了倾斜内锁型三维机织陶瓷基复合材料经向纤维束的弯曲和纬向纤维束的平直,纤维束的横截面形状尺寸和相邻纤维束之间的孔洞以及材料制造过程中碳纤维性能下降对弹性性能的影响。基于层合板理论,提出两种单胞应变状态假设分别对材料的九个弹性常数进行了推导计算,结果表明两种方法理论的预测值非常接近。计算结果与实验值比较吻合,表明所提出的细观力学模型是合理的,可以为纺织陶瓷基复合材料的优化设计提供有价值的参考。     

Z-pin复合材料细观模型固化残余应力研究

针对植入Z-pin后碳纤维增强复合材料的微观结构，通过施加Z-pin周期性边界约束条件，建立了Z-pin复合材料单层板单胞细观模型．考虑固化过程中树脂体积收缩、弹性模量随固化度变化和纤维因Z-pin进入偏转因素，运用有限单元法计算了单胞结构在固化成型工艺过程中树脂和纤维应力发展和分布，并研究了Z-pin直径和分布密度对单层板面内残余应力的影响．结果表明：凝胶点之前，树脂模量和残余应力很小，凝胶点之后，树脂模量和残余应力增加较快；残余应力分布与纤维偏转有关；Z-pin直径和分布密度增加会使固化残余应力增大．     

二维平纹编织复合材料压缩力学行为研究

通过对二维平纹编织复合材料压缩力学性能进行三维有限元计算,对复合材料在沿垂直于编织平面方向和纤维束方向压缩后的微观应力场进行了分析。对于沿垂直于编织平面方向的压缩(即横向压缩),发现在经纬向纤维束相互交替处基体材料的横向压应力最大,是其薄弱部位。分析表明,与0°/90°层合板相比,编织复合材料受到横向压缩载荷时,自由边界处层间剪应力较大,容易发生层间裂纹,与实验观察相吻合。通过改变纤维束截面尺寸和纤维束之间间隙大小研究了破坏面角度的变化。     

基于渐进均匀化的平纹编织复合材料低速冲击多尺度方法

针对平纹编织复合材料低速冲击响应和损伤问题,提出了一种多尺度分析方法. 首先, 建立微观尺度单胞模型,引入周期性边界条件,采用最大主应力失效准则和直接刚度退化模型表征纤维丝和基体的损伤起始与演化,预测了纤维束的弹性性能和强度性能. 其次,将这些性能参数代入介观尺度单胞模型,基于Hashin和Hou的混合失效准则以及连续介质损伤模型对介观尺度单胞进行6种边界条件下的渐进损伤模拟.然后采用渐进均匀化方法,以介观尺度单胞为媒介预测了0$^\circ$和90$^\circ$子胞的性能参数,并建立平纹编织复合材料的子胞模型,进而扩展成为材料的宏观尺度低速冲击模型. 在此基础上,研究了平纹编织复合材料低速冲击下的力学响应与损伤特征.结果表明:宏观冲击仿真和试验吻合较好, 验证了多尺度方法的正确性;最大接触力、材料吸能和分层面积均随冲击能量的增大而增大,分层损伤轮廓逐渐从椭圆形向圆形转化;基体拉伸和压缩损伤的长轴方向分别与子胞材料主方向正交和一致,损伤面积前者远大于后者.      

基于多尺度方法的三维编织C/C复合材料弯曲失效特性分析

为了预测三维编织C/C复合材料的弯曲失效行为,基于多尺度渐进展开理论,结合细观渐进损伤模型,建立了三维编织C/C复合材料宏细观多尺度分析模型。通过商业有限元软件ABAQUS用户子程序UMAT的二次开发,在宏观结构有限元分析中实时调用细观单胞模型进行细观渐进损伤分析,实现了宏细观尺度之间交互式信息传递和多尺度损伤模拟。利用上述模型对三点弯曲载荷下三维编织C/C复合材料梁的渐进损伤和失效过程进行了模拟,预测了梁的载荷-挠度曲线和弯曲强度,并与实验结果进行了对比分析,验证了基于多尺度方法的三维编织C/C复合材料弯曲强度预测模型的有效性,为此类材料及结构失效分析提供了一种手段。     

纤维增强复合材料弹性性能预测的域分解方法及应用

提出了新的有限元建模方法,即域分解方法,用于预测纤维增强复合材料单向带T300/BSL914C(环氧树脂)和AS4/3501-6(环氧树脂)的弹性性能。域分解方法基于区域叠合技术,分别建立单胞的整体域与纤维域模型用于代替传统有限元建模方法中单胞的基体域与纤维域模型。整体域是真实基体体积与纤维体积的叠加,两区域网格独立划分,互不影响。采用MSC.Nastran中的多节点约束Explicit单元,在整体域与纤维域节点之间建立位移连接属性模拟单胞基体域与纤维域之间的位移约束关系,从而实现两区域的耦合计算。计算结果表明:域分解方法单胞模型纤维增强方向弹性模量Ez预测值与试验值误差在7%以内,其余弹性常数也都与试验值吻合较好。域分解方法不仅可以大大简化纤维增强复合材料的细观力学建模,而且可以准确地预测纤维增强复合材料的弹性性能。     

混凝土多尺度应力响应方程及其数值模拟

针对混凝土的多相多尺度材料组成特征及其复杂力学响应问题, 首先, 根据混凝土中各组成材料的几何特征, 将C-S-H凝胶、硬化水泥浆体、砂浆及混凝土细观组成分别视为纳观、微观、亚细观和细观尺度上的复合材料, 并利用颗粒空间堆积方法, 重构了混凝土各尺度复合材料的简化几何模型; 其次, 基于重构的几何模型和等效夹杂理论, 通过等效刚度的升阶计算和应力响应的降阶计算, 建立各尺度复合材料应力响应之间的过渡关系, 推导混凝土多尺度应力响应方程, 并编制相应的计算程序; 最后, 以单轴压缩载荷作用为例, 数值计算载荷作用下混凝土各尺度复合材料中的应力响应, 分析骨料空间位置和相互作用以及水化产物刚度、几何形状和空间取向对其应力响应的影响规律. 结果表明, 单轴压缩载荷作用下, 混凝土细观组成中的应力分布并不均匀; 骨料颗粒之间的距离影响到混凝土中的应力分布, 其有效影响范围约为骨料粒径的6倍; 水泥水化产物的刚度、几何形状和空间取向是影响其应力分布的重要因素, 刚度越大, 所受应力越大, 与载荷作用方向的夹角越小, 长椭球形水化产物沿载荷作用方向的应力越大, 扁椭球形水化产物与之相反.      

三维编织复合材料几何建模及数值分析

论文根据三维编织复合材料的结构特点,把整体结构分为内部单胞、表面单胞和角单胞三种类型的子单胞,考虑空间纤维束的相互扭结和挤压所造成的纤维束的弯曲和截面变形,针对每种类型的子单胞,建立了相应的几何分析模型.引入周期性位移边界条件,建立了材料的弹性性能预报模型,得到了三维编织复合材料的工程弹性常数.通过数值比较可以看出,论文所给出的数值计算结果与实验结果吻合较好,从而验证了本文模型的有效性.     

基于均匀化理论的混凝土宏细观力学特性研究

在细观层次上将混凝土视为由砂浆基质、骨料及其界面组成的三相复合材料,在此基础上,采用有限元方法,对混凝土弹性本构关系进行数值模拟,利用位移渐近展开技术和均匀化理论建立了多尺度力学框架下的有限元平衡方程,重点考察了单胞尺寸对宏观力学性能的影响,得到了相对最大骨料粒径的最小单胞尺寸,即5～10倍最大骨料粒径.作为例证,对混凝土三点弯梁进行了宏细观尺度数值仿真研究,结果表明:基于本文方法可以较好地反映混凝土的宏观力学行为.此外,由于骨料、砂浆的相互作用,应力分布呈现出宏观渐变平滑与细观局部突变的特征.     

复合材料导向叶片的结构与材料一体化优化设计

将材料细观结构优化和宏观结构优化结合起来,从结构与材料两个尺度出发,发展复合材料导向叶片的结构与材料一体化优化设计方法。在建立平纹编织复合材料细观结构分析模型的基础上,采用细观力学有限元法进行材料刚度性能预测,热-固耦合分析方法进行涡轮导向叶片结构分析。同时以细观结构参数和宏观结构参数为设计变量,以叶片质量最小为优化目标,并要求满足应力和位移方面的约束,进行编织复合材料导向叶片的结构与材料一体化优化设计。     

Nonlinear viscoelastic multi-scale repetitive unit cell model of 3D woven composites with damage evolution

Three-dimensional (3D) textile composites have great potential applications to aircrafts and high speed vehicles because of the high strength/weight ratios and the capabilities of manufacturing complex, net-shape preforms. This paper reports the nonlinear viscoelastic responses and damage mechanisms of one kind of 3D textile composites, named as 3D orthogonal woven composite (3DOWC) under quasi-static tensile loading based on a micro/meso-scale repetitive unit cells (RUCs) model. In the RUCs model, the resin is described with a nonlinear viscoelastic material and the fibers/tows with an elastic material. The damage initiation and propagation in resin are simulated by the post-damage constitutive models with maximum principal theory failure criteria. The fibers/tows impregnated with resin are defined by elastic transverse-isotropic material model with ultimate strengths failure of ‘expanded smeared crack’ both along and perpendicular to fibers/tows axis direction. The engineering parameters and ultimate strengths of homogenized fibers/tows filled with matrix in meso-RUCs model are transferred from the numerical analysis of the micro-RUCs. The results are compared with experimental and theoretical values of RUC deformation and damage initiation and propagation under monotonic axial loading. The methodology of establishing the nonlinear visco-elastic multi-scale model of 3D textile composites without introducing the real fabric architecture in finite element analyses is explained. With the multi-scale RUCs model, the mechanical behaviors of other kinds of 3D textile composites can also be predicted.     

Two-scale analysis of a filament-wound cylindrical structure and application of periodic boundary conditions

A two-scale numerical approach to predict the effective in-plane properties of helical filament-wound thin-walled cylindrical tubes is provided. A meso-scale repeated unit cell (RUC) model for a filament-wound tube is established according to the manufacturing process, in which cross-overs, undulations and overlaps of fibre bundles are described using a bottom-up solid modelling technique. An approach to implement the general periodic boundary conditions in a finite element analysis scheme is also presented. As an application example, the effective in-plane elastic constants of glass fibre/epoxy filament-wound tubes are predicted and the influences of the number of the winding circuits and the shape of the fibre bundles are analyzed and discussed. In addition, the stress/strain distribution in the RUC is obtained which provides the essential information on the stress/strain concentration due to fibre cross-over/undulation/overlap. This then indicates the location where damage will initiate.     

Moving interfacial crack between piezoelectric ceramic and elastic layers

An analysis is performed for the problem of a finite Griffith crack moving with constant velocity along the interface of a two-layered strip composed of a piezoelectric ceramic and an elastic layers. The combined out-of-plane mechanical and in-plane electrical loads are applied to the strip. Fourier transforms are used to reduce the problem to a pair of dual integral equations, which is then expressed in terms of a Fredholm integral equation of the second kind. The dynamic stress intensity factor(DSIF) is determined, and numerical results show that DSIF depends on the crack length, the ratio of stiffness and thickness, and the magnitude and direction of electrical loads as well as the crack speed. In case that the crack moves along the interface of piezoelectric and elastic half planes, DSIF is independent of the crack speed.     

二维三轴编织复合材料的弹性性能分析

提出了二维三轴编织复合材料的几何模型,模型考虑了纤维束的弯曲扭转状态及空间交错特性等几何元素。基于体积平均法,建立了预测二维三轴编织复合材料弹性性能的理论分析模型;通过引入更普遍适用的周期性位移边界条件,结合二维三轴编织复合材料的细观实体结构,建立了分析其力学性能的有限元模型。两种模型预测结果均与试验结果吻合,证明了方法的合理有效性。分析了材料受载下的细观应力分布,并讨论了编织参数对材料性能的影响。研究表明,二维三轴编织复合材料轴向性能得到了增强,应力分布更均匀,编织角以及纤维体积含量对材料弹性性能影响较大。     

Elastic stability of inhomogeneous thin plates on an elastic foundation

We study the elastic stability of infinite inhomogeneous thin plates on an elastic foundation under in-plane compression. The elastic stiffness constants depend on the coordinate variable in the thickness direction of the plate. The elastic foundation is represented as a Winkler-type model characterized by linear and nonlinear spring constants. First we derive the Föppl–von Kármán equations by taking variations of the elastic strain energy. Next we develop the linear stability analysis of the plate under uniform in-plane compression and explicitly derive the critical loads and wave numbers for particular three cases. The effects of the material inhomogeneity, material orthotropy and loading orthotropy on the critical states are examined independently. Finally, we perform a weakly nonlinear analysis of the plate at the onset of the buckling instability. With the multiple scales method, the amplitude equations for the unstable modes that provide insight into the mode type and its amplitude are derived and then the effect of the material inhomogeneity on buckling modes are evaluated qualitatively.     

Hamilton体系下矩形薄板受抛物线压力载荷的屈曲分析

针对四边简支矩形薄板在两对边相向的非线性分布压力下的面内应力分布以及屈曲问题,应用弹性力学的Hamilton体系和Galerkin法进行了研究.基于弹性力学的平面矩形域Hamilton体系,根据辛本征向量展开解法,得到了对应于零本征值和非零本征值的含待定常数的实数型面内应力分布通解.依据必须满足的应力边界条件,导出了矩形薄板在抛物线分布载荷下的面内应力分布.考虑到应力分布表达式的复杂性,用完全的解析方法得到屈曲载荷是不可能的.因此,运用基于虚功原理的Galerkin法,根据四边简支矩形薄板弯曲的位移边界条件,给出了不同长宽比矩形薄板受抛物线分布载荷的屈曲临界载荷.通过与已有文献中DQ法给出的数值计算结果比较,表明了本文求解方法的有效性和正确性.基于所给出的结果,可望为解决矩形薄板在非线性分布载荷下的面内应力分布以及屈曲问题提供一种新的研究方法.     

In-plane modal testing of a free isotropic rectangular plate

In-plane vibration modes of an aluminum panel were experimentally identified from frequency response tests. Responses were measured on the panel edges and at selected locations on the panel surface. The measurements on the surface were made by attaching accelerometers oriented parallel to the panel plane. Resonance frequencies, relative damping ratios and mode shapes were established for the lowest 12 in-plane modes found in the frequency range between 1600 and 7000 Hz. A damping ratio of less than 0.05 percent of critical damping is proved to be valid for the aluminum panel. A finite element software was used to calculate 12 corresponding theoretical in-plane eigenfrequencies and mode shapes. An outline for a nondestructive procedure is suggested to estimate the input data for the elastic constants of an isotropic plate model. Two of the modes were used in analogy with the flexural vibration of beams and plates. The modes illustrate the deformation pattern including shear deformations, through the thickness, for the bending modes of thick beams or plates. The Rayleigh-Timoshenko theory also was used for the calculation of these two eigenfrequencies.     

用均匀化理论分析蜂窝结构的等效弹性参数

在线弹性范围内,根据均匀化理论,并结合有限元方法推导出适用于二维周期性结构的均匀化的有限元格式（Homo FEM）,计算出不同相对密度下的规则蜂窝结构的等效弹性模量Ee和泊松比νe.同时,利用蜂窝结构的代表胞元模型,用常规的有限元方法（FEM）模拟计算出相应的等效弹性参数.最后将两种数值计算结果与己有的理论公式进行了比较和分析讨论.结果表明:在考察的相对密度全场范围内（0～0.4）,HOmO FEM得到的蜂窝结构的 Ee和νe 与 FEM使用平面实体单元模拟计算得到的结果一致吻合,反映出 Homo FEM数值方法的客观准确性和可行性.而 Gibson公式和 W-K得到的等效弹性模量值 Ee只是在较小相对密度的情况下（小于0.15）与数值计算结果吻合.当结构的相对密度较大时,必须考虑胞棱附近区域由应力集中导致的复杂的应力和应变分布的影响.     

Non-linear dynamic response of a thin laminate subject to non-uniform thermal field

Non-linear dynamics behavior of a thin isotropic laminate in a simply supported boundary condition is studied for its response with both mechanical and thermal loads in effect. The thermal effects of both the in-plane and transverse non-uniform temperature variations in steady-state are considered. The equation of motion for the laminate deflection is reduced to the Duffing equation in a decoupled modal form by means of a generalized Galerkin's method. The stress field as a function of deflection and temperature variation is also obtained in a plane stress condition for its non-linear elastic behavior with von Karman strain field.For an exemplary laminated microstructure used as a printed wiring board, it is found that a high rise of the in-plane temperature increases the resonance frequency and could significantly increase the stresses of the lamina. The through thickness temperature variation has no significant effect on the deflection. Failure analysis is also made based on the composite failure criteria for a laminate to identify the critical mechanical and thermal loads.     

A crack in a piezoelectric layer bonded to a dissimilar elastic layer under transient load

Summary  A piezoelectric layer bonded to the surface of an elastic structure is considered. The piezoelectric and the elastic layers are infinite along the x-axis and have finite thickness in the y-direction. The polarization direction of the piezoelectric material is along the y-axis. By means of the method of singular integral equations, the solution in a Laplace transform plane is demonstrated. Laplace inversion yields the results in the time domain. Numerical values of the crack tip fields under in-plane transient electromechanical loading are obtained. The influence of layers thickness on stress and electric displacement intensity factors is investigated. Received 16 March 2000; accepted for publication 16 August 2000