当前位置:文档之家› 【达尔整理】ANSYS流固耦合分析实例命令流

【达尔整理】ANSYS流固耦合分析实例命令流

【达尔整理】ANSYS流固耦合分析实例命令流
【达尔整理】ANSYS流固耦合分析实例命令流

达尔文档DareDoc

分享知识传播快乐

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

本资料来源于网络,仅供学习交流

2015年10月达尔文档|DareDoc整理

目录

ANSYS流固耦合例子命令流.......................................................................... 错误!未定义书签。ANSYS流固耦合的方式 (3)

一个流固耦合模态分析的例子1 (3)

一个流固耦合模态分析的例子2 (4)

一个流固耦合建模的例子 (7)

一加筋板在水中的模态分析 (8)

一圆环在水中的模态分析 (10)

接触分析实例---包含初始间隙 (14)

耦合小程序 (19)

流固耦合练习 (21)

一个流固耦合的例子 (22)

使用物理环境法进行流固耦合的实例及讲解 (23)

针对液面晃动问题,ANSYS/LS-DYNA提供三种方法 (30)

1、流固耦合 (30)

2、SPH算法 (34)

3、ALE(接触算法) (38)

脱硫塔于浆液耦合的分析 (42)

ANSYS坝-库水流固耦合自振特性的例子 (47)

空库时的INP文件 (47)

满库时的INP文件 (49)

计算结果 (52)

ANSYS流固耦合的方式

一般说来,ANSYS的流固耦合主要有4种方式:

1,sequential

这需要用户进行APDL编程进行流固耦合

sequentia指的是顺序耦合

以采用MpCCI为例,你可以利用ANSYS和一个第三方CFD产品执行流固耦合分析。在这个方法中,基于网格的平行代码耦合界面(MpCCI) 将ANSYS和CFD程序耦合起来。即使网格上存在差别,MpCCI也能够实现流固界面的数据转换。ANSYS CD中包含有MpCCI库和一个相关实例。关于该方法的详细信息,参见ANSYS Coupled-Field Analysis Guide中的Sequential Couplin

2,FSI solver

流固耦合的设置过程非常简单,推荐你使用这种方式

3,multi-field solver

这是FSI solver的扩展,你可以使用它实现流体,结构,热,电磁等的耦合

4,直接采用特殊的单元进行直接耦合,耦合计算直接发生在单元刚度矩阵

一个流固耦合模态分析的例子1

这是一个流固耦合模态分析的典型事例,采用ANSYS/MECHANICAL可以完成。处理过程中需要注意以下几个方面的问题:

1、单元的选择;

2、流体材料模式;

3、流固耦合关系的定义;

4、模态提取方法。

length=2

width=3

height=2

/prep7

et,1,63

et,2,30 !选用FLUID30单元,用于流固耦合问题

r,1,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

一个流固耦合模态分析的例子2

一实例,水箱采用SHELL63单元,水箱中的水采用FLUID30单元,以下即为整个流固耦合模态计算的命令流文件:

length=1

width=0.6

height=0.8

/prep7

et,1,63

et,2,30 !选用FLUID30单元,用于流固耦合问题

r,1,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,2,

!

block,,length,,width,,height

esize,0.1

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

我得出的结果是:

SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE

1 29.195 1 1 1

2 0.0000 1 1 1

3 33.202 1 2 2

4 0.0000 1 2 2

5 37.598 1 3 3

6 0.0000 1 3 3

7 44.592 1 4 4

8 0.0000 1 4 4

9 47.917 1 5 5

10 0.0000 1 5 5

11 88.015 1 6 6

12 0.0000 1 6 6

13 88.565 1 7 7

14 0.0000 1 7 7

15 97.133 1 8 8

16 0.0000 1 8 8

17 109.49 1 9 9

18 0.0000 1 9 9

19 109.81 1 10 10

20 0.0000 1 10 10

为了比较我还计算了水箱不装水的情况,(水箱体模型的建立和约束与上面的相同,只是把流体部分去掉),计算结果是:

SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE

1 45.108 1 1 1

2 52.591 1 2 2

3 68.922 1 3 3

4 76.373 1 4 4

5 91.398 1 5 5

6 126.49 1 6 6

7 140.43 1 7 7

8 152.15 1 8 8

9 156.77 1 9 9

10 160.79 1 10 10

我想请问:

1。水箱装水时,怎么计算结果中有0频率啊?其中不等于0的频率是水箱体的频率么?2。水箱不装水时的结果和装水时的相差怎么这么大啊?还是我理解错了么?

企盼指教

一个流固耦合建模的例子

/prep7 !进行预处理模块

et,1,30, !定义1号单元为Fluid30 流固耦合单元

et,2,29 !定义2号单元为Fluid29平面流体单元

et,3,30 , ,1 !定义3号单元为Fluid30流体介质单元

et,4,63 !定义4号单元为Shell63壳体单元

et,5,188 !定义5号单元为Beam188梁单元

r,4,0.002 !定义4号单元的厚度为2㎝

mp,dens,4,7800 !定义4号物理属性包括有密度

mp,ex,4,2.1e11 !杨氏模量、

mp,nuxy,4,0.3 !泊松比

mp,sonc,1,1460 !设置水中声速

mp,dens,1,1000 !设置流体密度

sectype, 1, beam, T, ! 选取T型梁

secoffset,,orig !设置梁的方向

secdata,0.04,0.05,0.002,0.02,0,0,0,0,0,0

!所建立的圆柱壳体的参数:圆柱长为50 ㎝,半径为25 ㎝,壳体的壁厚为2 ㎝,

cyl4,0,0,0.25,,5 !形成圆面

k,9,0,0,0 ! 定义原点

k,10,0,0,0.5

lstr,9,10 !通过原点作直线

adrag,5,6,7,8,,,9 !通过放样形成圆柱

wpoff,0,0,0.1

asel,s,,,2,5

asbw,all,,, !移动工作平面与选取的侧面相切

!重复上面操作,形成四个环肋面

wpoff,0,0,-0.4 !工作平面回到原点位置上

k,31,0.2,0,0.1 !定义环肋的方向点

lsel,s,,,20 !选择要划分为环肋的线段

latt,4,5,5,,31,40,1 !定义线段物理属性

lesize,20,,,6 !划分数目

secnum,1

lmesh,20 !划分线段

type,2 ! 选取第二种单元

lsel,s,,,1,4 !选取线段

lesize,all,,,10 !线段划分数目

lesize,all,,,6

amesh,1 !将通过工作平面的面1 进行划分

esize,0.1,0 !选取沿面放样的网格大小

mat,1 !定义放样的形成单元的物理属性

vdrag,1,,,,,,9 !进行放样形成流体介质

一加筋板在水中的模态分析

