Users:General FEM Analysis/syntax

From Carat++ Public Wiki
Jump to: navigation, search
!===================================================================

!===================================================================
PC-ANALYSIS 1: STA_GEO_LIN
 SOLVER = PC-SOLVER 1
 OUTPUT = PC-OUT 1
 COMPCASE = LD-COM 1
 DOMAIN = EL-DOMAIN 1
!-------
! SOLVE_EIGEN_PROBLEM = 0
! EIGEN_SOLVER = PC-SOLVER 2
!===================================================================
!===================================================================
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
!==================================================================
!*******************************
!*      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
! DIAGNOSTICS = 1
! BANDWITH = CUTHILL_MCKEE
!===================================================================
PC-SOLVER 7: MKL_PARDISO
 BANDWITH = CUTHILL_MCKEE
 MATRIX_TYPE = FULL   !SYMMETRIC or FULL -> Default Vaule 'SYMMETRIC'
 DATA_STORAGE = IN_CORE   ! IN_CORE  OUT_OF_CORE  SWAP
!===================================================================
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
!===================================================================
ND-SET 2 NAME = set2
  PART 1   NODE_TYPE = ALL   ! ALL  CORNER  LINE  SURFACE  VOLUME   
!===================================================================
!
!
!###################################################################
!####                          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
!===================================================================
!***********************
!*      MATERIAL       *
!***********************
!===================================================================
EL-MAT 2 : ST_VENANT_KIRCHHOFF_ISOTROPIC
EMOD=20000  DENS=0.0 NUE=0.2
!===================================================================
!***********************
!*     TOPOLOGY        *
!***********************
EL-SET 1 NAME=ALL_ELEMS_OF_PART_1
  PART 1
!==================================================================  
EL-SET 1 NAME=ELEMENTS_BY_ID
  NEL 1 2 3
!==================================================================  
!
!
!###################################################################
!####                          LD-BLOCK                         ####
!###################################################################

BC-DIRICHLET 2 : MPC-COUPLING
 DOFS = DISP_Z   NODES = 550,551,553,554,555,556,557,558,559,560,561
!===================================================================
BC-DIRICHLET 2: SPC-NONZERO
 NODE    1002  DISP_Z=0.1  CURVE=LD-CURVE 1
 NODE    1003  DISP_Z=0.1  CURVE=LD-CURVE 1
!===================================================================
LD-NODE 1 TYPE=FORCE
  NODE 481  DIRECTION = NODE_DIR_REF   VAL=-1
!===================================================================
LD-ELEM 1 PART=1   
! CURVE=LD-CURVE 1
 TYPE=DEAD      D1=0.0   D2=0.0   D3=1.0  VAL=100.0
!===================================================================
LD-ELEM 2
 TYPE = TEMPERATURE
!  CURVE=LD-CURVE 1
 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-BOUND 1 : AVOID_SELF_PENETRATION
  MIN_JACOBIAN =  1e-3  ! 3e-4 ! 1e-4 ! 5e-6
  MAX_ITER = 30
!===================================================================
OPT-BOUND 1 : SOLID_THICKNESS
 MAX_THICKNESS = 1000
 MIN_THICKNESS = 0.1
!===================================================================
!
!
!
!===================================================================
OPT-RESPONSE_FCT 1 : STRAIN_ENERGY_NONLIN
! -- basic stuff
 WEIGHT=1.0 ANALYSIS=PC-ANALYSIS 2 ETA=1e-06
 GRAD=ADJOINT SA=SEMI_ANALYTIC FDA=FOREWARD
 LOAD_CONSTANT = 1
 DESVAR=OPT-VAR 1
 !===================================================================
OPT-RESPONSE_FCT 1 : STRAIN_ENERGY_LIN
! -- basic stuff
 WEIGHT=1.0 ANALYSIS=PC-ANALYSIS 1 ETA=1e-06
 GRAD=ADJOINT SA=SEMI_ANALYTIC FDA=FOREWARD
 DESVAR=OPT-VAR 1
 !===================================================================
OPT-RESPONSE_FCT 1 : DISPLACEMENT_LIN
 WEIGHT=1.0 ANALYSIS=PC-ANALYSIS 1 ETA=1e-06
 GRAD=ADJOINT SA=SEMI_ANALYTIC 
 FDA= CENTRAL ! FOREWARD
 DESVAR=OPT-VAR 2
 NODE = 552
 DOF = DISP_Y
 GOAL = MAX
 USE_ABS_VALUE = 1  !   1: YES   0: NO 
!===================================================================
OPT-RESPONSE_FCT 2 : DISPLACEMENT_NONLIN
! -- basic stuff
 WEIGHT=1.0 ANALYSIS=PC-ANALYSIS 2 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
 USE_ABS_VALUE = 1  !   1: YES   0: NO 
