Users:General FEM Analysis/syntax

!###################################################################
!####                          PC-BLOCK                         ####
!###################################################################
!************************
!*       INCLUDES       *
!************************
!===================================================================
PC-INCLUDE
  salve.txt
!===================================================================
!***************************
!*       MASTER JOBS       *
!***************************
!===================================================================
PC-PROBLEM
  MASTERJOB = PC-ANALYSIS 1
 ! MASTERJOB = OPT-CTR 1
!===================================================================
!===================================================================
PC-ANALYSIS 1: STA_GEO_LIN
SOLVER = PC-SOLVER 2
OUTPUT = PC-OUT 1
COMPCASE = LD-COM 1
DOMAIN = EL-DOMAIN 1
!===================================================================
!===================================================================
PC-ANALYSIS 2: STA_GEO_NONLIN
  PATHCONTROL = FORCE ! or DISPLACEMENT or ARCLENGTH
  SOLVER = PC-SOLVER 3
  OUTPUT = PC-OUT 1
  COMPCASE = LD-COM 1
  DOMAIN = EL-DOMAIN 1
  NUM_STEP = 3
  MAX_ITER_EQUILIBRIUM =30
  EQUILIBRIUM_ACCURACY = 1e-12
  CURVE=LD-CURVE 1
  TRACED_NODE=3246 
  TRACED_NODAL_DOF=DISP_Y
!
  STEP_LENGTH_CONTROL = FIXED
  STEP_LENGTH_CONTROL_REALS = 1.0
!
 ! SIMULTANEOUS_EIGENVALUE_ANALYSIS = 1
 ! EIGEN_SOLVER = PC-SOLVER 2
!
  REALLOC_MTX = 0  !  1  0
!***********
LD-CURVE 1 TYPE=DISCRETE
 TIME=0.000   VAL=0.000 
 TIME=3.000   VAL=1.000
!===================================================================
PC-ANALYSIS 3: EIGENFREQUENCY
  EIGEN_SOLVER = PC-SOLVER 5
  OUTPUT = PC-OUT 1
  COMPCASE = LD-COM 1
  DOMAIN = EL-DOMAIN 1
  NUM_ROOT = 4
  MASS_LUMPED = 0
!===================================================================
PC-ANALYSIS 4: LINEAR_BUCKLING
  EIGEN_SOLVER = PC-SOLVER 1
  LINEAR_SOLVER = PC-SOLVER 2
  OUTPUT = PC-OUT 1
  COMPCASE = LD-COM 1
  DOMAIN = EL-DOMAIN 1
  NUM_ROOT = 3
!==================================================================
!******************************
!*        GID                 *
!******************************
!===================================================================
PC-OUT 1 : GID
   GEOM=1
   DISP=1
   STRESS=1 
   PREC=10
   DOMAIN_DECOMP = 1
!===================================================================
!       GID FOR OPTIMIZATION
PC-OUT 2 : GID
   GEOM=1
   DISP=1
   STRESS=0
   PREC=8
   FPN=1
   PURE_GRAD = 1
   SMOOTH_GRAD = 1
   DESIGN_UPDATE = 1
   DESIGN_VARIABLE = 1
   FILTER_FUNCTION = 1
   SPREAD_FILES = 1
   DOMAIN_DECOMP = 1
   NODAL_DIRECTOR_REF = 1
   TEMPERATURE = 1
   REGULARIZATION_OBJECTS = 1
!===================================================================
!       GID FOR OPTIMIZATION'S ANALYSIS
PC-OUT 3 : GID
   GEOM=0
!===================================================================
!*******************************
!*      LINEAR SOLVERS         *
!*******************************
!===================================================================
PC-SOLVER 1: SANDIA_LDL
 BANDWITH = CUTHILL_MCKEE
!===================================================================
PC-SOLVER 2: TRILINOS_AMESOS
 TYPE=UMFPACK    !    SUPER_LU_DIST  KLU  UMFPACK
 BANDWITH = CUTHILL_MCKEE
!===================================================================
PC-SOLVER 3 : TRILINOS_AZTECOO
 SCALING = DIAG   !  DIAG  NONE
 PRECOND=DOMAIN_DECOMP  !NONE  JACOBI  NEUMANN  LS  IFPACK_SYMM_GS  DOMAIN_DECOMP
 SUBDOMAIN_SOLVE=ICC   !ICC LU
 REL_TOLERANCE=1E-7  MAX_NUM_ITER=2000000  
 REORDER=1  TYPE_OVERLAP=SYMMETRIC OVERLAP=1
 GRAPH_FILL=0  KEEP_INFO=0  REUSE=0 OUTPUT=0 DIAGNOSTICS=0
 BANDWITH = CUTHILL_MCKEE