再给大家一个实例!

考虑结构在水中的自振频率:例子是一加筋板在水中的模态分析。命令流如下:

FINISH

/CLEAR

/FILENAME,plane

/UNITS,SI

/TITLE,plane

/PREP7

!*********ELEMENT DEFINE********

ET,63,63

ET,4,beam4

et,30,fluid30

!****MATERIAL DEFINE*********

MP,EX,1,2.10E11

MP,DENS,1,7850

MP,NUXY,1,0.3

mp,dens,30,1025

mp,sonc,30,1500

mp,mu,30,0.5

!*******REAL CONSTANT***********

r,30,1e-06

r,50,0.05

r,75,0.375e-02,0.78125e-06,0.000016406

k,1

k,4,1

kfill,1,4,2,,1

kgen,4,1,4,1,,1/3,,10

a,1,2,12,11

*do,i,0,2

*do,j,0,2*10,10

a,1+i+j,2+i+j,12+i+j,11+i+j

*enddo

*enddo

!***************************fluid element****************

k,100,-14.5,-14.5

k,101,-14.5,15.5

k,102,15.5,15.5

k,103,15.5,-14.5

k,140,-14.5,-14.5,30

k,141,-14.5,15.5,30

k,142,15.5,15.5,30

k,143,15.5,-14.5,30

a,100,101,102,103,4,14,24,34,33,32,31,21,11,1 a,1,2,3,4,103,100

a,140,141,142,143

a,100,101,141,140

a,101,102,142,141

a,142,143,103,102

a,140,143,103,100

a,14,24,34,33,32,31,21,11,1,2,3,4

asel,u,,,1,

FLST,2,8,5,ORDE,2

FITEM,2,10

FITEM,2,-17

V A,P51X

nummrg,all

alls

MSHKEY,0

MSHAPE,0

esize,1

lsel,s,loc,y,1/3

lsel,r,loc,x,0,1

lsel,r,loc,z,0

latt,1,75,4

lmesh,all

lsel,s,loc,y,2/3

lsel,r,loc,x,0,1

lsel,r,loc,z,0

latt,1,75,4

lmesh,all

lsel,s,loc,x,1/3

lsel,r,loc,y,0,1

lsel,r,loc,z,0

latt,1,75,4

lmesh,all

lsel,s,loc,x,2/3

lsel,r,loc,y,0,1

lsel,r,loc,z,0

latt,1,75,4

lmesh,all

asel,s,,,1,9

aatt,1,50,63

amesh,all

alls

MSHAPE,1,3d

esize,3

vsel,s,,,1

type,30 $mat,30 $real,30

vmesh,all

alls

FINISH

/solu

alls

!**** 求解***********

!*********************

ANTYPE,MODAL

!MODOPT,lanb,25,0

MODOPT,UNSYMM,25,0

SOLVE

FINISH

总是出现error 说矩阵不对称,不可以用lanb计算。

总结:流体单元不能用对称的解法

应该采用非对称解法。

一圆环在水中的模态分析

finish

/clear

/PREP7

!定义单元类型

ET,1,PLANE42 ! structural element

ET,2,FLUID29 ! acoustic fluid element with ux & uy

ET,3,129 ! acoustic infinite line element

r,3,0.31242,0,0

ET,4,FLUID29,,1,0 ! acoustic fluid element without ux & uy !材料属性

MP,EX,1,2.068e11

MP,DENS,1,7929

MP,NUXY,1,0

MP,DENS,2,1030

MP,SONC,2,1460

! 创建四分之一模型

CYL4,0,0,0.254,0,0.26035,90

CYL4,0,0,0.26035,0,0.31242,90

! 选择属性,网格划分

ASEL,S,AREA,,1

AA TT,1,1,1,0

LESIZE,1,,,16,1

LESIZE,3,,,16,1

LESIZE,2,,,1,1

LESIZE,4,,,1,1

MSHKEY,1

MSHAPE,0,2D ! mapped quad mesh AMESH,1

ASEL,S,AREA,,2

AA TT,2,1,2,0

LESIZE,5,,,16,1

LESIZE,7,,,16,1

LESIZE,6,,,5

LESIZE,8,,,5

MSHKEY,0

MSHAPE,0,2D ! mapped quad mesh AMESH,2

! 关于Y轴镜像

nsym,x,1000,all ! offset node number by 1000 esym,,1000,all

! 关于y轴镜像

nsym,y,2000,all ! offset node number by 2000 esym,,2000,all

NUMMRG,ALL ! merge all quantities esel,s,type,,1

nsle,s

esln,s,0

nsle,s

esel,inve

nsle,s

emodif,all,type,4

esel,all

nsel,all

! 指定无限吸收边界

csys,1

nsel,s,loc,x,0.31242

type,3

real,3

mat,2

esurf

esel,all

nsel,all

! 标识流固交接面

nsel,s,loc,x,0.26035

esel,s,type,,2

sf,all,fsi,1

nsel,all

esel,all

FINISH

/solu

antype,modal

modopt,damp,10

mxpand,10,,,yes

solve

finish

为了便于对比,也对圆环在空气中做了模态分析finish

/clear

/PREP7

!定义单元类型

ET,1,PLANE42 ! structural element

!材料属性

MP,EX,1,2.068e11

MP,DENS,1,7929

MP,NUXY,1,0

! 创建四分之一模型

CYL4,0,0,0.254,0,0.26035,90

! 选择属性,网格划分

ASEL,S,AREA,,1

AA TT,1,1,1,0

LESIZE,1,,,16,1

LESIZE,3,,,16,1

LESIZE,2,,,1,1

LESIZE,4,,,1,1

MSHKEY,1

MSHAPE,0,2D ! mapped quad mesh

AMESH,1

! 关于Y轴镜像

nsym,x,1000,all ! offset node number by 1000

esym,,1000,all

! 关于y轴镜像

nsym,y,2000,all ! offset node number by 2000

esym,,2000,all

NUMMRG,ALL

/solu

antype,modal

modopt,lanb,10

mxpand,10,,,yes

solve

finish

在水中的自振频率为

SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE 1-0.19544E-10 1 1 1

2 0.29640E-0

3 1 1 1

3-0.21663E-10 1 2 2

4-0.29640E-03 1 2 2

5 0.30870E-03 1 3 3

6 0.0000 1 3 3

7-0.30870E-03 1 4 4

8 0.0000 1 4 4

9-0.53726E-03 1 5 5

10 0.57522E-11 1 5 5

11 0.53726E-03 1 6 6

12-0.89057E-11 1 6 6

13 0.98059E-01 1 7 7

14 35.232 1 7 7

15 0.98059E-01 1 8 8

16 -35.232 1 8 8

17 0.98061E-01 1 9 9

18 35.233 1 9 9

