How to Plot Stress Contour??

I am a beginner , My doubt is that I have to Plot Stress Contour of tapered cantilever beam analysed using FEA now I have the value of element Stresses like sigmax , sigmay , tauxy ? How to proceed ?? Please Help..

[darova]

Can you show/attach the data?

[Shubhang Jain]

Thankyou for responding first .

I have divided the beam in 10 elements in plane X-Y .Where X=(0:0.1:2) and Y=(0:0.025:0.5) . I have done programming to formulate Elemental stresses i.e 30*1 matrix having 10 sigmax, sigmay and tauxy each .

 Sigma_X=1.0e+03 * [-4.2996; -3.6313; -2.7780; -1.8074; -0.6545; 4.2996; 3.6313; 2.7780; 1.8074; 0.6545],
   Sigma_y =[ -321.9020; -20.6500;  -136.9501;  -65.3855;  -151.7759;  235.7197;  -103.1941;  -126.590; -42.9998;  -418.4672],
    Tau_xy =[ -795.8185;  -899.8615;  -878.6821;  -706.5482; -693.8920;  -543.8945;  -383.5610;  -575.8634; -412.3329;  -628.4221],   

Now I have to plot contour of Sigma_X , Sigma_y , Tau_xy .

I am attaching an image of beam so that you can understand the situation .

Please help me with the solution.

Thanks in Advance...

[darova]

DO you have deformations or stresses in each point?

[Shubhang Jain]

No not at nodal points , I have used matlab functions to determine elemental stiffness then assembled global stiffness to determine elemental stresses . Hence I have stress value at each element overall. If you have any idea to convert element stress value to nodal points then please suggest.

[darova]

Can you show what elements are?

Did you try patch?

[Shubhang Jain]

I have found out the stresses sigmax, sigmay, tauxy values at each nodes point.

Sigma_X = 1.0e+03 * [-4299.57587958638 ;	-3965.44553517593;   -3204.63754334495;	-2292.69609509300	-1230.97384069952;	-654.515387137445;	7.73070496506989e-12;  8.64019966684282e-12;	1.11413100967184e-11;	  1.02318153949454e-11;	2.84217094304040e-12; 	-1.93267624126747e-12;	4299.57587958640;	3965.44553517595;	3204.63754334497;	2292.69609509302;	1230.97384069952;	654.515387137441]

Sigma_Y = [-321.901962451787;	-171.276004860075;	-78.8000793342068;	-101.167806049306;	-108.580704946648;	-151.775909194735;	-43.0911561819484;	-52.5066078518920;	-96.8460482931672;	-92.9813362559606;	-169.657107326625;	-285.121579205828;	235.719650087890;	66.2627891562907;	-114.892017252128;	-84.7948664626148;	-230.733509706602;	-418.467249216922]

Tau_XY =[ -795.818460340878;	-847.839977192489;	-889.271783271534;	-792.615142694200;	-700.220084886960;	-693.891956884487;	-669.856459330285;	-655.783844638485;	-684.491978609799;	-643.356643356871;	-610.298792117163;	-661.157024793497;	-543.894458319692;	-463.727712084480;	-479.712173948064;	-494.098144019542;	-520.377499347365;	-628.422092702508] 

Can you please help me to find the stress contour , My submission is on 3 JULY !!

Thanks in advance.

[darova]

did you try surf?

[Shubhang Jain]

Sorry but I am not familiar with surf . I am coding for first time in Matlab . Can you please explain the solution of my problem. ??

[darova]

Can you explain what those values represent?

[darova]

Here is an example using pcolor

[x,y] = meshgrid(-2:0.1:2);
stress = x.^2+y.^2;
pcolor(x,y,stress)

[Shubhang Jain]

Sigma_X, Sigma_Y , Tau_XY are the stress components values at each nodes (18) in this case . X and Y axis are the boundary of the problem .

Answer 1

[darova]

Reshape your data to create a matrix

stress = reshape(Sigma_X,6,3);

And use pcolor

[x,y] = meshgrid(1:3,1:6);
Sigma_X = 1.0e+03 * [-4299.57587958638 	-3965.44553517593   -3204.63754334495	-2292.69609509300	-1230.97384069952	-654.515387137445	7.73070496506989e-12  8.64019966684282e-12	1.11413100967184e-11	  1.02318153949454e-11	2.84217094304040e-12 	-1.93267624126747e-12	4299.57587958640	3965.44553517595	3204.63754334497	2292.69609509302	1230.97384069952	654.515387137441];
Sigma_Y = [-321.901962451787	-171.276004860075	-78.8000793342068	-101.167806049306	-108.580704946648	-151.775909194735	-43.0911561819484	-52.5066078518920	-96.8460482931672	-92.9813362559606	-169.657107326625	-285.121579205828	235.719650087890	66.2627891562907	-114.892017252128	-84.7948664626148	-230.733509706602	-418.467249216922];
Tau_XY =[ -795.818460340878	-847.839977192489	-889.271783271534	-792.615142694200	-700.220084886960	-693.891956884487	-669.856459330285	-655.783844638485	-684.491978609799	-643.356643356871	-610.298792117163	-661.157024793497	-543.894458319692	-463.727712084480	-479.712173948064	-494.098144019542	-520.377499347365	-628.422092702508] ;

stress = reshape(Sigma_X,6,3);
pcolor(x,y,stress)

[Shubhang Jain]

Thanks for the tip . It really helped .

[darova]