中级:FSI-2FlapsLS-DYNA模拟显示双片接触通道流场的流体结构交互单用LS-DYNA隐式不可压缩流体求解器计算,并用它强大的双向FSI联通LS-DYNA隐式结构求解器//www.cannycloudtech.com/icfd/intermediate-examples/intermediate-fsi-2-flaps//www.cannycloudtech.com/@@site-logo/LS-DYNA-Examples-Logo480x80.png
中级:FSI-2Flaps
LS-DYNA模拟显示双片接触通道流场的流体结构交互单用LS-DYNA隐式不可压缩流体求解器计算,并用它强大的双向FSI联通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
流速精简
$----------------------------------------------------------------------------- $ $ Example provided by Dr.18luck快乐彩Nils Karajan(DYNAMore)$E-Mail:fo@dyamore.de$We:http://www.dyamore.de$Flights,2015DYNAMoreGmbH$ $X------------------------------------------------------------------------------ $X $X 1.运行文件原封不动$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 & *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