suzanne.soy/FreeCAD
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

FEM: transient thermomech analysis (steadystate == true) needs density

Browse Source
This commit is contained in:
Bernd Hahnebach 2016-09-25 19:36:10 +01:00 committed by Yorik van Havre
parent 005f7f6bc8
commit 80c95aa555
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/Mod/Fem/FemInputWriterCcx.py
View File

@ -299,7 +299,7 @@ class FemInputWriterCcx(FemInputWriter.FemInputWriter):
f.write('** Materials\n') f.write('** Materials\n')
f.write('** written by {} function\n'.format(sys._getframe().f_code.co_name)) f.write('** written by {} function\n'.format(sys._getframe().f_code.co_name))
f.write('** Young\'s modulus unit is MPa = N/mm2\n') f.write('** Young\'s modulus unit is MPa = N/mm2\n')
if self.analysis_type == "frequency" or self.selfweight_objects: if self.analysis_type == "frequency" or self.selfweight_objects or (self.analysis_type == "thermomech" and not self.solver_obj.SteadyState):
f.write('** Density\'s unit is t/mm^3\n') f.write('** Density\'s unit is t/mm^3\n')
if self.analysis_type == "thermomech": if self.analysis_type == "thermomech":
f.write('** Thermal conductivity unit is kW/mm/K = t*mm/K*s^3\n') f.write('** Thermal conductivity unit is kW/mm/K = t*mm/K*s^3\n')
@ -312,7 +312,7 @@ class FemInputWriterCcx(FemInputWriter.FemInputWriter):
YM = FreeCAD.Units.Quantity(mat_obj.Material['YoungsModulus']) YM = FreeCAD.Units.Quantity(mat_obj.Material['YoungsModulus'])
YM_in_MPa = float(YM.getValueAs('MPa')) YM_in_MPa = float(YM.getValueAs('MPa'))
PR = float(mat_obj.Material['PoissonRatio']) PR = float(mat_obj.Material['PoissonRatio'])
if self.analysis_type == "frequency" or self.selfweight_objects: if self.analysis_type == "frequency" or self.selfweight_objects or (self.analysis_type == "thermomech" and not self.solver_obj.SteadyState):
density = FreeCAD.Units.Quantity(mat_obj.Material['Density']) density = FreeCAD.Units.Quantity(mat_obj.Material['Density'])
density_in_tonne_per_mm3 = float(density.getValueAs('t/mm^3')) density_in_tonne_per_mm3 = float(density.getValueAs('t/mm^3'))
if self.analysis_type == "thermomech": if self.analysis_type == "thermomech":
@ -327,7 +327,7 @@ class FemInputWriterCcx(FemInputWriter.FemInputWriter):
f.write('*MATERIAL, NAME=' + mat_name + '\n') f.write('*MATERIAL, NAME=' + mat_name + '\n')
f.write('*ELASTIC \n') f.write('*ELASTIC \n')
f.write('{0:.0f}, {1:.3f}\n'.format(YM_in_MPa, PR)) f.write('{0:.0f}, {1:.3f}\n'.format(YM_in_MPa, PR))
if self.analysis_type == "frequency" or self.selfweight_objects: if self.analysis_type == "frequency" or self.selfweight_objects or (self.analysis_type == "thermomech" and not self.solver_obj.SteadyState):
f.write('*DENSITY \n') f.write('*DENSITY \n')
f.write('{0:.3e}, \n'.format(density_in_tonne_per_mm3)) f.write('{0:.3e}, \n'.format(density_in_tonne_per_mm3))
if self.analysis_type == "thermomech": if self.analysis_type == "thermomech":