; # PLATE ; [InputMesh] ; TriangleType = PlateElement ; ClampedBoundaryCondition = 1 ; chartXAType = id ; ; ShellType = 100 ; file = ../../../../data/VTK/unitSquare/unitSquareSubdiv4.vtk ; tangentSpaceType = 2 ; parserFileVTKPlotChart = ../../../../parameters/dkt/ExampleConfigurationsPlate/vtkPlotInterface_Plate_Undeformed.ini ; parserFileVTKPlotUndeformed = ../../../../parameters/dkt/ExampleConfigurationsPlate/vtkPlotInterface_Plate_Undeformed.ini ; parserFileVTKPlotDeformed = ../../../../parameters/dkt/ExampleConfigurationsPlate/vtkPlotInterface_Plate_Bdr-Left.ini ; parserFileVTKPlotStress = ../../../../parameters/dkt/ExampleConfigurationsPlate/vtkPlotInterface_Plate_Stress.ini ; ; [ConstraintProblemAdaptiveProlongationTypes] ; ProlongationMeshType = 1 ; ProlongationPfType = 1 ; ProlongationDispType = 1 ; ProlongationMeshBoundaryType = 1 ; constructProlongatedBoundaryByShellHandler = 0 ; ; [Force] ; numLoads = 1 ; Type1 = 11 ; ; Factor1 = -1. ; ; StartRange1 = -100.,-100.,-100. ; ; EndRange1 = 100.,100.,100. ; Factor1 = -250. ; StartRange1 = 0.45,0.45,-100. ; EndRange1 = 0.55,0.55,100. ; scaleForceWithThickness = 0 ; scaleForceWithThicknessExponent = 0 ; factor_force = 1. #HALF CYLINDER ; [InputMesh] ; TriangleType = PlateElementMappedToShell ; ClampedBoundaryCondition = 1 ; chartXAType = PlateToCylinder ; ; ShellType = 5 ; file = ../../../../data/VTK/unitSquare/unitSquareSubdiv4.vtk ; tangentSpaceType = 2 ; parserFileVTKPlotChart = ../../../../parameters/dkt/ExampleConfigurationsPlate/vtkPlotInterface_Plate_Undeformed.ini ; parserFileVTKPlotUndeformed = ../../../../parameters/dkt/ExampleConfigurationsCylinder/vtkPlotInterface_Cylinder_Undeformed.ini ; parserFileVTKPlotDeformed = ../../../../parameters/dkt/ExampleConfigurationsCylinder/vtkPlotInterface_Cylinder.ini ; parserFileVTKPlotStress = ../../../../parameters/dkt/ExampleConfigurationsCylinder/vtkPlotInterface_Cylinder_Stress.ini ; ; [ConstraintProblemAdaptiveProlongationTypes] ; ProlongationMeshType = 1 ; ProlongationPfType = 1 ; ProlongationDispType = 1 ; ProlongationMeshBoundaryType = 1 ; constructProlongatedBoundaryByShellHandler = 0 ; ; [Force] ; numLoads = 1 ; Type1 = 2 ; RadiusBall1 = 0.01 ; CenterBall1 = 0.31830988618379067153776752674502872,0.5,0.0 ; Load1 = -1.,0.0,0.0 ; scaleForceWithThickness = 1 ; scaleForceWithThicknessExponent = 1.0 #HALF CYLINDER, TOPBOTTOM [InputMesh] TriangleType = PlateElementMappedToShell ClampedBoundaryCondition = 1 chartXAType = PlateToCylinder ShellType = 6 file = ../../../../data/VTK/unitSquare/unitSquareSubdiv4.vtk tangentSpaceType = 2 parserFileVTKPlotChart = ../../../../parameters/dkt/examplePlate/vtkPlotInterface_Plate_Undeformed.ini parserFileVTKPlotUndeformed = ../../../../parameters/dkt/exampleCylinder/vtkPlotInterface_Cylinder_Undeformed.ini parserFileVTKPlotDeformed = ../../../../parameters/dkt/exampleCylinder/vtkPlotInterface_Cylinder.ini [ConstraintProblemAdaptiveProlongationTypes] ProlongationMeshType = 1 ProlongationPfType = 1 ProlongationDispType = 1 ProlongationMeshBoundaryType = 1 constructProlongatedBoundaryByShellHandler = 0 [Force] numLoads = 1 Type1 = 1 Load1 = -100.,0.0,0.0 StartRange1 = -100,-100,-100. EndRange1 = 100.,100.,100. scaleForceWithAreaOfHardMaterial = 1 scaleForceWithThickness = 0 scaleForceWithThicknessExponent = 1. factor_force = 1. #HALF SPHERE ; [InputMesh] ; TriangleType = PlateElementMappedToShell ; ClampedBoundaryCondition = 1 ; chartXAType = CircToSphere ; ; ShellType = 100 ; file = ../../../../data/VTK/circle/circ8_radius1_LoopSubdiv3.vtk ; tangentSpaceType = 2 ; parserFileVTKPlotChart = ../../../../parameters/dkt/ExampleConfigurationsPlate/vtkPlotInterface_Plate_Undeformed.ini ; parserFileVTKPlotUndeformed = ../../../../parameters/dkt/ExampleConfigurationsHalfSphere/vtkPlotInterface_HalfSphere_Undeformed.ini ; parserFileVTKPlotDeformed = ../../../../parameters/dkt/ExampleConfigurationsHalfSphere/vtkPlotInterface_HalfSphere.ini ; parserFileVTKPlotStress = ../../../../parameters/dkt/ExampleConfigurationsHalfSphere/vtkPlotInterface_HalfSphere_Stress.ini ; ; [ConstraintProblemAdaptiveProlongationTypes] ; ProlongationMeshType = 12 ; ProlongationPfType = 1 ; ProlongationDispType = 1 ; ProlongationMeshBoundaryType = 1 ; constructProlongatedBoundaryByShellHandler = 0 ; ; [Force] ; numLoads = 1 ; Type1 = 2 ; ; RadiusBall1 = 0.05 ; CenterBall1 = 0.0,0.0,1.0 ; Load1 = 0.,0.,-4.0 ; scaleForceWithThicknessExponent = 1.0 ; scaleForceWithThickness = 1 ; ; ; RadiusBall1 = 0.1 ; ; CenterBall1 = 0.0,0.0,1.0 ; ; Load1 = 0.,0.,-1.0 ; ; scaleForceWithThicknessExponent = 1.2 ; ; scaleForceWithThickness = 1 ; ; factor_force = 1. ;######################## [Material] ; #------------------------------------------------------------------------ ; ##### elastic energy parameter ##### ; #------------------------------------------------------------------------ MaterialFromDictionary = 0 # elastic materials: #if = 0 ElastModulus_hard = 1.e+2 ElastModulus_soft = 1.e+0 PoissonRatio_hard = 0. PoissonRatio_soft = 0. thickness_hard = 1. thickness_soft = 1. factor_membraneEnergy = 1.e-0 ; factor_bendingEnergy = 1.e-2 ; factor_bendingEnergy = 1.e-3 factor_bendingEnergy = 1.e-4 ; factor_bendingEnergy = 1.e-6 ; factor_bendingEnergy = 1.e-8 #chi(v) = ...: -1 - harmonic, 1,2,4 - polyonmial pfFunctionMaterialType = 4 #double well: 2,4 - polyonmial pfFunctionDoubleWellType = 4 [InitMaterial] #initMaterialType: # 1 - constant, # 2 - layer # 3 - ray # 4 - triangle # 10 - holes # 1000 Random ;######################## [saving] saveDirectory = ../../../../results/resultsDKT saveInitialization = 1 saveSolution = 1 plotResults = 0 removeOldRefinementSteps = 1 plotPointCloud = 0 VTKFileType = vtk ; VTKFileType = vtp LATEXType = pdflatex ; LATEXType = lualatex plotInterface = 1 plotStress = 1 thresholdInterface = 0.5 thresholdBdr = 1. thresholdHardRegion = 0.5 thresholdSoftRegion = 0.5 ;######################## [DerivativeTest] # order: 0 - no, 1 - FirstDerivativeTest, 2 - SecondDerivativeTest order = 0 stepSize = 1.e-7 tolerance = 1.e-5 ;######################## [DerivativeTestConstraint] # order: 0 - no, 1 - FirstDerivativeTest, 2 - SecondDerivativeTest, 12 - both order = 0 stepSize = 1.e-7 tolerance = 1.e-5 ;######################## [ConstraintProblem] ; #------------------------------------------------------------------------ ; ##### constraint problem (linelast ) ##### ; #------------------------------------------------------------------------ solveWithDirectSolver = 1 #if iterative solver: #tolerance = sqrt( tol^2 * systemmat.col() ) toleranceLinearSystem = 1.e-16 #maxIters = maxItersFac * systemmat.cols() maxItersFacLinearSystem = 0.25 # linear: 0 - Dirichlet, 1 - Laplace, 2 - DirichletLaplace, 3 - KirchhoffLove, # nonlinear:11 -NonlinMembraneHinge ElastEnergy = 31 #nonlinear solver: NEWTON or IPOPT or FirstTryNEWTONelseIPOPT ; nonlinearSolver = NEWTON nonlinearSolver = IPOPT ; nonlinearSolver = FirstTryNEWTONelseIPOPT ############### for ipopt # 0: MUMPS, x HSL_x x = 27,57,77,86,97 ; linearSolverTypeIpopt = 0 linearSolverTypeIpopt = 86 maxIterationsIpopt = 200 breakConditionIpopt = 1.e-10 ; breakConditionIpopt = 1.e-12 outputLevelIpopt = 0 ############### for newton # 2 - Cholmod, 1 Umfpack, 24 BICGSTAB linearSolver = 1 maxIterations = 200 breakCondition = 1.e-14 #set initialization of displacement and multiplier to solution for current material updateInitialization = 1 # 0 - no output , 1 - output if Newton fails, 2 - , 3 - Newton steps, 4 - always plot deformation when solve #recommended: 1 for material optimization, 3 for optimal deformation outputLevel = 1 # initType : 0 - zero, # Plate: 1 - simple, 2 - twist # general: 101 : h * 10^-3 (1,1,1) on interior nodes # 102 : h * 10^-3 f on interior nodes initType = 0 # MarkingType: # 1 - all # 2 - all twice # 11 - L2Gradient of material # 12 - L2Gradient and clamped elements (ie if all nodes on boundary) # 31 - L2Hessian of disp # 32 - L2Hessian of disp and clamped elements (ie if all nodes on boundary) ; adaptiveMarkingType = 32 adaptiveMarkingType = 2 adaptiveMarkingPercentil = 0.85 ; adaptiveMarkingPercentil = 0.75 ; adaptiveMarkingPercentil = 0.5 ; adaptiveMarkingPercentil = 0.15 adaptiveMarkingUseTotalStressVec = 1 ;######################## [MaterialOptimization] factorComplianceCost = 1.0 factorAreaCost = 0.0 factorInterfaceCost = 0.0 factorDoubleWell = 1. [Compare] thicknessVec = 0.01,0.003,0.001 ; thicknessVec = 0.01 factorForceVec = 1.,2.,4.,8. ; thicknessVec = 0.1,0.03,0.01,0.003,0.001 ; factorForceVec = 1.,2.,3.,4.,5.,6.,7.,8.,9.,10.,11.,12.,13.,14.,15.,16.,17.,18.,19.,20.,21.,22.,23.,24. ; factorForceVec = 1.,2.,3.,4.,5.,6.,7.,8.,9.,10.,11.,12.,13.,14.,15.,16.,17.,18.,19.,20.,21.,22.,23.,24.,25.,26.,27.,28.,29.,30.,31.,32. ; factorForceVec = 1.,2.,4.,8.,16.,32.,64.,128. ; checkScaledForce = 0 ; checkScaledForceForFinalSolution = 1 ; checkScaledForceForFinalSolutionVec = 2.,3.,4. ; checkScaledForceForFinalSolutionVec = 0.5,1.,1.5,2.,3.,4.,5.,6.,7. ; checkScaledForceForFinalSolutionVec = 1.,1.5,2.,3.,4.,6.,8.,12.,16. ; checkScaledForceForFinalSolutionVec = 1.,1.5,2.,3.,4.,6.,8.,12.,16.,24.,32.,48.,64.,96.,128.,192.,256.,382.,512. ; thicknessVec = 0.1,0.09,0.08,0.07,0.06,0.05,0.04,0.03,0.02,0.01,0.009,0.008,0.007,0.006,0.005,0.004,0.003,0.002,0.001 ############################## #for compareBound numBisections = 10 #possible quantities to compare: LInf, Estored valueToCompare = LInf LInfBound = 0.3 EstoredBound = 1.e-2 ;######################## [BASH] numThreads = 4 renewKerberosTicket = 0