Info

This question is closed. Reopen it to edit or answer.

Need help writing m file?

2 views (last 30 days)
Steven P.
Steven P. on 1 May 2021
Closed: John D'Errico on 1 May 2021
Hello,
I need assistance (Any I can get) with writing an m file for the following problem:
I have also attached some of the formulas and steps we were given to attempt solving this. I am not very familiar with writing code. So, the engineering properties are used to determine the S matrix which for orthotropic under plane stress is equal to: [1/E1 -v12/E1 0; -v12/E1 1/E2 0; 0 0 1/G12]. Q matrix is determined by the inverse of S : Q=inv(S). This however is in the local coordinate system (where theta or the angle is equal to 0). We can transform the q matrix of each ply according to the formula shown below for Q bar in the attached pictures. Q bar represents the global stiffness matrix in each ply which is transformed via its orientation using the below formulas where Qbar=Tsigma^-1 * Q* Tepsilon of each ply from the neutral axis (so n*thickness of each ply divided by 2 would represent the neutral axis). This is shown below as well. we need to define an array of h values which will represent the distance from the neutral axis to each interface. I am struggling the most here in writing code for this. Then, we need to determine the Q bar for each ply, which I am also having trouble with. Then, we are supposed to define the A, B and D matrices as a function of the h values we calculated (or their location) and the local Qbar values for each ply. Once that is done, we need to define a force N and moment M which we input and need to caluclate the global strains for each ply by the below formula of {N; M}=[A B; B D]{; K} where epsilon represents the global strain for each ply. Then we need to determine this by mulitplying the global strains in each ply times the Qbar in each ply which will give globals stresses. Then he wants us to plot the global stresses and strains.
I know that is a lot of infromation and can be quite confusing. Any help or direction would be greatly appreciated. Any comments in general would also be helpful. If anything needs to be clarified I can do that as well. I would usually not ask for this much help but we never really taught how to use matlab in class so I am not sure where to go from here.

Answers (0)

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!