示例#1
0
    def simulationCicle(self, **inputs):
#==============================================================================
#       Thermal model implementation (Every model should run on a loop)
#==============================================================================
        max_TC = max(self.mastic.thermal_conductivity,
                     self.airvoid.thermal_conductivity,
                     self.aggregate.thermal_conductivity)

        self.chemical = ChemicalModel(self.matrix_materials)
        self.chemical.applySimulationConditions(74.47)
        self.matrix_materials = self.chemical.simulate()

        self.mechanics = FEMMechanics(self.matrix_materials)
        self.mechanics.applySimulationConditions(inputs['force_input'])
        self.matrix_materials = self.mechanics.simulate()

        data1 = self.matrix_materials.copy()

        self.thermal = ThermalModel(self.matrix_materials, max_TC)
        self.thermal.applySimulationConditions()
        self.matrix_materials = self.thermal.simulate(inputs['thermal_steps'])

        self._calcNewModules(self.matrix_materials)
        # Change the EA in order to run the second chemical model

        self.chemical = ChemicalModel(self.matrix_materials)
        #A little change in the energy activation(EA) is emplemented
        # from a increase of 3.13 of the rca in three moths
        # it's obtained the following change in EA
        # 3.9700659563673262303399289700659563673262303399289700659*10e-7 each second
        # in 4 hours, that is, 144000 seconds, rca would increase 313/5475
        self.chemical.applySimulationConditions(313/5475.)
        self.matrix_materials = self.chemical.simulate()

        self.mechanics = FEMMechanics(self.matrix_materials)
        self.mechanics.applySimulationConditions(inputs['force_input'])
        self.matrix_materials = self.mechanics.simulate()

        data2 = self.matrix_materials.copy()

        return data1, data2
示例#2
0
class SimulationEngine(object):
    """
    This class configures the sample as a materials array in order to run the
    simulations defined in the thermal, mechanical and chemical models.
    """
    def __init__(self, collection, slice_id, **physical_cons):

        self.collection = collection

        # materials creation
        self.mastic = Material('mastic',
                               physical_cons['mastic_YM'], # young modulus
                               physical_cons['mastic_TC'], # thermal conducticity
                               physical_cons['mastic_CH']) # chemical

        self.aggregate = Material('aggregate',
                                  physical_cons['aggregate_YM'],
                                  physical_cons['aggregate_TC'],
                                  physical_cons['aggregate_CH'])

        self.airvoid = Material('airvoid',
                                physical_cons['air_YM'],
                                physical_cons['air_TC'],
                                physical_cons['air_CH'])

        vertical_slice = self._loadVerticalSlice(slice_id)

        # Structure data where the simulation takes place
        self.matrix_materials = self._getMatrixMaterials(vertical_slice)
        self.printToTxt()

    def printToTxt(self):
        r"""This function saves in a txt and jpg file the values of the Young's
        modulus for each material on the vertical slice selected by the user
        when the simulation is configured."""

        print "Exporting matrix materials..."
        # matrix_modulus = np.empty([self.matrix_materials.shape[0], self.matrix_materials.shape[1]], dtype=int)
        matrix_modulus = np.empty(self.matrix_materials.shape, dtype=int)

        # Creates the folder for the output
        if not os.path.exists("plain_images"):
            os.makedirs("plain_images")

        # Iterate over the slice to obtain the modulus from each node
        for i in xrange(self.matrix_materials.shape[0]):
            for j in xrange(self.matrix_materials.shape[1]):
                matrix_modulus[i, j] = self.matrix_materials[i, j].young_modulus

        # Saves the files to the plain_images folder
        np.savetxt('plain_images/materials.txt', matrix_modulus, delimiter=',', fmt='%i')
        imsave('plain_images/materials.jpg', matrix_modulus)


    def _loadVerticalSlice(self, slice_id):
        """Cut the slice of the collection in the position id"""
        vertical_slice = self.collection[:, :, slice_id]
        vertical_slice = vertical_slice.copy() # avoid side effects
        return vertical_slice.transpose()

    def _getMatrixMaterials(self, vertical_slice):
        """Create the matrix material from a vertical slice"""
        material_matrix = np.empty(vertical_slice.shape, dtype=object)

        for (x,y), _ in np.ndenumerate(vertical_slice):
            if vertical_slice[x, y] == 2:
                material_matrix[x,y] = copy.deepcopy(self.aggregate)
            elif vertical_slice[x,y] == 1:
                material_matrix[x,y] = copy.deepcopy(self.mastic)
            elif vertical_slice[x,y] == 0 or  vertical_slice[x,y] == -1:
                material_matrix[x,y] = copy.deepcopy(self.airvoid)

        print "Materials matrix created, size:", material_matrix.shape
        return material_matrix

    def _calcNewModules(self, MM):
        print "Recalculating Young's modules..."
        for i in xrange(MM.shape[0]):
            for j in xrange(MM.shape[1]):
                if MM[i,j].phase == 'mastic':
                    if MM[i,j].temperature <= 20:
                        MM[i,j].young_modulus = 16030

                    elif MM[i,j].temperature <= 35:
                        MM[i,j].young_modulus = 5148

                    else:
                        MM[i,j].young_modulus = 1527

    def simulationCicle(self, **inputs):
#==============================================================================
#       Thermal model implementation (Every model should run on a loop)
#==============================================================================
        max_TC = max(self.mastic.thermal_conductivity,
                     self.airvoid.thermal_conductivity,
                     self.aggregate.thermal_conductivity)

        self.chemical = ChemicalModel(self.matrix_materials)
        self.chemical.applySimulationConditions(74.47)
        self.matrix_materials = self.chemical.simulate()

        self.mechanics = FEMMechanics(self.matrix_materials)
        self.mechanics.applySimulationConditions(inputs['force_input'])
        self.matrix_materials = self.mechanics.simulate()

        data1 = self.matrix_materials.copy()

        self.thermal = ThermalModel(self.matrix_materials, max_TC)
        self.thermal.applySimulationConditions()
        self.matrix_materials = self.thermal.simulate(inputs['thermal_steps'])

        self._calcNewModules(self.matrix_materials)
        # Change the EA in order to run the second chemical model

        self.chemical = ChemicalModel(self.matrix_materials)
        #A little change in the energy activation(EA) is emplemented
        # from a increase of 3.13 of the rca in three moths
        # it's obtained the following change in EA
        # 3.9700659563673262303399289700659563673262303399289700659*10e-7 each second
        # in 4 hours, that is, 144000 seconds, rca would increase 313/5475
        self.chemical.applySimulationConditions(313/5475.)
        self.matrix_materials = self.chemical.simulate()

        self.mechanics = FEMMechanics(self.matrix_materials)
        self.mechanics.applySimulationConditions(inputs['force_input'])
        self.matrix_materials = self.mechanics.simulate()

        data2 = self.matrix_materials.copy()

        return data1, data2