Problem 11.6 - By integration of Moment equation
[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="399"
b=reactive.Value("__")
a=reactive.Value("__")
F=reactive.Value("__")
I=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 ft", 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 beam of lengths a = {a()} ft and b = {b()} ft is loaded with concentrated load F = {F()} kips at the free end. Determine the magnitude of the deflection at the free end if E = 29,000 ksi and I = {I()} in.<sup>4</sup>.")]
@reactive.Effect
@reactive.event(input.generate_problem)
def randomize_vars():
random.seed(input.ID())
b.set(random.randrange(10,50,1)/10)
a.set(round(b())*1.5)
F.set(random.randrange(10,100,1)/10)
I.set(random.randrange(20,100,1))
@reactive.Effect
@reactive.event(input.submit)
def _():
attempt_counter.set(attempt_counter() + 1) # Increment the attempt counter on each submission.
IE = I()*4.82253086*10**-5*29000/0.00694444
x = a()+b()
Ay = -F()*b()/a()
By = F()+Ay
C3 = -Ay*a()**2/6+By*a()**2/2
C4 = -By*a()**3/6
instr= abs(By*x**3/6-By*a()*x**2/2+Ay*x**3/6+C3*x+C4)/IE
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}. For reference in debugging this, the calculated instructor answer is {instr}")
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)