1、1、怎样查看模型? 答:plot grid 可以查看网格,plot grid num 可以查看节点号。2、请问在圆柱体四周如何施加约束条件? 答:可以用 fix . ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中 r2gp0_disp gp0_disp=gp_disp endif p_gp = gp_next(p_gp) endloop ;找最大值的坐标 p_gp=gp_head ss=0 loop while p_gp # null gp_disp = gp_xdisp(p_gp)*gp_xdisp(p
2、_gp) gp_disp = gp_disp + gp_ydisp(p_gp)*gp_ydisp(p_gp) gp_disp = gp_disp + gp_zdisp(p_gp)*gp_zdisp(p_gp) gp_disp = sqrt(gp_disp) if gp_disp=gp0_disp ss=ss+1 xtable(1,ss)=gp_xpos(p_gp) ytable(1,ss)=gp_ypos(p_gp) endif p_gp = gp_next(p_gp) endloop end get_gp_maxdisp 18、如何显示计算后模型的变形状态?答:pl sk magfac *(
3、一个放大的数字)19、如何查看剖面上任一点的位移值 答:print gp disp range x() y( ) z( )或 print gp disp range id plot gp disp range x() y( ) z( )或 plot gp disp range id20、group 命令划分斜面 答:group groupname colour range plane above dd * dip * origin * plane below dd * dip * origin * dd 和 dip 就是两个斜面的倾向和倾角 ,origin 是这个斜面通过的某个点,这个命令就可
4、以 group 两个斜面间的区域 21、ansys 导入 flac 的时候,如何划分一下 group 答:在 ansys 中对不同 group 用不同的材料号(mat)定义,再用ansys2flac 之类的程序转换成 flac 的.dat 文件就可,group 自动生成。22、 fix x y z range x -0.1 0.1 a fix x range x -0.1 0.1b 请问 a 式与 b 式有何区别? 答:a 表示 在 x=0 这个平面 是固定的 b 表示在 x=0 这个平面 只 x 方向不能发生运动。 23、如何知道每个单元的 ID? 答:用鼠标双击单元的表面,可以知道单元的
5、ID 和坐标。 24、如何进行切片? 答:plo set plane ori (点坐标 ) norm (法向矢量) plo con sz plane (显示 z 方向应力的切片 ) 这两行命令,用来显示指定面上的位移变化云图(把第二行改改,还可以显示应力云图,等其他) ;plot con syy y 向应力云图plot con szz outline on z 向应力云图加显示网格plot con szz effective on z 向有效应力plot con szz inter 1e4 云图显示增量(后面的值越小,图上颜色条带就越细, )plot con szz min -10e4 显示最
6、小值,图上小于这个数的值都用同一色表示plot con szz max 10e4 显示最大值,图上大于这个数的值都用同一色表示plot fip outside on 显示力的箭头 现在看边界,如果是不动面,用 fix x range如果边界带速率,或位移用 iniset grav 0 0,-9.81ini dens 1000设置重力就一定要设置密度,不考虑重力就不要设计密度。两者共存的,设计程序的时候一定要记的。set mech force 50solve用 SET 命令设定了最大不平衡力的上限值为 50N,即:计算的最大不平衡力一旦小于此值,我们便人为模型达到了平衡状态。求解过程终止。设置这
7、个目的是在一定精度要求下,加快计算(值越大,算的量越少,花的时间也越少) 。25、如何在程序中进行暂停,并可恢复计算? 答:在命令中加入 pause 命令,用 continue 进行继续。在我们分步求解中想得到某一个过程中的结果,不用等到全求完,还可以在分布求解错误的时候就进行改正,而不是等到结果出来。26、如何跳过某个计算步? 答:在计算中按空格键跳过本次计算,自动进入下一步 27、Fish 是什么东西?Fish 是否一定要学?答:是 FLAC3D 的内置语言,可以用来进行参数化模型、完成命令本身不能进行的功能。Fish 可以不用学,需要的时候查 Mannual 获得需要的变量就可以了。 2
8、8、如何调用一些可选模块? 答:config dyn (fluid, creep, cppudm) 29、range 的命令: 取圆球面的范围,就是同一个中心,不同半径之间夹的那一部分,用在给球面施加力,约束时用该命令,命令:range annulus center xyz(中心坐标)radius r1 r2(两半径) ,想要表达的是球面时,r1r2 取值相差小点即可,如 radius 2.9 3.1 这时范围就是半径为 3 的面。圆柱的环面施加力或约束时也用到该命令。range cylinder end1 xyz end2 xyz radius r(end 是表示柱体断面的中心点坐标)ran
9、ge plane above 表示的是 z 大于 0 的上部结构range plane blow 表示的是 z 小于 0 的下部结构range dd a dip b origin xyz 表示的是走向角度 a,倾向角度b 的面dd 表示,当面与 xoy 面相交时的垂线与 y 轴正向的夹角,顺为正,说白话,就是面的指向方向面是朝那个方向的意思。dip 面与 xoy 面的夹角,与 z 轴负向为正,origin 表示面经过 xyz 点;range normal xyz 法线方向为 xyz 的面range not 没试用过(范围外的东西,想想用到圆上取范围非常好) 。30、如何调用 FISH 文件?
10、答:新手调用 FISH 文件出错的原因:是:文件名为 fishcall.fis.txt(因为在记事本里这个 txt 的后缀是默认的) 。修改方法是:将 fishcall.fis.txt 里面的内容全部复制到新的记事本里,再按保存,保存时候将保存类型改为所有文件(*,*),这时你便可在文件名栏内输入 fishcall.fis,类型也就正确了。 31、问 a 式与 b 式有何区别?答:fix x y z range x -0.1 0.1 a fix x range x -0.1 0.1b a 表示 在 x=0 这个平面 是固定的; b 表示在 x=0 这个平面 只 x 方向不能发生运动。32、如何
11、查看各个时段不平衡力的具体数值? 答:采用 his 来记录计算,包括位移应力等 命令 his unbal;his gp(zone) zdis range (0 0 0) 或者 id=?; 导出数据命令 his write n vs m begin 时步 end 时步+ file filename.his n 表示纪录的 id m 表示时步; 要导出不平衡力的具体数值 his unbal step 100000 or solve his write 1 vs step begin 1 end 1000 file 123.his 使用上述命令就可以查看各个时步下的不平衡力的具体数值。33、如何显示
12、变形轮廓线的命令? plo ske magf 10 其中 10 为放大系数34、怎样查看模型的矢量图?pl grid vel (disp) redplo blo gro 可以看到不同的 group 的颜色分布(看模型的样子) ;plo gpfix red sk 看模型的边界情况plo fap red sk 看模型的体力分布plo con dis (xdis, ydis, zdis) 看模型的云图:位移plo con sz (sy, sx, sxy, syz, sxz) 看模型的云图:应力plo dis (xdis, ydis, zdis) 看模型的矢量图plo info 看模型有多少单元、节点
13、plot grid 可以查看网格,plot grid num 可以查看节点号。plot fix 查看边界条件plot grid mark 查看预先留设的开挖边界pl syy (sxx xdis ydis sig1) fill;查看模型的云图。注意,带 fill 关键字就是查看云图,不带 fill 关键字就是查看等直线图。plo add ske 剖面的沉降云图上增加网格plo add dis plane 剖面上的矢量plo add axe 剖面上坐标值35、怎样查看 INTERFACE?答:pl grid iface id redprint iface id ; id 是 interface 的
14、编号36、FLAC 怎么调整视图?答:ctrl+z 或 shift+z ,然后点击要放大部位的左上角和右下角就可以放大某个部位,如果要缩小,双击屏幕。37、判断模型中网格点是一个还是两个的方法。答:显示模型网格点的 ID,如果有两个 ID 号,说明是两个。节理面上的网格点应该有两个 ID 号。38、flac 中影响运行时间的因素答:(1)单元材料,结构面之间的刚度对比很大;(2)划分单元尺寸对比很大,尽量避免长细比大于 5:1 的细长单元。39、简单命令解释:答:1.plot add cont disp out on 时时显示云图2.plot block group gpnum on 显示节点
15、编号3.plot inter red ske 看接触面的效果.4.del range group 1 not 不删除 group 信息 5.del range group 2 删除 group 信息6.ini z add * range group * (name) 7.interface 1 face range x(*,*)y(*,*)z (*,*)或者 interface 1 face range plane orig(x y z 坐标)norm (法线方向) 建立接触面8.water table clear 清除水位线9.pl set back white 设置背景10.set plo
16、t bitmap size (800,510) 设置图片大小11.plot set caption size 37 设置图例数字大小12.plot set center 100 0 40 设置图片中心13.fix x y z range x -0.1 0.1 表示在 x=0 这个平面是固定的14.fix x range x -0.1 0.1 表示这个平面只 x 方向不能发生运动40、在固定面上施加位移时,fix 与 apply 不能同时使用,应该用 ini和 fix,才能施加上。41、在进行三轴压缩时候,如何施加围压?答:使用 apply sxx =-1e5 szz=-1e5 range cy
17、l end1 0,0,0 end2 0,2,0 radius 1,其中 X 和 Z 轴是在我们眼睛正对的面。42、在进行三轴压缩时候,如何对轴向力进行位移控制?答:使用 ini yvel 1e-7 range y -0.1 0.1 ini yvel 1e-7 range y 1.9 2.1 即可,数值可以变化。43、在 FLAC3D 计算中,如何记录程序中的值,如果不想时时记录该怎么弄?答:用 hist wtite n1 vs n2 file.txt, 想间隔记录可用 set hist_rep n44、想记录接触面上的相关力学参数该怎么办?答:hist inter i nstress /ndi
18、s/sstress/sdisp45、接触面设置参数包括哪几个?答:coh dial fric kn ks ten.46、如何输出单元格的数据?答;expgrid .flacsd47、建模的过程中,或同一个文件中有重复应用的命令,就可以用宏命令了。答:宏命令: macro 名称 命令的内容例:macro shiping plot con zdip outline on用 shiping 代替后面的引号里的东西。在建模中只要有相同的都可以用 shiping 代替,缩减了书写命令的时间,整体上也能让人看的更明白。注意,你命名的符号最好跟 falc 自带的命令不要重复,处理方法是,名称后带数字。如 s
19、hiping1;如果你想修改宏对象,改变命令的内容,在名称上加单引号就是替换的意思。 例 macro shiping1 新的命令 。48、group 的用法:group 直接加到模型的后面相当于给整个模型命名;group name range x xyy zz 给一个范围的体命名 。例:在 gen zone brick - size - group name49、plot 中 add 与 set 的区别add 原来有的加新的用 addset 原来没有的,设置用 setplot create GravV plot set plane dip=90 dd=0 origin=3,4,0 plot s
20、et rot 15 0 20 plot set center 2.5 4.2 4.0 plot add bound behind plot add bcont szz plane plot add axes plot show 。我们想显示别的图像,又不想把这个弄好的图像丢掉,则用 plot create name 建立一个图像,在建立新的图像时,plot create name1(新的名)显示以前的图像,那么输入如下命令:plot current name。另一种方法,输入: plot show name 50、不同的区域定义不同的本构模型model mohr range z=0,5 pro
21、p ClayeyGravel range z=0,5 在定义参数后面加个 range 范围,为了避免在输入命令时总要通过range全局坐标域的方法来限定该命令执行的区域,我们可以先通过 range 命令来给不同的区域命名,然后就可以用 range区域名的方法来限定范围了。这样就避免了每次输入命令要限定在该区域执行时,都要输入一遍坐标范围的麻烦。具体如何给区域命名,以及日后如何用区域名来限定该区域,可以看下面这个例子: range name trench x=-1,1 y=0,4 z=-2,2 model null range trench51、右手法则:建立模型时,各点的排放时有规律的。方法如
22、下:拿出右手伸直了 4 指合并指向就是 p0 到 p1 的方向,掌心指向就是 p0 到 p2 的方向,拇指与 4 指垂直的方向就是 p0 到 p3的方向。很多人说网格后的云图有带齿什么的,改进的方法就是重复用下面的命令Gen ADJust 用此命令来微调网格, 使之离散平滑。该命令可以连续使用以增进平滑效果的。52、如果想看模型在计算过程中的变化情况Set movie avi step 1 file 1.avi;定义动画,动画记录的频率为 1 时步一频率(计算量大的时候就没必要 1 时步的记录,可以 10,100 时步的记录,减少电脑的使用):Movie start;动画开始附件里这个命令“h
23、ist write 5 vs 2 file dibiao.txt” 能将结果按时步讲就算结果输出 hist n= 5 :N(Nstep)=5 指定了每迭代计算 5 次记录一次相关的值(如果不指定,则系统默认值为 10) ;hist unbal 命令的第二行指定记录最大不平衡力(the maximum unbalanced force) ;hist gp zdisp 4,4,8 记录坐标值为(4,4,8)的节点 Z 方向的位移(这个用的多,查看收敛量时就是这个命令,上面两个个人认为用处不大,也很少用过) 。53、应变硬化/软化模型 SS,如何选择参数表Pro ftab 1 ctab 2 dtab
24、 3Table 1 0 45 .05 42 .1 40 1 40Table 2 0 2e6 .05 1e6 .1 5e5 1 5e5Table 3 0 10 .05 3 .1 0答:对于 SS 模型来说,材料需要硬化或者软化,也就是体现在内聚力和摩察角的不断变大或变小。table 1 就是在 0, 0.05, 0.1, 1 的应变/位移值上体现出内聚力的变化。table 2、3 同理。点的选取决定了应力应变曲线的形状,她的数量可以是任意的。应力与应变之间的关系塑性应变 plastic strain 与参数之间的关系,自己随便定义gen zone cyl p0 0 0 0 p1 1 0 0 p2
25、 0 4 0 p3 0 0 1 size 12 30 12gen zone reflect norm 1,0,0gen zone reflect norm 0,0,1model sspro den 2500 bulk 2e8 she 1e8 co 2e6 fric 45 ten 1e6 dil 10pro ftab 1 ctab 2 dtab 3table 1 0 45 .05 42 .1 40 1 40table 2 0 2e6 .05 1e6 .1 5e5 1 5e5table 3 0 10 .05 3 .1 0fix x y z range y -.1 .1fix x y z range
26、 y 3.9 4.1ini yvel 2.5e-5 range y -.1 .1 ;施加恒定速度测试反力(应变控制)ini yvel -2.5e-5 range y 3.9 4.1def ax_strstr = 0pnt = gp_headloop while pnt # nullif gp_ypos(pnt) 0.2def ShowPlasticZoneszp = zone_headloop while zp # nullif z_prop(zp,es_plastic) 0.2 ;塑性剪应变 es_plastic; 塑性拉应变 et_plasticz_group(zp) = yieldels
27、ez_group(zp) = otherendifzp = z_next(zp)endLoopendShowPlasticZonesplo crea qqqplo add surf red range group yieldplo add axes greenplo set rot 123 313 3plo set mag 1.5ret54、查看网格塑性区Pl block state55、如何建立马蹄形隧道newgen zone radcyl size 5 8 8 10 第一部分gen zone radcyl size 5 8 8 10 第二部分gen zone brick size 5 8 10 第三部分gen zone reflect orig 0 0 0 norm 1 0 0 ; 镜像plo crea qqq plo add surf red plo sho 56、如何显示锚杆Plot sel geom node on cid on