Appendix B: MATLAB Script For Beltrami-Michell Equations: All All

Download as pdf or txt
Download as pdf or txt
You are on page 1of 1

Appendix B: MATLAB Script for Beltrami-Michell Equations

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)

Defining stiffness matrix for a linearly isotropic material.

G = E / (2 * (1 + nu));
epsilon_11 = (1 / E) * (sigma_11 - nu * sigma_22 - nu * sigma_33)

epsilon_11(x_1, x_2, x_3) =

epsilon_22 = (1 / E) * (sigma_22 - nu * sigma_11 - nu * sigma_33)

epsilon_22(x_1, x_2, x_3) =

epsilon_33 = (1 / E) * (sigma_33 - nu * sigma_11 - nu * sigma_22)

epsilon_33(x_1, x_2, x_3) =

gamma_12 = sigma_12 / G

gamma_12(x_1, x_2, x_3) =

gamma_13 = sigma_13 / G

gamma_13(x_1, x_2, x_3) =

gamma_23 = sigma_23

gamma_23(x_1, x_2, x_3) =

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy