基础知识:气流FSI

LS-DYNA模拟扩展基础程序:圆柱流添加结构部件并解决FSI问题结构部件除垂直方向外所有方向固定不变,流体甲板不变,唯一修改是时间步省时计算时间

$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $  FSI deck $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *KEYWORD *INCLUDE *ICFD_BOUNDARY_FSI *ICFD_CONTROL_FSI *END $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $  Solid deck $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *KEYWORD *PARAMETER *TITLE *BOUNDARY_SPC_SET *CONTROL_IMPLICIT_DYNAMICS *CONTROL_IMPLICIT_GENERAL *CONTROL_IMPLICIT_SOLUTION *CONTROL_TERMINATION *ELEMENT_SHELL *INCLUDE *PART *SECTION_SHELL *SET_NODE_LIST *MAT_RIGID *NODE *END $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $ $  Fluid deck $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *KEYWORD *TITLE *INCLUDE *DATABASE_BINARY_D3PLOT *DEFINE_CURVE_TITLE *DEFINE_CURVE_TITLE *ICFD_BOUNDARY_NONSLIP *ICFD_BOUNDARY_PRESCRIBED_VEL *ICFD_BOUNDARY_PRESCRIBED_PRE *ICFD_BOUNDARY_PRESCRIBED_TEMP *ICFD_CONTROL_OUTPUT *ICFD_CONTROL_TIME *ICFD_DATABASE_TIMESTEP *ICFD_DATABASE_TEMP *ICFD_INITIAL *ICFD_MAT *ICFD_PART *ICFD_PART_VOL *ICFD_SECTION *MESH_BL *MESH_BL_SYM *MESH_VOLUME *MESH_NODE *MESH_SURFACE_ELEMENT *END

18luck快乐彩-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------$ $X------------------------------------------------------------------------------ $X $X 1.fsi.k运行文件并选择内存=5mUnitedsLSDYNAMPR8.0$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