三个流固耦合分析实例

合集下载

ANSYS Workbench轴流叶轮机械流固耦合分析实例

ANSYS 14.0中Workbench提供了进行流固耦合（FSI）分析的模块，可以十分方便的对轴流叶轮机械进行气动载荷分析，包括最大变形量和等效应力分布。

1.进入ANSYS14.0 Workbench界面。

2.在左下角中的custom system模块中选择第一个流固耦合模块FSI：Fluid Flow(CFX)-staticstructural，双击。

3.屏幕中出现了FSI模块。

4.右击A5(solution)选择import solution，导入已经计算完毕的CFX结果.res文件。

5.导入结果后的界面如下图所示。

CFX部分已经完成了计算，所以不需要额外的设置。

6.双击B3（Geometry）进入结构分析的几何单元，初始单位选择meter。

7.导入一个叶片的几何实体，可以选择的几何文件类型很多，x_t、iges等等都可以。

在CFX中，我们通常计算的都是多个转子，多个叶片，但是在分析流固耦合时，只需导入自己关心的那个叶片就可以了。

8.然后点击Generate，就可以看到生成的叶片实体了。

8.关闭Geometry窗口回到Workbench截面，可以看到此时B3(Geometry)后已经变成了绿色的√，说明生成正确。

9.双击B4（model）进入。

可以看到Geometry、coordinate system、connections等项目前面已经是绿色的对号，不需要再进行设置。

10.单击mesh，在左下角的Details of mesh，如图进行设置。

10.右击mesh，选择generate mesh生成网格。

11.生成的叶片网格如图所示。

12.点击static structural ，选择工具栏中的support 下的fixed support，为叶片根部添加约束。

13.选中叶根面，点击左下角中的Apply，完成约束添加。

14.点击上工具栏中units，选择转速单位为RPM.15.如图所示添加转速16.按自己的算例输入转速。

第九讲流固耦合

刚体欧拉耦合
全耦合– 水下爆炸

水下爆炸对舰艇的影响

刚体欧拉耦合
刚体壳
水和空气使用3D多物质欧拉
全耦合 – 水面爆炸对舰艇的影响

靠近舰艇的空气中爆炸

刚体欧拉耦合刚体壳水和空气使用3D多物质欧拉
全耦合 – 爆炸对舰艇壳体的冲击变形

靠近舰艇的空气中爆炸

刚体欧拉耦合刚体壳水和空气使用3D多物质欧拉

覆盖的体和面在每个时间不均参与计算

效率没有刚性耦合高为有效的载荷周期进行耦合计算

不需要欧拉网格时，将其删除
全耦合 - 建筑物中爆炸
砖混结构的炸药库中27 吨炸药爆炸 (ESH)
欧拉冲击波求解器
全耦合 - 内弹道
气体在膛内流动，出口处膨胀泄压欧拉冲击波求解器
全耦合- 水下爆炸

炸药和水用多物质欧拉求解器钢质筒使用拉格朗日求解器
全耦合 - 爆炸侵彻钢筋混凝土
Euler Blast
Lagrange Concrete
Beam Reinforcements
全耦合- 地雷爆炸

空气爆炸

全耦合

破片碰撞

接触

侵蚀

余留的惯性
全耦合 – 爆炸侵彻 RPG

RPG爆炸冲击波和破片对 CFRP翼箱的破坏

空气中爆炸采用冲击波求解器 RPG壳体 (破片) 和翼箱
全耦合 – 玻璃碎片

有效减低飞散玻璃碎片的危险，德国国防部对各种汽车玻璃进行了安全评估 :

Test in Large Blast Simulator

流固耦合

流固耦合定义：它是研究变形固体在流场作用下的各种行为以及固体位形对流场影响这二者相互作用的一门科学。

流固耦合力学的重要特征是两相介质之间的相互作用,变形固体在流体载荷作用下会产生变形或运动。

变形或运动又反过来影响流，从而改变流体载荷的分布和大小，正是这种相互作用将在不同条件下产生形形色色的流固耦合现象。

(一)流固耦合动力学：求解方法与基本理论---张阿漫，戴绍仕●有限元法●边界元法●SPH法与谱单元法●瞬态载荷作用下流固耦合分析方法●小尺度物体的流固耦合振动●水下气泡与边界的耦合效应按耦合机理分两大类：1 耦合作用只发生在两相交界面---界面耦合（场间不相互重叠与渗透），耦合作用通过界面力（包括多相流的相间作用力等）起作用。

