銆愯揪灏旀暣鐞嗐�慉NSYS娴佸浐鑰﹀悎鍒嗘瀽瀹炰緥鍛戒护娴?- 鐧惧害鏂囧簱

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

$

*PART

Part 2 for Mat 2 and Elem Type 1 2 2 2 0 0 0 0 $

*PART

Part 3 for Mat 3 and Elem Type 1 3 2 3 0 0 0 0 *PART

Part 3 for Mat 3 and Elem Type 1 4 2 4 0 0 0 0 $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$

$ RIGID BOUNDRIES $

$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $

*CONTACT_SURFACE_TO_SURFACE 1 2 3 3 0 0 0 0

0.0000 0.0000 0.0000 0.0000 0.0000 0 0.00000.1000E+08 1.0000 0.0001 0.0000 0.0000 1.0000 1.0000 1.0000 1.0000 *CONTACT_SURFACE_TO_SURFACE 1 3 3 3 0 0 0 0

0.0000 0.0000 0.0000 0.0000 0.0000 0 0.00000.1000E+08 1.0000 0.0001 0.0000 0.0000 1.0000 1.0000 1.0000 1.0000 *CONTACT_SURFACE_TO_SURFACE 1 4 3 3 0 0 0 0

0.0000 0.0000 0.0000 0.0000 0.0000 0 0.00000.1000E+08 1.0000 0.0001 0.0000 0.0000 1.0000 1.0000 1.0000 1.0000 $

*DEFINE_CURVE

1 0 1.000 1.000 0.000 0.000

0.000000000000E+00 1.000000000000E-04 1.000000000000E+05 1.000000000000E-04

*BOUNDARY_PRESCRIBED_MOTION_RIGID 3 2 0 1 -1.00 0 0.000 5.00e3

*BOUNDARY_PRESCRIBED_MOTION_RIGID 4 2 0 1 -1.00 0 5.00e3 0.000 $

通常我们处理液面晃动采用上述方法，LS-DYNA求解器在下一个版本LS970中在MESHFREE功能上增加了更为稳定和高效的EFG

技术，在解决此类问题上将给予我们更大的灵活性和更多的选择。

41

DareDoc

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

我是通过设置流体单元的关键字还实现接触与非接触的，边界条件的话，我是用的底部固结，当然也可以在地脚锚栓处固结。这个差别都不是很大的

有人说ansys公司技术人员建议对流固耦合模态分析使用fluid80单元，然后看了以下ansys帮助中说 The fluid80 element is particularly well suited for calculating hydrostatic pressures and fluid/solid interactions

脱硫塔于浆液耦合的分析

上次做的脱硫塔于浆液耦合的分析中间有事就一直停下来没做，前两天又拾起来接着做。现在遇到了一些问题，我也查不出是什么原因造成的：就是我在用sf家族命令标记耦合界面时，出现了这样的错误

*** ERROR *** CP= 29.297 TIME= 22:15:44 Number of FLUID/SOLID interfaces do not match for FSIN boundary condition.

*** ERROR *** CP= 29.641 TIME= 22:15:48 Error during fsi setup - Skipping the remaining FSI analysis. 下面是我标记耦合界面的命令流： !对塔体底部加约束 lsel,s,,,54,63 nsll,s,1 d,all,all !求解

!*** FSI interfaces !!fluid side

asel,s,,,16,37,3 allsel,below,area nsla,s,1 sf,all,fsi,1

!structural side

asel,s,,,39,78 allsel,below,area nsla,s,1 sf,all,fsin,1

allsel /solu

fldata,solu,flow,1

42

DareDoc

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

fldata,solu,ale,1 fldata,solu,turb,1 fldata,solu,tran,1 fldata,iter,exec,50 !Fluid Properties

fldata7,prot,dens,constant fldata7,prot,visc,constant

fldata8,nomi,dens,1000. ! 1000 kg/m3 for density - water fldata8,nomi,visc,4.6E-4 ! 4.6E-4 kg-s/m (viscosity of water) fldata8,conv,pres,1.E-8 ! Tighten pressure equation convergence

fsan,on !FSI analysis on

fsco,all,1.0e-1 !FSI convergence criteria fsre,all,0.5 !FSI relaxation

fsou,1 !FSI output frequency fsit,5 !Number of stagger loops fsti,0.05 !FSI end time fsdt,0.01 !FSI delta time

fsor,fluid !FSI order of solution

fstr,solid, !Transient fluid for FSI analysis fstr,fluid, !Transient structural for FSI analysis deltim,0.01

fldata,time,step,0.01 save

流固耦合瞬态分析的例子

/batch a=1.0 b=1.0

gap=0.8*a/2 !gap=0.4 thick=0.01*a !thick=0.01 c=0.5*a !c=0.5

/PREP7

smrt,off !! 关闭智能化网格功能 et,1,142 !142 三维流体单元 et,2,45 et,3,141 et,4,42

keyopt,1,4,1 keyopt,3,4,1 msha,0

DareDoc

43

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

mshk,1

rect,0,gap,0,b/2-thick

rect,0,gap,b/2-thick,b/2+thick rect,0,gap,b/2+thick,b rect,gap,a/2,0,b/2-thick

rect,gap,a/2,b/2-thick,b/2+thick rect,gap,a/2,b/2+thick,b nummrg,kp numcmp,all

nx1=16-4

lesi,1,,,nx1,-1,1 lesi,3,,,nx1,-1,1 lesi,6,,,nx1,-1,1 lesi,9,,,nx1,-1,1

nx2=6-2

lesi,11,,,nx2,1,1 lesi,13,,,nx2,1/1,1 lesi,15,,,nx2,1/1,1 lesi,17,,,nx2,1/1,1

ny1=2

lesi,7,,,ny1,1,1 lesi,5,,,ny1,1,1 lesi,14,,,ny1,1,1

ny2=20-4

lesi,4,,,ny2,-2,1 lesi,2,,,ny2,-2,1 lesi,12,,,ny2,-2,1 lesi,10,,,ny2,-2,1 lesi,8,,,ny2,-2,1 lesi,16,,,ny2,-2,1

asel,s,,,2 aatt,,,4 asel,inve aatt,,,3 alls

amesh,all agen,2,2 aclear,2 adel,2,,,1

DareDoc

44