# Input File for FEM Simulation with Mesh Reference

## Mesh mesh_file <name.msh> <physical_entity_final_mesh_name>
mesh_file Benchmarks/Elements/Benchmark_TO/Benchmark2_SERENDIPITYHEX_PARAM_M.inp Volume
rst_folder Benchmarks/Elements/Benchmark_TO/Results/
physics decoupledthermoelectromechanical
rst_folder Benchmarks/Elements/Benchmark_TO/Results/
solver NR
InitialTemperature 298.15

## MATERIAL PROPERTIES: material <volumename> \n <propertyname> <value> \n ... 
##
material Volume
Seebeck 1 8.2121E-05 3.5980E-07
Penalty_Seebeck 0 1
ThermalConductivity 1 2.9176 -0.0043
Penalty_ThermalConductivity 0 3
ElectricalConductivity 2 1.028E+5 -5.369E2 1.824
Penalty_ElectricalConductivity 0 3
ThermalExpansionCoefficient_x 0 15E-6
ThermalExpansionCoefficient_y 0 15E-6
ThermalExpansionCoefficient_z 0 15E-6
YoungModulus 0 61.6E9
Penalty_YoungModulus 0 3
PoissonRatio 0 0.241


## Boundary Conditions: bc <boundaryname> <surfacename> <value>
bc Voltage Ymax 0.0
bc Temperature Ymin 298
bc Displacement_x Ymin 0
bc Displacement_y Ymin 0
bc Displacement_z Ymin 0
bc Pressure_y Ymax -100000

## Topopoly Optimization Objective: TopOpt_Objective <ObjectiveName> <SurfaceName> <DesignElements>
TopOpt_Objective AverageTemperature Ymax Volume
## Topopoly Optimization Extra Variable (bc): TopOpt_Variable <bctype> <SurfaceName> <minval> <initialval> <maxval> 
TopOpt_bc Voltage Ymin -0.06 -0.058 -0.056
## Topopoly Optimization Intial densities: TopOpt_Initial_x <value(0,1](if<0, random values )>
TopOpt_Initial_x 0.8
## Topopoly Optimization NR start point: TopOpt_Initial_x <value(0,1](if<0, random values )>
TopOpt_NR T previous  
TopOpt_NR V proporcional
## Topopoly Optimization Constraints
TopOpt_Constraint Power 0.3
TopOpt_Constraint Volume 0.8
TopOpt_Constraint Stress_Pnorm 10000
KSU 3

## Output Options
## output <data> <path/outputname>
##output T dataT.vtk