它的计算只要满足耦合界面力平衡，界面相容就可以了（其耦合效应是通过在方程中引入两相耦合面边界条件的平衡及协调关系来实现的）。

如气动弹性，水动弹性等。

按照两相间相对运动的大小及相互作用分为三类：(1)流体和固体结构之间有大的相对运动问题"最典型的例子是飞机机翼颤振和悬索桥振荡中存在的气固相互作用问题,一般习惯称为气动弹性力学问题"(2)具有流体有限位移的短期问题"这类问题由引起位形变化的流体中的爆炸或冲击引起"其特点是:我们极其关心的相互作用是在瞬间完成的,总位移是有限的,但流体的压缩性是十分重要的"(3)具有流体有限位移的长期问题"如近海结构对波或地震的响应!噪声振动的响应!充液容器的液固耦合振动!船水响应等都是这类问题的典型例子"对这类问题,主要关心的是耦合系统对外加动力荷载的动态响应"2 两域部分或全部重叠在一起，难以明显的分开，使描述物理现象的方程，特别是本构方程需要针对具体的物理现象来建立，其耦合效应应通过建立与不同单相介质的本构方程等微分方程来体现。

按耦合求解方法分两大类：1 直接耦合求解：直接耦合是在一个求解器中同时求解不同物理场的所有变量，需要针对具体的物理现象来建立本构方程，其耦合效应通过描述问题的微分方程来体现。

滚筒洗衣机ABAQUS流固耦合实例分析(步骤)

例子的来源是Abaqus CLE的官方教程，可是写的太粗线条，我还是搞了两天才做出了这个例子。

其实就是个滚筒洗衣机带着洗衣机里的水一起转的问题。

1. 分别为Eulerian domain和Lagrangian domain建立两个part建立Lagrangian domain的Part，类型设置为Discrete rigid，并设置Reference Point。

建立Eulerian domain的Part，类型设置为Eulerian，要注意Eulerian domain 和Lagrangian domain要保证有重叠的部分，这是一种弱耦合，数据在两个区域间抛来抛去，所以网格要有重叠部分。

这导致在Eulerian domain里有的部分是有材料的，有的地方是没有材料的。

为了之后设置材料分布时候方便，要把part实现划出几个辅助的partition。

黄色虚线是在划分partition时，为了指明Extrude/Sweep方向用到的辅助坐标轴。

2. 定义水的材料属性选择状态方程模型EOS中Us-Up，设置声速c0=1483m/s；密度为1000kg/m3；粘度为0.001kg/ms。

并把截面属性赋给Eulerian domain。

3. 把两个Part组装起来4. 新建一个Step-15. 为Eulerian domain和Lagrangian domain划分网格6. 设置接触新建一个Contact Property ，因为不是普通的面和面的接触，水中的任何的一个部分可能在流动区域里的任何一个地方和Lagrangian domain接触，设置Tangential Behavior为Rough，赋给水和洗衣机之间的关系。

新建一个Interaction，把刚才的Contact Property赋给它。

更重要的是设置接触的两个Surface。

其中一个Surface是Lagrangian domain 部分的内侧面，为Geometry类型，另一个Surface是Eulerian domain的全部网格，为Mesh类型。

流固耦合模态分析实例1

