Page 5 - Demo

Planar TrussExamplefor ComrelAdd-onRCPConsult, 2022-2025Page 5// Plot a main title of plot_pmsgctl={8.08.0.120105};_pmsgstr=\draw;// Plot initial and deformed states of truss with loadsb=reshape(a,nnp,2);// deflection vectorDF=XY+b*2;// deformed coordinates with scale factor 2.fori (1,nfe,1);_pline=1~6~XY[IEN[i,1],1]~XY[IEN[i,1],2]~XY[IEN[i,2],1]~XY[IEN[i,2],2]~1~11~1;draw;// cyan color for initial stateendfor;fori (1,nfe,1);_pline=1~6~DF[IEN[i,1],1]~DF[IEN[i,1],2]~DF[IEN[i,2],1]~DF[IEN[i,2],2]~1~12~1;draw;// red color for initial stateendfor;fori (1,nfe,1);x=(XY[IEN[i,1],1]+XY[IEN[i,2],1])/2;y=(XY[IEN[i,1],2]+XY[IEN[i,2],2])/2;_pmsgctl=x~y~.07~0~1~15~1;_pmsgstr=sprintf(\draw;endfor;fori (1,nnp,1);_pmsgctl=XY[i,1]~XY[i,2]~.07~0~1~0~1;_pmsgstr=sprintf(\draw;endfor;_pmsgctl=DF[4,1]~DF[4,2]0.2~.07~0~1~0~1;_pmsgstr=sprintf(\draw;// value of vertical deflection in node N04fd=zeros(ndof,1);fori (1,rows(P),1);fd[14+2*i]=P[i];// vertical forces [N]endfor;fd=reshape(fd,nnp,ned)*0.000025;// forces with scale factor 0.000025n=0;fori (1,nnp,1);iffd[i,2]>0.1;n=n+1;_parrow=XY[i,1]+fd[i,1]~XY[i,2]+fd[i,2]~XY[i,1]~XY[i,2]~3~0.1~1~9~1~6~2;draw;_pmsgctl=XY[i,1]+fd[i,1]+0.1~XY[i,2]+fd[i,2]~.07~0~1~0~1;_pmsgstr=sprintf(\draw;endif;endfor;elseifnf ==2;//! Internal member forces N// Plot a main title of plot_pmsgctl={8.08.0.120105};_pmsgstr=\draw;// Plot internal member forces Nll=ones(4,1);cc={1,1,1,1};ww={1,1,1,1};tt={6,6,6,6};fori (1,nfe,1);l=leng[i]/2;n=abs(NEL[i])*0.0000005;X=-l~l~l~-l;Y=-n~-n~n~n;co=cos(theta[i]);// cosine inclination anglesi=sin(theta[i]);// sine inclination angleXrot=X*coY*si;Yrot=X*si+Y*co;x=(XY[IEN[i,1],1]+XY[IEN[i,2],1])/2;y=(XY[IEN[i,1],2]+XY[IEN[i,2],2])/2;xo=Xrot[1]|Xrot[2]|Xrot[3]|Xrot[4];xe=Xrot[2]|Xrot[3]|Xrot[4]|Xrot[1];yo=Yrot[1]|Yrot[2]|Yrot[3]|Yrot[4];ye=Yrot[2]|Yrot[3]|Yrot[4]|Yrot[1];ifNEL[i]>=0;cn=12;else;cn=9;endif;_pline=ll~tt~xo +x~yo +y~xe +x~ye +y~ll~cc*cn~ww;draw;_pmsgctl=x-0.5~y~.05~0~1~0~1;_pmsgstr=sprintf(\draw;endfor;elseifnf ==3;//! Internal member Tensile stress// Plot a main title of plot_pmsgctl={6.08.0.120105};_pmsgstr=\draw;minA=minc(abs(NEL./reshape(area,nfe,1)*1e6));maxA=maxc(abs(NEL./reshape(area,nfe,1)*1e6));dA=maxA-minA;// Plot internal member von Mises stressll=ones(4,1);cc={1,1,1,1};ww={1,1,1,1};tt={6,6,6,6};fori (1,nfe,1);Su=abs(NEL[i]/area[i]*1e6);l=leng[i]/2;n=Su/dA*0.2;X=l~l~l~l;Y=n~n~n~n;co=cos(theta[i]);// cosine inclination anglesi=sin(theta[i]);// sine inclination angleXrot=X*co-Y*si;Yrot=X*si+Y*co;x=(XY[IEN[i,1],1]+XY[IEN[i,2],1])/2;y=(XY[IEN[i,1],2]+XY[IEN[i,2],2])/2;xo=Xrot[1]|Xrot[2]|Xrot[3]|Xrot[4];xe=Xrot[2]|Xrot[3]|Xrot[4]|Xrot[1];yo=Yrot[1]|Yrot[2]|Yrot[3]|Yrot[4];ye=Yrot[2]|Yrot[3]|Yrot[4]|Yrot[1];_pline=ll~tt~xo +x~yo +y~xe +x~ye +y~ll~cc*14~ww;draw;_pmsgctl=x0.5~y~.05~0~1~0~1;_pmsgstr=sprintf(\draw;endfor;endif;// Invoked version of Gauss Engine_pmsgctl={-1.9-1.52.080103};_pmsgstr=StrurelEngine;draw;// Show/Export a plotifStrurelPlotMode==2orStrurelPlotMode==3;ifStrurelPlotMode==2;q=\q\else;q=\endif;StrurelPlotCount=StrurelPlotCount+1;graphprt(\StrurelPlotName,StrurelPlotCount,StrurelPlotType)$+q);endif;endwind;endfor;endif;endif;// begin of plot blockretp(uout);// vertical deflection at bottom center (N04)endp;// end of procedure