(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 4.2' Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. *******************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 179228, 7962]*) (*NotebookOutlinePosition[ 180326, 7997]*) (* CellTagsIndexPosition[ 180282, 7993]*) (*WindowFrame->Normal*) Notebook[{ Cell[TextData[{ StyleBox["Here are the instructions:\n\n1. Complete the missing code in \ Cells 1 and 3 where you see blue text enclosed \n in (* ... *). Red \ statements are for testing. Once these cells work go to step 2.\n2. Execute \ Cells 1 through 5 in any order (twice each if you want to get rid of \n \ annoying warning messages) to initialize all necessary modules. \n \ Shortcut for ", FontFamily->"Palatino"], StyleBox["Mathematica 4.0 ", FontFamily->"Palatino", FontSlant->"Italic"], StyleBox[" users: select Kernel ->Evaluation->Evaluate\n \ Initialization (In ", FontFamily->"Palatino"], StyleBox["Mathematica", FontFamily->"Palatino", FontSlant->"Italic"], StyleBox[" 2.2, Action->Evaluate Initialization).\n3. Cell 6 contains the \ driver program to run the plane frame of Ex 22.3. Run by\n selecting the \ cell and executing.\n4. Execute Cell 6A to generate deformed shape plot \ frames. Animate by double\n clicking one of the frames. Animation may \ be slowed down by the playback\n button controls on the left side of the \ bottom window bar.\n5. If the plots produced by running a main program \ appear too small, click \n on the top one (only) with the mouse, \ grap a corner \"handle\" and enlarge it. \n Then rerun the program.", FontFamily->"Palatino"] }], "Text", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0.463401, 1, 0.300023]], Cell[TextData[{ "Cell 1. Assembly module for plane frame structure. The element stiffness \ module \n", StyleBox["PlaneBeamColumn2Stiffness ", FontFamily->"Courier", FontWeight->"Bold"], "that supports the assembler has been \ndescribed in Chapter 21." }], "Text", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1]], Cell[CellGroupData[{ Cell[TextData[{ "PlaneBeamColumn2Stiffness[ncoor_,mprop_,fprop_,opt_]:= Module[\n \ {x1,x2,y1,y2,x21,y21,Em,Gm,rho,alpha,A,Izz,num,c,s,L,LL,\n \ LLL,ra,rb,T,Kebar,Ke}, \n {{x1,y1},{x2,y2}}=ncoor; {x21,y21}={x2-x1,y2-y1}; \ \n {Em,Gm,rho,alpha}=mprop; {A,Izz}=fprop; {num}=opt;\n \ LL=Simplify[x21^2+y21^2]; L=PowerExpand[Sqrt[LL]]; LLL=L*LL; \n c=x21/L; \ s=y21/L; ra=Em*A/L; rb= 2*Em*Izz/LLL; \n Kebar= ra*{\n { \ 1,0,0,-1,0,0},{0,0,0,0,0,0},{0,0,0,0,0,0},\n {-1,0,0, \ 1,0,0},{0,0,0,0,0,0},{0,0,0,0,0,0}} +\n rb*{\n { 0,0,0,0,0,0},{0, \ 6, 3*L,0,-6, 3*L},{0,3*L,2*LL,0,-3*L, LL},\n { 0,0,0,0,0,0},{0,-6,-3*L,0, \ 6,-3*L},{0,3*L, LL,0,-3*L,2*LL}};\n \ T={{c,s,0,0,0,0},{-s,c,0,0,0,0},{0,0,1,0,0,0},\n \ {0,0,0,c,s,0},{0,0,0,-s,c,0},{0,0,0,0,0,1}};\n Ke=Transpose[T].Kebar.T; If \ [num,Ke=N[Ke]];\n Return[Ke]\n ];\n\ PlaneFrameMasterStiffness[nodcoor_,elenod_,\n \ elemat_,elefab_,eleopt_]:=Module[\n \ {numele=Length[elenod],numnod=Length[nodcoor],\n", " e,eNL,eftab,ni,nj,i,j,ncoor,mprop,fprop,opt,Ke,K},\n \ K=Table[0,{3*numnod},{3*numnod}];\n For [e=1, e<=numele, e++,\n", StyleBox[" eNL=elenod[[e]]; {ni,nj}=eNL; \n \ eftab={3*ni-2,3*ni-1,3*ni,3*nj-2,3*nj-1,3*nj}; \n \ ncoor={nodcoor[[ni]],nodcoor[[nj]]}; \n mprop=elemat[[e]]; \ fprop=elefab[[e]]; opt=eleopt; \n \ Ke=PlaneBeamColumn2Stiffness[ncoor,mprop,fprop,opt];\n Print[\"Ke for \ element \",e,\" is \",Ke//MatrixForm];\n neldof=Length[Ke];\n For \ [i=1, i<=neldof, i++, ii=eftab[[i]];\n For [j=i, j<=neldof, j++, \ jj=eftab[[j]];\n K[[jj,ii]]=K[[ii,jj]]+=Ke[[i,j]] \n ];\ \n ];", FontColor->RGBColor[0, 0, 1]], "\n ];\n Return[K]", "\n ];\n", StyleBox["nodcoor={{0,0},{10,0},{10,10}};\nelenod= {{1,2},{2,3},{1,3}}; \ elemat= Table[{100,0,0,0},{3}];\nelefab= {{1,10},{1/2,10},{2*Sqrt[2],10}}; \ eleopt= {True};\n\ K=PlaneFrameMasterStiffness[nodcoor,elenod,elemat,elefab,eleopt];\n\ Print[\"Master Stiffness of Example Frame:\"];\nPrint[K//MatrixForm]; \ Print[Chop[Eigenvalues[K]]]; ", FontColor->RGBColor[1, 0, 0]] }], "Input", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0.699214, 1, 0.0500191]], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Ke for element \"\>", "\[InvisibleSpace]", "1", "\[InvisibleSpace]", "\<\" is \"\>", "\[InvisibleSpace]", TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"10.`", "0.`", "0.`", \(-10.`\), "0.`", "0.`"}, {"0.`", "12.`", "60.`", "0.`", \(-12.`\), "60.`"}, {"0.`", "60.`", "400.`", "0.`", \(-60.`\), "200.`"}, {\(-10.`\), "0.`", "0.`", "10.`", "0.`", "0.`"}, {"0.`", \(-12.`\), \(-60.`\), "0.`", "12.`", \(-60.`\)}, {"0.`", "60.`", "200.`", "0.`", \(-60.`\), "400.`"} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]}], SequenceForm[ "Ke for element ", 1, " is ", MatrixForm[ {{10.0, 0.0, 0.0, -10.0, 0.0, 0.0}, {0.0, 12.0, 60.0, 0.0, -12.0, 60.0}, {0.0, 60.0, 400.0, 0.0, -60.0, 200.0}, {-10.0, 0.0, 0.0, 10.0, 0.0, 0.0}, {0.0, -12.0, -60.0, 0.0, 12.0, -60.0}, { 0.0, 60.0, 200.0, 0.0, -60.0, 400.0}}]], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Ke for element \"\>", "\[InvisibleSpace]", "2", "\[InvisibleSpace]", "\<\" is \"\>", "\[InvisibleSpace]", TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"12.`", "0.`", \(-60.`\), \(-12.`\), "0.`", \(-60.`\)}, {"0.`", "5.`", "0.`", "0.`", \(-5.`\), "0.`"}, {\(-60.`\), "0.`", "400.`", "60.`", "0.`", "200.`"}, {\(-12.`\), "0.`", "60.`", "12.`", "0.`", "60.`"}, {"0.`", \(-5.`\), "0.`", "0.`", "5.`", "0.`"}, {\(-60.`\), "0.`", "200.`", "60.`", "0.`", "400.`"} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]}], SequenceForm[ "Ke for element ", 2, " is ", MatrixForm[ {{12.0, 0.0, -60.0, -12.0, 0.0, -60.0}, {0.0, 5.0, 0.0, 0.0, -5.0, 0.0}, {-60.0, 0.0, 400.0, 60.0, 0.0, 200.0}, {-12.0, 0.0, 60.0, 12.0, 0.0, 60.0}, {0.0, -5.0, 0.0, 0.0, 5.0, 0.0}, {-60.0, 0.0, 200.0, 60.0, 0.0, 400.0}}]], Editable->False]], "Print"], Cell[BoxData[ InterpretationBox[ RowBox[{"\<\"Ke for element \"\>", "\[InvisibleSpace]", "3", "\[InvisibleSpace]", "\<\" is \"\>", "\[InvisibleSpace]", TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"12.121320343559642`", "7.878679656440358`", \(-21.213203435596427`\), \ \(-12.121320343559642`\), \(-7.878679656440358`\), \(-21.213203435596427`\)}, {"7.878679656440358`", "12.121320343559642`", "21.213203435596427`", \(-7.878679656440358`\), \ \(-12.121320343559642`\), "21.213203435596427`"}, {\(-21.213203435596427`\), "21.213203435596427`", "282.842712474619`", "21.213203435596427`", \(-21.213203435596427`\), "141.4213562373095`"}, {\(-12.121320343559642`\), \(-7.878679656440358`\), "21.213203435596427`", "12.121320343559642`", "7.878679656440358`", "21.213203435596427`"}, {\(-7.878679656440358`\), \(-12.121320343559642`\), \ \(-21.213203435596427`\), "7.878679656440358`", "12.121320343559642`", \(-21.213203435596427`\)}, {\(-21.213203435596427`\), "21.213203435596427`", "141.4213562373095`", "21.213203435596427`", \(-21.213203435596427`\), "282.842712474619`"} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]}], SequenceForm[ "Ke for element ", 3, " is ", MatrixForm[ {{12.121320343559642, 7.8786796564403581, -21.213203435596427, -12.121320343559642, \ -7.8786796564403581, -21.213203435596427}, {7.8786796564403581, 12.121320343559642, 21.213203435596427, -7.8786796564403581, -12.121320343559642, 21.213203435596427}, {-21.213203435596427, 21.213203435596427, 282.84271247461902, 21.213203435596427, -21.213203435596427, 141.42135623730951}, {-12.121320343559642, -7.8786796564403581, 21.213203435596427, 12.121320343559642, 7.8786796564403581, 21.213203435596427}, {-7.8786796564403581, -12.121320343559642, \ -21.213203435596427, 7.8786796564403581, 12.121320343559642, -21.213203435596427}, {-21.213203435596427, 21.213203435596427, 141.42135623730951, 21.213203435596427, -21.213203435596427, 282.84271247461902}}]], Editable->False]], "Print"], Cell[BoxData[ \("Master Stiffness of Example Frame:"\)], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"22.121320343559642`", "7.878679656440358`", \(-21.213203435596427`\), \(-10.`\), "0.`", "0.`", \(-12.121320343559642`\), \ \(-7.878679656440358`\), \(-21.213203435596427`\)}, {"7.878679656440358`", "24.121320343559642`", "81.21320343559643`", "0.`", \(-12.`\), "60.`", \(-7.878679656440358`\), \(-12.121320343559642`\), "21.213203435596427`"}, {\(-21.213203435596427`\), "81.21320343559643`", "682.842712474619`", "0.`", \(-60.`\), "200.`", "21.213203435596427`", \(-21.213203435596427`\), "141.4213562373095`"}, {\(-10.`\), "0.`", "0.`", "22.`", "0.`", \(-60.`\), \(-12.`\), "0.`", \(-60.`\)}, {"0.`", \(-12.`\), \(-60.`\), "0.`", "17.`", \(-60.`\), "0.`", \(-5.`\), "0.`"}, {"0.`", "60.`", "200.`", \(-60.`\), \(-60.`\), "800.`", "60.`", "0.`", "200.`"}, {\(-12.121320343559642`\), \(-7.878679656440358`\), "21.213203435596427`", \(-12.`\), "0.`", "60.`", "24.121320343559642`", "7.878679656440358`", "81.21320343559643`"}, {\(-7.878679656440358`\), \(-12.121320343559642`\), \ \(-21.213203435596427`\), "0.`", \(-5.`\), "0.`", "7.878679656440358`", "17.121320343559642`", \(-21.213203435596427`\)}, {\(-21.213203435596427`\), "21.213203435596427`", "141.4213562373095`", \(-60.`\), "0.`", "200.`", "81.21320343559643`", \(-21.213203435596427`\), "682.842712474619`"} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Print"], Cell[BoxData[ \({1121.7041041422942`, 555.3384692425468`, 531.6521065780339`, 46.61291206166076`, 22.497103073280275`, 14.366011225661897`, 0, 0, 0}\)], "Print"] }, Open ]], Cell["\<\ Cell 2. These modules apply the displacement boundary conditions \ by modifying the force vector and the master stiffness matrix. This implementation of \ ModifyNodeForcesFor DBC cannot handle nonzero prescribed displacements.\ \>", "Text", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1]], Cell[CellGroupData[{ Cell[TextData[{ "ModifiedMasterStiffness[pdof_,K_] := Module[\n \ {i,j,k,n=Length[K],np=Length[pdof],Kii,Kmod}, Kmod=K; \n For [k=1,k<=np,k++, \ i=pdof[[k]]; Kii=Kmod[[i,i]]; \n For [j=1,j<=n,j++, \ Kmod[[i,j]]=Kmod[[j,i]]=0];\n If [Kii==0, Kii=1]; Kmod[[i,i]]=Kii\n \ ]; \n Return[Kmod]\n];\n\nModifyNodeForces[pdof_,f_] := Module[\n \ {i,k,np=Length[pdof],fmod}, fmod=f; \n For [k=1,k<=np,k++, i=pdof[[k]]; \ fmod[[i]]=0];\n Return[fmod]\n];\n\n", StyleBox["K=Array[Kij,{6,6}];\nPrint[\"Assembled Master \ Stiffness:\"];Print[K//MatrixForm];\nK=ModifiedMasterStiffness[{1,2,4},K];\n\ Print[\"Master Stiffness Modified For Displacement B.C.:\"];\n\ Print[K//MatrixForm];\nf=Array[fi,{6}];\nPrint[\"Node Force Vector:\"]; \ Print[f];\nf=ModifyNodeForces[{1,2,4},f];\nPrint[\"Node Force Vector Modified \ For Displacement B.C.:\"];\nPrint[f];", FontColor->RGBColor[1, 0, 0]] }], "Input", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0.699214, 1, 0.0500191]], Cell[BoxData[ \("Assembled Master Stiffness:"\)], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {\(Kij[1, 1]\), \(Kij[1, 2]\), \(Kij[1, 3]\), \(Kij[1, 4]\), \(Kij[1, 5]\), \(Kij[1, 6]\)}, {\(Kij[2, 1]\), \(Kij[2, 2]\), \(Kij[2, 3]\), \(Kij[2, 4]\), \(Kij[2, 5]\), \(Kij[2, 6]\)}, {\(Kij[3, 1]\), \(Kij[3, 2]\), \(Kij[3, 3]\), \(Kij[3, 4]\), \(Kij[3, 5]\), \(Kij[3, 6]\)}, {\(Kij[4, 1]\), \(Kij[4, 2]\), \(Kij[4, 3]\), \(Kij[4, 4]\), \(Kij[4, 5]\), \(Kij[4, 6]\)}, {\(Kij[5, 1]\), \(Kij[5, 2]\), \(Kij[5, 3]\), \(Kij[5, 4]\), \(Kij[5, 5]\), \(Kij[5, 6]\)}, {\(Kij[6, 1]\), \(Kij[6, 2]\), \(Kij[6, 3]\), \(Kij[6, 4]\), \(Kij[6, 5]\), \(Kij[6, 6]\)} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Print"], Cell[BoxData[ \("Master Stiffness Modified For Displacement B.C.:"\)], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {\(Kij[1, 1]\), "0", "0", "0", "0", "0"}, {"0", \(Kij[2, 2]\), "0", "0", "0", "0"}, {"0", "0", \(Kij[3, 3]\), "0", \(Kij[3, 5]\), \(Kij[3, 6]\)}, {"0", "0", "0", \(Kij[4, 4]\), "0", "0"}, {"0", "0", \(Kij[5, 3]\), "0", \(Kij[5, 5]\), \(Kij[5, 6]\)}, {"0", "0", \(Kij[6, 3]\), "0", \(Kij[6, 5]\), \(Kij[6, 6]\)} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Print"], Cell[BoxData[ \("Node Force Vector:"\)], "Print"], Cell[BoxData[ \({fi[1], fi[2], fi[3], fi[4], fi[5], fi[6]}\)], "Print"], Cell[BoxData[ \("Node Force Vector Modified For Displacement B.C.:"\)], "Print"], Cell[BoxData[ \({0, 0, fi[3], 0, fi[5], fi[6]}\)], "Print"] }, Open ]], Cell["\<\ Cell 3. Computation of internal forces (axial forces) for all \ elements of a plane truss, from the computed nodal displacements.\ \>", "Text", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1]], Cell[CellGroupData[{ Cell[TextData[{ "PlaneFrameIntForces[nodcoor_,elenod_,elemat_,elefab_,\n eleopt_,u_]:= \ Module[{numele=Length[elenod],\n numnod=Length[nodcoor],e,eNL,eftab,ni,nj,i,\ \n ncoor,mprop,fprop,opt,ue,p},\n p=Table[0,{numele}]; ue=Table[0,{6}];\n \ For [e=1, e<=numele, e++, \n eNL=elenod[[e]]; {ni,nj}=eNL;\n \ ncoor={nodcoor[[ni]],nodcoor[[nj]]}; ", StyleBox[" ", FontColor->RGBColor[1, 0, 0]], " \n mprop=elemat[[e]]; fprop=elefab[[e]]; opt=eleopt;\n \ eftab={3*ni-2,3*ni-1,3*ni,3*nj-2,3*nj-1,3*nj}; \n For [i=1,i<=6,i++, \ ii=eftab[[i]]; ue[[i]]=u[[ii]]];\n \ p[[e]]=PlaneBeamColumn2IntForces[ncoor,mprop,fprop,opt,ue]\n ]; \n \ Return[p]\n];\nPlaneBeamColumn2IntForces[ncoor_,mprop_,fprop_,opt_,ue_]:=\n \ Module[{x1,x2,y1,y2,x21,y21,Em,Gm,rho,alpha,A,Izz,num,LL,L,\n \ dvy,cv1,cv2,pe=Table[0,{3}]},\n", StyleBox[" {{x1,y1},{x2,y2}}=ncoor; {x21,y21}={x2-x1,y2-y1};\n \ {Em,Gm,rho,alpha}=mprop; {A,Izz}=fprop; {num}=opt;\n If \ [num,{x21,y21,Em,A,Izz}=N[{x21,y21,Em,A,Izz}]]; \n LL=x21^2+y21^2; \ L=PowerExpand[Sqrt[LL]];\n \ pe[[1]]=Em*A*(x21*(ue[[4]]-ue[[1]])+y21*(ue[[5]]-ue[[2]]))/LL;\n \ duy=(x21*(ue[[5]]-ue[[2]])-y21*(ue[[4]]-ue[[1]]))/L;\n cv1= \ 6*duy/LL+2*(-2*ue[[3]]-ue[[6]] )/L;\n cv2=-6*duy/LL+2*( \ ue[[3]]+2*ue[[6]])/L;\n pe[[2]]=Em*Izz*cv1; pe[[3]]=Em*Izz*cv2;", FontColor->RGBColor[0, 0, 1]], "\n Return[pe]\n]; \n", StyleBox["ClearAll[L,Em,A1,A2,Izz1,Izz2,Izz3,uz1,uz2,uz2,uz3]; \n\ nodcoor={{0,0},{L,0},{L,L}}; elenod= {{1,2},{2,3},{1,3}};\nelemat= \ Table[{Em,0,0,0},{3}];\nelefab= {{A1,Izz1},{A2,Izz2},{A3,Izz3}}; eleopt= \ {False};\nu={0,uy1,0, 0,uy2,0, 0,uy3,0};\n\ p=PlaneFrameIntForces[nodcoor,elenod,elemat,elefab,eleopt,u];\nPrint[\"Int \ Forces of Example Frame:\"]; Print[p//MatrixForm]; ", FontColor->RGBColor[1, 0, 0]] }], "Input", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0.699214, 1, 0.0500191]], Cell[BoxData[ \("Int Forces of Example Frame:"\)], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ { "0", \(\(6\ Em\ Izz1\ \((\(-uy1\) + uy2)\)\)\/L\^2\), \(-\(\(6\ Em\ Izz1\ \((\(-uy1\) + uy2)\)\)\/L\^2\)\)}, {\(\(A2\ Em\ \((\(-uy2\) + uy3)\)\)\/L\), "0", "0"}, {\(\(A3\ Em\ \((\(-uy1\) + uy3)\)\)\/\(2\ L\)\), \(\(3\ Em\ Izz3\ \((\(-uy1\) + uy3)\)\)\/\(\@2\ L\^2\)\), \(-\(\(3\ Em\ Izz3\ \ \((\(-uy1\) + uy3)\)\)\/\(\@2\ L\^2\)\)\)} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Print"], Cell["\<\ Cell 4 collectively groups plot modules that support preprocessing, \ that is, draw pictures of the FEM model and BCs. It is broken down inro three cells: 4A, 4B, 4C because of the length of \ these modules.\ \>", "Text", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1]], Cell["\<\ Cell 4A. PlotLineElements draws only the element (bars) without \ any labels.\ \>", "Text", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1]], Cell[CellGroupData[{ Cell[TextData[{ "PlotLineElements[nodcoor_,elenod_,aspect_,title_]:= Module[\n \ {eNL,e,n,ni,nj,nc,xyc,sides={},numele=Length[elenod],\n \ numnod=Length[nodcoor],x,y,xmin,xmax,ymin,ymax,aratio},\n \ x=Table[nodcoor[[i,1]],{i,1,numnod}]; \n \ y=Table[nodcoor[[i,2]],{i,1,numnod}];\n \ {xmin,xmax,ymin,ymax}=N[{Min[x],Max[x],Min[y],Max[y]}];\n For \ [e=1,e<=numele,e++, eNL=elenod[[e]]; nc=Length[eNL]; \n If [nc!=2, \ Continue[]]; {ni,nj}=eNL; \n xyc={{x[[ni]],y[[ni]]},{x[[nj]],y[[nj]]}}; \ \n sides=AppendTo[sides,Graphics[Line[xyc]]]\n ];\n If [aspect>0, \ aratio=aspect, aratio=(ymax-ymin)/(xmax-xmin)];\n If [aspect<0, \ aratio=Automatic];\n Show[Graphics[AbsoluteThickness[2]],\n \ Graphics[RGBColor[0,0,0]],\n sides, Background->RGBColor[1,1,0],\n \ AspectRatio->aratio,PlotLabel->title ];\n ];\n", StyleBox[" \nClearAll[L,H];\nL=10000; H=6000; \n\ NodeCoordinates={{-L,0},{-L,H/2},{-L,H},{-L/2,H},{0,H},{L/2,H},\n \ {L,H},{L,H/2},{L,0},{0,H/2},{0,0}};\nElemNodeLists= \ {{1,2},{2,3},{3,4},{4,5},{5,6},{6,7},{7,8},\n \ {8,9},{5,10},{10,11}};\nnumnod=Length[nodcoor]; numele=Length[elenod];\n\ aspect=0; elabinfo={True,0.12}; nlabinfo={True,0.06};\n\n\ PlotLineElements[NodeCoordinates,ElemNodeLists,\n aspect,\"test mesh\"];", FontColor->RGBColor[1, 0, 0]] }], "Input", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0.699214, 1, 0.0500191]], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .3 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Background color 1 1 0 r MFill % Scaling calculations 0.5 4.7619e-05 0.00714286 4.7619e-05 [ [.5 .3125 -29 0 ] [.5 .3125 29 12 ] [ 0 0 0 0 ] [ 1 .3 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 0 1 r gsave .5 .3125 -90 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 0.000 0.000 rmoveto 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor (test) show 93.000 13.000 moveto (mesh) show 117.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore 0 0 m 1 0 L 1 .3 L 0 .3 L closepath clip newpath 0 0 0 r 2 Mabswid [ ] 0 setdash .02381 .00714 m .02381 .15 L s .02381 .15 m .02381 .29286 L s .02381 .29286 m .2619 .29286 L s .2619 .29286 m .5 .29286 L s .5 .29286 m .7381 .29286 L s .7381 .29286 m .97619 .29286 L s .97619 .29286 m .97619 .15 L s .97619 .15 m .97619 .00714 L s .5 .29286 m .5 .15 L s .5 .15 m .5 .00714 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageSize->{527.688, 158.25}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{-0, 1}, {0, 1}}, ImageCacheValid->False] }, Open ]], Cell["\<\ Cell 4B. PlotLineElementsAndNodes drwas element and nodes, and \ optionally labels them.\ \>", "Text", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1]], Cell[CellGroupData[{ Cell[TextData[{ "PlotLineElementsAndNodes[nodcoor_,elenod_,aspect_,title_,\n \ elabinfo_,nlabinfo_]:= Module[{i,j,k,e,eNL,n,ni,nj,x,y,\n \ delx,dely,xmin,xmax,ymin,ymax,xy0,epsx,epsy,ex,ey,\n \ numele=Length[elenod],elabels,nlabels,frade,fradn,\n \ numnod=Length[nodcoor],rade,radn,elab,nlab,pnod=pnlab={},\n \ pesid=pelab=pecirc=pedisk={},aratio},\n \ x=Table[nodcoor[[i,1]],{i,1,numnod}]; \n \ y=Table[nodcoor[[i,2]],{i,1,numnod}];\n \ {xmin,xmax,ymin,ymax}=N[{Min[x],Max[x],Min[y],Max[y]}];\n \ {delx,dely}={xmax-xmin,ymax-ymin};\n {elabels,frade}=elabinfo; \ {nlabels,fradn}=nlabinfo;\n radn=fradn*Min[delx,dely]; \ {ex,ey}={2*radn,1.6*radn}; \n rade=frade*Min[delx,dely]; \n For \ [e=1,e<=numele,e++, eNL=elenod[[e]]; n=Length[eNL];\n If [n!=2, \ Continue[]]; {ni,nj}=eNL;\n gside={nodcoor[[ni]],nodcoor[[nj]]}; \n \ xy0=(gside[[1]]+gside[[2]])/2; \n pesid= AppendTo[pesid, \ Graphics[Line[gside]]];\n If [elabels, \n \ pedisk=AppendTo[pedisk,Graphics[Disk[xy0,rade]]]; \n \ pecirc=AppendTo[pecirc, Graphics[Circle[xy0,rade]]];\n \ elab=FontForm[e,{\"Times\",11}];\n pelab= AppendTo[pelab, \ Graphics[Text[elab,xy0]]]]\n ];\n For [n=1,n<=numnod,n++, \ nlab=FontForm[n,{\"Times\",12}];\n If [nlabels, \n \ pnlab=AppendTo[pnlab, \n \ Graphics[Text[nlab,nodcoor[[n]]+{ex,ey}]]]];\n pnod= AppendTo[pnod, \n \ Graphics[Disk[nodcoor[[n]],radn]]]\n ]; \ epsx=0.04*(xmax-xmin); epsy=0.14*(ymax-ymin);\n If [aspect>0, \ aratio=aspect, aratio=(ymax-ymin)/(xmax-xmin)];\n If [aspect<0, \ aratio=Automatic];\n Show[Graphics[GrayLevel[0]],\n \ Graphics[AbsoluteThickness[3]],pesid,\n Graphics[GrayLevel[0]],pnod,\n \ Graphics[AbsoluteThickness[1]],\n \ Graphics[RGBColor[.99,.99,.99]],pedisk,\n \ Graphics[GrayLevel[0]],pecirc,pelab,\n Graphics[GrayLevel[0]],pnlab,\n \ Background->RGBColor[1,1,0],\n PlotLabel->title,\n \ PlotRange->{{xmin-epsx,xmax+epsx},{ymin-epsy,ymax+epsy}},\n \ Axes->False,AspectRatio->aratio];\n];\n\n", StyleBox["ClearAll[L,H];\nL=10000; H=6000; \n\ NodeCoordinates={{-L,0},{-L,H/2},{-L,H},{-L/2,H},{0,H},{L/2,H},\n \ {L,H},{L,H/2},{L,0},{0,H/2},{0,0}};\nElemNodeLists= \ {{1,2},{2,3},{3,4},{4,5},{5,6},{6,7},{7,8},\n \ {8,9},{5,10},{10,11}};\nnumnod=Length[nodcoor]; numele=Length[elenod];\n\ aspect=-1; elabinfo={True,0.10}; nlabinfo={True,0.04};\n\n\ PlotLineElementsAndNodes[NodeCoordinates,ElemNodeLists,aspect,\n \"test \ mesh with elem & node labels\",elabinfo,nlabinfo];", FontColor->RGBColor[1, 0, 0]] }], "Input", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0.699214, 1, 0.0500191]], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .35556 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Background color 1 1 0 r MFill % Scaling calculations 0.5 4.62963e-05 0.0388889 4.62963e-05 [ [.5 .36806 -101.125 0 ] [.5 .36806 101.125 12 ] [ 0 0 0 0 ] [ 1 .35556 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 0 1 r gsave .5 .36806 -162.125 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 0.000 0.000 rmoveto 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor (test) show 93.000 13.000 moveto (mesh) show 123.125 13.000 moveto (with) show 153.125 13.000 moveto (elem) show 183.250 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor (&) show 195.250 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor (node) show 225.250 13.000 moveto (labels) show 261.250 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore 0 0 m 1 0 L 1 .35556 L 0 .35556 L closepath clip newpath 0 g 3 Mabswid [ ] 0 setdash .03704 .03889 m .03704 .17778 L s .03704 .17778 m .03704 .31667 L s .03704 .31667 m .26852 .31667 L s .26852 .31667 m .5 .31667 L s .5 .31667 m .73148 .31667 L s .73148 .31667 m .96296 .31667 L s .96296 .31667 m .96296 .17778 L s .96296 .17778 m .96296 .03889 L s .5 .31667 m .5 .17778 L s .5 .17778 m .5 .03889 L s .03704 .03889 m .03704 .03889 .01111 0 365.73 arc F .03704 .17778 m .03704 .17778 .01111 0 365.73 arc F .03704 .31667 m .03704 .31667 .01111 0 365.73 arc F .26852 .31667 m .26852 .31667 .01111 0 365.73 arc F .5 .31667 m .5 .31667 .01111 0 365.73 arc F .73148 .31667 m .73148 .31667 .01111 0 365.73 arc F .96296 .31667 m .96296 .31667 .01111 0 365.73 arc F .96296 .17778 m .96296 .17778 .01111 0 365.73 arc F .96296 .03889 m .96296 .03889 .01111 0 365.73 arc F .5 .17778 m .5 .17778 .01111 0 365.73 arc F .5 .03889 m .5 .03889 .01111 0 365.73 arc F .99 .99 .99 r .03704 .10833 m .03704 .10833 .02778 0 365.73 arc F .03704 .24722 m .03704 .24722 .02778 0 365.73 arc F .15278 .31667 m .15278 .31667 .02778 0 365.73 arc F .38426 .31667 m .38426 .31667 .02778 0 365.73 arc F .61574 .31667 m .61574 .31667 .02778 0 365.73 arc F .84722 .31667 m .84722 .31667 .02778 0 365.73 arc F .96296 .24722 m .96296 .24722 .02778 0 365.73 arc F .96296 .10833 m .96296 .10833 .02778 0 365.73 arc F .5 .24722 m .5 .24722 .02778 0 365.73 arc F .5 .10833 m .5 .10833 .02778 0 365.73 arc F 0 g 1 Mabswid newpath .03704 .10833 .02778 0 365.73 arc s newpath .03704 .24722 .02778 0 365.73 arc s newpath .15278 .31667 .02778 0 365.73 arc s newpath .38426 .31667 .02778 0 365.73 arc s newpath .61574 .31667 .02778 0 365.73 arc s newpath .84722 .31667 .02778 0 365.73 arc s newpath .96296 .24722 .02778 0 365.73 arc s newpath .96296 .10833 .02778 0 365.73 arc s newpath .5 .24722 .02778 0 365.73 arc s newpath .5 .10833 .02778 0 365.73 arc s gsave .03704 .10833 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (1) show 1.000 setlinewidth grestore gsave .03704 .24722 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (2) show 1.000 setlinewidth grestore gsave .15278 .31667 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (3) show 1.000 setlinewidth grestore gsave .38426 .31667 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (4) show 1.000 setlinewidth grestore gsave .61574 .31667 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (5) show 1.000 setlinewidth grestore gsave .84722 .31667 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (6) show 1.000 setlinewidth grestore gsave .96296 .24722 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (7) show 1.000 setlinewidth grestore gsave .96296 .10833 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (8) show 1.000 setlinewidth grestore gsave .5 .24722 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (9) show 1.000 setlinewidth grestore gsave .5 .10833 -68.5 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (10) show 1.000 setlinewidth grestore gsave .05926 .05667 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (1) show 1.000 setlinewidth grestore gsave .05926 .19556 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (2) show 1.000 setlinewidth grestore gsave .05926 .33444 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (3) show 1.000 setlinewidth grestore gsave .29074 .33444 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (4) show 1.000 setlinewidth grestore gsave .52222 .33444 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (5) show 1.000 setlinewidth grestore gsave .7537 .33444 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (6) show 1.000 setlinewidth grestore gsave .98519 .33444 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (7) show 1.000 setlinewidth grestore gsave .98519 .19556 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (8) show 1.000 setlinewidth grestore gsave .98519 .05667 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (9) show 1.000 setlinewidth grestore gsave .52222 .19556 -69 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (10) show 1.000 setlinewidth grestore gsave .52222 .05667 -69 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (11) show 1.000 setlinewidth grestore % End of Graphics MathPictureEnd \ \>"], "Graphics", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageSize->{529.438, 188.188}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{-0, 1}, {0, 1}}, ImageCacheValid->False] }, Open ]], Cell[TextData[{ "Cell 4C. PlotBoundaryConditions tries to do like the title says: draw \ applies forces and\nsupport conditions. Not very successful yet, because ", StyleBox["Mathematica", FontSlant->"Italic"], " plotting capabilities are not\nexactly state of the art." }], "Text", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1]], Cell[TextData[StyleBox["(* Removed, it does not work properly and there \n \ is no time to fix*)", FontColor->RGBColor[1, 0, 0]]], "Input", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0.241825, 0.918593, 0.929992]], Cell["\<\ Cell 5. Plot modules that support postprocessing. As in the case \ of Calle 4, these are subdivided into Cells 5A, 5B and 5C to keep them short.\ \>", "Text", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1]], Cell["Cell 5A. PlotDeforcedShape does exactly what it says.", "Text", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1]], Cell[CellGroupData[{ Cell[TextData[{ "PlotPlaneFrameDeformedShape[nodcoor_,elenod_,u_,pcont_,\n \ aspect_,title_]:= Module[\n \ {eNL,e,i,ii,n,ni,nj,nc,xyc,pmem={},numele=Length[elenod],\n \ numnod=Length[nodcoor],x,y,eftab,ue=Table[0,{6}],\n \ xmin,xmax,ymin,ymax,xam,xax,yam,yax,xbm,xbx,ybm,ybx,\n \ amp,ampmax,ampmin,ns,aratio},\n {amp,ampmax,ampmin,ns}=pcont;\n \ x=Table[nodcoor[[i,1]],{i,1,numnod}]; \n \ y=Table[nodcoor[[i,2]],{i,1,numnod}];\n \ {xmin,xmax,ymin,ymax}=N[{Min[x],Max[x],Min[y],Max[y]}];\n For \ [e=1,e<=numele,e++, eNL=elenod[[e]]; nc=Length[eNL]; \n If [nc!=2, \ Continue[]]; {ni,nj}=eNL; \n xyc={{x[[ni]],y[[ni]]},{x[[nj]],y[[nj]]}}; \n\ eftab={3*ni-2,3*ni-1,3*ni,3*nj-2,3*nj-1,3*nj};\n For [i=1,i<=6,i++, \ ii=eftab[[i]]; ue[[i]]=u[[ii]]];\n \ pmem=Append[pmem,PlotBeamShape[xyc,ue,amp,ns]];\n \ {xam,xax,yam,yax}=BeamShapePlotBounds[xyc,ue,ampmax,ns];\n \ {xbm,xbx,ybm,ybx}=BeamShapePlotBounds[xyc,ue,ampmin,ns]; \n \ xmin=Min[xmin,xam,xbm]; ymin=Min[ymin,yam,ybm];\n xmax=Max[xmax,xax,xbx]; \ ymax=Max[ymax,yax,ybx];\n ];\n If [aspect>0, aratio=aspect, \ aratio=(ymax-ymin)/(xmax-xmin)];\n If [aspect<0, aratio=Automatic];\n \ pminmax={Graphics[RGBColor[1,1,1]],\n \ Graphics[AbsolutePointSize[1]],\n Graphics[Point[{xmin,ymin}]],\n \ Graphics[Point[{xmax,ymax}]]};\n Show[Graphics[AbsoluteThickness[2]],\n \ Graphics[RGBColor[0,0,0]],\n pmem, pminmax, \n \ Background->RGBColor[1,1,0],\n AspectRatio->aratio,PlotLabel->title ];\n \ ];\n \nPlotBeamShape[ncoor_,ue_,amp_,nsub_]:=Module[\n \ {x1,y1,x2,y2,x21,y21,L,ux1,uy1,theta1,ux2,uy2,theta2,\n \ uxbar1,uybar1,uxbar2,uybar2,k,xi,x,y,x0,y0,\n \ uxbar,uybar,xmin,xmax,ymin,ymax,xold,yold,p={}},\n {{x1,y1},{x2,y2}}=ncoor; \ x21=x2-x1; y21=y2-y1; \n {ux1,uy1,theta1,ux2,uy2,theta2}=amp*ue;\n \ L=N[Sqrt[x21^2+y21^2]];\n uxbar1= (x21*ux1+y21*uy1)/L; uxbar2= \ (x21*ux2+y21*uy2)/L; \n uybar1=(-y21*ux1+x21*uy1)/L; \ uybar2=(-y21*ux2+x21*uy2)/L; ", StyleBox[" ", FontColor->RGBColor[1, 0, 0]], "\n {xold,yold}={x1+ux1,y1+uy1}; \n For [k=1, k<=nsub, k++,\n \ xi=N[(2*k-nsub)/nsub];\n x0=0.5*(x1+x2+x21*xi); y0=0.5*(y1+y2+y21*xi);\n \ uxbar=0.5*(uxbar1+uxbar2+(uxbar2-uxbar1)*xi);\n \ uybar=0.125*(4*(uybar1+uybar2)+2*(uybar1-uybar2)*(xi^2-3)*\n \ xi+L*(xi^2-1)*(theta2-theta1+(theta1+theta2)*xi));\n \ x=x0+(uxbar*x21-uybar*y21)/L; y=y0+(uxbar*y21+uybar*x21)/L; \n \ AppendTo[p,Graphics[Line[{{xold,yold},{x,y}}]]];\n xold=x; yold=y;\n \ ]; Return[p];\n ];\n \n\ BeamShapePlotBounds[ncoor_,ue_,ampmax_,nsub_]:=Module[\n \ {x1,y1,x2,y2,x21,y21,L,ux1,uy1,theta1,ux2,uy2,theta2,\n \ uxbar1,uybar1,uxbar2,uybar2,k,xi,ux,uy,x0,y0,\n \ uxbar,uybar,xmin,xmax,ymin,ymax},\n {{x1,y1},{x2,y2}}=ncoor; x21=x2-x1; \ y21=y2-y1; \n {ux1,uy1,theta1,ux2,uy2,theta2}=ampmax*ue;\n \ L=N[Sqrt[x21^2+y21^2]];\n uxbar1= (x21*ux1+y21*uy1)/L; uxbar2= \ (x21*ux2+y21*uy2)/L; \n uybar1=(-y21*ux1+x21*uy1)/L; \ uybar2=(-y21*ux2+x21*uy2)/L; ", StyleBox[" ", FontColor->RGBColor[1, 0, 0]], "\n xmin=xmax=x1+ux1; ymin=ymax=y1+uy1;\n For [k=1, k<=nsub, k++,\n \ xi=N[(2*k-nsub)/nsub];\n x0=0.5*(x1+x2+x21*xi); y0=0.5*(y1+y2+y21*xi);\n \ uxbar=0.5*(uxbar1+uxbar2+(uxbar2-uxbar1)*xi);\n \ uybar=0.125*(4*(uybar1+uybar2)+2*(uybar1-uybar2)*(xi^2-3)*\n \ xi+L*(xi^2-1)*(theta2-theta1+(theta1+theta2)*xi));\n \ ux=(uxbar*x21-uybar*y21)/L; uy=(uxbar*y21+uybar*x21)/L; \n \ xmin=Min[xmin,x0+ux,x0-ux]; ymin=Max[ymin,y0+uy,y0-uy];\n \ xmax=Max[xmax,x0+ux,x0-ux]; ymax=Max[ymax,y0+uy,y0-uy];\n ]; \ Return[{xmin,xmax,ymin,ymax}];\n ]; \n", StyleBox[" \nClearAll[L,H];\nL=10000; H=6000; \n\ NodeCoordinates={{-L,0},{-L,H/2},{-L,H},{-L/2,H},{0,H},{L/2,H},\n \ {L,H},{L,H/2},{L,0},{0,H/2},{0,0}};\nElemNodeLists= \ {{1,2},{2,3},{3,4},{4,5},{5,6},{6,7},{7,8},\n \ {8,9},{5,10},{10,11}};\nnumnod=Length[NodeCoordinates]; \ numele=Length[ElemNodeLists];\naspect=0; \nu=Table[0,{3*numnod}]; u[[6]]=1; \ u[[9]]=-1; u[[12]]=1;\nu[[15]]=-1; u[[18]]=1; u[[21]]=-1; u[[24]]=1; \ u[[27]]=-1;\nu[[33]]=2;\npcont={0.5,1.0,0.0,8};\n\n\ PlotPlaneFrameDeformedShape[NodeCoordinates,ElemNodeLists,\n \ u,pcont,aspect,\"a crazy deformed shape\"];", FontColor->RGBColor[1, 0, 0]] }], "Input", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0.699214, 1, 0.0500191]], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .33721 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Background color 1 1 0 r MFill % Scaling calculations 0.5 4.42968e-05 0.00802879 4.42968e-05 [ [.5 .34971 -68 0 ] [.5 .34971 68 12 ] [ 0 0 0 0 ] [ 1 .33721 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 0 1 r gsave .5 .34971 -129 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 0.000 0.000 rmoveto 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor (a) show 75.000 13.000 moveto (crazy) show 111.000 13.000 moveto (deformed) show 165.000 13.000 moveto (shape) show 195.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore 0 0 m 1 0 L 1 .33721 L 0 .33721 L closepath clip newpath 0 0 0 r 2 Mabswid [ ] 0 setdash .05703 .00803 m .05794 .02464 L s .05794 .02464 m .06015 .04125 L s .06015 .04125 m .06287 .05786 L s .06287 .05786 m .06534 .07447 L s .06534 .07447 m .06677 .09109 L s .06677 .09109 m .06638 .1077 L s .06638 .1077 m .06339 .12431 L s .06339 .12431 m .05703 .14092 L s .05703 .14092 m .04976 .15753 L s .04976 .15753 m .04457 .17414 L s .04457 .17414 m .04146 .19075 L s .04146 .19075 m .04042 .20736 L s .04042 .20736 m .04146 .22398 L s .04146 .22398 m .04457 .24059 L s .04457 .24059 m .04976 .2572 L s .04976 .2572 m .05703 .27381 L s .05703 .27381 m .08472 .2617 L s .08472 .2617 m .1124 .25305 L s .1124 .25305 m .14009 .24785 L s .14009 .24785 m .16777 .24612 L s .16777 .24612 m .19546 .24785 L s .19546 .24785 m .22315 .25305 L s .22315 .25305 m .25083 .2617 L s .25083 .2617 m .27852 .27381 L s .27852 .27381 m .3062 .28592 L s .3062 .28592 m .33389 .29457 L s .33389 .29457 m .36157 .29976 L s .36157 .29976 m .38926 .3015 L s .38926 .3015 m .41694 .29976 L s .41694 .29976 m .44463 .29457 L s .44463 .29457 m .47231 .28592 L s .47231 .28592 m .5 .27381 L s .5 .27381 m .52769 .2617 L s .52769 .2617 m .55537 .25305 L s .55537 .25305 m .58306 .24785 L s .58306 .24785 m .61074 .24612 L s .61074 .24612 m .63843 .24785 L s .63843 .24785 m .66611 .25305 L s .66611 .25305 m .6938 .2617 L s .6938 .2617 m .72148 .27381 L s .72148 .27381 m .74917 .28592 L s .74917 .28592 m .77685 .29457 L s .77685 .29457 m .80454 .29976 L s .80454 .29976 m .83223 .3015 L s .83223 .3015 m .85991 .29976 L s .85991 .29976 m .8876 .29457 L s .8876 .29457 m .91528 .28592 L s .91528 .28592 m .94297 .27381 L s .94297 .27381 m .9357 .2572 L s .9357 .2572 m .93051 .24059 L s .93051 .24059 m .92739 .22398 L s .92739 .22398 m .92636 .20736 L s .92636 .20736 m .92739 .19075 L s .92739 .19075 m .93051 .17414 L s .93051 .17414 m .9357 .15753 L s .9357 .15753 m .94297 .14092 L s .94297 .14092 m .95024 .12431 L s .95024 .12431 m .95543 .1077 L s .95543 .1077 m .95854 .09109 L s .95854 .09109 m .95958 .07447 L s .95958 .07447 m .95854 .05786 L s .95854 .05786 m .95543 .04125 L s .95543 .04125 m .95024 .02464 L s .95024 .02464 m .94297 .00803 L s .5 .27381 m .49364 .2572 L s .49364 .2572 m .49066 .24059 L s .49066 .24059 m .49027 .22398 L s .49027 .22398 m .49169 .20736 L s .49169 .20736 m .49416 .19075 L s .49416 .19075 m .49689 .17414 L s .49689 .17414 m .49909 .15753 L s .49909 .15753 m .5 .14092 L s .5 .14092 m .49818 .12431 L s .49818 .12431 m .49377 .1077 L s .49377 .1077 m .48832 .09109 L s .48832 .09109 m .48339 .07447 L s .48339 .07447 m .48053 .05786 L s .48053 .05786 m .48131 .04125 L s .48131 .04125 m .48728 .02464 L s .48728 .02464 m .5 .00803 L s 1 1 1 r 1 Mabswid .02381 .00803 Mdot .97619 .32918 Mdot % End of Graphics MathPictureEnd \ \>"], "Graphics", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageSize->{472.438, 159.312}, ImageMargins->{{0, 0}, {0, 0}}, ImageRegion->{{-0, 1}, {0, 1}}, ImageCacheValid->False] }, Open ]], Cell[CellGroupData[{ Cell[TextData[{ "PlotPlaneFrameStress[nodcoor_,elenod_,elefab_,p_,fmax_,Nsub_,\n \ aspect_,title_]:= Module[\n \ {eNL,e,i,ii,n,ni,nj,nc,xyc,psig={},numele=Length[elenod],\n \ numnod=Length[nodcoor],x,y,pe,\n \ xmin,xmax,ymin,ymax,xam,xax,yam,yax,xbm,xbx,ybm,ybx,\n ns,aratio},\n \ x=Table[nodcoor[[i,1]],{i,1,numnod}]; \n \ y=Table[nodcoor[[i,2]],{i,1,numnod}];\n \ {xmin,xmax,ymin,ymax}=N[{Min[x],Max[x],Min[y],Max[y]}];\n For \ [e=1,e<=numele,e++, eNL=elenod[[e]]; nc=Length[eNL]; \n If [nc!=2, \ Continue[]]; {ni,nj}=eNL;\n fprop=elefab[[e]]; pe=p[[e]];\n \ xyc={{x[[ni]],y[[ni]]},{x[[nj]],y[[nj]]}}; \n \ eftab={3*ni-2,3*ni-1,3*ni,3*nj-2,3*nj-1,3*nj};\n \ psig=Append[psig,PlotBeamColumn2Stress[xyc,fprop,pe,fmax,Nsub]];\n ];\n \ If [aspect>0, aratio=aspect, aratio=(ymax-ymin)/(xmax-xmin)];\n If [aspect<0, \ aratio=Automatic];\n Show[Graphics[AbsoluteThickness[2]],\n psig,\n \ Background->RGBColor[1,1,0],\n AspectRatio->aratio,PlotLabel->title ];\n \ ];\n\nPlotBeamColumn2Stress[ncoor_,fprop_,pe_,fmax_,Nsub_]:=Module[\n \ {x1,y1,x2,y2,x21,y21,L,c,s,A,Izz,H,F,m1,m2,Hx,Hy,\n xyc,fc,p1,p2,psig={}},\n\ {{x1,y1},{x2,y2}}=ncoor; x21=x2-x1; y21=y2-y1;\n L=N[Sqrt[x21^2+y21^2]]; \ c=x21/L; s=y21/L; {A,Izz,H}=fprop;\n {F,m1,m2}=pe; Hx=H*s/2; Hy=H*c/2; \n \ xyc={{x1+Hx,y1-Hy},{x2+Hx,y2-Hy},{x2-Hx,y2+Hy},{x1-Hx,y1+Hy}};\n \ xyc=N[xyc]; p1=-m1*H/(2*Izz); p2=-m2*H/(2*Izz);\n \ fc=N[{F/A-p1,F/A-p2,F/A+p2,F/A+p1}];\n \ psig=PlotFunctionOverQuad[xyc,fc,fmax,Nsub];\n Return[psig];\n ];\n \n\ PlotFunctionOverQuad[xyc_,fc_,fmax_,Nsub_]:=Module[\n \ {Ne,Nev,xy1,xy2,xy3,i,j,ixi,ieta,xi,eta,nxi,neta,x1,x2,\n \ x3,x4,y1,y2,y3,y4,xc,yc,c1,c2,c3,d,f1,f2,f3,f4,f,poly={}},\n \ {{x1,y1},{x2,y2},{x3,y3},{x4,y4}}=xyc;\n xc={x1,x2,x3,x4}; \ yc={y1,y2,y3,y4};{f1,f2,f3,f4}=fc;\n \ Ne[xi_,eta_]:=N[{(1-xi)*(1-eta),(1+xi)*(1-eta),\n \ (1+xi)*(1+eta),(1-xi)*(1+eta)}/4];\n If [Length[Nsub]==0, nxi=neta= Nsub];\n\ If [Length[Nsub]>0, {nxi,neta}=Nsub];", StyleBox[" ", FontColor->RGBColor[1, 0, 0]], "\n Do [Do [ixi=(2*i-nxi-1)/nxi; ieta=(2*j-neta-1)/neta;\n \ {xi,eta}=N[{ixi-1/nxi,ieta-1/neta}];Nev=Ne[xi,eta];\n \ xy1={xc.Nev,yc.Nev};\n \ {xi,eta}=N[{ixi+1/nxi,ieta-1/neta}];Nev=Ne[xi,eta];\n \ xy2={xc.Nev,yc.Nev};\n \ {xi,eta}=N[{ixi+1/nxi,ieta+1/neta}];Nev=Ne[xi,eta];\n \ xy3={xc.Nev,yc.Nev};\n \ {xi,eta}=N[{ixi-1/nxi,ieta+1/neta}];Nev=Ne[xi,eta];\n \ xy4={xc.Nev,yc.Nev};\n Nev=Ne[N[ixi],N[ieta]]; \n \ {c1,c2,c3}=ContourPolyColor[fc.Nev,fmax];\n \ AppendTo[poly,Graphics[RGBColor[c1,c2,c3]]];\n \ AppendTo[poly,Graphics[Polygon[{xy1,xy2,xy3,xy4}]]],\n \ {j,1,neta}],{i,1,nxi}];\n Return[poly];\n ]; \n\n\ ContourPolyColor[f_,fmax_]:= Module[{r,RGBmax={1,0,0}, \n RGBmin={0,0,1}, \ RGBzero={1,1,1}, RGBout={0,0,0}},\n If [f==0 || fmax==0, \n \ Return[RGBzero]]; (* White if f=0 *)\n If [f>fmax || f<-fmax, \n \ Return[RGBout ]]; (* Black if outside range *)\n If [f>0, r= N[f/fmax]; \n \ Return[r*RGBmax+(1-r)*RGBzero]]; (* positive *)\n If [f<0, \ r=-N[f/fmax]; \n Return[r*RGBmin+(1-r)*RGBzero]]; (* negative *)\n]; \n\ \n", StyleBox["(*p=PlotBeamColumn2Stress[{{0,0},{10,0}},{1,6,2},{0,30,-30},5,{8,\ 4}];\nShow[p,AspectRatio->Automatic]; *)\n\nnodcoor={{0,0},{10,0},{10,10}};\n\ elenod={{1,2},{2,3},{1,3}};\nelefab=Table[{2,4,1},{3}];\n\ p={{0,8,-8},{0,12,-12},{0,10,-10}};\nfmax=1.5; Nsub={32,6}; aspect=-1; \n\ title=\"Axial stress in frame members\";\n\ PlotPlaneFrameStress[nodcoor,elenod,elefab,p,fmax,Nsub,\n \ aspect,title];", FontColor->RGBColor[1, 0, 0]], " \n" }], "Input", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, InitializationCell->True, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0.699214, 1, 0.0500191]], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: 1 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Background color 1 1 0 r MFill % Scaling calculations 0.0548332 0.0877483 0.0676837 0.0877483 [ [.5 1.0125 -89 0 ] [.5 1.0125 89 12 ] [ 0 0 0 0 ] [ 1 1 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 0 1 r gsave .5 1.0125 -150 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale /g { setgray} bind def /k { setcmykcolor} bind def /p { gsave} bind def /r { setrgbcolor} bind def /w { setlinewidth} bind def /C { curveto} bind def /F { fill} bind def /L { lineto} bind def /rL { rlineto} bind def /P { grestore} bind def /s { stroke} bind def /S { show} bind def /N {currentpoint 3 -1 roll show moveto} bind def /Msf { findfont exch scalefont [1 0 0 -1 0 0 ] makefont setfont} bind def /m { moveto} bind def /Mr { rmoveto} bind def /Mx {currentpoint exch pop moveto} bind def /My {currentpoint pop exch moveto} bind def /X {0 rmoveto} bind def /Y {0 exch rmoveto} bind def 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 0.000 0.000 rmoveto 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor (Axial) show 99.000 13.000 moveto (stress) show 141.000 13.000 moveto (in) show 159.000 13.000 moveto (frame) show 195.000 13.000 moveto (members) show 237.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore 0 0 m 1 0 L 1 1 L 0 1 L closepath clip newpath 1 .462 .462 r .05483 .02381 m .08225 .02381 L .08225 .03843 L .05483 .03843 L F 1 .677 .677 r .05483 .03843 m .08225 .03843 L .08225 .05306 L .05483 .05306 L F 1 .892 .892 r .05483 .05306 m .08225 .05306 L .08225 .06768 L .05483 .06768 L F .892 .892 1 r .05483 .06768 m .08225 .06768 L .08225 .08231 L .05483 .08231 L F .677 .677 1 r .05483 .08231 m .08225 .08231 L .08225 .09693 L .05483 .09693 L F .462 .462 1 r .05483 .09693 m .08225 .09693 L .08225 .11156 L .05483 .11156 L F 1 .497 .497 r .08225 .02381 m .10968 .02381 L .10968 .03843 L .08225 .03843 L F 1 .698 .698 r .08225 .03843 m .10968 .03843 L .10968 .05306 L .08225 .05306 L F 1 .899 .899 r .08225 .05306 m .10968 .05306 L .10968 .06768 L .08225 .06768 L F .899 .899 1 r .08225 .06768 m .10968 .06768 L .10968 .08231 L .08225 .08231 L F .698 .698 1 r .08225 .08231 m .10968 .08231 L .10968 .09693 L .08225 .09693 L F .497 .497 1 r .08225 .09693 m .10968 .09693 L .10968 .11156 L .08225 .11156 L F 1 .531 .531 r .10968 .02381 m .1371 .02381 L .1371 .03843 L .10968 .03843 L F 1 .719 .719 r .10968 .03843 m .1371 .03843 L .1371 .05306 L .10968 .05306 L F 1 .906 .906 r .10968 .05306 m .1371 .05306 L .1371 .06768 L .10968 .06768 L F .906 .906 1 r .10968 .06768 m .1371 .06768 L .1371 .08231 L .10968 .08231 L F .719 .719 1 r .10968 .08231 m .1371 .08231 L .1371 .09693 L .10968 .09693 L F .531 .531 1 r .10968 .09693 m .1371 .09693 L .1371 .11156 L .10968 .11156 L F 1 .566 .566 r .1371 .02381 m .16452 .02381 L .16452 .03843 L .1371 .03843 L F 1 .74 .74 r .1371 .03843 m .16452 .03843 L .16452 .05306 L .1371 .05306 L F 1 .913 .913 r .1371 .05306 m .16452 .05306 L .16452 .06768 L .1371 .06768 L F .913 .913 1 r .1371 .06768 m .16452 .06768 L .16452 .08231 L .1371 .08231 L F .74 .74 1 r .1371 .08231 m .16452 .08231 L .16452 .09693 L .1371 .09693 L F .566 .566 1 r .1371 .09693 m .16452 .09693 L .16452 .11156 L .1371 .11156 L F 1 .601 .601 r .16452 .02381 m .19194 .02381 L .19194 .03843 L .16452 .03843 L F 1 .76 .76 r .16452 .03843 m .19194 .03843 L .19194 .05306 L .16452 .05306 L F 1 .92 .92 r .16452 .05306 m .19194 .05306 L .19194 .06768 L .16452 .06768 L F .92 .92 1 r .16452 .06768 m .19194 .06768 L .19194 .08231 L .16452 .08231 L F .76 .76 1 r .16452 .08231 m .19194 .08231 L .19194 .09693 L .16452 .09693 L F .601 .601 1 r .16452 .09693 m .19194 .09693 L .19194 .11156 L .16452 .11156 L F 1 .635 .635 r .19194 .02381 m .21936 .02381 L .21936 .03843 L .19194 .03843 L F 1 .781 .781 r .19194 .03843 m .21936 .03843 L .21936 .05306 L .19194 .05306 L F 1 .927 .927 r .19194 .05306 m .21936 .05306 L .21936 .06768 L .19194 .06768 L F .927 .927 1 r .19194 .06768 m .21936 .06768 L .21936 .08231 L .19194 .08231 L F .781 .781 1 r .19194 .08231 m .21936 .08231 L .21936 .09693 L .19194 .09693 L F .635 .635 1 r .19194 .09693 m .21936 .09693 L .21936 .11156 L .19194 .11156 L F 1 .67 .67 r .21936 .02381 m .24678 .02381 L .24678 .03843 L .21936 .03843 L F 1 .802 .802 r .21936 .03843 m .24678 .03843 L .24678 .05306 L .21936 .05306 L F 1 .934 .934 r .21936 .05306 m .24678 .05306 L .24678 .06768 L .21936 .06768 L F .934 .934 1 r .21936 .06768 m .24678 .06768 L .24678 .08231 L .21936 .08231 L F .802 .802 1 r .21936 .08231 m .24678 .08231 L .24678 .09693 L .21936 .09693 L F .67 .67 1 r .21936 .09693 m .24678 .09693 L .24678 .11156 L .21936 .11156 L F 1 .705 .705 r .24678 .02381 m .2742 .02381 L .2742 .03843 L .24678 .03843 L F 1 .823 .823 r .24678 .03843 m .2742 .03843 L .2742 .05306 L .24678 .05306 L F 1 .941 .941 r .24678 .05306 m .2742 .05306 L .2742 .06768 L .24678 .06768 L F .941 .941 1 r .24678 .06768 m .2742 .06768 L .2742 .08231 L .24678 .08231 L F .823 .823 1 r .24678 .08231 m .2742 .08231 L .2742 .09693 L .24678 .09693 L F .705 .705 1 r .24678 .09693 m .2742 .09693 L .2742 .11156 L .24678 .11156 L F 1 .74 .74 r .2742 .02381 m .30163 .02381 L .30163 .03843 L .2742 .03843 L F 1 .844 .844 r .2742 .03843 m .30163 .03843 L .30163 .05306 L .2742 .05306 L F 1 .948 .948 r .2742 .05306 m .30163 .05306 L .30163 .06768 L .2742 .06768 L F .948 .948 1 r .2742 .06768 m .30163 .06768 L .30163 .08231 L .2742 .08231 L F .844 .844 1 r .2742 .08231 m .30163 .08231 L .30163 .09693 L .2742 .09693 L F .74 .74 1 r .2742 .09693 m .30163 .09693 L .30163 .11156 L .2742 .11156 L F 1 .774 .774 r .30163 .02381 m .32905 .02381 L .32905 .03843 L .30163 .03843 L F 1 .865 .865 r .30163 .03843 m .32905 .03843 L .32905 .05306 L .30163 .05306 L F 1 .955 .955 r .30163 .05306 m .32905 .05306 L .32905 .06768 L .30163 .06768 L F .955 .955 1 r .30163 .06768 m .32905 .06768 L .32905 .08231 L .30163 .08231 L F .865 .865 1 r .30163 .08231 m .32905 .08231 L .32905 .09693 L .30163 .09693 L F .774 .774 1 r .30163 .09693 m .32905 .09693 L .32905 .11156 L .30163 .11156 L F 1 .809 .809 r .32905 .02381 m .35647 .02381 L .35647 .03843 L .32905 .03843 L F 1 .885 .885 r .32905 .03843 m .35647 .03843 L .35647 .05306 L .32905 .05306 L F 1 .962 .962 r .32905 .05306 m .35647 .05306 L .35647 .06768 L .32905 .06768 L F .962 .962 1 r .32905 .06768 m .35647 .06768 L .35647 .08231 L .32905 .08231 L F .885 .885 1 r .32905 .08231 m .35647 .08231 L .35647 .09693 L .32905 .09693 L F .809 .809 1 r .32905 .09693 m .35647 .09693 L .35647 .11156 L .32905 .11156 L F 1 .844 .844 r .35647 .02381 m .38389 .02381 L .38389 .03843 L .35647 .03843 L F 1 .906 .906 r .35647 .03843 m .38389 .03843 L .38389 .05306 L .35647 .05306 L F 1 .969 .969 r .35647 .05306 m .38389 .05306 L .38389 .06768 L .35647 .06768 L F .969 .969 1 r .35647 .06768 m .38389 .06768 L .38389 .08231 L .35647 .08231 L F .906 .906 1 r .35647 .08231 m .38389 .08231 L .38389 .09693 L .35647 .09693 L F .844 .844 1 r .35647 .09693 m .38389 .09693 L .38389 .11156 L .35647 .11156 L F 1 .878 .878 r .38389 .02381 m .41131 .02381 L .41131 .03843 L .38389 .03843 L F 1 .927 .927 r .38389 .03843 m .41131 .03843 L .41131 .05306 L .38389 .05306 L F 1 .976 .976 r .38389 .05306 m .41131 .05306 L .41131 .06768 L .38389 .06768 L F .976 .976 1 r .38389 .06768 m .41131 .06768 L .41131 .08231 L .38389 .08231 L F .927 .927 1 r .38389 .08231 m .41131 .08231 L .41131 .09693 L .38389 .09693 L F .878 .878 1 r .38389 .09693 m .41131 .09693 L .41131 .11156 L .38389 .11156 L F 1 .913 .913 r .41131 .02381 m .43873 .02381 L .43873 .03843 L .41131 .03843 L F 1 .948 .948 r .41131 .03843 m .43873 .03843 L .43873 .05306 L .41131 .05306 L F 1 .983 .983 r .41131 .05306 m .43873 .05306 L .43873 .06768 L .41131 .06768 L F .983 .983 1 r .41131 .06768 m .43873 .06768 L .43873 .08231 L .41131 .08231 L F .948 .948 1 r .41131 .08231 m .43873 .08231 L .43873 .09693 L .41131 .09693 L F .913 .913 1 r .41131 .09693 m .43873 .09693 L .43873 .11156 L .41131 .11156 L F 1 .948 .948 r .43873 .02381 m .46615 .02381 L .46615 .03843 L .43873 .03843 L F 1 .969 .969 r .43873 .03843 m .46615 .03843 L .46615 .05306 L .43873 .05306 L F 1 .99 .99 r .43873 .05306 m .46615 .05306 L .46615 .06768 L .43873 .06768 L F .99 .99 1 r .43873 .06768 m .46615 .06768 L .46615 .08231 L .43873 .08231 L F .969 .969 1 r .43873 .08231 m .46615 .08231 L .46615 .09693 L .43873 .09693 L F .948 .948 1 r .43873 .09693 m .46615 .09693 L .46615 .11156 L .43873 .11156 L F 1 .983 .983 r .46615 .02381 m .49357 .02381 L .49357 .03843 L .46615 .03843 L F 1 .99 .99 r .46615 .03843 m .49357 .03843 L .49357 .05306 L .46615 .05306 L F 1 .997 .997 r .46615 .05306 m .49357 .05306 L .49357 .06768 L .46615 .06768 L F .997 .997 1 r .46615 .06768 m .49357 .06768 L .49357 .08231 L .46615 .08231 L F .99 .99 1 r .46615 .08231 m .49357 .08231 L .49357 .09693 L .46615 .09693 L F .983 .983 1 r .46615 .09693 m .49357 .09693 L .49357 .11156 L .46615 .11156 L F .49357 .02381 m .521 .02381 L .521 .03843 L .49357 .03843 L F .99 .99 1 r .49357 .03843 m .521 .03843 L .521 .05306 L .49357 .05306 L F .997 .997 1 r .49357 .05306 m .521 .05306 L .521 .06768 L .49357 .06768 L F 1 .997 .997 r .49357 .06768 m .521 .06768 L .521 .08231 L .49357 .08231 L F 1 .99 .99 r .49357 .08231 m .521 .08231 L .521 .09693 L .49357 .09693 L F 1 .983 .983 r .49357 .09693 m .521 .09693 L .521 .11156 L .49357 .11156 L F .948 .948 1 r .521 .02381 m .54842 .02381 L .54842 .03843 L .521 .03843 L F .969 .969 1 r .521 .03843 m .54842 .03843 L .54842 .05306 L .521 .05306 L F .99 .99 1 r .521 .05306 m .54842 .05306 L .54842 .06768 L .521 .06768 L F 1 .99 .99 r .521 .06768 m .54842 .06768 L .54842 .08231 L .521 .08231 L F 1 .969 .969 r .521 .08231 m .54842 .08231 L .54842 .09693 L .521 .09693 L F 1 .948 .948 r .521 .09693 m .54842 .09693 L .54842 .11156 L .521 .11156 L F .913 .913 1 r .54842 .02381 m .57584 .02381 L .57584 .03843 L .54842 .03843 L F .948 .948 1 r .54842 .03843 m .57584 .03843 L .57584 .05306 L .54842 .05306 L F .983 .983 1 r .54842 .05306 m .57584 .05306 L .57584 .06768 L .54842 .06768 L F 1 .983 .983 r .54842 .06768 m .57584 .06768 L .57584 .08231 L .54842 .08231 L F 1 .948 .948 r .54842 .08231 m .57584 .08231 L .57584 .09693 L .54842 .09693 L F 1 .913 .913 r .54842 .09693 m .57584 .09693 L .57584 .11156 L .54842 .11156 L F .878 .878 1 r .57584 .02381 m .60326 .02381 L .60326 .03843 L .57584 .03843 L F .927 .927 1 r .57584 .03843 m .60326 .03843 L .60326 .05306 L .57584 .05306 L F .976 .976 1 r .57584 .05306 m .60326 .05306 L .60326 .06768 L .57584 .06768 L F 1 .976 .976 r .57584 .06768 m .60326 .06768 L .60326 .08231 L .57584 .08231 L F 1 .927 .927 r .57584 .08231 m .60326 .08231 L .60326 .09693 L .57584 .09693 L F 1 .878 .878 r .57584 .09693 m .60326 .09693 L .60326 .11156 L .57584 .11156 L F .844 .844 1 r .60326 .02381 m .63068 .02381 L .63068 .03843 L .60326 .03843 L F .906 .906 1 r .60326 .03843 m .63068 .03843 L .63068 .05306 L .60326 .05306 L F .969 .969 1 r .60326 .05306 m .63068 .05306 L .63068 .06768 L .60326 .06768 L F 1 .969 .969 r .60326 .06768 m .63068 .06768 L .63068 .08231 L .60326 .08231 L F 1 .906 .906 r .60326 .08231 m .63068 .08231 L .63068 .09693 L .60326 .09693 L F 1 .844 .844 r .60326 .09693 m .63068 .09693 L .63068 .11156 L .60326 .11156 L F .809 .809 1 r .63068 .02381 m .6581 .02381 L .6581 .03843 L .63068 .03843 L F .885 .885 1 r .63068 .03843 m .6581 .03843 L .6581 .05306 L .63068 .05306 L F .962 .962 1 r .63068 .05306 m .6581 .05306 L .6581 .06768 L .63068 .06768 L F 1 .962 .962 r .63068 .06768 m .6581 .06768 L .6581 .08231 L .63068 .08231 L F 1 .885 .885 r .63068 .08231 m .6581 .08231 L .6581 .09693 L .63068 .09693 L F 1 .809 .809 r .63068 .09693 m .6581 .09693 L .6581 .11156 L .63068 .11156 L F .774 .774 1 r .6581 .02381 m .68552 .02381 L .68552 .03843 L .6581 .03843 L F .865 .865 1 r .6581 .03843 m .68552 .03843 L .68552 .05306 L .6581 .05306 L F .955 .955 1 r .6581 .05306 m .68552 .05306 L .68552 .06768 L .6581 .06768 L F 1 .955 .955 r .6581 .06768 m .68552 .06768 L .68552 .08231 L .6581 .08231 L F 1 .865 .865 r .6581 .08231 m .68552 .08231 L .68552 .09693 L .6581 .09693 L F 1 .774 .774 r .6581 .09693 m .68552 .09693 L .68552 .11156 L .6581 .11156 L F .74 .74 1 r .68552 .02381 m .71295 .02381 L .71295 .03843 L .68552 .03843 L F .844 .844 1 r .68552 .03843 m .71295 .03843 L .71295 .05306 L .68552 .05306 L F .948 .948 1 r .68552 .05306 m .71295 .05306 L .71295 .06768 L .68552 .06768 L F 1 .948 .948 r .68552 .06768 m .71295 .06768 L .71295 .08231 L .68552 .08231 L F 1 .844 .844 r .68552 .08231 m .71295 .08231 L .71295 .09693 L .68552 .09693 L F 1 .74 .74 r .68552 .09693 m .71295 .09693 L .71295 .11156 L .68552 .11156 L F .705 .705 1 r .71295 .02381 m .74037 .02381 L .74037 .03843 L .71295 .03843 L F .823 .823 1 r .71295 .03843 m .74037 .03843 L .74037 .05306 L .71295 .05306 L F .941 .941 1 r .71295 .05306 m .74037 .05306 L .74037 .06768 L .71295 .06768 L F 1 .941 .941 r .71295 .06768 m .74037 .06768 L .74037 .08231 L .71295 .08231 L F 1 .823 .823 r .71295 .08231 m .74037 .08231 L .74037 .09693 L .71295 .09693 L F 1 .705 .705 r .71295 .09693 m .74037 .09693 L .74037 .11156 L .71295 .11156 L F .67 .67 1 r .74037 .02381 m .76779 .02381 L .76779 .03843 L .74037 .03843 L F .802 .802 1 r .74037 .03843 m .76779 .03843 L .76779 .05306 L .74037 .05306 L F .934 .934 1 r .74037 .05306 m .76779 .05306 L .76779 .06768 L .74037 .06768 L F 1 .934 .934 r .74037 .06768 m .76779 .06768 L .76779 .08231 L .74037 .08231 L F 1 .802 .802 r .74037 .08231 m .76779 .08231 L .76779 .09693 L .74037 .09693 L F 1 .67 .67 r .74037 .09693 m .76779 .09693 L .76779 .11156 L .74037 .11156 L F .635 .635 1 r .76779 .02381 m .79521 .02381 L .79521 .03843 L .76779 .03843 L F .781 .781 1 r .76779 .03843 m .79521 .03843 L .79521 .05306 L .76779 .05306 L F .927 .927 1 r .76779 .05306 m .79521 .05306 L .79521 .06768 L .76779 .06768 L F 1 .927 .927 r .76779 .06768 m .79521 .06768 L .79521 .08231 L .76779 .08231 L F 1 .781 .781 r .76779 .08231 m .79521 .08231 L .79521 .09693 L .76779 .09693 L F 1 .635 .635 r .76779 .09693 m .79521 .09693 L .79521 .11156 L .76779 .11156 L F .601 .601 1 r .79521 .02381 m .82263 .02381 L .82263 .03843 L .79521 .03843 L F .76 .76 1 r .79521 .03843 m .82263 .03843 L .82263 .05306 L .79521 .05306 L F .92 .92 1 r .79521 .05306 m .82263 .05306 L .82263 .06768 L .79521 .06768 L F 1 .92 .92 r .79521 .06768 m .82263 .06768 L .82263 .08231 L .79521 .08231 L F 1 .76 .76 r .79521 .08231 m .82263 .08231 L .82263 .09693 L .79521 .09693 L F 1 .601 .601 r .79521 .09693 m .82263 .09693 L .82263 .11156 L .79521 .11156 L F .566 .566 1 r .82263 .02381 m .85005 .02381 L .85005 .03843 L .82263 .03843 L F .74 .74 1 r .82263 .03843 m .85005 .03843 L .85005 .05306 L .82263 .05306 L F .913 .913 1 r .82263 .05306 m .85005 .05306 L .85005 .06768 L .82263 .06768 L F 1 .913 .913 r .82263 .06768 m .85005 .06768 L .85005 .08231 L .82263 .08231 L F 1 .74 .74 r .82263 .08231 m .85005 .08231 L .85005 .09693 L .82263 .09693 L F 1 .566 .566 r .82263 .09693 m .85005 .09693 L .85005 .11156 L .82263 .11156 L F .531 .531 1 r .85005 .02381 m .87747 .02381 L .87747 .03843 L .85005 .03843 L F .719 .719 1 r .85005 .03843 m .87747 .03843 L .87747 .05306 L .85005 .05306 L F .906 .906 1 r .85005 .05306 m .87747 .05306 L .87747 .06768 L .85005 .06768 L F 1 .906 .906 r .85005 .06768 m .87747 .06768 L .87747 .08231 L .85005 .08231 L F 1 .719 .719 r .85005 .08231 m .87747 .08231 L .87747 .09693 L .85005 .09693 L F 1 .531 .531 r .85005 .09693 m .87747 .09693 L .87747 .11156 L .85005 .11156 L F .497 .497 1 r .87747 .02381 m .90489 .02381 L .90489 .03843 L .87747 .03843 L F .698 .698 1 r .87747 .03843 m .90489 .03843 L .90489 .05306 L .87747 .05306 L F .899 .899 1 r .87747 .05306 m .90489 .05306 L .90489 .06768 L .87747 .06768 L F 1 .899 .899 r .87747 .06768 m .90489 .06768 L .90489 .08231 L .87747 .08231 L F 1 .698 .698 r .87747 .08231 m .90489 .08231 L .90489 .09693 L .87747 .09693 L F 1 .497 .497 r .87747 .09693 m .90489 .09693 L .90489 .11156 L .87747 .11156 L F .462 .462 1 r .90489 .02381 m .93232 .02381 L .93232 .03843 L .90489 .03843 L F .677 .677 1 r .90489 .03843 m .93232 .03843 L .93232 .05306 L .90489 .05306 L F .892 .892 1 r .90489 .05306 m .93232 .05306 L .93232 .06768 L .90489 .06768 L F 1 .892 .892 r .90489 .06768 m .93232 .06768 L .93232 .08231 L .90489 .08231 L F 1 .677 .677 r .90489 .08231 m .93232 .08231 L .93232 .09693 L .90489 .09693 L F 1 .462 .462 r .90489 .09693 m .93232 .09693 L .93232 .11156 L .90489 .11156 L F 1 .193 .193 r .97619 .06768 m .97619 .09511 L .96157 .09511 L .96157 .06768 L F 1 .516 .516 r .96157 .06768 m .96157 .09511 L .94694 .09511 L .94694 .06768 L F 1 .839 .839 r .94694 .06768 m .94694 .09511 L .93232 .09511 L .93232 .06768 L F .839 .839 1 r .93232 .06768 m .93232 .09511 L .91769 .09511 L .91769 .06768 L F .516 .516 1 r .91769 .06768 m .91769 .09511 L .90307 .09511 L .90307 .06768 L F .193 .193 1 r .90307 .06768 m .90307 .09511 L .88844 .09511 L .88844 .06768 L F 1 .245 .245 r .97619 .09511 m .97619 .12253 L .96157 .12253 L .96157 .09511 L F 1 .547 .547 r .96157 .09511 m .96157 .12253 L .94694 .12253 L .94694 .09511 L F 1 .849 .849 r .94694 .09511 m .94694 .12253 L .93232 .12253 L .93232 .09511 L F .849 .849 1 r .93232 .09511 m .93232 .12253 L .91769 .12253 L .91769 .09511 L F .547 .547 1 r .91769 .09511 m .91769 .12253 L .90307 .12253 L .90307 .09511 L F .245 .245 1 r .90307 .09511 m .90307 .12253 L .88844 .12253 L .88844 .09511 L F 1 .297 .297 r .97619 .12253 m .97619 .14995 L .96157 .14995 L .96157 .12253 L F 1 .578 .578 r .96157 .12253 m .96157 .14995 L .94694 .14995 L .94694 .12253 L F 1 .859 .859 r .94694 .12253 m .94694 .14995 L .93232 .14995 L .93232 .12253 L F .859 .859 1 r .93232 .12253 m .93232 .14995 L .91769 .14995 L .91769 .12253 L F .578 .578 1 r .91769 .12253 m .91769 .14995 L .90307 .14995 L .90307 .12253 L F .297 .297 1 r .90307 .12253 m .90307 .14995 L .88844 .14995 L .88844 .12253 L F 1 .349 .349 r .97619 .14995 m .97619 .17737 L .96157 .17737 L .96157 .14995 L F 1 .609 .609 r .96157 .14995 m .96157 .17737 L .94694 .17737 L .94694 .14995 L F 1 .87 .87 r .94694 .14995 m .94694 .17737 L .93232 .17737 L .93232 .14995 L F .87 .87 1 r .93232 .14995 m .93232 .17737 L .91769 .17737 L .91769 .14995 L F .609 .609 1 r .91769 .14995 m .91769 .17737 L .90307 .17737 L .90307 .14995 L F .349 .349 1 r .90307 .14995 m .90307 .17737 L .88844 .17737 L .88844 .14995 L F 1 .401 .401 r .97619 .17737 m .97619 .20479 L .96157 .20479 L .96157 .17737 L F 1 .641 .641 r .96157 .17737 m .96157 .20479 L .94694 .20479 L .94694 .17737 L F 1 .88 .88 r .94694 .17737 m .94694 .20479 L .93232 .20479 L .93232 .17737 L F .88 .88 1 r .93232 .17737 m .93232 .20479 L .91769 .20479 L .91769 .17737 L F .641 .641 1 r .91769 .17737 m .91769 .20479 L .90307 .20479 L .90307 .17737 L F .401 .401 1 r .90307 .17737 m .90307 .20479 L .88844 .20479 L .88844 .17737 L F 1 .453 .453 r .97619 .20479 m .97619 .23221 L .96157 .23221 L .96157 .20479 L F 1 .672 .672 r .96157 .20479 m .96157 .23221 L .94694 .23221 L .94694 .20479 L F 1 .891 .891 r .94694 .20479 m .94694 .23221 L .93232 .23221 L .93232 .20479 L F .891 .891 1 r .93232 .20479 m .93232 .23221 L .91769 .23221 L .91769 .20479 L F .672 .672 1 r .91769 .20479 m .91769 .23221 L .90307 .23221 L .90307 .20479 L F .453 .453 1 r .90307 .20479 m .90307 .23221 L .88844 .23221 L .88844 .20479 L F 1 .505 .505 r .97619 .23221 m .97619 .25963 L .96157 .25963 L .96157 .23221 L F 1 .703 .703 r .96157 .23221 m .96157 .25963 L .94694 .25963 L .94694 .23221 L F 1 .901 .901 r .94694 .23221 m .94694 .25963 L .93232 .25963 L .93232 .23221 L F .901 .901 1 r .93232 .23221 m .93232 .25963 L .91769 .25963 L .91769 .23221 L F .703 .703 1 r .91769 .23221 m .91769 .25963 L .90307 .25963 L .90307 .23221 L F .505 .505 1 r .90307 .23221 m .90307 .25963 L .88844 .25963 L .88844 .23221 L F 1 .557 .557 r .97619 .25963 m .97619 .28705 L .96157 .28705 L .96157 .25963 L F 1 .734 .734 r .96157 .25963 m .96157 .28705 L .94694 .28705 L .94694 .25963 L F 1 .911 .911 r .94694 .25963 m .94694 .28705 L .93232 .28705 L .93232 .25963 L F .911 .911 1 r .93232 .25963 m .93232 .28705 L .91769 .28705 L .91769 .25963 L F .734 .734 1 r .91769 .25963 m .91769 .28705 L .90307 .28705 L .90307 .25963 L F .557 .557 1 r .90307 .25963 m .90307 .28705 L .88844 .28705 L .88844 .25963 L F 1 .609 .609 r .97619 .28705 m .97619 .31448 L .96157 .31448 L .96157 .28705 L F 1 .766 .766 r .96157 .28705 m .96157 .31448 L .94694 .31448 L .94694 .28705 L F 1 .922 .922 r .94694 .28705 m .94694 .31448 L .93232 .31448 L .93232 .28705 L F .922 .922 1 r .93232 .28705 m .93232 .31448 L .91769 .31448 L .91769 .28705 L F .766 .766 1 r .91769 .28705 m .91769 .31448 L .90307 .31448 L .90307 .28705 L F .609 .609 1 r .90307 .28705 m .90307 .31448 L .88844 .31448 L .88844 .28705 L F 1 .661 .661 r .97619 .31448 m .97619 .3419 L .96157 .3419 L .96157 .31448 L F 1 .797 .797 r .96157 .31448 m .96157 .3419 L .94694 .3419 L .94694 .31448 L F 1 .932 .932 r .94694 .31448 m .94694 .3419 L .93232 .3419 L .93232 .31448 L F .932 .932 1 r .93232 .31448 m .93232 .3419 L .91769 .3419 L .91769 .31448 L F .797 .797 1 r .91769 .31448 m .91769 .3419 L .90307 .3419 L .90307 .31448 L F .661 .661 1 r .90307 .31448 m .90307 .3419 L .88844 .3419 L .88844 .31448 L F 1 .714 .714 r .97619 .3419 m .97619 .36932 L .96157 .36932 L .96157 .3419 L F 1 .828 .828 r .96157 .3419 m .96157 .36932 L .94694 .36932 L .94694 .3419 L F 1 .943 .943 r .94694 .3419 m .94694 .36932 L .93232 .36932 L .93232 .3419 L F .943 .943 1 r .93232 .3419 m .93232 .36932 L .91769 .36932 L .91769 .3419 L F .828 .828 1 r .91769 .3419 m .91769 .36932 L .90307 .36932 L .90307 .3419 L F .714 .714 1 r .90307 .3419 m .90307 .36932 L .88844 .36932 L .88844 .3419 L F 1 .766 .766 r .97619 .36932 m .97619 .39674 L .96157 .39674 L .96157 .36932 L F 1 .859 .859 r .96157 .36932 m .96157 .39674 L .94694 .39674 L .94694 .36932 L F 1 .953 .953 r .94694 .36932 m .94694 .39674 L .93232 .39674 L .93232 .36932 L F .953 .953 1 r .93232 .36932 m .93232 .39674 L .91769 .39674 L .91769 .36932 L F .859 .859 1 r .91769 .36932 m .91769 .39674 L .90307 .39674 L .90307 .36932 L F .766 .766 1 r .90307 .36932 m .90307 .39674 L .88844 .39674 L .88844 .36932 L F 1 .818 .818 r .97619 .39674 m .97619 .42416 L .96157 .42416 L .96157 .39674 L F 1 .891 .891 r .96157 .39674 m .96157 .42416 L .94694 .42416 L .94694 .39674 L F 1 .964 .964 r .94694 .39674 m .94694 .42416 L .93232 .42416 L .93232 .39674 L F .964 .964 1 r .93232 .39674 m .93232 .42416 L .91769 .42416 L .91769 .39674 L F .891 .891 1 r .91769 .39674 m .91769 .42416 L .90307 .42416 L .90307 .39674 L F .818 .818 1 r .90307 .39674 m .90307 .42416 L .88844 .42416 L .88844 .39674 L F 1 .87 .87 r .97619 .42416 m .97619 .45158 L .96157 .45158 L .96157 .42416 L F 1 .922 .922 r .96157 .42416 m .96157 .45158 L .94694 .45158 L .94694 .42416 L F 1 .974 .974 r .94694 .42416 m .94694 .45158 L .93232 .45158 L .93232 .42416 L F .974 .974 1 r .93232 .42416 m .93232 .45158 L .91769 .45158 L .91769 .42416 L F .922 .922 1 r .91769 .42416 m .91769 .45158 L .90307 .45158 L .90307 .42416 L F .87 .87 1 r .90307 .42416 m .90307 .45158 L .88844 .45158 L .88844 .42416 L F 1 .922 .922 r .97619 .45158 m .97619 .479 L .96157 .479 L .96157 .45158 L F 1 .953 .953 r .96157 .45158 m .96157 .479 L .94694 .479 L .94694 .45158 L F 1 .984 .984 r .94694 .45158 m .94694 .479 L .93232 .479 L .93232 .45158 L F .984 .984 1 r .93232 .45158 m .93232 .479 L .91769 .479 L .91769 .45158 L F .953 .953 1 r .91769 .45158 m .91769 .479 L .90307 .479 L .90307 .45158 L F .922 .922 1 r .90307 .45158 m .90307 .479 L .88844 .479 L .88844 .45158 L F 1 .974 .974 r .97619 .479 m .97619 .50643 L .96157 .50643 L .96157 .479 L F 1 .984 .984 r .96157 .479 m .96157 .50643 L .94694 .50643 L .94694 .479 L F 1 .995 .995 r .94694 .479 m .94694 .50643 L .93232 .50643 L .93232 .479 L F .995 .995 1 r .93232 .479 m .93232 .50643 L .91769 .50643 L .91769 .479 L F .984 .984 1 r .91769 .479 m .91769 .50643 L .90307 .50643 L .90307 .479 L F .974 .974 1 r .90307 .479 m .90307 .50643 L .88844 .50643 L .88844 .479 L F .97619 .50643 m .97619 .53385 L .96157 .53385 L .96157 .50643 L F .984 .984 1 r .96157 .50643 m .96157 .53385 L .94694 .53385 L .94694 .50643 L F .995 .995 1 r .94694 .50643 m .94694 .53385 L .93232 .53385 L .93232 .50643 L F 1 .995 .995 r .93232 .50643 m .93232 .53385 L .91769 .53385 L .91769 .50643 L F 1 .984 .984 r .91769 .50643 m .91769 .53385 L .90307 .53385 L .90307 .50643 L F 1 .974 .974 r .90307 .50643 m .90307 .53385 L .88844 .53385 L .88844 .50643 L F .922 .922 1 r .97619 .53385 m .97619 .56127 L .96157 .56127 L .96157 .53385 L F .953 .953 1 r .96157 .53385 m .96157 .56127 L .94694 .56127 L .94694 .53385 L F .984 .984 1 r .94694 .53385 m .94694 .56127 L .93232 .56127 L .93232 .53385 L F 1 .984 .984 r .93232 .53385 m .93232 .56127 L .91769 .56127 L .91769 .53385 L F 1 .953 .953 r .91769 .53385 m .91769 .56127 L .90307 .56127 L .90307 .53385 L F 1 .922 .922 r .90307 .53385 m .90307 .56127 L .88844 .56127 L .88844 .53385 L F .87 .87 1 r .97619 .56127 m .97619 .58869 L .96157 .58869 L .96157 .56127 L F .922 .922 1 r .96157 .56127 m .96157 .58869 L .94694 .58869 L .94694 .56127 L F .974 .974 1 r .94694 .56127 m .94694 .58869 L .93232 .58869 L .93232 .56127 L F 1 .974 .974 r .93232 .56127 m .93232 .58869 L .91769 .58869 L .91769 .56127 L F 1 .922 .922 r .91769 .56127 m .91769 .58869 L .90307 .58869 L .90307 .56127 L F 1 .87 .87 r .90307 .56127 m .90307 .58869 L .88844 .58869 L .88844 .56127 L F .818 .818 1 r .97619 .58869 m .97619 .61611 L .96157 .61611 L .96157 .58869 L F .891 .891 1 r .96157 .58869 m .96157 .61611 L .94694 .61611 L .94694 .58869 L F .964 .964 1 r .94694 .58869 m .94694 .61611 L .93232 .61611 L .93232 .58869 L F 1 .964 .964 r .93232 .58869 m .93232 .61611 L .91769 .61611 L .91769 .58869 L F 1 .891 .891 r .91769 .58869 m .91769 .61611 L .90307 .61611 L .90307 .58869 L F 1 .818 .818 r .90307 .58869 m .90307 .61611 L .88844 .61611 L .88844 .58869 L F .766 .766 1 r .97619 .61611 m .97619 .64353 L .96157 .64353 L .96157 .61611 L F .859 .859 1 r .96157 .61611 m .96157 .64353 L .94694 .64353 L .94694 .61611 L F .953 .953 1 r .94694 .61611 m .94694 .64353 L .93232 .64353 L .93232 .61611 L F 1 .953 .953 r .93232 .61611 m .93232 .64353 L .91769 .64353 L .91769 .61611 L F 1 .859 .859 r .91769 .61611 m .91769 .64353 L .90307 .64353 L .90307 .61611 L F 1 .766 .766 r .90307 .61611 m .90307 .64353 L .88844 .64353 L .88844 .61611 L F .714 .714 1 r .97619 .64353 m .97619 .67095 L .96157 .67095 L .96157 .64353 L F .828 .828 1 r .96157 .64353 m .96157 .67095 L .94694 .67095 L .94694 .64353 L F .943 .943 1 r .94694 .64353 m .94694 .67095 L .93232 .67095 L .93232 .64353 L F 1 .943 .943 r .93232 .64353 m .93232 .67095 L .91769 .67095 L .91769 .64353 L F 1 .828 .828 r .91769 .64353 m .91769 .67095 L .90307 .67095 L .90307 .64353 L F 1 .714 .714 r .90307 .64353 m .90307 .67095 L .88844 .67095 L .88844 .64353 L F .661 .661 1 r .97619 .67095 m .97619 .69837 L .96157 .69837 L .96157 .67095 L F .797 .797 1 r .96157 .67095 m .96157 .69837 L .94694 .69837 L .94694 .67095 L F .932 .932 1 r .94694 .67095 m .94694 .69837 L .93232 .69837 L .93232 .67095 L F 1 .932 .932 r .93232 .67095 m .93232 .69837 L .91769 .69837 L .91769 .67095 L F 1 .797 .797 r .91769 .67095 m .91769 .69837 L .90307 .69837 L .90307 .67095 L F 1 .661 .661 r .90307 .67095 m .90307 .69837 L .88844 .69837 L .88844 .67095 L F .609 .609 1 r .97619 .69837 m .97619 .7258 L .96157 .7258 L .96157 .69837 L F .766 .766 1 r .96157 .69837 m .96157 .7258 L .94694 .7258 L .94694 .69837 L F .922 .922 1 r .94694 .69837 m .94694 .7258 L .93232 .7258 L .93232 .69837 L F 1 .922 .922 r .93232 .69837 m .93232 .7258 L .91769 .7258 L .91769 .69837 L F 1 .766 .766 r .91769 .69837 m .91769 .7258 L .90307 .7258 L .90307 .69837 L F 1 .609 .609 r .90307 .69837 m .90307 .7258 L .88844 .7258 L .88844 .69837 L F .557 .557 1 r .97619 .7258 m .97619 .75322 L .96157 .75322 L .96157 .7258 L F .734 .734 1 r .96157 .7258 m .96157 .75322 L .94694 .75322 L .94694 .7258 L F .911 .911 1 r .94694 .7258 m .94694 .75322 L .93232 .75322 L .93232 .7258 L F 1 .911 .911 r .93232 .7258 m .93232 .75322 L .91769 .75322 L .91769 .7258 L F 1 .734 .734 r .91769 .7258 m .91769 .75322 L .90307 .75322 L .90307 .7258 L F 1 .557 .557 r .90307 .7258 m .90307 .75322 L .88844 .75322 L .88844 .7258 L F .505 .505 1 r .97619 .75322 m .97619 .78064 L .96157 .78064 L .96157 .75322 L F .703 .703 1 r .96157 .75322 m .96157 .78064 L .94694 .78064 L .94694 .75322 L F .901 .901 1 r .94694 .75322 m .94694 .78064 L .93232 .78064 L .93232 .75322 L F 1 .901 .901 r .93232 .75322 m .93232 .78064 L .91769 .78064 L .91769 .75322 L F 1 .703 .703 r .91769 .75322 m .91769 .78064 L .90307 .78064 L .90307 .75322 L F 1 .505 .505 r .90307 .75322 m .90307 .78064 L .88844 .78064 L .88844 .75322 L F .453 .453 1 r .97619 .78064 m .97619 .80806 L .96157 .80806 L .96157 .78064 L F .672 .672 1 r .96157 .78064 m .96157 .80806 L .94694 .80806 L .94694 .78064 L F .891 .891 1 r .94694 .78064 m .94694 .80806 L .93232 .80806 L .93232 .78064 L F 1 .891 .891 r .93232 .78064 m .93232 .80806 L .91769 .80806 L .91769 .78064 L F 1 .672 .672 r .91769 .78064 m .91769 .80806 L .90307 .80806 L .90307 .78064 L F 1 .453 .453 r .90307 .78064 m .90307 .80806 L .88844 .80806 L .88844 .78064 L F .401 .401 1 r .97619 .80806 m .97619 .83548 L .96157 .83548 L .96157 .80806 L F .641 .641 1 r .96157 .80806 m .96157 .83548 L .94694 .83548 L .94694 .80806 L F .88 .88 1 r .94694 .80806 m .94694 .83548 L .93232 .83548 L .93232 .80806 L F 1 .88 .88 r .93232 .80806 m .93232 .83548 L .91769 .83548 L .91769 .80806 L F 1 .641 .641 r .91769 .80806 m .91769 .83548 L .90307 .83548 L .90307 .80806 L F 1 .401 .401 r .90307 .80806 m .90307 .83548 L .88844 .83548 L .88844 .80806 L F .349 .349 1 r .97619 .83548 m .97619 .8629 L .96157 .8629 L .96157 .83548 L F .609 .609 1 r .96157 .83548 m .96157 .8629 L .94694 .8629 L .94694 .83548 L F .87 .87 1 r .94694 .83548 m .94694 .8629 L .93232 .8629 L .93232 .83548 L F 1 .87 .87 r .93232 .83548 m .93232 .8629 L .91769 .8629 L .91769 .83548 L F 1 .609 .609 r .91769 .83548 m .91769 .8629 L .90307 .8629 L .90307 .83548 L F 1 .349 .349 r .90307 .83548 m .90307 .8629 L .88844 .8629 L .88844 .83548 L F .297 .297 1 r .97619 .8629 m .97619 .89032 L .96157 .89032 L .96157 .8629 L F .578 .578 1 r .96157 .8629 m .96157 .89032 L .94694 .89032 L .94694 .8629 L F .859 .859 1 r .94694 .8629 m .94694 .89032 L .93232 .89032 L .93232 .8629 L F 1 .859 .859 r .93232 .8629 m .93232 .89032 L .91769 .89032 L .91769 .8629 L F 1 .578 .578 r .91769 .8629 m .91769 .89032 L .90307 .89032 L .90307 .8629 L F 1 .297 .297 r .90307 .8629 m .90307 .89032 L .88844 .89032 L .88844 .8629 L F .245 .245 1 r .97619 .89032 m .97619 .91775 L .96157 .91775 L .96157 .89032 L F .547 .547 1 r .96157 .89032 m .96157 .91775 L .94694 .91775 L .94694 .89032 L F .849 .849 1 r .94694 .89032 m .94694 .91775 L .93232 .91775 L .93232 .89032 L F 1 .849 .849 r .93232 .89032 m .93232 .91775 L .91769 .91775 L .91769 .89032 L F 1 .547 .547 r .91769 .89032 m .91769 .91775 L .90307 .91775 L .90307 .89032 L F 1 .245 .245 r .90307 .89032 m .90307 .91775 L .88844 .91775 L .88844 .89032 L F .193 .193 1 r .97619 .91775 m .97619 .94517 L .96157 .94517 L .96157 .91775 L F .516 .516 1 r .96157 .91775 m .96157 .94517 L .94694 .94517 L .94694 .91775 L F .839 .839 1 r .94694 .91775 m .94694 .94517 L .93232 .94517 L .93232 .91775 L F 1 .839 .839 r .93232 .91775 m .93232 .94517 L .91769 .94517 L .91769 .91775 L F 1 .516 .516 r .91769 .91775 m .91769 .94517 L .90307 .94517 L .90307 .91775 L F 1 .193 .193 r .90307 .91775 m .90307 .94517 L .88844 .94517 L .88844 .91775 L F 1 .327 .327 r .08586 .03666 m .11328 .06408 L .10294 .07442 L .07552 .047 L F 1 .596 .596 r .07552 .047 m .10294 .07442 L .0926 .08476 L .06517 .05734 L F 1 .865 .865 r .06517 .05734 m .0926 .08476 L .08225 .09511 L .05483 .06768 L F .865 .865 1 r .05483 .06768 m .08225 .09511 L .07191 .10545 L .04449 .07802 L F .596 .596 1 r .04449 .07802 m .07191 .10545 L .06157 .11579 L .03415 .08837 L F .327 .327 1 r .03415 .08837 m .06157 .11579 L .05123 .12613 L .02381 .09871 L F 1 .371 .371 r .11328 .06408 m .1407 .0915 L .13036 .10184 L .10294 .07442 L F 1 .622 .622 r .10294 .07442 m .13036 .10184 L .12002 .11219 L .0926 .08476 L F 1 .874 .874 r .0926 .08476 m .12002 .11219 L .10968 .12253 L .08225 .09511 L F .874 .874 1 r .08225 .09511 m .10968 .12253 L .09933 .13287 L .07191 .10545 L F .622 .622 1 r .07191 .10545 m .09933 .13287 L .08899 .14321 L .06157 .11579 L F .371 .371 1 r .06157 .11579 m .08899 .14321 L .07865 .15355 L .05123 .12613 L F 1 .414 .414 r .1407 .0915 m .16812 .11892 L .15778 .12927 L .13036 .10184 L F 1 .648 .648 r .13036 .10184 m .15778 .12927 L .14744 .13961 L .12002 .11219 L F 1 .883 .883 r .12002 .11219 m .14744 .13961 L .1371 .14995 L .10968 .12253 L F .883 .883 1 r .10968 .12253 m .1371 .14995 L .12676 .16029 L .09933 .13287 L F .648 .648 1 r .09933 .13287 m .12676 .16029 L .11641 .17063 L .08899 .14321 L F .414 .414 1 r .08899 .14321 m .11641 .17063 L .10607 .18097 L .07865 .15355 L F 1 .457 .457 r .16812 .11892 m .19554 .14635 L .1852 .15669 L .15778 .12927 L F 1 .674 .674 r .15778 .12927 m .1852 .15669 L .17486 .16703 L .14744 .13961 L F 1 .891 .891 r .14744 .13961 m .17486 .16703 L .16452 .17737 L .1371 .14995 L F .891 .891 1 r .1371 .14995 m .16452 .17737 L .15418 .18771 L .12676 .16029 L F .674 .674 1 r .12676 .16029 m .15418 .18771 L .14384 .19805 L .11641 .17063 L F .457 .457 1 r .11641 .17063 m .14384 .19805 L .13349 .20839 L .10607 .18097 L F 1 .501 .501 r .19554 .14635 m .22296 .17377 L .21262 .18411 L .1852 .15669 L F 1 .701 .701 r .1852 .15669 m .21262 .18411 L .20228 .19445 L .17486 .16703 L F 1 .9 .9 r .17486 .16703 m .20228 .19445 L .19194 .20479 L .16452 .17737 L F .9 .9 1 r .16452 .17737 m .19194 .20479 L .1816 .21513 L .15418 .18771 L F .701 .701 1 r .15418 .18771 m .1816 .21513 L .17126 .22547 L .14384 .19805 L F .501 .501 1 r .14384 .19805 m .17126 .22547 L .16092 .23581 L .13349 .20839 L F 1 .544 .544 r .22296 .17377 m .25039 .20119 L .24004 .21153 L .21262 .18411 L F 1 .727 .727 r .21262 .18411 m .24004 .21153 L .2297 .22187 L .20228 .19445 L F 1 .909 .909 r .20228 .19445 m .2297 .22187 L .21936 .23221 L .19194 .20479 L F .909 .909 1 r .19194 .20479 m .21936 .23221 L .20902 .24255 L .1816 .21513 L F .727 .727 1 r .1816 .21513 m .20902 .24255 L .19868 .25289 L .17126 .22547 L F .544 .544 1 r .17126 .22547 m .19868 .25289 L .18834 .26324 L .16092 .23581 L F 1 .588 .588 r .25039 .20119 m .27781 .22861 L .26747 .23895 L .24004 .21153 L F 1 .753 .753 r .24004 .21153 m .26747 .23895 L .25712 .24929 L .2297 .22187 L F 1 .918 .918 r .2297 .22187 m .25712 .24929 L .24678 .25963 L .21936 .23221 L F .918 .918 1 r .21936 .23221 m .24678 .25963 L .23644 .26997 L .20902 .24255 L F .753 .753 1 r .20902 .24255 m .23644 .26997 L .2261 .28032 L .19868 .25289 L F .588 .588 1 r .19868 .25289 m .2261 .28032 L .21576 .29066 L .18834 .26324 L F 1 .631 .631 r .27781 .22861 m .30523 .25603 L .29489 .26637 L .26747 .23895 L F 1 .779 .779 r .26747 .23895 m .29489 .26637 L .28455 .27671 L .25712 .24929 L F 1 .926 .926 r .25712 .24929 m .28455 .27671 L .2742 .28705 L .24678 .25963 L F .926 .926 1 r .24678 .25963 m .2742 .28705 L .26386 .2974 L .23644 .26997 L F .779 .779 1 r .23644 .26997 m .26386 .2974 L .25352 .30774 L .2261 .28032 L F .631 .631 1 r .2261 .28032 m .25352 .30774 L .24318 .31808 L .21576 .29066 L F 1 .674 .674 r .30523 .25603 m .33265 .28345 L .32231 .29379 L .29489 .26637 L F 1 .805 .805 r .29489 .26637 m .32231 .29379 L .31197 .30413 L .28455 .27671 L F 1 .935 .935 r .28455 .27671 m .31197 .30413 L .30163 .31448 L .2742 .28705 L F .935 .935 1 r .2742 .28705 m .30163 .31448 L .29128 .32482 L .26386 .2974 L F .805 .805 1 r .26386 .2974 m .29128 .32482 L .28094 .33516 L .25352 .30774 L F .674 .674 1 r .25352 .30774 m .28094 .33516 L .2706 .3455 L .24318 .31808 L F 1 .718 .718 r .33265 .28345 m .36007 .31087 L .34973 .32121 L .32231 .29379 L F 1 .831 .831 r .32231 .29379 m .34973 .32121 L .33939 .33156 L .31197 .30413 L F 1 .944 .944 r .31197 .30413 m .33939 .33156 L .32905 .3419 L .30163 .31448 L F .944 .944 1 r .30163 .31448 m .32905 .3419 L .31871 .35224 L .29128 .32482 L F .831 .831 1 r .29128 .32482 m .31871 .35224 L .30836 .36258 L .28094 .33516 L F .718 .718 1 r .28094 .33516 m .30836 .36258 L .29802 .37292 L .2706 .3455 L F 1 .761 .761 r .36007 .31087 m .38749 .33829 L .37715 .34864 L .34973 .32121 L F 1 .857 .857 r .34973 .32121 m .37715 .34864 L .36681 .35898 L .33939 .33156 L F 1 .952 .952 r .33939 .33156 m .36681 .35898 L .35647 .36932 L .32905 .3419 L F .952 .952 1 r .32905 .3419 m .35647 .36932 L .34613 .37966 L .31871 .35224 L F .857 .857 1 r .31871 .35224 m .34613 .37966 L .33579 .39 L .30836 .36258 L F .761 .761 1 r .30836 .36258 m .33579 .39 L .32544 .40034 L .29802 .37292 L F 1 .805 .805 r .38749 .33829 m .41491 .36572 L .40457 .37606 L .37715 .34864 L F 1 .883 .883 r .37715 .34864 m .40457 .37606 L .39423 .3864 L .36681 .35898 L F 1 .961 .961 r .36681 .35898 m .39423 .3864 L .38389 .39674 L .35647 .36932 L F .961 .961 1 r .35647 .36932 m .38389 .39674 L .37355 .40708 L .34613 .37966 L F .883 .883 1 r .34613 .37966 m .37355 .40708 L .36321 .41742 L .33579 .39 L F .805 .805 1 r .33579 .39 m .36321 .41742 L .35287 .42776 L .32544 .40034 L F 1 .848 .848 r .41491 .36572 m .44233 .39314 L .43199 .40348 L .40457 .37606 L F 1 .909 .909 r .40457 .37606 m .43199 .40348 L .42165 .41382 L .39423 .3864 L F 1 .97 .97 r .39423 .3864 m .42165 .41382 L .41131 .42416 L .38389 .39674 L F .97 .97 1 r .38389 .39674 m .41131 .42416 L .40097 .4345 L .37355 .40708 L F .909 .909 1 r .37355 .40708 m .40097 .4345 L .39063 .44484 L .36321 .41742 L F .848 .848 1 r .36321 .41742 m .39063 .44484 L .38029 .45518 L .35287 .42776 L F 1 .891 .891 r .44233 .39314 m .46976 .42056 L .45941 .4309 L .43199 .40348 L F 1 .935 .935 r .43199 .40348 m .45941 .4309 L .44907 .44124 L .42165 .41382 L F 1 .978 .978 r .42165 .41382 m .44907 .44124 L .43873 .45158 L .41131 .42416 L F .978 .978 1 r .41131 .42416 m .43873 .45158 L .42839 .46192 L .40097 .4345 L F .935 .935 1 r .40097 .4345 m .42839 .46192 L .41805 .47227 L .39063 .44484 L F .891 .891 1 r .39063 .44484 m .41805 .47227 L .40771 .48261 L .38029 .45518 L F 1 .935 .935 r .46976 .42056 m .49718 .44798 L .48684 .45832 L .45941 .4309 L F 1 .961 .961 r .45941 .4309 m .48684 .45832 L .47649 .46866 L .44907 .44124 L F 1 .987 .987 r .44907 .44124 m .47649 .46866 L .46615 .479 L .43873 .45158 L F .987 .987 1 r .43873 .45158 m .46615 .479 L .45581 .48935 L .42839 .46192 L F .961 .961 1 r .42839 .46192 m .45581 .48935 L .44547 .49969 L .41805 .47227 L F .935 .935 1 r .41805 .47227 m .44547 .49969 L .43513 .51003 L .40771 .48261 L F 1 .978 .978 r .49718 .44798 m .5246 .4754 L .51426 .48574 L .48684 .45832 L F 1 .987 .987 r .48684 .45832 m .51426 .48574 L .50392 .49608 L .47649 .46866 L F 1 .996 .996 r .47649 .46866 m .50392 .49608 L .49357 .50643 L .46615 .479 L F .996 .996 1 r .46615 .479 m .49357 .50643 L .48323 .51677 L .45581 .48935 L F .987 .987 1 r .45581 .48935 m .48323 .51677 L .47289 .52711 L .44547 .49969 L F .978 .978 1 r .44547 .49969 m .47289 .52711 L .46255 .53745 L .43513 .51003 L F .5246 .4754 m .55202 .50282 L .54168 .51316 L .51426 .48574 L F .987 .987 1 r .51426 .48574 m .54168 .51316 L .53134 .52351 L .50392 .49608 L F .996 .996 1 r .50392 .49608 m .53134 .52351 L .521 .53385 L .49357 .50643 L F 1 .996 .996 r .49357 .50643 m .521 .53385 L .51065 .54419 L .48323 .51677 L F 1 .987 .987 r .48323 .51677 m .51065 .54419 L .50031 .55453 L .47289 .52711 L F 1 .978 .978 r .47289 .52711 m .50031 .55453 L .48997 .56487 L .46255 .53745 L F .935 .935 1 r .55202 .50282 m .57944 .53024 L .5691 .54059 L .54168 .51316 L F .961 .961 1 r .54168 .51316 m .5691 .54059 L .55876 .55093 L .53134 .52351 L F .987 .987 1 r .53134 .52351 m .55876 .55093 L .54842 .56127 L .521 .53385 L F 1 .987 .987 r .521 .53385 m .54842 .56127 L .53808 .57161 L .51065 .54419 L F 1 .961 .961 r .51065 .54419 m .53808 .57161 L .52773 .58195 L .50031 .55453 L F 1 .935 .935 r .50031 .55453 m .52773 .58195 L .51739 .59229 L .48997 .56487 L F .891 .891 1 r .57944 .53024 m .60686 .55767 L .59652 .56801 L .5691 .54059 L F .935 .935 1 r .5691 .54059 m .59652 .56801 L .58618 .57835 L .55876 .55093 L F .978 .978 1 r .55876 .55093 m .58618 .57835 L .57584 .58869 L .54842 .56127 L F 1 .978 .978 r .54842 .56127 m .57584 .58869 L .5655 .59903 L .53808 .57161 L F 1 .935 .935 r .53808 .57161 m .5655 .59903 L .55516 .60937 L .52773 .58195 L F 1 .891 .891 r .52773 .58195 m .55516 .60937 L .54482 .61971 L .51739 .59229 L F .848 .848 1 r .60686 .55767 m .63428 .58509 L .62394 .59543 L .59652 .56801 L F .909 .909 1 r .59652 .56801 m .62394 .59543 L .6136 .60577 L .58618 .57835 L F .97 .97 1 r .58618 .57835 m .6136 .60577 L .60326 .61611 L .57584 .58869 L F 1 .97 .97 r .57584 .58869 m .60326 .61611 L .59292 .62645 L .5655 .59903 L F 1 .909 .909 r .5655 .59903 m .59292 .62645 L .58258 .63679 L .55516 .60937 L F 1 .848 .848 r .55516 .60937 m .58258 .63679 L .57224 .64713 L .54482 .61971 L F .805 .805 1 r .63428 .58509 m .66171 .61251 L .65136 .62285 L .62394 .59543 L F .883 .883 1 r .62394 .59543 m .65136 .62285 L .64102 .63319 L .6136 .60577 L F .961 .961 1 r .6136 .60577 m .64102 .63319 L .63068 .64353 L .60326 .61611 L F 1 .961 .961 r .60326 .61611 m .63068 .64353 L .62034 .65387 L .59292 .62645 L F 1 .883 .883 r .59292 .62645 m .62034 .65387 L .61 .66421 L .58258 .63679 L F 1 .805 .805 r .58258 .63679 m .61 .66421 L .59966 .67456 L .57224 .64713 L F .761 .761 1 r .66171 .61251 m .68913 .63993 L .67879 .65027 L .65136 .62285 L F .857 .857 1 r .65136 .62285 m .67879 .65027 L .66844 .66061 L .64102 .63319 L F .952 .952 1 r .64102 .63319 m .66844 .66061 L .6581 .67095 L .63068 .64353 L F 1 .952 .952 r .63068 .64353 m .6581 .67095 L .64776 .68129 L .62034 .65387 L F 1 .857 .857 r .62034 .65387 m .64776 .68129 L .63742 .69164 L .61 .66421 L F 1 .761 .761 r .61 .66421 m .63742 .69164 L .62708 .70198 L .59966 .67456 L F .718 .718 1 r .68913 .63993 m .71655 .66735 L .70621 .67769 L .67879 .65027 L F .831 .831 1 r .67879 .65027 m .70621 .67769 L .69587 .68803 L .66844 .66061 L F .944 .944 1 r .66844 .66061 m .69587 .68803 L .68552 .69837 L .6581 .67095 L F 1 .944 .944 r .6581 .67095 m .68552 .69837 L .67518 .70872 L .64776 .68129 L F 1 .831 .831 r .64776 .68129 m .67518 .70872 L .66484 .71906 L .63742 .69164 L F 1 .718 .718 r .63742 .69164 m .66484 .71906 L .6545 .7294 L .62708 .70198 L F .674 .674 1 r .71655 .66735 m .74397 .69477 L .73363 .70511 L .70621 .67769 L F .805 .805 1 r .70621 .67769 m .73363 .70511 L .72329 .71545 L .69587 .68803 L F .935 .935 1 r .69587 .68803 m .72329 .71545 L .71295 .7258 L .68552 .69837 L F 1 .935 .935 r .68552 .69837 m .71295 .7258 L .7026 .73614 L .67518 .70872 L F 1 .805 .805 r .67518 .70872 m .7026 .73614 L .69226 .74648 L .66484 .71906 L F 1 .674 .674 r .66484 .71906 m .69226 .74648 L .68192 .75682 L .6545 .7294 L F .631 .631 1 r .74397 .69477 m .77139 .72219 L .76105 .73253 L .73363 .70511 L F .779 .779 1 r .73363 .70511 m .76105 .73253 L .75071 .74288 L .72329 .71545 L F .926 .926 1 r .72329 .71545 m .75071 .74288 L .74037 .75322 L .71295 .7258 L F 1 .926 .926 r .71295 .7258 m .74037 .75322 L .73003 .76356 L .7026 .73614 L F 1 .779 .779 r .7026 .73614 m .73003 .76356 L .71968 .7739 L .69226 .74648 L F 1 .631 .631 r .69226 .74648 m .71968 .7739 L .70934 .78424 L .68192 .75682 L F .588 .588 1 r .77139 .72219 m .79881 .74961 L .78847 .75996 L .76105 .73253 L F .753 .753 1 r .76105 .73253 m .78847 .75996 L .77813 .7703 L .75071 .74288 L F .918 .918 1 r .75071 .74288 m .77813 .7703 L .76779 .78064 L .74037 .75322 L F 1 .918 .918 r .74037 .75322 m .76779 .78064 L .75745 .79098 L .73003 .76356 L F 1 .753 .753 r .73003 .76356 m .75745 .79098 L .74711 .80132 L .71968 .7739 L F 1 .588 .588 r .71968 .7739 m .74711 .80132 L .73676 .81166 L .70934 .78424 L F .544 .544 1 r .79881 .74961 m .82623 .77704 L .81589 .78738 L .78847 .75996 L F .727 .727 1 r .78847 .75996 m .81589 .78738 L .80555 .79772 L .77813 .7703 L F .909 .909 1 r .77813 .7703 m .80555 .79772 L .79521 .80806 L .76779 .78064 L F 1 .909 .909 r .76779 .78064 m .79521 .80806 L .78487 .8184 L .75745 .79098 L F 1 .727 .727 r .75745 .79098 m .78487 .8184 L .77453 .82874 L .74711 .80132 L F 1 .544 .544 r .74711 .80132 m .77453 .82874 L .76419 .83908 L .73676 .81166 L F .501 .501 1 r .82623 .77704 m .85365 .80446 L .84331 .8148 L .81589 .78738 L F .701 .701 1 r .81589 .78738 m .84331 .8148 L .83297 .82514 L .80555 .79772 L F .9 .9 1 r .80555 .79772 m .83297 .82514 L .82263 .83548 L .79521 .80806 L F 1 .9 .9 r .79521 .80806 m .82263 .83548 L .81229 .84582 L .78487 .8184 L F 1 .701 .701 r .78487 .8184 m .81229 .84582 L .80195 .85616 L .77453 .82874 L F 1 .501 .501 r .77453 .82874 m .80195 .85616 L .79161 .86651 L .76419 .83908 L F .457 .457 1 r .85365 .80446 m .88108 .83188 L .87073 .84222 L .84331 .8148 L F .674 .674 1 r .84331 .8148 m .87073 .84222 L .86039 .85256 L .83297 .82514 L F .891 .891 1 r .83297 .82514 m .86039 .85256 L .85005 .8629 L .82263 .83548 L F 1 .891 .891 r .82263 .83548 m .85005 .8629 L .83971 .87324 L .81229 .84582 L F 1 .674 .674 r .81229 .84582 m .83971 .87324 L .82937 .88359 L .80195 .85616 L F 1 .457 .457 r .80195 .85616 m .82937 .88359 L .81903 .89393 L .79161 .86651 L F .414 .414 1 r .88108 .83188 m .9085 .8593 L .89816 .86964 L .87073 .84222 L F .648 .648 1 r .87073 .84222 m .89816 .86964 L .88781 .87998 L .86039 .85256 L F .883 .883 1 r .86039 .85256 m .88781 .87998 L .87747 .89032 L .85005 .8629 L F 1 .883 .883 r .85005 .8629 m .87747 .89032 L .86713 .90067 L .83971 .87324 L F 1 .648 .648 r .83971 .87324 m .86713 .90067 L .85679 .91101 L .82937 .88359 L F 1 .414 .414 r .82937 .88359 m .85679 .91101 L .84645 .92135 L .81903 .89393 L F .371 .371 1 r .9085 .8593 m .93592 .88672 L .92558 .89706 L .89816 .86964 L F .622 .622 1 r .89816 .86964 m .92558 .89706 L .91524 .9074 L .88781 .87998 L F .874 .874 1 r .88781 .87998 m .91524 .9074 L .90489 .91775 L .87747 .89032 L F 1 .874 .874 r .87747 .89032 m .90489 .91775 L .89455 .92809 L .86713 .90067 L F 1 .622 .622 r .86713 .90067 m .89455 .92809 L .88421 .93843 L .85679 .91101 L F 1 .371 .371 r .85679 .91101 m .88421 .93843 L .87387 .94877 L .84645 .92135 L F .327 .327 1 r .93592 .88672 m .96334 .91414 L .953 .92448 L .92558 .89706 L F .596 .596 1 r .92558 .89706 m .953 .92448 L .94266 .93483 L .91524 .9074 L F .865 .865 1 r .91524 .9074 m .94266 .93483 L .93232 .94517 L .90489 .91775 L F 1 .865 .865 r .90489 .91775 m .93232 .94517 L .92198 .95551 L .89455 .92809 L F 1 .596 .596 r .89455 .92809 m .92198 .95551 L .91163 .96585 L .88421 .93843 L F 1 .327 .327 r .88421 .93843 m .91163 .96585 L .90129 .97619 L .87387 .94877 L F % End of Graphics MathPictureEnd \ \>"], "Graphics", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageSize->{246, 246}, ImageMargins->{{7, 0}, {0, 0}}, ImageRegion->{{-0, 1}, {0, 1}}, ImageCacheValid->False] }, Open ]], Cell["\<\ Cell 5B. Plot Axial Stress level displays stress level and sign \ using color: red for tension, blue for compression, white for zero stress. A black background is \ used.\ \>", "Text", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1]], Cell[CellGroupData[{ Cell[TextData[{ "PlotAxialStressLevel[nodcoor_,elenod_,s_,sfac_,aspect_,title_]:= \n \ Module[{eNL,e,n,ni,nj,nc,xyc,sides={},numele=Length[elenod],\n \ numnod=Length[nodcoor],x,y,xmin,xmax,ymin,ymax,\n \ smin,smax,fmax,sval,c1,c2,c3,aratio},\n \ x=Table[nodcoor[[i,1]],{i,1,numnod}]; \n \ y=Table[nodcoor[[i,2]],{i,1,numnod}];\n \ {xmin,xmax,ymin,ymax}=N[{Min[x],Max[x],Min[y],Max[y]}];\n smin=Min[s]; \ smax=Max[s]; fmax=Max[Abs[smax],Abs[smin]];\n For [e=1,e<=numele,e++, \ eNL=elenod[[e]]; nc=Length[eNL];\n sval=N[s[[e]]*sfac];\n \ {c1,c2,c3}=LineColor[sval,fmax];", StyleBox[" ", FontColor->RGBColor[1, 0, 0]], "\n If [nc!=2, Continue[]]; {ni,nj}=eNL; \n \ xyc={{x[[ni]],y[[ni]]},{x[[nj]],y[[nj]]}}; \n \ sides=AppendTo[sides,Graphics[RGBColor[c1,c2,c3]]];\n \ sides=AppendTo[sides,Graphics[Line[xyc]]]\n ];\n If [aspect>0, \ aratio=aspect, aratio=(ymax-ymin)/(xmax-xmin)];\n \ Show[Graphics[AbsoluteThickness[3]],\n Graphics[RGBColor[0,0,0]], sides, \ \n Background->GrayLevel[0],\n AspectRatio->aratio, \ PlotLabel->title ];\n ];\n\nLineColor[f_,fmax_]:= Module[{r,RGBmax={1,0,0}, \ \n RGBmin={0,0,1}, RGBzero={1,1,1}, RGBout={0,0,0}},\n If [f==0 || \ fmax==0, \n Return[RGBzero]]; (* White if f=0 *)\n If [f>fmax || \ f<-fmax, \n Return[RGBout ]]; (* Black if outside range *)\n If \ [f>0, r= N[f/fmax]; \n Return[r*RGBmax+(1-r)*RGBzero]]; (* positive *)\n\ If [f<0, r=-N[f/fmax]; \n Return[r*RGBmin+(1-r)*RGBzero]]; (* \ negative *)\n];\n", StyleBox["\nnodcoor={{0,0},{10,5},{10,0},{20,8},{20,0},{30,9},{30,0},\n \ {40,8},{40,0},{50,5},{50,0},{60,0}};\nelenod= \ {{1,3},{3,5},{5,7},{7,9},{9,11},{11,12},\n \ {1,2},{2,4},{4,6},{6,8},{8,10},{10,12},\n \ {2,3},{4,5},{6,7},{8,9},{10,11},\n {2,5},{4,7},{7,8},{9,10}};\n\ numele=Length[elenod];\nu={0,0, 0,-0.10, 0,-0.10, 0,-0.15, 0,-0.15, 0,-0.20, \ \n 0,-0.20, 0,-0.15, 0,-0.15, 0,-0.10, 0,-0.10, 0,0};\n\ sigma=Table[ne-10,{ne,1,numele}];\naspect=0;\n\ PlotAxialStressLevel[nodcoor,elenod,sigma,1.0,aspect,\n \"test axial stress \ level plot\"];\n \n(* \nu={0,0, 0,-0.10, 0,-0.10, 0,-0.15, 0,-0.15, \ 0,-0.20, \n 0,-0.20, 0,-0.15, 0,-0.15, 0,-0.10, 0,-0.10, 0,0};\nFor \ [t=0.,t<=N[Pi],t=t+N[Pi/6], amp=10*Sin[t];\n \ PlotDeformedShape[nodcoor,elenod,u,amp,aspect,\n \"deformed shape\"]]; *)", FontColor->RGBColor[1, 0, 0]] }], "Input", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0.699214, 1, 0.0500191]], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .15 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Background color 0 g MFill % Scaling calculations 0.0238095 0.015873 0.00357143 0.015873 [ [(test axial stress level plot)] .5 .15 0 -2 Msboxa [ 0 0 0 0 ] [ 1 .15 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 1 g p [(test axial stress level plot)] .5 .15 0 -2 Mshowa P 0 0 m 1 0 L 1 .15 L 0 .15 L closepath clip newpath p 0 0 0 r .182 .182 1 r 3 Mabswid .02381 .00357 m .18254 .00357 L s .273 .273 1 r .18254 .00357 m .34127 .00357 L s .364 .364 1 r .34127 .00357 m .5 .00357 L s .455 .455 1 r .5 .00357 m .65873 .00357 L s .545 .545 1 r .65873 .00357 m .81746 .00357 L s .636 .636 1 r .81746 .00357 m .97619 .00357 L s .727 .727 1 r .02381 .00357 m .18254 .08294 L s .818 .818 1 r .18254 .08294 m .34127 .13056 L s .909 .909 1 r .34127 .13056 m .5 .14643 L s 1 1 1 r .5 .14643 m .65873 .13056 L s 1 .909 .909 r .65873 .13056 m .81746 .08294 L s 1 .818 .818 r .81746 .08294 m .97619 .00357 L s 1 .727 .727 r .18254 .08294 m .18254 .00357 L s 1 .636 .636 r .34127 .13056 m .34127 .00357 L s 1 .545 .545 r .5 .14643 m .5 .00357 L s 1 .455 .455 r .65873 .13056 m .65873 .00357 L s 1 .364 .364 r .81746 .08294 m .81746 .00357 L s 1 .273 .273 r .18254 .08294 m .34127 .00357 L s 1 .182 .182 r .34127 .13056 m .5 .00357 L s 1 .091 .091 r .5 .00357 m .65873 .13056 L s 1 0 0 r .65873 .00357 m .81746 .08294 L s P % End of Graphics MathPictureEnd \ \>"], "Graphics", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageSize->{361, 54}, ImageMargins->{{34, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgMl6O00QASl3OF/PDYl6OeXo0001Ool003l000QNoahD_`=nmalQ7`Qm224]O`=l Y@Il@R0ig`Ql01m6?`=mJaiBW`QoFSl0005oo`00@@0025ko70000gkg7@0027d80WjD80001GbU00=n 27`07@000ge[700027mJ@@000Goo00130008G_lJ0003O_LK0008O@P01000 OYAnU7jD7P001WbU0gha1Wa27@000gg>00<007`0O0001W`06`000ge[6P0027mJ@`000Goo00150007 G_lI0003O_LI0008O@P30003OY@M0006O:D01@00OS5n0`0027`06@000ge[ 6@001gmJA@000Goo00170007G_lG0003O_LG0008O@P50003OY@L0006O:D20003OS420006O48K0003 OLh50008O00G0003OF/G0007OeY70001Ool004T000MNoaD000=nmaD000Qm20L000=nU1/000ElY@@0 00=n<@@000El@QX000=mcPL000Ql01D000=mJaD000MoFTT0005oo`00BP0025ko4`000gkg4`0027d8 2@000gjD6@001WbU1@000gha1@001Ga26@000gg>2@0027`04`000ge[4`0027mJBP000Goo001<0008 G_lA0003O_LA0008O@P;0003OY@H0006O:D60003OS460006O48G0003OLh;0008O00A0003OF/A0008 OeY<0001Ool004h000QNo`l000=nma0000Mm20d000=nU1L000IlY@L000=n<@L000Il@QH000=mcPd0 00Ml010000=mJ`l000QoFTh0005oo`00D00025ko3@000gkg3P001gd83`000gjD5P001GbU2@000gha 20001Wa25@000gg>3`001g`03P000ge[3@0027mJD0000Goo001B0008G_l;0003O_L<0007O@PA0003 OY@D0006O:D:0003OS4:0005O48D0003OLhA0007O00<0003OF/;0008OeYB0001Ool005@000QNo`T0 00=nm`X000Qm218000=nU1<000IlY@/000=n<@/000Il@Q8000=mcQ8000Ql00X000=mJ`T000QoFU@0 005oo`00EP0025ko1`000gkg200027d850000gjD4P001WbU30000gha30001Wa24@000gg>500027`0 20000ge[1`0027mJEP000Goo001H0007G_l60003O_L60008O@PF0003OY@@0006O:D>0003OS4=0006 O48@0003OLhF0008O0060003OF/60007OeYH0001Ool005X000MNo`@000=nm`@000Qm21P000=nU0l0 00IlY@l000=n<@l000El@Pl000=mcQP000Ql00@000=mJ`@000MoFUX0005oo`00G0001eko0P000gkg 0P0027d86P000gjD3P001WbU40000gha40001Ga23P000gg>6P0027`00P000ge[0P001gmJG0000Goo 001M0008G_l3O_L8O@PL0003OY@=0005O:DB0003OS4A0006O48<0003OLhL0008O003OF/8OeYM0001 Ool005l000INo`03O_Mm27d800Im21h000=nU0/000IlYA<000=n8000 1g`017mJH@000Goo001S0002G_l5O@P3JelO0003OY@90006O:DE0003OS4E0005O4890003OLhO0003 Okd5O002OeYS0001Ool006D000=m20U[Ga/000=nU0P000ElYAL000=n600037nmJ0000Goo 001[000=JelD0003OY@50006O:DI0003OS4H0006O4850003OLhD000=Oke[0001Ool006l000e[Ga00 00=nU0@000IlYAX000=n3@003GnmLP000Goo001f000=Jel90003OY@00`00O:ElY@04O:DM0003 OS4L0006O4801000OLimcWg>2@003GnmMP000Goo001j000"], ImageRangeCache->{{{0, 360}, {53, 0}} -> {-16.6032, -0.225048, 0.258907, \ 0.258907}}] }, Open ]], Cell["\<\ Cell 6. The main driver program. Defines the frame problem and \ runs the analysis.\ \>", "Text", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1]], Cell["\<\ ClearAll[L,H,a,Em,P]; L=10000; H=6000; Em=35000; a=500; A=a^2; Izz=a^4/12; P=4800; NodeCoordinates={{-L,0},{-L,H/2},{-L,H},{-L/2,H},{0,H},{L/2,H}, {L,H},{L,H/2},{L,0},{0,H/2},{0,0}}; ElemNodeLists= {{1,2},{2,3},{3,4},{4,5},{5,6},{6,7},{7,8}, {8,9},{5,10},{10,11}}; numnod=Length[NodeCoordinates]; numele=Length[ElemNodeLists]; numdof=3*numnod; ElemMaterial= Table[{Em,0,0,0},{numele}]; ElemFabrication= Table[{A,Izz},{numele}]; ProcessOptions= {True}; aspect=0; PlotLineElements[NodeCoordinates,ElemNodeLists,aspect, \"frame mesh\"]; PlotLineElementsAndNodes[NodeCoordinates,ElemNodeLists,aspect, \"frame mesh with elem & node labels\",{True,0.10},{True,0.04}]; FreedomTag=FreedomValue=Table[{0,0,0},{numnod}]; FreedomValue[[4]]={0,-P,0}; FreedomValue[[3]]={P/2,0,0}; Print[\"Applied node forces=\"]; Print[FreedomValue]; FreedomTag[[1]]=FreedomTag[[9]]=FreedomTag[[11]]={1,1,1}; (* fixed *) f=Flatten[FreedomValue]; K=PlaneFrameMasterStiffness[NodeCoordinates, ElemNodeLists,ElemMaterial,ElemFabrication,ProcessOptions]; pdof={}; For[n=1,n<=numnod,n++, For[j=1,j<=3,j++, If [FreedomTag[[n,j]]>0, AppendTo[pdof,3*(n-1)+j]]]]; Print[\"pdof=\",pdof]; Kmod=ModifiedMasterStiffness[pdof,K]; fmod=ModifyNodeForces [pdof,f]; u=LinearSolve[Kmod,fmod]; u=Chop[u,.000001]; Print[\"Computed Nodal Displacements:\"]; NodeForces=NodeDisplacements=Table[{0,0,0},{numnod}]; For [n=1,n<=numnod,n++, For[j=1,j<=3,j++, NodeDisplacements[[n,j]]=u[[3*(n-1)+j]] ]]; Print[NodeDisplacements//MatrixForm] f=Simplify[K.u]; f=Chop[f,.000001]; Print[\"External Node Forces Including Reactions:\"]; For [n=1,n<=numnod,n++, For[j=1,j<=3,j++, NodeForces[[n,j]]=f[[3*(n-1)+j]] ]]; Print[NodeForces//MatrixForm]; p=PlaneFrameIntForces[NodeCoordinates,ElemNodeLists, ElemMaterial,ElemFabrication,eleopt,u]; p=Chop[p,.000001]; Print[\"Internal Member Forces:\"]; Print[p//MatrixForm]; amp=10000.; PlotPlaneFrameDeformedShape[NodeCoordinates,ElemNodeLists,u, {amp,amp,0,16},aspect,\"Deformed shape (magnified x10000)\"]; H=500; ElemFabrication= Table[{A,Izz,H},{numele}]; fmax=0.25; Nsub{16,8}; aspect=-1; title=\"Axial stress in frame members\"; PlotPlaneFrameStress[NodeCoordinates,ElemNodeLists, ElemFabrication,p,fmax,Nsub,aspect,title]; \ \>", "Input", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 0]], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .3 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Background color 1 1 0 r MFill % Scaling calculations 0.5 4.7619e-05 0.00714286 4.7619e-05 [ [.5 .3125 -32 0 ] [.5 .3125 32 12 ] [ 0 0 0 0 ] [ 1 .3 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 0 1 r gsave .5 .3125 -93 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 0.000 0.000 rmoveto 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor (frame) show 99.000 13.000 moveto (mesh) show 123.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore 0 0 m 1 0 L 1 .3 L 0 .3 L closepath clip newpath 0 0 0 r 2 Mabswid [ ] 0 setdash .02381 .00714 m .02381 .15 L s .02381 .15 m .02381 .29286 L s .02381 .29286 m .2619 .29286 L s .2619 .29286 m .5 .29286 L s .5 .29286 m .7381 .29286 L s .7381 .29286 m .97619 .29286 L s .97619 .29286 m .97619 .15 L s .97619 .15 m .97619 .00714 L s .5 .29286 m .5 .15 L s .5 .15 m .5 .00714 L s % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 86.375}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHg"], ImageRangeCache->{{{0, 287}, {85.375, 0}} -> {-12834.5, -150.033, 89.4388, \ 89.4388}}], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .3 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Background color 1 1 0 r MFill % Scaling calculations 0.5 4.62963e-05 0.0328125 3.90625e-05 [ [.5 .3125 -104.125 0 ] [.5 .3125 104.125 12 ] [ 0 0 0 0 ] [ 1 .3 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 0 1 r gsave .5 .3125 -165.125 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 0.000 0.000 rmoveto 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor (frame) show 99.000 13.000 moveto (mesh) show 129.125 13.000 moveto (with) show 159.125 13.000 moveto (elem) show 189.250 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor (&) show 201.250 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor (node) show 231.250 13.000 moveto (labels) show 267.250 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore 0 0 m 1 0 L 1 .3 L 0 .3 L closepath clip newpath 0 g 3 Mabswid [ ] 0 setdash .03704 .03281 m .03704 .15 L s .03704 .15 m .03704 .26719 L s .03704 .26719 m .26852 .26719 L s .26852 .26719 m .5 .26719 L s .5 .26719 m .73148 .26719 L s .73148 .26719 m .96296 .26719 L s .96296 .26719 m .96296 .15 L s .96296 .15 m .96296 .03281 L s .5 .26719 m .5 .15 L s .5 .15 m .5 .03281 L s .03704 .03281 m matrix currentmatrix 0.0111111 0.009375 scale 3.3336 3.49973 1 0 365.73 arc setmatrix F .03704 .15 m matrix currentmatrix 0.0111111 0.009375 scale 3.3336 16 1 0 365.73 arc setmatrix F .03704 .26719 m matrix currentmatrix 0.0111111 0.009375 scale 3.3336 28.5003 1 0 365.73 arc setmatrix F .26852 .26719 m matrix currentmatrix 0.0111111 0.009375 scale 24.1668 28.5003 1 0 365.73 arc setmatrix F .5 .26719 m matrix currentmatrix 0.0111111 0.009375 scale 45 28.5003 1 0 365.73 arc setmatrix F .73148 .26719 m matrix currentmatrix 0.0111111 0.009375 scale 65.8332 28.5003 1 0 365.73 arc setmatrix F .96296 .26719 m matrix currentmatrix 0.0111111 0.009375 scale 86.6664 28.5003 1 0 365.73 arc setmatrix F .96296 .15 m matrix currentmatrix 0.0111111 0.009375 scale 86.6664 16 1 0 365.73 arc setmatrix F .96296 .03281 m matrix currentmatrix 0.0111111 0.009375 scale 86.6664 3.49973 1 0 365.73 arc setmatrix F .5 .15 m matrix currentmatrix 0.0111111 0.009375 scale 45 16 1 0 365.73 arc setmatrix F .5 .03281 m matrix currentmatrix 0.0111111 0.009375 scale 45 3.49973 1 0 365.73 arc setmatrix F .99 .99 .99 r .03704 .09141 m matrix currentmatrix 0.0277778 0.0234375 scale 1.33344 3.90016 1 0 365.73 arc setmatrix F .03704 .20859 m matrix currentmatrix 0.0277778 0.0234375 scale 1.33344 8.89984 1 0 365.73 arc setmatrix F .15278 .26719 m matrix currentmatrix 0.0277778 0.0234375 scale 5.50008 11.4001 1 0 365.73 arc setmatrix F .38426 .26719 m matrix currentmatrix 0.0277778 0.0234375 scale 13.8334 11.4001 1 0 365.73 arc setmatrix F .61574 .26719 m matrix currentmatrix 0.0277778 0.0234375 scale 22.1666 11.4001 1 0 365.73 arc setmatrix F .84722 .26719 m matrix currentmatrix 0.0277778 0.0234375 scale 30.4999 11.4001 1 0 365.73 arc setmatrix F .96296 .20859 m matrix currentmatrix 0.0277778 0.0234375 scale 34.6666 8.89984 1 0 365.73 arc setmatrix F .96296 .09141 m matrix currentmatrix 0.0277778 0.0234375 scale 34.6666 3.90016 1 0 365.73 arc setmatrix F .5 .20859 m matrix currentmatrix 0.0277778 0.0234375 scale 18 8.89984 1 0 365.73 arc setmatrix F .5 .09141 m matrix currentmatrix 0.0277778 0.0234375 scale 18 3.90016 1 0 365.73 arc setmatrix F 0 g 1 Mabswid newpath matrix currentmatrix 0.0277778 0.0234375 scale 1.33344 3.90016 1 0 365.73 arc setmatrix s newpath matrix currentmatrix 0.0277778 0.0234375 scale 1.33344 8.89984 1 0 365.73 arc setmatrix s newpath matrix currentmatrix 0.0277778 0.0234375 scale 5.50008 11.4001 1 0 365.73 arc setmatrix s newpath matrix currentmatrix 0.0277778 0.0234375 scale 13.8334 11.4001 1 0 365.73 arc setmatrix s newpath matrix currentmatrix 0.0277778 0.0234375 scale 22.1666 11.4001 1 0 365.73 arc setmatrix s newpath matrix currentmatrix 0.0277778 0.0234375 scale 30.4999 11.4001 1 0 365.73 arc setmatrix s newpath matrix currentmatrix 0.0277778 0.0234375 scale 34.6666 8.89984 1 0 365.73 arc setmatrix s newpath matrix currentmatrix 0.0277778 0.0234375 scale 34.6666 3.90016 1 0 365.73 arc setmatrix s newpath matrix currentmatrix 0.0277778 0.0234375 scale 18 8.89984 1 0 365.73 arc setmatrix s newpath matrix currentmatrix 0.0277778 0.0234375 scale 18 3.90016 1 0 365.73 arc setmatrix s gsave .03704 .09141 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (1) show 1.000 setlinewidth grestore gsave .03704 .20859 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (2) show 1.000 setlinewidth grestore gsave .15278 .26719 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (3) show 1.000 setlinewidth grestore gsave .38426 .26719 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (4) show 1.000 setlinewidth grestore gsave .61574 .26719 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (5) show 1.000 setlinewidth grestore gsave .84722 .26719 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (6) show 1.000 setlinewidth grestore gsave .96296 .20859 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (7) show 1.000 setlinewidth grestore gsave .96296 .09141 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (8) show 1.000 setlinewidth grestore gsave .5 .20859 -65.75 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (9) show 1.000 setlinewidth grestore gsave .5 .09141 -68.5 -10.5 Mabsadd m 1 1 Mabs scale currentpoint translate 0 21 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 11.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (10) show 1.000 setlinewidth grestore gsave .05926 .04781 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (1) show 1.000 setlinewidth grestore gsave .05926 .165 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (2) show 1.000 setlinewidth grestore gsave .05926 .28219 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (3) show 1.000 setlinewidth grestore gsave .29074 .28219 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (4) show 1.000 setlinewidth grestore gsave .52222 .28219 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (5) show 1.000 setlinewidth grestore gsave .7537 .28219 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (6) show 1.000 setlinewidth grestore gsave .98519 .28219 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (7) show 1.000 setlinewidth grestore gsave .98519 .165 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (8) show 1.000 setlinewidth grestore gsave .98519 .04781 -66 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (9) show 1.000 setlinewidth grestore gsave .52222 .165 -69 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (10) show 1.000 setlinewidth grestore gsave .52222 .04781 -69 -11 Mabsadd m 1 1 Mabs scale currentpoint translate 0 22 translate 1 -1 scale 63.000 14.000 moveto %%IncludeResource: font Times %%IncludeFont: Times /Times findfont 12.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 0.000 setrgbcolor (11) show 1.000 setlinewidth grestore % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 86.375}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgOn001@00Oomoogoo00001Woo00<007oPOn005goP00<0 07oPOn0077oP0Goo000YOn000`00On1oh00=On030003Ool01000Ool0000067oP00D007oPOn1oh000 00moh080009oo`80009oo`8001aoh003001oh7oP00eoh0800003Ool007oo00Aoo`8001Qoh005001o h7oPOn00000?On0200000goo00000003Ool2000KOn000`00On1oh00KOn01Ool002Moh080019oh0H0 01]oh004001oh7oP000AOn06000JOn04000BOn07000JOn04000BOn06000MOn000`00On1oh00KOn01 Ool002Uoh003001oh7oP035oh003001oh000035oh003001oh7oP03Aoh003001oh7oP03Moh003001o h7oP01Yoh05oo`009GoP00D007oPOn1oh00003Aoh080035oh003001oh7oP03Eoh003001oh7oP039o h005001oh7oPOn00000LOn01Ool002Ioh0<003Ioh003001oh7oP02moh0D003Aoh080039oh0D001ao h05oo`00ogoP87oP0Goo003oOn0POn01Ool00?moh21oh05oo`00ogoP87oP0Goo003oOn0POn01Ool0 02Uoh0D07`04On007`0O01l4On0401l01WoP01moh00OOn007`9oh0<07`Qoh00501moh00OOn007`02 On0301l2On0401l2On0601l7On000`0OOn007`02On0501l3On0201l2On0601l7On0301l2On0501l2 On0301l2On001@0OOn007goP01l027oP0P0O0WoP00<07goPOn0017oP0`0O00Aoh00O01moh0<07`=o h0@07`9oh0<07`Qoh0D07`9oh0@07`03On007`0O00807`=oh0<07`9oh0D07`03On007`0O00807bao h05oo`00:goP00<07goPOn000WoP00<07goPOn000WoP00@07goPOn007`9oh00701moh00OOn007goP 01l02goP00L07goP01moh00OOn007`09On00100OOn1oh00O0WoP00<07goPOn001WoP00<07goP01l0 17oP00<07goPOn000WoP00@07goPOn007`9oh00401moh7oP01l7On000`0OOn1oh005On001@0OOn1o h7oP01l01GoP00D07goP01moh00O00Moh00301moh7oP00807`Qoh00601moh7oP01moh00O0goP00<0 7goP01l00WoP00@07goPOn007`eoh00501moh7oPOn007`02On001@0OOn1oh7oP01l00WoP00<07goP 01l01goP00<07goPOn001GoP00<07goPOn00:GoP0Goo000[On000`0OOn1oh002On000`0OOn1oh003 On0301l2On001P0OOn007goP01moh0D07`Moh00601moh00OOn007goP1@0O0WoP0`0O0goP00@07goP On007`Moh00501moh00OOn007`03On000`0OOn1oh002On000`0OOn1oh003On00100OOn1oh00O1goP 1@0O0goP00@07goPOn1oh0D07`06On007goP01moh00O1goP00D07goPOn1oh00O00Qoh00601moh7oP 01moh00O0goP00<07goP01l00WoP00<07goPOn001@0O2GoP00<07goPOn000WoP0`0O0goP00D07goP On007goP00D07`=oh00301moh7oP009oh0<07baoh05oo`00:goP00<07goPOn000WoP0P0O0goP00@0 7goPOn007`9oh00701moh00OOn007goP01l00goP00<07goPOn001GoP00L07goP01moh00OOn007`03 On000`0OOn007`06On00100OOn1oh00O1goP00D07goP01moh00O00=oh00301moh7oP009oh00301mo h7oP00=oh00401moh7oP01l7On001@0OOn1oh7oP01l00goP00D07goPOn1oh00O00=oh00701moh00O On007goP01l027oP00D07goP01moh00O00Moh00601moh7oP01moh00O0goP00<07goP01l00WoP00@0 7goPOn007`=oh00301moh7oP00Moh00501moh7oPOn007`02On001@0OOn1oh7oP01l00WoP00<07goP 01l00goP00D07goPOn1oh00O00=oh00301moh7oP02eoh05oo`00:WoP100O00Aoh00O01moh0807`9o h0807`=oh0807`03On007goP009oh0<07`Qoh0807`03On007goP009oh0<07`=oh0@07`9oh0<07`Qo h0T07`=oh0@07`=oh0<07`Uoh0<07`Aoh00301moh7oP009oh0<07`9oh0807`03On007goP00Uoh003 01moh7oP00Moh0@07`=oh0<07`=oh0<07`=oh0<07`Yoh00301moh7oP009oh0807`Aoh0<07`=oh0<0 7`Aoh00301moh7oP009oh0@07b]oh05oo`00:goP00<07goPOn00On0301l9On0201l:On0301lcOn01Ool0 0?moh21oh05oo`00ogoP87oP0Goo0000\ \>"], ImageRangeCache->{{{0, 287}, {85.375, 0}} -> {-13201.2, -840.04, 91.9941, \ 109.03}}], Cell[BoxData[ \("Applied node forces="\)], "Print"], Cell[BoxData[ \({{0, 0, 0}, {0, 0, 0}, {2400, 0, 0}, {0, \(-4800\), 0}, {0, 0, 0}, {0, 0, 0}, {0, 0, 0}, {0, 0, 0}, {0, 0, 0}, {0, 0, 0}, {0, 0, 0}}\)], "Print"], Cell[BoxData[ InterpretationBox[\("pdof="\[InvisibleSpace]{1, 2, 3, 25, 26, 27, 31, 32, 33}\), SequenceForm[ "pdof=", {1, 2, 3, 25, 26, 27, 31, 32, 33}], Editable->False]], "Print"], Cell[BoxData[ \("Computed Nodal Displacements:"\)], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"0", "0", "0"}, { "0.0303661299298521`", \(-0.0006749364398764822`\), \ \(-0.000022309647797797032`\)}, { "0.14625125385415613`", \(-0.0013498728797529645`\), \ \(-0.000057012662662968005`\)}, {"0.14473638353426252`", \(-0.23184978289767522`\), "9.80846131423104`*^-6"}, {"0.1432215132143689`", \(-0.0018790617245506927`\), "0.000017461304099165202`"}, {"0.14283861873655787`", "0.052644001650203116`", "2.264789013239436`*^-6"}, { "0.1424557242587468`", \(-0.00006249396712490459`\), \ \(-0.000025430519497667474`\)}, { "0.05215497250612281`", \(-0.000031246983562452294`\), \ \(-0.000029256301190269837`\)}, {"0", "0", "0"}, { "0.08470673468155837`", \(-0.0009395308622753466`\), \ \(-0.00004017070432838353`\)}, {"0", "0", "0"} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Print"], Cell[BoxData[ \("External Node Forces Including Reactions:"\)], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {"251.02305981370523`", "1968.5646163064064`", "979086.3702150256`"}, {"0", "0", "0"}, {"2400.000000000063`", "0", "0"}, {"0", \(-4799.999999999991`\), "0"}, {"0", "0", "0"}, {"0", "0", "0"}, {"0", "0", "0"}, {"0", "0", "0"}, {\(-670.0653361692871`\), "91.13703539048585`", "2.7828246390793542`*^6"}, {"0", "0", "0"}, {\(-1980.9577236444627`\), "2740.298348303094`", "5.412364799864999`*^6"} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Print"], Cell[BoxData[ \("Internal Member Forces:"\)], "Print"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", GridBox[{ {\(-1968.5646163064066`\), \(-979086.3702150262`\), \ \(-1.7321555496561409`*^6\)}, {\(-1968.5646163064066`\), \(-1.7321555496561404`*^6\), \ \(-2.4852247290972746`*^6\)}, {\(-2651.023059813823`\), \(-2.4852247290972765`*^6\), "7.357598352434706`*^6"}, {\(-2651.023059813823`\), "7.35759835243471`*^6", \(-6.79957856603326`*^6\)}, {\(-670.0653361693121`\), \(-326197.02403148904`\), \ \(-781882.2009839318`\)}, {\(-670.0653361693606`\), \(-781882.2009839318`\), \ \(-1.2375673779363635`*^6\)}, {\(-91.13703539048586`\), \(-1.237567377936362`*^6\), "772628.6305714928`"}, {\(-91.13703539048586`\), "772628.630571494`", "2.7828246390793542`*^6"}, {\(-2740.2983483030926`\), \(-6.473381542001768`*^6\), \ \(-530508.3710683899`\)}, {\(-2740.2983483030935`\), \(-530508.3710683886`\), "5.412364799864998`*^6"} }], "\[NoBreak]", ")"}], Function[ BoxForm`e$, MatrixForm[ BoxForm`e$]]]], "Print"], Cell[GraphicsData["PostScript", "\<\ %! %%Creator: Mathematica %%AspectRatio: .36396 MathPictureStart /Mabs { Mgmatrix idtransform Mtmatrix dtransform } bind def /Mabsadd { Mabs 3 -1 roll add 3 1 roll add exch } bind def %% Graphics %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10 scalefont setfont % Background color 1 1 0 r MFill % Scaling calculations 0.50079 4.16122e-05 0.00866573 4.16122e-05 [ [.5 .37646 -101 0 ] [.5 .37646 101 12 ] [ 0 0 0 0 ] [ 1 .36396 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath 0 0 1 r gsave .5 .37646 -162 -4 Mabsadd m 1 1 Mabs scale currentpoint translate 0 20 translate 1 -1 scale 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 0.000 0.000 rmoveto 63.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor (Deformed) show 117.000 13.000 moveto (shape) show %%IncludeResource: font Math2Mono %%IncludeFont: Math2Mono /Math2Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 153.000 13.000 moveto (H) show 159.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor (magnified) show 219.000 13.000 moveto (x10000) show %%IncludeResource: font Math2Mono %%IncludeFont: Math2Mono /Math2Mono findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 255.000 13.000 moveto (L) show 261.000 13.000 moveto %%IncludeResource: font Courier %%IncludeFont: Courier /Courier findfont 10.000 scalefont [1 0 0 -1 0 0 ] makefont setfont 0.000 0.000 1.000 setrgbcolor 0.000 0.000 rmoveto 1.000 setlinewidth grestore 0 0 m 1 0 L 1 .36396 L 0 .36396 L closepath clip newpath 0 0 0 r 2 Mabswid [ ] 0 setdash .08467 .00867 m .08471 .01645 L s .08471 .01645 m .08483 .02424 L s .08483 .02424 m .08504 .03202 L s .08504 .03202 m .08534 .0398 L s .08534 .0398 m .08573 .04759 L s .08573 .04759 m .08622 .05537 L s .08622 .05537 m .08681 .06316 L s .08681 .06316 m .0875 .07094 L s .0875 .07094 m .08831 .07873 L s .08831 .07873 m .08923 .08651 L s .08923 .08651 m .09026 .0943 L s .09026 .0943 m .09141 .10208 L s .09141 .10208 m .09269 .10987 L s .09269 .10987 m .0941 .11765 L s .0941 .11765 m .09563 .12544 L s .09563 .12544 m .0973 .13322 L s .0973 .13322 m .09911 .14101 L s .09911 .14101 m .10107 .14879 L s .10107 .14879 m .10317 .15658 L s .10317 .15658 m .10542 .16436 L s .10542 .16436 m .10782 .17215 L s .10782 .17215 m .11039 .17993 L s .11039 .17993 m .11311 .18771 L s .11311 .18771 m .116 .1955 L s .116 .1955 m .11906 .20328 L s .11906 .20328 m .12229 .21107 L s .12229 .21107 m .1257 .21885 L s .1257 .21885 m .12929 .22664 L s .12929 .22664 m .13306 .23442 L s .13306 .23442 m .13702 .24221 L s .13702 .24221 m .14118 .24999 L s .14118 .24999 m .14553 .25778 L s .14553 .25778 m .15849 .25011 L s .15849 .25011 m .17146 .24202 L s .17146 .24202 m .18442 .23366 L s .18442 .23366 m .19738 .22515 L s .19738 .22515 m .21035 .21664 L s .21035 .21664 m .22331 .20826 L s .22331 .20826 m .23628 .20015 L s .23628 .20015 m .24924 .19244 L s .24924 .19244 m .26221 .18528 L s .26221 .18528 m .27517 .17879 L s .27517 .17879 m .28813 .17313 L s .28813 .17313 m .3011 .16842 L s .3011 .16842 m .31406 .1648 L s .31406 .1648 m .32703 .16241 L s .32703 .16241 m .33999 .16138 L s .33999 .16138 m .35296 .16186 L s .35296 .16186 m .36592 .16392 L s .36592 .16392 m .37889 .16743 L s .37889 .16743 m .39185 .17218 L s .39185 .17218 m .40481 .17798 L s .40481 .17798 m .41778 .18463 L s .41778 .18463 m .43074 .19194 L s .43074 .19194 m .44371 .1997 L s .44371 .1997 m .45667 .20772 L s .45667 .20772 m .46964 .2158 L s .46964 .2158 m .4826 .22375 L s .4826 .22375 m .49557 .23137 L s .49557 .23137 m .50853 .23845 L s .50853 .23845 m .52149 .24481 L s .52149 .24481 m .53446 .25025 L s .53446 .25025 m .54742 .25456 L s .54742 .25456 m .56039 .25756 L s .56039 .25756 m .57338 .25979 L s .57338 .25979 m .58637 .26194 L s .58637 .26194 m .59937 .26401 L s .59937 .26401 m .61236 .26599 L s .61236 .26599 m .62536 .26787 L s .62536 .26787 m .63835 .26964 L s .63835 .26964 m .65134 .27131 L s .65134 .27131 m .66434 .27285 L s .66434 .27285 m .67733 .27428 L s .67733 .27428 m .69033 .27557 L s .69033 .27557 m .70332 .27673 L s .70332 .27673 m .71631 .27774 L s .71631 .27774 m .72931 .27861 L s .72931 .27861 m .7423 .27932 L s .7423 .27932 m .7553 .27986 L s .7553 .27986 m .76829 .28025 L s .76829 .28025 m .78128 .28045 L s .78128 .28045 m .79428 .28048 L s .79428 .28048 m .80727 .28032 L s .80727 .28032 m .82026 .27996 L s .82026 .27996 m .83326 .27941 L s .83326 .27941 m .84625 .27865 L s .84625 .27865 m .85925 .27767 L s .85925 .27767 m .87224 .27648 L s .87224 .27648 m .88523 .27507 L s .88523 .27507 m .89823 .27342 L s .89823 .27342 m .91122 .27153 L s .91122 .27153 m .92422 .2694 L s .92422 .2694 m .93721 .26702 L s .93721 .26702 m .9502 .26438 L s .9502 .26438 m .9632 .26148 L s .9632 .26148 m .97619 .25831 L s .97619 .25831 m .97416 .25051 L s .97416 .25051 m .97204 .24271 L s .97204 .24271 m .96984 .23491 L s .96984 .23491 m .96757 .22711 L s .96757 .22711 m .96524 .21931 L s .96524 .21931 m .96286 .2115 L s .96286 .2115 m .96044 .2037 L s .96044 .2037 m .958 .1959 L s .958 .1959 m .95554 .1881 L s .95554 .1881 m .95306 .1803 L s .95306 .1803 m .95059 .1725 L s .95059 .1725 m .94813 .1647 L s .94813 .1647 m .9457 .15689 L s .9457 .15689 m .94329 .14909 L s .94329 .14909 m .94093 .14129 L s .94093 .14129 m .93861 .13349 L s .93861 .13349 m .93636 .12569 L s .93636 .12569 m .93419 .11789 L s .93419 .11789 m .93209 .11008 L s .93209 .11008 m .93009 .10228 L s .93009 .10228 m .92819 .09448 L s .92819 .09448 m .9264 .08668 L s .9264 .08668 m .92473 .07888 L s .92473 .07888 m .9232 .07108 L s .9232 .07108 m .92181 .06328 L s .92181 .06328 m .92057 .05547 L s .92057 .05547 m .91949 .04767 L s .91949 .04767 m .91859 .03987 L s .91859 .03987 m .91787 .03207 L s .91787 .03207 m .91734 .02427 L s .91734 .02427 m .91702 .01647 L s .91702 .01647 m .91691 .00867 L s .56039 .25756 m .56149 .24978 L s .56149 .24978 m .56211 .242 L s .56211 .242 m .56227 .23422 L s .56227 .23422 m .562 .22645 L s .562 .22645 m .56133 .21867 L s .56133 .21867 m .56028 .21089 L s .56028 .21089 m .5589 .20311 L s .5589 .20311 m .55721 .19533 L s .55721 .19533 m .55523 .18756 L s .55523 .18756 m .553 .17978 L s .553 .17978 m .55056 .172 L s .55056 .172 m .54792 .16422 L s .54792 .16422 m .54512 .15644 L s .54512 .15644 m .54218 .14867 L s .54218 .14867 m .53915 .14089 L s .53915 .14089 m .53604 .13311 L s .53604 .13311 m .53289 .12533 L s .53289 .12533 m .52972 .11756 L s .52972 .11756 m .52658 .10978 L s .52658 .10978 m .52348 .102 L s .52348 .102 m .52046 .09422 L s .52046 .09422 m .51754 .08644 L s .51754 .08644 m .51476 .07867 L s .51476 .07867 m .51215 .07089 L s .51215 .07089 m .50973 .06311 L s .50973 .06311 m .50753 .05533 L s .50753 .05533 m .5056 .04755 L s .5056 .04755 m .50395 .03978 L s .50395 .03978 m .50261 .032 L s .50261 .032 m .50162 .02422 L s .50162 .02422 m .501 .01644 L s .501 .01644 m .50079 .00867 L s 1 1 1 r 1 Mabswid .02381 .00867 Mdot .97619 .35529 Mdot % End of Graphics MathPictureEnd \ \>"], "Graphics", ImageSize->{288, 104.812}, ImageMargins->{{43, 0}, {0, 0}}, ImageRegion->{{0, 1}, {0, 1}}, ImageCache->GraphicsData["Bitmap", "\<\ CF5dJ6E]HGAYHf4PAg9QL6QYHgWoP0Goo003oOn0POn01Ool00?moh21oh05oo`00ogoP87oP 0Goo003oOn0POn01Ool00?moh21oh05oo`00ogoP87oP0Goo003oOn0POn01Ool00?moh21oh05oo`00 ogoP87oP0Goo003oOn0POn01Ool00?moh21oh05oo`00ogoP87oP0Goo003oOn0POn01Ool00?moh21o h05oo`00ogoP87oP0Goo003oOn0POn01Ool00?moh0Aoh003Oomoh7oP01Uoh05oo`00ogoP87oP0Goo 003oOn0POn01Ool00?moh21oh05oo`00ogoP87oP0Goo003oOn0POn01Ool00?moh21oh05oo`00M7oP 0`0O4goP00<07goPOn0037oP0`0ODGoP00<07goPOn00;goP0Goo001eOn000`0OOn1oh00AOn000`0O On1oh00@On000`0OOn1oh01?On000`0OOn1oh00^On01Ool002aoh0@07`=oh0<07`9oh0D07`9oh0<0 7`9oh0<07`=oh00501moh00OOn007`02On0301l3On0401l7On0401l2On0601l00goP01l07`0201l2 On0301l3On0301l:On002@0OOn1oh7oP01moh00OOn007`02On0401l2On0401l01GoP01l07`0OOn00 1`0O00=oh00O01l00`0O00=oh00O01l00`0O0WoP0`0O0goP100O1goP0P0O00=oh00O01l00WoP0`0O 17oP0P0O17oP0P0O17oP0P0O17oP0P0O1GoP00<07goPOn00;GoP0Goo000]On001P0OOn1oh00OOn00 7`Moh00501moh7oPOn007`03On00100OOn1oh00O17oP00L07goP01moh00OOn007`05On00100OOn1o h00O37oP00@07goPOn007`9oh00301moh00O009oh00501moh7oPOn007`02On000`0OOn007`0=On00 2`0OOn1oh7oP01moh00OOn007goP01l00WoP00@07goPOn007`=oh00401moh7oP01l2On001@0OOn1o h7oP01l01GoP00<07goPOn000goP00D07goPOn1oh00O00Eoh00401moh7oP01l9On000`0OOn007`04 On000`0OOn1oh002On00100OOn1oh00O0WoP00@07goPOn007`9oh00401moh7oP01l2On00100OOn1o h00O17oP00<07goPOn00;GoP0Goo000]On001@0OOn1oh00OOn001@0O0goP00D07goPOn1oh00O00=o h00401moh7oP01l4On001P0OOn007goP01moh0D07`05On007goPOn007`09On0301l3On00100OOn1o h00O0WoP0`0O0goP00D07goPOn007goP00D07`Uoh00901moh7oPOn007goP01moh00O009oh0<07`9o h00501moh7oPOn007`02On00100OOn1oh00O0goP00<07goPOn000goP00<07goPOn000goP00@07goP On1oh0D07`05On007goPOn007`0:On000`0OOn1oh003On000`0OOn1oh002On00100OOn1oh00O0WoP 00@07goPOn007`9oh00401moh7oP01l2On00100OOn1oh00O17oP00<07goPOn00;GoP0Goo000]On00 1P0OOn1oh00OOn007`=oh00501moh7oPOn007`03On001@0OOn1oh7oP01l00WoP0P0O0goP00L07goP 01moh00OOn007`03On000`0OOn007`02On000`0OOn1oh006On000`0OOn1oh004On001P0OOn1oh00O On007`9oh00501moh7oPOn007`02On000`0OOn007`03On000`0OOn1oh007On002`0OOn1oh7oP01mo h00OOn007goP01l00WoP00@07goPOn007`=oh00401moh7oP01l2On001@0OOn1oh7oP01l01GoP00<0 7goPOn000goP00D07goPOn1oh00O00=oh00301moh00O009oh00301moh7oP00Moh00301moh00O00Ao h00301moh7oP009oh00401moh7oP01l2On00100OOn1oh00O0WoP00@07goPOn007`9oh00401moh7oP 01l4On000`0OOn1oh00]On01Ool002eoh00401moh7oP01l2On0301l3On0401l2On0301l2On0201l0 17oP01l07goP0P0O00=oh00OOn000WoP0`0O0goP0`0O2GoP100O0WoP0`0O0goP0P0O0goP100O0goP 0`0O2goP00<07goPOn000P0O00=oh00OOn000WoP0P0O17oP2@0O0WoP0`0O17oP100O00Aoh00O01l0 7`Aoh0<07`=oh0<07`Qoh0807`03On007`0O00=oh00301moh7oP009oh00401moh7oP01l2On00100O On1oh00O0WoP00@07goPOn007`9oh00401moh7oP01l3On000`0OOn1oh00^On01Ool002eoh00401mo h7oP01l9On000`0OOn1oh00LOn000`0OOn1oh00=On000`0OOn1oh00NOn000`0OOn1oh00OOn000`0O On1oh00@On000`0OOn1oh00=On0201l4On00100OOn1oh00O0WoP00@07goPOn007`9oh00401moh7oP 01l2On00100OOn1oh00O0WoP00<07goPOn00;goP0Goo000/On0401l;On0201lKOn0201l>On0201ll On000`0OOn1oh004On0201l3On000`0OOn1oh009On0201l@On000`0OOn1oh003On0201l4On0201l4 On0201l4On0201leOn01Ool00?moh21oh05oo`00ogoP87oP0Goo0000\ \>"], ImageRangeCache->{{{0, 287}, {103.812, 0}} -> {-14159.1, -208.295, 98.5373, \ 98.5373}}], Cell["\<\ Cell 6A. Animation of deflections. Uses the displacements u \ computed by driver program to plot a sequence of deformed shapes with amplitude varying as sine of \ pseudotime t. To animate, click twice on a plot, then use the VCR buttons that appear on \ the bottom of the window to regulate frame speed.\ \>", "Text", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 1]], Cell["\<\ ampmax=10000.; ampmin=0; For [t=0.,t<=N[Pi],t=t+N[Pi/6], amp=ampmax*Sin[t]; PlotPlaneFrameDeformedShape[NodeCoordinates,ElemNodeLists, u,{amp,ampmax,ampmin,16},-1,\"deformed shape\"]]; \ \>", "Input", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 0]], Cell["\<\ ClearAll[L,H,a,Em,P]; R=10; Em=10000; b=1; h=5/10; A=b*h; Izz=b*h^3/12; P=4; ne=32; nnod=ne+1; NodeCoordinates=Table[{0,0},{nnod}]; For [n=1,n<=nnod,n++, arg=Pi/(2*ne)*(n-1); NodeCoordinates[[n]]=R*{Sin[arg],Cos[arg]}]; Print[\"NodeCoordinates=\",NodeCoordinates]; ElemNodeLists=Table[{0,0},{ne}]; For [e=1,e<=ne,e++,ElemNodeLists[[e]]={e,e+1}]; numnod=Length[NodeCoordinates]; numele=Length[ElemNodeLists]; numdof=3*numnod; ElemMaterial= Table[{Em,0,0,0},{numele}]; ElemFabrication= Table[{A,Izz},{numele}]; ProcessOptions= {True}; aspect=0; PlotLineElements[NodeCoordinates,ElemNodeLists,aspect, \"frame mesh\"]; PlotLineElementsAndNodes[NodeCoordinates,ElemNodeLists,aspect, \"frame mesh with elem & node labels\",{True,0.03},{True,0.02}]; FreedomTag=FreedomValue=Table[{0,0,0},{numnod}]; FreedomValue[[1]]={0,-P,0}; Print[\"Applied node forces=\"]; Print[FreedomValue]; FreedomTag[[1]]={1,0,1}; FreedomTag[[numnod]]={0,1,1}; f=Flatten[FreedomValue]; K=PlaneFrameMasterStiffness[NodeCoordinates, ElemNodeLists,ElemMaterial,ElemFabrication,ProcessOptions]; Print[\"K=\",K//MatrixForm]; pdof={}; For[n=1,n<=numnod,n++, For[j=1,j<=3,j++, If [FreedomTag[[n,j]]>0, AppendTo[pdof,3*(n-1)+j]]]]; Print[\"pdof=\",pdof]; Kmod=ModifiedMasterStiffness[pdof,K]; Print[\"Kmod=\",Kmod//MatrixForm]; Print[\"eigs of Kmod=\",Chop[Eigenvalues[N[Kmod]]]]; fmod=ModifyNodeForces [pdof,f]; u=LinearSolve[Kmod,fmod]; u=Chop[u,.000001]; Print[\"Computed Nodal Displacements:\"]; NodeForces=NodeDisplacements=Table[{0,0,0},{numnod}]; For [n=1,n<=numnod,n++, For[j=1,j<=3,j++, NodeDisplacements[[n,j]]=u[[3*(n-1)+j]] ]]; Print[NodeDisplacements//MatrixForm]; f=Simplify[K.u]; f=Chop[f,.000001]; Print[\"External Node Forces Including Reactions:\"]; For [n=1,n<=numnod,n++, For[j=1,j<=3,j++, NodeForces[[n,j]]=f[[3*(n-1)+j]] ]]; Print[NodeForces//MatrixForm]; p=PlaneFrameIntForces[NodeCoordinates,ElemNodeLists, ElemMaterial,ElemFabrication,eleopt,u]; p=Chop[p,.000001]; Print[\"Internal Member Forces:\"]; Print[p//MatrixForm]; amp=1.; PlotPlaneFrameDeformedShape[NodeCoordinates,ElemNodeLists,u, {amp,amp,0,16},aspect,\"Deformed shape (magnified x10000)\"]; \ \>", "Input", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 0]] }, Open ]], Cell["\<\ ClearAll[L,H,a,Em,P]; NodeCoordinates={{0,a},{a,0}}; numnod=2; numele=1; Print[\"NodeCoordinates=\",NodeCoordinates]; ElemNodeLists={{1,2}}; ElemMaterial= {{Em,0,0,0}}; ElemFabrication= {{A,\[Alpha]*a^2*A}}; ProcessOptions= {False}; (*aspect=0; PlotLineElements[NodeCoordinates,ElemNodeLists,aspect, \"frame mesh\"]; PlotLineElementsAndNodes[NodeCoordinates,ElemNodeLists,aspect, \"frame mesh with elem & node labels\",{True,0.03},{True,0.02}];*) FreedomTag={{1,0,1},{0,1,1}}; FreedomValue={{0,P,0},{0,0,0}}; f=Flatten[FreedomValue]; K=PlaneFrameMasterStiffness[NodeCoordinates, ElemNodeLists,ElemMaterial,ElemFabrication,ProcessOptions]; Print[\"K=\",K//MatrixForm]; pdof={}; For[n=1,n<=numnod,n++, For[j=1,j<=3,j++, If [FreedomTag[[n,j]]>0, AppendTo[pdof,3*(n-1)+j]]]]; Print[\"pdof=\",pdof]; Kmod=ModifiedMasterStiffness[pdof,K]; Print[\"Kmod=\",Kmod//MatrixForm]; (*Print[\"eigs of Kmod=\",Chop[Eigenvalues[N[Kmod]]]];*) fmod=ModifyNodeForces [pdof,f]; u=LinearSolve[Kmod,fmod]; u=Simplify[u]; Print[\"Computed Nodal Displacements:\"]; NodeForces=NodeDisplacements=Table[{0,0,0},{numnod}]; For [n=1,n<=numnod,n++, For[j=1,j<=3,j++, NodeDisplacements[[n,j]]=u[[3*(n-1)+j]] ]]; Print[NodeDisplacements//MatrixForm]; f=Simplify[K.u]; f=Chop[f,.000001]; Print[\"External Node Forces Including Reactions:\"]; For [n=1,n<=numnod,n++, For[j=1,j<=3,j++, NodeForces[[n,j]]=f[[3*(n-1)+j]] ]]; Print[NodeForces//MatrixForm]; p=PlaneFrameIntForces[NodeCoordinates,ElemNodeLists, ElemMaterial,ElemFabrication,eleopt,u]; p=Simplify[p]; Print[\"Internal Member Forces:\"]; Print[p//MatrixForm];\ \>", "Input", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 0]], Cell["UVA plane - 3-member plane frame model of wings", "Text", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[1, 1, 0]], Cell["\<\ ClearAll[L,H,a,b,Em,P,A,Izz]; Em=185648; a=30; b=9; A=0.7449; \ Izz=0.4434; P=2*4.4723; NodeCoordinates={{0,b},{a,b},{a,0},{0,0} }; numnod=Length[NodeCoordinates]; Print[\"NodeCoordinates=\",NodeCoordinates]; ElemNodeLists={{1,2},{2,3},{3,4}}; numele=Length[ElemNodeLists]; \ Print[numele]; ElemMaterial= Table[{Em,0,0,0},{numele}]; ElemFabrication= Table[{A,Izz},{numele}]; ProcessOptions= {False}; aspect=0; PlotLineElements[NodeCoordinates,ElemNodeLists,aspect, \"UVA mesh\"]; PlotLineElementsAndNodes[NodeCoordinates,ElemNodeLists,aspect, \"UVA mesh with elem & node labels\",{True,0.12},{True,0.08}]; FreedomTag={{1,1,1},{0,0,0},{0,0,0},{1,1,1}}; FreedomValue={{0,0,0},{0,P/2,0},{0,P/2,0},{0,0,0}}; f=Flatten[FreedomValue]; K=PlaneFrameMasterStiffness[NodeCoordinates, ElemNodeLists,ElemMaterial,ElemFabrication,ProcessOptions]; K=N[K]; Print[\"K=\",K//MatrixForm]; Print[\"eigs of K=\",Chop[Eigenvalues[N[K]]]]; pdof={}; For[n=1,n<=numnod,n++, For[j=1,j<=3,j++, If [FreedomTag[[n,j]]>0, AppendTo[pdof,3*(n-1)+j]]]]; Print[\"pdof=\",pdof]; Kmod=ModifiedMasterStiffness[pdof,K]; Print[\"Kmod=\",Kmod//MatrixForm]; Print[\"eigs of Kmod=\",Chop[Eigenvalues[N[Kmod]]]]; fmod=ModifyNodeForces [pdof,f]; u=LinearSolve[Kmod,fmod]; u=Simplify[u]; Print[\"Computed Nodal Displacements:\"]; NodeForces=NodeDisplacements=Table[{0,0,0},{numnod}]; For [n=1,n<=numnod,n++, For[j=1,j<=3,j++, NodeDisplacements[[n,j]]=u[[3*(n-1)+j]] ]]; Print[NodeDisplacements//MatrixForm]; f=Simplify[K.u]; f=Chop[f,.000001]; Print[\"External Node Forces Including Reactions:\"]; For [n=1,n<=numnod,n++, For[j=1,j<=3,j++, NodeForces[[n,j]]=f[[3*(n-1)+j]] ]]; Print[NodeForces//MatrixForm]; p=PlaneFrameIntForces[NodeCoordinates,ElemNodeLists, ElemMaterial,ElemFabrication,eleopt,u]; p=Simplify[p]; Print[\"Internal Member Forces:\"]; Print[p//MatrixForm]; amp=20.; PlotPlaneFrameDeformedShape[NodeCoordinates,ElemNodeLists,u, {amp,amp,0,16},aspect,\"Deformed shape (magnified x1)\"]; H=2; ElemFabrication= Table[{A,Izz,H},{numele}]; fmax=200; Nsub={16,8}; aspect=-1; title=\"Axial stress in frame members\"; PlotPlaneFrameStress[NodeCoordinates,ElemNodeLists, ElemFabrication,p,fmax,Nsub,aspect,title]; \ \>", "Input", CellFrame->True, CellMargins->{{12, 51}, {Inherited, Inherited}}, CellLabelMargins->{{6, Inherited}, {Inherited, Inherited}}, ImageRegion->{{-0, 1}, {0, 1}}, Background->RGBColor[0, 1, 0]] }, FrontEndVersion->"4.2 for Macintosh", ScreenRectangle->{{0, 1920}, {0, 1180}}, AutoGeneratedPackage->None, WindowToolbars->{}, CellGrouping->Manual, WindowSize->{1436, 1119}, WindowMargins->{{72, Automatic}, {Automatic, 5}}, PrivateNotebookOptions->{"ColorPalette"->{RGBColor, -1}}, ShowCellLabel->True, ShowCellTags->False, RenderingOptions->{"ObjectDithering"->True, "RasterDithering"->False}, Magnification->1.5, MacintoshSystemPageSetup->"\<\ 00<0001804P000000]P2:?oQon82n@960dL5:0?l0080001804P000000]P2:001 0000I00000400`<300000BL?00400@00000000000000060801T1T00000000000 00000000000000000000000000000000\>" ] (******************************************************************* Cached data follows. If you edit this Notebook file directly, not using Mathematica, you must remove the line containing CacheID at the top of the file. The cache data will then be recreated when you save this file from within Mathematica. *******************************************************************) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[1754, 51, 1601, 31, 463, "Text"], Cell[3358, 84, 479, 12, 117, "Text"], Cell[CellGroupData[{ Cell[3862, 100, 2392, 44, 942, "Input", InitializationCell->True], Cell[6257, 146, 1122, 20, 146, "Print"], Cell[7382, 168, 1114, 20, 146, "Print"], Cell[8499, 190, 2526, 45, 146, "Print"], Cell[11028, 237, 69, 1, 32, "Print"], Cell[11100, 240, 1829, 34, 214, "Print"], Cell[12932, 276, 179, 3, 32, "Print"] }, Open ]], Cell[13126, 282, 451, 11, 117, "Text"], Cell[CellGroupData[{ Cell[13602, 297, 1160, 21, 542, "Input", InitializationCell->True], Cell[14765, 320, 62, 1, 32, "Print"], Cell[14830, 323, 904, 17, 146, "Print"], Cell[15737, 342, 83, 1, 32, "Print"], Cell[15823, 345, 562, 11, 146, "Print"], Cell[16388, 358, 53, 1, 32, "Print"], Cell[16444, 361, 75, 1, 32, "Print"], Cell[16522, 364, 84, 1, 32, "Print"], Cell[16609, 367, 63, 1, 32, "Print"] }, Open ]], Cell[16687, 371, 353, 9, 93, "Text"], Cell[CellGroupData[{ Cell[17065, 384, 2076, 38, 762, "Input", InitializationCell->True], Cell[19144, 424, 63, 1, 32, "Print"], Cell[19210, 427, 667, 14, 120, "Print"], Cell[19880, 443, 432, 11, 117, "Text"], Cell[20315, 456, 301, 8, 69, "Text"], Cell[CellGroupData[{ Cell[20641, 468, 1590, 28, 642, "Input", InitializationCell->True], Cell[22234, 498, 2566, 132, 274, 2308, 125, "GraphicsData", "PostScript", \ "Graphics", ImageCacheValid->False] }, Open ]], Cell[24815, 633, 311, 8, 69, "Text"], Cell[CellGroupData[{ Cell[25151, 645, 2941, 50, 1182, "Input", InitializationCell->True], Cell[28095, 697, 23529, 962, 319, 23270, 955, "GraphicsData", "PostScript", \ "Graphics", ImageCacheValid->False] }, Open ]], Cell[51639, 1662, 492, 11, 119, "Text"], Cell[52134, 1675, 393, 8, 82, "Input", InitializationCell->True], Cell[52530, 1685, 370, 9, 93, "Text"], Cell[52903, 1696, 268, 5, 69, "Text"], Cell[CellGroupData[{ Cell[53196, 1705, 4601, 79, 1902, "Input", InitializationCell->True], Cell[57800, 1786, 5082, 350, 275, 4823, 343, "GraphicsData", "PostScript", \ "Graphics", ImageCacheValid->False] }, Open ]], Cell[CellGroupData[{ Cell[62919, 2141, 3953, 70, 1782, "Input", InitializationCell->True], Cell[66875, 2213, 47873, 3558, 405, 47622, 3551, "GraphicsData", \ "PostScript", "Graphics", ImageCacheValid->False] }, Open ]], Cell[114763, 5774, 394, 10, 93, "Text"], Cell[CellGroupData[{ Cell[115182, 5788, 2666, 47, 1142, "Input"], Cell[117851, 5837, 6968, 196, 117, 1554, 123, "GraphicsData", "PostScript", \ "Graphics"] }, Open ]], Cell[124834, 6036, 308, 8, 69, "Text"], Cell[125145, 6046, 2538, 69, 1302, "Input"], Cell[127686, 6117, 5290, 154, 142, 1717, 105, "GraphicsData", "PostScript", \ "Graphics"], Cell[132979, 6273, 21009, 695, 142, 11833, 577, "GraphicsData", "PostScript", \ "Graphics"], Cell[153991, 6970, 55, 1, 32, "Print"], Cell[154049, 6973, 183, 3, 32, "Print"], Cell[154235, 6978, 208, 4, 32, "Print"], Cell[154446, 6984, 64, 1, 32, "Print"], Cell[154513, 6987, 1158, 29, 270, "Print"], Cell[155674, 7018, 76, 1, 32, "Print"], Cell[155753, 7021, 748, 19, 270, "Print"], Cell[156504, 7042, 58, 1, 32, "Print"], Cell[156565, 7045, 1179, 25, 282, "Print"], Cell[157747, 7072, 13511, 661, 170, 7466, 582, "GraphicsData", "PostScript", \ "Graphics"], Cell[171261, 7735, 529, 13, 141, "Text"], Cell[171793, 7750, 421, 11, 142, "Input"], Cell[172217, 7763, 2419, 66, 1242, "Input"] }, Open ]], Cell[174651, 7832, 1846, 52, 962, "Input"], Cell[176500, 7886, 262, 5, 69, "Text"], Cell[176765, 7893, 2459, 67, 1222, "Input"] } ] *) (******************************************************************* End of Mathematica Notebook file. *******************************************************************)