Appendix B: MATLAB Script For Beltrami-Michell Equations: All All
Appendix B: MATLAB Script For Beltrami-Michell Equations: All All
Appendix B: MATLAB Script For Beltrami-Michell Equations: All All
clear all
close all
clc
Defining symbolic variables for displacement field, body forces and material properties.
syms sigma_11(x_1, x_2, x_3) sigma_22(x_1, x_2, x_3) sigma_33(x_1, x_2, x_3)
syms sigma_12(x_1, x_2, x_3) sigma_13(x_1, x_2, x_3) sigma_23(x_1, x_2, x_3)
syms E G nu
assume(E > 0)
assume(nu > 0 & nu~= 0.5)
G = E / (2 * (1 + nu));
epsilon_11 = (1 / E) * (sigma_11 - nu * sigma_22 - nu * sigma_33)
gamma_12 = sigma_12 / G
gamma_13 = sigma_13 / G
gamma_23 = sigma_23