应力与变形的数值模型方法_数值模拟软件FLAC介绍_龚纪文

FLAC-3D(ThreeDimensionalFastLagrangianAnalysisofContinua)是美国ItascaConsultingGouplnc开发的三维快速拉格朗日分析程序,该程序能较好地模拟地质材料在达到强度极限或屈服极限时,发生的破坏或塑性流动的力学行为,特别适用于分析渐进破坏和失稳以及模拟大变形.FLAC3D分析的使用领域根据手册总结如下:(1)承受荷载能力与变形分析:用于边坡稳定和基础设计(2)渐进破坏与坍塌反演:用于硬岩采矿和隧道设计(3)断层构造的影响研究:用于采矿设计(4)施加于地质体锚索支护所提供的支护力研究:岩锚和土钉的设计(5)排水和不排水加载条件下全饱和流体流动和孔隙压力扩散研究:挡土墙结构的地下水流动,和土体固结研究(6)粘性材料的蠕变特性:用于碳酸钾盐矿设计(7)陡滑面地质结构的动态加载:用于地震工程和矿山岩爆研究(8)爆炸荷载和振动的动态响应:用于隧道开挖和采矿活动(9)结构的地震感应:用于土坝设计(10)由于温度诱发荷载所导致的变形和结构的不稳定(11)大变形材料分析:用于研究粮仓谷物流动和放矿的矿石流动10种材料本构模型Flac3D中为岩土工程问题的求解开发了特有的本构模型,总共包含了10种材料模型:(1)开挖模型null(2)3个弹性模型(各向同性,横观各向同性和正交各向同性弹性模型)(3)6个塑性模型(Drucker-Prager模型、Morh-Coulomb模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型和修正的cam粘土模型).Flac3D网格中的每个区域可以给以不同的材料模型,并且还允许指定材料参数的统计分布和变化梯度.还包含了节理单元,也称为界面单元,能够模拟两种或多种材料界面不同材料性质的间断特性.节理允许发生滑动或分离,因此可以用来模拟岩体中的断层、节理或摩擦边界.FLAC3D中的网格生成器gen,通过匹配、连接由网格生成器生成局部网格,能够方便地生成所需要的三维结构网格.还可以自动产生交岔结构网格(比如说相交的巷道),三维网格由整体坐标系x,y,z系统所确定,这就提供了比较灵活的产生和定义三维空间参数.五种计算模式(l)静力模式:这是FLAC-3D默认模式,通过动态松弛方法得静态解.(2)动力模式:用户可以直接输人加速度、速度或应力波作为系统的边界条件或初始条件,边界可以固定边界和自由边界.动力计算可以与渗流问题相藕合.(3)蠕变模式:有五种蠕变本构模型可供选择以模拟材料的应力-应变-时间关系:Maxwell模型、双指数模型、参考蠕变模型、粘塑性模型、脆盐模型.(4)渗流模式:可以模拟地下水流、孔隙压力耗散以及可变形孔隙介质与其间的粘性流体的耦合.渗流服从各向同性达西定律,流体和孔隙介质均被看作可变形体.考虑非稳定流,将稳定流看作是非稳定流的特例.边界条件可以是固定孔隙压力或恒定流,可以模拟水源或深井.渗流计算可以与静力、动力或温度计算耦合,也可以单独计算.(5)温度模式:可以模拟材料中的瞬态热传导以及温度应力.温度计算可以与静力、动力或渗流计算藕合,也可单独计算.模拟多种结构形式(l)对于通常的岩体、土体或其他材料实体,用八节点六面体单元模拟.(2)FIAC-3D包含有四种结构单元:梁单元、锚单元、桩单元、壳单元.可用来模拟岩土工程中的人工结构如支护、衬砌、锚索、岩栓、土工织物、摩擦桩、板桩等.(3)FLAC-3D的网格中可以有界面,这种界面将计算网格分割为若干部分,界面两边的网格可以分离,也可以发生滑动,因此,界面可以模拟节理、断层或虚拟的物理边界.有多种边界条件边界方位可以任意变化,边界条件可以是速度边界、应力边界,单元内部可以给定初始应力,节点可以给定初始位移、速度等,还可以给定地下水位以计算有效应力、所有给定量都可以具有空间梯度分布.FLAC-3D内嵌语言FISHFLAC-3D具有强大内嵌语言FISH,使得用户可以定义新的变量或函数,以适应用户的特殊需要,例如,利用HSH做以下事情:(l)用户可以自定义材料的空间分布规律,如非线性分布等.(2)用户可以定义变量,追踪其变化规律并绘图表示或打印输出.(3)用户可以自己设计FLAC-3D内部没有的单元形态.(4)在数值试验中可以进行伺服控制.(5)用户可以指定特殊的边界条件.(6)自动进行参数分析(7)利用FLAC-3D内部定义的Fish变量或函数,用户可以获得计算过程中节点、单元参数,如坐标、位移、速度、材料参数、应力、应变、不平衡力等.FLAC-3D前后处理功能FLAC-3D具有强大的自动三维网格生成器,内部定义了多种单元形态,用户还可以利用FISH自定义单元形态,通过组合基本单元,可以生成非常复杂的三维网格,比如交叉隧洞等.在计算过程中的任何时刻用户都可以用高分辨率的彩色或灰度图或数据文件输出结果,以对结果进行实时分析,图形可以表示网格、结构以及有关变量的等值线图、矢量图、曲线图等,可以给出计算域的任意截面上的变量图或等直线图,计算域可以旋转以从不同的角度观测计算结果.FLAC3D计算分析一般步骤与大多数程序采用数据输入方式不同,FLAC采用的是命令驱动方式.命令字控制着程序的运行.在必要时,尤其是绘图,还可以启动FLAc用户交互式图形界面.为了建立FLAC计算模型,必须进行以下三个方面的工作:(1)有限差分网格(2)本构特性与材料性质(3)边界条件与初始条件完成上述工作后,可以获得模型的初始平衡状态,也就是模拟开挖前的原岩应力状态.然后,进行工程开挖或改变边界条件来进行工程的响应分析,类似于FLAC的显式有限差分程序的问题求解.与传统的隐式求解程序不同,FLAC采用一种显式的时间步来求解代数方程.进行一系列计算步后达到问题的解.在FLAC中,达到问题所需的计算步能够通过程序或用户加以控制,但是,用户必须确定计算步是否已经达到问题的最终的解.后处理(一)用tecplot绘制曲线(1)第一主应力(2)xdisp、ydisp、zdisp、disp(二)用excel做曲线隧道(1)做地表沉降槽(zdisp)(2)地表横向位移(xdisp)(3)隧道中线竖向沉降曲线(zdisp)(4)提取位移矢量图,(5)显示初期支护结构内力(6)显示state(找塑性区)基坑(1)做地表沉降槽(zdisp)(2)提取位移矢量图,(3)显示初期支护结构内力(4)显示state(找塑性区)边坡(1)做安全系数和应变图模型最优化用FLAC3D解决问题时,为了得到最有效的分析使模型最优化是很重要的.(1)检查模型运行时间:一个FLAC3D例子的运行时间是区域数的4/3倍.这个规则适用于平衡条件下的弹性问题.对于塑性问题,运行时间会有点改变,但是不会很大,但是如果发生塑性流动,这个时间将会大的多.对一个具体模型检查自己机子的计算速度很重要.一个简单的方法就是运行基准测试.然后基于区域数的改变,用这个速度评估具体模型的计算速度.(2)影响运行时间的因素:FLAC3D有时会需要较长时间才可以收敛主要发生在下列情况下:(a)材料本身刚度变异或材料与结构及接触面之间的刚度差异很大.(b)划分的区域尺寸相差很大.这些尺寸差异越大编码就越无效.在做详细分析前应该研究刚度差异的影响.例如,一个荷载作用下的刚性板,可以用一系列顶点固定的网格代替,并施以等速度.(记住FIX命令确定速度,而不是位移.)地下水的出现将使体积模量发生明显的增加(流体-固体相互作用).(3)考虑网格划分的密度:FLAC3D使用常应变单元.如果应力/应变曲线倾斜度比较高,那么你将需要许多区域来代表多变的分区.通过运行划分密度不同的同一个问题来检查影响.FLAC3D应用常应变区域,因为当用多的少节点单元与用比较少的多节点单元模拟塑性流动时相比更准确.应尽可能保持网格,尤其是重要区域网格的统一.避免长细比大于5:1的细长单元,并避免单元尺寸跳跃式变化(即应使用平滑的网格).应用GENERATE命令中的比率关键词,使细划分区域平滑过渡到粗划分区域.(4)自动发现平衡状态:默认情况下,当执行SOLVE命令时,系统将自动发现力的平衡.当模型中所有网格顶点中所有力的平均量级与其中最大的不平衡力的量级的比率小于1*10时,认为达到了平衡状态.注意一个网格顶点的力由内力(例如,由于重力)和外力(例如,由于所加的应力边界条件)共同引起.因为比率是没有尺寸的,所以对于有不同的单元体系的模型,在大多数情况下,不平衡力和所加力比率的限制给静力平衡提供了一个精确的限制.同时还提供了其他的比率限制;可以用SETratio命令施加.如果默认的比率限制不能为静力平衡提供一个足够精确的限制,那么应考虑可供选择的比率限制.默认的比率限制同样可用于热分析和流体分析的稳定状态求解.对于热分析,是对不平衡热流量和所加的热流量量级进行评估,而不是力.对于流体分析,对不平衡流度和所加流度量级进行评估.(5)考虑选择阻尼:对于静力分析,默认的阻尼是局部阻尼,对于消除大多数网格顶点的速度分量周期性为零时的动能很有效.这是因为质量的调节过程依赖于速度的改变.局部阻尼对于求解静力平衡是一个非常有效的计算法则且不会引入错误的阻尼力(见Cundall1987).如果在求解最后状态,重要区域的网格海域的速度分量不为零,那么说明默认的阻尼对于达到平衡状态是不够的.有另外一种形式的阻尼,叫组合阻尼,相比局部阻尼可以使稳定状态达到更好的收敛,这时网格将发生明显的刚性移动.例如,求解轴向荷载作用下桩的承载力或模拟蠕变时都可能发生.使用SETmechanicaldampcombined命令来调用组合阻尼.组合阻尼对于减小动能方面不如局部阻尼有效,所以应注意使系统的动力激发最小化.可以用SETmechanicaldamplocal命令转换到默认阻尼.(6)检查模型反应:FLAC3D显示了一个相试的物理系统是怎样变化的.做一个简单的试验证明你在做你认为你在做的事情.例如,如果荷载和实体在几何尺寸上都是对称的,当然反应也是对称的.改变了模型以后,执行几个时步(假如,5或10步),证明初始反应是正确的,并且发生的位置是正确的.对应力或位移的期望值做一个估计,与FLAC3D的输出结果作比较.如果你对模型施加了一个猛烈的冲击,你将会得到猛烈的反应.如果你对模型作了一些看起来不合理的事情,你一定要等待奇怪的结果.如果在分析的一个给定阶段,得到了意外值,那么回顾到这个阶段所用的时步.在进行模拟前很关键的是检查输出结果.例如,除了一个角点速度很大外,一切都很合理,那么在你理解原因前不要继续下去.这种情况下,你可能没有给定适当的网格边界.(7)初始化变量:在模拟基坑开挖过程时,在达到目的前通常要初始化网格顶点位移.因为计算次序法则不要求位移,所以可以初始化位移,这只是由网格顶点的速度决定,并有益于用户初始化速度却是一件难事.如果设定网格顶点的速度为一常数,那么这些点在设置否则前保持不变.所以,不要为了清除这些网格的速度而简单的初始化它们为零...这将影响模拟结果.然而,有时设定速度为零是有用的(例如,消除所有的动能).(8)最小化静力分析的瞬时效应:对于连续性静力分析,经过许多阶段逐步接近结果是很重要的...即,当问题条件突然改变时,通过最小化瞬时波的影响,使结果更加“静力”.使FLAC3D解决办法更加静态的方法有两种.(a)当突然发生一个变化时(例如,通过使区域值为零模拟开挖),设定强度性能为很高的值以得到静力平衡.然后为了确保不平衡力很低,设定性能为真实值,再计算,这样,由瞬时现象引起的失败就不会发生了.(b)当移动材料时,用FISH函数或表格记录来逐步减少荷载.(9)改变模型材料:FLAC3D对一个模拟中所用的材料数没有限制.这个准则已经尺寸化,允许用户在自己所用版本的FLAC3D中最大尺寸网格的每个区域(假如设定的)使用不同的材料.(10)运行在现场原位应力和重力作用下的问题:有很多问题在建模时需要考虑现场原位应力和重力的作用.这种问题的一个例子是深层矿业开挖:回填.此时大多数岩石受很高的原位应力区的影响(即,自重应力由于网孔尺寸的限制可以忽略不计),但是回填桩的放置使自重应力发展导致岩石在荷载作用下可能坍塌.在这些模拟中要注意的重点(因为任何一种模拟都有重力的作用)是网格的至少三个点在空间上应固定...否则,整个网格在重力作用下将转动.如果你曾经注意到整个网格在重力加速度矢量方向发生转动,那么你可能忘记在空间上固定网格了.FLAC3D主要适明模拟计算地质材料和岩土上程的力学行为。

