非连续变形计算力学模型基本原理及其线性规划解.pdfVIP

第40卷第3期 大 连 理 工 大 学 学 报 Vol 4O．No 3 2 0 0 0年 5月 JournM of Dalian University of Technology Mly 2 0 0 0 一 习 文章编号：1000-8608(2000)03—0351—07 非连续变形计算力学模型基本原理及其线性规划解 黎 勇 ， 栾 茂 田“ (1．大连理工大学 土术工程系，辽宁 大连 116024； -}U 2．大连理工大学 海岸和近海工程国家重点实验室，辽宁大连 116024) 丁U 摘要：采用Coulomb摩擦定律考虑界面的摩擦接烛案件，从一般的平面问题出发，根据其边值 方程的野形式，利用加报采量法．在分别建立广义有限单元的位穆插值模式和接触力元(印接触 面上由力分布模式)的基础上。推导了多体系统非连续变形计算力学模型的基丰控制方程t给出 了其墁性规划的教学列竞和求梓试问恿的数值算法．算帝1分析表明．谈方法能够模报多体相互 作用系统的变形和应力及接触内力． 关键词：有限元}线性规划；接触面／非连续变形；接舷力元；广义有限元 中图分类号：TU43；0343．3；0242．21 文献标识码：A 勃字援型 O 引 言 题，根据系统的平衡方程和摩擦接触、滑移和分离 等边值条件，给出了其平衡方程弱形式．进而根据 自从 2O世纪 6O年代有限元法产生与发展以 数学加权余量法，推导出了非连续变形计算力学模 来，相继建立了许多不同的离散模型与数值方法来 型的动力控制方程及其离散形式的数学规划列式． 分析岩土工程中的非连续变形问题，例如基于 1．1 平面问题边值方程的一般形式 Goodman无厚度接触单元和Desai薄层单元的界 假定多体系统由m个物体构成．已知t时刻系 面单元模型一有限元方法 ]，从Kawai的刚体一弹 统构形为n，每个物体构形为n．，边界为S(观)一 簧模型(RBSM)出发的刚体有限元方法【4 ]，到 U U ，其中 为位移边界． 为外力边界， Cundall的个别单元法0 和石根华的非连续变形分 析(DDA)0 所形成的各种离散单元法，以及将连续 为接触边界，n= ：n．，s(n)= ：s(皿)．边值 百 有限元与非连续变形分析相统一的流形方法 增 ． 问题表述为，在给定的边界条件下，求解t+At时 这些方法对接触界面的处理大都是通过施加刚性 刻多体系统的构形，即确定应力函数向量增量 法向和切向弹簧来实现的，通过改变弹簧的刚度来 ( ( ， ) ( ，Y) r( 。 )) 、位移函数向量增 模拟接触界面的力学特性，如接触界面两侧物体的 量(“ ， ) ， )) 和界面接触力向量增量 滑移、分离、闭合、再滑移和再分离现象，所给出