Verification, Validation, Uncertainty and Reporting
Comprehensive university textbook unit covering fundamental principles, governing relationships, and engineering applications of Verification, Validation, Uncertainty and Reporting.
§8.1 Code and solution verification
**Target competency.** Derive the appropriate approximate relation, apply it with coherent units, and identify one situation in which it would produce a misleading engineering result. **Primary relation:** `r=Ku-F`. **Field study:** Compare numerical displacement against a known axial-bar solution.
§8.2 Mesh convergence and uncertainty budget
**Target competency.** Derive the appropriate approximate relation, apply it with coherent units, and identify one situation in which it would produce a misleading engineering result. **Primary relation:** `U_total^2≈∑ U_i^2 (independent)`. **Field study:** Compare two thickness alternatives under uncertain loads.
§8.3 Experimental validation planning
**Target competency.** Derive the appropriate approximate relation, apply it with coherent units, and identify one situation in which it would produce a misleading engineering result. **Primary relation:** `e=|Q_FE-Q_test|/|Q_test|`. **Field study:** Design a deflection test for a laboratory cantilever.
§8.4 Safety margin and decision criteria
**Target competency.** Derive the appropriate approximate relation, apply it with coherent units, and identify one situation in which it would produce a misleading engineering result. **Primary relation:** `MS=allowable/demand-1`. **Field study:** Review a bracket thickness reduction with multiple failure modes.
§8.5 Reproducible simulation reports
**Target competency.** Derive the appropriate approximate relation, apply it with coherent units, and identify one situation in which it would produce a misleading engineering result. **Primary relation:** `decision=claims+evidence+assumptions`. **Field study:** Write a one-page model audit for a design review committee.
Step-by-Step Finite Element Analysis and Computational Mechanics Solutions
Three original questions connect calculations, definitions and interpretation. Open each solution to follow the reasoning.
Problem Statement:
Worked Engineering Derivation:
A small hand-checkable model associated with this unit uses **normalized equilibrium residual**, based on `r_norm=abs(R)/max(abs(F),1)`. Illustrative inputs are **Unbalanced nodal force=15 N, Applied force scale=7000 N**. Converting the displayed units into SI inside the calculation gives **0.2142857 %**. This number comes from the stated elementary relation and is *not* a full FE simulation or a measured response. Recreate the substitution independently, confirm the dimensional units, and explain the physical meaning of the output sign and magnitude. The associated browser model provides bounded parameter changes and a sensitivity plot. For negative outcomes such as compressive stress or reverse heat flow, retain the sign in the balance even if a subsequent design comparison uses absolute magnitude. Compute a second scenario by changing the dominant parameter by ten percent while holding all others fixed. If the governing relation is nonlinear in that parameter, do not assume the output changes by ten percent. For example, beam deflection scales with length cubed under Euler–Bernoulli assumptions; this makes a modest length error important. An area or thickness variation can also change local stress and global stiffness differently. When the relevant approximation instead assumes constant geometry, explain why replacing the original problem by a new geometry requires re-analysis rather than simply scaling the final result.
Problem Statement:
Worked Engineering Derivation:
A small hand-checkable model associated with this unit uses **normalized equilibrium residual**, based on `r_norm=abs(R)/max(abs(F),1)`. Illustrative inputs are **Unbalanced nodal force=15 N, Applied force scale=7000 N**. Converting the displayed units into SI inside the calculation gives **0.2142857 %**. This number comes from the stated elementary relation and is *not* a full FE simulation or a measured response. Recreate the substitution independently, confirm the dimensional units, and explain the physical meaning of the output sign and magnitude. The associated browser model provides bounded parameter changes and a sensitivity plot. For negative outcomes such as compressive stress or reverse heat flow, retain the sign in the balance even if a subsequent design comparison uses absolute magnitude. Compute a second scenario by changing the dominant parameter by ten percent while holding all others fixed. If the governing relation is nonlinear in that parameter, do not assume the output changes by ten percent. For example, beam deflection scales with length cubed under Euler–Bernoulli assumptions; this makes a modest length error important. An area or thickness variation can also change local stress and global stiffness differently. When the relevant approximation instead assumes constant geometry, explain why replacing the original problem by a new geometry requires re-analysis rather than simply scaling the final result.
Problem Statement:
Worked Engineering Derivation:
A small hand-checkable model associated with this unit uses **normalized equilibrium residual**, based on `r_norm=abs(R)/max(abs(F),1)`. Illustrative inputs are **Unbalanced nodal force=15 N, Applied force scale=7000 N**. Converting the displayed units into SI inside the calculation gives **0.2142857 %**. This number comes from the stated elementary relation and is *not* a full FE simulation or a measured response. Recreate the substitution independently, confirm the dimensional units, and explain the physical meaning of the output sign and magnitude. The associated browser model provides bounded parameter changes and a sensitivity plot. For negative outcomes such as compressive stress or reverse heat flow, retain the sign in the balance even if a subsequent design comparison uses absolute magnitude. Compute a second scenario by changing the dominant parameter by ten percent while holding all others fixed. If the governing relation is nonlinear in that parameter, do not assume the output changes by ten percent. For example, beam deflection scales with length cubed under Euler–Bernoulli assumptions; this makes a modest length error important. An area or thickness variation can also change local stress and global stiffness differently. When the relevant approximation instead assumes constant geometry, explain why replacing the original problem by a new geometry requires re-analysis rather than simply scaling the final result.
Problem Statement:
Worked Engineering Derivation:
A small hand-checkable model associated with this unit uses **model-to-test relative discrepancy**, based on `e=abs(Q_FE-Q_test)/abs(Q_test)`. Illustrative inputs are **FE prediction=2.45 mm, Measured response=2.38 mm**. Converting the displayed units into SI inside the calculation gives **2.941176 %**. This number comes from the stated elementary relation and is *not* a full FE simulation or a measured response. Recreate the substitution independently, confirm the dimensional units, and explain the physical meaning of the output sign and magnitude. The associated browser model provides bounded parameter changes and a sensitivity plot. For negative outcomes such as compressive stress or reverse heat flow, retain the sign in the balance even if a subsequent design comparison uses absolute magnitude. Compute a second scenario by changing the dominant parameter by ten percent while holding all others fixed. If the governing relation is nonlinear in that parameter, do not assume the output changes by ten percent. For example, beam deflection scales with length cubed under Euler–Bernoulli assumptions; this makes a modest length error important. An area or thickness variation can also change local stress and global stiffness differently. When the relevant approximation instead assumes constant geometry, explain why replacing the original problem by a new geometry requires re-analysis rather than simply scaling the final result.
Problem Statement:
Worked Engineering Derivation:
A small hand-checkable model associated with this unit uses **model-to-test relative discrepancy**, based on `e=abs(Q_FE-Q_test)/abs(Q_test)`. Illustrative inputs are **FE prediction=2.45 mm, Measured response=2.38 mm**. Converting the displayed units into SI inside the calculation gives **2.941176 %**. This number comes from the stated elementary relation and is *not* a full FE simulation or a measured response. Recreate the substitution independently, confirm the dimensional units, and explain the physical meaning of the output sign and magnitude. The associated browser model provides bounded parameter changes and a sensitivity plot. For negative outcomes such as compressive stress or reverse heat flow, retain the sign in the balance even if a subsequent design comparison uses absolute magnitude. Compute a second scenario by changing the dominant parameter by ten percent while holding all others fixed. If the governing relation is nonlinear in that parameter, do not assume the output changes by ten percent. For example, beam deflection scales with length cubed under Euler–Bernoulli assumptions; this makes a modest length error important. An area or thickness variation can also change local stress and global stiffness differently. When the relevant approximation instead assumes constant geometry, explain why replacing the original problem by a new geometry requires re-analysis rather than simply scaling the final result.