并从 中总结 线路 交 叉对桥 墩 沉降 的影 响 。 关键 词 京 沪 高铁 桥墩沉降 Ｆ Ｌ Ａ Ｃ ３ Ｄ 数值 分析 中图分 类号 ： Ｔ Ｕ ４ ３ ３ 文 献标 识码 ： Ａ
京 沪 高铁 是我 国中长期 铁路 网规 划 中投 资规 模最
大、 技术含量最高的一项工程 ， 全线优先采用以桥代路 的方式 ， 最 大 限度节 约 宝贵 的土 地资 源 ， 桥梁 长度 约 占 正线长度的 ８ ６ ． ５ ％， 故桥墩沉降问题是关系到线路和 列 车安 全 运行 的重 大技 术 问题 。通过 数值 计算 软 件对 京沪高速铁路某桥墩沉降量进行了数值建模分析。
点力取 ５ ３ ０ ｋ Ｎ， 计算沿线路方 向线荷载为 １ ７ ６ ｋ Ｎ ／ ｍ，
地 表压力 为 １ ７ ． ６ ｋ Ｐ ａ 。
可以模拟岩土或其他材料 的三维力学行 为。Ｆ Ｌ Ａ Ｃ ３ Ｄ
收稿 日期 ： ２ ０ １ ３— ０ ５—３ １ 第一作者简介 ： 闫东旭 （ １ ９ ８ ４ 一） ， 男， ２ ０ １ ２年毕业 于河海 大学岩土 工程 专业 ， 工 学博士 ， 工程师。
由于新建德州南 联络线需进 行下挖 和增 建雨棚 柱， 数值计算分三种工况进行 ： ①无新建联络线影响工 况； ②新建联络线开挖工况 ； ③新建联络线开挖后加载 工况 。联络线与京沪高速铁路交叉状况如 图 １ 所示 ，
线路交叉对京沪高铁桥墩沉降影响数值分析 ： 闰东旭
５
测工程地质报告［ Ｒ ］ ． 天津 ： 铁三院 ， ２ ０ １ ２
７ ． ２ １
４ ． ８ ５
２ ． ４ ６
１ ． ６ ６
１ ６ ． １
１ ３ ． ５
６ ７
５ ３

