def __init__(self, beamMesh, backgroundMesh, beamThickness, initialGap):
## Read in 2d Solid Mesh
self.beam_thickness = beamThickness
self.Gap = initialGap
beamReader = FluentCase(beamMesh)
beamReader.read()
print "read solid mesh"
self.solidMeshes = beamReader.getMeshList()
self.geomFields = models.GeomFields('geom')
self.solidMetricsCalculator = models.MeshMetricsCalculatorA(
self.geomFields, self.solidMeshes)
self.solidMetricsCalculator.init()
## Define plate and deformation model
self.plateFields = models.PlateFields('plate')
self.pmodel = models.PlateModelA(self.geomFields, self.plateFields,
self.solidMeshes)
self.dmodel = models.PlateDeformationModelA(self.geomFields,
self.plateFields,
self.solidMeshes)
bcMap = self.pmodel.getBCMap()
## Apply a default Boundary Condition
#for i, bc in bcMap.iteritems():
#bc.bcType = 'SpecifiedTraction'
## Read in 3d Background Mesh
fluidReader = FluentCase(backgroundMesh)
fluidReader.read()
self.fluidMeshes = fluidReader.getMeshList()
self.fluidMetricsCalculator = models.MeshMetricsCalculatorA(
self.geomFields, self.fluidMeshes)
self.fluidMetricsCalculator.init()
## Define electric model
self.elecFields = models.ElectricFields('elec')
self.emodel = models.ElectricModelA(self.geomFields, self.elecFields,
self.fluidMeshes)
bcMap = self.emodel.getBCMap()
## Apply Default boundary conditions
for i, bc in bcMap.iteritems():
bc.bcType = "Symmetry"
self.solidBoundaryMeshes = [
m.extrude(1, beamThickness, True) for m in self.solidMeshes
]
self.solidBoundaryMetricsCalculator = models.MeshMetricsCalculatorA(
self.geomFields, self.solidBoundaryMeshes)
self.solidBoundaryMetricsCalculator.init()
def createModels(self):
if self.enablePlateModel == True:
self.plateFields = models.PlateFields('plate')
self.pmodel = models.PlateModelA(self.geomFields, self.plateFields,
self.solidMeshes)
self.dmodel = models.PlateDeformationModelA(
self.geomFields, self.plateFields, self.solidMeshes)
if self.enableElecModel == True:
self.elecFields = models.ElectricFields('elec')
self.emodel = models.ElectricModelA(self.geomFields,
self.elecFields,
self.fluidMeshesNew)
if self.enableFlowModel == True:
self.flowFields = models.FlowFields('flow')
self.fmodel = models.FlowModelA(self.geomFields, self.flowFields,
self.fluidMeshes)
#currently contact model is implemented in ComputeForce.py
#if self.enableContactModel == True:
# self.contactFields = models.ContactFields('contact')
# self.cmodel = models.ContactModelA(self.geomFields,self.contactFields,self.fluidMeshesNew)
print 'models are created'
beamReader = FluentCase(fileBase_input + "beam_2D.cas")
beamReader.read()
solidMeshes = beamReader.getMeshList()
solidBoundaryMeshes = [m.extrude(1, beam_thickness, True) for m in solidMeshes]
solidBoundaryMetricsCalculator = models.MeshMetricsCalculatorA(
geomFields, solidBoundaryMeshes)
solidBoundaryMetricsCalculator.init()
### output files
probeFile = open(fileBase_output + "centerDisplacement.dat", "w")
forceFile = open(fileBase_output + "force.dat", "w")
### electric model and boundary condition ###
elecFields = models.ElectricFields('elec')
emodel = models.ElectricModelA(geomFields, elecFields, fluidMeshes)
bcMap = emodel.getBCMap()
### specify a potential field; use it as boundary condition
pf = fvmbaseExt.Field('potential')
fgs = fluidMeshes[0].getBoundaryFaceGroups()
for fg in fgs:
if fg.id == fluidBot:
faces = fg.site
nFaces = faces.getCount()
areaMag = geomFields.areaMag[faces]
xf = geomFields.coordinate[faces].asNumPyArray()
pot = areaMag.newSizedClone(nFaces)
