三维摩擦接触问题不动点算法.pdf

第 28卷 第 1期 力 学 季 刊 Vo1．28No．1 2007年 3月 CHINESEQUARTERLYOFMECHANICS March2007 三维摩擦接触问题不动点算法 祁德庆 ，高云开 ，徐连 民 (同济大学 建筑工程系 ，上海 200092) 摘要 ：三维摩擦接触 问题具有多重非线性性质 ，使求解变得 比较困难 ，为 了解决此问题 ，建立 了三维摩擦接触 问 题的模型，介绍了一种非光滑混合不动点算法。该算法克服了在接触面上 由于可能的滑动状态有无穷多个而难 以确定的难点 ，算法未引人人工变量 ，计算量较小，计算结果精确满足接触状态条件 ，收敛性得到保证。根据此 算法编制程序，将叠合悬臂梁算例数值计算结果与商用有限元软件进行比较，也表明了不动点算法的有效性。 关键词 ：摩擦接触 问题 ；不动点算法 ；收敛性 中图分类号：0343．3 文献标识码 ：A 文章编号：0254．0053(2007l01．149．4 FixedPointMethodonThreeDimensionalFrictionaIContactProblem De-qing，GAO Yun-kai，XU 。竹-min (DepartmentofConstructionEngineering，T0ngjiUniversity，Shanghai200092，China) Abstract：Themultiplenonlinearpropertiesofthreedimensionalfrictionalcontactproblemsmakeitsca1． culationdifficult．Modelofconatctproblemsandafixedpointmethodwereintroducedtosolvetheprob． 1ems．Thismethodovercomedthedifficultyofdeterminingtheslippingangleatcontactproblems，which wascausedbytheunceratindirectionsoftherelativeslipping．Noadditionalvariableisimported．Thenu． mericalresultsaccordwiththecontactstatesexactlyandtheconvergenceisguaranteed．Thismethodwas appliedtocomputethecantileverbeam．Thecomparisonoffixedpointmethodandfiniteelementsoftware suggeststhatthismethodishighlyefficient． Keywords：frictionalcontactproblems；fixedpointmethod；convergence 工程实际中存在着大量的摩擦接触问题 ，对于二维问题，由于接触面上只能沿着接触点的切线有2个 可能的滑动方向，已有了如线性互补口 等较有效的算法。但是对于三维问题 ，由于在接触面上存在无穷 多个可能的滑动方向，这为问题求解及算法收敛性分析带来很大麻烦，实际计算中得到大量应用的是迭代 法 ，该方法过程简单 、易于实施 ，但是其收敛性不能保证 。 文献E4]中提出了三维摩擦接触问题的一种混合不动点算法，本文结合该算法提出了三维接触问题的 模型和不动点算法求解过程 ，并在具体算例中将该算法与商用程序ABAQUS的接触单元做 出比较。 三维摩擦接触问题模型 假定 B 和 B。为 2个接触体 ，可能接触的边界为 r 和 r ，假定两个接触体只有小变形和小转动 ，且 接触边界充分靠近，U (a=1，2)为 I1边界上的位移，d。为已知 的初始间距，对于微小变形问题，依据距 离最近 的原则假定接触点对后 ，d。可由距离公式得出。则两接触体的相对间距 收稿 日期 2006—11—20 作者简介 ：祁德庆 (1960