Multi-Body Contact Dynamics in Zebulon

This is an amusing calculation showing multi-body contact with the implicit dynamic solution. The letters are made of a fairly massive soft elastic material and were meshed using a script which takes vector boundary input from external files (Z-program interface with Zmaster which will be available in the next release). The letters were a vectorized translation of Helvetica font.

****calcul dynamic
  ***mesh updated_lagrangian_plane_strain
    ***resolution newton
      **sequence
        *time      1.1 2.
        *increment 10  200
        *iteration 6
        *algorithm p1p2p3
   **automatic_time automatic_by_sequence ETO 0.1
    *divergence 2. 10
    *security  1.25
  % *max_dtime 0.1
  ***bc
    **gravity
      ZeBuLoN d2 -9.8 tab1
    **impose_nodal_dof
      fix U1 0.0
      fix U2 0.0
  ***table
   **name tab1
    *time  0. 0.01 1.e9
    *value 0. 1.   1.
  ***contact soft
    **conv 1.e-8 100
    **soft_param 1.e3 6.0 1.e-8
    %
    % ZeBuLoN to box contact
    %
    **zone
     *impactor ZeBuLoN
     *target   box
     *friction 0.5
     *gap      1.e-6
     *warning_distance 0.4
    %
    %
    %
    **zone
     *impactor eBuLoN
     *target   Z
     *friction 0.5
     *gap      1.e-6
     *warning_distance 0.1
    **zone
     *impactor BuLoN
     *target   e
     *friction 0.5
     *gap      1.e-6
     *warning_distance 0.1
    **zone
     *impactor uLoN
     *target   B
     *friction 0.5
     *gap      1.e-6
     *warning_distance 0.1

    **zone
     *impactor LoN
     *target   u
     *friction 0.5
     *gap      1.e-6
     *warning_distance 0.1

    **zone
     *impactor oN
     *target   L
     *friction 0.5
     *gap      1.e-6
     *warning_distance 0.1

    **zone
     *impactor N
     *target   o
     *friction 0.5
     *gap      1.e-6
     *warning_distance 0.1

  ***material
   **elset ZeBuLoN
    *file calc3.inp 1
   **elset A15
    *file calc3.inp 2
****return

***behavior linear_elastic lagrange_polar
 **elasticity
     young 2000.0
     poisson 0.3
 ***coefficient
   masvol 1.e-2
***return

***behavior linear_elastic lagrange_polar
 **elasticity
     young 200000.0
     poisson 0.3
 ***coefficient
   masvol 1.
***return