第2425集流固耦合的理论与计算上甲丙dxf建立三维潜水面

|（E）Section24-25：Principleandcalculationoffluid-solidcouplingDXFfileforestablishingadivingsurfacein©：都 水粒气水粒气principleprincipleofeffectivePorewater

EffectiveEffectiveBiot系数定义为当孔隙Biot系数定义为当孔隙BiotBiotprincipleofeffective

PorewaterEffectiveEffectivewaterpressure

modelconfigurefluidmodelfluidactiveon

neor学 回路

环 环路In的区域用流体的null材料定义（与Null本构模型区别）。In的区域用流体的null材料定义（与Null本构模型区别）。GridNotConfiguredforFluid没有打开没有打开Modelconfigurefluid或者：modelfluidactive

zonezonegridpointinitializepore-pressureS0gradient（Gz*z）rangeposition-Z-orwiththezonewater ne(orzonewatersetcommand.此时必须设置流体密度zonewaterdensity与重力值modelgravity。打开Modelconfigurefluid或者：modelfluidactive打开Modelconfigurefluid或者：modelfluidactive

配流体模型zonefluidcmodelassignisotropic/anisotropicrange….其中若想设置隔水边界模型，用zonefluidcmodelassignnullrange…（注意与null本构模型的区别） Zonewater命令(orzonegridpointZonewater命令(orzonegridpointinitializepore-pressure) 如果FLAC3D打开渗流计算板块(即设置了modelconfigurefluid),则固体饱和密度设置：zoneinitializedensityf固体干密度设置：zonePropertydensity流体密度设置（全局设置）：zonewater流体密度范围设置：zoneinitializefluid-densityBiotCoefficient&Biot

Tip:BiotTip:Biot系数是zone属性，见Page10；用zonefluidpropertybiotfbiot模量是节点变量，由zonegridpointinitializefluid-moduluszonefluidbioton 量导致更快的收敛速度。无论如何，从数值计算的观点，没有必要大于20倍的土体的体积模量。zonefluidpropertyporosityfrangezonezonegridpointinitializesaturationfrange...

细粒土中，孔隙水可以承受明显的拉力（负超孔隙水压力）。在FLAC3D中，流出，土体不再饱和）。负超孔隙水压的极限用INITIALftens命令设置，程 zonegridpointinitializefluid-tensionfrange...fluid-modulus使用：zone pore-pressurerange…进行设孔隙水压力固定（Pore-pressurefixed）：这是一个可渗透边界，沿着使用：zonegridpointfixpore-pressurefrangepore-pressure固定。或者，用zonefaceapplypore-pressure命令设置外部或者内部的边界条件。需 ，如果不是设定外边界面上的孔隙水压力，必须用interior 。用zonefaceapply命令的一个好处就是：可以用history控制孔隙水压力（使用fish功能）。zonezonegridpointfixwell将规定的流入或流出应用于边界网格点。如果添加modelmodelconfigurefluid;modelmodelconfigurefluid;geometry--zonecreatecylindersize201010point1(020)point2(001)point3(20;fluidflowmodel--zonefluidcmodelassignzonefluidpropertypermeability0.5e-4porosity0.5biot1zonefluidbiotonzonegridpointinitializebiot;iffluidmodulusis;setflbiot;inifmoduluszoneinitializefluid-densityzonegridpointinitializefluid-tension;initialconditions--modelrangecreate'mound'cylinderend-100-1end-2002radius1.01zonegridpointinitializesaturation0.0zonegridpointinitializesaturation1.0rangenamed-rangemoundzonegridpointinitializepore-pressure1e4grad00-1e4...rangenamed-range;settingsmodelgravity00-10modelmechanicalactiveoffmodelfluidactiveonmodelsave;test--modelsolvefluidtime-totalmodeltitle"Spreadingofagroundwatermound(t=0.35)"modelsave'mound1'modelsolvefluidtime-totalmodeltitle"Spreadingofagroundwatermound(t=modelsavemodelsolvefluidtime-totalmodeltitle"Spreadingofagroundwatermound(t=0.65)"modelsave'mound3'modelsolvefluidtime-totalmodeltitle"Spreadingofagroundwatermound(t=0.85)"modelsave'mound4'modelmodelmodelconfigurefluidzonecreatebricksize40120point0000point12000point2010point300;mechanicalmodel--zonepropertybulk5e7shear3e7fric25cohesion1e5tens1e10zonegridpointfixvelocity-xrangeunionposition-x0position-x20zonegridpointfixvelocity-y;applyloadslowlyfishdefinerampzonefaceapplystress-normal=-0.3e6fish@ramp...rangeposition-x03position-z10;fluidflowmodel--;porepressurefixedatzeroatthesur