膜-基复合材料界面剪应力分析的影响系数法

提出了一种计算膜 基复合材料界面剪应力的新方法。将薄膜、基体分别作为二维和三维问题进行有限元计算 ,利用影响系数 ,以剪应力为未知量 ,求出同一点处薄膜、基体的位移。通过薄膜、基体在界面上位移相等的条件建立线性方程组 ,联立求解 ,计算出界面剪应力。结果表明 ,本文提出的方法思路清晰 ,易于实现 ,精度令人满意 ,为膜 基复合材料结构强度分析提供了一种新的途径     

含裂粘补结构胶层剪应力的计算

提出了一种计算含裂粘补结构胶层剪应力的新方法。采用有限元计算补片的影响系数，使用已有解析解计算蒙皮的影响系数。利用位移协调条件，列出有限多组联立方程。解出含裂粘补结构胶层剪应力。计算结果较好的反映了含裂粘补结构胶层剪应力的实际变化趋势，该方法同样可应用于其他结构如高梯度基膜复合材料剪应力的计算。     

涂层厚度对硬脆涂层结合强度的影响

应用有限元分析的方法分析二维模型涂层层数以及涂层厚度对涂层/基体结合强度的影响已得到广泛研究。然而二维模型相比三维模型未能考虑宽度方向应力应变的影响,导致求解精度低。因此对三维模型涂层层数以及厚度对涂层/基体结合强度影响的研究显得实为重要,但这方面研究比较少。本文建立有限元三维模型,通过比较有限元模拟结果与赫兹接触精确解析解验证了模型的可靠性,讨论了涂层厚度对涂层/基体界面最大剪切力影响关系。结果表明,剪应力最小值随着涂层厚度的增加逐渐减小,然而剪应力最大值却随着厚度的增加先减小后增大,在0.4mm时取得最小值;随着涂层厚度的增加,涂层/基体结合处剪应力Z方向突变减小,涂层内部剪应力Z方向突变增大,但是突变位置向涂层表面靠近,而且涂层/基体结合处X方向的剪应力变化随着涂层厚度增加趋于平缓。     

复合材料层合壳有限元分析的预测-修正法

对于复合材料层合壳的有限元分析，本文根据Reissner-Mindlin型的全局位移场给出了一个预测一修正法。首先按照一般的有限元分析过程(没有引入剪切修正系数)计算出全局响应(如挠度，频率和屈曲载荷等)的预测值；然后利用Lagrange插值构造横向剪应力的一般形式，使得满足层间连续和表面上为零的条件，通过最小二乘法拟合三维应力平衡方程获得横向剪应力；最后在单元上计算和引入剪切修正系数，再经过有限元分析计算出全局响应的修正值。     

复合材料层合圆柱壳的有限元分析

本文根据Reissner-Mindlin型的全局位移场(一阶和三阶)，应用有限元预测一修正法，数值计算和分析了机械载荷作用下复合材料层合圆柱壳的挠度和横向剪应力。首先按照一般的有限元分析过程(没有引入剪切修正系数)计算出层合圆柱壳的挠度预测值；然后利用Lagrange插值构造横向剪应力的一般形式，使得满足层间连续和表面上为零的条件，通过最小二乘法拟合三维应力平衡方程获得横向剪应力；最后在单元上计算和引入剪切修正系数，再经过有限元分析计算出层合圆柱壳的挠度修正值。数值计算结果与三维线弹性解的比较表明，挠度修正值和横向剪应力的精度是十分满意的。     

基于"梁段"跨中FRP-混凝土界面剪应力数值计算

纤维复合材料(FRP)加固混凝土梁跨中可能发生剥离破坏,研究者已提出几种计算FRP-混凝土粘结界面剪应力τ的方法,但所得结果有较大的差异.本文先通过典型的FRP-混凝土界面的粘结-滑移本构关系和FRP加固梁实验结果,评估既有计算方法存在的问题.进而提出一个简单的"梁段"有限元计算模型,并证实其合理性和可靠性.再用此模型计算讨论了多种参数影响下的FRP-混凝土界面剪应力τ分布特征.基于实验研究和有限元计算结果,指出FRP-混凝土界面剪应力最大值τu不适合作为剥离破坏准则,并建议考虑将滑移量最大值δ作为FRP-混凝土界面的剥离破坏准则.     

共晶基陶瓷复合材料的强度模型

根据细观结构内界面的强约束特性,通过纤维-基体内界面切应力确定了共晶陶瓷棒体的细观应力场.然后分析了两相界面处位错塞积产生的应力集中,获得基体内的最大应力,当最大拉应力等于基体理论断裂强度时,得到共晶棒体的断裂强度的解析表达式.考虑共晶陶瓷棒体长度和方位的随机性,根据概率理论得到共晶陶瓷基复合材料的宏观强度的理论模型.结果表明复合材料的宏观强度与亚微米纤维的直径和长度、以及亚微米纤维、基体、共晶陶瓷棒体的弹性常数有关.理论与实验结果十分接近,说明文中理论模型是合理的,同时证明了共晶界面对陶瓷复合材料的重要影响.     

