Problem 8.7 - Area Moment of Inertia
[Problem adapted from © Kurt Gramoll CC BY NC-SA 4.0]
#| standalone: true
#| viewerHeight: 600
#| components: [viewer]
from shiny import App, render, ui, reactive
import random
import asyncio
import io
import math
import string
from datetime import datetime
from pathlib import Path
def generate_random_letters(length):
# Generate a random string of letters of specified length
return ''.join(random.choice(string.ascii_lowercase) for _ in range(length))
problem_ID="347"
b=reactive.Value("__")
c=reactive.Value("__")
d=reactive.Value("__")
a=reactive.Value("__")
d1=reactive.Value("__")
d2=reactive.Value("__")
attempts=["Timestamp,Attempt,Answer,Feedback\n"]
app_ui = ui.page_fluid(
ui.markdown("**Please enter your ID number from your instructor and click to generate your problem**"),
ui.input_text("ID","", placeholder="Enter ID Number Here"),
ui.input_action_button("generate_problem", "Generate Problem", class_="btn-primary"),
ui.markdown("**Problem Statement**"),
ui.output_ui("ui_problem_statement"),
ui.input_text("answer","Your Answer in units of in\u2074", placeholder="Please enter your answer"),
ui.input_action_button("submit", "Submit Answer", class_="btn-primary"),
ui.download_button("download", "Download File to Submit", class_="btn-success"),
)
def server(input, output, session):
# Initialize a counter for attempts
attempt_counter = reactive.Value(0)
@output
@render.ui
def ui_problem_statement():
return[ui.markdown(f"The cross-section shows a concrete beam with two hollow round holes. Determine the area moment of inertia about the beam's centroid. Assume lengths a = {a()} in., b = {b()} in., c = {c()} in., d = {d()} in., d<sub>1</sub> = {d1()} in., and d<sub>2</sub> = {d2()} in.")]
@reactive.Effect
@reactive.event(input.generate_problem)
def randomize_vars():
random.seed(input.ID())
b.set(random.randrange(20, 60, 1)/10)
c.set(round(b()*3,1))
d.set(b()*2)
a.set(b()*2.5)
d1.set(round(b()*1.5,1))
d2.set(d1()*2)
@reactive.Effect
@reactive.event(input.submit)
def _():
attempt_counter.set(attempt_counter() + 1) # Increment the attempt counter on each submission.
A1 = (b()+c()+d())*2*a()
h1 = (b()+c()+d())/2
A2 = (d2()/2)**2*math.pi
h2 = d()
A3 = (d1()/2)**2*math.pi
h3 = d()+c()
h = (A1*h1-A2*h2-A3*h3)/(A1-A2-A3)
y1 = h-h1
y2 = h-h2
y3 = h-h3
I1 = 2*a()*(b()+c()+d())**3/12
I2 = (d2()/2)**4*math.pi/4
I3 = (d1()/2)**4*math.pi/4
instr = I1+A1*y1**2-(I2+A2*y2**2)-(I3+A3*y3**2)
if math.isclose(float(input.answer()), instr, rel_tol=0.01):
check = "*Correct*"
correct_indicator = "JL"
else:
check = "*Not Correct.*"
correct_indicator = "JG"
# Generate random parts for the encoded attempt.
random_start = generate_random_letters(4)
random_middle = generate_random_letters(4)
random_end = generate_random_letters(4)
encoded_attempt = f"{random_start}{problem_ID}-{random_middle}{attempt_counter()}{correct_indicator}-{random_end}{input.ID()}"
# Store the most recent encoded attempt in a reactive value so it persists across submissions
session.encoded_attempt = reactive.Value(encoded_attempt)
# Append the attempt data to the attempts list without the encoded attempt
attempts.append(f"{datetime.now()}, {attempt_counter()}, {input.answer()}, {check}\n")
# Show feedback to the user.
feedback = ui.markdown(f"Your answer of {input.answer()} is {check}.")
m = ui.modal(
feedback,
title="Feedback",
easy_close=True
)
ui.modal_show(m)
@session.download(filename=lambda: f"Problem_Log-{problem_ID}-{input.ID()}.csv")
async def download():
# Start the CSV with the encoded attempt (without label)
final_encoded = session.encoded_attempt() if session.encoded_attempt is not None else "No attempts"
yield f"{final_encoded}\n\n"
# Write the header for the remaining CSV data once
yield "Timestamp,Attempt,Answer,Feedback\n"
# Write the attempts data, ensure that the header from the attempts list is not written again
for attempt in attempts[1:]: # Skip the first element which is the header
await asyncio.sleep(0.25) # This delay may not be necessary; adjust as needed
yield attempt
# App installation
app = App(app_ui, server)