选择底面的线将底面的线分成4份然后在弹出的对话框中设为20再选择侧面的线把侧面的线分成203选择划分单元的类型如下图所示4点击meshinggenerate划分出图示网格
3D 非线性简单算例
1、问题描述 模拟土体在逐渐被压的过程中，应力应变变化过程。在本算例中，本构模型 采用非线性弹性 Duncan-Chang 模型。 其中，Duncan-Chang 模型计算参数取值，见下表： K n
（４）点击“Meshing”—〉 “Generate” ，划分出图示网格。
（5）把前处理数据转化为SciFEA 所需计算格式。点击【Files】-【Save】 菜单，保存为GID 文件。点击菜单【Calculate】-【Calculate】，弹出“process info”消息框，如图所示，点击“ok”按钮，然后退出Gid，至此前处理工作结 束。
(。 )
c
Rf
G
F
D
2.0e2 0.2
32 5e4 0.8 0.3 0.15 3.0
2、求解步骤 2.1 选择项目 （1）启动 SciFEA，选择“项目”-> “新建项目” 菜单或者选择新建项目按 钮“ ” ，弹出如下对话框。
(2)在“项目名称”后输入该项目的名称，再点击“问题类型”栏中的“二 维 Duncan”选项。如图所示。 (3)点击“OK”按钮完成项目类型的选择。
第二步，点击鼠标右键，再点击“Contextual”—〉 “Join C-a”
第三步，在图示的平面内点击一点，就可以显示出这一点的应力应变变 化曲线图，如下图所示。
第四步，同时画出第二主应力Ｓ２和第三主应力Ｓ３的曲线，如图所示。
本例中，Ｓ２的曲线和Ｓ３的曲线重合。 第五步，画出剪应力的变化曲线，点击“View results”—〉 “Graphs” —〉 “Point Analysis”—〉 “unodm0”—〉 “tr”,如图所示。