!===================================================================
PC-SOLVER 4: TRILINOS_CG
 PRECOND = IFPACK_SYMM_GS ! NONE NEUMANN IFPACK_IC IFPACK_SYMM_GS
 SCALING = DIAG   ! NONE DIAG
 REL_TOLERANCE=1E-8
 MAX_NUM_ITER=10000
 OUTPUT = 100
 BANDWITH = CUTHILL_MCKEE
!===================================================================
PC-SOLVER 5: CROUT_SKYLINE
 BANDWITH = CUTHILL_MCKEE
!===================================================================
PC-SOLVER 6: CHOLMOD
 BANDWITH = CUTHILL_MCKEE
!===================================================================
PC-SOLVER 1: EIGEN_SUBSPACE_IT
 SUBSOLVE = CROUT_SKYLINE ! or CHOLMOD
 SUB_SOLVER = PC-SOLVER 1 
 SUBSOLVE_EIG = LAPACK ! or JACOBI or QR
 BANDWITH = CUTHILL_MCKEE
 REL_TOLERANCE = 1E-12
 MAX_NUM_ITER = 5000
!=====================================================================
!*******************************
!*       EIGEN SOLVERS         *
!*******************************
!===================================================================
PC-SOLVER 6: TRILINOS_EIGEN_SUBSPACE_IT
 BANDWITH = CUTHILL_MCKEE
 REL_TOLERANCE = 1E-7
 MAX_NUM_ITER = 1000
 SUB_SOLVER = PC-SOLVER 2
 OUTPUT = WARNINGS
!===================================================================
!===================================================================
PC-SOLVER 7: TRILINOS_ANASAZI
 BANDWITH = CUTHILL_MCKEE
 REL_TOLERANCE = 1E-7
 TYPE=BLOCK_KRYLOV_SCHUR
 MAX_NUM_ITER = 1000
 SUB_SOLVER = PC-SOLVER 2
 OUTPUT = WARNINGS
!===================================================================
PC-SOLVER 8: EIGEN_SUBSPACE_IT
 SUBSOLVE = CHOLMOD   ! CHOLMOD  or  CROUT_SKYLINE
 SUB_SOLVER = PC-SOLVER 1
 SUBSOLVE_EIG = LAPACK
 BANDWITH = CUTHILL_MCKEE
 REL_TOLERANCE = 1E-9
 MAX_NUM_ITER = 30
!===================================================================
!
!###################################################################
!####                          ND-BLOCK                         ####
!###################################################################
!===================================================================
ND-COOR
  NODE  1  X  0  Y  0  Z  0
  NODE  2  X  0.125  Y  0  Z  0
!===================================================================
ND-SET 1 NAME = set1
  NODE  1  2  3
!===================================================================
!
!
!###################################################################
!####                          EL-BLOCK                         ####
!###################################################################
!***************************
!*       ELEMENT PROPS     *
!***************************
!=================================================================== 
EL-PROP 1: SOLIDHEXA1
INT_TYPE_HEXA = FULL  !  REDUCED  or FULL
EAS = 0
MAT = EL-MAT 1
!===================================================================
!=================================================================== 
EL-PROP 2: SOLIDTETRA1
INT_TYPE_TETRA = FULL
MAT = EL-MAT 1
!===================================================================
!=================================================================== 
EL-PROP 3: SOLIDPENTA1
INT_TYPE_PENTA = FULL
EAS = 0
MAT = EL-MAT 1
!===================================================================
!=================================================================== 
EL-PROP 4 : SHELL8
MAT=EL-MAT 1
THICKNESS=0.1 
INT_TYPE_SHELL8 = FULL
SDC=1.0
EAS =0,0,0,0,0 
ANS =NONE
FORCE = RST_ORTHO
!===================================================================
!=================================================================== 
EL-PROP 5 : QUAD1
MAT=EL-MAT 1
THICKNESS=0.1 
INT_TYPE_QUAD = FULL
EAS=0
MODE=PSTRESS
!=================================================================== 
EL-PROP 6 : TRUSS1
MAT= EL-MAT 1        AREA=1.0
PRESTRESS       SIG11=0.0
LAGRANGE=TOTAL
!===================================================================
EL-PROP 7 : SHELL9_COMPOSITE
INT_TYPE_SHELL9 = FULL
SDC= 1.0
EAS = 4,4,4,4,4
!EAS = 0,0,0,0,0
ANS = NONE
FORCE = XYZ        !  RST_ORTHO  XYZ
!---------------------
ROTATORIC_DOFS = DIR_DIFF  ! DIR_DIFF  GLOB_6
STRESS_OUTPUT = SIG_EPS ! PRIN_VM  SIG_EPS
STRESS_EVAL   = SURF        ! GP SURF
!---------------------
MATERIAL_ANGLE = 0          ! in degrees
MATERIAL_DIRECTION = 1, 0, 0
PROJECTION = X-AXIS  !  Y-AXIS   Z-AXIS   DIRECTOR
!---------------------
  MAT=EL-MAT 1  THICK=0.25  THETA=0        !  PLY 1        bottom
  MAT=EL-MAT 1  THICK=0.25  THETA=90       !  PLY 2          ||
  MAT=EL-MAT 1  THICK=0.25  THETA=90       !  PLY 3          \/
  MAT=EL-MAT 1  THICK=0.25  THETA=0        !  PLY 4         top
