钢-混凝土组合箱梁梁段有限元法研究

组合梁界面滑移将减小组合梁刚度,增大变形,影响构件性能;剪应力沿截面横向分布不均匀,造成其弯曲正应力的横向分布呈曲线形状,即剪力滞效应;同时组合梁往往重载,具有较小的跨高比,剪切变形不可忽略.根据虚功原理,建立了同时考虑滑移效应、剪力滞及剪切变形效应的组合梁单元刚度矩阵及等效节点力向量,并在此基础上编制了组合梁梁段有限元程序.利用本文程序对现有组合梁试件的混凝土顶板应力、钢梁底板应力、跨中挠度和梁端滑移进行了计算,并将本文计算结果与解析法计算结果及试验结果进行了比较.结果表明,本文计算结果与解析法计算结果及试验结果吻合良好;同时,本文计算结果具有较好的稳定性,验证了本文计算方法的正确性.本文所建立的梁单元刚度矩阵同时考虑了剪切变形、剪力滞及滑移效应的影响,符合工程实际,为有限梁段法分析组合箱梁提供了理论基础.     

A variational method for unidirectional fiber-reinforced composites with matrix creep

A variational method is developed for analyzing the matrix creep induced time-dependent change in fiber stress profiles in unidirectional composites. A functional of admissible profiles of fiber stress rate is presented by supposing a fiber broken in matrix as well as a fiber pulled out from matrix. The functional is shown to have the stationary function satisfying an incremental differential equation based on the shear lag assumption. Then, the stationary function is approximately determined by assuming bilinear profiles of fiber stress and a power law of matrix creep, leading to analytical solutions for the time-dependent change in fiber stress profiles. The solutions are verified on the basis of an energy balance equation and a finite difference computation. Moreover, it is shown that the solution for the fiber pull-out model agrees well with an experiment on a single carbon fiber/acrylic model composite if the initial slip at fiber/matrix interface is taken into account. In addition, the solution for the fiber breakage model is used for evaluating the characteristic time in long-term creep rupture of unidirectional composite.     

Thin interphase/imperfect interface in elasticity with application to coated fiber composites

The imperfect interface conditions which are equivalent to the effect of a thin elastic interphase are derived by a Taylor expansion method in terms of interface displacement and traction jumps. Plane and cylindrical interfaces are analyzed as special cases. The effective elastic moduli of a unidirectional coated fiber composite are obtained on the basis of the derived imperfect interface conditions. High accuracy of the method is demonstrated by comparison of solutions of several problems in terms of the imperfect interface conditions or explicit presence of interphase as a third phase. The problems considered are transverse shear of a coated infinite fiber in infinite matrix and effective transverse bulk and shear moduli and effective axial shear modulus of a coated fiber composite. Unlike previous elastic imperfect interface conditions in the literature, the present ones are valid for the entire range of interphase stiffness, from very small to very large.     

Multiple cracking of fiber/matrix composites—Analysis of normal extension

This paper is concerned with the fracture of a fiber embedded in a matrix of finite radius. There is a periodic array of cracks in the fiber along the central axis of the medium. The paper accounts for the cases of axial extension and residual temperature change of the composite medium. Fourier and Hankel transforms are used to reduce the problem to the solution of a system of dual integral equations, which are solved by the singular integral equation technique. Rigorous fracture mechanics analysis, which exactly satisfies all boundary conditions of the problem, is conducted. Numerical solutions for the crack tip field and the stress in the fiber are obtained for various values such as crack radius, crack spacing and fiber volume fraction.     

纤维增强橡胶复合材料界面力学性能有限元分析

基于剪滞理论,引入双线性内聚力模型研究了纤维与基体界面应力传递机理.采用ABAQUS模拟了非理想界面在单纤维拔出过程中的脱粘失效,分析了不同脱粘阶段界面剪应力分布情况,以及界面刚度和纤维长径比对界面应力传递和拔出载荷的影响规律.结果表明,在纤维受载失效过程中,纤维的拔出过程可分为4个阶段,即界面的完全粘结、损伤演化、逐渐脱粘、完全脱粘.界面的刚度和纤维长径比对界面应力传递与最大拔出力均有一定的影响.界面刚度、纤维长径比主要影响纤维的最大拔出载荷以及界面脱粘失效位移.     

Dynamic fiber debonding and frictional push-out in mode composite systems: Experimental observations

