木结构燕尾榫节点应力分析研究.docx

? ? 木结构燕尾榫节点应力分析研究* ? ? 张 俊 李源河 李英洁 徐伟涛 姚利宏 （1.内蒙古农业大学材料科学与艺术设计学院，内蒙古 呼和浩特 010018; 2.国家林业和草原局林产工业规划设计院，北京 100010) 木结构古建筑是一种以木构架为承重系统，墙体作围护结构的建筑[1]。榫卯节点是木构架的主要连接方式，其柱端与梁端分别制成卯口和榫头形式，将梁端的榫头插入柱顶预留卯口中[2]，形成一种凹、凸结合的特殊结构系统。在地震荷载下，榫卯节点多通过节点间摩擦滑移与挤压进行消震减震，是木结构耗能的主要部位[3-4]，具有一定的拉伸、压缩、弯曲和抗扭强度和半刚性特点[5-7]，其对木结构建筑整体稳定性起着至关重要的作用。 燕尾榫有“万榫之母”之称，是梁柱构件典型的节点连接形式之一。李义柱[8]试验发现燕尾榫节点的刚度与节点转角成反比，弯矩与节点转角成正比，其滞回环形状均呈反Z形，具有显著的“捏拢”和“滑移”效应。周乾等[9]对比燕尾榫与直榫发现，两者卯口处都有不同程度拉、压、剪破坏，且燕尾榫更为严重，榫头更易破坏。张令心等[10]发现燕尾榫节点榫颈与卯口处最易发生破坏，木材的材性和节点的形状均会影响破坏的形式。杨娜[11]首先构建了燕尾榫节点初始刚度与极限弯矩的回归模型，并依此推导出古建筑木结构整体结构计算的燕尾榫节点三参数幂函数模型。吴洋等[12]发现燕尾榫底部受力最为集中，额枋两端受力相对薄弱，拔榫是燕尾榫主要的破坏形式。罗勇[13]得到了准确的燕尾榫木构架弯矩-转角滞回曲线。本文在已有燕尾榫节点受力研究基础上，从模型比例与等材方面展开研究，利用ABAQUS对其进行拔榫量与应力分析，以期为后续燕尾榫节点受力分析研究提供参考。 1 燕尾榫节点受力机理 在燕尾榫节点连接梁枋时，因其限制了沿榫宽和榫深方向的自由度，只允许沿榫厚方向的移动，因此具有良好的接合强度。本文利用ABAQUS有限元分析软件模拟不同等材与不同比例的燕尾榫节点，分析其在单调加载和低周循环往复加载下的力学性能。 如图1 所示，假定燕尾榫节点向下转动为正方向[14]。榫头沿着榫颈转动时，由于上部与卯口挤压，受到来自卯口的压应力N2，同时还受到来自卯口与榫头上表面间的反向摩擦力f。榫头下部正向转动时，相对卯口滑移量较少，摩擦力可忽略不计，主要受与卯口之间(垂直于枋长方向)的斜向合力N3[15]。 图1 燕尾榫节点加载受力示意图Fig.1 Schematic diagram of the stress on the dovetail joint 2 有限元模型建立 2.1 木材材性定义 木材顺纹本构模型和横纹本构模型分别如图2a、b所示，其中横纹本构模型考虑了弱强化段[16]。图中εcu为木材顺纹极限压应变，εco为木材顺纹受压屈服应变，εtu为木材顺纹受拉屈服应变，Ecs为木材顺纹受压弹性模量，Ets为木材顺纹受拉弹性模量。εL为木材横纹受压屈服应变，εN为木材横纹受压极限应变，E1为木材横纹受压弹性段模量，E2为木材横纹强化阶段的受压弹性模量，σL为木材的横纹受压强度，σN为木材横纹强化段结束时的受压强度。 图2 木材本构模型Fig.2 Constitutive model for wood ABAQUS软件中木材弹性特征用工程弹性常数来定义。本次模拟试验选择樟子松材（Pinus sylvestris）为研究对象，其弹性工程常数见表1[17]。 表1 樟子松弹性工程常数Tab.1 Engineering constants of Pinus sylvestris var. mongolica Litv. 2.2 节点模型建立 依照《营造法式》“材分制”标准[19-20]，设计了模型比例1∶3.2、1∶2.2、1∶4.4的二等材燕尾榫节点（依次命名为YS-1、YS-2、YS-3）与模型比例1∶3.2三等材燕尾榫节点（YS-4）。燕尾榫节点缩尺模型尺寸通过现代国际基本单位与宋尺转换得出，节点尺寸示意图如图3所示。 图3 燕尾榫节点尺寸示意图Fig.3 Dimensional diagram of dovetail joint 2.3 摩擦参数定义 燕尾榫节点的榫头与卯口接触面上存在切向和法向作用。接触面之间为法向作用时，定义接触方式为“硬接触”[21]，不限制榫卯节点接触面间的压力。当两接触面相互分离压力为零时，作用其面上的接触约束被解除。接触面之间为切向作用时，榫卯节点接触面产生摩擦力。由于在柱顶局部平板枋施加了轴向均布荷载，为了保持局部平板枋与柱接触，因此两者之间采用绑定连接[22]。榫头与卯口部位，选择柱端卯口作为从面，榫头作为主面；柱顶端面与局部平板枋接触部位，主面选择柱顶端面，从面选择局部平板枋接触底面。 2.4 网格与单元定义 采用C3D8