!===================================================================
EL-PROP 8 : SHELL9
INT_TYPE_SHELL9 = FULL
MAT=EL-MAT 1
THICKNESS=0.1 
SDC= 1.0
EAS = 4,4,4,4,4
!EAS = 0,0,0,0,0
ANS = NONE
FORCE = RST_ORTHO
!---------------------
ROTATORIC_DOFS = DIR_DIFF  ! DIR_DIFF  GLOB_6
MATERIAL_ANGLE = 0    ! in degrees
STRESS_OUTPUT = SIG_EPS ! PRIN_VM  SIG_EPS
STRESS_EVAL   = SURF        ! GP SURF
!===================================================================
!***********************
!*      MATERIAL       *
!***********************
!===================================================================
EL-MAT 1 : LIN_ELAST_ISOTROPIC
EMOD=20000 ALPHAT=0.0 DENS=0.0 NUE=0.2
!===================================================================
!===================================================================
EL-MAT 2 : ST_VENANT_KIRCHHOFF_ISOTROPIC
EMOD=20000  DENS=0.0 NUE=0.2
!===================================================================
EL-MAT 3 : LIN_ELAST_ORTHOTROPIC_3D
!
 EMOD1=4.32E8  EMOD2=4.32E8
 G12=2.16E8     G13=2.16E8     G23=2.16E8
!
 NUE12=0.0
 ALPHAT1=0.0  ALPHAT2=0.0  ALPHAT3=0.0
 DENS=4
 XSI=1.2
!
!===================================================================
!***********************
!*     TOPOLOGY        *
!***********************
!===================================================================
EL-DOMAIN 1
 ELEMENTS = EL-TOP 1
!===================================================================
!===================================================================
EL-PART 1 NAME=p1
!===================================================================
EL-TOP 1
  NEL 1  1  1  1  2  3  4  5  6  7  8
  NEL 2  1  1  5  6  7  8  9  10  11  12
!===================================================================
!
!
!###################################################################
!####                          LD-BLOCK                         ####
!###################################################################
!===================================================================
LD-COM 1 
 TYPE=LD-NODE 1 FAC= 1.0
 TYPE=LD-ELEM 1 FAC= 1.0
 TYPE=BC-DIRICHLET 1
!=================================================================== 
!===================================================================
BC-DIRICHLET 1: SPC-ZERO
  NODE 1  DISP_X, DISP_Y, DISP_Z
  NODE 4  DISP_X, DISP_Y, DISP_Z
!===================================================================
LD-NODE 1 TYPE=FORCE
  NODE 481  D1=0.0  D2=0.0  D3=1.0   VAL=-1
!===================================================================
LD-ELEM 1 PART=1
 TYPE=DEAD      D1=0.0   D2=0.0   D3=1.0  VAL=100.0
!===================================================================
LD-ELEM 2
 TYPE = TEMPERATURE
 NODES  = 10,11,12  LAYER = 1  VAL = +10
 ND-SET = 2         LAYER = 1  VAL = -10