圆环在水中（模态分析）【ANSYS】默认分类2009-11-18 16:48:28 阅读31 评论0 字号：大中小finish/clear/PREP7!定义单元类型ET,1,PLANE42 ! structural elementET,2,FLUID29 ! acoustic fluid element with ux & uyET,3,129 ! acoustic infinite line elementr,3,0.31242,0,0ET,4,FLUID29,,1,0 ! acoustic fluid element without ux & uy!材料属性MP,EX,1,2.068e11MP,DENS,1,7929MP,NUXY,1,0MP,DENS,2,1030MP,SONC,2,1460! 创建四分之一模型CYL4,0,0,0.254,0,0.26035,90CYL4,0,0,0.26035,0,0.31242,90! 选择属性,网格划分ASEL,S,AREA,,1AATT,1,1,1,0LESIZE,1,,,16,1LESIZE,3,,,16,1LESIZE,2,,,1,1LESIZE,4,,,1,1MSHKEY,1MSHAPE,0,2D ! mapped quad meshAMESH,1ASEL,S,AREA,,2AATT,2,1,2,0LESIZE,5,,,16,1LESIZE,7,,,16,1LESIZE,6,,,5LESIZE,8,,,5MSHKEY,0MSHAPE,0,2D ! mapped quad meshAMESH,2! 关于Y轴镜像nsym,x,1000,all ! offset node number by 1000esym,,1000,all! 关于y轴镜像nsym,y,2000,all ! offset node number by 2000esym,,2000,allNUMMRG,ALL ! merge all quantitiesesel,s,type,,1nsle,sesln,s,0 选择附在节点上的单元（0表示单元有一个节点被选中即选中单元）nsle,sesel,inve 反选单元nsle,semodif,all,type,4esel,allnsel,all! 指定无限吸收边界csys,1nsel,s,loc,x,0.31242type,3real,3mat,2esurfesel,allnsel,all! 标识流固交接面nsel,s,loc,x,0.26035esel,s,type,,2sf,all,fsi,1nsel,allesel,allFINISH/soluantype,modalmodopt,damp,10mxpand,10,,,yessolvefinish水箱盛水【ANSYS】默认分类2009-11-18 16:48:38 阅读61 评论0 字号：大中小length=2width=3height=2/prep7et,1,63et,2,30 !选用FLUID30单元，用于流固耦合问题r,1,0.01mp,ex,1,2e11mp,nuxy,1,0.3mp,dens,1,7800mp,dens,2,1000 !定义Acoustics材料来描述流体材料-水mp,sonc,2,1400mp,mu,0,!block,,length,,width,,heightesize,0.5mshkey,1!type,1mat,1real,1asel,u,loc,y,widthamesh,allalls!type,2mat,2vmesh,allfini/soluantype,2modopt,unsym,10 !非对称模态提取方法处理流固耦合问题eqslv,frontmxpand,10,,,1nsel,s,loc,x,nsel,a,loc,x,lengthnsel,r,loc,yd,all,,,,,,ux,uy,uz,nsel,s,loc,y,width,d,all,pres,0allsasel,u,loc,y,width,sfa,all,,fsi !定义流固耦合界面allssolvfini/post1set,firstplnsol,u,sum,2,1finiANSYS 声学计算算例水下圆柱壳体的建模与声学分析使用有限元软件ANSYS进行计算和分析时水下环肋圆柱壳体有限元模型的建立及结构声学分析主要分为以下一些步骤：1.建立壳体的实体模型（包括有圆柱壳体的建立，给圆柱壳体加环肋）；2.圆柱壳体外部流体介质的生成；3.对圆柱壳体和流体介质进行有限元4.设置流固耦合单元，并设置外部声场边界条件；5.在求解器中进行振动模态求解和受激励的谐响应求6.求解结果进行后处理分析。

fluent流固耦合传热算例

fluent流固耦合传热算例【原创实用版】目录1.Fluent 流固耦合传热简介2.Fluent 软件的应用范围3.流固耦合传热的算例分析4.Fluent 软件在流固耦合传热中的应用技巧5.总结正文一、Fluent 流固耦合传热简介流固耦合传热是一种复杂的热传递过程，涉及到流体和固体之间的相互作用。

在这种过程中，流体与固体之间的热传递机制和热流动特性都需要考虑。

Fluent 是一款强大的计算流体力学（CFD）软件，可以模拟流固耦合传热过程，为研究人员和工程师提供可靠的解决方案。

二、Fluent 软件的应用范围Fluent 软件广泛应用于各种流体动力学问题的仿真和分析中，包括流固耦合传热问题。

它可以模拟多种流体流动和传热模式，如强制对流、自然对流和湍流等。

同时，Fluent 也可以考虑固体的热传导和热膨胀等特性，为研究者提供全面的热传递分析手段。

三、流固耦合传热的算例分析在 Fluent 中，可以通过设置耦合界面和热流边界条件来模拟流固耦合传热问题。

例如，可以考虑一个流体与固体相接触的系统，通过调整流体和固体的热传导系数、对流换热系数等参数，观察不同条件下的热传递特性。

