def drawElementValues(self, figElementValues):
"""Draws element values in specified figure"""
cfv.figure(figElementValues.nr)
cfv.clf()
cfv.draw_element_values(self.vonMises,
self.coords,
self.edof,
self.mesh.dofs_per_node,
self.mesh.el_type,
self.a,
draw_elements=True,
draw_undisplaced_mesh=False,
title="Example 06 effective stress")
def drawElementValues(self, figElementValues):
"""Draw element values in provided figure"""
cfv.figure(figElementValues.nr)
cfv.clf()
cfv.draw_element_values(self.results.elForces,
self.mesh.coords,
self.mesh.edof,
self.rect.dofsPerNode,
self.rect.elementType,
self.results.a,
draw_elements=True,
draw_undisplaced_mesh=False,
title="Effective Stress",
magnfac=1)
solver.addBC(rect.left_id, 0.0)
solver.addForceTotal(rect.top_id, -10e5, dimension=2)
results = solver.execute()
# ---- Visualise results ----------------------------------------------------
cfu.info("Drawing results...")
cfv.figure()
cfv.draw_geometry(rect.geometry(), title="Geometry")
cfv.figure()
cfv.draw_mesh(mesh.coords, mesh.edof, rect.dofs_per_node, rect.element_type,
filled=True, title="Mesh") #Draws the mesh.
cfv.figure()
cfv.draw_displacements(results.a, mesh.coords, mesh.edof, rect.dofs_per_node, rect.element_type,
draw_undisplaced_mesh=False, title="Displacements",
magnfac=1)
cfv.figure()
cfv.draw_element_values(results.el_forces, mesh.coords, mesh.edof, rect.dofs_per_node, rect.element_type, results.a,
draw_elements=True, draw_undisplaced_mesh=False,
title="Effective Stress", magnfac=1)
#cfv.colorBar().SetLabel("Effective stress")
cfu.info("Done drawing...")
cfv.show_and_wait()
es, et = cfc.planqs(ex[i,:], ey[i,:], ep, D, ed[i,:])
# Calc and append effective stress to list.
vonMises.append( sqrt( pow(es[0],2) - es[0]*es[1] + pow(es[1],2) + 3*es[2] ) )
## es: [sigx sigy tauxy]
# ---- Visualise results ----------------------------------------------------
cfu.info("Visualising...")
cfv.drawGeometry(g, draw_points=False, label_curves=True)
cfv.figure()
cfv.draw_element_values(vonMises, coords, edof, mesh.dofs_per_node, mesh.el_type, a, draw_elements=True, draw_undisplaced_mesh=False, title="Example 06 effective stress")
cfv.figure()
cfv.draw_displacements(a, coords, edof, mesh.dofs_per_node, mesh.el_type, draw_undisplaced_mesh=True, title="Example 06")
# Make use of attribute 'nodesOnCurve' in GmshMesher to draw some arrows on
# the right hand side of the mesh:
rightSideNodes = set()
# 4 and 5 are the IDs of the curves where we applied the forces.
for curveID in [4,5]:
# Get the nodes, without duplicates.
el_type=el_type,
filled=True,
title="Mesh") #Draws the mesh.
cfv.figure()
cfv.draw_displacements(a,
coords,
edof,
dofs_per_node,
el_type,
draw_undisplaced_mesh=False,
title="Displacements",
magnfac=25.0)
cfv.figure()
cfv.draw_element_values(von_mises,
coords,
edof,
dofs_per_node,
el_type,
a,
draw_elements=True,
draw_undisplaced_mesh=False,
title="Effective Stress",
magnfac=25.0)
cfv.get_color_bar().SetLabel("Effective Stress")
cfu.info("Done drawing...")
cfv.showAndWait()