用线性互站方程求解空间弹性接触问题.pdfVIP

第一章基本理论 本章简要回顾了本文所需应用的基本知识。变分不等方程部分以一简单问 题为例，引出了变分不等方程，并介绍了与本文有关的变分不等方程理论；接 触问题部分建立了两个物体接触情况下，有无摩擦时的接触互补形式；线性互 补方程部分则给出线性互补方程的定义和本文采用的解法(Lemke算法)。 第一节 变分不等方程简介 变分不等方程是近代应用数学的新成果之一，近年来在力学和物理学中得 到广泛的应用。变分不等方程经数值方法离散，导出相应线性互补方程，面线 性互补方程已有一些较成熟的解法，如Lemke算法。 下面以一典型自由边界问题——～维障碍问题为例，介绍变分不等方程 的形成。 P L X D 『 』 —- T— 图1-1．1 例：一等截面柱体，截面积为A，弹性模量为E，长为L，上部固定，下 端自由，见(图1-1．1)。下端受有合力为P的均布法向力，距下端D远处有刚 性障碍，D((L，当P大到一定程度，即拉伸量大于或等于D之后，柱体伸 长将受到障碍的限制，同时障碍将对柱体下端有作用力。求得柱体下端的位移 6 ——————————点丛塑巡塑 和障碍阻力，就可以求出该问题的应力场和位移场。 由弹性力学知，柱内各点必须满足： 平衡方程 Gx．：=0 (1．1．1) 物理方程 o。=Es。 ¨小2 几何方程 占z=“1 ㈤H) 边界条件 u(o)=0． u{L1≤D “山4 且如果 “(上)D 则口，=PIA (1-1．5) “(三)=D贝0 Ox≤P／A 另外，障碍对柱体的作用力为： Q2P。仃x(L)A (1一l_61 为了导出所需要的变分不等方程，我们对方程(1-1．1)式运用虚位移原理 得 o=J[肛．如)旷(工)一“扛)-k幽 =f【甄巾)旷(小小)协出 =EA毛(“)【y(z)一H(x)】|：一f￡如，(“)【毛(矿)-sx(H)幽出 =(P—Q)【y(三)一(训一fl瞻(“)【‘(矿)一＆(越)】幽出 (】-1．7) 其中， “∈U耐， VV∈U“，u_为可能位移的集合． 定义：c(虬矿)=￡如，(Ⅳ)q(V)dn O-l-s) 户，u2PU(L) (1-l忉 其中Q表示全域，并注意到： 当Q0时．u(L)=D，V(L)sD，于是有 0 Q【V(L)一u(L)】S (I-I—10a) 当Q=0时．U(L)=D，V(L)sD，于是有 L