薄层单元求解带初间隙摩擦接触问题研究.pdfVIP

·550· 薄层单元求解带初间隙摩擦接触问题 姚纬明，张春华，李同春 (河海大学水利水电工程学院．江苏南京210098) 摘要： 针对带初闻隙摩擦接触问题，本文提出了基于矩形或柱状常规薄层单元并以高斯点来描述接触状态的求解方法。 运用本文方法可以有效地模拟施工过程带来的接触缝问题，比常规方法计算精确度高。 关键词： 薄层单元：摩擦接触；劲度矩阵：赢斯点：接触半宽 Thin Elementfor ContactProblemswith Layer Solving FrictionandInitial Gaps Yao Chunhua,Li Weiming，ZhangTongchun and fCollegeofWaterConservancyHydropowerEng．，HohaiUniv，Nanjin9210098，China) atcontact withfrictionandinitial contactstatesimulationthe Abstract：Aimingproblems gaps，the byusing Gauss inthe orthecolumnelementis inthis contactduetOthe points rectangle presentedpaper．The joints and has thanthe methods． constructioncanbe simulated．it easilyaccurately highercomputationprecision regular words：thinelement：frictioncontact：stiffnessmatrix：Gauss half-width Key layer point：contact 1 前言 不同材料的迭合面(如坝与地基的接触面等)、结构开裂后的缝面等均属于接触问题。接触面上应力 条件一般需满足摩尔库仑屈服准则。若在接触面人为设置初始缝隙，其问题的求解归属于带初间隙的摩擦 接触问题。由于接触面的状态(闭合或张开)直接影响到结构整体应力分布，进而影响到结构整体安全度， 因此越来越受到人们的重视。 过研究结点与结点之间的接触状态来求解。当非线性接触变形较大时。必须转而采用结点与线段(二维) 或结点与面(三维)之间的接触形式进行研究。 接触问题的求解通常采用的方法主要有三类：第一类是罚有限元法，该方法的基本思想是将接触面的 非嵌入条件作为罚项吸收到系统的总势能泛函中，使原先的约束优化问题转化为罚优化问题，并通过调整 罚因子寻求真解：第二类方法为拉格朗日乘子法，这两类方法的特点均是通过附加约束方程用优化方法求 解：第三类方法称为变分不等式法，这类方法主要通过对结点接触状态进行迭代，将不等式约束转化等式 约束求解，由于这类方法不需增加求解变量，不改变刚度矩阵的对称性因而有一定的应用价值。 上述方法对接触面上接触状态的判断均是对结点进行的。虽然对一般接触问题的求解有效．但存在以 下三大缺点： (a)当接触区域为不连续曲面时，采用结点与结点或结点与线(面)接触模型无法模拟结点处的接 cn， 作者简介：媾纬明(1957)．男．江芬入