def displayBlocks(self, setToDisplay= None, caption= None, fileName=None):
'''Display the blocks (points, lines, surfaces and volumes)
of the set.
:param setToDisplay: set to display.
:param caption: title of the graphic.
:param fileName: name of the file to plot the graphic. Defaults to
None, in that case an screen display is generated
'''
if(setToDisplay==None):
setToDisplay= self.modelSpace.getTotalSet()
if(caption==None):
caption= setToDisplay.name+' set; blocks'
displaySettings= vtk_CAD_graphic.DisplaySettingsBlockTopo()
displaySettings.cameraParameters= self.getCameraParameters()
displaySettings.displayBlocks(setToDisplay,caption= caption, fileName= fileName)
CooMaxY= 1 old_stderr= sys.stderr # Problem type feProblem= xc.FEProblem() preprocessor= feProblem.getPreprocessor nodes= preprocessor.getNodeHandler modelSpace= predefined_spaces.SolidMechanics2D(nodes)# Model space definition # Definition of points. points= preprocessor.getMultiBlockTopology.getPoints pt1= points.newPoint(geom.Pos3d(0.0,0.0,0.0)) pt2= points.newPoint(geom.Pos3d(CooMaxX/2,0.0,0.0)) pt3= points.newPoint(geom.Pos3d(CooMaxX,0.0,0.0)) pt4= points.newPoint(geom.Pos3d(0.0,CooMaxY,0.0)) pt5= points.newPoint(geom.Pos3d(CooMaxX/2,CooMaxY,0.0)) pt6= points.newPoint(geom.Pos3d(CooMaxX,CooMaxY,0.0)) # Definition of quadrilateral surfaces. surfaces= preprocessor.getMultiBlockTopology.getSurfaces s1= surfaces.newQuadSurfacePts(pt1.tag,pt2.tag,pt5.tag,pt4.tag) s2= surfaces.newQuadSurfacePts(pt2.tag,pt3.tag,pt6.tag,pt5.tag) # Display surfaces. from postprocess.xcVtk.CAD_model import vtk_CAD_graphic displaySettings= vtk_CAD_graphic.DisplaySettingsBlockTopo() setDisp= feProblem.getPreprocessor.getSets['total'] displaySettings.displayBlocks(setDisp,caption= setDisp.name+' set')
