Excel学习―快速学习Excel插入图表

1、当执行这个FISH函数后，生成(nptab1+10,nprof组数据，这些数据为Aux2调用。注意：本段程序产生的数据是为那个特定的曲面边坡制作的，不是一个general purpose的LIB5）FISH函数aux2的作用是使用WATER TABLE face命令产生水面，把table 1的ytable值作为z坐标，aux1产生的数据作为x,y坐标。def aux2command  set grav 0 0 -10  water dens 1000end_commandloop n(1,nptab1-1  loop k (1,npr

2、of-1     xx1=xtable(n+10,k     yy1=ytable(n+10,k     zz1=ytable(1,n     xx2=xtable(n+10,k+1     yy2=ytable(n+10,k+1     zz2=zz1     xx3=xtable(n+11,k+1 

3、60;   yy3=ytable(n+11,k+1     zz3=ytable(1,n+1     xx4=xtable(n+11,k     yy4=ytable(n+11,k     zz4=zz3     command        WATER TABLE face xx1,yy1,zz1 xx2,yy2,z

4、z2 xx3,yy3,zz3        WATER TABLE face xx1,yy1,zz1 xx3,yy3,zz3 xx4,yy4,zz4      end_command   end_loopend_loopendaux224.产生结构元的几何形状(Geometry CreationFLAC3D有6种类型的结构单元，它们分别是梁单元beamSELs,锚杆单元cableSELs,桩单元pileSELs,壳单元shellSELs,地质网格单元geogridSELs和衬

5、砌单元linerSELs.每一种类型的结构由相应分量对象的集合来表示。例如锚杆结构(cable由一系列的cableSELs组成，而衬砌结构(liner则由一系列的linerSELs组成。每种结构元都由两个不同的辨识号码来表示，以便与其它结构元相区别。用Seltype命令和相应的关键字来定义结构单元:Seltype beam/cable/geogrid/liner/pile/shell,每一种支护类型的性质用SEL 命令的property关键字定义.History SEL keywordsID：用来标识一个结构对象，比如第一根锚杆的ID1，第二根锚杆的ID=2，它是对整个结构进行编号；CID:用

6、来表示一个结构对象内部的分量，比如一根锚杆划分为3段，则这根锚杆的CID分别是1,2,3,用来表示这根锚杆的不同部分。ID是为结构进行编号，CID是为结构的分量进行编号。比如下面的例子： sel pile id=1 begin=(10.0, 1.0, 0.0 end=(10.0, 1.0, -10.0 nseg=5sel pile id=2 begin=(10.0, 3.0, 0.0 end=(10.0, 3.0, -10.0 nseg=5在这个例子中，设置了2根桩，第一根桩以id=1来表示，第二根桩以id=2来表示；第一根桩分为5段，则cid=15，第二根桩也分为5段，则cid=610。 可

7、以用plot sel geom 命令来查看设置是否正确。 问题: 我设了三根桩，sel pile id=1 begin=(10.0, 1.0, 0.0 end=(10.0, 1.0, -3.0 nseg=3sel pile id=2 begin=(10.0, 3.0, 0.0 end=(10.0, 3.0, -3.0 nseg=3sel pile id=3 begin=(10.0, 5.0, 0.0 end=(10.0, 5.0, -3.0 nseg=3图中红色的数字1-3-4-2代表什么？是node？还是link？生成node的时候，如果周围存在zone，会自动生成link，而且link和n

8、ode位置总是重合的25. 基于ansys复杂地质体flac模型自动生成首先将复杂的三维地质模型建于ansys，然后导出结点坐标信息和单元信息，最后运用此程序就可以生成复杂地质体的flac命令，用flac程序运行此命令就可以生成复杂地质体的flac模型。复杂的三维地质模型一般是在cad中建立的,然后导入ansys中，通过此它建立复杂地形的flac模型。用autocad建好模型后，输出为.sat文件，就可以导入到ansys，不用软件转换。26. FLAC结果数据的导出pl set back white ;设置背景set plot bitmap size (800,510 ;设置图片大小plot

9、set caption size 37 ; 设置图例数字大小 (default=35 范围10-50mainwin size 1 1 position 0 0.8 ; command windowplot set window size 0.8 0.8 position 0.1 0.0 ; graphics windowplot set center 100 0 40 根据情况更改。27.一段fish命令的解释def aux1 定义函数loop n(1,nptab1 开始一个循环，从1循环到变量nptab1rr=xtable(1,n 从table 1 的x列中读第n个数付给变量rrzz=yta

10、ble(1,n 从table 1 的y列中读第n个数付给变量zzloop k(1,nprof 开始第二个循环，从1循环到变量nprofcase_of k 下面是一个选择结构，利用变量k判断alfa=0.5*(pi+.1*float(k-2/float(nprof-2xtable(n+10,k=rr*cos(alfaytable(n+10,k=rr*sin(alfacase 1 当k1的时候执行下面语句xtable(n+10,k=rr table n+10的x列第k个值等于变量rrytable(n+10,k=-40.0 table n+10的y列第k个值等于变量-40case 2 当k2的时候执

11、行下面语句xtable(n+10,k=rr table n+10的x列第k个值等于变量rrytable(n+10,k=0.0 table n+10的y列第k个值等于变量0end_case 结束选择结构end_loop   结束第二个循环end_loop 结束第一个循环end 函数定义结束aux1 执行什么定义的函数28. 如何删除指定区域的塑性区？可以通过分组塑性区来查看def plasticarrangezp=zone_headloop while zp # nullcurr_state=z_state(zp,0 if curr_state # 0 z_group(zp=

12、'yield'（这以后的什么意思？）else z_group(zp='other'endifzp=z_next(zpendloopendplasticarrangepr curr_state是将塑性区分组为yield，再用model null ran gr yield29.有关命令的解释1） apply nvel 0.1 plane dip 60 dd 270 range left_boun apply 施加边界条件，plane法向速度为0.1，plane 面的位置：倾角60（在xoy平面从-Z轴方向逆时针旋转），dd 走向270度，从+y方向绕z轴顺时针旋转；

13、2） ini szz -5e4 grad 0 0 -1e4ini 初始边界条件，szz=-5e4+x*0+y*0+-z*-1e4,注：-5e4是range范围起点坐标的SZZ值，x,y,z所施加应力网格点的坐标。30. 问：做的边坡开挖是一个四面体或是契体，除了在建模过程中单独建模再开挖以外，还有什么方法？单独建模的话，模型比较复杂，而且网格不太好统一。回答：用group定义,然后在挖掉即可,斜面你可以采用wedge来建立,或用brick建立,但是设置斜面要用:norm (X X X ori (X,X,X above x x1 x2 y y1 y2 z z1 z231. 渗流示例;-参数部分b

14、ulk modulus, K 390 MPa 体积模量（土质比较硬）shear modulus, G 280 MPa 剪切模量soil dry density, d 1200 kg/m3 土的密度water density, w 1000 kg/m3 水的密度wall density, wal 1500 kg/m3 档土墙的密度permeability, k 1012 m2/Pa-s 渗透系数porosity, n 0.3 孔隙率fluid bulk modulus, Kf 2.0 GPa 流体的体积模量;-;-以下是命令流;-; Excavation in a saturated soil;

15、-config fluid ;设置流体算法; - geometrical model - 建立模型gen zone brick p1 12 0 0 p2 0 12 0 p3 0 0 12 size 12 12 12 rat 1 1 1group soil ;设置土 组group excavate range x 0 4 y 0 4 z 0 5 ;设置开挖范围组group wal1 range x 4 5 y 0 5 z 0 7 ; 设置档土墙 1 组group wal2 range x 0 4 y 4 5 z 0 7 ;设置档土墙 2 组group wall range group wal1

16、any group wal2 any ;档土墙1 和2 合并为 一个档土墙组; - fluid flow model - 流体模型model fl_iso ;各项同性流体prop perm 1e-12 poro 0.3 ;设置渗透系数和？比 or ？率ini fdensity 1e3 ; 初始化 水密度ini fmod 2.0e9 ftens -1e-3 ; 初始化水的体积模量 及 ？model fl_null range group excavate ; 开挖部分没有水流model fl_null range group wall ;档土墙 不渗水ini pp 0 grad 0 0 1e4

17、;初始化孔隙水压力fix pp range z -0.1 0.1 ;约束 应该是整个模型的底部吧？fix pp range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ;开挖部分也给约束住了，而其他的地方不管; - mechanical model - 力学模型model elas ; 弹塑性prop bul 3.9e6 shea 2.8e6 ;设置体积模量和 剪切模量，这两个必须，有一次我只设置了体积模量，就错的一塌糊涂model null range group excavate ;力学模型的开挖ini density 1.2e3 ;初始化 密度 ，应该是土的密度ini

18、 density 1.5e3 range group wall ;初始化墙（档土墙）的密度，fix x range x -.1 .1 ;施加约束，据说差分法约束0 的时候，必须取一个-0.1<0<0.1,范围fix x range x 11.9 12.1 ;施加约束，或者叫做边界条件fix y range y -.1 .1 ;施加约束fix y range y 11.9 12.1fix z range z 11.9 12.1; initial total stresses 出世总应力ini szz 0 grad 0 0 -1.5e4 ;z方向的应力大一些ini sxx 0 grad

19、 0 0 -1.2e4 ;水平方向的，小一些ini syy 0 grad 0 0 -1.2e4 ;apply nstress 0 grad 0 0 -1.2e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0 apply nstress 0 grad 0 0 -1.2e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0apply nstress -7.5e4 range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1; - setting - 设置set gravity 0 0 10 ;设置重力加速度; - initial state

20、 - ;初始化状态solve force 1 ; check initial equilibrium 求解平衡; - histories - 记录set hist_rep 40 ;记录步为40hist fltime ;记录时间hist gp pp 0 0 6 ;记录某点的孔隙水压力hist gp xdis 4 0 0 ;水平位移hist gp xdis 4 0 2 ;记录水平位移hist gp xdis 4 2 0hist gp zdis 0 0 5 ;记录竖向位移hist gp zdis 2 0 5hist gp zdis 4 0 5hist gp zdis 2 2 5hist gp zdi

21、s 4 2 5hist gp zdis 4 4 5hist gp zdis 10 0 1hist gp zdis 10 0 2; - excavation - 开挖set fluid off ;关闭水; apply pore pressure at walls 在墙上施加孔隙水压力apply nstress 0 grad 0 0 -1.e4 range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0apply nstress 0 grad 0 0 -1.e4 range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0apply nstress -5.e4 range x

22、 0.0 4.0 y 0.0 4.0 z 4.9 5.1solve ;force 1 ;求解save exc1.sav ; 保存; - drainage - 排水apply remove nstress ;撤掉刚才的压力def relaxsetup ;定义一个 释放 函数的参数step0 = stependrelaxsetup ;调用def relax ;定义一个释放函数rstep = step - step0 if rstep < ncyc thenrelax=1.0-(float(rstep/float(ncycelserelax = 0.0endifend ;结束set ncyc

23、= 1000 ;先赋值，随后调用apply nstress 0 grad 0 0 -1.e4 hist relax &range x 0.0 4.0 y 3.9 4.1 z 0.0 5.0apply nstress 0 grad 0 0 -1.e4 hist relax &range x 3.9 4.1 y 0.0 4.0 z 0.0 5.0apply nstress -5.e4 hist relax &range x 0.0 4.0 y 0.0 4.0 z 4.9 5.1cyc ncycsolve ;求解save exc2.sav ;保存; - percolation

24、 - 渗透set fluid on ;设置水算法fix pp 0 range x -0.1 4.1 y -0.1 4.1 z 4.9 5.1 ;施加孔隙水压力cyc 9000 ;循环save exc3.sav ;保存plot create excav ;建立一个显示的 视图 一下内容都是显示内容plot set rot 200 0 195plot set cent 6 6 6plot set dist 39.18plot set magn 0.8plot set plane ori 0 0 0 normal 0 -1 0plot add cont pp out onplot add block

25、 group lgra range group wallplot add cont pp int 10000 max 110000 range x 0 4 y 0 4 z 5 7plot add flow planeplot showret32.VB中调用FLAC与FLAC调用VB程序在VB中调用Flac3d.exe的源代码：Private Sub Command1_Click('定义打开Flac时返回的值Dim Flac As Long'在VB中调用Flac3dFlac = Shell("E:工具箱FLAC 3Dflac3d v3.00-251Flac3d v3.0

26、0f3d300.exe", vbNormalFocus' 激活Flac,使其窗口获得焦点'AppActivate Flac'使用代码模拟键盘输入' 按下回车键给计算器SendKeys "Enter", True'SendKeys "gen zone brick", True'SendKeys "Enter", True'SendKeys "plot surface", True'SendKeys "Enter", True

27、'SendKeys "Enter", True'调用命令流SendKeys "call F:VBworktest.dat", TrueSendKeys "Enter", TrueEnd Subtest.dat：restore tcav3.savcall stable.datstable.dat：;显示支撑轴力得变化plot hist 7;excavate the fourth layermodel null ran group pcav4;支撑稳定性判断程序;读出数据，将其输到文件c:pandat1.txt中;read

28、 the force from flacdef stableloop t(1,200 ;每500步判断一次 ;计算200步command step 200endcommandcommand;删除旧有文件Sys del F:VBworkt1.txt;打开文件开关,并将支撑轴力录入到c:pandat1.txt中。set log onset logfile F:VBworkt1.txtprint sel beam forceset log off;调用vb程序Sys start /wait F:VBworkinputdata.exeendcommand;11-1a_size = 1IO_READ = 0IO_WRITE = 1IO_FISH = 0IO_ASCII =1 ;当需要读入，需使用ASCII点击“下一步”，如下图所示：第二步：1、添加系列1（11-1） 点击“系列”“添加”将光标移至按住;11-