$-----------------------------------------------------------------------------
$
$ Example provided by Marcus Timgren (DYNAmore Nordic)
$
$ E-Mail: info@dynamore.de
$ Web: http://www.dynamore.de
$
$ Copyright, 2015 DYNAmore GmbH
$ Copying for non-commercial usage allowed if
$ copy bears this notice completely.
$
$X------------------------------------------------------------------------------
$X
$X 1. Run file main_fsi.k and select memory=5m.
$X Requires LS-DYNA MPP R8.0.0 (or higher) with double precision.
$X------------------------------------------------------------------------------
$X------------------------------------------------------------------------------
$X
$X
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ main_fsi.k $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$-------------------------------------------------------------------------------
$
*KEYWORD
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ INCLUDE SOLID AND FLUID DECK
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*INCLUDE
struc.k
*INCLUDE
i.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ BOUNDARY CONDITIONS
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_BOUNDARY_FSI
$# pid
4
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_CONTROL_FSI
$# owc bt dt idc lcidsf
0 0.01.00000E28 0.25 0
*END
$
$
$-------------------------------------------------------------------------------
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ Solid Deck $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$-------------------------------------------------------------------------------
*KEYWORD
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARAMETERS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
$# prmr1 val1 prmr2 val2 prmr3 val3 prmr4 val4
R t_ends100.0
$--- Struc
$
Rrho_struc1.20
R e_struc2e11
R nu_struc0.3
Rdtm_struc0.05
$--- Struc
$
$--- Struc
$
$--- Struc
$
*INCLUDE
mesh_struc.k
*INCLUDE
boundary_condition.k
*INCLUDE
implicit.k
$--- Struc
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*CONTROL_TERMINATION
$# endtim endcyc dtmin endeng endmas
&t_ends 0 0.0 0.01.000000E8
*MAT_RIGID
$# mid ro e pr n couple m alias
1&rho_struc&e_struc &nu_struc 0.0 0.0 0.0
$# cmo con1 con2
0.0 0 0
$#lco or a1 a2 a3 v1 v2 v3
0.0 0.0 0.0 0.0 0.0 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARTS/ SECTION/ MATERIAL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PART
$# title
$# pid secid mid eosid hgid grav adpopt tmid
1 1 1 0 0 0 0 0
*SECTION_SHELL
$# secid elform shrf nip propt qr/irid icomp setyp
1 12 1.0 2 1.0 0 0 1
$# t1 t2 t3 t4 nloc marea idof edgset
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0
*END
$
$-------------------------------------------------------------------------------
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ Fluid Deck $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$-------------------------------------------------------------------------------
*KEYWORD
*TITLE
$# title
ICFD Cylindrical Flow
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ PARAMETERS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
R t_end100.0
R dt_plot1.00
$
$--- Fluid
$
R v_inlet1.0
Rrho_fluid1.0
R mu_fluid0.005
R dt_fluid0.050
$
$--- Fluid
*TITLE
$# title
ICFD Cylinder Flow
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ INCLUDE SURFACE MESH
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*INCLUDE
mesh.k
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ DATABASE CARD FOR RESULTS OUTPUT
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BINARY_D3PLOT
$# dt lcdt beam npltc psetid
&dt_plot 0 0 0 0
$# ioopt
0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ LOAD CURVES FOR THE BOUNDARY CONDITIONS
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DEFINE_CURVE_TITLE
Velocity inlet
$# lcid sidr sfa sfo offa offo dattyp lcint
1 0 1.0&v_inlet 0.0 0.0 0 0
$# a1 o1
0.0 1.0
10000.0 1.0
*DEFINE_CURVE_TITLE
Pressure outlet
$# lcid sidr sfa sfo offa offo dattyp lcint
2 0 1.0 1.0 0.0 0.0 0 0
$# a1 o1
0.0 0.0
10000.0 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD BOUNDARY/INITIAL CONDITIONS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_BOUNDARY_FREESLIP
$# pid
3
*ICFD_BOUNDARY_NONSLIP
$# pid
4
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# pid dof vad lcid sf vid death birth
1 1 1 1 1.0 01.00000E28 0.0
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# pid dof vad lcid sf vid death birth
1 2 1 2 1.0 01.00000E28 0.0
*ICFD_BOUNDARY_PRESCRIBED_PRE
$# pid lcid sf death birth
2 2 1.01.00000E28 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD CONTROL CARDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_CONTROL_TIME
$# ttm dt cfl lcidsf dtmin dtmax
&t_end &dt_fluid 1.0 0 0.0 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ DATABASE (OUTPUT) $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_DATABASE_DRAG
$# pid cpid dtout
4 0 0.0
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD PARTS/ SECTION/ MATERIAL $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_MAT
$# mid flg ro vis st
1 1&rho_fluid&mu_fluid 0.0
*ICFD_PART
$# pid secid mid
1 1 1
*ICFD_PART
$# pid secid mid
2 1 1
*ICFD_PART
$# pid secid mid
3 1 1
*ICFD_PART
$# pid secid mid
4 1 1
*ICFD_PART_VOL
$# pid secid mid
10 1 1
$# spid1 spid2 spid3 spid4 spid5 spid6 spid7 spid8
1 2 3 4 0 0 0 0
*ICFD_SECTION
$# sid
1
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$ $
$ ICFD MESH KEYWORDS $
$ $
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_BL
$# pid nelth blth blfe blst
4 2 0.0 0.0 0
*MESH_VOLUME
$# volid
1
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
1 2 3 4 0 0 0 0
*END