flac3d流固耦合命令流

第1部分 命令流按照顺序进行

2-1定义一个FISH函数 new def abc

abc = 25 * 3 + 5 End print abc

2-2使用一个变量 new def abc hh = 25

abc = hh * 3 + 5 End Print hh Print abc

2-3对变量和函数的理解 new def abc hh = 25

abc = hh * 3 + 5 End

set abc=0 hh=0 print hh print abc print hh new def abc

abc = hh * 3 + 5 end

set hh=25 print abc

set abc=0 hh=0 print hh print abc print hh

2-4获取变量的历史记录 new

gen zone brick size 1 2 1 model mohr

prop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10

fix x y z range y -0.1 0.1

apply yvel -1e-5 range y 1.9 2.1 plot set rotation 0 0 45 plot blo

new

gen zone radtun p0 0,0,0 p1 25,0,0 p2 0,50, 0 p3 0,0,25 size 4 25 4 10 dim 4 4 4 4 ratio 1 1 1 1.1 fill

gen zone reflect normal 1 0 0 ori 0 0 0 gen zone reflect normal 0 0 1 ori 0 0 0 mo mohr

pro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5 ini dens 2000

fix x range x -25.1,-24.9 fix x range x 24.9 25.1 fix y range y 49.9 50.1 fix z range z -25.1 -24.9 fix z range z 24.9 25.1 PLOT SURF PAUSE K

sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10 sel cable id=1 beg 0,29,0 end 0,3

1、怎样查看模型？

答：plot grid 可以查看网格，plot grid num 可以查看节点号。

2、请问在圆柱体四周如何施加约束条件？

答：可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not，其中r2

fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 1 0 0 end2 1 4 0 rad 1

3、怎么能把一个PLOT的图像数据导出来以便用其他软件绘图？ 答：用set log on 命令，把数据导出来，转到excel里处理一下，然后用surfer或者什么作图软件绘制就行了。

4、用命令建立模型后，如何显示点的坐标？ 答：使用plo blo gro gpnum on命令

5、关于gauss_dev对性质进行高斯正态分布的问题？

答：根据手册上的说明：下面的命令设定一个平均摩擦角为40度，标准方差是±5％。则命令如下：prop friction 40 gauss_dev 2

问题：请问gauss_dev 2中的2是如何计算的？如果把±5％改为±10％，则命令应如何写？40×

ANSYS流固耦合分析实例命令流

达尔文档 DareDoc 分享知识 传播快乐 ANSYS流固耦合分析实例命令流

本资料来源于网络，仅供学习交流 2015年10月 达尔文档|DareDoc 整理

1

DareDoc

ANSYS流固耦合分析实例命令流

目 录

ANSYS流固耦合例子命令流 .......................................................................... 错误！未定义书签。 ANSYS流固耦合的方式 ..................................................................................................................... 3 一个流固耦合模态分析的例子1 ....................................................................................................... 3 一个流固耦合模态分析的例子2 ...............................

hist reset ;清空里面所有的监测点

nstress=-40e3 hist ramp的意思，很简单 nstress = -40e3 * ramp hist相当于乘法，为了方便控制动荷载的幅值。

plotitems——Add——strucktural elements——geometry ——cables 显示锚杆的几何形状

Pl-add-vector-apply body force 显示施加的力，不过是在求解之后显示 Creat 创建一新视图

plotitems------add-----vetor----apply body force 显示荷载的施加位置

Array var(4,3) 二维数组 array var(5)一维5个元素的数组 Var(m,n)=urand ;随机 var—数组的名字 Pnt=z_next(pnt) 查找下一个单元体表 Pl plas 显示塑性区

Array x_b(x,y) z_b(x,y) 创建两个x行y列的表格，可以给其赋值，以便调用。Array就是数组的意思。

Hist id=1 gp xd 7 0 0 监测7 0 0 点处的x方向位移，ID=1 材料参数dilation 剪胀角

hist reset ;清空里面所有的监测点

nstress=-40e3 hist ramp的意思，很简单 nstress = -40e3 * ramp hist相当于乘法，为了方便控制动荷载的幅值。

plotitems——Add——strucktural elements——geometry ——cables 显示锚杆的几何形状

Pl-add-vector-apply body force 显示施加的力，不过是在求解之后显示 Creat 创建一新视图

plotitems------add-----vetor----apply body force 显示荷载的施加位置

Array var(4,3) 二维数组 array var(5)一维5个元素的数组 Var(m,n)=urand ;随机 var—数组的名字 Pnt=z_next(pnt) 查找下一个单元体表 Pl plas 显示塑性区

Array x_b(x,y) z_b(x,y) 创建两个x行y列的表格，可以给其赋值，以便调用。Array就是数组的意思。

Hist id=1 gp xd 7 0 0 监测7 0 0 点处的x方向位移，ID=1 材料参数dilation 剪胀角

hist reset ;清空里面所有的监测点

nstress=-40e3 hist ramp的意思，很简单 nstress = -40e3 * ramp hist相当于乘法，为了方便控制动荷载的幅值。

