首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
林姗  李春光  孙冠华  王水林  杨永涛 《岩土力学》2018,39(10):3863-3874
接触问题是一种常见的非线性问题,如何能够很好地模拟接触面的变形及受力特性,以及实现对变形体间的接触问题的真实模拟是该领域研究的难点问题。基于二维接触问题的实际物理意义,分别在法向和切向建立等价的互补模型。用非线性互补函数(NCP)中的Fischer-Burmeister(FB)函数将互补函数模型转化为非光滑方程组表达,用常规的Newton法求解。同时,基于高斯积分法可以用较少的积分点达到较高的精度,为了进一步提高求解精度,改善不连续的通病,对面-面接触模型在高斯点上对接触面进行处理,可通过调节积分点数目对求解精度进行控制,方法易于理解,实现方便。在此基础上建立二维接触有限元模型,通过一系列工程算例验证该方法的可行性与有效性。结果表明,与ABAQUS有限元的计算结果相比,该方法有着较高的精度,更真实地反映问题的实际。  相似文献   

2.
马文涛  师俊平  李宁 《岩土力学》2012,33(10):3145-3150
针对摩擦接触问题,给出一种新型的无网格数值方法。该方法基于单位分解思想,在标准无网格Galerkin法的位移模式中嵌入不连续函数和裂尖奇异函数,分别反映接触面的不连续性和接触面端点的奇异性;结合接触摩擦定律,继而构造出求解摩擦接触问题的无网格线性互补模型。在该方法中,接触面方程使用接触点对离散,全局离散系统方程则转化为标准的线性互补问题,可以很方便地使用Lemke算法求解。算例分析证明,本文方法的正确性和有效性。  相似文献   

3.
成层地基轴对称问题的基本解   总被引:1,自引:2,他引:1  
基于各向同性弹性体基本方程的转换及Hankel变换,得到了轴对称荷载作用下,无限弹性层轴向不同深度经Hankel变换的位移应力向量间的传递矩阵,运用该传递矩阵可求解成层的、层间完全接触情形和各向同性空间轴对称问题。  相似文献   

4.
模拟三维裂纹问题的扩展有限元法   总被引:4,自引:1,他引:3  
余天堂 《岩土力学》2010,31(10):3280-3285
扩展有限元法是一种在常规有限元框架内求解强和弱不连续问题的新型数值方法,其计算网格与不连续面相互独立,因此模拟移动不连续面时无需对网格进行重新剖分。给出了模拟三维裂纹问题的扩展有限元法。在常规有限元位移模式中,基于单位分解的思想加进一个阶跃函数和二维渐近裂尖位移场,反映裂纹处位移的不连续性。用两个水平集函数表示裂纹。采用线性互补法求解裂纹面非线性接触条件,不需要迭代,提高了计算效率。采用两点位移外推法计算裂纹前缘应力强度因子。给出了3个三维弹性静力问题算例,其结果显示了所提方法能获得高精度的应力强度因子,并能有效地处理裂纹面间的接触问题,同时表明扩展有限元结合线性互补法求解不连续问题具有较好的前景。  相似文献   

5.
姜育松  苏超 《岩土力学》2011,32(11):3497-3502
给出了一种求解多块体接触的有限元解法,即将接触面条件精确引入并以节点接触力为基本未知量,把高度接触非线性问题凝缩在可能接触面上进行。首先从两块体三维接触力学模型着手,从整体有限元方程导出增量形式的协调方程,在此基础上详细讨论了多块体接触协调方程的建立方法。随后以广义Mohr-Coulomb准则为摩擦力条件,给出了增量形式的定解和判定条件。程序实施时,针对三维接触问题滑动方向无穷的问题,同时对接触状态及滑动摩擦力分量进行判定,使得迭代收敛快速准确;对于多块体常见的非光滑接触,即角点问题,将角点分类,并根据角点实际的接触状态分情况给出了更为合理的处理方法。通过算例验证后的程序应用于不连续面地下洞室开挖模拟,取得了满意的结果。  相似文献   