四、Fluent 软件在流固耦合传热中的应用技巧为了获得准确的仿真结果，需要注意以下几点：1.网格划分：在仿真中，需要对流体和固体部分进行适当的网格划分，以确保计算精度。

2.耦合设置：在设置耦合界面时，需要选择正确的耦合方式，如耦合热流或耦合应力等。

3.边界条件：在设置热流边界条件时，需要考虑流体与固体之间的热交换方式，如对流换热或传导换热等。

4.物质属性：需要正确设置流体和固体的物质属性，如比热容、密度和热传导系数等。

五、总结Fluent 软件在流固耦合传热方面的应用具有广泛的实用性，可以模拟各种复杂的热传递过程。

流固耦合的定义

流固耦合的定义《流固耦合那些事儿》嘿，大家好！今天咱来唠唠“流固耦合”这个听起来有点高大上的玩意儿。

流固耦合嘛，简单说就是流体和固体之间相互影响、相互作用的一种关系。

你想想看哈，水在河里哗哗流，那河边的石头是不是老被水冲来冲去？这就是流体（水）对固体（石头）的作用。

反过来呢，石头要是长得奇形怪状的，也会影响水流的走向吧？这就是固体对流体的影响。

这就是流固耦合在我们生活中常见的例子。

流固耦合就像是一场“拔河比赛”。

流体和固体就像是两队选手，互相拉扯、较劲。

有时候流体力量大，就会把固体冲得摇摇晃晃；有时候固体很顽固，就会让流体拐个弯儿走。

就拿风吹大树来说吧。

风就是那流体，大树就是那固体。

大风呼呼一吹，大树就得跟着晃悠，要是风再大点儿，说不定大树就被吹倒啦！这就是风（流体）和大树（固体）之间的流固耦合。

再比如，咱们盖房子的时候。

房子是固体吧，可要是遇到地震了，那地面晃来晃去的，就跟流体似的。

这时候房子就得承受这种晃动感，要是房子不够结实，那就得遭殃咯！流固耦合的世界啊，充满了各种有趣又奇妙的现象。

有时候让我们头疼，有时候又让我们惊叹。

搞科研的人呢，就喜欢研究这些，想搞清楚流体和固体之间到底是咋互动的。

不过研究流固耦合可不是件容易的事儿啊！就像看一场复杂的“武林争霸赛”，得时刻盯着双方的招式，还得分析出其中的门道。

有时候一个不小心，就可能被复杂的现象给绕晕咯！但这也是它的魅力所在，让人越研究越觉得有意思。

总之呢，流固耦合就是这么个神奇的东西。

它在我们生活中无处不在，影响着我们的方方面面。

无论是大自然中的现象，还是我们人类制造的各种工程，都少不了流固耦合的身影。

所以啊，咱们可得好好认识它、了解它，这样才能更好地和这个奇妙的世界相处嘛！你说是不是这个理儿？。

多孔介质流固耦合模型

多孔介质流固耦合模型全文共四篇示例，供读者参考第一篇示例：多孔介质是一种具有孔隙结构的介质，其内部空间可以充满气体或液体。

在工程应用中，多孔介质广泛存在于地下水层、土壤、岩石等领域，其流体运动受到物理和化学环境的复杂影响，因此需要建立多孔介质流固耦合模型来研究其流动特性。

多孔介质流固耦合模型是将多孔介质内部的固相颗粒和流体相耦合起来，通过数学方程描述多孔介质内部的流体流动和固体颗粒的运动规律。

多孔介质流固耦合模型的构建涉及多孔介质内部流动的描述、固相颗粒的运动规律、力学和流体动力学方程的建立等多个方面。

在多孔介质内部的流动描述方面，常用的方法包括达西定律和布里渊方程。

达西定律是描述多孔介质内部流体流动速度与渗透率之间的关系的经典定律，其表示为：\[ v = -K \frac{∇p}{μ} \]v为流速，K为渗透率，p为压力，μ为粘度。

布里渊方程则是描述多孔介质内部流体流动速度随位置变化的变化规律，其表示为：ρ为流体密度，g为重力加速度，z为高度。

达西定律和布里渊方程可以结合使用，建立多孔介质内部流动的数学模型。

固相颗粒的运动规律是多孔介质流固耦合模型的另一个关键方面。

通常情况下，多孔介质内部的固相颗粒会受到流体作用力和固相颗粒之间的相互作用力的影响，其运动规律可以用牛顿第二定律和达西定律描述。

