fino_reltol = $1 fino_penalty_weight = $2 SHELL "gmsh -v 0 -3 -clmax 0.03 cylinder.geo -o cylinder-fine.msh" MESH FILE_PATH cylinder-fine.msh E = 1 nu = 0.3 PHYSICAL_GROUP bottom BC tangential radial PHYSICAL_GROUP top BC P=1 FINO_STEP # compute min/max stresses MESH_FIND_MINMAX FUNCTION sigmax MIN sigmax_min MAX sigmax_max MESH_FIND_MINMAX FUNCTION sigmay MIN sigmay_min MAX sigmay_max MESH_FIND_MINMAX FUNCTION sigmaz MIN sigmaz_min MAX sigmaz_max MESH_FIND_MINMAX FUNCTION tauxy MIN tauxy_min MAX tauxy_max MESH_FIND_MINMAX FUNCTION tauyz MIN tauyz_min MAX tauyz_max MESH_FIND_MINMAX FUNCTION tauzx MIN tauzx_min MAX tauzx_max PRINT "stress\tmin\t\tmax" PRINT SEP "\t" "% +5e" "sigmax" sigmax_min sigmax_max PRINT SEP "\t" "% .5e" "sigmay" sigmay_min sigmay_max PRINT SEP "\t" "% .5e" "sigmaz" sigmaz_min sigmaz_max PRINT SEP "\t" "% .5e" "tauxy" tauxy_min tauxy_max PRINT SEP "\t" "% .5e" "tauyz" tauyz_min tauyz_max PRINT SEP "\t" "% .5e" "tauzx" tauzx_min tauzx_max