6.
对对流占优的三维溶质运移问题提出了分步广义迎风解法,首先利用N.N.Ya-nenko对水动力弥散方程分步求解的思想,将原来的一个定解问题分解为两个定解问题即对流定解问题和扩散定解问题,对对流定解问题采用广义迎风对偶单元均衡法求解,对扩散定解问题采用一般的Galerkin有限元法求解,不仅避免了用一般有限元法和有限差分法求解对流占优的地下水水质数学模型时常出现数值弥散和过量问题,而且避免了求节点速度这一步,节省运算步骤,对井点的浓度变化给出了更合适的求解方法。  相似文献   

7.
刘俊卿  曹书文  孙莹  王炜 《岩土力学》2006,27(Z1):993-996
求解了无拉力双参数地基上厚、薄板受任意移动荷载作用下的动力响应问题和板的脱离问题。从板-地基接触系统的分区广义势能泛函出发,通过建立和分析板单元区、地基单元区和连接板-地基的接触单元区的势能泛函,给出求解无拉力双参数地基板受任意移动荷载作用下的动力响应问题的分区广义协调元法。利用纽马克(Newmark)法经MATLAB编程叠代数值计算,验证了文中方法分析无拉力双参数地基板的动力响应问题是有效的。  相似文献   

8.
赵强  焦玉勇  张秀丽  谢壁婷  王龙  黄刚海 《岩土力学》2019,40(11):4515-4522
非连续变形分析方法(DDA)是一种平行于有限元法的新型数值计算方法,该方法基于最小势能原理,把每个离散块体的变形、运动和块体之间的接触统一到平衡方程中进行隐式求解。然而,传统DDA方法在计算过程中需组装整体刚度矩阵并联立求解方程组,在用于大型岩土工程问题的三维数值模拟时占用内存较大、耗时较长、计算效率极低。因此,提出一种基于显式时间积分的三维球颗粒DDA方法。该方法在求解过程中不需要组装整体刚度矩阵,在求解加速度时,由于质量矩阵为对角矩阵,可存储为一维向量占用内存较少,且可分块逐自由度求解,效率较高,在接触判断上采用最大位移准则简化了接触算法,采用较小的时步,保证了计算的精确性;通过几个典型算例验证了该方法的准确性及计算效率。  相似文献   

9.
三维溶质运移问题的分步广义迎风解法   总被引:1,自引:0,他引:1  
对对流占优的三维溶质运移问题提出了分步广义的迎风解法,首先利用N,N,Yanenko对水动力弥散方程分步求解的思想,将原来的一个定解问题分解为两个定解问题即对流定解问题和扩散定解问题,对对流定解问题采用广义迎风对偶单元均衡法求解,对扩散定解问题采用一般的Galerkin有限元法求解,不仅避免了用一般有限元法和有限差分法求解对流占优的地下水水质数学模型时常出现数值弥散和过量问题,而且避免了求节眯速度  相似文献   

10.
毛坚强 《岩土力学》2004,25(10):1592-1598
大多数的岩土工程问题都会涉及到岩土体-结构之间的相互作用,它们之间界面的合理模拟对计算结果有重要的影响。目前常用的Goodman节理单元、Desai薄层单元等虽然简单好用,但也存在着明显的缺陷。相比之下,作为接触问题计算,可更为合理地模拟接触面的力学行为。针对变形体-刚体这种特殊的接触问题,应用变形体的虚功原理建立了相应的虚功方程,在此基础上得出有限元计算公式,并建立了弹塑性接触问题有限元计算的实现方法。通过对刚性条形基础基底压力分布和刚性挡墙土压力分布等算例的计算,验证了该法的正确性和求解工程问题的有效性。  相似文献   

11.
叶照桂 《探矿工程》2004,31(3):29-30
通过岩庄水库涵洞顶管施工的接触灌浆事例,阐述了有关接触灌浆工艺的几个基本问题,对讨论的问题进行了重点分析。  相似文献   

12.
基于遗传算法最短路径问题的探讨   总被引:9,自引:0,他引:9  
对用遗传算法求解最短路径问题作了有益的尝试,详细分析了求解最佳路径的遗传算法的构成要素,提出了一种新的交叉算子,并且论证了算法参数对结果的影响。通过仿真实验,给出了算法的主要性能参数,证明了算法的可行性。并指出了遗传算法求解最短路径问题的不足之处。  相似文献   