19 0.98061E-01 1 10 10

20 -35.233 1 10 10

在空气中的自振频率为

SET TIME/FREQ LOAD STEP SUBSTEP CUMULATIVE

1 0.0000 1 1 1

2 0.0000 1 2 2

3 0.73609E-03 1 3 3

4 60.80

5 1 4 4

5 60.805 1 5 5

6 172.9

7 1 6 6

7 172.97 1 7 7

8 334.40 1 8 8

9 334.40 1 9 9

10 546.59 1 10 10

主要有以下疑问:

1)考虑流固耦合,做模态分析时流体单元是否只能用fluid29(2d)和fluid30(3d),对于fluid129和fluid130在耦合中具体起到什么作用,能不能不设,而用边界约束条件代替?

2)流体范围怎样确定,如本例中(CYL4,0,0,0.26035,0,0.31242,90),外半径为0.31242。如果不是环形的,如一块当水板,该怎样考虑?

3)如果不考虑流体的压缩性,把声速设的很大,MP,SONC,2,1e20,就可以了。

4)从自振频率可以看出,在水中和在空气中,圆环的自振频率差别特别大,且振型也大相径庭,为什么?在水中时,模态提取方法用damp(为什么不能用unsym),特征值的虚部代表角频率,为什么第一阶为正,第二阶为负,而第三阶和第四阶都为0,第六阶、八阶、十阶都为负。应该是从小到大才对?

5)在空气中时,模态提取方法用lanb,为什么第一阶第二阶的频率都为0。

请高手指点迷津,急盼中

对以上问题的解答:

频率为零,一般是发生了刚体位移,估计你是把水抽走,而没有限制圆环。

1。圆环在水中振动必然导致波动(其实就是声波)在水中传播,当声波到达水的另一个界面时就会发生反射(除非水和另一个相邻体的声阻抗是匹配的)。水和金属中的声速相差不大,即可压缩性相差不大。两种可压缩性相差不大的物质的相互作用对两者影响都很大。圆环在水中振动,水对圆环的反作用是由于反射波引起的,流固耦合中采用fluid129和130就是最大程度的减弱反射波。

2。声波从圆环开始传播,随着传播距离的增加,波阵面不断增大,振幅不断减小。同时由于水的衰减,声波也不断减弱。如果水的空间越大,则反射波返回圆环的路径越长,衰减也就越多,影响也就越小。fluid129和130对反射波的衰减(通过很小的反射实现)有限,因此还需要水要有足够的空间。fluid129和130离结构应该大于0.2λ(λ=c/f,c为水中声速)。以上的做法在误差允许的情况下等效于水在无限大水空间中的情况。如果是挡水板,水就是有限空间了,情况也不一样。

3。声速加大情况也不一样,就是不知是不是你所要的情况?

4。空气作为介质,由于其声速比金属小很多,可压缩性大很多,影响可以忽略不计。而水的影响就不同了。这可能就是频率和振型不同的原因吧?我试了你的例子,各种提取方法都可以。

5。空气的影响忽略不计,因此需要对圆环进行约束。你没有约束,那么就会发生静态位移即频率为零。圆环有两个对称轴,因此会发生频率成对出现的情况。也就是说,两个方向上有同样的振型。

接触分析实例---包含初始间隙

fini

/clear,nostart

/prep7

et,1,82

KEYOPT,1,3,3

r,1,0.5

mp,ex,1,1e9

mp,prxy,1,0.3

k,1,0,0

k,2,10,0

k,3,10,5

k,4,6.2,5

k,5,7.5,3.4

k,6,2.5,3.4

k,7,3.8,5

k,8,0,5

a,1,2,3,4,5,6,7,8 LFILLT,6,5,0.18, , LFILLT,5,4,0.18, , FLST,2,3,4

FITEM,2,9

FITEM,2,11

FITEM,2,10

AL,P51X

FLST,2,3,4

FITEM,2,13

FITEM,2,14

FITEM,2,12

AL,P51X

FLST,2,3,5,ORDE,2 FITEM,2,1

FITEM,2,-3

AADD,P51X

rect,0,10,4.8,5

ASBA, 4, 1

gap=0.02

k,24,6.2-gap,5

k,25,7.5-gap,3.4

k,26,2.5+gap,3.4

k,27,3.8+gap,5

a,24,25,26,27

LFILLT,4,3,0.2, ,

LFILLT,3,2,0.2, ,

FLST,2,3,4

FITEM,2,7

FITEM,2,10

FITEM,2,8

AL,P51X

FLST,2,3,4

FITEM,2,13

FITEM,2,14

FITEM,2,11

AL,P51X

FLST,3,2,5,ORDE,2

FITEM,3,3

FITEM,3,-4

ASBA, 1,P51X

rect,3.8+gap,6.2-gap,5,10

rect,3.8+gap,3.8+gap+8,10,12 FLST,2,3,5,ORDE,3

FITEM,2,1

FITEM,2,3

FITEM,2,5

AADD,P51X

rect,3.8+gap+8,3.8+gap+8+2,10,12 FLST,2,2,5,ORDE,2

FITEM,2,1

FITEM,2,4

AGLUE,P51X

CYL4,2.0,1.8,0.6

CYL4,7.0,1.8,0.6

FLST,2,3,5,ORDE,3

FITEM,2,2

FITEM,2,4

FITEM,2,-5

AOVLAP,P51X

esize,0.2

amesh,all

FLST,5,135,2,ORDE,32 FITEM,5,485

FITEM,5,576

FITEM,5,-577 FITEM,5,621 FITEM,5,-625 FITEM,5,707 FITEM,5,-711 FITEM,5,716 FITEM,5,741 FITEM,5,-745 FITEM,5,750 FITEM,5,-751 FITEM,5,766 FITEM,5,797 FITEM,5,-798 FITEM,5,854 FITEM,5,888 FITEM,5,-938 FITEM,5,1101 FITEM,5,1103 FITEM,5,1420 FITEM,5,1628 FITEM,5,1653 FITEM,5,1696 FITEM,5,1699 FITEM,5,-1702 FITEM,5,1726 FITEM,5,-1728 FITEM,5,1852 FITEM,5,-1874 FITEM,5,2044 FITEM,5,-2066 CM,_Y,ELEM ESEL, , , ,P51X CM,_Y1,ELEM CMSEL,S,_Y CMDELE,_Y EREF,_Y1, , ,1,0,1,1 CMDELE,_Y1

ET,2,TARGE169 ET,3,CONTA172 R,3, , ,

R,3,0,0,0.1, 10,0,0 R,4, , ,

R,4,0,0,0.1, 10,-0.02,0

lsel,s,,,9

lsel,a,,,5

lsel,a,,,12

nsll,s,1

type,3

