API tools faq
paste
Advertisement
Guest User

SenasteProgrammeeettt

a guest
Nov 28th, 2014
177
0
Never
Add comment
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 12.22 KB | None | 0 0
raw download clone embed print report
  1. %EquilibriumSolverV0.1
  2. clear all, close all, clc
  3. %% Värden
  4. v=70;
  5. R=20;
  6. m=150;
  7. A=0.6;
  8. df=800;
  9. db=200;
  10. h=20;
  11. h1=30;
  12. b=1000;
  13. b1=100;
  14. bd=100;
  15. bb=100;
  16. rd=80;
  17. rb=80;
  18. a1=5;
  19. a2=-8;
  20. c=0.25;
  21. dh=120;
  22. ns=1;
  23. sigmas=(200e6);
  24. g=9.81;
  25. rho=1.21;
  26. v=v*1000/3600;
  27. df=df*1e-3;
  28. db=db*1e-3;
  29. h=h*1e-2;
  30. h1=h1*1e-2;
  31. b=b*1e-3;
  32. b1=b1*1e-3;
  33. bd=bd*1e-3;
  34. bb=bb*1e-3;
  35. rd=rd*1e-3;
  36. rb=rb*1e-3;
  37. dh=dh*1e-3;
  38. d1=b1;
  39. d2=b/2+bd;
  40. d3=b-b1;
  41. d4=b;
  42.  
  43. %% Jämvikt Fall D
  44. FL=0.5*rho*c*A*(0.7*v)^2;
  45. FD=FL;
  46. Nb=(FD*(h+h1)+m*g*df)/(db+df);
  47. Nf=m*g-Nb;
  48. Vi=m*g*0.5-(h*(m*(0.7*v)^2)/(R*b));
  49. Vy=m*g*0.5+(h*(m*(0.7*v)^2)/(R*b));
  50. Hb=(m*(0.7*v)^2)/(R*(df+db))*df;
  51. Vbi=Vi*(1/(1+(Nf/Nb)));
  52. Vby=Vy*(1/(1+(Nf/Nb)));
  53. Hby=Hb*(1/(1+(Vbi/Vby)));
  54. Hbi=Hb-Hby;
  55. R1x=Hbi+Hby;
  56. Fk=FD*dh*0.5*(1/rd);
  57. R1y=-Fk*(b/2-b1-bd)/(b-2*b1)+(FD*(2*b1-b))/(2*(b-2*b1));
  58. R2y=-Fk*(b/2-b1+bd)/(b-2*b1)-(FD*(b-b1))/(2*(b-2*b1))+(FD*b1*0.5)/(b-2*b1);
  59. R2z=(Vbi*b1+Hbi*0.5*dh+Hby*0.5*dh-Vby*(b-b1))/(b-2*b1);
  60. R1z=(Vby*b1-Vbi*(b-b1)-(Hbi+Hby)*0.5*dh)/(b-2*b1);
  61. %% Snittkrafter
  62. MHbi=Hbi*dh/2;
  63. MHby=Hby*dh/2;
  64.  
  65. X1y=0:0.001:b1;
  66. X2y=b1+0.001:0.001:(0.5*b+bd);
  67. X3y=(0.5*b+bd+0.001):0.001:(b-b1);
  68. X4y=(b-b1+0.001):0.001:b;
  69. Xiy=[X1y X2y X3y X4y];
  70. %http://pastebin.com/5RgCS36F
  71. %% Fall D
  72.  
  73. %Z-led
  74. Tz1D=-Vbi*[1:length(X1y)].^0; My1D=(FD/2)*X1y; %från 0 till <b1
  75. Tz2D=(-Vbi-R1z)*[1:length(X2y)].^0; My2D=(FD/2)*X2y+R1y*(X2y-b1);  %från b1 till <0.5*b+bd
  76. Tz3D=(-Vbi-R1z)*[1:length(X3y)].^0; My3D=(FD/2)*X3y+R1y*(X3y-b1)+Fk*(X3y-0.5*b-bd); %Från 0.5*b+bd till <b-b1
  77. Tz4D=(-Vbi-R1z-R2z)*[1:length(X4y)].^0; My4D=0.5*FD*X4y+R1y*(X4y-b1)+Fk*(X4y-0.5*b-bd)+R2y*(X4y-b+b1); %från b-b1 till <b
  78.  
  79. TzD=[Tz1D Tz2D Tz3D Tz4D]; %dessa ska vi plotta mha av
  80. MyD=[My1D My2D My3D My4D]; %same goes for this badboy
  81.  
  82.  
  83. %y-led
  84. Ty1D=-FD/2*([1:length(X1y)].^0); Mz1D=-MHbi-Vbi*X1y; %0 till <b1
  85. Ty2D=((-FD/2)-R1y)*([1:length(X2y)].^0);  Mz2D=-MHbi-Vbi*X2y-R1z*(X2y-b1); %b1 till <b-b1
  86. Ty3D=((-FD/2)-R1y-Fk)*([1:length(X3y)].^0); Mz3D=-MHbi-Vbi*X3y-R1z*(X3y-b1);%Momentberäkning bör göras om %b-b1 till <L
  87. Ty4D=(-FD*0.5-R1y-Fk-R2y)*([1:length(X4y)].^0);  Mz4D=-MHbi-Vbi*X4y-R1z*(X4y-b1)-R2z*(X4y-b+b1);%Educational purpose only
  88.  
  89. TyD=[Ty1D Ty2D Ty3D Ty4D];
  90. MzD=[Mz1D Mz2D Mz3D Mz4D];
  91.  
  92. MtotD=sqrt(MyD.^2+MzD.^2);
  93. %%
  94.  
  95. % Fall C
  96. %% Jämvikt Fall C
  97. vc=0;
  98. FLC=0.5*rho*c*A*(vc)^2;
  99. FDC=FLC+m*a1;
  100. NbC=(FDC*(h+h1)+m*(g*df-a1*h1))/(db+df);
  101. NfC=m*g-NbC;
  102. FkC=(FDC*dh)/(2*rd);
  103. R2yC=(FkC*(d2-d1)+FDC*(d4/2-d1))/(d3-d1);
  104. R1yC=FkC+FDC-R2yC;
  105. R1zC=NbC*(1-(d4/2-d1)/(d3-d1));
  106. R2zC=R1zC;
  107.  
  108.  
  109.  
  110. %z-led
  111. Tz1C=(-FDC/2)*([1:length(X1y)].^0); Mz1C=(FDC/2)*X1y;
  112. Tz2C= (R1yC-FDC/2)*([1:length(X2y)].^0); Mz2C=(R1yC*d1)-(R1yC-FDC/2)*X2y;
  113. Tz3C=(R1yC-FkC-FDC/2)*([1:length(X3y)].^0); Mz3C=R1yC*d1-FkC*d2-(R1yC-FkC-FDC/2)*X3y;
  114. Tz4C=(R1yC+R2yC-FkC-FDC/2)*([1:length(X4y)].^0); Mz4C=R1yC*d1+R2yC*d3-FkC*d2-(R1yC+R2yC-FkC-FDC/2)*X4y;
  115.  
  116. %y-led
  117. Ty1C=(-NbC/2)*([1:length(X1y)].^0); My1C=(-NbC/2)*X1y;
  118. Ty2C=(R1zC-NbC/2)*([1:length(X2y)].^0); My2C=(R1zC-NbC/2)*X2y-R1zC*d1;
  119. Ty3C=(R1zC-NbC/2)*([1:length(X3y)].^0); My3C=(R1zC-NbC/2)*X3y-R1zC*d1;
  120. Ty4C=(R1zC+R2zC-NbC/2)*([1:length(X4y)].^0); My4C=(R1zC+R2zC-NbC/2)*X4y-R1zC*d1-R2zC*d3;
  121.  
  122. TzC=[Tz1C Tz2C Tz3C Tz4C];
  123. MzC=[Mz1C Mz2C Mz3C Mz4C];
  124.  
  125. TyC=[Ty1C Ty2C Ty3C Ty4C];
  126. MyC=[My1C My2C My3C My4C];
  127.  
  128. MtotC=sqrt(MyC.^2+MzC.^2);
  129.  
  130.  
  131. %%
  132. %% FALL B
  133. a1=0;
  134.  
  135. FLB=0.5*rho*c*A*(v)^2;
  136. FDB=FLB+m*a1;
  137. NbB=(FDB*(h+h1)+m*(g*df-a1*h1))/(db+df);
  138. NfB=m*g-NbB;
  139. FkB=(FDB*dh)/(2*rd);
  140. R2yB=(FkB*(d2-d1)+FDB*((d4/2)-d1))/(d3-d1);
  141. R1yB=FkB+FDB-R2yB;
  142. R1zB=NbB*(1-((d4/2)-d1)/(d3-d1));
  143. R2zB=R1zB;
  144.  
  145. %z-led
  146. Tz1B=(-FDB/2)*([1:length(X1y)].^0); Mz1B=(FDB/2)*X1y;
  147. Tz2B= (R1yB-FDB/2)*([1:length(X2y)].^0); Mz2B=(R1yB*d1)-(R1yB-FDB/2)*X2y;
  148. Tz3B=(R1yB-FkB-FDB/2)*([1:length(X3y)].^0); Mz3B=R1yB*d1-FkB*d2-(R1yB-FkB-FDB/2)*X3y;
  149. Tz4B=(R1yB+R2yB-FkB-FDB/2)*([1:length(X4y)].^0); Mz4B=R1yB*d1+R2yB*d3-FkB*d2-(R1yB+R2yB-FkB-FDB/2)*X4y;
  150.  
  151. %y-led
  152. Ty1B=(-NbB/2)*([1:length(X1y)].^0); My1B=(-NbB/2)*X1y;
  153. Ty2B=(R1zB-NbB/2)*([1:length(X2y)].^0); My2B=(R1zB-NbB/2)*X2y-R1zB*d1;
  154. Ty3B=(R1zB-NbB/2)*([1:length(X3y)].^0); My3B=(R1zB-NbB/2)*X3y-R1zB*d1;
  155. Ty4B=(R1zB+R2zB-NbB/2)*([1:length(X4y)].^0); My4B=(R1zB+R2zB-NbB/2)*X4y-R1zB*d1-R2zB*d3;
  156.  
  157. TzB=[Tz1B Tz2B Tz3B Tz4B];
  158. MzB=[Mz1B Mz2B Mz3B Mz4B];
  159.  
  160. TyB=[Ty1B Ty2B Ty3B Ty4B];
  161. MyB=[My1B My2B My3B My4B];
  162.  
  163. MtotB=sqrt(MyB.^2+MzB.^2);
  164.  
  165.  
  166.  
  167. % Fall A
  168. %% Jämvikt Fall A
  169. v=0;
  170. a1=0;
  171. FLA=0.5*rho*c*A*(v)^2;
  172. FDA=FLA+m*a1;
  173. NbA=(FDA*(h+h1)+m*(g*df-a1*h1))/(db+df);
  174. NfA=m*g-NbA;
  175. FkA=(FDA*dh)/(2*rd);
  176. R2yA=(FkA*(d2-d1)+FDA*(d4/2-d1))/(d3-d1);
  177. R1yA=FkA+FDA-R2yA;
  178. R1zA=NbA*(1-(d4/2-d1)/(d3-d1));
  179. R2zA=R1zA;
  180.  
  181. %z-led
  182. Tz1A=(-FDA/2)*([1:length(X1y)].^0); Mz1A=(FDA/2)*X1y;
  183. Tz2A= (R1yA-FDA/2)*([1:length(X2y)].^0); Mz2A=(R1yA*d1)-(R1yA-FDA/2)*X2y;
  184. Tz3A=(R1yA-FkA-FDA/2)*([1:length(X3y)].^0); Mz3A=R1yA*d1-FkA*d2-(R1yA-FkA-FDA/2)*X3y;
  185. Tz4A=(R1yA+R2yA-FkA-FDA/2)*([1:length(X4y)].^0); Mz4A=R1yA*d1+R2yA*d3-FkA*d2-(R1yA+R2yA-FkA-FDA/2)*X4y;
  186.  
  187. %y-led
  188. Ty1A=(-NbA/2)*([1:length(X1y)].^0); My1A=(-NbA/2)*X1y;
  189. Ty2A=(R1zA-NbA/2)*([1:length(X2y)].^0); My2A=(R1zA-NbA/2)*X2y-R1zA*d1;
  190. Ty3A=(R1zA-NbA/2)*([1:length(X3y)].^0); My3A=(R1zA-NbA/2)*X3y-R1zA*d1;
  191. Ty4A=(R1zA+R2zA-NbA/2)*([1:length(X4y)].^0); My4A=(R1zA+R2zA-NbA/2)*X4y-R1zA*d1-R2zA*d3;
  192.  
  193. TzA=[Tz1A Tz2A Tz3A Tz4A];
  194. MzA=[Mz1A Mz2A Mz3A Mz4A];
  195.  
  196. TyA=[Ty1A Ty2A Ty3A Ty4A];
  197. MyA=[My1A My2A My3A My4A];
  198.  
  199. MtotA=sqrt(MyA.^2+MzA.^2);
  200.  
  201.  
  202. % %% Figure 2 PLOTS
  203. % h2=figure('name','Mtot');
  204. % plot(Xiy,MtotA,'b')
  205. % hold on
  206. % plot(Xiy,MtotB,'r')
  207. % hold on
  208. % plot(Xiy,MtotC,'g')
  209. % hold on
  210. % plot(Xiy,MtotD,'k')
  211. %
  212. % legend('Total moment Fall A', 'Total moment Fall B', 'Total moment Fall C', 'Total moment Fall D')
  213. % %% Figure 1 %PLOTS
  214. % h=figure('name','T och M-plots');
  215. % %Plot 1
  216. % subplot(2,2,1)
  217. % plot(Xiy,TzA,'b')
  218. % hold on
  219. % plot(Xiy,TzB,'r')
  220. % hold on
  221. % plot(Xiy,TzC,'g')
  222. % hold on
  223. % plot(Xiy,TyD,'k')
  224. % legend('Tvärkraft, z-rikt. Fall A', 'Tvärkraft, z-rikt. Fall B', 'Tvärkraft, z-rikt. Fall C', 'Tvärkraft, z-rikt. Fall D')
  225. % %--------------------
  226. % %Plot 2
  227. % subplot(2,2,2)
  228. % plot(Xiy,TyA,'b')
  229. % hold on
  230. % plot(Xiy,TyB,'r')
  231. % hold on
  232. % plot(Xiy,TyC,'g')
  233. % hold on
  234. % plot(Xiy,TzD,'k')
  235. % legend('Tvärkraft, y-rikt. Fall A', 'Tvärkraft, y-rikt. Fall B', 'Tvärkraft, y-rikt. Fall C', 'Tvärkraft, y-rikt. Fall D')
  236. % %--------------------
  237. % %Plot 3
  238. % subplot(2,2,3)
  239. % plot(Xiy,MyA,'b')
  240. % hold on
  241. % plot(Xiy,MyB,'r')
  242. % hold on
  243. % plot(Xiy,MyC,'g')
  244. % hold on
  245. % plot(Xiy,MzD,'k')
  246. % legend('Moment, y-rikt. Fall A', 'Moment, y-rikt. Fall B', 'Moment, y-rikt. Fall C', 'Moment, y-rikt. Fall D')
  247. % %--------------------
  248. % %Plot 4
  249. % subplot(2,2,4)
  250. % plot(Xiy,MzA,'b')
  251. % hold on
  252. % plot(Xiy,MzB,'r')
  253. % hold on
  254. % plot(Xiy,MzC,'g')
  255. % hold on
  256. % plot(Xiy,MyD,'k')
  257. % legend('Moment, z-rikt. Fall A', 'Moment, z-rikt. Fall B', 'Moment, z-rikt. Fall C', 'Moment, z-rikt. Fall D')
  258.  
  259. %--------------------
  260. %% Vridmoment och normalkrafter
  261.  
  262. %FALL D
  263. N1D=Hbi*([1:length(X1y)].^0);
  264. N2D=(Hbi-R1x)*([1:length(X2y)].^0);
  265. N3D=(Hbi-R1x)*([1:length(X3y)].^0);
  266. N4D=(Hby+Hbi-R1x)*([1:length(X4y)].^0);
  267. ND=[N1D N2D N3D N4D];
  268. NtotD=[N1D N2D N3D N4D];
  269. NmaxD=max(NtotD);
  270. Mv1D=((FD*dh)/4)*X1y;
  271. Mv2D=((FD*dh)/4)*X2y;
  272. Mv3D=(-(Fk*rd)+((FD*dh)/4))*X3y;
  273. Mv4D=-(Fk*rd)*X4y;
  274. MvTotD=[Mv1D Mv2D Mv3D Mv4D];
  275. TmaxD=sqrt(TzD.^2+TyD.^2);
  276.  
  277. %FALL C
  278. Mv1C=((FDC*dh)/4)*([1:length(X1y)].^0);
  279. Mv2C=((FDC*dh)/4)*([1:length(X2y)].^0);
  280. Mv3C=(((FDC*dh)/4)+(FkC*rd))*([1:length(X3y)].^0);
  281. Mv4C=(((FDC*dh)/2)-(FkC*rd))*([1:length(X4y)].^0);
  282. MvTotC=[Mv1C Mv2C Mv3C Mv4C];
  283. TmaxC=sqrt(TzC.^2+TyC.^2);
  284. NC=0*([1:length(Xiy)].^0);
  285.  
  286. %FALL B
  287. Mv1B=((FDB*dh)/4)*([1:length(X1y)].^0);
  288. Mv2B=((FDB*dh)/4)*([1:length(X2y)].^0);
  289. Mv3B=(((FDB*dh)/4)+(FkB*rd))*([1:length(X3y)].^0);
  290. Mv4B=(FkB*rd)*([1:length(X4y)].^0);
  291. MvTotB=[Mv1B Mv2B Mv3B Mv4B];
  292. TmaxB=sqrt(TzB.^2+TyB.^2);
  293. NB=0*([1:length(Xiy)].^0);
  294.  
  295. %FALL A
  296. Mv1A=((FDA*dh)/4)*([1:length(X1y)].^0);
  297. Mv2A=((FDA*dh)/4)*([1:length(X2y)].^0);
  298. Mv3A=(((FDA*dh)/4)+(FkA*rd))*([1:length(X3y)].^0);
  299. Mv4A=(FkA*rd)*([1:length(X4y)].^0);
  300. MvTotA=[Mv1A Mv2A Mv3A Mv4A];
  301. TmaxA=sqrt(TzA.^2+TyA.^2);
  302. NA=0*([1:length(Xiy)].^0);
  303.  
  304. %skit i detta tillsvidare
  305.  
  306. %Dimensionering
  307.  
  308. axelnsDelar=[0.6 1 1 0.6];
  309. diameterD=[];
  310. slutDiameterD=[];
  311. diameterC=[];
  312. slutDiameterC=[];
  313. diameterB=[];
  314. slutDiameterB=[];
  315. diameterA=[];
  316. slutDiameterA=[];
  317.  
  318.  
  319.  
  320. ghettoMatris=[101 500 300 100];
  321.  
  322. %ta fram max värden för de olika grejerna
  323. counterMin=0;
  324. counterMax=0;
  325. for i = 1:4
  326.     rotterD=[];
  327.     rotterC=[];
  328.     rotterB=[];
  329.     rotterA=[];
  330.     counterMin=counterMax+1;
  331.     counterMax=counterMax+ghettoMatris(i);
  332.    
  333.     for k = counterMin:counterMax %axelns alla delar (aka x axeln)
  334.         %Fall D
  335.         rD=roots([1 0 0 0 -1/(pi^2*sigmas^2)*(16*ND(k).^2+48*TmaxD(k).^2)/(axelnsDelar(i)^4) -(256*MtotD(k).*ND(k)+384*MvTotD(k).*TmaxD(k))/(pi^2*sigmas^2)/(axelnsDelar(i)^5) -(1024*MtotD(k).^2+768*MvTotD(k).^2)/(pi^2*sigmas^2)/(axelnsDelar(i)^6)]);
  336.         rD=rD(imag(rD)==0);
  337.         rotterD=[rotterD rD'];
  338.        
  339.         %Fall C
  340.         rC=roots([1 0 0 0 -1/(pi^2*sigmas^2)*(16*NC(k).^2+48*TmaxC(k).^2)/(axelnsDelar(i)^4) -(256*MtotC(k).*NC(k)+384*MvTotC(k).*TmaxC(k))/(pi^2*sigmas^2)/(axelnsDelar(i)^5) -(1024*MtotC(k).^2+768*MvTotC(k).^2)/(pi^2*sigmas^2)/(axelnsDelar(i)^6)]);
  341.         rC=rC(imag(rC)==0);
  342.         rotterC=[rotterC rC'];
  343.        
  344.         %Fall B
  345.         rB=roots([1 0 0 0 -1/(pi^2*sigmas^2)*(16*NB(k).^2+48*TmaxB(k).^2)/(axelnsDelar(i)^4) -(256*MtotB(k).*NB(k)+384*MvTotB(k).*TmaxB(k))/(pi^2*sigmas^2)/(axelnsDelar(i)^5) -(1024*MtotB(k).^2+768*MvTotB(k).^2)/(pi^2*sigmas^2)/(axelnsDelar(i)^6)]);
  346.         rB=rB(imag(rB)==0);
  347.         rotterB=[rotterB rB'];
  348.        
  349.         %Fall A
  350.         rA=roots([1 0 0 0 -1/(pi^2*sigmas^2)*(16*NA(k).^2+48*TmaxA(k).^2)/(axelnsDelar(i)^4) -(256*MtotA(k).*NA(k)+384*MvTotA(k).*TmaxA(k))/(pi^2*sigmas^2)/(axelnsDelar(i)^5) -(1024*MtotA(k).^2+768*MvTotA(k).^2)/(pi^2*sigmas^2)/(axelnsDelar(i)^6)]);
  351.         rA=rA(imag(rA)==0);
  352.         rotterA=[rotterA rA'];
  353.     end
  354.     %Fall D
  355.     diameterD=max(abs(rotterD));
  356.     diameterD=max(diameterD);
  357.     slutDiameterD=[slutDiameterD diameterD];
  358.    
  359.     %Fall C
  360.     diameterC=max(abs(rotterC));
  361.     diameterC=max(diameterC);
  362.     slutDiameterC=[slutDiameterC diameterC];
  363.    
  364.      %Fall B
  365.     diameterB=max(abs(rotterB));
  366.     diameterB=max(diameterB);
  367.     slutDiameterB=[slutDiameterB diameterB];
  368.    
  369.      %Fall C
  370.     diameterA=max(abs(rotterA));
  371.     diameterA=max(diameterA);
  372.     slutDiameterA=[slutDiameterA diameterA];
  373. end
  374. %--------------------
  375. %DiameterMax
  376. dMax=max(slutDiameterD);
  377. cMax=max(slutDiameterC);
  378. bMax=max(slutDiameterB);
  379. aMax=max(slutDiameterA);
  380.  
  381. %NormalkraftsMax
  382. N_MaxD=max(ND);
  383. N_MaxC=max(NC);
  384. N_MaxB=max(NB);
  385. N_MaxA=max(NA);
  386.  
  387. %MaxMoment
  388. MMaxD=max(MtotD);
  389.  
  390. %Klar med sånt skit, dags för plots
  391.  
  392. %% Von mises plots stuff
  393.  
  394. tauMaxD=((MvTotD*16)/(pi*dMax^3)+((4*TmaxD)/(pi*dMax^2)));
  395. tauMaxC=((MvTotC*16)/(pi*cMax^3)+((4*TmaxC)/(pi*cMax^2)));
  396. tauMaxB=((MvTotB*16)/(pi*bMax^3)+((4*TmaxB)/(pi*bMax^2)));
  397. tauMaxA=((MvTotA*16)/(pi*aMax^3)+((4*TmaxA)/(pi*aMax^2)));
  398. h3=figure('name','TauMax-plots');
  399. plot(Xiy, tauMaxD, 'b')
  400. hold on
  401. plot(Xiy, tauMaxC, 'r')
  402. hold on
  403. plot(Xiy, tauMaxB, 'g')
  404. hold on
  405. plot(Xiy, tauMaxA, 'k')
  406. legend('Max Tau. Fall D', 'Max Tau. Fall C', 'Max Tau. Fall B', 'Max Tau. Fall A')
  407.  
  408.  
  409. NtotC=0;
  410. NtotB=0;
  411. NtotA=0;
  412. h4=figure('name','SigmaMax-plots');
  413. sigmaDMax = ((4*NtotD)/(pi*dMax^2) + (32*MtotD)/(pi*dMax^3));
  414. sigmaCMax = ((4*NtotC)/(pi*cMax^2) + (32*MtotC)/(pi*cMax^3));
  415. sigmaBMax = ((4*NtotB)/(pi*bMax^2) + (32*MtotB)/(pi*bMax^3));
  416. sigmaAMax = ((4*NtotA)/(pi*aMax^2) + (32*MtotA)/(pi*aMax^3));
  417. plot(Xiy, sigmaDMax, 'b');
  418. hold on
  419. plot(Xiy, sigmaCMax, 'r');
  420. hold on
  421. plot(Xiy, sigmaBMax, 'g');
  422. hold on
  423. plot(Xiy, sigmaAMax, 'k');
  424. legend('Max Sigma. Fall D', 'Max Sigma. Fall C', 'Max Sigma. Fall B', 'Max Sigma. Fall A')
  425.  
  426. h5=figure('name', 'VonMises-plot');
  427.  
  428. sigmaMisesD=sqrt(sigmaDMax.^2+3*tauMaxD.^2);
  429. sigmaMisesC=sqrt(sigmaCMax.^2+3*tauMaxC.^2);
  430. sigmaMisesB=sqrt(sigmaBMax.^2+3*tauMaxB.^2);
  431. sigmaMisesA=sqrt(sigmaAMax.^2+3*tauMaxA.^2);
  432.  
  433. plot(Xiy, sigmaMisesD, 'b');
  434. hold on
  435. plot(Xiy, sigmaMisesC, 'r');
  436. hold on
  437. plot(Xiy, sigmaMisesB, 'g');
  438. hold on
  439. plot(Xiy, sigmaMisesA, 'k');
  440.  
  441. legend('Von Mises. Fall D', 'Von Mises. Fall C', 'Von Mises. Fall B', 'Von Mises. Fall A')
  442.  
  443. %Skriv ut övergångs spänningen för alla falll.....
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement
Public Pastes
Advertisement
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy.  OK, I Understand
Not a member of Pastebin yet?
Sign Up, it unlocks many cool features!