# USE UMATS8R5.FOR
# When using UMAT, must enter *TRANSVERSE SHEAR STIFFNESS

# open Ex_3.8.cae and set the WD to have all files local
Menu:   File, Open, c:\Simulia\User\Ex_3.8\Ex_3.8.cae
        Save As, c:\Simulia\User\Ex_3.14\Ex_3.14.cae
        Set Work Directory, c:\Simulia\User\Ex_3.14

# only edit what needs to be changed
Module: Property
# unidirectional material from Ex. 3.8 used to calculate H44, H45, H55
Menu:   Material, Name [unidirectional], Mechanical, Elasticity, Elastic, Lamina,
        [133860     7706     0.301     4306     4306     2760], OK
        # or right click on data field, read from File: [uni_lam_prop.txt]
        # add Fail Stress values
        Suboptions, Fail stress,
        [1830    -1096    57    -228    71    -1]
        # or right click on data field, read from File: [uni_fail_stress.txt]
        OK, OK
            
Module: Job
Menu:   Job, Manager, Edit, 
        (tab) General, [checkmark] Print model definition data,
        Submit, OK
        Results    # max. U3 = 8.704 mm
        # open Job-1.dat file to read 
        # Kts11=35833.0=K11=H55, Kts12=0.0=K12=H45, Kts22=23000.0=K22=H44

# run with UMAT
# create user material to be used with UMAT
Module: Property
Menu:   Material, Name [ud-user], 
        General, User material, Mechanical, Mechanical constants:
        [133860., 7706., 0.301, 0.396, 
        4306., 1830., 1096., 57.,  228.,   71.,   -1.]
        # or right click on data field, read from file: [props.txt], OK
        Depvar, 2, OK, Dis  # see umats8r5.for for interpretation

Menu:   Section, Edit, Section-1, 
        (tab) Basic, 
        Section integration, During analysis, 
        Thickness, Shell thickness, Value, [10]
        Material, [ud-user]
        Thickness integration rule, Simpson
        Thickness integration points, [3]
        (tab) Advanced, 
        Transverse Shear Stiffness, # checkmark Specify Values
        K11: [35833], K12: [0], K22: [23000]

Module: Step
Menu:   Output, Field output requests, Edit, F-Output-1,
        Output variables: Edit variables: [S,E,U,RF,SDV], OK
            
Module: Job
Menu:   Job, Manager,
        (tab) Advanced
        User subroutine file: [umats8r5.for], OK 
        Submit
        Results
            
        # The following must appear in the .inp file if using UMAT
        # ** Section: Section-1
        # *Shell Section, elset=_PickedSet2, material=uni-user, orientation=Ori-1
        # 10., 3
        # *Transverse Shear
        # 35883., 23000., 0.
        # and
        # *Material, name=uni-user
        # *Depvar
        # 2,
        # *User Material, constants=11
        # 133860., 7706., 0.301, 0.396, 4306., 1830., 1096.,   57.
        # 228.,   71.,   -1.
        
Module: Visualization
Menu:   Options, Contour, Quilt
Menu:   Plot, Contours on Deformed Shape, SDV1
Menu:   Result, Selection method: Plies, Section points: Bottommost
        # the bottom side reaches F_I~1.767E-2
Menu:   Result, Selection method: Plies, Section points: Topmost
        # the top side reaches F_I~9.444E-3
        # or use script TsaiWu.py to get almost the same results