文章编号:1009-6825(2011)07-0236-02基于F L A C 3D软件的土石坝应力应变分析收稿日期:2010-11-13作者简介:孙展杰(1983-),男,助理工程师,塔里木河流域管理局,新疆库尔勒　841000孙　倩(1987-),女,塔里木河流域管理局,新疆库尔勒　841000孙展杰　孙　倩摘　要:利用有限差分计算分析软件F L A C 3D对某土石坝坝体进行了应力与变形非线性有限差分计算,得到坝体在施工期与稳定蓄水期的应力与变形值,计算结果表明,在现有可能不利条件下,该坝坝体、坝基及复合土工膜的应力变形均在合理可控范围之内,大坝结构安全能够满足要求。

关键词:F L A C 3D软件,土石坝,非线性有限差分,应力应变分析中图分类号:T V 312文献标识码:A1　概述土石坝是世界上最古老、应用最广的一种坝型,具有就地取材、结构简单、施工方便等优点。

但是由于土石料性质极其复杂,使土石坝在计算中的应力、位移和实际中的应力、位移相符合尚有一定的难度。

随着计算机技术的发展和广泛应用,各种数值分析方法开始应用到岩土工程中来,其中有限单元法因为其对边界条件的适应性以及能很好的模拟各种土的本构关系,使其在岩土工程计算中得到广泛的应用。

F L A C 3D(F a s t L a g r a n g i a n A n a l y s i s o f C o n t i n u a ,连续介质快速拉格朗日分析)是由C u n d a l l 和美国I T A S C A 公司开发出的有限差分数值计算程序,主要适用于地质和岩土工程的力学分析。

该程序自1986年问世后,经不断改版,已经日趋完善。

前国际岩石力学学会主席C .F a i r h u r s t 评价它:“现在它是国际上广泛应用的可靠程序”[1]。

F L A C 程序中提供了由空模型、弹性模型和塑性模型组成的10种基本的本构关系模型,所有模型都能通过相同的迭代数值计算格式得到解决:给定前一步的应力条件和当前步的整体应变增量,能够计算出对应的应变增量和新的应力条件。

应力时空演化模型全文共四篇示例，供读者参考第一篇示例：应力时空演化模型是一种用来描述和分析物质在时空中受力的变化和演化规律的理论模型。

在现实世界中，物质受到各种各样的力的作用，有时候这些力可能是由外部施加的，比如重力、电磁力等，有时候则是由物质内部的相互作用引起的，比如分子之间的相互作用、原子核的相互作用等。

这些力会使得物质的形状、结构、性质等发生变化，而应力时空演化模型正是用来描述和预测这种变化的模型。

应力时空演化模型是一个综合了物理力学、力学、热力学等多个领域知识的理论模型，它不仅能够描述物质受力的变化规律，还可以分析物质的形变、应力分布、应变分布等变化过程。

在应力时空演化模型中，物质被看作是由许多微观粒子组成的，每个微观粒子都受到各种各样的力的作用，这些力的大小和方向会导致物质的形变和应力分布的变化。

在应力时空演化模型中，时空的概念非常重要。

时空是描述事物的位置和状态的基本要素，而应力时空演化模型正是用来描述物质的时空变化规律的。

在这个模型中，时空被看作是一个连续不断的整体，物质在时空中的运动和演化是由一系列微分方程来描述的，这些微分方程以时间和空间为自变量，表示了物质的运动和演化规律。

应力时空演化模型的基本假设是物质是连续的，力的作用是均匀的，物质是各向同性的等。

在这些基本假设下，可以建立起物质的应力时空演化模型，用来描述和分析物质在时空中受力的变化和演化规律。

这个模型可以应用在各种各样的领域，比如材料科学、土壤力学、岩土工程等等。

应力时空演化模型的建立和应用是一个复杂而繁琐的过程，需要大量的实验数据和数学计算来支撑。

在建立模型的过程中，往往需要利用现代科学技术的手段，比如计算机模拟、数值计算等，来对物质的力学性质进行研究和分析。

在模型应用的过程中，需要对模型的参数进行优化和调整，以使模型能够更好地描述和预测物质的力学行为。

应力时空演化模型是一种用来描述和分析物质受力的变化和演化规律的理论模型，它能够帮助我们更好地理解和预测物质的形变和应力分布等性质。

j-c本构公式参数
J-C本构公式是一种描述材料变形行为的数学模型，其参数是根据材料的实验数据进行拟合得到的。

J-C本构模型常用的参数包括：
1. 弹性模量（Young's modulus）：描述材料在受力时的刚度，代表了单位应变下的应力变化率。

2. 屈服应力（Yield stress）：材料在开始塑性变形之前所能承受的最大应力。

3. 塑性流动指数（Plastic flow exponent）：描述材料的塑性变形特性，一般取值在0-1之间。

4. 塑性应力指数（Plastic stress exponent）：描述材料的塑性变形特性，一般取值在1-10之间。

5. 蠕变指数（Creep exponent）：描述材料在长时间持续应力下发生的蠕变变形，一般取值在1-10之间。

6. 洛伦兹-杨比（Lorenz-Yang ratio）：描述材料的半非晶性，一般取值在1-10之间。

7. 晶粒尺度指数（Grain scale exponent）：描述材料的晶粒尺度对材料力学性能的影响，一般取值在0-1之间。

8. 温度指数（Temperature exponent）：描述材料的温度对材料力学性能的影响，一般取值在0-1之间。

以上参数是针对J-C本构模型的常用参数，具体的参数取值需要根据实验数据进行拟合。

应力应变公式曲线方程应力应变公式是描述材料在受力作用下产生的变形的数学表达式。

它是材料力学中最基本且重要的方程之一，可以用来研究材料的力学性质和预测材料的变形行为。

应力应变公式的研究在工程设计、材料科学、结构力学等领域具有重要的理论和应用价值。

首先，我们来了解应力应变公式的基本概念和意义。

应力是指材料单位面积上承受的力，通常用σ表示，单位是帕斯卡（Pa）。

