Problem 6.29 - Statically Indeterminate Torsion
[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 = "289"
T = reactive.Value("__")
L = reactive.Value("__")
ro = reactive.Value("__")
ri = reactive.Value("__")
GA = 3800000
GS = 11000000
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 ksi", 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"A composite circular rod is made from aluminum (G = 3,800 ksi) and steel (G = 11,000 ksi) as shown. A torque T = {T()} lb-ft is applied to the free end. What is the maximum stress in either material? Assume length L = {L()} in., outer radius r<sub>o</sub> = {ro()} in., and inner radius r<sub>i</sub> = {ri()} in."
)
]
@reactive.Effect
@reactive.event(input.generate_problem)
def randomize_vars():
random.seed(input.ID())
T.set(random.randrange(300, 2000, 100))
L.set(random.randrange(10, 60, 1))
ro.set(random.randrange(10, 60, 1) / 10)
ri.set(round(ro() / 1.5, 1))
@reactive.Effect
@reactive.event(input.submit)
def _():
attempt_counter.set(
attempt_counter() + 1
) # Increment the attempt counter on each submission.
JS = math.pi * ri() ** 4 / 2
JA = math.pi / 2 * (ro() ** 4 - ri() ** 4)
TA = T() * 12 * GA * JA / (GS * JS + GA * JA)
TS = T() * 12 - TA
instr = TS * ri() / JS / 1000
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)