!===================================================================
!
!===================================================================
!###################################################################
!####                         OPT-BLOCK                         ####
!###################################################################
!===================================================================
OPT-BOUND 1 : ABSOLUTE_BOX
  CEN_X=4500.0 CEN_Y=-480.0 CEN_Z=964.2 DIM_X=500.0 DIM_Y=960.0 DIM_Z=10   
!===================================================================
OPT-BOUND 1 : ABSOLUTE_PLANE
  P_X = 0        P_Y = 0        P_Z = 10
  NORMAL_X = 0     NORMAL_Y = 0     NORMAL_Z = 1
  HALFSPACE =  POSITIVE  !  POSITIVE or NEGATIVE
!===================================================================
OPT-BOUND 1 : SCALAR_BOUND
 LOWER = 0.05       UPPER = 0.15
!===================================================================
OPT-BOUND 1 : CURVATURE
 MAX_MEAN_CURVATURE = 0.1
!===================================================================
!===================================================================
OPT-RESPONSE_FCT 1 : DISPLACEMENT_LIN
 WEIGHT=1.0 ANALYSIS=PC-ANALYSIS 1 ETA=1e-06
 GRAD=ADJOINT SA=SEMI_ANALYTIC FDA=FOREWARD
 DESVAR=OPT-VAR 2
 NODE = 552
 DOF = DISP_Y
 GOAL = MAX
!===================================================================
OPT-RESPONSE_FCT 2 : DISPLACEMENT_LIN
! -- basic stuff
 WEIGHT=1.0 ANALYSIS=PC-ANALYSIS 1 ETA=1e-06
 GRAD=ADJOINT SA=SEMI_ANALYTIC FDA=FOREWARD
 DESVAR=OPT-VAR 2
! -- response function dependant parameters
 LOAD_CONSTANT = 1
 NODE = 552
 DOF = DISP_Z
 GOAL = MIN
! -- constraint parameters
 REL_LIMIT = 1.00
 ! ABS_LIMIT = 1-e3
 REL_TOLERANCE = 0.05
 CONSTRAINT_TYPE = INEQUALITY_LT   ! 'EQUALITY' or  'INEQUALITY_LT'
 LAMBDA_ABS_MAX = 40
!===================================================================
OPT-RESPONSE_FCT 3 : MASS
 WEIGHT=1.0 DOMAIN=EL-DOMAIN 1 ETA=1e-06
 GRAD=ADJOINT SA=SEMI_ANALYTIC FDA=FOREWARD
 DESVAR=OPT-VAR 1
! -- constraint parameters
 REL_LIMIT = 1.00
 REL_TOLERANCE = 0.01
 CONSTRAINT_TYPE = EQUALITY   ! 'EQUALITY' or  'INEQUALITY_LT'
 LAMBDA_ABS_MAX = 40
!===================================================================
!            BUCKLING CONSTRAINT
OPT-RESPONSE_FCT 4 : EIGENVALUE_NLN_BUCKLING
 WEIGHT=1.0 ANALYSIS=PC-ANALYSIS 3 ETA=1e-06
 GRAD=DIRECT 
 SA=SEMI_ANALYTIC ! EXACT_SEMI_ANALYTIC  GLOBAL_FD
 FDA=FOREWARD
 DESVAR=OPT-VAR 1
!----------
 GOAL = MAX
 NUM_ROOTS = 1
 RHO = 1.0
! -- constraint parameters
 ABS_LIMIT = 79.5812   !  - sqrt(min_eigenvalue * 0.5 / initial_eigenvalue) + 80
 REL_TOLERANCE = 0.01
 CONSTRAINT_TYPE = INEQUALITY_LT   ! 'EQUALITY' or  'INEQUALITY_LT'
 LAMBDA_ABS_MAX = 40
!===================================================================
OPT-RESPONSE_FCT 2 :  CURVATURE_PATCH_KS   ! -- LIMITATION OF SURFACE CURVATURE --
 WEIGHT=1.0  ETA=1e-06
 GRAD=ADJOINT SA=SEMI_ANALYTIC FDA=FOREWARD
 DESVAR=OPT-VAR 1
! -- specific parameters 
 NODES=ND-SET 1
 PATCH_SIZE = 3
 RHO = 25
! -- constraint parameters
 ABS_LIMIT = 0.04         ! MAX MEAN CURVATURE
 REL_TOLERANCE = 0.05
 CONSTRAINT_TYPE = INEQUALITY_LT   ! 'EQUALITY' or  'INEQUALITY_LT'
 LAMBDA_ABS_MAX = 6