而应变是指材料在受力作用下的变形程度，通常用ε表示，它是一个无量纲的比值。

应力和应变之间的关系可以通过应力应变公式来表达。

应力应变公式一般可以表示为σ=Eε，其中E是材料的弹性模量，代表材料的刚度和弹性性能。

弹性模量越大，材料的刚度越高，变形程度越小；弹性模量越小，材料的变形程度越大。

这个公式告诉我们应力和应变之间的关系是线性的，材料在弹性范围内可以按照线性关系变形。

然而，事实上，材料在受力作用下，并不总是按照线性关系变形。

很多材料在受力后会出现变形的非线性现象，这时候就需要引入非线性应力应变公式来描述材料的变形行为。

一般来说，非线性应力应变关系可以表示为σ=σ0+Kε^n，其中σ0代表应力偏移量，K代表应力与应变之间的系数，n代表非线性指数。

在实际应用中，根据不同材料的力学性质和应变特点，可以选择不同的应力应变公式来描述材料的变形行为。

例如，对于弹性材料来说，可以选择线性应力应变公式；对于塑性材料来说，可以选择非线性应力应变公式。

在材料设计和结构分析中，正确选择并应用适合的应力应变公式，可以更准确地预测和分析材料的变形行为，为工程设计提供可靠的依据。

除了应力应变公式，还有一些与之相关的概念和重要参数需要考虑。

例如，屈服强度是指材料在允许的变形范围内承受的最大应力；断裂强度是指材料在断裂前能承受的最大应力；刚度是指材料在受力下的抵抗能力；蠕变是指材料长时间作用下的变形现象等等。

这些概念和参数可以从不同角度对材料的力学性能进行研究和评价。

在工程实践中，应力应变公式的研究和应用可以用于材料的选取、结构的设计和分析以及性能的评估等方面。

材料变形与力学性能分析的数值模拟研究引言：在材料科学领域，对材料的变形及其力学性能进行分析和研究是至关重要的。

通过数值模拟方法可以模拟和预测材料在不同加载条件下的力学行为，为工程设计和材料优化提供指导。

本文将重点探讨材料变形与力学性能分析的数值模拟研究，并讨论其应用于实际工程中的意义与价值。

一、材料变形的数值模拟方法材料变形的数值模拟主要通过有限元法来实现。

有限元法是一种基于数值分析的方法，通过将大型物体或结构分解为有限个小的单元，再根据单元间的相互关系进行计算，从而得到整体的变形和力学性能。

常见的有限元法包括线性有限元法、非线性有限元法和动态有限元法。

其中，非线性有限元法可以更准确地模拟材料在高应变下的行为。

二、材料变形过程的数值模拟研究材料的变形过程通常涉及应力、应变和位移等参数，这些参数的变化对材料的力学性能产生重要影响。

数值模拟研究可以通过对材料的载荷施加和边界条件设置来模拟材料在外力作用下的变形行为。

利用有限元法的数学模型，可以预测和分析材料在不同工况下的应力分布、应变变化和位移分布。

这种数值模拟方法可以帮助工程师和科学家理解材料的变形机制，并优化设计，降低结构或零件的失效风险。

三、数值模拟在研究材料力学性能中的应用1. 材料强度分析：数值模拟可以通过载荷施加和边界条件的设定，模拟材料在不同应力下的应变变化。

通过分析应力-应变曲线，可以得到材料的弹性模量、屈服强度、抗拉强度等力学性能指标。

这些性能指标的研究可以帮助工程师选择最合适的材料，合理设计结构，确保产品的可靠性和安全性。

2. 材料疲劳分析：材料在长时间的循环加载下会出现疲劳失效。

数值模拟可以模拟材料在循环加载下的变形过程，预测疲劳寿命，并通过优化设计和材料选择来延长材料的使用寿命。

疲劳分析的数值模拟在航空航天、汽车、桥梁等领域具有重要应用价值。

3. 材料塑性变形分析：在一些特定加载条件下，材料会发生塑性变形，从而改变材料的形状和性能。

β= σ
fc
σ
fc
二维非线性指标
β=
OP σ2 σ = 1 = σ 2 f σ 1 f OF
!
三维非线性指标：Ottosen法
保持σ1, σ2不变，改变σ3 直至与破坏面相交得到 交点(σ1, σ2, σ3f)
β= σ3 σ3f
!
β = 1 处于破坏状态
ε
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》
混凝土三维本构模型的核心
!
各种本构模型的本质差别
!
!
!
混凝土的应力应变关系，主要是建立在 一维情况下 寻找合适的指标，将一维的应力应变关 系拓展到三维 非线性指标：
! ! !
!
K=
E 3(1 − 2ν )
E 2(1 + ν )
E=
ν=
9 KG 2K + G
G=
3K − 2G 2(3K + G )
非线性弹性本构 β ε pl 弹塑性本构 损伤力学本构 D
清华大学研究生课程——《钢筋混凝土有限元》
清华大学研究生课程——《钢筋混凝土有限元》
全量模型
!
Cedolin 模型
σ oct = 3K s ε oct τ oct = 3Gs γ oct
Ks = ab −ε oct / c + d K0 Gs = pq −γ oct / v + sγ oct + t G0
β=
J2 J2f
=
5.292 = 0 .587 9.018
β < 0.8
ν s = ν 0 = 0.18
清华大学研究生课程——《钢筋混凝土有限元》

焊接应力和变形的数值模拟研究焊接应力和变形的数值模拟研究摘要：焊接是一种重要的金属连接方式，在工程领域得到广泛应用。

然而，焊接过程中产生的应力和变形可能会导致焊接接头的质量不达标或功能失效。

为了实现高质量的焊接，数值模拟成为了研究焊接应力和变形的重要工具。

本文通过对焊接过程中应力和变形的数值模拟研究进行总结和分析，旨在深入探讨焊接过程中应力和变形的产生机制和影响因素，为实现高质量的焊接提供理论指导和工程应用建议。