real,3

esurf,all

alls,

lsel,s,,,19

lsel,a,,,20

nsll,s,1

type,3

real,4

esurf,all

alls,

lsel,s,,,7

lsel,a,,,3

lsel,a,,,11

nsll,s,1

type,2

real,3

esurf,all

alls,

lsel,s,,,25

lsel,a,,,26

nsll,s,1

type,2

real,4

esurf,all

alls,

FLST,2,2,5,ORDE,2 FITEM,2,4 FITEM,2,-5

DA,P51X,ALL, FLST,2,1,4,ORDE,1 FITEM,2,6

SFL,P51X,PRES,500,

/solu

antype,0

nlgeom,on

outres,all,all

nsubst,200,200,2

neqit,1000

solve

耦合小程序

最近用到耦合,写了一段小程序,奉献出来,与大家共享。

如果有很多节点,每两个节点位置相同,如果将这些杂乱无章的节点耦合,是件很麻烦的事,可用这段程序,轻松解决。

cpnum=0

cmsel,s,n-zhong !需要耦合的节点

*GET,n_num,NODE,,COUNT, , , , !节点总数

*do,i,1,n_num

cmsel,s,n-zhong

*GET,n_lowest,NODE,,NUM,MIN, , , , !号码最小的节点

*GET,n_x,NODE,n_lowest,LOC,X !该节点坐标

*GET,n_y,NODE,n_lowest,LOC,Y

*GET,n_z,NODE,n_lowest,LOC,Z

NSEL,s,LOC,X,n_x-0.3,n_x+0.3 !寻找与该节点位置相同的节点NSEL,R,LOC,Y,n_y-0.3,n_y+0.3

NSEL,R,LOC,z,n_z-0.3,n_z+0.3

cm,n_cp_cp,node !位置相同的节点形成一个组

cmsel,s,n-zhong

cmsel,u,n_cp_cp

cm,n-zhong,node !取消这些点后剩余的点形成组*GET,n_num_1,NODE,,COUNT, , , , !节点总数*if,n_num_1,lt,2,exit !如果节点数小于二则退出

cmsel,s,n_cp_cp

*GET,n_num,NODE,,COUNT, , , ,

*if,n_num,gt,1,then

CP,cpnum+1,ux,all

CP,cpnum+2,uy,all

CP,cpnum+3,uz,all

cpnum=cpnum+3

*else

*endif

*enddo

该段程序可用

CPINTF,UX,0.001

CPINTF,UY,0.001

CPINTF,UZ,0.001

代替

*DO,I,2,296,3

CP,I,UX,I,I+2

*ENDDO

*DO,I,2,296,3

CP,I,UY,I,I+2

*ENDDO

*DO,I,2,296,3

Ansys常见命令流

Ansys命令流 第一天 目标:熟悉ANSYS基本关键字的含义 k --> Keypoints 关键点 l --> Lines 线 a --> Area 面 v --> Volumes 体 e --> Elements 单元 n --> Nodes 节点 cm --> component 组元 et --> element type 单元类型 mp --> material property 材料属性 r --> real constant 实常数 d --> DOF constraint 约束 f --> Force Load 集中力 sf --> Surface load on nodes 表面载荷 bf --> Body Force on Nodes 体载荷 ic --> Initial Conditions 初始条件 第二天 目标:了解命令流的整体结构,掌握每个模块的标识 !文件说明段 /BATCH /TITILE,test analysis !定义工作标题 /FILENAME,test !定义工作文件名 /PREP7 !进入前处理模块标识 !定义单元,材料属性,实常数段 ET,1,SHELL63 !指定单元类型 ET,2,SOLID45 !指定体单元 MP,EX,1,2E8 !指定弹性模量 MP,PRXY,1,0.3 !输入泊松比 MP,DENS,1,7.8E3 !输入材料密度 R,1,0.001 !指定壳单元实常数-厚度...... !建立模型 K,1,0,0,, !定义关键点 K,2,50,0,,

K,3,50,10,, K,4,10,10,, K,5,10,50,, K,6,0,50,, A,1,2,3,4,5,6, !由关键点生成面 ...... !划分网格 ESIZE,1,0, AMESH,1 ...... FINISH !前处理结束标识 /SOLU !进入求解模块标识 !施加约束和载荷 DL,5,,ALL SFL,3,PRES,1000 SFL,2,PRES,1000 ...... SOLVE !求解标识 FINISH !求解模块结束标识 /POST1 !进入通用后处理器标识 ...... /POST26 !进入时间历程后处理器 …… /EXIT,SAVE !退出并存盘 以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助 /ANGLE !指定绕轴旋转视图 /DIST !说明对视图进行缩放 /DEVICE !设置图例的显示,如:风格,字体等 /REPLOT !重新显示当前图例 /RESET !恢复缺省的图形设置 /VIEW !设置观察方向 /ZOOM !对图形显示窗口的某一区域进行缩放

ansys旋转经典命令流

1 旋转摩擦 (1) 2. 电磁三d命令流实例(论坛看到) (11) 3. 帮助感应加热例子induction heating of a solid cylinder billet (15) 4. 感应加热温度场的数值模拟(论文)inducheat30命令流 (19) 5. 如何施加恒定的角速度?Simwe仿真论坛 (24) 6. 旋转一个已经生成好的物体 (27) 7. 产生这样的磁力线 (28) 8. 旋转摩擦生热简单例子(二维旋转) (32) 8.1. 原版 (32) 8.2. 部分gui操作 (35) 9. VM229 Input Listing (39) 10 轴承---耦合+接触分析 (47) 11. 板的冲压仿真 (52) 1 旋转摩擦 FINISH /FILNAME,Exercise24 !定义隐式热分析文件名 /PREP7 !进入前处理器 ET,1,SOLID5 !选择单元类型 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7800 !定义材料1的密度 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,C,1,,460 !定义材料1的比热 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,KXX,1,,66.6 !定义材料1的热传导系数 MPTEMP,,,,,,,, MPTEMP,1,0 UIMP,1,REFT,,,30 !定义材料1的热膨胀系数的参考温度 MPDATA,ALPX,1,,1.06e-5 !定义材料1的热膨胀系数MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,206e9 !定义材料1的弹性模量 MPDATA,PRXY,1,,0.3 !定义材料1的泊松比 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,2,,8900 !定义材料2的密度 MPTEMP,,,,,,,, MPTEMP,1,0

ansys workbench 流固耦合计算实例