!===================================================================
OPT-RESPONSE_FCT 1 : EIGENFREQUENCY_KS
 WEIGHT=1.0 ANALYSIS=PC-ANALYSIS 1 ETA=1e-06
 GRAD=ADJOINT SA=SEMI_ANALYTIC FDA=FOREWARD
 RHO=100
 DESVAR=OPT-VAR 1 
!===================================================================
!===================================================================
OPT-FILTER 1 : TOPO   !   GLOBAL   TOPO
 FUNCTIONS = OPT-LINKFUNC 1
!===================================================================
OPT-LINKFUNC 1 : CIRC 
 SHAPE=LIN   ! CONST  LIN QUAD  CUB
 RADIUS          = 1.1
 BOUNDARY_RADIUS = 40
 MINIMUM_RADIUS = 1.1
 PERFORM_DIRECTION_WEIGHTING = 1
 CONCIDER_BLOCK_INTERACTION = 1
!===================================================================
!===================================================================
OPT-VAR 1
 TYPE=ND-SET 2 SUBTYPE=VEC  D1=0.0  D2=0.0 D3=1.0      BOUND=OPT-BOUND 1       STEPLENGTH_FACTOR = 0.01 
 TYPE=ND-SET 1 SUBTYPE=VEC  D1=0.0  D2=0.0 D3=-1.0         BOUNDVAR = 1
!
 TYPE=EL-SET 1 SUBTYPE=THICKNESS 
 TYPE=EL-SET 1 SUBTYPE=THICKNESS LAYER_ID = 1
!
 TYPE=EL-SET 1 SUBTYPE=MATERIAL_ANGLES
 TYPE=EL-SET 1 SUBTYPE=MATERIAL_ANGLE        LAYER_IDS = 1,2,3
!===================================================================
!===================================================================
OPT-CTR 1 : ALM_NAND      !ALM_NAND   CG_NAND    STEEP_DESC_NAND
 FILTER = OPT-FILTER 1
!
 LINE_SEARCH = OPT-LINE_SEARCH 1
!
 DESIGN_SPACE_BOUNDS = ND-SET 1
!
 OBJ=OPT-RESPONSE_FCT 1
 CON=OPT-RESPONSE_FCT 3,4
!
 OUTPUT=PC-OUT 1
 DOMAIN = EL-DOMAIN 1
 REGULARIZATION = EL-REGULARIZATION 1
!
 CONVERGENCE_CONTROL = OPT-CONVERGENCE 1
!
 RESTART_DATA_FREQ = -1  ! 0: no restart  -1: write resart data at the end of opt  n>1: write data every n steps
!********** ALM-PARAMETERS ******
 PENALTY = 1.0
 PENALTY_MAX = 10.0
 GAMMA = 1.1
!===================================================================
OPT-LINE_SEARCH 1 : FIXED
!OPT-LINE_SEARCH 1 : 3-POINT
!OPT-LINE_SEARCH 1 : 3-POINT_ADAPTIVE
!OPT-LINE_SEARCH 1 : ARMIJO
!OPT-LINE_SEARCH 1 : ARMIJO_3-POINT_ADAPTIVE
!
MAX_STEP_LENGTH = 0.05
MIN_STEP_LENGTH = 0.01
ARMIJO_PARAMETER = 0.2
!
MAX_NUM_RESTART   = 3
!===================================================================
OPT-CONVERGENCE 1
  CHECK_CONVERGENCE = 1
  OBJECTIVE_TOL     = 1e-2
  MAX_STEP          = 40
  CONVERGENCE_STEPS = 3
!===================================================================
!###################################################################
!####                         REG-BLOCK                         ####
!###################################################################
!===================================================================
EL-REGULARIZATION 1 : LAPLACE
 DOMAIN = EL-DOMAIN 1
!
! AUTOMATIC_NODE_DETECTION = 1
!
 LINE_NODES = ND-SET 7,8,9,10,11,12,13,14,15,16,17,18
 SURFACE_NODES = ND-SET 1,2,3,4,5,6
 VOLUME_NODES = ND-SET 19
 NMB_STEPS = 200
 REL_TOLERANCE = 0.1
 USE_SIMPLIFIED_APPROACH = 0
 PERFORM_INEXACT_PARALLEL_REG = 0
!===================================================================