In the present work a modified Split Hopkinson Pressure Bar (SHPB) system is adopted to perform dynamic fiber push-out experiments on model single fiber composite systems. A tapered punch and a support connect a monofilament composite with the incident and transmitted bars of the SHPB. The tapered punch is used to apply compressive loading to a single fiber (either steel or aluminum) embedded in a surrounding matrix material (EPON 862). The SHPB allows real time measurement of relative fiber/matrix displacement and push-out force, as the debonding and push-out event progresses. Using this technique we have studied the effect of loading rate, material mismatch, fiber length, and surface roughness on the push-out event. It was seen that maximum push-out force increases with increasing loading rate. In addition dynamic interfacial strength and toughness is highly dependent on fiber surface roughness. Results from a finite element analysis incorporating a cohesive failure model were used to extract interface strength and toughness values. It was found that the particular aluminum/EPON interface used is characterized by a dynamic shear failure strength of 48±8 MPa, a mode II fracture toughness of 160±40 N/m, and a friction coefficient of 0.2 at a sliding rate of 6 m/s. For the rates tested here these quantities were found to be approximately constant.     

Stiffness properties of particulate composites containing debonded particles

This paper considers the problem of determining the nonlinear bimodular stiffness properties, i.e., the tensile and compressive Young’s moduli and Poisson’s ratios, and the shear modulus, of particulate composite materials with particle–matrix interfacial debonding. It treats the general case in which some of the particles are debonded while the others remain intact. The Mori–Tanaka approach is extended to formulate the method of solution for the present problem. The resulting auxiliary problem of a single debonded particle in an infinite matrix subjected to a remote stress equal to the average matrix stress, for which Eshelby’s solution does not exist, is solved by the finite element method accounting for the particle–matrix separation and contact at the debonded particle–matrix interface. Because of the nonlinear nature of the problem, an iterative process is employed in calculating the stiffness properties. The predicted stiffness properties are compared to the exact solutions of the stiffness properties of particulate composites with body-centered cubic packing arrangement.     

纤维增强韧性基体界面力学行为

分析了纤维增强韧性基体的界面力学行为及其失效机理.按剪滞理论和应变硬化规律研究微复合材料的弹塑性变形和应力状态.讨论了幂硬化和线性硬化基体的弹塑性变形和界面应力分布,并给出纤维应力和位移的表达式.按最大剪应力强度理论建立了纤维/基体界面失效准则,推导出弹塑性界面失效的平均剪应力随纤维埋入长度的变化关系.     

Debonding and fiber pull-out in reinforced composites

In order to evaluate the strength of fiber-reinforced composites, there is first the need to investigate the interfacial debonding and the pull-out of fibers in a fractured composite with intact fibers. This type of problem in crack bridging has been investigated by several authors based on different models and assumptions [1–7]. In this study, we will consider a three-dimensional model of a single fiber of finite length bonded by a finite cylindrical matrix with an initial crack existing in a portion of the interface. In the model, one end of the cylinder is so constrained that the axial component of displacement vanishes. A tensile stress is applied to the fiber at the other end. The aim is to determine the pull-out of the fiber and the critical condition for interfacial debonding. Both the fiber and the matrix are treated as elastic materials. Analysis is made based on a method using Papkovich-Neuber displacement potential functions for the problem of an elastic solid subjected to axisymmetrical boundary conditions. Solutions are found by means of the technique of trigonometrical series. Effects of initial misfit strains and frictional sliding between the fiber and the matrix over the interfacial crack are also included in the study.     

肋板结构受压构件的剪力滞效应分析

建立了肋板结构受压构件考虑剪力滞效应时的微分方程,并以该微分方程的解作为位移插值函数,得到了较为精确的有限段单元刚度矩阵,数值算例证明本方法具有好的精度,且具有一维有限元法的便利和广泛适应性。     

Cohesive modeling of transverse cracking in laminates under in-plane loading with a single layer of elements per ply

This study aims to bridge the gap between classical understanding of transverse cracking in cross-ply laminates and recent computational methods for the modeling of progressive laminate failure. Specifically, the study investigates under what conditions a finite element model with cohesive X-FEM cracks can reproduce the in situ effect for the ply strength. It is shown that it is possible to do so with a single element across the thickness of the ply, provided that the interface stiffness is properly selected. The optimal value for this interface stiffness is derived with an analytical shear lag model. It is also shown that, when the appropriate statistical variation of properties has been applied, models with a single element through the thickness of a ply can predict the density of transverse matrix cracks.     