Oscillating Plate with Two-Way Fluid-Structure Interaction Introduction This tutorial includes: ?Features ?Overview of the Problem to Solve ?Setting up the Solid Physics in Simulation (ANSYS Workbench) ?Setting up the Fluid Physics and ANSYS Multi-field Settings in ANSYS CFX-Pre ?Obtaining a Solution using ANSYS CFX-Solver Manager ?Viewing Results in ANSYS CFX-Post If this is the first tutorial you are working with, it is important to review the following topics before beginning: ?Setting the Working Directory ?Changing the Display Colors Unless you plan on running a session file, you should copy the sample files used in this tutorial from the installation folder for your software (/examples/) to your working directory. This prevents you from overwriting source files provided with your installation. If you plan to use a session file, please refer to Playing a Session File. Sample files referenced by this tutorial include: ?OscillatingPlate.pre ?OscillatingPlate.agdb ?OscillatingPlate.gtm ?OscillatingPlate.inp 1.Features This tutorial addresses the following features of ANSYS CFX.

ANSYS命令流实例

/PREP7 !进入前处理 ANTYPE,STATIC !设置分析类型为静力结构分析 PSTRES,ON !用于后面的模态分析中考虑预应力(该开关不影响静力分析) ET,1,LINK10 !选取单元类型1(单向杆单元) KEYOPT,1,3,0 !设置仅承受拉应力,KEYOPT(3)=0 R,1,306796E-8,543248E-8 !设置实常数,包括绳索截面积(306796E-8),初始应变(543248E-8) MP,EX,1,30E6 !定义材料的弹性模量(1号材料) MP,DENS,1,73E-5 !定义材料的密度(1号材料) N,1 ! 定义第1号节点 N,14,100 ! 定义第14号节点 FILL ! 均分填满第2号至第13号节点 E,1,2 !由节点1及节点2生成单元 EGEN,13,1,1 !依序复制生成13个单元 D,ALL,ALL ! 对所有节点施加固定约束 FINISH ! 前处理结束 /SOLU ! 进入求解模块,求解预应力引起的应力状态 SOLVE ! 求解 FINISH ! 退出求解模块 /POST1 ! 进入一般的后处理 ETABLE,STRS,LS,1 !针对LINK10单元,建立单元列表STRS,通过LS及特征号1来获得单元的轴向应力 *GET,STRSS,ELEM,13,ETAB,STRS !针对单元列表STRS, 提取13号单元的应力 FINISH ! 后处理结束 /POST26 ! 进入时间历程后处理,处理支反力 RFORCE,2,1,F,X !将1号节点上的x方向支反力提取,并存储到2号变量中 STORE ! 存储 *GET,FORCE,V ARI,2,EXTREM,VMAX !将2号变量的最大值赋给参数FORCE /SOLU ! 再次进入求解模块,模态分析 ANTYPE,MODAL ! 模态分析 MODOPT,SUBSP,3 ! 选择子空间迭代法,求3阶模态 MXPAND,3 ! 设定3阶模态扩展 PSTRES,ON ! 用于在模态分析中考虑预应力(还需在前面的静力分析中也同时打开) DDELE,2,UX,13 ! 删除从2号节点到13号节点上的UX约束 DDELE,2,UY,13 !删除从2号节点到13号节点上的UY约束 SOLVE !求解 *GET,FREQ1,MODE,1,FREQ ! 提取第1阶模态共振频率,并赋值给参数FREQ1 *GET,FREQ2,MODE,2,FREQ ! 提取第2阶模态共振频率,并赋值给参数FREQ2 *GET,FREQ3,MODE,3,FREQ ! 提取第3阶模态共振频率,并赋值给参数FREQ3 *STATUS !列出所有参数的实际内容

ansys命令流解释

对ansys主要命令的解释 本文给出了ansys主要命令的一些解释。 1, /PREP7 ! 加载前处理模块 2, /CLEAR,NOSTART ! 清除已有的数据, 不读入启动文件的设置(不加载初始化文件)初始化文件是用于记录用户和系统选项设置的文本文件 /CLEAR, START !清除系统中的所有数据,读入启动文件的设置/FILENAME, EX10.5 ! 定义工程文件名称 /TITLE, EX10.5 SOLID MODEL OF AN AXIAL BEARING ! 指定标题 4, F,2,FY,-1000 ! 在2号节点上施加沿着-Y方向大小为1000N 的集中力 6, FINISH ! 退出模块命令 7, /POST1 ! 加载后处理模块 8, PLDISP,2 ! 显示结构变形图,参数“2”表示用虚线绘制出原来结构的轮廓 9, ETABLE,STRS,LS,1 ! 用轴向应力SAXL的编号”LS,1”定义单元表STRS ETABLE, MFORX,SMISC,1 ! 以杆单元的轴力为内容, 建立单元表MFORX

ETABLE, SAXL, LS, 1 ! 以杆单元的轴向应力为内容, 建立单元表SAXL ETABLE, EPELAXL, LEPEL, 1 ! 以杆单元的轴向应变为内容, 建立单元表EPELAXL ETABLE,STRS_ST,LS,1 !以杆件的轴向应力“LS,1”为内容定义单元表STRS_ST ETABLE, STRS_CO, LS,1 !以杆件的轴向应力“LS,1”定义单元表STRS_CO ETABLE,STRSX,S,X ! 定义X方向的应力为单元表STRSX ETABLE,STRSY,S,Y ! 定义Y方向的应力为单元表STRSY *GET,STRSS_ST,ELEM,STEEL_E, ETAB, STRS_ST !从单元表STRS_ST中提取STEEL_E单元的应力结果,存入变量STRSS_ST; *GET, STRSS_CO,ELEM,COPPER_E,ETAB,STRS_CO”从单元表STRS_CO中提取COPPER_E单元的应力结果,存入变量STRSS_CO 10 FINISH !退出以前的模块 11, /CLEAR, START ! 清除系统中的所有数据,读入启动文件的设置 12 /UNITS, SI !申明采用国际单位制 14 /NUMBER, 2 !只显示编号, 不使用彩色 /NUMBER, 0 ! 显示编号, 并使用彩色 15 /SOLU ! 进入求解模块:定义力和位移边界条件,并求解 ANTYPE, STATIC ! 申明分析类型是静力分析(STATIC或者0)

流固耦合ansys命令流

