- 3.66 MB
- 2022-04-29 14:28:47 发布
- 1、本文档共5页，可阅读全部内容。
- 2、本文档内容版权归属内容提供方，所产生的收益全部归内容提供方所有。如果您对本文有版权争议，可选择认领，认领后既往收益都归您。
- 3、本文档由用户上传，本站不保证质量和数量令人满意，可能有诸多瑕疵，付费之前，请仔细先通过免费阅读内容等途径辨别内容交易风险。如存在严重挂羊头卖狗肉之情形，可联系本站下载客服投诉处理。
- 文档侵权举报电话：19940600175。
'FLAC软件教程
第一讲FLAC简介
XDispx-位移YDispy-位移XVelX-速度YVely-速度还有几个可选用的关键词以协助实现INITIAL命令:Mark只有标记的结点值将被初始化
Unglued胶结命令GLUE将使交界面的上下两部分胶结在一起,使之不能滑动和分离,命令UNGLUE则为命令GLUE之逆。MArk<范围>在给定范围内的结点将被标记。计算的进行与结点的标记与否无关,但标记的结点可以限定一个区域作为命令INITIAL,
PROP和MODEL的作用域。命令GEN将自动对结点进行标记。<范围>的形式可以是I=i1,i2,J=j1,j2,i和j的顺序无关。Model关键词<范围>该命令对于给定的区域或范围赋予有关的本构模型。关键词有:
ANisotropic横观各向同性弹性模型Elastic弹性,各同性本构模型Mohr-Coulomb摩尔-库伦塑性模型Null零模型,用于开挖掉的单元SS应变软化Ubiquitous彻体节理模型NEW该命令可以不用退出FLAC重新开始一个新问题。
Plot关键词<关键词>...<范围>用缺省值。可以在一行上连写几个关键词以在一个图上画出几个变量。在绘图前先要用MODEL命令赋予网格以一定的本构模型。关键词及其意义为:BEam画出结构单元的几何形状Boundary画出网格的内外边界Cable画出锚束的几何形状
Disp位移矢量E_p塑性应变的轮廓线(只限于应变软化材料)Grid画出画出纪录在第nhis个历史上的变量值PP孔隙压力的等值线RF以矢量形式画出的固定结点的反力STAte画出单元中心当前的屈服状态
STress主应力矢量SXXxx-应力的等值线(全应力)SYYyy-应力的等值线(全应力)SXYxy-应力的等值线Velocity以箭头画出的速度矢量WAter水位线XDispx-位移的等值线YDispy-位移的等值线
XVelx-速度YVely-速度所谓开关其本身也是关键词,用于设置图形的的某些特征,有:color如红为Red,绿为GREen,黄为Yellow等。
Interval=c将等值线的间距置为c.在PLOT命令前发扫描线命令SCLIN可以画出等值线的值。Max=v在画矢量时,置箭头的最大长度为v.Noh不画标题Zero不画零值线
如果要放大或缩小图形,可以在PLOT命令前用WINDOW命令来设置窗口的大小。Print关键词<关键词>...<范围>和绘图命令PLOT一样,结点变量的值只有在给定的材料模型和性质后才可以打印出来。关键词有:
Apply打印出所施加力或压力的大小和范围Fix打印出固定x或y的结点Limits对SOLVE命令的限制值MArk打印出标记结点MEm占用的内存STruct打印出结构单元上有关结点的力,力矩和位移
Interface打印出交界面的数据,包括结点力和单位法向矢量主要的网格关键词Bulk体积模量COhesion粘结力Den质量密度Dilation剪胀角E_p全塑性应变(应变软化模型)
Friction摩擦系数SHear剪切模量XX-坐标YY-坐标PP孔隙压力JFric节理摩擦力JCoh节理粘结力JAngle节理角度
NUYxy-x泊松比NUZxz-x泊松比SIG1最大主应力SIG2最小主应力SXXXX-应力SYYYY-应力SXYXY-应力Theta最小主应力与x-轴的夹角
XDisX-位移YDisY-位移XModX-模量YModY-模量XVelX-速度YVelY-速度ASXXXX-应力(三角形A)BSXXXX-应力(三角形B)
CSXXXX-应力(三角形C)DSXXXX-应力(三角形D)ASYYYY-应力(三角形A)BSYYYY-应力(三角形B)CSYYYY-应力(三角形C)DSYYYY-应力(三角形D)ASXYXY-应力(三角形A)BSXYXY-应力(三角形B)
CSXYXY-应力(三角形C)DSXYXY-应力(三角形D)State塑性状态0弹性1正在屈服中2曾经屈服,现为弹性3已超过单轴抗张力4屈服并超过单轴抗张力
5已经超过抗张力6彻体节理正在屈服7彻体节理过去屈服过,现为弹性状态Tables打印所存贮的表格XReactionX-反力YReactionY-反力
PROp关键词=值<...><范围>本命令为MODEL命令赋予材料的性质。下面的关键词对PROPERTY命令加以补充:Regioni,j所有为标记结点所连续包围的区域将赋予该性质。i,j单元可为标记区域中的任一单元
Varxvyv性质可以有变差,xv和yv分别为x-方向和y-方向的变差值<范围>可为I=i1,i2,J=j1,j2,i和j孰者在前无妨,但<范围>必需在输入行的最后。各种模型所需输入的性质为:
弹性(1)剪切模量(2)体积模量(3)密度
摩尔-库伦(1)剪切模量(2)体积模量(3)密度(4)摩擦角(5)粘结力(6)剪胀角
横观各向同性(1)剪切模量(2)x-模量(3)y-模量(4)密度(5)NUYx(6)NUYz
彻体节理(1)剪切模量(2)体积模量(3)密度(4)粘结力(整体材料)(5)摩擦力(整体材料)(6)节理粘结力(7)节理摩擦力(8)节理角(9)整体材料的剪胀角
应变-硬化/软化(1)剪切模量(2)体积模量(3)密度(4)初始摩擦角(5)初始粘结力(6)摩擦角和塑性应变的关系表的表号(7)粘结力和塑性应变的关系表的表号(8)剪胀角和塑性应变的关系表的表号如果表号给0,则取关键词COHESION,DILATION或FRICTION所给出的值。
Quit该命令同STOP,FLAC终止运行。REstore文件名把以前用SAVE命令所存的文件读入内存恢复现场。RETurn该命令应是输入数据文件中的最后一个命令。SAve文件名将内存中问题的现场存入文件中,如果已经有同名的文件,则该文件将被复盖。
SClinnx1,y1x2,y2该命令在屏幕上作一个扫描线,以切割等值线,交点以A到Z的字母来表示。在一张图上最多可以有5条扫描线。扫描线的参数为:n=扫描线的编号(必需为1,2,3,4或5)x1,y1=扫描线的起点坐标x2,y2=扫描线的终点坐标
SEt<关键词>该命令设置整体条件,其关键词有:Aspecta调整屏幕的高宽比,以使输出不变形Colnn为打印输出的最大列数。用户必需在运行FLAC前设置好打印机的列数。
EGA设置640x350的EGA图形模式Forcef设置SOLVE命令时的不平衡力的限值。Gravityg重力加速度g及其与负y-轴的夹角Large大变形(坐标更新)
LogON在当前盘上打开名为FLAC.LOG的文件,以纪录FLAC的运行进程。如果FLAC.LOG文件已经存在,则将被复盖。OFF关闭FLAC.LOG文件。如果在稍后再SETLOGON,则屏幕显示将继续纪录在FLAC.LOG文件上。
Small小变形(坐标不更新)STepnSOLVE命令进行计算循环的限制数目。TimetSOLVE命令进行计算的限制时间的分钟数。
SOLve<关键词=值...>用本命令进行计算循环,关键词为:Step时步数目Time运行时间(分钟数)Force不平衡力缺省值为:S=500时步,T=5分钟,F=100
最长的时限为1440分钟(24小时),如果要求运算超过一天,则可以用几个SOLVE的命令。FLAC在运行中如要中止,可以按退出键Esc退出循环,将控制转给用户,进行交互操作。STEpn执行n个时步的循环运算。
STopFLAC终止运行。STRuct关键词...STRUCT命令用于定义结构单元的几何形状,性质等条件。STRUCT命令要求说明结构单元的类型和形状以及它与FLAC网格的连接情况。定义结构单元的关键词有:
(1)结构类型关键词梁一个梁单元系在关键词BEAM后用梁的几何形状和性质的关键词及其值来表征。锚束一个锚束单元系在关键词CABLE后以其几何形状和性质的关键词及其值来表征。
(2)结构单元的几何形状及其结点联结关键词BEGIN表征梁或锚束的起始点关键词END表征梁或锚束的终了点结构单元的端点可以在BEGIN和END关键词后直接用下面三个之一的关键词来表征:
Gridi,j结构单元的起点或终点为结点i,jNodenFLAC将自动按照构件单元输入的顺序为结构单元的结点编号,此处n为所编的第n号结点。x,y直接给端点以x和y坐标。
(3)结构单元的性质类型关键词Propn结构单元的类型要赋予性质数n,用structprop=n再加上关键词及其值来说明。(4)结构单元分段的关键词Segk将结构单元等分成k个分段。
(5)锚束预加张力的关键词TensiontSTRuctNode=n关键词每个结构结点用下列的关键词可以有不同的边界或几何条件。FFIx可以将n结点的x-和/或y-位移或转动固定。
FRee此命令为命令FIx之逆。Pi使结点n为铰接,即放松力矩。Loadfx,fy,m可对结点n施加x-和/或y-方向的力或力矩。x,yn结点的x-坐标和y-坐标。
STRuctProp=n关键词=值性质的关键词有:E=值I=值(梁的惯性矩)Area=断面值KBond=灰浆的粘结刚度SBond=灰浆的粘结强度Yield=锚束的屈服强度(力)
SYstem不退出FLAC暂时进入DOS系统,可以执行的DOS命令有:DIR,COPY,REN,DEL和TYPE.用回车键可以回到FLAC.Tablenx1y1建成一个点的坐标的表,以供FLAC引用。n是表号。
TItle在标题命令TITLE的下一行输入标题的名称。Unmark<范围>去掉在给定范围内结点的标记。WAter引入水面下的孔隙压力,孔隙压力为正。命令WATERTABLE=0将撤去所有结点的孔隙压力。
WIndow该命令按照xlo(x-低坐标),xhi(x-高坐标)和ylo(y-低坐标),yhi(y-高坐标)建立一个窗口,否则为自动建立(AUTO).
第二讲FLAC实例
FLAC基本单元
分析模式大小与RAM之关系
Flac3.3（解密）计算规模（不限）内存64MB弹性,各同性本构模型模型360万个单元摩尔-库伦塑性模型302万个单元
Flac3.3（解密）计算规模内存128MB弹性,各同性本构模型模型722万个单元摩尔-库伦塑性模型562万个单元内存180MB弹性,各同性本构模型模型1000万个单元摩尔-库伦塑性模型810万个单元
计算机计算规模的确定
flacw_sp180Gr1000,1000MePrmem
Flac3.3（解密）使用方法首次运行时，点击flac33.bat，按提示输入：96－05－10，就ok了
改批处理文件flac33.bat:@echooffsetbackdate5-10-96copyflaccfg.orgflac.cfgflac64setoday
Gr900,900MmPrmem
单位
正负号方向
(1)应力－正号代表张力，负号代表压力(2)剪应力详见下图，图中所示剪应力为正号(3)应变－正应变表示伸长，负应变代表压缩(4)剪应变－剪应变的正负号与剪应力相同(5)孔隙压力－孔隙压力永远为正(6)重力－正号的重力物质往下拉，负号的重力将物质往上提。
FLAC模型
1．空洞模式(nullmodel)使用于土壤被移除或开挖2．弹性模式3．塑性模式，包括a.Drucker－Pragermodelb.Mohr－Coulombmodelc.ubiquitous－jointmodeld.strain－hardening/softeningmodele.bilinearstrain－hardening/softeningmodelf.double－yieldmodelgmodifiedcam－claymodel
此外，另有选购(option)模式，包括：1.动力模式(DynamicOption)2.热力模式(ThermalOption)3.潜变模式(CreepOption)使用者另可使用FISH语言去建构独特需要的模型。
FLAC使用前准备步骤
步骤1：依比例画出所欲分析之资料于纸上画出地点之位置、地层资料、并简标示距离及深度资料。步骤2：换算输入资料成同一单位将现有地层资料，如Density,Bulkmodulus,Young`smodulus,tension,cohesion,frictionAngle等资料，换算成同一单位。附注：需谨慎检查输入资料之单位，如因单位不同而造成过大或过小的值，将会造成FLAC无法计算，而产生ERROR讯息。
步骤3：应用公式简略计算应用公式或依据经验，简略算出FLAC输出资料的范围。步骤4：建立x,y座标与nodei,j之关系于图上距离及深度之关系，建立x,y座标系统，再由x,y座标系统，转换与网格间系，为了便于以后输出资料的分析，故应确实掌握网格之位置及其相对应的x,y座标。建议在敏感区域使用较密之网格，其它地区则使用较疏之网格。
FLAC输入数据文件编写步骤FLAC只认文本文件可以用记事本、EDIT、PE2等编辑，但不能用word文件
FLAC数据文件编写顺序需依实际工程进行之逻辑步骤建立。如欲分析开挖题目，则需先求出未挖前之应力分布，再以此应力分布求开挖后之土壤变形及是否会崩毁，如因开挖后造成崩毁，则FLAC程式将无法继续，萤幕将出现ERROR讯息，如BadGeometryZone，---------"-----------．
FLAC程式前几行之顺序为oConfig ________oGrid ________oModel ________以后各行编写顺序，则建议以下之步骤：1.求起始之应力平衡(1)建立座标x,y与网格关系，使用Gen命令Genx1,y1x2,y2,x3,y3x4,y4i=i0,i1j=j0,j1FLAC可自动产生x,y座标与网格之关系，将网格依其疏密程度需要之不同，实际控制网格之座标。
(2)设定材料性质：prop(3)设定外力：SetGravity,ApplyPressure,inisxx,Syy(4)设定边界条件：fix,free(5)求起始之应力平衡：solve(6)储存：Save
２求工程之影响求出区域内之应力分布情况后，再依工程之流程及步骤阶段执行各工程进行过程之影响，建议使用以下之步骤：(1)调出起初之应力平衡状态：re_____.sav(2)设定新的材料性质：model,prop(3)设定新的支撑性质：struct(4)设定新的外力(5)设定边界条件(6)求工程时之应力平衡(7)储存
分析结果输出及图形绘制
分析成果绘图a.直接绘图（不存图档）flac：setplotwindowsflac：plotpenGRYD(欲画出格网及y向变位)说明:输入上述两行指令后，连结的打印机会直接印出图形，约占半页A4之纸张，如欲印出全页，则须设定打印机横向打印。b.绘图(存图档)flac：setoutYD.EMF(设以下要画之内容档名为.emf)
flac：setplotemfcolor(设彩色印制)flac：plotpenGRYD(画出之内容为格网及Y向变位)说明：YD.EMF可用Word软体叫出并绘图印出分析内容或成果数据flac：setlogYD.TXT(YD.TXT为所要列印内容之档名)flac：printydi=1,10j=1,5(设要印出y向位移量)flac：setlogoff说明：输入上述三个指令后，可用Word，Nodepad等软体调出YD.TXT并打印。
例一
;callex1.txtgrid6,5mmplotgridnumberplotgridgnum
;callex2.txt;FunctiontocalculateKandGfromEandnu.;Remembertosete_modandp_ratio!!!defDeriveKGs_mod=e_mod/(2.0*(1.0+p_ratio))b_mod=e_mod/(3.0*(1.0-2.0*p_ratio))end
;callex3.txtsete_mod=5e8p_ratio=0.25;callex4.txtDeriveKG;callex5.txtpropertydens=1000bulk=b_modshear=s_mod
;callex6.txt;sxx=0.5*syydefk0_setloopi(1,izones)loopj(1,jzones)sxx(i,j)=0.5*syy(x,y)end_loopend_loopEnd
例二
;ASimpleAnalysis--FailureofaSandySlope;calls1.txtgrid4020modelmohrgen-15,-16-15,2060,2060,-16;calls2.txttable1-15,00,034.641,2060,20gentable1
;calls3.txtmodelnullreg1,15;calls4.txttitleSlopeStabilityAnalysisplotgrid
;calls5.txtfixxi=1fixxyj=1fixxi=41propbulk=2e8she=1e8den=2500propfric=45coh=1e6ten=1e6setgrav=9.81histunbalhistxdi=25j=19histydi=25j=19solve
;calls6.txtpropcoh=0ten=0solvesavesavesl_ini.sav
;calls7.txt;Equation1:FoS=tan(actualfrictionangle)/tan(frictionangleatfailure)defssolvefric_act=friction(1,1)loopk(1,int_fric)fname="sl_"+string(k)+".sav"fric_fac=fric_inc*(float(k)-1.0)fric_cal=fric_ref-fric_facFoS=tan(fric_act*degrad)/tan(fric_cal*degrad)
commandinixdis=0ydis=0inixvel=0yvel=0printfric_calpropfric=fric_calsolvestep4000end_commandxtable(1,k)=fric_calytable(1,k)=xvel(25,19)xtable(2,k)=fric_calytable(2,k)=yvel(25,19)
commandsave@fnameend_command;stopseriesofrunsifunbalanceforce>100ifunbal>100thencommandprintfric_calFoSquitend_commandend_ifend_loopend
;calls8.txtsetfric_ref=31fric_inc=0.2int_fric=12ssolve;calls9.txtrestsl_9.savprinttable1printtable2window10,50-1040plotbounyvelzero
例题例3;Modelingcollapseofanunsupportedtrench;createmodelnewgrid5,5
;assignmaterialmodelandpropertiesmodelmohrpropbulk=1e8shear=.3e8fric=35propdens=1000coh=1e10ten=1e10;printandplotgridprintxy
plotgrid
;assignboundaryconditionsandapplygravityfixyj=1fixxi=1fixxi=6setgrav=9.81Plotgrfix
;monitorhistoryofmovementhisnstep=5hisydisi=2j=6;solvefortheinitialequilibriumstatesolve;plotandprinttheresultsplothis1
plothis1
titleasimpletrenchexcavationexamplesclin1(1,0)(1,5)plotholdsyyyelbougreplotholdsxxyelbougreprintsxxsyy
savetrench.sav;savetheinitialequilibriumstate;changethematerialstrengthandexcavatethetrenchpropcoh=0setlargeinitxdis=0ydis=0
modelnulli=3j=3,5
;calculatetheprogressivefailureofthetrench;andplottheresultspausestep100plotholdplasticboundary
plotholdxvyellint=5e-6disredmax=1e-2bougreen
step400 plotholdgrid
plotholdplasbou
plotholdsyyzeroint=2500dispmax=0.2magebougreen
nowwemakeamovieoftheprogressivefailure;startfromtheinitialequilibriumstate;andexcavatethetrenchasbeforepauseresttrench.sav
propcoh=0setlargeInixdis=0ydis=0modelnulli=3j=3,5;useaFISHfunctiontocapturethemovieplotsframesautomatically
defshowifn_step=0thenstep_int=25elsestep_int=n_stependifloopk(1,snap_shot)command
stepstep_intplotxdzint=.03max.25min-.25filldisredmax=0.5bouiwhendcommandendloopend;make30plotsata30stepintervalforthemovie
setsnap_shot=30n_step=30movieonpauseshowmovieoff;viewthemovie;pressthespacebartopausethemovie
;pressthespacebaragaintocontinuethemovie;presstheEnterkeytoexitthemoviepausemovieview;nowinstallhorizontalstrutstosupporttheexcavation
;startfromtheinitialequilibriumstate;andexcavatethetrenchasbeforeresttrench.savpropcoh=0setlargeinitxdis=0ydis=0histreset
hisxdisi=3j=4modelnulli=3j=3,5;installthreestrutsacrossthetrenchstructbeambegingrid3,4endgrid4,4seg=2structbeambegingrid3,5endgrid4,5seg=2
structbeambegingrid3,6endgrid4,6seg=2structprop1e=200e9i2.3e-5area4.8e-3step400;plottheresults:theexcavationisnowstable
plotstrdispblue25structbeam
plotholdhist1
分析范例－边坡稳定分析
分析网格之建立如下二图所示
01title02SLOPEUNDERGRAVITATIONALLOAD03grid20,1004;Mohr-Coulombmodel05mm06;soilproperties-notelargecohesiontoforceinitialelastic07;behaviorfordetermininginitialstresssate.Thiswillprevent08;slopefailurewheninitializingthegravitystresses09props=.3e8b=1e8d=1500fri=20coh=1e10ten=1e1010;warpgridtoformaslope:
11gen0,00,320,320,0j1412gensame9,1020,10samei621j41113marki=1,6j=414marki=6,j=4,1115modelnullregion1,1016;displacementboundaryconditions17fixxi=118fixxi=2119fixxyj=120;applygravity21setgrav=9.81
22;displacementhistoryofslope23hisydisi=10j=1024;solveforinitialgravitystresses25slove26;saveinitialstae27savesll.sav28;resetdisplacementcomponentstozero29inixdis=0ydis=030;setcohesionto031propcoh=0
32;uselargestrainlogic33setlarge34step20035plotholdbodisxvelmin-2.25e-4max0int2.5e-5zero36saves12.sav37step80038plotholdbodisxvelmin━1e-3max0int2e-4zero39saves13.sav
40restsll.sav41inixdis=0ydis=042propcoh=le4tens0.043setlarge44solve45savesl4.sav46inixdis=0.0ydis=0.047;installphreaticsurfaceinslope48watertable1den100049table1(0,5)(6.11,5)(20,9)
50defwetden51 loopi(1,izones)52 loopj(1,jzones)53 ifmode(i,j)>1 then54 xa=(x(i,j)+x(i+1,j)+(i+1,j+1)+x(i,j+1)55 xc=0.25*xa56 ya=(y(i,j)+y(i+1,j)+y(i+1,j+1)+y(i,j+1)57 yc=0.25*ya58 ifyc：strucnode1fixr定义loadfx,fy,m：strucnode4load0-1e60注意beggr511、begnode1、begx,y区别
例题14（liner）gr1515mmgen000303030300gen1414146060606014rat1.21.2i816j816gen14-301414601460-30rat1.2.833i816j18gen-30-30-3014141414-30rat.833.833i18j18
gen-3014-306014601414rat.8331.2i18j816gencircle14144genadjustprops5.75e9b6.6e9d2000coh1e7fr35fixxyi1fixxyi16fixxyj1fixxyj16inisxx-60e6syy-30e6solve
mnreg88strucprop1e13.8e9i2.8e-4a.15strucbeambeggr79endgr89strucbeambeggr89endgr99strucbeambeggr99endgr98strucbeambeggr98endgr97strucbeambeggr97endgr87strucbeambeggr87endgr77strucbeambeggr77endgr78strucbeambeggr78endgr79return
分析与思考复杂网格的建立巷道的支护gen000303030300可以不要后面的gen可改变前面gen定义的大小。思考：弱面的设置、锚网支护
第四讲FLAC实例
例题（x）newtitleHoleinMohr-Coulombmedium(dilatancy=0)-iebboundary-76x76elementsg7676mmgen-2.5-2.5-2.52.52.52.52.5-2.5gencirc0,0,2.5
mnmmreg38,38gencirc0,0,1.0mnreg38,38propshear=2.8e9bulk=3.9e9dens=2500coh=3.45e6fric=30dil=0ten=1e10inisxx=-30e6syy=-30e6szz=-30e6;iebxcycrbulkshearieb0,0,2.53.9e92.8e9
hisyvi10j15hisydi10j15hisxdi4j10step20000savex.savret
分析与思考ieb0,0,2.53.9e92.8e9相当于外圆边界固定命令格式：iebxcycrbulkshear思考：实现外圆边界固定是否还有其它办法？
例题（x1）grid100100mmohrpropd2000shear.25e8bu.5e8propfric40rdev=2tens=0dil=0gen-50-50-5050505050-50gencirc0050
mnreg1,1mnreg1,100mnreg100,1mnreg100,100fixxymarkinixv-.25e-3var.5e-30iniyv.25e-3var0-.5e-3inisxx-1e5syy-1e5szz-1e5wind-5555-5555step1000plotholdboussrint1.0e-5ret
分析与思考gencirc0050与fixxymark联合实现了外圆边界固定inixv-.25e-3var.5e-30表示xv在x方向线性变化，命令格式为：inixv=vx0varxy其含义为在j=1时其值由vx0线性变化为vx0+x在i=1时其值由vx0线性变化为vx0+y思考iniyv.25e-3var0-.5e-3的含义是什么？
例题（x2）*createawaveysurfacetostudyinitialstressconditionsg2420megenarc0,173,1790genarc6,173,17180genarc12,1715,17180genarc18,1715,17180genarc24,1724,2090
mnreg=6,20mnreg=18,20fixxi=1fixxi=25fixyj=1propd2000sh1e9bulk2e9setgrav=10*sethorizontalstress=2*verticalinisyy-3.4e5var=0,3.4e5inisxx-6.8e5var=0,6.8e5
*notein-situstressesgotozeroatmeandurface*----solvef=1e2,step=1000sclin15,05,20sclin212.5,012.5,20plotholdgridplotholdbousyyint2.5e4plotholdbousxxint5.0e4ret
分析与思考波浪形边界的形成
例题（wall）titleTESTWALL--;establishproblemgeometrygrid35,35modelmohrgen-3068-301035103568gen-5071-50103-20103-2071ratio.91i=1,11j=4,36gen-5065-5071-2071-2065ratio.9.9i=1,11j=1,4
gen-2065-2071071065ratio1.9i=11,31j=1,4gen06507110711065ratio1.2.9i=31,36j=1,4gen0710103101031071ratio1.21i=31,36j=4,36table1-60.0,100-44.8,100-40.3,99-36.6,98.5table1-33.6,99-30.6,100.2-27.3,102.1-25.2,102.8-24,103gentable1modnullregi=2j=35
inix-29y100.9i7j34inix-27y102.1i8j35;boundaryconditionsforconsolidationfixxi1fixxi36fixyj1;initialstressesassumecoefficientofearthpressureatrest=0.45inisyy-4446var04446inisxx-2000var02000iniszz-2000var02000
;conditionsduringexecutionsetlargesetgrav32.2;propertiesformediumdensewellgradedsandpropdens3.63fric36coh0dil7.5bulk1.33e6shear0.8e6;displacementmonitoringpointshistn100histunbalhistxdispi31j36histydispi31j36
histxdispi31j30histydispi31j30histxdispi31j24histydispi31j24histxdispi31j18histydispi31j18histxdispi31j12histydispi31j12histxdispi31j6histydispi31j6
;consolidatesolvefor=10;saveconsolidatedstatesaven0.sav;re-initializesomeparameterspriortostartingprobleminixdisp=0iniydisp=0inixvel0i136j136iniyvel0i136j136
fixxyi1fixxyj1fixxyi36;;STAGEI;modnulli3135j3135titleEXCAVATION--STAGEI
;shotcretepropertiesstructprop1e=5e8i=3e-3a=.333structbeambeggrid3136endgrid3135seg1prop1structbeambeggrid3135endgrid3134seg1prop1structbeambeggrid3134endgrid3133seg1prop1structbeambeggrid3133endgrid3132seg1prop1structbeambeggrid3132endgrid3131seg1prop1
;nailpropertiestoprowstructprop2a8.5e-3e9.3e8y=16360kbond=1.4e7sbond220struccablebegin-11.9,98.65end-0.1,100.75seg6prop2solvefor=10savenI.sav;;STAGEII;titleEXCAVATION--STAGEII
modnuli3135j2930structbeambeggrid3131endgrid3130seg1prop1structbeambeggrid3130endgrid3129seg1prop1;nailproperties--row2structprop3a8.5e-3e1.2e9y21040kbond=1.8e7sbond1430struccablebegin-11.9,94.15end-0.1,96.25seg6prop3solvefor=10savenII.sav
;STAGEIIItitleEXCAVATION--STAGEIIImodelnulli3135j2428structbeambeggrid3129endgrid3128seg1prop1structbeambeggrid3128endgrid3127seg1prop1structbeambeggrid3127endgrid3126seg1prop1structbeambeggrid3126endgrid3125seg1prop1
structbeambeggrid3125endgrid3124seg1prop1;nails--row3struccablebegin-11.9,89.65end-0.1,91.75seg6prop4struprop4a8.5e-3e4.65e8kbond7e6sbond667y8180;tieback-anchorstruccablebegin-44.3,83.9end-20.788.1seg6prop5struprop5a.0103e4.65e8y24750kbond7e6sbond1000
;tieback-ungroutedportionstruccablebeggrid3125endnode41seg1prop6ten20689struprop6a.0103e4.65e8y24750kbond1sbond1solvefor=10saveniii.sav;STAGEIVtitleEXCAVATION--STAGEIVmodelnulli3135j2023
structbeambeggrid3124endgrid3123seg1prop1structbeambeggrid3123endgrid3122seg1prop1structbeambeggrid3122endgrid3121seg1prop1structbeambeggrid3121endgrid3120seg1prop1;nailproperties-row4to7structprop7a8.5e-3e9.3e8sbond1333yi16360kbond1.4e7
structcablebeg-11.985.17end-0.187.25seg6prop7solvefor=10saveniv.sav;;STAGEV;titleEXCAVATION--STAGEVmodnulli3135j1519
structbeambeggrid3120endgrid3119seg1prop1structbeambeggrid3119endgrid3118seg1prop1structbeambeggrid3118endgrid3117seg1prop1structbeambeggrid3117endgrid3116seg1prop1structbeambeggrid3116endgrid3115seg1prop1structcablebeg-11.980.67end-0.182.75seg6prop7
solvefor=10savenv.sav;STAGEVItitleEXCAVATION--STAGEVImodnulli3135j1114structbeambeggrid3115endgrid3114seg1prop1structbeambeggrid3114endgrid3113seg1prop1structbeambeggrid3113endgrid3112seg1prop1
structbeambeggrid3112endgrid3111seg1prop1structcablebeg-11.976.17end-0.178.25seg6prop7solvefor=10savenvi.sav;STAGEVIItitleEXCAVATION--STAGEVIImodnulli3135j610structbeambeggrid3111endgrid3110seg1prop1
structbeambeggrid3110endgrid319seg1prop1structbeambeggrid319endgrid318seg1prop1structbeambeggrid318endgrid317seg1prop1structbeambeggrid317endgrid316seg1prop1structcablebeg-11.971.67end-0.173.75seg6prop7solvefor=10savenvii.sav
;STAGEVIII;titleEXCAVATION--STAGEVIIImodnulli3135j5structbeambeggrid316endgrid315seg1prop1solvefor=10savenviii.savret
例题（a1a）目的：模拟复杂情况下的煤层开采时的应力、变形。一、计算模型取计算模型大小：长100米，高60米。计算模型岩性：十层岩层（包括煤层），九层弱面。岩层厚度、力学性能与柱状图相同。计算模型网格大小：0.5米0.5米。计算规模：2.5９３万节点，2.56万个单元。十层岩层（包括煤层），九层弱面，如图
计算模型网格
边界条件
煤层开采前的应力ｙ分布
煤层开采前的应力Ｘ分布
煤层开采前的应力ｚ分布
二、计算方案分层开采煤层位于计算模型左侧的五十米，综采的煤层位于计算模型右侧的五十米。1．由于煤层上分层早已开采，而且时间达一年，所以先模拟上分层开采。2．模拟上分层开采稳定后，模拟下分层开采对应力、变形影响。3．模拟上分层开采稳定后，模拟综采对应力、变形影响。
三、计算结果分析1．上分层开挖
上分层开挖后的变形
上分层开挖后的应力y在上分层右下角应力集中最明显，压应力y由开挖前13.08MPa，开挖后升高为79.98MPa，是开挖前的6.115倍。
2．下分层开挖后的变形
下分层开挖后的应力y
在上分层右上角应力集中最明显，压应力y由开挖前13.08MPa，开挖后升高为74.93MPa，是开挖前的5.729倍。下分层右下角应力集中，压应力y由开挖前13.08MPa，开挖后升高为64.28MPa，是开挖前的4.914倍。在上分层开挖后压应力y为79.98MPa，下分层后降为64.28MPa。
3．模拟综采对应力、变形影响综采区开挖（留八米煤柱无支架）后模拟变形：
应力y压应力y最大值在上分层边界右侧2—3米处，应力值高达111.5MPa。
从上分层边界至综采区留八米煤柱（无支架）的附近局部图可以看出，煤柱已有较大区域破坏。
从上分层边界至综采区留八米煤柱（无支架）的附近局部图可以看出，煤柱有较大的侧向位移
从上分层边界至综采区留八米煤柱（无支架）的附近局部图可以看出，煤柱已有较大区域破坏，但仍有较大的速度，说明变形还将继续，即没有稳定。综上所述，上分层边界至综采区留八米煤柱（无支架），将导致煤柱整体破坏。要保证开采安全，必须采取支护措施。
4．模拟综采对应力、变形影响（1）综采区开挖（综采区开挖下半部分留三米煤柱，综采区开挖上半部分留九米加支架）后模拟
变形图
应力y图压应力y最大值发生在综采区边界左侧1.0—1.5米处，应力值高达117.5MPa。
局部应力y图
局部变形图
局部位移图
电动吸引器洗胃护理教学查房
洗胃目的清除胃内毒物和刺激物,挽救病人的生命减轻胃黏膜水肿为某些手术和检查做准备。
评估1、病人的病情、意识状态、生命体征、瞳孔、合作的程度、及心理健康状况。 2、病人中毒的情况，中毒的时间和途径、毒物的种类、性质、有无腐蚀性。 3、病人的疾病史，是否患有肝硬化伴食道静脉曲张、近期曾发生过上消化道出血或胃穿孔。 4、病人的口鼻腔黏膜的情况，有无义齿和其他疾患。
提问：为什么要对病人进行评估？答：可以了解病人中毒的程度：毒物的性质；以便采取相应的抢救程序和选择相应的洗胃溶液。对病人健康史的评估，可以判断病人有无洗胃的禁忌症，和预防并发症的发生。
准备1、护士的准备：核对医嘱，穿戴整齐、洗手、戴口罩。 2、用物的准备：自动洗胃机一台（处于功能位）治疗巾内放置洗胃管、镊子、50ml针管、纱布、石蜡油、水温计。治疗巾外放橡皮单、手套、弯盘、治疗碗内盛漱口水、标本容器、手消毒液、记录本、笔、水桶（桶内盛洗胃液5000—10000ml温度25—38度为宜）、污水桶、必要时准备开口器、舌钳、压舌板、电动吸引器一台。
操作步骤1、根据病人的情况、洗胃的目的、或者毒物的性质选择、准备洗胃的溶液：根据毒物的性质选择相对抗的洗胃溶液，当毒物性质不明时选用温开水或者生理盐水洗胃最适宜。可避免配制洗胃液而延误抢救时机，并可防止洗胃液选择不当而导致不良后果。
提问：1）、乐果、1605、1059等禁忌用什么洗胃溶液洗胃？ 答：高锰酸钾。 （2）、敌百虫禁用什么洗胃溶液？为什么？ 答：禁用碱性溶液，如2%碳酸氢钠溶液，因为遇碱性药物可分解出毒性更强的敌敌畏。
操作步骤2、把准备好的用物推至床旁，认真查对病人姓名，向别人解释洗胃的目的和方法。要点说明：对清醒的病人要消除病人的紧张情绪，使病人更好的配合操作；对自服毒物的病人应该注意病人的心理变化，耐心而有效的作好心理护理。
操作步骤3、连接好电动洗胃机的管道，打开电源开关和吸引器装置，处于备用状态。 4、协助病人取左侧卧位，昏迷者取去枕平卧位，头偏向一侧，把橡胶单围于病人的胸前，避免水和呕吐物弄湿、弄脏病人的衣服。如果病人有义齿，应先取下，弯盘放在病人的口鼻处。要点说明：注意保护病人，防止意外发生：避免义齿滑脱进入食道或气道内，造成对病人的伤害；烦躁的病人防止发生摔伤。
操作步骤5、戴好手套，取出胃管，测量应该插入的长度，从病人的发际到剑突的距离，一般为45—55厘米（小儿从到耳垂到鼻尖到剑突）。提问：如果胃管过长会导致怎样的结果？答：胃管插入过长会使液体直接进入十二指肠，不能吸出洗胃溶液，造成腹胀。也可以使胃管在胃内盘曲，使液体的进入不畅，影响洗胃效果。
操作步骤6、将润滑油倒在纱布上，润滑胃管的前端，放置好口圈，从口中缓缓插入，当胃管插入约10-15厘米的时候，胃管前端到达咽喉部，这时请病人深呼吸或做吞咽动作，昏迷病人当胃管插至咽喉部时可抬起病人头部，使其下颌靠近胸骨柄，加大咽部通道弧度，将胃管徐徐插入。这样可减轻插管时的阻力，减轻病人的不适。
操作步骤要点说明：插管时动作要轻、稳、勿损伤食道粘膜，尽量减少对病人的刺激；插管过程中如果病人出现呛咳、呼吸困难、(提问）应该立即拔除胃管，待病人休息后再插，避免误入气管；为昏迷病人插管采取平卧位，头偏向一则，用开口器从臼齿处打开病人的口腔，垫口圈如果舌后坠，用舌钳将舌拉出，再插管；如果有气道阻塞则先气管插管再插胃管。证实胃管进入胃内后用胶布固定好胃管。
提问：如何证实胃管进入胃内？回答： （1）注射器能够抽出胃内容物；或者进入胃内后有胃内容物从胃管内流出。 （2）用注射器抽吸10毫升空气迅速从胃管末端注入，用听诊器在胃部能够听到气过水声 （3）把胃管放入装水的治疗碗里面，如果没有气泡冒出，说明胃管在胃内
操作步骤7、遵医嘱留取毒物标本送检，抽净胃内容物，与洗胃机连接，调节参数，每次灌洗量300—500ml，反复冲洗直至洗出液澄清无味为止。 要点说明：要注意观察病人的病情变化（提问洗胃过程中观察内容），特别是血压、脉搏、呼吸、意识、表情、瞳孔、有无腹痛、腹胀及洗出胃液的颜色。每次灌入量和吸出量应该保持基本相等，以免造成胃储留；吸出的胃内容物为血性，或病人出现虚脱现象等异常情况时应停止洗胃，通知医生，采取相应的急救措施。并及时详细准确地做好抢救记录。
操作步骤8、准备拔管，洗胃完毕，分离胃管与洗胃机的接头，让胃内液体从胃管内留出完毕，先反折胃管，再用纱布包住胃管迅速拔出放在弯盘内。协助患者漱口，用纱布擦净病人的口角及面部，，取下橡胶单，脱手套，安置病人，整理用物做好终末消毒，清洗各种管道，及储液瓶，放入含有效氯500mg/l消毒液中浸泡2小时然后用蒸馏水冲洗干净待用。提问：为什么拔管时要反折胃管？（纪元） 答：避免胃管内的液体误入气管内，导致吸入性肺炎或窒息。
操作步骤9、消毒手，记录洗胃液的量及名称，呕吐物的颜色和气味，病人的情况。做好仪器使用维护登记。提问：如果幽门梗阻的病人洗胃，可在什么时间进行？ 答：饭后4-6小时或者空腹进行，并记录储留量，了解梗阻的情况。
提问：洗胃的注意事项有哪些？1）、急性中毒的病人应立即催吐或洗胃以减轻毒物的吸收，洗胃时行左侧卧位，插管动作轻柔勿损伤食道粘膜，中毒原因不明者，先抽出胃内容物送检，可先注入温开水或生理盐水。（2）在洗胃过程中要注意观察病人的病情变化，特别是血压、脉搏、呼吸、意识和洗出液的性质|；每次灌入量和吸出量应该保持基本相等，以免造成胃储留。如果发现异常，应停止洗胃，通知医生采取相应的急救措施并做好记录。 （3）、吞服强酸、强碱等腐蚀性的药物，禁忌洗胃。以免造成穿孔，给予牛奶、豆浆、蛋清、米汤等。以保护胃粘膜。
提问：洗胃的注意事项有哪些？（4）、消化道溃疡、食管阻塞、食管静脉曲张、胃癌等一般不洗胃。昏迷病人宜谨慎，且应去枕平卧，头偏向一侧，幽门梗阻者宜在饭后4-6小时洗胃，且应记录潴留量。 （5）、每次灌入量以300—500ml为宜，如灌入量过多，液体可从鼻腔内涌出而引起窒息，同时还易产生急性胃扩张，使胃内压上升，增加毒物吸收，突然的胃扩张又易兴奋迷走神经，引起反射性心跳骤停，对心肺疾病患者更应慎重。 （6）、用电动吸引器或自动洗胃机洗胃，应先检查机器功能是否正常，以确保安全。
谢谢大家！'
您可能关注的文档
- 最新ERP实施与应用PPT课件PPT课件.ppt
- 最新ERP入门基本知识课件PPT.ppt
- 最新ERwin使用教程课件PPT.ppt
- 最新EXCEL常用图表课件PPT.ppt
- 最新Eye & Vision Screening：眼睛和视力筛查课件PPT.ppt
- 最新F-信息技术促进了交流课件PPT.ppt
- 最新fact-or-opinion从阅读到写作PPT课件PPT.ppt
- 最新FFC vs VC PRODUCTS FFC(纯耐火土)和VC(玻璃瓷)产品课件PPT.ppt
- 最新flash-第三四课课件PPT.ppt
- 最新Flash工作界面介绍及制作小球课件PPT.ppt
- 最新Flash动画制作之—制作形状补间动画 05课件PPT.ppt
- 最新Flash动画设计与制作(9)导出和发布动画课件PPT.ppt
- 最新flashcs6基础教程之FlashCS6简介课件PPT.ppt
- 最新fragstats-操作课件PPT.ppt
- 最新Freemind制作思维导图教程84993788课件PPT.ppt
- 最新gambit例子课件PPT.ppt
- 最新GB预包装食品营养标签通则月-(2)课件PPT.ppt
- 最新GE-CT原理机构与维护__training课件PPT.ppt
相关文档
最近下载
微信/支付宝扫码支付下载-
关注微信公众号售出明细实时看