1. 引言焊接是一种通过热和压力将金属材料连接在一起的加工工艺。

焊接过程中，由于高温和热应力的作用，焊接接头容易产生应力和变形。

这些应力和变形可能会影响焊接接头的力学性能和寿命，甚至导致焊接接头的破裂。

因此，研究焊接应力和变形的产生机制和特性，对于确保焊接接头的质量和可靠性具有重要意义。

2. 焊接应力和变形的数值模拟方法数值模拟是研究焊接应力和变形的有效手段之一。

常用的数值模拟方法包括有限元方法、有限差分法、边界元方法等。

其中，有限元方法是最常用的数值模拟方法之一，具有计算精度高、计算效率高等优点。

本文将主要以有限元方法为基础，进行焊接应力和变形的数值模拟研究。

3. 焊接应力和变形的产生机理焊接过程中，应力和变形的产生主要是由于以下几个因素的综合作用：焊接热源的热输入、材料的热物理性质、焊接接头的几何形状等。

具体而言，焊接热源的热输入会引起焊接接头的温度场分布的变化，从而导致焊接接头内部发生热膨胀或热收缩，产生应力和变形。

另外，焊接接头在冷却过程中也会发生冷却收缩，进一步增加焊接接头的应力和变形。

4. 焊接应力和变形的影响因素焊接应力和变形的程度和特性受多种因素的影响，包括焊接接头的几何形状、材料的物理性质、焊接过程的参数等。

焊接接头的几何形状对于焊接应力和变形的分布和大小具有重要影响。

一般来说，焊接接头的几何形状越复杂，焊接应力和变形的程度越大。

材料的物理性质也会影响焊接应力和变形的特性，包括热膨胀系数、热导率等。

FLAC数值模拟软件及其在地学中的应用
龚纪文;崔建军;席先武;林舸
【期刊名称】《大地构造与成矿学》
【年(卷),期】2002(26)3
【摘要】FLAC是以拉格朗日差分法为基础的一种数值模拟应用软件,主要用于岩土力学工程.近年地质学家将其应用到地质构造和块体运动学、动力学研究中,以证实或证伪地球学家们提出的相关地球动力学模型.本文初步介绍了FLAC应用软件的基本原理、使用方法以及它在国内外地质学上的几个应用实例.
【总页数】5页(P321-325)
【作者】龚纪文;崔建军;席先武;林舸
【作者单位】中国科学院广州地球化学研究所,广东,广州,510640;华北石油地质局第九普查大队,内蒙古,榆次,030600;中国科学院广州地球化学研究所,广东,广州,510640;中国科学院广州地球化学研究所,广东,广州,510640
【正文语种】中文
【中图分类】TP392
【相关文献】
1.拉格朗日元法及其应用软件FLAC在边坡稳定分析中的应用 [J], 杨丽萍;宋子岭;吴野
2.FLAC基本原理及其在地学中的应用 [J], 杨立强;张中杰;林舸;王岳军
3.浅谈数值模拟软件FLAC3D在某边坡工程中的应用 [J], 苏恩华;王胜
4.强度折减法在FLAC 3D数值模拟中的应用 [J], 章显斌;余以强
5.数值模拟软件FLAC及其在地学应用简介 [J], 谢建华;夏斌;徐振华;张宴华
因版权原因，仅展示原文概要，查看原文内容请购买。

焊接是一种常见的金属加工工艺，它通常用于将两个或更多金属件连接起来。

在焊接过程中，金属会受到热量的影响，从而产生应力和变形。

为了更好地理解焊接过程中的应力和变形机理，以及预测焊接件的性能和寿命，数值模拟研究成为了焊接领域的研究热点之一。

1. 应力分布的数值模拟焊接过程中，焊缝和母材会受到热量的影响，产生应力。

通过有限元分析等数值模拟方法，可以准确地预测焊接件中的应力分布。

这对于避免焊接件的变形、裂纹和疲劳寿命的提高至关重要。

数值模拟可以帮助工程师优化焊接工艺参数，减少焊接应力，提高焊接件的质量和性能。

2. 变形控制的数值模拟除了应力之外，焊接过程中还伴随着焊接件的变形。

焊接变形可能导致产品尺寸的偏差，从而影响其装配质量和外观。

数值模拟可以帮助工程师预测焊接件的变形情况，优化焊接设计，减少变形产生的影响。

通过数值模拟，工程师可以选择合适的焊接序列、布局和残余应力的预处理方法，从而有效地控制焊接变形。

3. 焊接残余应力的数值模拟焊接过程中产生的应力不仅会影响焊接件的性能和寿命，还会导致焊接残余应力的存在。

焊接残余应力可能导致产品的破坏和失效，因此需要对其进行有效的控制。

数值模拟可以模拟焊接残余应力的分布和大小，帮助工程师选择合适的残余应力消除方法，如后续热处理、切割释应力等，从而提高焊接件的质量和可靠性。

总结数值模拟是研究焊接应力和变形的重要手段，通过数值模拟，工程师可以更好地理解焊接过程中的物理现象，预测和控制焊接件的应力和变形，提高焊接件的质量和性能。

相信随着数值模拟技术的不断发展和完善，焊接应力和变形的研究将会取得更加深入和全面的成果，为焊接工艺的改进和创新提供更可靠的技朧支撑。

在焊接工艺中，焊接应力和变形的研究一直是焊接工程领域的一个重要课题。

在实际工程中，焊接应力和变形的控制对于确保焊接件的质量、性能及使用寿命至关重要。

针对焊接过程中产生的应力和变形问题，数值模拟成为了研究人员以及工程师进行预测和优化的重要工具。