牛顿第二定律表示为：\[ F = ma \]F为作用力，m为质量，a为加速度。

达西定律描述了固相颗粒受到的流体作用力与流速之间的关系。

通过以上分析，我们可以看出，多孔介质流固耦合模型是描述多孔介质内部流动和固相颗粒运动的重要数学模型，其建立需要考虑多个方面的因素。

在实际工程应用中，多孔介质流固耦合模型可以用于研究地下水流动、土壤固结、岩石力学等问题，为工程设计和科学研究提供重要的参考依据。

希望通过本文的介绍，读者能够对多孔介质流固耦合模型有更深入的了解，并在实际应用中取得更好的研究成果。

第二篇示例：多孔介质是一种具有孔隙结构的材料，它可以允许流体通过或在其中存储流体。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

length=2 !定义体各种变量参数，长宽高 width=3 height=2 /prep7 et,1,63 !选用壳模型 et,2,30 !选用FLUID30单元，用于流固耦合问题 r,1,0.01 增加实常数，壳厚为0.01 mp,ex,1,2e11 mp,nuxy,1,0.3 mp,dens,1,7800 !定义壳单元的各种单元属性 mp,dens,2,1000 !定义Acoustics材料来描述流体材料-水 mp,sonc,2,1400 !定义声单元声速 mp,mu,0, !定义吸声系数 ! block,,length,,width,,height !建立长方体 esize,0.5 mshkey,1 ! type,1 !选择壳单元 mat,1 real,1 asel,u,loc,y,width !选择面 amesh,all !划分面单元 alls !选择所有项 ! type,2 !选择声单元 mat,2 vmesh,all !划分体单元 fini

/solu antype,2 modopt,unsym,10 !非对称模态提取方法处理流固耦合问题 eqslv,front mxpand,10,,,1 nsel,s,loc,x, nsel,a,loc,x,length nsel,r,loc,y d,all,,,,,,ux,uy,uz, nsel,s,loc,y,width, d,all,pres,0 !上面几步为定义边界条件和约束 alls asel,u,loc,y,width, sfa,all,,fsi !定义流固耦合界面 alls !选择所有项 solv !求解 fini

/post1 !后处理 set,first plnsol,u,sum,2,1 !显示图形 fini

/PREP7 !定义壳材料与性质 !壳元素与材料 ET,1,shell63 $MP,EX,1,201E9 $MP,prxy,1,0.26 $MP,dens,1,7.85E3 $r,1,0.006 !流体元素与材料 ET,2,FLUID80 $MP,EX,2,1.5e9 $MP,DENS,2,0.84e3 $mp,visc,2,1.0e-10

!以下这个keyoption怎么用? 如过用1，就会显示[Element 877 may not have a positive Z coordinate IF KEYOPT(2) = 1.]，显示这个错误代表要做什么修正吗?所以我暂时用KEYOPT(2) = 0就可以跑。 KEYOPT,2,2,0

!建立壳关键点 K,1,10,0,0 $K,2,10,0,12

!建立中心线关键点 k,3,0,0,0 $k,4,0,0,20

!定义壳壁线 L,1,2 $L,1,3

!以关键点3,4为中心线旋转360度生成壳体 AROTAT,all,,,,,,3,4,360

!划分壳体网格 AATT,1,1,1 $esize,2 $mshape,0,3D $mshkey,2 $amesh,all $alls

!延伸出水位体积 VEXT,2,8,2,0,0,10,0,0,0 $vglue,all csys,1 !划分水位网格 type,2 $mat,2 $esize,2 $mshape,0,3D $mshkey,1 $vmesh,all alls !以上建模应该没太大问题

!以下是耦合，我在流固界面上的网格是重合节点，特别是下面这两段落我很不确定该怎么设定，感觉问题就出在这边了!这里解决了应该就可以。要怎么改?或是用CP? 或是NUMMRG? 重点是流体和固体要一起动，通常设定不好就流体自己动，或是流体都跑到壳体外面去了，流体跟壳不应该穿越，而是一起有行为。 csys,1 !将工作平面定义为柱坐标。 nsel,s,loc,x,10 nrotate,all !旋转节点坐标系。 CPINTF,UX,0.0001, !将径向约束（即X方向）加到节点上。

