Test Example

from planestress.analysis import PlaneStress
from planestress.pre import LoadCase, Material, circle, rectangle, steel_material

d = 300
b = 250
h_r = 25
t = 16
u = 0.05

steel = steel_material(thickness=t)
rect = rectangle(0.5 * d, 0.5 * b, steel)
hole = circle(2 * h_r, n=32)
geom = rect - hole
geom.plot_geometry()

<Axes: title={'center': 'Geometry'}>

lhs_support = geom.add_line_support((0, h_r), (0, 0.5 * d), "x", 0.0)
bot_support = geom.add_line_support((25, 0), (125, 0), "y", 0.0)
top_disp = geom.add_line_support((0, 0.5 * d), (0.5 * b, 0.5 * d), "y", u)
lc = LoadCase([lhs_support, bot_support, top_disp])

geom.create_mesh(mesh_sizes=10)
geom.plot_mesh()

<Axes: title={'center': 'Finite Element Mesh'}>

ps = PlaneStress(geom, [lc])
results_list = ps.solve()
res = results_list[0]

res.plot_deformed_shape(500)

<Axes: title={'center': 'Deformed Shape [ds = 500]'}>

res.plot_displacement_contour("x", normalize=False, colormap="viridis")
res.plot_displacement_contour("y", normalize=False, colormap="viridis", contours=True)
res.plot_displacement_contour("xy", normalize=False, colormap="viridis")

<Axes: title={'center': 'Displacement Contours [xy]'}>

res.plot_stress("xx")
res.plot_stress("yy", normalize=False, colormap="viridis")
res.plot_stress("xy", contours=True)

<Axes: title={'center': 'Stress Contour Plot - $\\sigma_{xy}$'}>