破裂过程的数值方法研究.pdfVIP

破裂过程的数值方法 赵纪生12刘艳琼1陶夏新1’3师黎静1 (1中国地震局工程力学研究所哈尔滨150080；2冻土工程国家重点实验室兰州730000； 3哈尔滨工业大学土木工程学院哈尔滨150080) 摘要 本文应用势能原理导出了位移间断萌生条件，联合物质导数描述的位移间断进化条件和间断连 接，对控制方程的弱形式进行了有限元数值离散。 关键词 位移间断；势能原理；离散；积分 1问题的描述 一个固体的分析域体为¨它的应力边界为，位移边界为瓯，平衡方程为： V盯+f=p／／ (1) 本构关系为： 仃=De或 寸=D毒 (2) 式中，弹性阶段时，D=D。；弹塑性阶段时，D=D。一Dp。Dp与材料的加载条件、塑性势函数 和塑性刚度相关。 分析固体破坏时的强度条件为 或 (3) 痧(盯，叹(鼋))=0 ．厂(盯)一嚷(g)=0 我们的目的是在拟静力加载条件下，建立满足式(1)～(3)和初边值条件的位移间断破坏过程的 数值分析方法(图1)。其中包括间断的萌生、生长，间断之间的连接和位移间断进化三个主要部分。 图1固体的破裂过程 ·-——161·-·—- 2位移间断条件 拟静力加载时，固体内部出现位移间断，即“=万+I“1日。，间断面两侧作用力均可视为外力， 相应的势能原理【11为， ll吐广删一budV—L矿郦一工，加+衄一∥I“+p 一点，以一dS一工，仆一p (4) 式中，W：广竺de，w为应变能密度。 ,Iba￡ 万H=0，即为： ^广8u(Vcr+功+L面(肌彩+工，抛+(n+O--m工广万¨(咒协加 +工．抛一(咒一仃一妒’)+工一万Pl(行一仃一∽=o (5) 式(5)左边前三项分别对应平衡方程、应力边界和间断面应力协调条件，我们关注的是第四、六 项，由于其必须符合强度条件，即：妒7=矿q和妒’=，l一吼。位移间断条件为 由于条件中没有对应力应变矩阵进行限定，因此上述条件可以是弹性筘阵、弹塑性矩阵，也 可以考虑应变硬化、应变软化情况。 当破坏强度吒随【H1或《五l变化(尤其是减小)的情况，拟静力平衡状态被打破，分析域出 现不平衡力，使静力平衡状态变化为动力平衡状态。由万FI=0，考虑到式(5)和分析体的运动 状态的改变，在某一时刻t，有 一JI。却’H一工，双∽H=一Ir砌·舀 (7) 当位移间断面考虑位移／速度弱化时，分析域将出现位移间断加速破坏现象。但对于界面为 位移／速度强化时，破坏进化不明显，除非加载继续增强。 当把分析域作为整体分析时，由关系矿7=矿呒、妒’=胛一叹、“++“一=o和I“+l+I“一l=《“l， 式(4)可以改写为 II-f、。P删y^。Ffi[dV—L≯邪+￡，+妒’删衄 (8) 式(8)的最后一项为位移间断界面一侧的作用力与另一侧表现出的滑移距离的积，显示了系统能 量的消耗。同时本文拟重点考察破坏强度吼随《“l变化引起的位移间断破坏过程，因此有： II=IV％SI腑y一工、町budV—L妒娜+￡，n一吼4“p (9a) A=工告∥2dy (9b) 式中，人为系统的动能。 考虑到分析城内不同时间阶段