共查询到20条相似文献，搜索用时 31 毫秒
1.
结合边界元法和离散元法,提出一种可以进行计算颗粒内部应力和破碎路径的方法。该方法利用离散元法求解颗粒的相互作用和每个颗粒上的荷载。然后利用边界元法计算颗粒的应力分布,为了实现动态平衡,将颗粒的加速度视为恒定大小的体力。但体力导致边界积分方程中出现域积分,故采用直线积分法将域积分转化为边界积分,以保证边界元法降维的优势。为了提高边界元的计算效率,对于几何形状相似的颗粒,以其中一个颗粒作为模板颗粒,只需要计算模板颗粒在局部坐标系中的系数矩阵,其他相似颗粒可以通过局部和全局坐标系之间的映射获得。在得到应力后,基于Hoek-Brown准则来判断颗粒是否破碎。此外,将破坏路径简化为直线,并采用最小二乘法拟合得到破坏路径。 相似文献
2.
介绍了联合运用积分方程法与有限元法(简称混合法)来计算三维非均匀介质中电磁场分布的理论方法,并进行了数值模拟实验。混合法的原理是引入一个包围非均匀目标体的虚构边界,在边界内部的场用有限元法模拟,在边界外部(包括边界)的场用积分方程表达,二者在边界上通过场的连续性耦合起来。数值实验结果表明,混合法既能显著地减小网格规模,又能灵活地模拟复杂的介质情况,且计算精度较高。 相似文献
3.
基于边界元法的边坡矢量和稳定分析 总被引：4，自引：0，他引：4
矢量和法物理力学意义明确,计算简单,且能根据边坡当前的应力分布状态合理地评价其整体稳定性状态。其中边坡的应力状态通常是采用有限元法来求解。由于边界元法具有研究问题降阶、离散化带来的误差值仅产生在边界以及计算量小等优点,在工程中得到了广泛应用;对于平面问题,以源点作为原点,以所积分单元的切向和法向为坐标轴建立局部坐标系,对于线性单元可以得到所有积分的解析解。因此,可以得到计算区域内部任意点的场变量的解析解,这就保证了位于边界附近区域场变量的精度。利用边界元法得到二维边坡体内连续的应力分布状态,使用矢量和法对该边坡进行稳定性分析,并且与基于有限元的矢量和法、极限平衡法进行对比分析。边坡圆弧滑面和折线滑面的计算结果表明,基于边界元法得到的矢量和安全系数和基于有限元的矢量和法、极限平衡法基本一致;边界元法对应的矢量和安全系数对边界单元尺寸不敏感。 相似文献
4.
碎石桩复合地基的作用机理分析及沉降计算 总被引：24，自引：2，他引：22
通过对规范方法的讨论及对碎石桩复合地基桩土竖向及径向应力应变关系分析,提出桩土变形模式,建立桩土应力、应力比表达式及沉降计算方法。 相似文献
5.
采用自制的压–拉转换装置,配合RMT 150C岩石力学试验系统及数字散斑相关方法,对双边非对称裂纹类岩石平板试样进行直接拉伸试验,得到类岩石试样的拉应力–应变曲线、试样表面应变场演化过程和裂纹扩展模式。研究发现,类岩石试样直接拉伸的拉应力–应变曲线大致可以分为4个阶段：（1）近似线性阶段,预制裂纹基本不起作用,应力随应变增加较快,试样表面应变场的分布主要受试样内部的孔隙及颗粒的影响;（2）整体缓慢增加阶段,两预制裂纹和试样内部的孔隙及颗粒共同影响试样表面应变场的分布,整体上应力随应变呈增加趋势;（3）短暂峰值过渡阶段,试样中某个预制裂纹对试样表面应变场的分布起决定性作用;（4）破坏阶段,裂纹起裂位置在应变场相对集中区域,并扩展导致试样破坏。对于直接拉伸条件下的双边非对称裂纹平板试样,其中某条预制裂纹会率先扩展,先向远离前方裂纹的方向扩展,再向靠近前方裂纹的方向扩展,对采用数值模拟方法研究张拉应力状态下裂纹相互作用扩展规律具有重要意义。 相似文献
6.
在平面应变压缩条件下,采用拉格朗日元法研究了材料缺陷对岩样应变局部化及宏观力学行为的影响。在数值计算中,采用了莫尔-库仑与拉破坏复合的破坏准则,峰后岩石的本构关系为线性应变软化。对于理想岩样（不含任何缺陷）,从剪切应变率的等值线图及变形网格图发现,变形场相对试样垂直对称轴对称。在试样的边界上设置材料缺陷后,变形场的对称性被打破。材料缺陷附近是局部化现象启动的主要位置。与不含材料缺陷试样相比,含缺陷试样的局部化提前启动。当缺陷位于试样左边界中部附近时,试样内部出现多条剪切带,发生韧性剪切破坏,峰后应力-轴向应变曲线和应力-侧向应变曲线均倾向于韧性,岩样的稳定性增强;当缺陷位于试样上、下端面附近时,仅出现一条贯穿试样左右边界的剪切带,发生脆性剪切破坏,峰后两种曲线倾向于脆性,易发生失稳破坏。从数值结果中还观测到了剪切带的跳跃或迁移现象及剪切带的相互竞争、此消彼涨的现象。当岩样中不包含任何缺陷时,试样应力-变形曲线的峰值强度最高。随着缺陷的上移,峰值强度下降,直到峰值强度基本保持不变。缺陷越接近于岩样的下端,对岩样应力-变形曲线的峰值强度的影响越小。 相似文献
7.
8.
数值模拟探讨青藏高原东部应力应变场及断裂活动性 总被引：1，自引：0，他引：1
青藏高原东部是我国活动构造非常强烈的区域,现今的应力应变状态对活动构造研究有重要参考意义。本文依据地质资料和岩石力学性质构建相应的有限元数值模型,以GPS观测为边界约束,计算得到了青藏高原东部应力应变场,在此基础上探讨了主要断裂的活动特征。结果显示,距离印度板块和欧亚大陆板块碰撞带较近的青藏高原,汇聚速度、最大压应变率和最大压应力均较高,而随着远离碰撞带各个数值相应降低;在青藏高原东部,以喜马拉雅东构造结为旋转中心的一系列弧形断裂滑动速率相对较高,如玉树—鲜水河—小江断裂带,巴塘—金沙江断裂带,嘉黎断裂带等,这些走滑断裂带控制着青藏高原东南缘物质旋转挤出的特征,而其它断裂带滑动速率则相对较低,在块体内部则更低;结果同时表明应变速率分配高的断裂达到屈服的时间短,对应的地震复发周期短。 相似文献
9.
本文讨论了在均匀谐变场激发下导电导磁体电磁响应的边界单元法计算问题。文中首先建立矢量位的边值问题,然后导出其边界积分方程,最后用边界元法进行求解。本方法的计算结果具有精度高、速度快的优点。 相似文献
10.
本程序采用非线性变模量模型,并应用了岩石的无侧向应变加卸载试验曲线和动抗拉强度。计算中用应力波传播的加卸载判别准则来决定单元的加卸载状态,对受拉破坏单元进行内部应力调整和应力转移。使用增量法求得每个时刻的位移、速度、加速度和应力、应变。结果表明计算的和实测的应力应变比较接近。 相似文献
11.
本文以二层介质中存在球体的模型为例，讨论了用边界元方法解决层状介质中存在三维不均匀体情况下点源场视电阻率的计算问题。文中通过对基本解的选取，避免了对水平界面的积分计算，并且在积分方程中消去了位的法向导数，使得积分方程的形式更为简单，也大大地减少了所需的边界剖分节点数和未知量的个数，节省了计算机内存，提高了计算精度。计算的结果表明，边界元法对本问题的研究是成功的。 相似文献
12.
构造应力场边界载荷反演的有限元逆逼近法 总被引：2，自引：0，他引：2
构造应力场问题的边界条件、加载方式是地应力场分析的难题,由于构造运动未知和地质构造复杂,无法直接求解。边界力反演是进行构造应力场研究的重要手段。从弹性力学基本方程出发,针对油气储层地应力场反问题,结合对区域构造应力场认识,利用阻尼最小二乘法建立了应力场反演的约束模型,提出了一种用于反演构造应力场边界力的反分析方法,该方法根据部分实测点的应力值进行边界力反演,并与有限元方法和有限差分法相结合,以边界力为参量逼近目标函数。根据研究区少量地应力实测资料采用该方法进行最优化计算,使得计算应力场与实测应力达到最优拟合,以弥补目前的回归反演方法和边界载荷调整法的不足来提高计算精度。算例表明了该方法的有效性,可用于解决构造应力场边界力的反演问题。 相似文献
13.
鄂尔多斯盆地奥陶纪“L”状边缘隆起演化过程及其构造背景 总被引：1，自引：6，他引：1
将鄂尔多斯盆地简化为近似直角梯形的等厚各向同性弹性薄板模型进行应力—应变场数值模拟,力图揭示鄂尔多斯地区奥陶纪发育的“L”状边缘隆起的形成背景及演化过程。模拟中Z轴方向应变(zε)正值区对应于盆地内部的隆起区,故zε正等值线形态对应于隆起形态。设定的多种不同的边界条件中,只有在南、西边界分别施以大小相同的边界应力时,εz正等值线形态才出现与实际地质情况相类似的“L”状,该边界条件下的模拟再现了盆地奥陶纪“L”状隆起的形成过程:在盆地南、西边界首先产生两个月牙状隆起,随着应力传递,南、西边界处的两个隆起扩展而连接,形成一个“L”状隆起,随后该隆起不断向盆地内部延伸。模拟结果证实由于奥陶纪鄂尔多斯盆地南缘秦岭海槽的向北俯冲及西缘贺兰海槽和秦祁海槽的东西向扩张,鄂尔多斯盆地南、西边界同时受到挤压作用,正是因为两边界所受挤压力大小相似,“L”状隆起才于盆地南西边界发育。 相似文献
14.
大型地下洞室群围岩应力-损伤-渗流耦合分析 总被引：5，自引：2，他引：3
以渗透体力来考虑渗流场的力学效应,建立了应力-损伤-渗透系数关系方程来考虑应力和损伤对渗流场的影响,结合岩体结构的三维弹塑性损伤有限元分析,建立了大型地下洞室开挖围岩应力-损伤-渗流耦合的计算模型。该计算模型求解的难度主要体现在岩体材料弹塑性、损伤、渗流自由面边界、渗流溢出边界、应力-损伤-渗流相互影响关系等。提出分步迭代法对以上因素进行归纳后按一定顺序分别进行迭代求解,取得了良好的计算效果。将该方法应用于某水电站大型地下洞室群分析,得出了一系列有意义的结论。 相似文献
15.
建立尾砂的本构模型是开展尾矿坝数值模拟和安全评价的重要基础,而目前尾砂的本构模型研究多集中于非线性弹性模型,如Duncan-Chang模型,关于弹塑性模型的研究较少。结合尾砂的应力-应变特性,提出了一种适用于描述尾砂力学特性的改进广义塑性模型。基于用户自定义材料子程序UMAT,将提出的模型在ABAQUS中二次开发实现,应力积分采用Runge-Kutta显式积分。通过三轴试验模拟验证,偏差应力曲线表明,有限元计算结果可以反映围压对应力-应变关系的影响,抗剪强度随应变逐渐硬化,达到峰值强度,随后发生应变软化,抗剪强度有所下降。体应变曲线结果表明,广义塑性模型可以很好地描述体应变曲线的体缩-体胀发展过程,与试验曲线吻合较好。同时有限元数值模拟结果和理论值误差很小,且和试验结果拟合较好。该研究成果可进一步用于尾矿坝的应力变形计算和安全评价。 相似文献
16.
利用格林公式及点源场的边界条件,将地下半空间位函数的体积积分变为地面位函教的积分方程以后,用边界元方法解此积分方程便可得到地面各节点的电位值。其计算精度和速度取决于剖分方式和节点密度。本文采用变网格的剖分方式,用较少的点实现了计算,并兼顾了精度和速度。 相似文献
17.
鄂尔多斯盆地石炭纪中央古隆起形成机制 总被引：1，自引：4，他引：1
将鄂尔多斯盆地简化为受南北挤压的等厚各向同性弹性薄板模型进行应力—应变场数值模拟,力图揭示鄂尔多斯地区石炭纪出现的细腰状中央古隆起的形成机制。模拟中Z轴方向应变(εz)正值区对应于盆地内部的隆起区,zε正等值线形态对应于隆起形态;单轴挤压条件下,εz正等值线总会呈现沿应力轴方向延伸的细腰状形态;点作用力产生的zε正等值线范围局限,而线作用力产生的εz正值区分布较广。模拟结果表明鄂尔多斯盆地石炭纪细腰状隆起是在南北边界受挤压条件下,应力—应变在盆地内部传递过程中所必然出现的结果,南北向点作用力比南北向线作用力产生的zε等值线形态更接近于鄂尔多斯盆地中央古隆起形态;点作用力可能代表了石炭纪微板块间的点碰撞或者弧—陆碰撞。 相似文献
18.
断裂问题的扩展有限元法研究 总被引：3，自引：0，他引：3
扩展有限元（extended finite element method,XFEM）是近年来发展起来的、在常规有限元框架内求解不连续问题的有效数值计算方法,其基于单位分解的思想,在常规有限元位移模式中加入能够反映裂纹面不连续性的跳跃函数及裂尖渐进位移场函数,避免了采用常规有限元计算断裂问题时需要对裂纹尖端重新加密网格造成的不便。在推导扩展有限元算法的基础上,分析了应力强度因子的J积分计算方法及积分区域的选取。采用XFEM对I型裂纹进行了计算,有限元网格独立于裂纹面,无需在裂纹尖端加密网格;分析了积分区域、网格密度对应力强度因子计算精度的影响,指出了计算应力强度因子的合适参数,验证了此方法的可靠性和准确性。 相似文献
19.
巷道开挖围岩能量释放与偏应力应变能生成的分析计算 总被引：1，自引：0，他引：1
巷道开挖,围岩能量释放,同时在围岩中产生偏应力。围岩应力是原岩应力与偏应力的叠加,偏应力或偏应力能控制岩体破坏。在静水压力 和岩体体积应变为0的条件下,利用文[1]在弹性、非线性软化本构模型导得的巷道围岩应力分布表达式,用重积分计算了围岩弹性区和软化区中的偏应力应变能 ,证明了 可以简捷地用地应力 关于巷壁位移 做一次积分再乘巷壁周长的途径来得到,阐述了该计算途径的原理。巷道开挖过程围岩释放的能量等于围岩压力 关于 的积分乘巷壁周长。由此,可通过 ~ 曲线、 ~ 曲线所围面积的几何形式,表示围岩偏应力能、围岩弹性能释放量随 变化的情况。所得研究结果可以深化围岩由于开挖产生的力学响应及挖成后围岩工况规律的认识。 相似文献
20.
青海玉树M_S7.1地震发震过程的数值模拟 总被引：1，自引：1，他引：0
根据玉树地区的地应力场、速度场和断层展布,对青海玉树2010年4月14日MS7.1级地震发震机理进行了数值模拟。将围岩看成弹性体,断层看成具有应变软化的弹塑性体,断层和围岩组成统一的地质介质系统。在地应力、孔隙压力及边界位移的作用下,应力逐渐积累,当达到断层摩擦破坏强度时,断层产生应变软化,断层突然滑动,能量突然释放,应力突然下降,形成地震。模拟结果表明:玉树7.1级地震是在印度板块向北推挤,青藏高原向东南侧向挤压,在玉树地区形成主压应力为北东80°方向的水平应力场,使甘孜-玉树断裂带产生左旋走滑错动形成的。计算结果给出了应力降、能量释放量、断层走滑错动量、地震复发周期、应力积累速度等重要参数,模拟结果与野外调查资料具有较好的一致性。 相似文献