结构塑性极限分析上限法数值计算方法研究

论文摘要

在土木工程中,结构物的极限承载力和破坏模式的确定是一项重要的研究题和工程问题。分析此类问题的方法大致分为两类:一类是弹塑性的增量分;另一类则是塑性极限分析方法。极限分析的上限和下限方法以塑性极限定为理论基础,是工程结构的设计和分析中直接而又严格的极限状态分析方法。工程的实际应用中通常采用极限分析的数值方法。其中,在工程结构的极限析中较为常用。另外,极限分析方法最终需要求解一个数学规划问题,根据体的情况大致可分为线性和非线性规划问题。随着问题维数的增加,数学规问题将可能成为大规模优化问题,其求解成为了一个难题。因此本文从经典性极限分析理论出发,进一步改进运动许可速度场的构造方法,并将数值优领域中提出的新算法应用于数值极限分析上限的数学规划问题的求解中,取的主要成果如下。在刚体有限元上限分析中,如果将安全系数定义为目标函数,则数学规划题就成为了带有约束的非线性规划问题。本文首次采用一种新型的优化算法P-free方法求解此非线性规划问题。求解非线性规划的常用算法为序列二次规（SQP）方法。然而,在初始点任意的情况下,传统的SQP在求解刚体有限上限分析法的非线性规划模型时出现了子问题不相容的问题而导致得不到最解,且在每个迭代步中都要花费大量计算来求解一个二次规划（QP）问题。对这一非线性规划问题,文中采用了一种新型非线性优化算法——QP-free法来求解刚体有限元上限分析法的非线性规划问题。该方法的转轴操作可以免子问题不相容的问题,并且在每个迭代步中将求解QP问题转化为求解三具有相同系数矩阵的线性方程组。而根据虚功率方程将安全系数表示为运动可速度场的函数,目的就是使得非线性规划问题的目标函数避免了其导数成常数向量,便于采用QP-free算法进行求解。通过两类算法对经典边坡稳定题的对比分析,QP-free算法比传统的SQP算法则更为有效。在上述的刚体有限元上限分析法中,数学规划模型的非线性是由于采用安系数作为评价指标而引起的。为了避免求解非线性规划问题,文中采用了加系数的定义来衡量结构的极限载荷,并且根据临界加速度的概念可由加载系计算安全系数,实现了利用线性规划模型来求解结构的安全系数。在此基础,通过在刚体单元之间的界面上的两个端点施加相关联流动法则,使得在构造运动许可速度场时可以考虑刚体单元形心速度的转动分量,因而得到了考虑转动破坏模式的刚体有限元上限法的线性规划模型。进而将该方法应用于边坡稳定分析和地基承载力计算中。根据边坡稳定和地基承载力的经典问题的分析验证了本章中所提出的刚体有限元上限法及其线性规划模型的正确性。另外,通过单纯形方法和原对偶内点算法对线性规划问题求解的对比发现,后者比前者的计算效率要高,尤其是针对大规模问题的情况。然而,对于大多数的岩土材料而言,相关联流动法则并不能反映其真实的力学特性。需要在塑性极限分析中考虑非相关联流动法则。因此,在考虑相关联流动法则的刚体有限元上限法中,应当将该方法推广到考虑非相关联流动法则的情况。在构造考虑转动模式运动许可速度场时,在刚体单元之间界面的两端点处施加非相关联流动法则,使得刚体有限元上限法可以在非相关联流动法则的情况下建立求解结构极限承载力的线性规划模型。并将考虑非相关联流动法则的刚体有限元上限法应用于膨胀系数对边坡稳定安全系数的影响研究中。另外,还可以将这一方法推广到研究土体剪胀性对土钉极限抗拔力的影响。为了克服调整网格划分带了的困难,在径向基点插值法构造运动许可速度场的基础上,并结合Cartesian积分变换数值积分法对能量耗散功率的近似计算,推导了无需任何网格的数值极限分析上限法的非线性数学规划模型。为了求解此类无网格上限法,采用了一种逐步区分刚性和塑性区域的直接迭代算法。针对摩擦型材料而言,无论强度参数取何值时,直接迭代算法的迭代控制因子都等于零,因此该方法不能应用于求解服从Mohr-Coulomb或Drucker-Prager破坏准则的摩擦型材料的极限承载力问题,而只能求解服从椭圆型屈服函数的非模型材料的极限承载力的问题。通过对垂直边坡的临界高度、带孔板的极限承载力以及厚壁圆筒的极限扩张压力的计算和分析,验证了本章所提出的无网格上限法及其无搜索迭代算法。

