高级 : 流血

举个例子血液流出动脉Carrreau非牛顿模型已被选为血液模型windecsel类型边界条件已加诸不同容器插口,象征外部验血负载LS-DYNA版本R10.1或后使用

*KEYWORD *TITLE *DATABASE_BINARY_D3PLOT *DEFINE_CURVE *ICFD_BOUNDARY_NONSLIP *ICFD_BOUNDARY_PRESCRIBED_PRE *ICFD_BOUNDARY_WINDKESSEL *ICFD_CONTROL_MESH *ICFD_CONTROL_TIME *ICFD_DATABASE_FLUX *ICFD_DATABASE_TIMESTEP *ICFD_MAT *ICFD_MODEL_NONNEWT *ICFD_PART *ICFD_PART_VOL *ICFD_SECTION *INCLUDE *MESH_SURFACE_ELEMENT *MESH_SURFACE_NODE *MESH_VOLUME *PARAMETERS *END

18luck快乐彩LSSTC和Hossein Mohammadi大学提供$ $X------------------------------------------------------------------------------ $X $X 1.运行文件原封不动$X    Requires LS-DYNA MPP Dev 117000 (or higher) with double precision  $X $X------------------------------------------------------------------------------ $# UNITS: (g/cm/s) $X------------------------------------------------------------------------------ $X *KEYWORD *TITLE ICFD Blood flow *INCLUDE mesh.k $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $                                                                              $ $                             PARAMETERS                                       $ $                                                                              $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *PARAMETER R    t_end       0.5 R dt_fluid      5e-5 R  dt_plot     0.005 Rp_init       126656 $ $--- Fluid $ $ Carreau model for blood flow Rrho_fluid     1.059 R mu_fluid     0.035 R mu_o          0.56 R n           0.3568 R lambda       3.313 $ Windkessel parameters R ra4          45000 R ca4         8.e-08 R ram4        75000.Rv43.e-6Rv470000R5100000R51.0E-07Rl5160000Rv54.e-6Rv540000Rra625000RCA66.0E-07RRcv61.8e-5Rrv69500$---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $                                                                              $ $                           ICFD CONTROL CARDS                                 $ $                                                                              $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *ICFD_CONTROL_MESH       1.05 *ICFD_CONTROL_TIME $#     ttm        dt       cfl     &t_end &dt_fluid $---+----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                1         1&rho_fluid &mu_fluid               2 *ICFD_MODEL_NONNEWT          2        2      &mu_o       &n  &mu_fluid   &lambda *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 $#     pid     secid       mid          6         1         1 *ICFD_PART_VOL $#     pid     secid       mid          7         1         1 $#   spid1     spid2     spid3     spid4     spid5     spid6     spid7     spid8          1         2         3         4         5         6  *ICFD_SECTION $#     sid          1 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $                                                                              $ $                    ICFD BOUNDARY/INITIAL/LOAD CONDITIONS                     $ $                                                                              $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *ICFD_BOUNDARY_PRESCRIBED_PRE          1         2 *ICFD_BOUNDARY_PRESCRIBED_PRE          2         1 *ICFD_BOUNDARY_NONSLIP $#     pid          3 *ICFD_BOUNDARY_WINDKESSEL $#     pid     wtype        r1        c1        r2          4         3      &ra4      &ca4     &ram4 $#  p2lcid        c2        r3                 &cv4      &rv4 *ICFD_BOUNDARY_WINDKESSEL $#     pid     wtype        r1        c1        r2          5         3      &ra5      &ca5     &ram5 $#  p2lcid        c2        r3                 &cv5      &rv5 *ICFD_BOUNDARY_WINDKESSEL $#     pid     wtype        r1        c1        r2          6         3      &ra6      &ca6     &ram6 $#  p2lcid        c2        r3                 &cv6      &rv6 *ICFD_INITIAL          0         0         0         0         0   &p_init *DEFINE_CURVE 1,,,1,-0.0243, 0,12239 0.0243,126656.2705 0.0328,129456.0407 0.0412,133189.0676 0.0496,137322.0617 0.058,141321.7334 0.0665,146654.629 0.0749,150654.3007 0.0833,153320.7485 0.0917,155987.1963 0.1,157320.4202 0.109,158653.6441 0.117,159986.868 0.125,159986.868 0.134,159986.868 0.142,159986.868 0.151,159986.868 0.159,158653.6441 0.168,155987.1963 0.176,153320.7485 0.184,150654.3007 0.193,147987.8529 0.201,145321.4051 0.21,143988.1812 0.218,141321.7334 0.227,138655.2856 0.235,135988.8378 0.243,134655.6139 0.252,134655.6139 0.26,135988.8378 0.269,138655.2856 0.277,139988.5095 0.286,142654.9573 0.294,143988.1812 0.302,143988.1812 0.311,142654.9573 0.319,141321.7334 0.328,138655.2856 0.336,137322.0617 0.344,135988.8378 0.353,134655.6139 0.361,133322.39 0.37,133055.7452 0.378,132922.4228 0.387,132655.7781 0.395,131589.1989 0.403,130122.6526 0.412,129189.3959 0.42,128656.1064 0.429,128389.4616 0.437,127989.4944 0.446,127589.5272 0.454,127189.5601 0.462,126789.5929 0.471,125989.6586 0.479,125189.7242 0.488,124389.7899 0.496,123589.8555 0.505,123189.8884 *DEFINE_CURVE 2 0,0 1000,0 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $                                                                              $ $                            ICFD MESH KEYWORDS                                $ $                                                                              $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *MESH_BL          3         1 *MESH_BL_SYM          1 *MESH_BL_SYM          2 *MESH_BL_SYM          4 *MESH_BL_SYM          5 *MESH_BL_SYM          6 *MESH_VOLUME $#   volid         10 $#    pid1      pid2      pid3      pid4      pid5      pid6      pid7      pid8          1         2         3         4         5         6 $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 $                                                                              $ $                             DATABASE (OUTPUT)                                $ $                                                                              $ $---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8 *ICFD_DATABASE_FLUX          1 *ICFD_DATABASE_FLUX          2 *ICFD_DATABASE_FLUX          4 *ICFD_DATABASE_FLUX          5 *ICFD_DATABASE_FLUX          6 *ICFD_DATABASE_TIMESTEP          1 *DATABASE_BINARY_D3PLOT $#      dt      lcdt      beam     npltc    psetid &dt_plot *END