polycalc APIpolycalc.online ↗

"""Walkthrough: PolycalcOnline SDK — structural engineering calculations."""
import json
from parse_apis.polycalc_online_api import PolycalcOnline, InputFormatInvalid

client = PolycalcOnline()

# --- Beam analysis: simply supported beam with a point load ---
supports = json.dumps([
    {"type": "pin", "position": 0},
    {"type": "roller", "position": 10},
])
point_loads = json.dumps([{"position": 5, "magnitude": 1000}])

beam = client.beam_analyses.calculate(
    beam_length=10,
    supports=supports,
    point_loads=point_loads,
    query_x=5,
)
print(f"Beam span: {beam.beam_length}")
for r in beam.reactions:
    print(f"  {r.support_type} @ {r.position}: Fy={r.force_y}")

# Query point gives shear & moment at the requested location
if beam.query_point is not None:
    qp = beam.query_point
    print(f"  At x={qp.x}: shear={qp.shear}, moment={qp.moment}")

# --- Mohr's circle: plane stress transformation ---
mohr = client.mohrs_circles.calculate(sigma_x=100, sigma_y=-50, tau_xy=30)
print(f"Principal stresses: σ1={mohr.principal_stresses.sigma_1}, σ2={mohr.principal_stresses.sigma_2}")
print(f"Max shear: τ_max={mohr.maximum_shear.tau_max}")
print(f"Circle center={mohr.circle_properties.center}, radius={mohr.circle_properties.radius}")

# --- Centroid of a composite section ---
shapes = json.dumps([
    {"type": "rectangle", "width": 12, "height": 8, "center_x": 6, "center_y": 4},
    {"type": "circle", "radius": 2, "center_x": 6, "center_y": 4, "hole": True},
])

try:
    section = client.cross_sections.calculate(shapes=shapes)
except InputFormatInvalid as e:
    print(f"Invalid input: {e.message}")
else:
    print(f"Net area: {section.net_area}")
    print(f"Centroid: ({section.centroid_x}, {section.centroid_y})")
    print(f"Ix={section.ix}, Iy={section.iy}")
    for part in section.parts:
        label = "(hole)" if part.hole else "(solid)"
        print(f"  {part.type} {label}: area={part.area}")

# --- Truss analysis: simple triangle truss ---
nodes = json.dumps([
    {"x": 0, "y": 0, "support": "pin"},
    {"x": 4, "y": 0, "support": "roller"},
    {"x": 2, "y": 3},
])
members = json.dumps([
    {"start_node": 0, "end_node": 1, "E": 200e9, "A": 0.01},
    {"start_node": 0, "end_node": 2, "E": 200e9, "A": 0.01},
    {"start_node": 1, "end_node": 2, "E": 200e9, "A": 0.01},
])
loads = json.dumps([{"node": 2, "fx": 0, "fy": -10000}])

truss = client.truss_analyses.calculate(nodes=nodes, members=members, loads=loads)
print(f"Truss: {truss.num_nodes} nodes, {truss.num_members} members")
for mf in truss.member_forces:
    print(f"  Member {mf.member}: {mf.force:.1f} N ({mf.state})")
for rx in truss.reactions:
    print(f"  Node {rx.node} ({rx.support_type}): Fx={rx.force_x}, Fy={rx.force_y}")

print("exercised: beam_analyses.calculate / mohrs_circles.calculate / cross_sections.calculate / truss_analyses.calculate")