new config fluid title 基于流固耦合作用下的双龙富水隧道稳定性研宄 set fluid off set log on set logfile yang 1 .log genzonradcylpOOOOpl 9.00 0p2 0 50 0 p3 0 0 8 size4 20 64 dim6 5 6 5 rat 1 1 I 1 group 围岩 gen zon cshell pOOOOpl 6.0 0 0 p2 0 50 0 p3 0 0 5.0 size 4 2064 dim 5.6 4.6 5.6 4.6 rat 1 1 1 1 group初期支护 gen zon cshell pO 0 0 0 pi 5.6 0 0 p2 0 50 0 p3 0 0 4.6 size 4 2064 dim 5.0 4.0 5.0 4.0 rat 1 1 1 1 group 二次衬砲fill group原岩 gen zon radcyl pOOOOpl 0 0 ?8.0 p2 0 50 0 p3 9.0 0 0 size 4 20 64 dim 3 6 3 6 rat 1 1 1 1 group围岩2 gen zon cshell pOOO Opl 0 0 -3.0 p2 0 50 0 p3 6.0 0 0 size 4 2064 dim 2.6 5.6 2.6 5.6 rat 1 1 1 1 group仰拱初期支护 gen zon cshell pO 0 0 Opl 0 0 -2.6 p2 0 50 0 p3 5.6 0 0 size 4 2064 dim2 5 2 5 rat 1 1 1 1 group仰拱二次衬砲fill group仰拱原岩 gen zone reflect normal -10 0 gen zone radtun pO 0 0 0 pi 45 0 0 p2 0 50 0 p3 0 0 20 size 3 20 3 12 dim 9 8 9 8 rat 1 1 1 1.1 group围岩3 gen zon reflect dip 0 ori 0 0 0 range x09y0 50z8 20 gen zon reflect dip 0 ori 0 0 0 range x 9 45 y 0 50 z 0 20 gen zon reflect dip 90 dd 270 ori 0 0 0 range x09y0 50z8 20 gen zon reflect dip 90 dd 270 ori 0 0 0 range x 0 9 y 0 50 z -8 -20 gen zon reflect dip 90 dd 270 ori 0 0 0 range x 9 45 y 0 50 z -20 20 save shuitun一model.sav model fl_iso prop perm 1.23e-9 poro 0.45 range z 4.5 20 prop perm 4.70e-10 poro 0.4 range z -20 4.5 set fl biot off ini fdensity le3 ini sat 1.0 ini food 2.0e9 ftens -le-3 ;假设围岩岩体符合mohr-coulomb本构模型,给围岩陚参数命令流如下, ;mohr-coulomb model model mohr def derive s_modl=E一modl/(2.0*(l .0+p—ratio 1)) b_modl=E_modl/(3.0*(1.0-2.0*p_ratiol)) s_mod2=E—mod2/(2.0*(l .0+p_ratio2)) bjmod2~E_mod2/(3.0*( 1.0-2.0*p_ratio2)) end

ansys命令流

第一天目标: 熟悉ANSYS基本关键字的含义k --> Keypoints关键点l --> Lines线a --> Area 面v --> Volumes体e --> Elements单元n --> Nodes节点cm --> component组元et --> element type单元类型mp --> material property材料属性r --> real constant实常数d --> DOF constraint约束f --> Force Load集中力sf --> Surface load on nodes 表面载荷bf --> Body Force on Nodes体载荷ic --> Initial Conditions初始条件第二天目标: 了解命令流的整体结构,掌握每个模块的标识!文件说明段/BATCH/TILE,test analysis!定义工作标题/FILENAME,test!定义工作文件名/PREP7!进入前处理模块标识!定义单元,材料属性,实常数段ET,1,SHELL63!指定单元类型ET,2,SOLID45!指定体单元MP,EX,1,2E8!指定弹性模量MP,PRXY,1, 0.3!输入泊松比MP,DENS,1, 7.8E3!输入材料密度R,1, 0.001!指定壳单元实常数-厚度......!建立模型K,1,0,0,,!定义关键点 K,2,50,0,,K,3,50,10,,K,4,10,10,,K,5,10,50,,K,6,0,50,,A,1,2,3,4,5,6,!由关键点生成面......!划分网格ESIZE,1,0,AMESH, 1......FINISH!前处理结束标识/SOLU!进入求解模块标识!施加约束和载荷DL,5,,ALLSFL,3,PRES,1000SFL,2,PRES, 1000......SOLVE!求解标识FINISH!求解模块结束标识/POST1!进入通用后处理器标识....../POST26!进入时间历程后处理器……/EXIT,SAVE!退出并存盘以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助/ANGLE!指定绕轴旋转视图/DIST!说明对视图进行缩放/DEVICE!设置图例的显示,如: 风格,字体等/REPLOT!重新显示当前图例/RESET!恢复缺省的图形设置/VIEW!设置观察方向/ZOOM!对图形显示窗口的某一区域进行缩放第三天生成关键点和线部分 1.生成关键点K,关键点编号,X坐标,Y坐标,Z坐标例:

ANSYS 结构稳态(静力)分析之经典实例-命令流格式

ANSYS 结构稳态(静力)分析之经典实例-命令流格式.txt两人之间的感情就像织毛衣,建立 的时候一针一线,小心而漫长,拆除的时候只要轻轻一拉。。。。/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元;Eight-node brick element ET,2,PLANE42 ! 2维平面结构;Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1;杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1;泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生 一个倒角圆,并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT

ansys实例命令流-弹塑性分析命令流

/FILNAME,Elastic-Plasitc,1 /TITLE, Elastic-Plasitc Analysis !前处理。 /PREP7 !**定义梁单元189。 ET,1,BEAM189 !定义单元。 !**梁截面1。 SECTYPE, 1, BEAM, HREC, , 0 !定义梁截面。SECOFFSET, CENT SECDATA,50,100,6,6,6,6,0,0,0,0 !定义梁截面完成。 !**定义材料。 MPTEMP,,,,,,,, !定义弹塑性材料模型。MPTEMP,1,0 MPDATA,EX,1,,2.05e5 MPDATA,PRXY,1,,0.3 TB,BISO,1,1,2, TBTEMP,0 TBDATA,,150,18600,,,, !定义弹塑性材料模型。!**建立几何模型。 K,1, , , , K,2 ,900, K,3 ,,50 LSTR, 1, 2 !**网格划分。 FLST,5,1,4,ORDE,1 !定义网格密度。FITEM,5,1 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,50, , , , ,1 !定义网格密度完成。CM,_Y,LINE !网格划分。 LSEL, , , , 1 CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,1, ,1, , 3, ,1 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 LMESH, 1 !网格划分完成。 !施加载荷及求解。 FINISH /SOL

!**施加约束。 FLST,2,1,3,ORDE,1 !施加约束。FITEM,2,1 /GO DK,P51X, , , ,0,UX,UY,UZ,ROTX, , , FLST,2,1,3,ORDE,1 FITEM,2,2 /GO DK,P51X, , , ,0,UY,UZ,ROTX, , , , !施加约束完成。 !**加载。 FLST,2,50,2,ORDE,2 FITEM,2,1 FITEM,2,-50 SFBEAM,P51X,1,PRES,100, , , , , , LSWRITE,1, !定义载荷步1完成。FLST,2,50,2,ORDE,2 !定义载荷步2。FITEM,2,1 FITEM,2,-50 SFEDELE,P51X,1,PRES LSWRITE,2, !定义载荷步2完成。!设定求解步并求解。 LSSOLVE,1,2,1,

