FLAC3D高手总结.doc

?;模型镜像 gen zone radcylinder size 25 1 25 25gen zone reflect normal -1 0 0 origin x y z(面上一点）;沿X轴镜像,通过对称平面法线向量确定对称面gen zone reflect normal 0 0 -1 ;沿z轴镜像 ;绘图控制pl add contour szz outline on ;在模型中显示位移-应变曲线hist gp ydisp 0,0,0hist zone syy 0,1,0hist zone syy 1,1,0pl his -2 -3 vs 1 ;在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负); -表示对其值作mirror ;对模型进行压缩实验的方法;即在模型两侧施加相反方向的速度ini yvel 1e-7 range y -.1 .1ini yvel -1e-7 range y 1.9 2.1 ;修改模型的坐标值ini x add -100 y add -100 z add -100 ;显示云图的同时也显示模型网格轮廓plot add cont disp outline on gradient/average;gradient更精确 ;输入角度、弧度方法 pi=π，90°为90.0*degraddef set_valsptA = 25.0 * sin(pi/2);ptA=25.0ptB = 25.0 * cos( 60.0*degrad );ptB=12.5ptC=pi;ptC=3.1415926endset_valsprint ptA ptB ptC ;施加结构单元方法sel shell id=5 range cylinder end1=(0.0, 0.0,0.0) end2=(0.0,25.0,0.0) radius=24.5 notplot add sel geom black black cid on scale=0.03sel node init zpos add -25.0 ;如何显示某一平面plot create name_planeplot set plane? origin 3 4 0 normal 1 0 0plot add cont disp plane behind shade onplot add sel geom blackplot add axes redplot current name_planeplot show ;如何施加法向应力并显示所施加的力apply nstress -3.0E7 range cylinder end1 0 -1 0 end2010 radius 39.9 not;柱体外表面施加法向力apply nstress -5.0E6 range cylinder end1 0 -1 0 end2010 radius 2.1;柱体内表面施加法向力plot grid green fap black;显示施加的力 ;获取某网格点的ID，坐标值def aaget_id=gp_id(gp_near(3,0,0))endaaprint get_idprint gp position range id 3 ;保存命令流文件set log onset logfile hua.log;文件存放的路径与调用的文本路径相同......set log off ;fish函数的说明与解释def find_add? ;定义fish函数find_addhead = null?? ;给head赋值p_gp = gp_head? ;第一个网格结点的指针赋给p_gploop_while p_gp # null? ;当p_gp值不为null时作循环x_pos = gp_xpos(p_gp)? ;将指针为p_gp的结点的x坐标值赋给x_posif x_pos = 10.0 then? ;如果x_pos = 10.0 则(执行)new = get_mem(2)? ;从主内存空间里得到2个fish变量对象并返回第一个对象的开始地址mem(new) = head? ;将head类型和数值置于地址为new的fish变量mem(new+1)=p_gp? ;第二个变量存储符合条件的zone地址head = new?????? ;将new值赋给headendifp_gp = gp_next(p_gp)? ;将结点指针为p_gp的下一个结点的指针赋给p_gpendloop;结束循环end? ;结束fish函数;实际上这个fish函数为满足条件（x坐标为10的）的zone的地址