! -- 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
!===================================================================
OPT-RESPONSE_FCT 10 :  NLN_BUCKLING
  WEIGHT=1.0 ANALYSIS=PC-ANALYSIS 3 ETA=1e-06
 GRAD=ADJOINT 
 SA=SEMI_ANALYTIC ! EXACT_SEMI_ANALYTIC  GLOBAL_FD
 FDA=FOREWARD
 DESVAR=OPT-VAR 1
!----------
 GOAL = MIN
!===================================================================
!            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-RESPONSE_FCT 2 : SUM  !  response computes sum of design variables
 WEIGHT=1.0   DESVAR=OPT-VAR 1 
! -- dummys -- 
 ETA=1e-05
 GRAD=ADJOINT SA=SEMI_ANALYTIC FDA=FOREWARD 
!===================================================================
OPT-RESPONSE_FCT 3 : COST_BENEFIT  !  cost-benefit of two pre-defined reponses 
 COST_RESPONSE    =  OPT-RESPONSE_FCT 2
 BENEFIT_RESPONSE =  OPT-RESPONSE_FCT 1
 WEIGHT=1.0  
! -- dummys -- 
 DESVAR=OPT-VAR 1   ETA=1e-05
 GRAD=ADJOINT SA=SEMI_ANALYTIC FDA=FOREWARD 
!===================================================================
OPT-RESPONSE_FCT 10 : LEAST_SQUARES
 SUB_RESPONSE = OPT-RESPONSE_FCT     1,  2,    3 !, 4
 TARGET_VALUE =                    -10, -3,  -30 !, -0.5
 GOAL = MAX ! MIN MAX
! dummy values ...
  DESVAR=OPT-VAR 1
  GRAD=ADJOINT 
  SA=EXACT_SEMI_ANALYTIC
  FDA=FOREWARD  
!===================================================================
OPT-RESPONSE_FCT 3 : LOAD_FACTOR_DISP_CTRL  !  STRAIN_ENERGY_NONLIN_DISP_CTRL
 WEIGHT=1.0 ANALYSIS=PC-ANALYSIS 1 ETA=1e-06
 GRAD=ADJOINT 
 SA= SEMI_ANALYTIC  !  SEMI_ANALYTIC    GLOBAL_FD
 FDA=FOREWARD
 DESVAR=OPT-VAR 1
 LOAD_CONSTANT = 1
 GOAL = MIN
 USE_ABS_VALUE = 0 !  1  0
! -- constraint parameters
! REL_LIMIT = 1.00
! REL_TOLERANCE = 0.025
! CONSTRAINT_TYPE = INEQUALITY_LT   ! 'EQUALITY' or  'INEQUALITY_LT'
! LAMBDA_ABS_MAX = 40
!===================================================================
OPT-RESPONSE_FCT 6 : RESPONSE_SUM
 SUB_RESPONSE = OPT-RESPONSE_FCT     4,  5
 SUB_WEIGHTING_FTRS =                1, -1
 GOAL = MIN ! MIN MAX
! dummy values ...
  DESVAR=OPT-VAR 1
  GRAD=ADJOINT 
  SA=EXACT_SEMI_ANALYTIC
  FDA=FOREWARD  
!===================================================================
!
!
!
!===================================================================
OPT-FILTER 1 : TOPO   !   GLOBAL   TOPO
 FUNCTIONS = OPT-LINKFUNC 1                             ! FILTER ALL DESIGN VARIABLES
! FUNCTION = OPT-LINKFUNC 1          DESVAR=OPT-VAR 1   ! FILTER ONLY STECIFIED DESVAR BLOCK
!===================================================================
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=ND-SET 1 SUBTYPE=DIR
 TYPE=ND-SET 1 SUBTYPE=FIXED_DIR  
 TYPE=ND-SET 2 SUBTYPE=MAPPED_FIXED_DIR    D1_1=0  D1_2=1  D1_3=0   D2_1=1  D2_2=0  D2_3=0
!
 TYPE=EL-SET 1 SUBTYPE=THICKNESS 
 TYPE=EL-SET 1 SUBTYPE=THICKNESS LAYER_ID = 1  MODE=SINGLE   ! SINGLE  PATCH
!
 TYPE=EL-SET 1 SUBTYPE=MATERIAL_ANGLES
 TYPE=EL-SET 1 SUBTYPE=MATERIAL_ANGLE        LAYER_IDS = 1,2,3
!
 TYPE=LD-NODE 1  SPLIT=1  ! 1: consider each node as a single variable  0:block mode     
!
 TYPE=LD-ELEM 1  SPLIT=1  ! 1: consider each node as a single variable  0:block mode  
!===================================================================
!
!
!
!===================================================================
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
!
 PERFORM_RESTART = 0 ! 0  1
!********** 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
!===================================================================
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
!
! OUTPUT = PC-OUT 1
! OUTPUT_FREQUENCY = 10
!
 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
 APPROACH = SIMPLE   !  EXACT  EQC
 PERFORM_INEXACT_PARALLEL_REG = 0
!===================================================================





Whos here now:   Members 0   Guests 1   Bots & Crawlers 0
 
Personal tools
Content for Developers