ansys命令流最全详细介绍二

三 生成关键点和线部分 1.生成关键点 K,关键点编号,X坐标,Y坐标,Z坐标 例:K,1,0,0,0 2.在激活坐标系生成直线 LSTR,关键点P1,关键点P2 例LSTR,1,2 3.在两个关键点之间连线 L,关键点P1,关键点P2 例L,1,2 注:此命令会随当前的激活坐标系不同而生成直线或弧线 4.由三个关键点生成弧线 LARC,关键点P1,关键点P2,关键点PC,半径RAD 例LARC,1,3,2,0.05 注:关键点PC是用来控制弧线的凹向 5.通过圆心半径生成圆弧

CIRCLE,关键点圆心,半径RAD,,,,圆弧段数NSEG 例:CIRCLE,1,0.05,,,,4 6.通过关键点生成样条线 BSPLIN,关键点P1,关键点P2,关键点P3,关键点P4,关键点P5,关键点P6 例:BSPLIN,1,2,3,4,5,6 7.生成倒角线 LFILLT,线NL1,线NL2,倒角半径RAD 例LFILLT,1,2,0.005 8.通过关键点生成面 A,关键点P1,关键点P2,关键点P3,关键点P4,关键点P5,关键点P6,P7,P8... 例:A,1,2,3,4 9.通过线生成面 AL,线L1,线L2,线L3,线L4,线L5,线L6,线L7,线L8,线L9,线L10 例:AL,5,6,7,8 10.通过线的滑移生成面

ASKIN,线NL1,线NL2,线NL3,线NL4,线NL5,线NL6,线NL7,线NL8,线NL9 例:ASKIN,1,4,5,6,7,8 注:线1为滑移的导向线 四 目标:掌握常用的实体-面的生成 生成矩形面 1.通过矩形角上定位点生成面 BLC4,定位点X方向坐标XCORNER,定位点Y方向坐标YCORNER,矩形宽度WIDTH,矩形高度HEIGHT,矩形深度DEPTH 例:BLC4,0,0,5,3,0 2.通过矩形中心定位点生成面 BLC5,定位点X方向坐标XCENTER,定位点Y方向坐标YCENTER,矩形宽度WIDTH,矩形高度HEIGHT,矩形深度DEPTH 注:与上条命令的不同就在于矩形的定位点不一样 例:BLC5,2.5,1.5,5,3,0 3.通过在工作平面定义矩形X.Y坐标生成面 RECTNG,矩形左边界X坐标X1,矩形右边界X坐标X2,矩形下边界Y

ansys14workbench血管流固耦合分析实例

Ansys14 workbench血管流固耦合实例 根据收集得一些资料,进行学习后,试着做了这个ansys14workbench得血管流固耦合模拟,感觉能够耦合上,仅就是熟悉流固耦合分析过程,不一定正确,仅供参考,希望大家多讨论。谢谢! 1、先在proe5中建立血管与血液流体区得模型(两者装配起来),或者直接在workbench中建模。 图1 模型图 2、新建工程。在workbench中toolbox中选custom system,双击FSI: FluidFlow(fluent)->static structure、 图2 计算工程 3、修改engineering data,因为系统缺省材料就是钢,需要构建血管材料,如图3所示。先复制steel,而后修改密度1150kg/m3,杨氏模量 4、5e8Pa,泊松比0、3,重新命名,最后在主菜单中点击“update project”保存、

图3 修改工程材料 4、模型导入,进入gemetry模块,import外部模型文件。 图4 模型导入图 5、进入FLUENT网格划分。 在workbench工程视图中得Mesh上点击右键,选择Edit…,如图5所示,进入网格划分meshing界面,如图6所示。我们这里需要去掉血管部分,只保留血液几何。

图5 进入网格划分

图6 禁用血管模型 6、设置网格方法。 默认就是采用ICEM CFD进行网格划分,设置方式如图7所示,截面圆弧边分为12份,纵截面得边均分为10份,网格结果如图8所示。另外在这个界面中要设置边界得几何面,如inlet、outlet、symmetry 图7 设置网格划分方式 图8 最终出网格

ANSYS命令流(入门必备)

ANSYS命令集 /EXIT,Slab,Fname,Ext,Dir Slab=ALL 保存所有资料 Slab=NOSA VE所有更改资料不保存 Slab=MODEL保存实体模型,有限元 模型,负载的资料(系统默认)例:/EXIT,ALL -------------------------------------------------------- /FILNAM,Fname Fname=工作文件名称,不要扩展名例:/FILNAM,Sanpangzi --------------------------------------------------------/SA VE,Fname,Ext,Dir 保存目前所有的Datebase资料,即 更新Jobname.db --------------------------------------------------------/RESUME,Fname,Ext,Dir,NOPAR 回到最后SA VE时的Datebase 状态 --------------------------------------------------------/CLEAR 清除所有Datebase资料 -------------------------------------------------------- LOCAL,KCN,KCS,XC,YC,ZC,THXY,THYZ,THZX,PAR1,PAR2 定义区域坐标系统 KCN 区域坐标系统代号,大于10的任何号码

KCS=0,1,20=笛卡儿坐标1=圆柱坐标2=球面坐标XC,YC,ZC 该区域坐标原点与整体坐标原点的关系 THXY,THYZ,THZX 该区域坐标与整体坐标XYZ轴的关系例:LOCAL,11,1,1,1,0 -------------------------------------------------------- CSYS,0,1,2声明当前坐标系统 例:CSYS,0 -------------------------------------------------------- /UNITS,LABEL 声明系统分析时所用的单位 LABEL=SI (米,千克,秒) LABEL=CGS (厘米,克,秒) LABEL=BFT (英尺) LABEL=BIN (英寸) 例:LABEL,SI -------------------------------------------------------- /PREP7进入通用前处理器-------------------------------------------------------- N,NODE,X,Y,Z,THXY,THYZ,THZX 定义节点NODE 节点号码X,Y,Z 节点在当前坐标系中位置 例:N,1,2,3,4 -------------------------------------------------------- NDELE,NODE1,NODE2,NINC 删除已建立的节点

几个ansys经典实例(长见识)