纤维增强混凝土材料的界面剪应力分布研究

针对纤维与混凝土界面的破坏过程,提出了几种简化的粘结-滑移本构模型,以双线性局部粘结-滑移本构模型为基础,在受力平衡和变形协调的基本原理基础上,推导了纤维脱粘过程中界面剪应力的解析解.采用弹簧粘结单元,通过数值方法模拟了纤维与混凝土之间的粘结-滑移过程,给出了纤维与混凝土界面脱粘过程中界面剪应力的分布、变化情况.对解析解、有限元计算结果和试验结果之间的差异进行了对比分析,验证了简化模型的合理性和有效性.     

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

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

对矩形截面杆弹性自由扭转的探讨

利用柱体扭转问题的经典弹性力学解析解，结合自编电算程序绘制了沿截面控制线 上剪应力的分布图,同时利用ANSYS有限元分析软件模拟了等直矩形截面杆的自由扭 转问题. 将各种材料力学教材中的截面剪应力分布图与准确计算结果进行了比较，并 对截面剪应力分布规律进行了总结.     

Stress-function variational method for interfacial stress analysis of adhesively bonded joints

High interfacial stresses at the free edges of adherends are responsible for the debonding failure of adhesively bonded joints (ABJs). In this paper, a general stress-function variational method is formulated to determinate the interfacial shear and normal (peeling) stresses in ABJs in high accuracy. By extending authors’ prior work in stress analysis of bonded joints (Wu and Jenson, 2011), all the planar stress components in the adherends and adhesive layer of an ABJ are expressed in terms of four unknown interfacial stress functions, which are introduced at the upper and lower surfaces of the adhesive layer. A set of governing ordinary differential equations (ODEs) of the four interfacial stress functions is obtained via minimizing the complimentary strain energy of the ABJ, which is further solved by using eigenfunctions. The obtained semi-analytic stress field can satisfy all the traction boundary conditions (BCs) of the ABJ, especially the stress continuity across the bonding lines and the shear-free condition at the ends of adherends and adhesive layer. As an example, the stress field in an adhesively single-sided strap joint is determined by the present method, whose numerical accuracy and reliability are validated by finite element method (FEM) and compared to existing models in the literature. Parameter studies are performed to examine the dependencies of the interfacial stresses of the exemplified ABJ upon the geometries, moduli and temperature change of the adherends and adhesive layer, respectively. The present method is applicable for scaling analysis of joint strength, optimal design of ABJs, etc.     

Finite element analysis of interfacial stresses in FRP–RC hybrid beams

The interfacial stresses in fiber reinforced plastic (FRP)–reinforced concrete (RC) hybrid beams were studied by the finite element method. The mesh sensitivity test shows that the finite element results for interfacial stresses are not sensitive to the finite element mesh. The finite element analysis then is used to calculate the interfacial stress distribution and evaluate the effect of the structural parameters on the interfacial behavior. It is shown that both the normal and shear stresses at the interface are influenced by the material and geometry parameters of the composite beam. This research is helpful for the understanding on mechanical behavior of the interface and design of the FRP–RC hybrid structures.     

Non-steady-state elastic–plastic behavior of metal forming (rolling) with a generalized friction law

An analytical and generalized friction law is formulated and a three-dimensional large-strain non-steady-state elastic–plastic finite element analysis has been performed for rolling process. The contact/friction problem at the interface between the workpiece and the rolls is treated rigorously by using this new friction law. The numerical results, including the evolutions of roll torque and roll forces, and the interfacial normal and shear stresses, are presented and discussed.     

An axisymmetric stress analysis in a single fibre composite of finite length under a thermal expansion mismatch

In this paper, an exact solution is derived for the characterization of thermal stresses in a single-fibre composite of finite length. All the required boundary and interfacial conditions of the thermo-elastic problem are thus satisfied exactly. The proposed method involves a particular solution that is added to a three-dimensional (3D) complementary displacement field which satisfies automatically the Navier's equations. Based on experimental data provided by fibre Bragg grating sensors, an axisymmetric analysis is used then to determine the residual stress field inside the composite due to matrix shrinkage. The numerical results clearly indicate that all stress components vary significantly near the ends. An abrupt change of the shear stresses is thus predicted close to the edges. The results of the model are also found to be in good agreement with those obtained from finite element simulations. A comparison of the proposed approach with three other published theoretical models is also presented.     

Cohesive zone modeling of interfacial stresses in plated beams