摘 要 ：矿 山压 力 是 影 响底 板 承 压水 突 出 的 重 要 因 素 ， 是底板破坏 的力源 ； 了 解 煤 层 底 板 的 应 力 分 布
是 研 究 煤 层 底 板 破 坏 深 度 的 重要 前 提 ， 有 利 于加 深 对 突水 机 理 的 认 识 。本 文 采 用 Ｆ Ｌ ＡＣ模 拟 软 件 求 解 采
应力逐步增加 ， 底板 垂直应力 ｄ ｙ为 底 板 深 度 Ｘ 的 二 次 函数 ， 而 垂 向位 移 逐 步 减 小 ， 底 板垂直方 向的应 变 ￡ 为底板深度 Ｘ的负指数函数 。 关 键 词 ：采 场 底 板 ；Ｆ Ｌ ＡＣ；应 力 一 应 变 ；本 构 关 系 中图 分 类 号 ：Ｔ Ｄ３ ５ ３ 文 献 标 识 码 ：Ａ 文 章 编 号 ：１ ０ ０ ４ — ４ ０ ５ １ （ ２ ０ Ｉ ３ ） １ １ － ０ ０ ６ ３ — ０ ３
（１ ． Ｋｅ ｙ Ｌａｂ ．ｏ ｆ Ｍｉ ｎｅ Ｄｉ ｓ ａ ｓ ｔ ｅ ｒ ｐ ｒ ｅ ｖ ｅ ｎｔ ｉ ｏｎ ａ ｎｄ ｃ ｏｎ ｔ ｒ ｏ ｌ ，Ｓ ａｎ ｈｅ ０ ６５ ２０ １， Ｃｈｉ ｎａ； ２ ． Ｓａ ｆ ｅ ｔ ｙ Ｅｎｇｉ ｎｅ ｅ ｒ ｉ ｎｇ Ｃｏ ｌ ｌ ｅ ｇｅ，
场 的 围岩 应 力 一 应变关 系， 分 析 了不 同 阶段 采 空 区 底 板 应 力 变 化 过 程 ， 并 根 据 底 板 不 同 深 度 的 垂 直 应 力 与
垂直位移值 ， 运用数 学方法建立应力 、 应 变 随 采 场 深 度 变 化 的 本 构 方 程 。结 果 表 明 ： 随 着 深 度 的增 加 ， 垂 直
Ｎｏ ｒ ｔ ｈ Ｃｈ ｉ ｎ ａ Ｉ ｎ ｓ ｔ ｉ ｔ ｕ ｔ ｅ ｏ ｆ Ｓ ｃ ｉ ｅ ｎ ｃ ｅ ａ ｎ ｄ Ｔｅ ｃ ｈ ｎ ｏ ｌ ｏ ｇ ｙ ，Ｓ ａ ｎ ｈ ｅ ０ ６ ５ ２ ０ １ ，Ｃ ｈ ｉ ｎ ａ ；３ ．Ｓ ｔ ａ ｔ ｅ Ｋｅ ｙ Ｌａ ｂ ｏ ｒ ａ ｔ ｏ ｒ ｙ ｏ ｆ Ｃｏ ａ ｌ Ｒｅ ｓ ｏ ｕ ｒ ｃ ｅ ｓ

ＦＩ ＡＣ（ ｓ ｇ ａ ｇａ ａ ｙ ｉ ｆ Ｃｏ ｔｎ ａ 连 Ｆａ ｔ Ｌａ ｒ ｎ ｉｎ Ａｎ ｌ ｓｓ ｏ ｎ ｉ ｕ ，
１ Ｆ ＡＣ 软 件 的 基 本 原 理 Ｌ
Ｆ ＡＣ软 件 的基 本 原 理 即 是 拉 格 朗 日差 分 法 。 Ｌ
拉 格 朗 日元 法 源 于 流 体 力 学 。在 流 体 力 学 中 有 两
工 程 问日分 析 ） 是 基 于 该 方 法 的显 式 正 有限差 分程 序 ， 不但 能处理一 般 的大 变形 问题 ， 它 而 且 可 以 在 模 型 中 加 人 节 理 、弱 面 等 地 质 构
造 ７ 。
在 国外 Ｆ Ｉ ＡＣ 软 件 得 到 了 较 大 的 应 用 ， 在 而
种 主 要 的 研 究 方 法 ， 种 是 定 点 观 察 法 ， 称 欧 拉 一 亦 法 ； 一 种 是 随 机 观 察 法 ， 为 拉 格 朗 日法 。后 者 另 称 是 研 究 每 个 流 体 质 点 随 时 间 而 变 化 的 状 态 ， 研 究 即 某一流体质点 在任 一段 时间 内的运 动轨迹 、 度 、 速 压 力 等 特 征 。把 拉 格 朗 日法 移 值 到 固体 力 学 中 ， 把 所 研 究 的 区域 划 分 成 网格 ， 结 点 就 相 当于 流 体 质 其 点 ， 后 按 时步 用 拉 格 朗 日法 来 研 究 网 格 结 点 的 运 然 动 ， 种 方 法 就 是 拉 格 朗 日元 法 。 它 的 优 点 是 占用 这 内存 少 ， 解 速 度 快 ， 于 用 微 机 求 解 较 大 规 模 的 求 便
拉 格 朗 日元 法 是 ～ 种 利 用 拖 带 坐 标 系 分 析 大
９ Ｏ年 代 初 ， 内才 引 进 该 软 件 ， 要 应 用 于 岩 土 力 国 主 学 分 析 ， 矿 体 滑 坡 、 矿 开 采 沉 陷 预 测 、 利 枢 纽 例 煤 水