Three-dimensional continuum analysis for a unidirectional composite with a broken fiber

The three-dimensional problem of a periodic unidirectional composite with a penny-shaped crack traversing one of the fibers is analyzed by the continuum equations of elasticity. The solution of the crack problem is represented by a superposition of weighted unit normal displacement jump solutions, everyone of which forms a Green’s function. The Green’s functions for the unbounded periodic composite are obtained by the combined use of the representative cell method and the higher-order theory. The representative cell method, based on the triple discrete Fourier transform, allows the reduction of the problem of an infinite domain to a problem of a finite one in the transform space. This problem is solved by the higher-order theory according to which the transformed displacement vector is expressed by a second order expansion in terms of local coordinates, in conjunction with the equilibrium equations and the relevant boundary conditions. The actual elastic field is obtained by a numerical evaluation of the inverse transform. The accuracy of the suggested approach is verified by a comparison with the exact analytical solution for a penny-shaped crack embedded in a homogeneous medium. Results for a unidirectional composite with a broken fiber are given for various fiber volume fractions and fiber-to-matrix stiffness ratios. It is shown that for certain parameter combinations the use of the average stress in the fiber, as it is employed in the framework of the shear lag approach, for the prediction of composite’s strength, leads to an over estimation. To this end, the concept of “point stress concentration factor” is introduced to characterize the strength of the composite with a broken fiber. Several generalizations of the proposed approach are offered.     

Transverse conductivity and longitudinal shear of elliptic fiber composite with imperfect interface

The paper addresses the problem of calculating the local fields and effective transport properties and longitudinal shear stiffness of elliptic fiber composite with imperfect interface. The Rayleigh type representative unit cell approach has been used. The micro geometry of composite is modeled by a periodic structure with a unit cell containing multiple elliptic inclusions. The developed method combines the superposition principle, the technique of complex potentials and certain new results in the theory of special functions. An appropriate choice of the potentials provides reducing the boundary-value problem to an ordinary, well-posed set of linear algebraic equations. The exact finite form expression of the effective stiffness tensor has been obtained by analytical averaging the local gradient and flux fields. The convergence of solution has been verified and the parametric study of the model has been performed. The obtained accurate, statistically meaningful results illustrate a substantial effect of imperfect interface on the effective behavior of composite.     

剪力滞效应对箱形梁挠度影响的研究

从剪力滞翘曲应力的轴向平衡条件出发,选取双室箱梁的合理翘曲位移函数,引入相应于剪力滞翘曲变形的惯性矩和惯性积等几何特性,用能量变分法建立薄壁箱梁剪力滞效应分析的控制微分方程。通过求解控制微分方程,导出集中荷载和均布荷载作用下简支箱梁和悬臂箱梁的挠度公式及有限梁段单元刚度矩阵,模型试验和ANSYS壳单元计算结果证实了其正确性。结合简支、悬臂和连续箱梁数值算例,具体分析剪力滞效应对箱梁挠度的提高程度。结果表明,无论在集中荷载还是均布荷载作用下,剪力滞效应对简支箱梁的挠度均有显著提高。在集中荷载作用下,剪力滞效应对连续箱梁挠度的提高可达14%;对于跨宽比约为4.0～6.0的简支箱梁,可将按初等梁计算的跨中挠度乘以提高系数1.05～1.11;计算悬臂箱梁的挠度时,一般可以忽略剪力滞效应的影响。     

Fiber push-out study of a copper matrix composite with an engineered interface: Experiments and cohesive element simulation

The fiber push-out test is a basic method to probe the mechanical properties of the fiber/matrix interface of fiber-reinforced metal matrix composites. In order to estimate the interfacial properties, parameters should be calibrated using the measured load–displacement data and theoretical models. In the case of a soft matrix composite, the possible plastic yield of the matrix has to be considered for the calibration. Since the conventional shear lag models are based on elastic behavior, a detailed assessment of the plastic effect is needed for accurate calibration. In this paper, experimental and simulation studies are presented regarding the effect of matrix plasticity on the push-out behavior of a copper matrix composite with strong interface bonding. Microscopic images exhibited significant local plastic deformation near the fibers leading to salient nonlinear response in the global load–displacement curve. For comparison, uncoated interface with no chemical bonding was also examined where the nonlinearity was not observed. A progressive FEM simulation was conducted for a complete push-out process using the cohesive zone model and inverse fitting. Excellent coincidence was achieved with the measured push-out curve. The predicted results confirmed the experimental observations.     

