Users:General FEM Analysis/syntax

From Carat++ Public Wiki
< Users:General FEM Analysis
Revision as of 11:43, 1 October 2012 by Helmut (Talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

!################################################################### !#### 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_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 

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 = 5
RHO = 25

! -- constraint parameters 

ABS_LIMIT = 12.0
REL_TOLERANCE = 0.05
CONSTRAINT_TYPE = INEQUALITY_LT   ! 'EQUALITY' or  'INEQUALITY_LT'
LAMBDA_ABS_MAX = 6

!=================================================================== !=================================================================== 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

!===================================================================





Whos here now:   Members 0   Guests 0   Bots & Crawlers 1
 
Retrieved from "http://carat.st.bv.tum.de/caratuserswiki/index.php/Users:General_FEM_Analysis/syntax"
Personal tools
Content for Users
Content for Developers