平面问题斜支座的处理 如图5-7所示,为一个带斜支座的平面应力结构,其中位置2及3处为固定约束,位置4处为一个45o的斜支座,试用一个4节点矩形单元分析该结构的位移场。 (a)平面结构(b)有限元分析模型 图5-7 带斜支座的平面结构 基于ANSYS平台,分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。 (7) 模型加约束 左边施加X,Y方向的位移约束 ANSYS Main Menu: Solution →Define Loads →Apply →-Structural→Displacement On Nodes →选取2,3号节点→OK →Lab2: All DOF(施加X,Y方向的位移约束) →OK 以下提供两种方法处理斜支座问题,使用时选择一种方法。 ?采用约束方程来处理斜支座 ANSYS Main Menu:Preprocessor →Coupling/ Ceqn →Constraint Eqn :Const :0, NODE1:4, Lab1: UX,C1:1,NODE2:4,Lab2:UY,C2:1→OK 或者?采用斜支座的局部坐标来施加位移约束 ANSYS Utility Menu:WorkPlane →Local Coordinate System →Create local system →At specified LOC + →单击图形中的任意一点→OK →XC、YC、ZC分别设定为2,0,0,THXY:45 →OK ANSYS Main Menu:Preprocessor →modeling →Move / Modify →Rotate Node CS →To active CS → 选择4号节点 ANSYS Main Menu:Solution →Define Loads →Apply →Structural →Displacement On Nodes →选取4号节点→OK →选择Lab2:UY(施加Y方向的位移约束) →OK 命令流; !---方法1 begin----以下的一条命令为采用约束方程的方式对斜支座进行处理 CE,1,0,4,UX,1,4,UY,-1 !建立约束方程(No.1): 0=node4_UX*1+node_UY*(-1) !---方法1 end --- !--- 方法2 begin --以下三条命令为定义局部坐标系,进行旋转,施加位移约束 !local,11,0,2,0,0,45 !在4号节点建立局部坐标系 !nrotat, 4 !将4号节点坐标系旋转为与局部坐标系相同 !D,4,UY !在局部坐标下添加位移约束 !--- 方法2 end

几个ansys流固耦合的例子

一般说来,ANSYS的流固耦合主要有4种方式: 1,sequential 这需要用户进行APDL编程进行流固耦合 sequentia指的是顺序耦合 以采用MpCCI为例,你可以利用ANSYS和一个第三方CFD产品执行流固耦合分析。在这个方法中,基于网格的平行代码耦合界面(MpCCI) 将ANSYS和CFD程序耦合起来。即使网格上存在差别,MpCCI也能够实现流固界面的数据转换。ANSYS CD中包含有MpCCI库和一个相关实例。关于该方法的详细信息,参见ANSYS Coupled-Field Analysis Guide中的Sequential Couplin 2,FSI solver 流固耦合的设置过程非常简单,推荐你使用这种方式 3,multi-field solver 这是FSI solver的扩展,你可以使用它实现流体,结构,热,电磁等的耦合 4,直接采用特殊的单元进行直接耦合,耦合计算直接发生在单元刚度矩阵 一个流固耦合的例子 length=2 width=3 height=2 /prep7 et,1,63 et,2,30 !选用FLUID30单元,用于流固耦合问题 r,1,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

ANSYS最常用命令流+中文注释(超级大全)

ANSYS最常用命令流+中文注释 VSBV, NV1, NV2, SEPO, KEEP1, KEEP2 —Subtracts volumes from volumes,用于2个solid相减操作,最终目的是要nv1-nv2=?通过后面的参数设置,可以得到很多种情况:sepo项是2个体的边界情况,当缺省的时候,是表示2个体相减后,其边界是公用的,当为sepo的时候,表示相减后,2个体有各自的独立边界。keep1与keep2是询问相减后,保留哪个体?当第一个为keep时,保留nv1,都缺省的时候,操作结果最终只有一个体,比如:vsbv,1,2,sepo,,keep,表示执行1-2的操作,结果是保留体2,体1被删除,还有一个1-2的结果体,现在一共是2个体(即1-2与2),且都各自有自己的边界。如vsbv,1,2,,keep,,则为1-2后,剩下体1和体1-2,且2个体在边界处公用。同理,将v换成a 及l是对面和线进行减操作! mp,lab, mat, co, c1,…….c4 定义材料号及特性 lab: 待定义的特性项目(ex,alpx,reft,prxy,nuxy,gxy,mu,dens) ex: 弹性模量 nuxy: 小泊松比 alpx: 热膨胀系数 reft: 参考温度 reft: 参考温度 prxy: 主泊松比 gxy: 剪切模量 mu: 摩擦系数 dens: 质量密度 mat: 材料编号(缺省为当前材料号) co: 材料特性值,或材料之特性,温度曲线中的常数项 c1-c4: 材料的特性-温度曲线中1次项,2次项,3次项,4次项的系数 定义DP材料: 首先要定义EX和泊松比:MP,EX,MA T,…… MP,NUXY,MAT,…… 定义DP材料单元表(这里不考虑温度):TB,DP,MA T 进入单元表并编辑添加单元表:TBDATA,1,C TBDATA,2,ψ TBDATA,3,…… 如定义:EX=1E8,NUXY=0.3,C=27,ψ=45的命令如下:MP,EX,1,1E8 MP,NUXY,1,0.3 TB,DP,1 TBDATA,1,27 TBDATA,2,45这里要注意的是,在前处理的最初,要将角度单位转化到“度”,即命令:*afun,deg VSEL, Type, Item, Comp, VMIN, VMAX, VINC, KSWP Type,是选择的方式,有选择(s),补选(a),不选(u),全选(all)、反选(inv)等,其余方式不常用 Item, Comp 是选取的原则以及下面的子项 如volu 就是根据实体编号选择, loc 就是根据坐标选取,它的comp就可以是实体的某方向坐标! 其余还有材料类型、实常数等 MIN, VMAX, VINC,这个就不必说了吧! ,例:vsel,s,volu,,14 vsel,a,volu,,17,23,2 上面的命令选中了实体编号为14,17,19,21,23的五个实体 VDELE, NV1, NV2, NINC, KSWP: 删除未分网格的体 nv1:初始体号 nv2:最终的体号 ninc:体号之间的间隔 kswp=0:只删除体 kswp=1:删除体及组成关键点,线面 如果nv1=all,则nv2,ninc不起作用 其后面常常跟着一条显示命令VPLO,或aplo,nplo,这个湿没有参数的命令,输入后直接回车,就可以显示刚刚选择了的体、面或节点,很实用的哦! Nsel, type, item, comp, vmin, vmax, vinc, kabs 选择一组节点为下一步做准备 Type: S: 选择一组新节点(缺省) R: 在当前组中再选择 A: 再选一组附加于当前组 U: 在当前组中不选一部分 All: 恢复为选中所有 None: 全不选 Inve: 反向选择 Stat: 显示当前选择状态 Item: loc: 坐标 node: 节点号

ANSYS_结构稳态(静力)分析之经典实例-命令流格式

/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元;Eight-node brick element ET,2,PLANE42 ! 2维平面结构;Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1;杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1;泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生一个倒角圆,并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT ! Line numbers off !

相关主题
相关文档 最新文档