flac--学习

FLAC原理实例与应用指南(附光盘一张)作者：刘波，韩彦辉编著 出版社：人民交通出版社 1.FLAC简介FLAC——FastLagrangianAnalysisofContinuaFLAC建立在拉格朗日算法基础上，采用有限差分显式算法来获得模型全部运动方程(包括内变量)的时间步长解，从而可以追踪材料的渐进破坏和垮落，这对研究岩土工程设计是非常重要的。FLAC适用模拟计算岩土材料力学行为，特别适合模拟大变形和扭曲，包括材料的高度非线性(应变硬化/软化)、不可逆剪切破坏和压密、粘弹(蠕变)、孔隙介质的应力—渗流耦合、热—力耦合以及动力学问 题 快递公司问题件快递公司问题件货款处理关于圆的周长面积重点题型关于解方程组的题及答案关于南海问题 等。2.FLAC软件应用2.1建模（1）设计计算模型的尺寸（2）规划计算网格数目和分布（3）安排工程对象（开挖、支护等）（4）给出材料的力学参数（5）确定边界条件2.2网格生成 Grid i,j例如：grid30,202.3网格规划：Genx1,y1x2,y2x3,y3x4,y4   例如：Gen0,00,1010,2020,0 2.4 赋给单元材料性质mode (弹性)propd1800e-6bu12.5sh5.77i=1,20j=1,10propd2400e-6bu1250sh577 i=1,20j=11,20modm（弹塑性Mohr-Coulumb准则）propd1800e-6bu12.5sh5.77c0fri20ten0.015regi,j2.5 赋给模型边界条件（1）固定边界（结点）Fixxi=1,j=1,21Fixyi=1,21j=1（2）施加边界力(结点)applyyf -10i=1,21j=21 或applysyy -10i=1,21j=21applyxf -5i=21,j=1,21 或 applysxx -5i=21,j=1,212.6 赋给模型初始条件赋单元内应力（单元）inisxx-10 i=1,20j=1,20inisyy -5var04 i=1,21 j=1,212.7 计算Setgrav9.81SetlargeStep1000Savetest.sav2.8 结果显示Plotgrid 显示网格Plotbo 显示边界Plotplas显示塑性区Plotsig1fi显示最大主应力1Plotsig2fi显示最小主应力2Plotsdiffi显示主应力差(1-2)Plotstr 显示主应力矢量场Plotxdisfi 显示X方向位移Plotydisfi 显示Y方向位移Plotdisp 显示位移矢量场Plot e_p--塑性应变的轮廓线（只限于应变软化材料）；ff—以矢量形式画出固定结点的反力；pp—孔隙压力的等值线；wate—水位线；FLAC使用步骤FLAC使用步骤1.FLAC程式使用前准备步骤步骤1：依比例画出所欲分析的资料在纸上画出地点的位置、地层资料、并简单标示距离及深度资料。步骤2：换算输入资料成同一单位将现有地层资料，如Density,Bulkmodulus,Young`smodulus,tension,cohesion,frictionAngle等资料，换算成同一单位。　附注：需谨慎检查输入资料的单位，如因单位不同而造成过大或过小的值，将会造成FLAC无法计算，而产生ERROR信息。　  步骤3：应用公式简略计算应用公式或依据经验，简略算出FLAC输出资料的范围，以作为Debug及输出资料分析时验证。步骤4：建立x,y坐标与nodei,j之间的关系在图上距离及深度的关系，建立x,y坐标系统，再由x,y坐标系统，转换与网格间系，为了便于以后输出资料的分析，故应确实掌握网格的位置及其相对应的x,y坐标。建议在敏感区域使用较密的网格，其它地方则使用较疏的网格，刚开始执行程序时，不宜使用网格太大的网格数目，因尽量使网格总数少于1000，以节省时间。2.FLAC3D程序的编写步骤：1Config　________2Grid ________3Model　________4求起始的应力平衡(1)建立x,y坐标与网格的关系，建议使用Gen指示：Genx1,y1x2,y2,x3,y3x4,y4i=i0,i1j=j0,j1详细指令参见使用手册，FLAC程式可自动产生x,y坐标与网格的关系，但由于产生的网格坐标不易控制，将对其它的工作产生负面影响，故依使用前步骤4所建立的关系，将网格依其疏密程度需要的不同，实际控制网格的坐标。(2)设定材料性质：prop(3)设定外力：SetGrav,ApplyPressure,inisxx,Syy(4)设定边界条件：fix,free(5)求起始的应力平衡：solve(6)储存：Save5求工程的影响求出区域内的应力分布情况后，再依工程的 流程 快递问题件怎么处理流程河南自建厂房流程下载关于规范招聘需求审批流程制作流程表下载邮件下载流程设计 及步骤阶段执行各工程进行过程的影响，建议使用以下的步骤：(1)调出起初的应力平衡：re_____.sav(2)设定新的材料性质：model,prop(3)设定新的支撑性质：struct(4)设定新的外力apply(5)设定边界条件apply,fix(6)求工程时的应力平衡solve,step,cycle(7)储存savefile.sav执行变更材料开挖；·增加或者删除网格点荷载或压力；·改变某个单元的材料模型或者性质参数；·指定或者解除某个网格点的速率用modelnull命令进行开挖；用applyforce和applyyforce命令可以给任意网格点指定网格点荷载；用apply命令可以在模型边界上修改压力或者应力；用model和property命令来改变材料模型和性质；用fix/free命令指定或者解除网格点速率。一些难点：关于gen命令的使用：gen里的same并不是代 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 上面一个语句相同位置的点坐标。same可以代表前面你用过的GEN命令中的指定的任意一个点，但是前提是这两个点有共同的grid(I,J)号same指的是前面的任意一点，不一定是前1个点或者是后1个点，只要是按照顺时针建模型即可！！ApplySOLVE命令受限定的平衡比率（10-3），限定的不平衡力（100个力单位），限定的步数（100,000步）和限定的计算机运行时间（1440分钟）的控制，括号里的值是默认值。当达到其中任何一个限制时，计算就会停止。在上面的例子中，首先达到平衡比率10-3。想要用不平衡力来控制计算，应该在SOLVE命令之前输入SETsratio=50命令。现在计算就会在力达到100的限值时停止。可以用SET命令改变任何一个求解限值，如SETforce=50会将不平衡限值改为50。另外，也可用SOLVEforce=50命令。限制会继续有效，直到再次改变或者执行一个新的命令，这样就会将限值重新设置为默认值。当使用SOLVE命令时，检查计算是否过早的停止是非常重要的，例如，如果计算多于100,000步才能达到平衡，就要用SETstep命令来增加步骤的限值。菜单驱动形式记录窗口记录显示当前模型项目状态所需命令。这个记录可以输出作为一个数据文件，其中是一系列显示所正在分析问题的FLAC命令属性的代码被分配给单元通过propertynp。如果关键词没有被指定，则默认值np=1。通过structproperty=np命令为每个不同类型的单元分配几何和材料性质。属性号被分配给一个结构单元组。组是有一样属性号的同一类型结构单元片段的一组相近的设置。如果同一类型结构单元的带有一样的属性号两个单独的组被连接，它们将变成有相同组号的一个组。如果带有不同的属性号两个单独的组被连接，每个组仍保留其原有的组号。组的号可通过plotstructnumber确定。分配给结构节点和耦合弹力一属性号，属性号和与结构节点相连的单元片段相协调。任何属性特征与连接点FLAC中的7种加固结构单元1.梁单元.允许弯曲,并可以通过节点或接触面与介质（岩石/土）的网格连接。接触面连接允许脱开和滑移.塑性屈服受轴向荷载控制，可指定塑性矩;Beam2.衬砌单元.* 与梁单元类似,屈服准则中包括弯曲应力;Liner3.锚索单元.不抗弯.锚索节点在法向方向上依附于网格的变位，顺锚索的切向方向通过剪切和滑动弹簧单元与网格相互作用.允许轴向屈服.锚索节点也可以刚性地与网格点相连接;Cable4.桩单元. 抗弯.在法向和切向方向通过具有屈服功能的弹簧与网格连接;Pile5.锚杆单元.* 与桩相似可以考虑围压的变化、粘结材料的应变软化、以及单元的张裂;Rockbolt6.条形锚单元.* 与锚杆单元类似,但不抗弯.条形锚/土界面服从非线性的剪切破坏包线，该包线与围压相关;Strip7.Supportmembers.使用简单的1维弹簧单元连接(两个)自由面(上的点)（例如,可以模拟矿井巷道中的木支架支护).Support.*FLAC5.0新开发特征梁的keywordArea 梁的代表性区域E杨氏模量 应力/位移弹性模量height梁元件的代表性高度i梁元件的惯性矩pmomvalue 梁的塑性矩也就是力矩的极限，Radius半径spacingvalue元件在平面外的方向上的间隔默认为1thexpvalue 梁的线性热膨胀系数widthvalue 梁元件代表性宽度