BEGIN{FS=",";i=0;} NR==2,NR==182{ i=i+1; x[i]=$2 y[i]=$3 dx[i]=$4 dy[i]=$5 N[i]=$6 S[i]=$7 M[i]=$8 s_in[i]=$10 s_ou[i]=$11 } END{ nd=i pi=3.141592654 ww=0.047 Eg=2.8 tauw=158 fnameW="fig_eps_disp.eps" xmax=100;xmin=0;dsx=20 printf "reset\n" printf "set terminal postscript eps enhanced font \"Helvetica\" 24\n" printf "set output \"%s\"\n",fnameW # printf "set datafile separator \",\"\n" printf "set xlabel \"Displacement (mm)\"\n" printf "set ylabel \"Location (y/R)\"\n" printf "set xtics %g\n",dsx printf "set format x \"%%.0f\"\n" printf "set ytics 0.25\n" printf "set format y \"%%.2f\"\n" printf "set xrange [%g:%g]\n",xmin,xmax printf "set yrange [-1:1]\n" printf "set grid xtics ytics mxtics mytics\n" printf "set xzeroaxis linetype 1\n" #Horizontal line printf "set label \"W=%.3fMPa R=1536mm t=22mm Eg=%.1fMPa\" at graph 0.0,graph 1.05 left\n",ww,Eg printf "plot \\\n" printf "\"-\" with lines linetype 1 linewidth 5 title \"{/Symbol D}X\", \\\n" printf "\"-\" with lines linetype 2 linewidth 5 title \"{/Symbol d}y\"\n" for(i=1;i<=nd;i++){ printf "%g %g\n",2*dx[i],y[i]/1536 } printf "e\n" for(i=1;i<=nd;i++){ printf "%g %g\n",-dy[i],y[i]/1536 } printf "e\n" fnameW="fig_eps_dist_SR.eps" xmin1=-15000;xmax1=15000;dsx1=5000 xmin2=-150;xmax2=150;dsx2=50 printf "reset\n" printf "set terminal postscript eps enhanced font \"Helvetica\" 24\n" printf "set output \"%s\"\n",fnameW # printf "set datafile separator \",\"\n" printf "set xlabel \"Moment (N-mm/mm)\"\n" printf "set x2label \"Axial force, Shear force (N/mm)\"\n" printf "set xtics %g\n",dsx1 printf "set x2tics %g\n",dsx2 printf "set format x \"%%.0f\"\n" printf "set ylabel \"Location (y/R)\"\n" printf "set ytics 0.25\n" printf "set format y \"%%.2f\"\n" printf "set xrange [%g:%g]\n",xmin1,xmax1 printf "set x2range [%g:%g]\n",xmin2,xmax2 printf "set yrange [-1:1]\n" printf "set grid xtics ytics mxtics mytics\n" printf "set xzeroaxis linetype 1\n" #Horizontal line printf "set yzeroaxis linetype 1\n" #Vertical line printf "plot \\\n" printf "\"-\" axis x2y1 with lines linetype 2 linewidth 5 title \"N\", \\\n" printf "\"-\" axis x2y1 with lines linetype 3 linewidth 5 title \"S\", \\\n" printf "\"-\" axis x1y1 with lines linetype 1 linewidth 5 title \"M\"\n" for(i=1;i<=nd;i++){printf "%g %g\n",N[i],y[i]/1536} printf "e\n" for(i=1;i<=nd;i++){printf "%g %g\n",S[i],y[i]/1536} printf "e\n" for(i=1;i<=nd;i++){printf "%g %g\n",M[i],y[i]/1536} printf "e\n" fnameW="fig_eps_stress.eps" xmax=200;xmin=-200;dsx=50 printf "reset\n" printf "set terminal postscript eps enhanced font \"Helvetica\" 24\n" printf "set output \"%s\"\n",fnameW # printf "set datafile separator \",\"\n" printf "set xlabel \"Normal stress (MPa)\"\n" printf "set ylabel \"Location (y/R)\"\n" printf "set xtics %g\n",dsx printf "set format x \"%%.0f\"\n" printf "set ytics 0.25\n" printf "set format y \"%%.2f\"\n" printf "set xrange [%g:%g]\n",xmin,xmax printf "set yrange [-1:1]\n" printf "set grid xtics ytics mxtics mytics\n" printf "set xzeroaxis linetype 1\n" #Horizontal line printf "set yzeroaxis linetype 1\n" #Vertical line printf "set label \"W=%.3fMPa R=1536mm t=22mm Eg=%.1fMPa\" at graph 0.0,graph 1.05 left\n",ww,Eg printf "plot \\\n" printf "\"-\" with lines linetype 1 linewidth 5 title \"Inside (with water)\", \\\n" printf "\"-\" with lines linetype 2 linewidth 5 title \"Outside (with water)\", \\\n" printf "\"-\" with lines linetype 3 linewidth 5 title \"Inside (no water)\", \\\n" printf "\"-\" with lines linetype 4 linewidth 5 title \"Outside (no water)\"\n" for(i=1;i<=nd;i++){ printf "%g %g\n",s_in[i]+0.9*1536/22,y[i]/1536 } printf "e\n" for(i=1;i<=nd;i++){ printf "%g %g\n",s_ou[i]+0.9*1536/22,y[i]/1536 } printf "e\n" for(i=1;i<=nd;i++){ printf "%g %g\n",s_in[i],y[i]/1536 } printf "e\n" for(i=1;i<=nd;i++){ printf "%g %g\n",s_ou[i],y[i]/1536 } printf "e\n" fnameW="fig_eps_dist_tau.eps" xmax=200;xmin=0;ds=20 printf "reset\n" printf "set terminal postscript eps enhanced font \"Helvetica\" 24\n" printf "set output \"%s\"\n",fnameW # printf "set datafile separator \",\"\n" printf "set xlabel \"Shearing stress of fillet welded joint (MPa)\"\n" printf "set ylabel \"Location (y/R)\"\n" printf "set xtics %g\n",ds printf "set format x \"%%.0f\"\n" printf "set ytics 0.25\n" printf "set format y \"%%.2f\"\n" printf "set xrange [%g:%g]\n",0,xmax printf "set yrange [-1:1]\n" printf "set grid xtics ytics mxtics mytics\n" printf "set label \"W=%.3fMPa R=1536mm t=22mm Eg=%.1fMPa\" at graph 0.0,graph 1.05 left\n",ww,Eg printf "set label \"{/Symbol t}w=%.0fMPa\" at first %g,first %g left rotate by 90\n",tauw,tauw+5,-0.7 printf "plot \\\n" printf "\"-\" with lines linetype 1 linewidth 5 title \"WL=95%% (with water)\", \\\n" printf "\"-\" with lines linetype 2 linewidth 5 title \"WL=65%% (no water)\", \\\n" printf "\"-\" with lines linetype 5 linewidth 5 title \"Shear strength of weld joint\"\n" sn=sin(81/180*pi) cs=cos(81/180*pi) s=3.2 L=0.95 for(i=1;i<=nd;i++){ NN1=(N[i]+1536*0.9)*sn*3.2/23.2 NM1=-M[i]*sn*13.32/1179.633*3.2 VV1=V[i]*sn*3.2/23.2 NN2=(N[i]+1536*0.9)*cs*3.2/23.2 NM2=-M[i]*cs*13.32/1179.633*3.2 VV2=V[i]*cs*3.2/23.2 tau1=(NN1+NM1)/(0.7*s*L*2) tau2=(NN2+NM2)/(0.7*s*L) tau3=(VV1+VV2)/(0.7*s*L) tau=sqrt(tau1^2+tau2^2+tau3^2) printf "%g %g\n",tau,y[i]/1536 } printf "e\n" L=0.65 for(i=1;i<=nd;i++){ NN1=(N[i])*sn*3.2/23.2 NM1=-M[i]*sn*13.32/1179.633*3.2 VV1=V[i]*sn*3.2/23.2 NN2=(N[i])*cs*3.2/23.2 NM2=-M[i]*cs*13.32/1179.633*3.2 VV2=V[i]*cs*3.2/23.2 tau1=(NN1+NM1)/(0.7*s*L*2) tau2=(NN2+NM2)/(0.7*s*L) tau3=(VV1+VV2)/(0.7*s*L) tau=sqrt(tau1^2+tau2^2+tau3^2) printf "%g %g\n",tau,y[i]/1536 } printf "e\n" printf "%g %g\n",tauw,y[1]/1536 printf "%g %g\n",tauw,y[nd]/1536 printf "e\n" printf "q\n" }