Intermediate: FSI - 2 Flaps
This LS-DYNA simulation shows the fluid-structure interaction (FSI) of two flaps exposed to a flow field in a channel. It was computed solely using the implicit incompressible fluid solver (ICFD) of LS-DYNA and its strong two-way FSI coupling possibilities to the implicit structural solver in LS-DYNA.
https://www.dynaexamples.com/icfd/intermediate-examples/intermediate-fsi-2-flaps
https://www.dynaexamples.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
Intermediate: FSI - 2 Flaps
This LS-DYNA simulation shows the fluid-structure interaction (FSI) of two flaps exposed to a flow field in a channel. It was computed solely using the implicit incompressible fluid solver (ICFD) of LS-DYNA and its strong two-way FSI coupling possibilities to the implicit structural solver in LS-DYNA.
*KEYWORD 100M *TITLE *BOUNDARY_SPC_SET *CONTROL_TERMINATION *CONTROL_IMPLICIT_GENERAL *CONTROL_IMPLICIT_DYNAMICS *CONTROL_IMPLICIT_SOLUTION *DATABASE_BINARY_D3PLOT *DEFINE_CURVE *ELEMENT_SHELL *ICFD_BOUNDARY_FREESLIP *ICFD_BOUNDARY_FSI *ICFD_BOUNDARY_PRESCRIBED_VEL *ICFD_BOUNDARY_PRESCRIBED_PRE *ICFD_BOUNDARY_NONSLIP *ICFD_CONTROL_FSI *ICFD_CONTROL_TIME *ICFD_CONTROL_OUTPUT *ICFD_DATABASE_DRAG *ICFD_DATABASE_FLUX *ICFD_MAT *ICFD_PART *ICFD_SECTION *INCLUDE *MAT_ELASTIC *MESH_EMBEDSHELL *MESH_SURFACE_ELEMENT *MESH_SURFACE_NODE *MESH_VOLUME *NODE *PARAMETER *PART *SECTION_SHELL *SET_NODE_LIST_TITLE *END
Fluid velocity streamlines
$-----------------------------------------------------------------------------
$
$ Example provided by Dr. Nils Karajan (DYNAmore)
$
$ 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 as is.
$X Requires LS-DYNA MPP R8.0.0 (or higher) with double precision
$X
$X------------------------------------------------------------------------------
$# UNITS: kg / m / s / N / Pa / Nm (J) / Pa*s
$X------------------------------------------------------------------------------
$X
*KEYWORD 100M
*TITLE
$# title
FSI of two flaps
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ GLOBAL PARAMETER DEFINITION
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PARAMETER
$# name val
r t_end 2.0
$
r dtmech 0.01
r dticfd 0.01
r dt_out 0.01
$
r ramp_vel 0.1
r velocity 5.00
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ FSI DEFINITION
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_BOUNDARY_FSI
$# pid
4
*ICFD_BOUNDARY_FSI
$# pid
5
*ICFD_CONTROL_FSI
$# owc bt dt idc
$===============================================================================
$
$ F L U I D P R O B L E M
$
$===============================================================================
$
$
$---+----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
&t_end &dticfd
*ICFD_CONTROL_OUTPUT
$# msgl outl dtout
0 &dt_out
$
$---
$
*ICFD_DATABASE_FLUX
$# pid
1
*ICFD_DATABASE_FLUX
$# pid
2
*ICFD_DATABASE_DRAG
$# pid cpid
4
*ICFD_DATABASE_DRAG
$# pid cpid
5
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ MESH
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*MESH_VOLUME
$# volid
1
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
1 2 3
*MESH_EMBEDSHELL
$# volid
1
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
4 5
$
$---
$
$*MESH_BL
$# pid nelth
$ 3 3
$*MESH_BL_SYM
$# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8
$ 1 2
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ BOUNDARY CONDITIONS
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# pid dof vad lcid sf vid death birth
1 1 1 2
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# pid dof vad lcid sf vid death birth
1 2 1 2
*ICFD_BOUNDARY_PRESCRIBED_VEL
$# pid dof vad lcid sf vid death birth
1 3 1 1
*ICFD_BOUNDARY_PRESCRIBED_PRE
$# pid lcid sf death birth
2 2
$
$---
$
*DEFINE_CURVE
$# lcid sidr sfa sfo offa offo dattyp
1
$# a1 o1
0.000 0.000
&ramp_vel &velocity
&t_end &velocity
*DEFINE_CURVE
$# lcid sidr sfa sfo offa offo dattyp
2
$# a1 o1
0.000 0.000
1000.000 0.000
$
$---
$
*ICFD_BOUNDARY_FREESLIP
$# pid
3
*ICFD_BOUNDARY_NONSLIP
$# pid
4
*ICFD_BOUNDARY_NONSLIP
$# pid
5
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ SECTION / MATERIAL
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*ICFD_SECTION
$# sid
1
*ICFD_MAT
$# http://www.engineeringtoolbox.com/air-properties-d_156.html
$# http://en.wikipedia.org/wiki/Viscosity#Dynamic_viscosity_.CE.BC
$# dynamic viscosity = (kinematic viscosity * density)
$# mid flg ro vis thd
1 1 1.205 18.21e-6 0.000
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ PART
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*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
$# pid secid mid
5 1 1
*ICFD_PART_VOL
$# pid secid mid
5 1 1
$# spid1 spid2 spid3 spid4 spid5 spid6 spid7 spid8
1 2 3 4 5
$===============================================================================
$
$ S T R U C T U R A L P R O B L E M
$
$===============================================================================
$
$
$---+----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_end
$
$---
$
*CONTROL_IMPLICIT_GENERAL
$# imflag dt0 imform nsbs igs cnstn form zero_v
1&dtmech
*CONTROL_IMPLICIT_DYNAMICS
$# imass gamma beta tdybir tdydth tdybur irate
1
*CONTROL_IMPLICIT_SOLUTION
$# nsolvr ilimit maxref dctol ectol rctol lstol abstol
12 35 50
$# dnorm diverg istif nlprint nlnorm d3itctl cpchk
1 1 1 1
$# arcctl arcdir arclen arcmth arcdmp
$# lsmtd lsdir irad srad awgt sred
5
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ DATABASE CARDS
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*DATABASE_BINARY_D3PLOT
$# dt lcdt beam npltc psetid
&dt_out
$# ioopt
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ PART / SECTION / MATERIAL
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*PART
$# title
flap 1
$# pid secid mid eosid hgid grav adpopt tmid
44 44 44
*PART
$# title
flap 2
$# pid secid mid eosid hgid grav adpopt tmid
55 44 44
*SECTION_SHELL
$# secid elform shrf nip propt qr/irid icomp setyp
44 16 3
$# t1 t2 t3 t4 nloc marea idof edgset
0.0015
*MAT_ELASTIC
$# mid ro e pr da db not used
44 7850.00 2.10e+11 0.30
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
$
$ BOUNDARY CONDITIONS
$
$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
*BOUNDARY_SPC_SET
$# nsid cid dofx dofy dofz dofrx dofry dofrz
44 1 1 1 1 1 1
*BOUNDARY_SPC_SET
$# nsid cid dofx dofy dofz dofrx dofry dofrz
55 1 1 1 1 1 1
*END