"""
This is an empty docstring
"""

import math
from oofem_classes import Constraint, Force, Node, Element, Structure, Visualizer


struct = Structure()

lb = 15.0
r  = 457.2/2000
t  = 10.0/1000
a  = math.pi * ((r ** 2) - ((r-t) ** 2))
e  = 2.1e11
c1 = Constraint([False, False, False])
c2 = Constraint([True, True, False])
f  = Force([0, -20e3, -100e3])

n1 = struct.add_node(Node([0.0, 0.0, lb * ((2.0/3.0) ** 0.5)]))
n2 = struct.add_node(Node([0.0, lb/(3 ** 0.5), 0.0]))
n3 = struct.add_node(Node([-lb/2, -lb/(12.0 ** 0.5), 0.0]))
n4 = struct.add_node(Node([lb/2, -lb/(12.0 ** 0.5), 0.0]))

n1.set_force(f)
n2.set_constraint(c1)
n3.set_constraint(c1)
n4.set_constraint(c2)

struct.add_element(e, a, 0, 1)
struct.add_element(e, a, 0, 2)
struct.add_element(e, a, 0, 3)
struct.add_element(e, a, 1, 2)
struct.add_element(e, a, 2, 3)
struct.add_element(e, a, 3, 1)