Files
jyvikoandClaude Opus 5.5 6d9452ad84 Verification oracle: tests/oracle cases, cargo xtask oracle, fixtures, docs/oracle.md, CI job
The oracle binary (modelcraft-testkit) runs each case's recipe through the
command engine, measures its quantities by JSON pointers into step outputs
or commands, compares them with the cited references and the frozen
fixtures, and writes docs/oracle.md. 20 theory cases (patch tests of every
solid element, Timoshenko cantilevers, Lame, Kirsch, beam frequencies,
torsion and bending through couplings, explicit waves and impact, CDT
invariants), 8 published benchmarks (MacNeal-Harder patch, straight and
curved cantilevers, thick cylinder, Scordelis-Lo with solids, NAFEMS LE1,
LE10, FV32) and 10 pending cases with citations.

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.9 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": "MacNeal–Harder thick-walled cylinder, nearly incompressible (hex8 with incompatible modes)",
"kind": "benchmark",
"source": "MacNeal and Harder, 'A proposed standard set of problems to test finite element accuracy', Finite Elements in Analysis and Design 1 (1985) 3–20: Fig. 10 (p. 9) geometry, material and radial mesh; Table 5(c) (p. 11) radial displacements at the inner radius",
"quantities": [
{
"id": "ur_1",
"quantity": "radial displacement at the inner radius, ν = 0.49",
"measure": {
"command": "solve.probe",
"params": {
"step": "Pressure",
"at": [
[
3.0,
0,
0.0
]
]
},
"pointer": "/nodes/0/displacement/0"
},
"value": 0.0050399,
"tolerance": {
"relative": 0.02
},
"reason": "MacNeal and Harder grade an element A within 2% of the theoretical value on this prescribed coarse mesh (p. 14)"
},
{
"id": "ur_2",
"quantity": "radial displacement at the inner radius, ν = 0.499",
"measure": {
"command": "solve.probe",
"params": {
"step": "Pressure",
"at": [
[
3.0,
0,
2.0
]
]
},
"pointer": "/nodes/0/displacement/0"
},
"value": 0.0050602,
"tolerance": {
"relative": 0.02
},
"reason": "MacNeal and Harder grade an element A within 2% of the theoretical value on this prescribed coarse mesh (p. 14)"
},
{
"id": "ur_3",
"quantity": "radial displacement at the inner radius, ν = 0.4999",
"measure": {
"command": "solve.probe",
"params": {
"step": "Pressure",
"at": [
[
3.0,
0,
4.0
]
]
},
"pointer": "/nodes/0/displacement/0"
},
"value": 0.0050623,
"tolerance": {
"relative": 0.02
},
"reason": "MacNeal and Harder grade an element A within 2% of the theoretical value on this prescribed coarse mesh (p. 14)"
}
]
}