Users:General FEM Analysis/syntax
From Carat++ Public Wiki
!=================================================================== !=================================================================== 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 |