flac3d 基坑降水 流固耦合命令流

1.ÉøÁ÷ʾÀý;------------------------------²ÎÊý²¿·Öbulkmodulus,K390MPaÌå»ýÄ£Á¿£¨ÍÁÖʱȽÏÓ²£©shearmodulus,G280MPa¼ôÇÐÄ£Á¿soildrydensity,¦Ñd1200kg/m3ÍÁµÄÃܶÈwaterdensity,¦Ñw1000kg/m3Ë®µÄÃܶÈwalldensity,¦Ñwal1500kg/m3µµÍÁǽµÄÃܶÈpermeability,k10?12m2/Pa-sÉø͸ϵÊýporosity,n0.3¿×϶ÂÊfluidbulkmodulus,Kf2.0GPa;-------------------------------------;-----ÒÔÏÂÊÇÃüÁîÁ÷;------------------------------------------------------;Excavationinasaturatedsoil;------------------------------------------------------configfluid;;---geometricalmodel---½¨Á¢Ä£ÐÍgenzonebrickp11800p20180p30018size181818rat111groupsoil;ÉèÖÃÍÁ×égroupexcavaterangex04y04z05;ÉèÖÿªÍÚ·¶Î§×égroupwal1rangex45y05z07;groupwal2rangex04y45z07;groupwallrangegroupwal1anygroupwal2any;µµÍÁǽ1ºÍ2ºÏ²¢ÎªÒ»¸öµµÍÁǽ×é;---fluidflowmodel---modelfl_iso;¸÷ÏîͬÐÔÁ÷Ìåpropperm1e-12poro0.3;ÉèÖÃÉø͸ϵÊýºÍ£¿£¿±Èor£¿£¿ÂÊinifdensity1e3;³õʼ»¯Ë®ÃܶÈinifmod2.0e9ftens-1e-3;³õʼ»¯Ë®µÄÌå»ýÄ£Á¿¼°£¿£¿modelfl_nullrangegroupexcavate;¿ªÍÚ²¿·ÖûÓÐË®Á÷modelfl_nullrangegroupwall;µµÍÁǽ²»ÉøË®inipp0grad001e4;³õʼ»¯¿×϶ˮѹÁ¦fixpprangez-0.10.1;Ô¼ÊøÓ¦¸ÃÊÇÕû¸öÄ£Ð͵ĵײ¿°É£¿fixpprangex-0.14.1y-0.14.1z4.95.1;¿ªÍÚ²¿·ÖÒ²¸øÔ¼ÊøסÁË£¬¶øÆäËûµÄµØ·½²»¹Ü;---mechanicalmodel---Á¦Ñ§Ä£ÐÍmodelelas;µ¯ËÜÐÔpropbul3.9e6shea2.8e6;ÉèÖÃÌå»ýÄ£Á¿ºÍ¼ôÇÐÄ£Á¿£¬ÕâÁ½¸ö±ØÐ룬ÓÐÒ»´ÎÎÒÖ»ÉèÖÃÁËÌå»ýÄ£Á¿£¬¾Í´íµÄÒ»ËúºýÍ¿modelnullrangegroupexcavate;Á¦Ñ§Ä£Ð͵ĿªÍÚinidensity1.2e3;³õʼ»¯Ãܶȣ¬Ó¦¸ÃÊÇÍÁµÄÃܶÈinidensity1.5e3rangegroupwall;³õʼ»¯Ç½£¨µµÍÁǽ£©µÄÃܶȣ¬fixxrangex-.1.1;Ê©¼ÓÔ¼Êø£¬¾Ý˵²î·Ö·¨Ô¼Êø0µÄʱºò£¬±ØÐëÈ¡Ò»¸ö-0.1<0<0.1,·¶Î§fixxrangex11.912.1;Ê©¼ÓÔ¼Êø£¬»òÕß½Ð×ö±ß½çÌõ¼þfixyrangey-.1.1;Ê©¼ÓÔ¼Êøfixyrangey11.912.1fixzrangez11.912.1;initialtotalstresses³öÊÀ×ÜÓ¦Á¦iniszz0grad00-1.5e4;z·½ÏòµÄÓ¦Á¦´óһЩinisxx0grad00-1.2e4;ˮƽ·½ÏòµÄ£¬Ð¡Ò»Ð©inisyy0grad00-1.2e4;applynstress0grad00-1.2e4rangex0.04.0y3.94.1z0.05.0applynstress0grad00-1.2e4rangex3.94.1y0.04.0z0.05.0applynstress-7.5e4rangex0.04.0y0.04.0z4.95.1;---setting---ÉèÖÃsetgravity0010;ÉèÖÃÖØÁ¦¼ÓËÙ¶È;---initialstate---;³õʼ»¯×´Ì¬solveforce1;checkinitialequilibriumÇó½âƽºâ;---histories---¼Ç¼sethist_rep40;¼Ç¼²½Îª40histfltime;¼Ç¼ʱ¼ähistgppp006;¼Ç¼ijµãµÄ¿×϶ˮѹÁ¦histgpxdis400;ˮƽλÒÆhistgpxdis402;¼Ç¼ˮƽλÒÆhistgpxdis420histgpzdis005;¼Ç¼ÊúÏòλÒÆhistgpzdis205histgpzdis405histgpzdis225histgpzdis425histgpzdis445histgpzdis1001histgpzdis1002;;---excavation---¿ªÍÚsetfluidoff;¹Ø±ÕË®;applyporepressureatwallsÔÚǽÉÏÊ©¼Ó¿×϶ˮѹÁ¦applynstress0grad00-1.e4rangex0.04.0y3.94.1z0.05.0applynstress0grad00-1.e4rangex3.94.1y0.04.0z0.05.0applynstress-5.e4rangex0.04.0y0.04.0z4.95.1solve;force1;Çó½âsaveexc1.sav;±£´æ;;---drainage---applyremovenstress;³·µô¸Õ²ÅµÄѹÁ¦defrelaxsetup;¶¨ÒåÒ»¸öÊͷź¯ÊýµÄ²ÎÊýstep0=stependrelaxsetup;µ÷ÓÃdefrelax;¶¨ÒåÒ»¸öÊͷź¯Êýrstep=step-step0ifrstep<ncycthenrelax=1.0-(float(rstep)/float(ncyc))elserelax=0.0endifend;½áÊøsetncyc=1000;Ïȸ³Öµ£¬Ëæºóµ÷ÓÃapplynstress0grad00-1.e4histrelax&rangex0.04.0y3.94.1z0.05.0applynstress0grad00-1.e4histrelax&rangex3.94.1y0.04.0z0.05.0applynstress-5.e4histrelax&rangex0.04.0y0.04.0z4.95.1cycncycsolve;Çó½âsaveexc2.sav;±£´æ;;---percolation---Éø͸setfluidon;ÉèÖÃË®Ëã·¨fixpp0rangex-0.14.1y-0.14.1z4.95.1;Ê©¼Ó¿×϶ˮѹÁ¦cyc9000;Ñ­»·saveexc3.sav;±£´æplotcreateexcav;½¨Á¢Ò»¸öÏÔʾµÄÊÓͼһÏÂÄÚÈݶ¼ÊÇÏÔʾÄÚÈÝplotsetrot2000195plotsetcent666plotsetdist39.18plotsetmagn0.8plotsetplaneori000normal0-10plotaddcontppoutonplotaddblockgrouplgrarangegroupwallplotaddcontppint10000max110000rangex04y04z57plotaddflowplaneplotshowret