dong 发表于 2007-7-23 14:49:11

流固耦合的例子！不运行，不知道原因

*
*****************************************************************************************
*本例题是一个流固耦合模态分析的算例；
*
*模型是一个液体中浮体结构的流固耦合频域分析；
*
*流固耦合模态分析的流体必须采用势流体，在建立势流体与结构的耦合模型时如果
*使势流体的节点与相邻结构的节点拟和，ADINA就可以自动建立流固耦合界面，如果节点
*不拟合则需要手动建立流固耦合界面；
*
*对于此类问题应提前根据流体与液体的密度计算好吃水线深度，在建立几何模型时把相互位置建好。
*******************************************************************************************
*--- by ADINA: AUI version 8.4.0 ---*
*
DATABASE NEW SAVE=NO PROMPT=NO
FEPROGRAM ADINA
CONTROL FILEVERSION=V84
*
FEPROGRAM PROGRAM=ADINA
***************************************************************
*建立几何模型
***************************************************************
*
COORDINATES POINT SYSTEM=0
@CLEAR
1 0.00000000000000 0.00000000000000 0.00000000000000 0
2 3.50000000000000 0.00000000000000 0.00000000000000 0
3 6.50000000000000 0.00000000000000 0.00000000000000 0
4 10.0000000000000 0.00000000000000 0.00000000000000 0
@
*
LINE STRAIGHT NAME=1 P1=1 P2=2
*
LINE STRAIGHT NAME=2 P1=2 P2=3
*
LINE STRAIGHT NAME=3 P1=3 P2=4
*
SUBDIVIDE LINE NAME=1 MODE=DIVISIONS NDIV=4 RATIO=1.00000000000000,
   PROGRESS=GEOMETRIC CBIAS=NO
@CLEAR
2
3
@
*
SURFACE EXTRUDED NAME=1 LINE=1 DX=0.00000000000000,
   DY=4.00000000000000 DZ=0.00000000000000 SYSTEM=0 PCOINCID=YES,
   PTOLERAN=1.00000000000000E-05 NDIV=8 OPTION=VECTOR
@CLEAR
2
3
@
*
SURFACE EXTRUDED NAME=4 LINE=6 DX=0.00000000000000,
   DY=2.00000000000000 DZ=0.00000000000000 SYSTEM=0 PCOINCID=YES,
   PTOLERAN=1.00000000000000E-05 NDIV=2 OPTION=VECTOR
@CLEAR
8
10
@
*
SURFACE EXTRUDED NAME=7 LINE=13 DX=0.00000000000000,
   DY=4.00000000000000 DZ=0.00000000000000 SYSTEM=0 PCOINCID=YES,
   PTOLERAN=1.00000000000000E-05 NDIV=4 OPTION=VECTOR
@CLEAR
15
17
@
*
VOLUME EXTRUDED NAME=1 SURFACE=1 DX=0.00000000000000,
   DY=0.00000000000000 DZ=4.00000000000000 SYSTEM=0 PCOINCID=YES,
   PTOLERAN=1.00000000000000E-05 NDIV=4 OPTION=VECTOR,
   RATIO=1.00000000000000 PROGRESS=GEOMETRIC CBIAS=NO
@CLEAR
2
3
4
5
6
7
8
9
@
*
VOLUME EXTRUDED NAME=10 SURFACE=14 DX=0.00000000000000,
   DY=0.00000000000000 DZ=1.00000000000000 SYSTEM=0 PCOINCID=YES,
   PTOLERAN=1.00000000000000E-05 NDIV=3 OPTION=VECTOR,
   RATIO=1.00000000000000 PROGRESS=GEOMETRIC CBIAS=NO
@CLEAR
18
22
32
29
26
36
39
42
@
*
VOLUME EXTRUDED NAME=19 SURFACE=62 DX=0.00000000000000,
   DY=0.00000000000000 DZ=1.00000000000000 SYSTEM=0 PCOINCID=YES,
   PTOLERAN=1.00000000000000E-05 NDIV=3 OPTION=VECTOR,
   RATIO=1.00000000000000 PROGRESS=GEOMETRIC CBIAS=NO
***************************************************************
*定义材料，对于本例题浮体的密度是液体密度的0.5倍，在做几何模型时
*浮体的吃水线恰好在高度1/2的位置。
***************************************************************
*
MATERIAL ELASTIC NAME=1 E=2.09000000000000E+09 NU=0.300000000000000,
   DENSITY=500.000000000000 ALPHA=0.00000000000000 MDESCRIP='NONE'
*
MATERIAL FLUID NAME=2 K=3.00000000000000E+09 DENSITY=1000.00000000000,
   GRAVITY=0.00000000000000 X0=0.00000000000000 Y0=0.00000000000000,
   Z0=0.00000000000000 MDESCRIP='NONE'
***************************************************************
*定义自由液面边界条件，对于浮体结构的模态分析必须用自由液面条件
***************************************************************
*
POTENTIAL-IN FREE-SURFACE NAME=1 GTYPE=SURFACES
@CLEAR
55 0
51 0
47 0
65 0
59 0
75 0
72 0
69 0
@
***************************************************************
*施加重力荷载
***************************************************************
*
LOAD MASS-PROPORTIONAL NAME=1 MAGNITUD=9.80000000000000,
   AX=0.00000000000000 AY=0.00000000000000 AZ=-1.00000000000000,
   INTERPRE=BODY-FORCE
