Files
jyvikoandClaude Opus 5.5 96c9992c92 oracle: 2D cases (Kirsch plane stress, Lamé plane strain and axisymmetric, rotating disc and cylinder, NAFEMS LE1 with quad8) and J2 uniaxial explicit
The pending 2D LE1 case becomes a real one; the plasticity cases stay
pending with what they still need (static J2 for the collapse pressure,
the transcribed measurements for the Taylor bar).

Co-Authored-By: Claude Opus 5.5 <[email protected]>
Claude-Session: https://claude.ai/code/session_01NmAttHhUZTokgPeEJia8sm
2026-10-13 21:42:19 -04:00

77 lines
1.6 KiB
JSON
Raw Permalink Blame History

This file contains ambiguous Unicode characters
This file contains Unicode characters that might be confused with other characters. If you think that this is intentional, you can safely ignore this warning. Use the Escape button to reveal them.
{
"title": "Thick-walled cylinder under internal pressure, 2D plane strain (Lamé)",
"kind": "theory",
"source": "Lamé's solution, Timoshenko and Goodier, Theory of Elasticity, 3rd ed., McGraw-Hill 1970, §28",
"quantities": [
{
"id": "ur_a",
"quantity": "radial displacement at the bore (10, 0)",
"measure": {
"command": "solve.probe",
"params": {
"step": "Pressure",
"at": [
[
10.0,
0,
0
]
]
},
"pointer": "/nodes/0/displacement/0"
},
"value": 0.00907936507936508,
"tolerance": {
"relative": 0.001
},
"reason": "quadratic elements on 8 radial divisions resolve u_r = A r + B/r to well under 0.1%"
},
{
"id": "hoop_a",
"quantity": "hoop stress σ_θ = σ_yy at the bore (10, 0)",
"measure": {
"command": "solve.probe",
"params": {
"step": "Pressure",
"at": [
[
10.0,
0,
0
]
]
},
"pointer": "/nodes/0/stress/1"
},
"value": 166.66666666666666,
"tolerance": {
"relative": 0.01
},
"reason": "stress extrapolated to the bore surface and averaged over two elements: about 0.5%"
},
{
"id": "szz",
"quantity": "axial stress σ_zz = 2νA at the bore (plane strain)",
"measure": {
"command": "solve.probe",
"params": {
"step": "Pressure",
"at": [
[
10.0,
0,
0
]
]
},
"pointer": "/nodes/0/stress/2"
},
"value": 20.0,
"tolerance": {
"relative": 0.01
},
"reason": "σ_zz = ν(σ_r + σ_θ) follows the in-plane stresses, extrapolated as they are"
}
]
}