Users:General FEM Analysis/syntax
From Carat++ Public Wiki
!################################################################### !#### 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 !=================================================================== 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 !=================================================================== !******************************* !* 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 ! ID PART PROPERTY NODES ... 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 !=================================================================== BC-DIRICHLET 1 : RBE MASTERNODE = NODE 1 ! MASTERNODE IS NOT IN NODE LIST DOFS = DISP_X DISP_Y DISP_Z NODES = 2,3,4,5,6,7,8,9,10 !=================================================================== 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_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 ! -- 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-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 ! TYPE=LD-NODE 1 SPLIT=1 ! TYPE=LD-ELEM 100 ! 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 !=================================================================== !=================================================================== 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 |