IJKF 软件由美国 R?57@5 公司于 (; 世纪 S; 年 代开发并于 (; 世纪 2; 年代得以广泛利用。目前该 软件在国外已被广泛应用于工程地质、岩土力学以 及构造地质学和成矿学等研究领域 T 82SH，UB@VQ W5&， !" #$% M 82SS，F%AX5$$，!" #$% M 822;，45?BA# Y#%@’，!" #$% M 8223，EW5A6，!" #$% M 822)，EW5A6， !" #$% M 822)，K,A#$X， !" #$% M 822S，EW5A6，!" #$% M 8222，Z$B[’,，!" #$% \ 。国内于 (; 世纪 2; 年代初才 引进该软件，主要应用于工程地质和岩土力学分析， 如矿体滑坡、煤矿开采沉陷预测、水利枢纽岩体稳定 性分析、采矿巷道稳定性研究等 T 822S，吴洪词M
8 概述
数值模拟是以计算机软件进行数值分析的一 种方法。经典的数值分析方法主要有拉格朗日法和 欧拉法。IJKF ! I57? J56,5A6B5A KA5$O7B7 #P F#AQ ?BA%5，连续介质快速拉格朗日分析1 即是基于拉格朗 日差分法的一种显式有限差分程序，它不仅适宜于 处理大尺度、大变形工程和地质问题，而且可以在初 始模型中加入诸如断裂、节理构造等地质因素。
#. ) $%&’ 软件概述 $%&’ 软件是以命令驱动的程序，命令驱动程序
使得 $%&’ 在工程分析方面更加通用。其最大的优 点在于在分析或运行过程中，可以随时查看运行过 程中某个状态的矢量分布，或查看某个矢量的演变 历程，同时亦可将该状态下的计算结果保存成 /&0 文件供后期处理。$%&’ 程序运行时，它输入的是命 令，其输入方法为：键盘交互式输入和数据文件输 入。在 $%&’ 中有 12 多个主命令和近 122 个命令变 量，但是常用的主命令和命令变量仅几十个。可以使 用或产生的文件类型包括：初始化文件 3 4567. 898 : 、数 据文件 3 ". ;6< : 、保存文件 3 ". =6> : 、日志文件 3 4567. 5?@ : 、历史文件 3 4567. A8= : 和动画文件 3 4567. B?> : 。

采矿工程数值计算方法——FLAC建模技巧与工程应用FLAC词汇讲解：Bulk 体积模量COhesion 粘结力Den 质量密度Dilation 剪胀角E_p 全塑性应变(应变软化模型)Friction 摩擦系数SHear 剪切模量X X-坐标Y Y-坐标PP 孔隙压力JFric 节理摩擦力JCoh 节理粘结力JAngle 节理角度NUYx y-x 泊松比NUZx z-x 泊松比SIG1 最大主应力SIG2 最小主应力SXX XX-应力SYY YY-应力SXY XY-应力Theta 最小主应力与x-轴的夹角XDis X-位移YDis Y-位移XMod X-模量YMod Y-模量XVel X-速度YVel Y-速度ASXX XX-应力(三角形A)BSXX XX-应力(三角形B)CSXX XX-应力(三角形C)DSXX XX-应力(三角形D)ASYY YY-应力(三角形A)BSYY YY-应力(三角形B)CSYY YY-应力(三角形C)DSYY YY-应力(三角形D)ASXY XY-应力(三角形A)BSXY XY-应力(三角形B)CSXY XY-应力(三角形C)DSXY XY-应力(三角形D)State 塑性状态0 弹性1 正在屈服中2 曾经屈服, 现为弹性3 已超过单轴抗张力4 屈服并超过单轴抗张力5 已经超过抗张力6 彻体节理正在屈服7 彻体节理过去屈服过, 现为弹性状态1 FLAC建模方法1.1 建模Flac命令汇总：横观各向同性弹性模型ANisotropic彻体节理模型Ubiquitous输出图形set plot jpgPl hard file 文件名（edit-copy to clipboard）create Trench（创建了一个叫做“TRENCH”的视图窗口）add surface yellow（表面显示黄色）显示锚杆pl add sel force在屏幕上显示：plot show显示坐标轴：plot add axes black设置模型的颜色：plot add surface yellow不同的组显示不同的颜色：plo blo group设置背景颜色：pl set back white设置图片大小：set plot bitmap size (800,510) ;设置图例数字大小(default=35 范围10-50)：plot set caption size 37command window ：mainwin size 1 1 position 0 0.8graphics window ：plot set window size 0.8 0.8 position 0.1 0.0plot set center 100 0 40如何进行切片：plo set plane ori (点坐标) norm (法向矢量)plo set plane ori 176.7621 426.9150 21.649 norm 1 0 0显示切片pl bl gr p ou off dis p scale 0.02 ;pl con disp plane outline off dis plane 显示变形轮廓线的命令plo ske magf 10显示模型：Pl bl gr显示点的编号：plo blo gro gpnum on显示网格：Plot grid显示边界：Plot bo查看模型的边界情况: plo gpfix red sk查看模型的体力分布：plo fap red sk求解边坡的安全系数solve fos显示安全系数plot fos/print fos查看剪应变：pl con ssi outline on vel red查看塑性区：pl bl st模型的云图：位移：plo con dis (xdis, ydis, zdis)plo bcon dis应力：call(sy, sx, sxy, syz, sxz)查看最小主应力pl con smin查看最大主应力pl con smax最大拉应力云图：plot cont smax最大压应力云图：plot cont smin查看szz的应力云图：plo con szz ou on模型的矢量图：plo dis (xdis, ydis, zdis)看模型有多少单元、节点：plo info查询节点的坐标pr gp p range id xxxx输出模型的后处理图：File/Print type/Jpg file，然后选择File/Print，将保存格式选择为jpe文件施加面力：app nstress调整视图的大小、角度：综合使用x, y, z, m, Shift键，配合使用Ctrl+R，Ctrl+Z等快捷键单元的ID：用鼠标双击单元的表面，可以知道单元的ID和坐标进行切片：plo set plane ori (点坐标) norm (法向矢量)plo con sz plane (显示z方向应力的切片)file里print setup选项可以设置输出图片的大小和质量调用已保存的结果：rest +文件名；或者File / Restor暂停计算：Esc程序中进行暂停：命令中加入pause命令，用continue进行继续把一个PLOT的图像数据导出来以便用其他软件绘图：用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。