The elastic analysis of interfacial stresses in plated beams has been the subject of several investigations. These studies provided both first-order and higher-order solutions for the distributions of interfacial shear and normal stresses close to the plate end in the elastic range. The notable attention devoted to this topic was driven by the need to develop predictive models for plate end debonding mechanisms, as the early models of this type adopted debonding criteria based on interfacial stresses. Currently, approaches based on fracture mechanics are becoming increasingly established. Cohesive zone modeling bridges the gap between the stress- and energy-based approaches. While several cohesive zone analyses of bonded joints subjected to mode-II loading are available, limited studies have been conducted on cohesive zone modeling of interfacial stresses in plated beams. Moreover, the few available studies present complex formulations for which no closed-form solutions can be found. This paper presents an analytical cohesive zone model for the determination of interfacial stresses in plated beams. A first-order analysis is conducted, leading to closed-form solutions for the interfacial shear stresses. The mode-II cohesive law is taken as bilinear, as this simple shape is able to capture the essential properties of the interface. A closed-form expression for the debonding load is proposed, and the comparison between cohesive zone modeling and linear-elastic fracture mechanics predictions is discussed. Analytical predictions are also compared with results of a numerical finite element model where the interface is described with zero-thickness contact elements, using the node-to-segment strategy and incorporating decohesion and contact within a unified framework.     

一种预测层合板层间剪应力的新后处理方法

层合板是航空航天领域典型的承力构件,过大的层间应力是导致其分层失效的主要原因.准确的层间应力预测往往依赖于三维平衡方程后处理方法(TPM).然而,该方法需要计算面内应力的一阶导,使得基于C0型板理论构造的线性单元无法使用TPM计算横向剪应力.本文在三维平衡方程后处理方法的基础上,提出了一种新后处理方法(NPM).新后处理方法通过虚功等效法消除了三维平衡方程后处理方法中产生的位移参数的高阶导.基于提出的新后处理方法和C0型板理论,仅需使用线性单元就可以预测层合板的横向剪应力.为了验证所提方法的有效性,本文基于修正锯齿理论(RZT)和所提方法构造了一种C0连续的三节点三角形线性板单元.数值算例表明,所提方法和三维平衡方程后处理方法具有相同的计算精度,提出的板单元能够准确高效地预测层合板的横向剪应力.此外,所提方法便于结合现有的有限元商用软件使用,基于商用软件中板壳单元获得的节点位移,使用新后处理方法极易获得准确的层间剪应力.     

Interfacial stresses in curved members bonded with a thin plate

The use of steel plates or externally bonded fibre-reinforced polymer laminates for flexural strengthening of concrete, masonry, timber or metallic structures is a technique that has become very popular. The effectiveness of this technique hinges heavily on the performance of the bond between the strengthening plate and the substrate, which has been the subject of many existing studies. In particular, the interfacial stresses between a beam and a soffit plate within the linear elastic range have been addressed by numerous analytical investigations. Surprisingly, none of these investigations has examined interfacial stresses in members with a curved soffit, despite that such members are often found in practice. This paper presents an analytical model for the interfacial stresses between a curved member of uniform section size and a thin plate bonded to its soffit. The governing differential equations for the interfacial shear and normal stresses are formulated and then solved with appropriate simplifying assumptions. Two numerical examples are presented to illustrate the effect of the curvature of the member on the interfacial stress distributions in a simply supported curved beam for the two cases of a point load and a uniformly distributed load. The analytical solution is verified by comparing its predictions with those from a finite element model.     

一类基于Hamilton体系的半解析法

利用偏微分方程在Hamilton体系中的表示和二类变量变分原理,结合有限元法,提出一类基于Hamilton体系的半解析法。本文以二阶非齐次椭圆型方程为例,给出了这类半解析法中的一种的有限元列式和算例,与解析解、Ritz法、有限条法和有限元法的结果比较表明,此法具有较高的精度,还可求解其它其些偏微分方程,有一定的普遍意义。     

A NEW SEMI-ANALYTIC ALGORITHM OF NEARLY SINGULAR INTEGRALS IN HIGH ORDER BOUNDARY ELEMENT ANALYSIS OF 2D ELASTICITY

针对边界元法中高阶单元中几乎奇异积分计算难题,解剖了二维边界元法高阶单元的几何特征,定义源点相对高阶单元的接近度。将高阶单元上奇异积分核函数用近似奇异函数逼近,从而分离出积分核中主导的奇异函数部分,其奇异积分核分解为规则核函 数和奇异核函数两项积分之和。规则核函数用常规高斯数值积分,再对奇异核函数积分导出解析公式,从而建立了一种新的半解析法,用于高阶边界单元上几乎强奇异和超奇异积分计算。给出3个算例,采用边界元法高阶单元的半解析法计算了弹性力学薄体结构和近边界点位移/应力,并与线性边界元正则化算法结果作了比较,结果表明提出的二次元的半解析算法更加有效。特别是分析薄体结构,采用正则化算法的线性边界元分析比有限元有显著优势,而用提出的二次边界元半解析算法分析比其线性元的有效接近度又减小了4个量级。