*
APPLY-LOAD BODY=0
@CLEAR
1'MASS-PROPORTIONAL' 1'MODEL' 0 0 1 0.00000000000000 0 -1 0 0 0,
'NO' 0.00000000000000 0.00000000000000 1 0
@
***************************************************************
*定义单元组
***************************************************************
*
EGROUP THREEDSOLID NAME=1 DISPLACE=DEFAULT STRAINS=DEFAULT MATERIAL=1,
   RSINT=DEFAULT TINT=DEFAULT RESULTS=STRESSES DEGEN=YES FORMULAT=0,
   STRESSRE=GLOBAL INITIALS=NONE FRACTUR=NO CMASS=DEFAULT,
   STRAIN-F=0 UL-FORMU=DEFAULT LVUS1=0 LVUS2=0 SED=NO RUPTURE=ADINA,
   INCOMPAT=DEFAULT TIME-OFF=0.00000000000000 POROUS=NO,
   WTMC=1.00000000000000 OPTION=NONE DESCRIPT='NONE' PRINT=DEFAULT,
   SAVE=DEFAULT TBIRTH=0.00000000000000 TDEATH=0.00000000000000
*
EGROUP FLUID3 NAME=2 DISPLACE=DEFAULT IPO=0 MATERIAL=2 RSINT=DEFAULT,
   TINT=DEFAULT RESULTS=PRESSURES DEGEN=DEFAUL FORMULAT=2,
   CMASS=DEFAULT DESCRIPT='NONE' PRINT=DEFAULT SAVE=DEFAULT,
   TBIRTH=0.00000000000000 TDEATH=0.00000000000000
***************************************************************
*划分网格
***************************************************************
*
GVOLUME NODES=8 PATTERN=0 NCOINCID=BOUNDARIES NCFACE=123456 NCEDGE=,
'123456789ABC' NCVERTEX=12345678 NCTOLERA=1.00000000000000E-05,
   SUBSTRUC=0 GROUP=1 MESHING=MAPPED PREFSHAP=AUTOMATIC,
   DEGENERA=YES COLLAPSE=NO MIDNODES=CURVED METHOD=DELAUNAY,
   BOUNDARY=ADVFRONT
@CLEAR
14
19
@
*
GVOLUME NODES=8 PATTERN=0 NCOINCID=BOUNDARIES NCFACE=123456 NCEDGE=,
'123456789ABC' NCVERTEX=12345678 NCTOLERA=1.00000000000000E-05,
   SUBSTRUC=0 GROUP=2 MESHING=MAPPED PREFSHAP=AUTOMATIC,
   DEGENERA=YES COLLAPSE=NO MIDNODES=CURVED METHOD=DELAUNAY,
   BOUNDARY=ADVFRONT
@CLEAR
1
2
3
4
5
6
7
8
9
10
11
12
13
15
16
17
18
@
***************************************************************
*求解控制
***************************************************************
*
MASTER ANALYSIS=FREQUENCIES MODEX=EXECUTE TSTART=0.00000000000000,
   IDOF=0 OVALIZAT=NONE FLUIDPOT=AUTOMATIC CYCLICPA=1 IPOSIT=STOP,
   REACTION=YES INITIALS=NO FSINTERA=NO IRINT=DEFAULT CMASS=NO,
   SHELLNDO=AUTOMATIC AUTOMATI=OFF SOLVER=SPARSE,
   CONTACT-=CONSTRAINT-FUNCTION TRELEASE=0.00000000000000,
   RESTART-=NO FRACTURE=NO LOAD-CAS=NO LOAD-PEN=NO MAXSOLME=0,
   MTOTM=2 RECL=3000 SINGULAR=YES STIFFNES=0.000100000000000000,
   MAP-OUTP=NONE MAP-FORM=NO NODAL-DE='' POROUS-C=NO ADAPTIVE=0,
   ZOOM-LAB=1 AXIS-CYC=0 PERIODIC=NO VECTOR-S=GEOMETRY EPSI-FIR=NO,
   STABILIZ=NO STABFACT=1.00000000000000E-12 RESULTS=PORTHOLE,
   FEFCORR=NO BOLTSTEP=1 EXTEND-S=YES CONVERT-=NO DEGEN=YES
*
FREQUENCIES METHOD=LANCZOS-ITERATION NEIGEN=20 NMODE=0 IPRINT=NO,
   RIGID-BO=NO RSHIFT=0.00000000000000 CUTOFF=1.00000000000000E+08,
   NITEMM=DEFAULT STURM-CH=NO TOLERANC=DEFAULT INTERVAL=NO,
   FMIN=0.00000000000000 FMAX=DEFAULT MODALSTR=NO STATIC=NO,
   NSHIFT=AUTO NSHIFT-B=50
这是adina网站提供的流固耦合的例子，不知道还缺什么东东，不能求解。我用的是adina8.2，望大虾给看看！万分感谢！

duanxc314 发表于 2007-7-23 16:54:21

"CONTROL FILEVERSION=V84"很显然是用adina8.4做的，这种软件一般很难向下兼容的。

huahua2005 发表于 2007-7-23 17:36:08

所以只能根据命令流自己来操作，或者试着把命令流中的与8.2版本不同的部分去掉

[ 本帖最后由 huahua2005 于 2007-7-23 17:38 编辑 ]

dong 发表于 2007-7-25 17:23:44

谢谢！

谢谢！本人是老菜，能否提供一点帮助？帮忙修改一下！哪怕给点详细建议也好啊!万分感谢！:handshake

武汉游侠 发表于 2007-7-26 19:48:51

很简单的一个例子，你照着思路做一遍就是的
页: [1]
查看完整版本: 流固耦合的例子！不运行，不知道原因