plotitems——Add——strucktural elements——geometry ——cables 显示锚杆的几何形状

Pl-add-vector-apply body force 显示施加的力，不过是在求解之后显示 Creat 创建一新视图

plotitems------add-----vetor----apply body force 显示荷载的施加位置

Array var(4,3) 二维数组 array var(5)一维5个元素的数组 Var(m,n)=urand ;随机 var—数组的名字 Pnt=z_next(pnt) 查找下一个单元体表 Pl plas 显示塑性区

Array x_b(x,y) z_b(x,y) 创建两个x行y列的表格，可以给其赋值，以便调用。Array就是数组的意思。

Hist id=1 gp xd 7 0 0 监测7 0 0 点处的x方向位移，ID=1 材料参数dilation 剪胀角

作为流流合版块的版主，我感到惭愧。因为我几乎就没真正应用流固耦合做过工程。第一次应用流固耦合还要追溯到做硕士毕业论文的时候，当时做的是高压水射流切割，属于一个大课题中的小项。主要用的软件是fluent。但是利用fluent是没办法计算射流的切割效果的，流体软件只能计算流场参数（压力、速度、温度等），对于应力计算实在是力不从心。我不知道导师是从哪里听来的风声，说让使用mpcci将fluent与abaqus耦合计算固体变形乃至断裂。当时也是初生牛犊不怕虎，老师说用那就用呗，于是开始关注固体计算，关注abaqus，关注mpcci。然而现实是残酷的，流体与固体采用不同的计算网格（流体用欧拉网格，固体采用拉格朗日网格），对于断裂的问题，单纯采用abaqus勉强可算，然而耦合上流体之后，通常计算会以出现负体积而告终。多次的失败终于磨灭了导师的耐心，于是项目转而采用LS-DYNA的ALE进行解决，而我的毕业论文，则彻底的舍弃了这一部分。搞射流的自然离不开喷嘴的设计，在研究射流喷嘴结构在高压流体作用下的材料行为，于是又涉及到了流固耦合问题，这次很幸运，虽然压力很高，然而压差并不大，喷嘴的变形处于弹性小变形阶段，我采用workbench中的CFX+ANSYS

作为流流合版块的版主，我感到惭愧。因为我几乎就没真正应用流固耦合做过工程。第一次应用流固耦合还要追溯到做硕士毕业论文的时候，当时做的是高压水射流切割，属于一个大课题中的小项。主要用的软件是fluent。但是利用fluent是没办法计算射流的切割效果的，流体软件只能计算流场参数（压力、速度、温度等），对于应力计算实在是力不从心。我不知道导师是从哪里听来的风声，说让使用mpcci将fluent与abaqus耦合计算固体变形乃至断裂。当时也是初生牛犊不怕虎，老师说用那就用呗，于是开始关注固体计算，关注abaqus，关注mpcci。然而现实是残酷的，流体与固体采用不同的计算网格（流体用欧拉网格，固体采用拉格朗日网格），对于断裂的问题，单纯采用abaqus勉强可算，然而耦合上流体之后，通常计算会以出现负体积而告终。多次的失败终于磨灭了导师的耐心，于是项目转而采用LS-DYNA的ALE进行解决，而我的毕业论文，则彻底的舍弃了这一部分。搞射流的自然离不开喷嘴的设计，在研究射流喷嘴结构在高压流体作用下的材料行为，于是又涉及到了流固耦合问题，这次很幸运，虽然压力很高，然而压差并不大，喷嘴的变形处于弹性小变形阶段，我采用workbench中的CFX+ANSYS

作为流流合版块的版主，我感到惭愧。因为我几乎就没真正应用流固耦合做过工程。第一次应用流固耦合还要追溯到做硕士毕业论文的时候，当时做的是高压水射流切割，属于一个大课题中的小项。主要用的软件是fluent。但是利用fluent是没办法计算射流的切割效果的，流体软件只能计算流场参数（压力、速度、温度等），对于应力计算实在是力不从心。我不知道导师是从哪里听来的风声，说让使用mpcci将fluent与abaqus耦合计算固体变形乃至断裂。当时也是初生牛犊不怕虎，老师说用那就用呗，于是开始关注固体计算，关注abaqus，关注mpcci。然而现实是残酷的，流体与固体采用不同的计算网格（流体用欧拉网格，固体采用拉格朗日网格），对于断裂的问题，单纯采用abaqus勉强可算，然而耦合上流体之后，通常计算会以出现负体积而告终。多次的失败终于磨灭了导师的耐心，于是项目转而采用LS-DYNA的ALE进行解决，而我的毕业论文，则彻底的舍弃了这一部分。搞射流的自然离不开喷嘴的设计，在研究射流喷嘴结构在高压流体作用下的材料行为，于是又涉及到了流固耦合问题，这次很幸运，虽然压力很高，然而压差并不大，喷嘴的变形处于弹性小变形阶段，我采用workbench中的CFX+ANSYS