弹塑性接触问题数值计算的二次规划算法研究.pdfVIP

第27卷专辑 固体力学学报 V01．27S．Issue 2006年12月 ACTAMECHANICASOLIDASINICA December2006 弹塑性接触问题数值计算的二次规划算法。 郭小明 (东南大学工卷力学系，南京，210096) 擅耍 针对弹塑性接触问题所推得的数值求解式子，运用二次规划法具体设计了算法．该算法采用有限的 基底交换运算就可得到收敛的数值解，具有较好的收敛性及较小的计算工作量．工程计算算倒结果表明文章所提 出的接触问题的求解方法是有效的，针对结构中接触问题所建立的效值计算模型能真实反映实际工作状况是可靠 的．文章中还详细给出了接触单元相关矩阵和向量的具体推导形式． 关t谓 接触问题．二次规划，有限元．算法 0 引育 口=1，2，…-m，卢=1，2} 接触状态方程的高度非线性，决定了接触问题 由文C3]的具体推导，可以得到解除了本构状 求解的难度，给工程应用带来了非常大的困难，如果 态约束方程及接触状态约束方程的弹塑性接触问题 再考虑其它非线性(如弹塑性接触问题)，在数值求 等价的变分不等式模式 解时会碰到许多困难II,S】． 求{口．A，页)∈K，使得 运用数学规划法是求解非线性问题的一种很有 效的工具，将接触问题考虑为有约束的最小化问题， 用序列线性规划或二次规划技术，则可认为是继经 c(r—A，．：I)一d(u，r—A)+．f(，．一A)≥ 典解析方法和迭代法之后的又一类方法．本文针对 L(口一口)， V{口，r，r)∈K (1) 弹塑性接触问题所建立的变分模式，通过有限元离 其中4，b，c，d．J，三，占，；，d一，厂的详细定义见文[3]． 散后得到了一组数值求解的式子，运用二次规划理 从所定义的解空间看，模型的实质就是：在所有 论，具体设计了数值计算的二次规划算法，该算法采 满足几何方程及位移边界条件的可能位移场中，真 用有限的基底交换运算就可得到收敛的数值解，具 实解{Ⅱ，A，夏)满足式(1)，该式子除了包含了原问题 有较好的收敛性及较小的计算工作量． 的平衡条件、协调条件和边界条件外，有效地解除了 1 弹塑性接触问题的变分不等式模式D1 弹塑性本构状态方程和接触状态方程的约束． 对于弹塑性接触问题，确定系统状态的物理量 2 有限元离散形式 有位移嘶，应变c#．应力％，接触位移昌，接触力n 在对式(1)进行有限元离散前，先引入相对位 等变量，由于其中的一些变量可通过定解方程从另 移形函数的概念【3]l 一些量导出，因此取位移嘶作状态变量；塑性流动 设接触单元由两个单元节点组成(如图1)，局 参数和滑动比例因子A，j作为控制变量进行讨论． 部坐标系下单元位移为 定义空间 矿={磷n，菇u，距’，砰’)T (2) Hl(n)为Sobolev空间，L2(n)为Hilbert空间 相对位移三可以用节点位移插值 H}(n)=(口I口∈H1(D)，“lt；矿) ；=丽 (3) 对为相对位移形函数． 网(n)一{翟l瞿∈Hl(D)，越lL；o} 日i(n)=[Hi(n)]’．霄2(D)=[H2(n)]。因此可得到局部接