13.
郑宏  谭国焕  刘德富 《岩土力学》2006,27(11):1880-1884
所谓近似对称矩阵是指矩阵中仅有极少一部分元素是非对称的,在将对角线以上的非对称元素用其对角线以下的对称元素替代后,该矩阵就变成了一个对称矩阵。在求解非线性连续介质力学问题时常常会遇到近似非对称矩阵。基于Sherman-Morrison定理,给出了一种新的近似非对称矩阵的分解算法。在确保数值稳定性的前提下,无论在求解效率还是在内存开销方面新算法都优于一般的非对称稀疏矩阵的求解器,且仅需对传统的基于LDLT分解的求解器略做修改,即可开发出适应于对称和非对称稀疏矩阵的求解器。最后用一个摩擦接触算例,显示了新算法的优越性。  相似文献   

14.
热传导有限元的神经计算原理及应用   总被引:1,自引:0,他引:1  
根据热传导有限元分析方法的特点,把有限元计算问题转化为带约束的非线性优化问题,并给出求解该问题的改进Hopfield(TH)网络,最后对一个简单温度场神经计算进行数值仿真,仿真结果表明,神经网络能完成有限元模型的求解。  相似文献   

15.
土体固结弹塑性分析的参数二次规划理论及有限元解   总被引:1,自引:0,他引:1  
基于广义Biot理论对土体弹塑性固结过程进行求解,建立了问题对应的参变量变分理论,并给出了数学证明,对此基础上推导了有限元分析列式,问题的求解最终化为参数二次规划问题,本文提出的方法适用于固结弹塑性分析的关联与非关流流动问题,也可处理各类软化问题。  相似文献   

16.
根据水质模型的具体特点,对不同的方程采用不同方法,水流问题用有限元法;对流弥散方程先用算子分裂的方法分解为两个方程,即对流方程和弥散方程,前者用高精度广义迎风格式求解,对弥散方程则采用多单元均衡格式法求解,最后合成为高精度广义迎风均衡格式求出溶质浓度。通过对数值实验例子的计算和实验溶质迁移的模拟,可以看出在求解对流弥散定解问题时,广义迎风均衡格式克服了有限元数值波动和浓度出现负值的问题,与有限元相比有较大改进。  相似文献   

17.
三维稳定渗流的 p 型自适应有限元分析   总被引:3,自引:2,他引:3  
费文平  陈胜宏 《岩土力学》2004,25(2):211-215
建立了三维稳定渗流的p型自适应有限元分析的误差估计器。充分利用p型有限元前处理少、计算精度高和收敛速度快的优点,分析了p型有限元法求解渗流问题的技术特点,给出了求解三维稳定渗流的自适应升阶策略。通过具体算例,研究了单元阶次及网格尺度对计算结果的影响,验证了三维p型有限元法在求解渗流问题时的可行性。  相似文献   

18.
工程体的稳定性问题最终要求确定的解。由概率模型模拟到简单耦合及复杂耦合求解逼近结构面确定性研究表明,求解耦合问题,一方面要求遵守确定问题求解的精确性约束,另一方面要求满足随机性能问题的概率特征约束,同时又要考虑结构面分布的差异性。  相似文献   

19.
石露  李小春  王伟  白冰 《岩土力学》2011,32(12):3805-3811
在处理含有裂隙这类不连续问题时,常规有限元方法需要对裂隙尖端部位进行局部网格加密,当裂隙扩展时还需要进行网络重构。基于单位分解思想的扩展有限元成功解决了常规有限元难以处理的裂隙类不连续问题。在研究复合裂隙时,通常需要考虑裂隙的接触问题。基于互补理论,建立了裂隙面上相对位移和接触力的互补方程,并采用牛顿法求解,无需开闭迭代,且能够快速收敛。最后,对含裂隙平板进行受压数值试验,计算结果表明,基于互补理论的扩展有限元接触算法能够有效地阻止裂隙两端网格的相互嵌入,且获得裂隙面上的应力分布与实际一致。  相似文献   

20.
Laplace方程Cauchy问题的一种数值解法   总被引:2,自引:1,他引:2  
研究了Laplace方程Cauchy问题的数值求解,该问题是一个典型的病态问题。利用格林(Green)公式将Laplace方程的Cauchy问题转化为Hausdorff矩问题。本文利用一种新方法,即矩问题的积分方程方法,求解矩问题,设计了二维Laplace方程Cauchy问题稳定的算法,给出了近似解的误差估计,并对二维Cauchy问题进行了数值模拟。  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号