def display_network_residual_forces(network,
display=True,
layer=None,
clear_layer=False,
scale=1.0,
color=(0, 255, 255),
attr_name='is_anchor'):
tol = rhino.get_tolerance()
guids = rhino.get_objects(name='{0}.force:residual.*'.format(network.attributes['name']))
rhino.delete_objects(guids)
if not display:
return
lines = []
for key, attr in network.vertices(True):
if attr[attr_name]:
continue
force = attr['rx'], attr['ry'], attr['rz']
start = network.vertex_coordinates(key)
end = [start[i] + scale * force[i] for i in range(3)]
length = distance_point_point(start, end)
arrow = 'end'
name = '{0}.force:residual.{1}'.format(network.attributes['name'], key)
if length < tol:
continue
lines.append({
'start' : start,
'end' : end,
'name' : name,
'color' : color,
'arrow' : arrow,
})
rhino.xdraw_lines(lines, layer=layer, clear=clear_layer)
def display_network_applied_loads(network,
display=True,
layer=None,
clear_layer=False,
scale=1.0,
color=(0, 0, 255)):
tol = rhino.get_tolerance()
guids = rhino.get_objects(name='{0}.force:load.*'.format(network.attributes['name']))
rhino.delete_objects(guids)
if not display:
return
lines = []
for key, attr in network.vertices(True):
load = attr['px'], attr['py'], attr['pz']
end = network.vertex_coordinates(key)
start = [end[i] - scale * load[i] for i in range(3)]
length = distance_point_point(start, end)
arrow = 'end'
name = '{0}.force:load.{1}'.format(network.attributes['name'], key)
if length < tol:
continue
lines.append({
'start': start,
'end' : end,
'name' : name,
'color': color,
'arrow': arrow,
})
rhino.xdraw_lines(lines, layer=layer, clear=clear_layer)
def display_network_axial_forces(network,
display=True,
layer=None,
clear_layer=False,
scale=1.0,
attr_name='f',
color_tension=(255, 0, 0),
color_compression=(0, 0, 255)):
"""Display the axial forces in the edges of a network.
Parameters:
network (compas.datastructures.network.Network):
The network object.
display (bool): Optional.
If ``True``, display the axial forces.
If ``False``, don't display the axial forces.
Default is ``True``.
layer (str): Optional.
The layer to draw in. Default is ``None``.
clear_layer (bool): Optional.
Flag for clearing the layer.
Default is ``False``.
scale (float): Optional.
The scale of the forces.
Default is ``1.0``.
attr_name (str): Optional.
The name of the edge attribute storing the force value.
Default is ``'f'``.
color_tension (tuple): Optional.
The color to use for tension forces.
Default is ``(255, 0, 0)``.
color_compression (tuple): Optional.
The color to use for compression forces.
Default is ``(0, 0, 255)``.
Example:
.. code-block:: python
import random
import compas
from compas.cad import rhino as compas_rhino
from compas.datastructures.network import Network
network = Network.from_obj(compas.get_data('lines.obj'))
for u, v, attr in network.edges(True):
attr['f'] = random.choice([-1.0, 1.0]) * random.randint(1, 10)
compas_rhino.display_network_axial_forces(network)
See Also:
* :func:`display_network_reaction_forces`
* :func:`display_network_residual_forces`
* :func:`display_network_selfweight`
"""
tol = rhino.get_tolerance()
objects = rhino.get_objects(name='{0}.force:axial.*'.format(network.attributes['name']))
rhino.delete_objects(objects)
if not display:
return
lines = []
for u, v, attr in network.edges(True):
start = network.vertex_coordinates(u)
end = network.vertex_coordinates(v)
force = attr['f']
color = color_tension if force > 0.0 else color_compression
radius = scale * ((force ** 2) ** 0.5 / 3.14159) ** 0.5
name = '{0}.force:axial.{1}-{2}'.format(network.attributes['name'], u, v)
if radius < tol:
continue
lines.append({
'start' : start,
'end' : end,
'name' : name,
'color' : color,
'radius': radius,
})
rhino.xdraw_cylinders(lines, layer=layer, clear=clear_layer)