Effective mechanical properties of unidirectional composites in the presence of imperfect interface

In this paper, the equivalent inclusion method is implemented to estimate the effective mechanical properties of unidirectional composites in the presence of an imperfect interface. For this purpose, a representative volume element containing three constituents, a matrix, and interface layer, and a fiber component, is considered. A periodic eigenstrain defined in terms of Fourier series is then employed to homogenize non-dilute multi-phase composites. In order to take into account the interphase imperfection effects on mechanical properties of composites, a stiffness parameter in terms of a matrix and interphase elastic modulus is introduced. Consistency conditions are also modified accordingly in such a way that only the part of the fiber lateral stiffness is to be effective in estimating the equivalent composite mechanical properties. Employing the modified consistency equations together with the energy equivalence relation leads to a set of linear equations that are consequently used to estimate the average values of eigenstrain in non-homogeneous phases. It is shown that for composites with both soft and hard reinforcements, largest stiffness parameter that indicates complete fiber–matrix interfacial debonding causes the same equivalent lateral properties.     

A time dependent micro-mechanical fiber composite model for inelastic zone growth in viscoelastic matrices

In this work, a fiber composite model is developed to predict the time dependent stress transfer behavior due to fiber fractures, as driven by the viscoelastic behavior of the polymer matrix, and the initiation and propagation of inelastic zones. We validate this model using in situ, room temperature, micro-Raman spectroscopy fiber strain measurements. Multifiber composites were placed under constant load creep tests and the fiber strains were evaluated with time after one fiber break occurred. These composite specimens ranged in fiber volume fraction and strain level. Comparison between prediction and MRS measurements allows us to characterize key in situ material parameters, the critical matrix shear strain for inelastic zones and interfacial frictional slip shear stress. We find that the inelastic zone is predominately either shear yielding or interfacial slipping, and the type depends on the local fiber spacing.     

短纤维金属基复合材料应力传递的有限元分析

基于短纤维增强金属基复合材料的单纤维轴对称和三维细观力学模型,利用弹塑性有限元分析方法对该复合材料中基体与纤维间的应力传递进行研究,研究中主要讨论了基体、纤维和界面的力学性能以及纤维位向的变化对应力传递和应力分布的影响。研究表明,复合材料微结构参数的变化将显著影响基体与纤维间的应力传递和复合材料中的应力分布,复合材料设计过程中必须考虑合理的微结构特征。     

单纤维压出实验模型的力学分析

单纤维压出实验是一种细观实验方法,用于测量复合材料的界面强度。本文采用子域法处理多域问题,利用轴对称边界元程序,对单纤维压出实验模型的应力传递进行了分析,并与常用的剪切滞后理论的分析结果进行了比较,发现在大部分纤维埋置区域,两种分析吻合很好,在纤维与基体的界面端部,边界元的分析结果显示界面应力在该区域有奇异现象,而剪滞理论则无法反映应力奇异现象。应力奇异现象的存在使得我们对该实验的界面强度判据需要     

单纤维段裂试验评述

单纤维段裂试验作为复合材料界面剪切强度的一种测试方法被沿用至今.但是, 这种方法的可信度已受到一些研究者的质疑.为了明确单纤维段裂试验的问题, 本文首先对试验技术、试验结果分析等方面作了概述, 并指出: 纤维段裂的饱和状态是单纤维段裂试验的终点标志,以及临界长度是由试验得到的唯一数据, 而这二点是这种试验方法独具的特点, 同时也是这种试验方法难以克服的缺陷.在单纤维段裂试验中, 按照纤维段界面端处的局部损伤模式, 有3种界面端应力奇异性分析的问题需要予以考虑:(1)纤维断裂, 基体没有开裂, 和界面没有脱粘;(2)纤维断裂, 基体开裂, 但界面没有脱粘;(3)纤维断裂, 界面脱粘, 基体已开裂或基体未开裂.在单纤维段裂试验的界面端应力奇异性分析的基础上, 本文对单纤维段裂试验的可靠性进行了研究.结论是: 任何纤维和基体组成的复合材料的单纤维段裂试验都存在界面端应力奇异性, 这就排除了用单纤维段裂试验测定界面剪切强度的可能性.