论文目录

• 摘要
• ABSTRACT
• 第1章 绪论
• 1.1 课题背景及研究的目的和意义
• 1.1.1 课题背景
• 1.1.2 课题的目的和意义
• 1.2 极限分析方法及其相关理论的研究现状
• 1.2.1 数值极限分析上限法计算理论的研究现状
• 1.2.2 极限分析方法的应用现状
• 1.3 本文主要研究内容
• 第2章 刚体有限元上限分析法的非线性规划模型及QP-FREE 算法求解
• 2.1 引言
• 2.2 刚体有限元上限分析法的非线性规划模型的改进
• 2.2.1 基于刚体有限元的运动许可速度场的构造
• 2.2.2 改进的非线性规划模型
• 2.3 QP-free 算法及其应用
• 2.3.1 QP-free 算法的具体步骤
• 2.3.2 垂直边坡临界高度的计算
• 2.3.3 承受垂直载荷的二维边坡稳定安全系数计算
• 2.4 本章小结
• 第3章 刚体有限元上限分析法的大规模线性规划模型及其应用
• 3.1 引言
• 3.2 考虑转动模式的刚体有限元上限法的线性规划模型
• 3.2.1 加载系数定义下极限分析上限法的一般形式
• 3.2.2 相关联流动法则情况下刚体有限元上限法的线性规划模型
• 3.2.3 非相关联流动法则情况下刚体有限元上限法的线性规划模型
• 3.3 岩土边坡稳定性分析
• 3.3.1 垂直边坡的临界高度问题
• 3.3.2 二维均质边坡的稳定性分析
• 3.3.3 二维含有软弱夹层的非均质边坡的稳定性分析
• 3.3.4 非相关流动法则对边坡稳定安全系数的影响分析
• 3.4 条形基础的极限承载力计算
• 3.4.1 不考虑土体自重的情况
• 3.4.2 考虑土体自重的情况
• 3.5 土钉极限抗拔力计算
• 3.5.1 土钉抗拔机理及其抗拔力计算的讨论
• 3.5.2 土钉极限抗拔力上限解
• 3.6 本章小结
• 第4章 无网格上限法的非线性规划模型及其直接迭代算法
• 4.1 引言
• 4.2 无网格上限法的非线性规划模型
• 4.2.1 移动最小二乘近似和径向基点插值法的对比
• 4.2.2 基于径向基点插值法的运动许可速度场的构造
• 4.2.3 基于Cartesian 积分变换的数值积分法
• 4.2.4 无网格上限法的非线性规划模型
• 4.3 无网格上限法非线性规划模型的直接迭代算法
• 4.4 垂直边坡临界高度的计算
• 4.5 板的极限承载力计算
• 4.5.1 von Mises 屈服准则下板的极限承载力分析
• 4.5.2 MDE 屈服准则下板的极限承载力分析
• 4.6 厚壁圆筒的极限扩张压力的计算
• 4.6.1 厚壁圆筒塑性极限内压力计算
• 4.6.2 组合厚壁圆筒塑性极限内压力计算
• 4.7 本章小结
• 结论
• 参考文献
• 附录
• 附录1 刚体有限元上限法线性规划模型建立的源程序
• 附录2 刚体有限元上限法的非线性规划问题的QP-free 算法的源程序
• 附录3 无网格上限法的非线性规划问题的建立及其直接迭代算法的源程序
• 攻读博士学位期间发表的论文
• 致谢
• 个人简历