nsel,s,loc,z,0 nrotate,all CPINTF,UZ,0.0001,

!边界条件，将底部固定，并给予Z方向加速度。 NSEL,S,LOC,Z,0 $D,ALL,ALL $acel,,,9.8

fini /solu antype,modal modopt,reduc,10,, mxpand,10,

csys,1 !Z上柱坐标系 !定义主自由度，由图显示感觉是没问题，但我也不太确定。 Esel,s,type,,1 !选择壳 Nsle,s,all !所有点 Nsel,u,loc,z,0 !排除边界条件 m,all,ux !(径向)x方向的主自由度

Esel,s,type,,2 !选择液体 Nsel,s,loc,z,10 !再选择液面表面 m,all,uz !(竖向)z方向的主自由度

alls solve fini 水坝空库 /BATCH KEYW,PR_SET,1 KEYW,PR_STRUC,1 KEYW,PR_THERM,0 KEYW,PR_FLUID,0 KEYW,PR_ELMAG,0 KEYW,MAGNOD,0 KEYW,MAGEDG,0 KEYW,MAGHFE,0 KEYW,MAGELC,0 KEYW,PR_MULTI,0 KEYW,PR_CFD,0 /GO !*

/prep7 !* define material proterties mp,dens,1,2650. !mat 1 for dam mp,ex,1,3.15e10 mp,prxy,1,.167

!* define element type et,1,PLANE42,,,2

!* define geometry lwater=618. k, 1,0.0,0.0,0.0 k, 2,70.2,0.,0. k, 3,0.,66.5,0. k, 4,21.9875,66.5,0. k, 5,0.,103.,0. k, 6,14.8,103.,0.

a, 1, 2, 4, 3 a, 3, 4, 6, 5

asel,s,loc,x,0.,100. aatt,1,,1 cm,adam,area allsel,all !* mesh geometry ESIZE,0,10 lsel,s,loc,y,0.1,66. !坝下部剖分分数 lesize,all,,,15

mshape,0,2D mshkey,1 allsel,all amesh,all

finish /solu

antype,modal MODOPT,LANB,30 MXPAND,30, , ,0

esel,s,mat,,1 !坝体约束 nsle,s nsel,r,loc,y,-1.0,1.0 d,all,ux,0. d,all,uy,0.

/pbc,all,,1 /pnum,type,1 /number,1 gplot

allsel,all save solve 满库 /BATCH KEYW,PR_SET,1 KEYW,PR_STRUC,1 KEYW,PR_THERM,0 KEYW,PR_FLUID,0 KEYW,PR_ELMAG,0 KEYW,MAGNOD,0 KEYW,MAGEDG,0 KEYW,MAGHFE,0 KEYW,MAGELC,0 KEYW,PR_MULTI,0 KEYW,PR_CFD,0 /GO !*

/prep7 !* define material proterties mp,dens,1,2650. !mat 1 for dam mp,ex,1,3.15e10 mp,prxy,1,.167 mp,dens,2,1000. !mat 2 for water mp,sonc,2,1440

!* define element type et,1,PLANE42,,,2 et,2,29 et,3,29,,1

!* define geometry lwater=618. k, 1,0.0,0.0,0.0 k, 2,70.2,0.,0. k, 3,0.,66.5,0. k, 4,21.9875,66.5,0. k, 5,0.,103.,0. k, 6,14.8,103.,0. k, 7,-1.*lwater,103.,0. k, 8,-1.*lwater,66.5,0. k, 9,-1.*lwater,0.,0.

a, 1, 2, 4, 3 a, 3, 4, 6, 5 a, 8, 3, 5, 7 a, 9,1, 3, 8

asel,s,loc,x,0.,100. aatt,1,,1 cm,adam,area asel,s,loc,x,-1*lwater,0. aatt,2,,3 cm,awater,area

allsel,all !* mesh geometry ESIZE,0,10 lsel,s,loc,y,0.1,66. !坝下部剖分分数 lesize,all,,,15 lsel,s,loc,x,-0.1,-1*lwater-1. !水体长度方向剖分分数 lsel,r,loc,y,-0.1,67. lesize,all,,,40,0.5 lsel,s,loc,x,-0.1,-1*lwater-1. !水体长度方向剖分分数 lsel,r,loc,y,100.,104. lesize,all,,,40,2.0