;---------------------------------------------------------------------------------------------------------
; FLAC3D Data File
;---------------------------------------------------------------------------------------------------------
; Name:         TestStoreShellStrains
; Author:       Roozbeh Geraili Mikola, PhD, PE
; Websites:     www.roozbehgm.com & www.geowizard.org 
; Last Modified: 05-09-2023
; Description:  Simple model to test the 'store_shell_strains()' FISH function
; Notes:        
;---------------------------------------------------------------------------------------------------------

model new
model title ...
    "Calculate Principal Strains in Shell Structural Elements"
; Create shells and assign properties
struct shell create by-quad (0,0,0) (10,0,0) (10,0,1) (0,0,1) size (50,5) ...
                    cross-diagonal element-type=dkt-csth
struct shell property isotropic=(2e11,0.25) thick=1.0
; Boundary conditions
struct node fix velocity-x range position-x=-0.01 0.01
struct node initialize velocity=(1e-6,0,0) range position-x=9.99 10.01
struct node fix velocity-x range position-x=9.99 10.01

; Load FISH functions
program call 'StoreShellStrains.fis'

; Create FISH function to record max. principal strain in the middle of shell
define max_prin_strain_mid
    local sp = struct.near((5.0,0.0,0.5))
    max_prin_strain_mid = struct.extra(sp,1)
end

; Assign history point in the tip and middle of the shell
history interval 100
fish history name='Max_Prin_Strain_Mid' max_prin_strain_mid
structure node history name='Tip_X_Disp' displacement-x position (10.0,0.0,0.5)

; Solve model
model large-strain off
model cycle 50000

model save 'ShellStrain'