$ Pleasea note: Use the DEV version with a timestamp later than $ "Tue Jul 7 19:50:22 2020" *KEYWORD *TITLE $# title Stiffener Coupled Eulerian/Lagrangian *CONTROL_MPP_DECOMPOSITION_DISTRIBUTE_ALE_ELEMENTS *INCLUDE mesh_sale.k *CONTROL_TERMINATION $# endtim endcyc dtmin endeng endmas 300.000000 *CONTROL_TIMESTEP $# dtinit tssfac isdo tslimt dt2ms lctm erode ms1st 0.000 0.600000 $# dt2msf dt2mslc imscl 0.000 0 0 *CONTROL_HOURGLASS $# ihq qh 1 1.0000E-5 *CONTROL_ENERGY $# hgen rwen slnten rylen 2 1 2 2 *DATABASE_ELOUT $# dt binary lcur ioopt 1.000000 *DATABASE_GLSTAT $# dt binary lcur ioopt 1.0000E-8 *DATABASE_NODOUT $# dt binary lcur ioopt dthf binhf 1.0000E-8 *DATABASE_BINARY_D3PLOT $# dt lcdt beam npltc psetid 10.000000 $# ioopt 0 *DATABASE_BINARY_D3THDT $# dt lcdt beam npltc psetid 1.000000 *DATABASE_EXTENT_BINARY $# neiph neips maxint strflg sigflg epsflg rltflg engflg 0 0 0 0 0 0 0 0 $# cmpflg ieverp beamip dcomp shge stssz n3thdt ialemat 0 0 0 0 0 0 0 1 $# nintsld pkp_sen sclp unused msscl therm 0 0 1.000000 *PART $# title material type # 3 (Kinematic/Isotropic Elastic-Plastic) $# pid secid mid eosid hgid grav adpopt tmid 1001 1 1 *PART $# title material type # 3 (Kinematic/Isotropic Elastic-Plastic) $# pid secid mid eosid hgid grav adpopt tmid 1 1 1 *SECTION_SHELL $# secid elform shrf nip propt qr/irid icomp setyp 1 2 0.000 3 $# t1 t2 t3 t4 nloc marea idof edgset 0.200000 0.200000 0.200000 0.200000 *MAT_PLASTIC_KINEMATIC $# mid ro e pr sigy etan beta 1 7.830000 2.070000 0.300000 0.008000 $# src srp fs vp 0.000 0.000 0.000 0.000 *DEFINE_CURVE $ lcid sidr sfa sfo offa offo dattyp 1 0 0.0 0.0 0.0 0.0 $ a1 o1 0.0 0.0 0.1 0.0010 *SET_PART_LIST $# sid da1 da2 da3 da4 1 $# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8 1 *SET_PART_LIST $# sid da1 da2 da3 da4 2 $# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8 9 *SET_PART_LIST $# sid da1 da2 da3 da4 1011 $# pid1 pid2 pid3 pid4 pid5 pid6 pid7 pid8 10 11 *SET_MULTI-MATERIAL_GROUP_LIST 1 1 2 *SET_MULTI-MATERIAL_GROUP_LIST 2 3 *ALE_STRUCTURED_FSI $# slave master sstyp mstyp mcoup 1 2 0 0 -1 $# start end pfac fric norm 0.1 1 *END *CONSTRAINED_LAGRANGE_IN_SOLID $# slave master sstyp mstyp nquad ctype direc mcoup 1 2 0 0 2 4 1 -1 $# start end pfac fric frcmin norm normtyp damp 0.1 1 *END $ $ ASF only does "compression only". This means we only push ALE fluild out if the $ fluid is penetrating into the structure. Fluids could move away freely from the $ structure. $ $ A separate ASF card is needed if we do not allow the moving away motion. In the $ extra card, we do not allow the fluid on the opposite side to penetration in. $ As shown in the following card. We flip the Lagrange structure normal (from 0 to $ 1), and we couple to air (-2) instead of HE and water (-1) $ $ This two cards together, is equivelant to DIREC=1 option in CONSTRAINED_LAG_IN $ _SOLID. $ $ DIREC=1 option is not supported by ASF. This is because ASF has automated leakage $ detect and control. So leakage control is not an option. It is always on. $ But CLIS's leakage control, when activated, only act on the "compression" motionr. $ Moving away motion, no leakage control, so this means potentially, the fluid $ on the other side could still leak in, without being detected and controlled. $ $ To make things simpler, the two cards approach is taken by *ALE_STRUCTURED_FSI $ This way, fluid from both side are properly constrained on their own sides. Not $ allowed to flow pass the structure. Leakage on both side are treated and controled. $ $ Below the first card is to couple HE and water to the structure (normal pointing $ dowwards). The second card is to couple air to the structure (normal pointing $ upwards) *ALE_STRUCTURED_FSI $# slave master sstyp mstyp mcoup 1 2 0 0 -1 $# start end pfac fric frcmin norm normtyp damp 0.1 1 *ALE_STRUCTURED_FSI $# slave master sstyp mstyp mcoup 1 2 0 0 -2 $# start end pfac fric frcmin norm normtyp damp 0.10 *END