(************** Content-type: application/mathematica **************

                    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[     41548,       1293]*)
(*NotebookOutlinePosition[     42555,       1324]*)
(*  CellTagsIndexPosition[     42511,       1320]*)
(*WindowFrame->Normal*)



Notebook[{
Cell[TextData[StyleBox["Problem 3. Using  thin-walled assumption.",
  FontFamily->"Times New Roman"]], "Title",
  CellFrame->{{0, 0}, {0, 3}},
  FontFamily->"Courier"],

Cell[BoxData[
    \(Off[General::spell1, General::spell, Solve::"\<ifun\>"]\)], "Input"],

Cell[CellGroupData[{

Cell["Cross section properties (all in mm)", "Section",
  FontFamily->"Times New Roman"],

Cell[BoxData[{
    \(\ \(\[Beta] = \ 
        t\_2\/t\_1\  /. \ {t\_2\  \[Rule] \ 12, \ \ t\_1\  \[Rule] \ 
              10};\)\), "\[IndentingNewLine]", 
    \(\(t\_1 = t\_2\/\[Beta];\)\), "\[IndentingNewLine]", 
    \(\(L = 3000;\)\), "\[IndentingNewLine]", 
    \(\(vars = {\ b \[Rule] \ 200, h \[Rule] \ 200, 
          t\_2\  \[Rule] \ 12};\)\)}], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Total Cross sectional Area of the beam", "Section",
  FontFamily->"Times New Roman"],

Cell["Section 1: Web ", "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(h\_1\  = h;\)\), "\[IndentingNewLine]", 
    \(\(A\_1\  = \ h\_1\ t\_1;\)\), "\[IndentingNewLine]", 
    \(Print[\*"\"\<The redefined web area is \!\(A\_1\) = \>\"", 
      A\_1, "\< = \>", \ 
      A\_1\  /. vars\  // N, \*"\"\< \!\(mm\^2\)\>\""]\)}], "Input"],

Cell[BoxData[
    InterpretationBox[\("The redefined web area is \!\(A\_1\) = "\
\[InvisibleSpace]\(5\ h\ t\_2\)\/6\[InvisibleSpace]" = \
"\[InvisibleSpace]2000.`\[InvisibleSpace]" \!\(mm\^2\)"\),
      SequenceForm[ "The redefined web area is \!\(A\_1\) = ", 
        Times[ 
          Rational[ 5, 6], h, 
          Subscript[ t, 2]], " = ", .2*^4, " \!\(mm\^2\)"],
      Editable->False]], "Print"]
}, Open  ]],

Cell["Section 2: Flange", "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(h\_2\  = t\_2;\)\), "\[IndentingNewLine]", 
    \(\(A\_2\  = \ b\ h\_2\ ;\)\), "\[IndentingNewLine]", 
    \(Print[\*"\"\<The redefined flange area is \!\(A\_2\) = \>\"", 
      A\_2, "\< = \>", \ 
      A\_2 /. vars\  // N, \*"\"\< \!\(mm\^2\)\>\""]\)}], "Input"],

Cell[BoxData[
    InterpretationBox[\("The redefined flange area is \!\(A\_2\) = "\
\[InvisibleSpace]\(b\ t\_2\)\[InvisibleSpace]" = "\[InvisibleSpace]2400.`\
\[InvisibleSpace]" \!\(mm\^2\)"\),
      SequenceForm[ "The redefined flange area is \!\(A\_2\) = ", 
        Times[ b, 
          Subscript[ t, 2]], " = ", .24*^4, " \!\(mm\^2\)"],
      Editable->False]], "Print"]
}, Open  ]],

Cell["Total area", "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(A\_\(\(Tp\)\(\ \)\) = \ 
      ExpandAll[\[Sum]\+\(j = 1\)\%2 A\_j]\ \), "\[IndentingNewLine]", 
    \(Print["\<For thin-walls neglect all higher order terms of thickness\>"]\
\), "\[IndentingNewLine]", 
    \(\(A\_\(\(T\)\(\ \)\) = 
        Normal[Series[
            A\_\(\(Tp\)\(\ \)\), {t\_2, 0, 1}]];\)\), "\[IndentingNewLine]", 
    \(Print[\*"\"\<The total area is \!\(A\_\(\(T\)\(\\\ \)\)\) = \>\"", 
      Expand[A\_\(\(T\)\(\ \)\)\ ], "\< = \>", \ 
      A\_\(\(T\)\(\ \)\)\  /. vars // 
        N, \*"\"\< \!\(mm\^2\)\>\""]\)}], "Input"],

Cell[BoxData[
    \(b\ t\_2 + \(5\ h\ t\_2\)\/6\)], "Output"],

Cell[BoxData[
    \("For thin-walls neglect all higher order terms of thickness"\)], "Print"],

Cell[BoxData[
    InterpretationBox[\("The total area is \!\(A\_\(\(T\)\(\\ \)\)\) = "\
\[InvisibleSpace]\(b\ t\_2 + \(5\ h\ t\_2\)\/6\)\[InvisibleSpace]" = "\
\[InvisibleSpace]4400.`\[InvisibleSpace]" \!\(mm\^2\)"\),
      SequenceForm[ "The total area is \!\(A\_\(\(T\)\(\\ \)\)\) = ", 
        Plus[ 
          Times[ b, 
            Subscript[ t, 2]], 
          Times[ 
            Rational[ 5, 6], h, 
            Subscript[ t, 2]]], " = ", .44*^4, " \!\(mm\^2\)"],
      Editable->False]], "Print"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
  "Centroid: ",
  Cell[BoxData[
      FormBox[
        SubscriptBox[
          FormBox[\(x\&_\),
            "TraditionalForm"], "G"], TraditionalForm]]],
  " (let's place the origin at point B)"
}], "Section",
  FontFamily->"Times New Roman"],

Cell[TextData[{
  "Centroid: ",
  Cell[BoxData[
      \(\(\(x\&_\_G\)\(=\)\(\(Sum\ of\ all\ \(A\_j\) 
              x\&_\_j\)\/\(Total\ Area\)\)\(\ \)\)\)]]
}], "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(x\&_\_1 = \ 0;\)\), "\[IndentingNewLine]", 
    \(\(x\&_\_2 = \ b\/2\ ;\)\), "\[IndentingNewLine]", 
    \(\(Qx = 
        ExpandAll[\ \((\[Sum]\+\(j = 1\)\%2\( A\_j\) 
                x\&_\_j\ )\)];\)\), "\[IndentingNewLine]", 
    \(Print["\<The first moment of area about the B'C' is:\>"]\), "\
\[IndentingNewLine]", 
    \(ExpandAll[\ \((Qx\ )\)]\), "\[IndentingNewLine]", 
    \(Print["\<The centroid is located at \>"]\), "\[IndentingNewLine]", 
    \(x\&_\_Gp = 
      FullSimplify[ExpandAll[\ \((Qx\/A\_T\ )\)]]\), "\[IndentingNewLine]", 
    \(Print["\<For thin-walls neglect all higher order terms of thickness\>"]\
\), "\[IndentingNewLine]", 
    \(\(x\&_\_G = 
        Normal[Series[x\&_\_Gp, {t\_2, 0, 1}]];\)\), "\[IndentingNewLine]", 
    \(Print[\*"\"\<\!\(x\&_\_G\) =\>\"", x\&_\_G, "\< = \>", 
      x\&_\_G /. vars\  // N, "\< mm\>"]\)}], "Input"],

Cell[BoxData[
    \("The first moment of area about the B'C' is:"\)], "Print"],

Cell[BoxData[
    \(\(b\^2\ t\_2\)\/2\)], "Output"],

Cell[BoxData[
    \("The centroid is located at "\)], "Print"],

Cell[BoxData[
    \(\(3\ b\^2\)\/\(6\ b + 5\ h\)\)], "Output"],

Cell[BoxData[
    \("For thin-walls neglect all higher order terms of thickness"\)], "Print"],

Cell[BoxData[
    InterpretationBox[\("\!\(x\&_\_G\) ="\[InvisibleSpace]\(3\ b\^2\)\/\(6\ b \
+ 5\ h\)\[InvisibleSpace]" = "\[InvisibleSpace]54.54545454545455`\
\[InvisibleSpace]" mm"\),
      SequenceForm[ "\!\(x\&_\_G\) =", 
        Times[ 3, 
          Power[ b, 2], 
          Power[ 
            Plus[ 
              Times[ 6, b], 
              Times[ 5, h]], -1]], " = ", 54.545454545454547, " mm"],
      Editable->False]], "Print"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
  "Centroid: ",
  Cell[BoxData[
      FormBox[
        SubscriptBox[
          FormBox[\(y\&_\),
            "TraditionalForm"], "G"], TraditionalForm]]],
  " (let's place the origin at point B)"
}], "Section",
  FontFamily->"Times New Roman"],

Cell[TextData[{
  "Centroid: ",
  Cell[BoxData[
      \(\(\(y\&_\_G\)\(=\)\(\(Sum\ of\ all\ \(A\_j\) 
              y\&_\_j\)\/\(Total\ Area\)\)\(\ \)\)\)]]
}], "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(y\&_\_1 = \ ExpandAll[\(-\((h\/2)\)\)];\)\), "\[IndentingNewLine]", 
    \(\(y\&_\_2 = \ 0\ ;\)\), "\[IndentingNewLine]", 
    \(\(Qy = 
        ExpandAll[\ \((\[Sum]\+\(j = 1\)\%2\( A\_j\) 
                y\&_\_j\ )\)];\)\), "\[IndentingNewLine]", 
    \(Print["\<The first moment of area about the B'A' is:\>"]\), "\
\[IndentingNewLine]", 
    \(\(-ExpandAll[\ \((Qy\ )\)]\)\ \), "\[IndentingNewLine]", 
    \(Print["\<The centroid is located at \>"]\), "\[IndentingNewLine]", 
    \(y\&_\_Gp = \(-FullSimplify[\ \((Qy\/A\_T\ )\)]\)\ \), "\
\[IndentingNewLine]", 
    \(Print["\<For thin-walls neglect all higher order terms of thickness\>"]\
\), "\[IndentingNewLine]", 
    \(\(y\&_\_G = 
        Normal[Series[y\&_\_Gp, {t\_2, 0, 1}]];\)\), "\[IndentingNewLine]", 
    \(Print["\<The centroid is located at \>"]\), "\[IndentingNewLine]", 
    \(Print[\*"\"\<\!\(y\&_\_G\) =\>\"", y\&_\_G, "\< = \>", 
      y\&_\_G /. vars\  // N, "\< mm\>"]\)}], "Input"],

Cell[BoxData[
    \("The first moment of area about the B'A' is:"\)], "Print"],

Cell[BoxData[
    \(\(5\ h\^2\ t\_2\)\/12\)], "Output"],

Cell[BoxData[
    \("The centroid is located at "\)], "Print"],

Cell[BoxData[
    \(\(5\ h\^2\)\/\(2\ \((6\ b + 5\ h)\)\)\)], "Output"],

Cell[BoxData[
    \("For thin-walls neglect all higher order terms of thickness"\)], "Print"],

Cell[BoxData[
    \("The centroid is located at "\)], "Print"],

Cell[BoxData[
    InterpretationBox[\("\!\(y\&_\_G\) ="\[InvisibleSpace]\(5\ h\^2\)\/\(2\ \
\((6\ b + 5\ h)\)\)\[InvisibleSpace]" = "\[InvisibleSpace]45.45454545454545`\
\[InvisibleSpace]" mm"\),
      SequenceForm[ "\!\(y\&_\_G\) =", 
        Times[ 
          Rational[ 5, 2], 
          Power[ h, 2], 
          Power[ 
            Plus[ 
              Times[ 6, b], 
              Times[ 5, h]], -1]], " = ", 45.454545454545453, " mm"],
      Editable->False]], "Print"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(200 - \[InvisibleSpace]45.45454545454545`\)], "Input"],

Cell[BoxData[
    \(154.54545454545456`\)], "Output"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
  "Second Moment of Area ",
  Cell[BoxData[
      \(\[CapitalIota]\_xx\)]],
  " about the neutral axis ( which is located at centroid C)"
}], "Section",
  FontFamily->"Times New Roman"],

Cell["Section 1: Web", "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(dy\_1 = \(-\((h\/2 - \((y\&_\_G\ )\))\)\)\), "\[IndentingNewLine]", 
    \(\(\[CapitalIota]\_xc1 = \(\(t\_1\) h\^3\)\/12;\)\), \
"\[IndentingNewLine]", 
    \(\[CapitalIota]\_xx1p = \ 
      ExpandAll[\ \[CapitalIota]\_xc1\  + \ \(A\_1\) 
            dy\_1\^2\ ]\ \), "\[IndentingNewLine]", 
    \(Print["\<For thin-walls neglect all higher order terms of thickness\>"]\
\), "\[IndentingNewLine]", 
    \(\[CapitalIota]\_xx1 = 
      Normal[Series[\[CapitalIota]\_xx1p, {t\_2, 0, 
            1}]]\), "\[IndentingNewLine]", 
    \(Print["\<The second moment of area about the neutral axis for the web \
is: \>"]\), "\[IndentingNewLine]", 
    \(Print[\*"\"\<\!\(\[CapitalIota]\_xx1\) =\>\"", \[CapitalIota]\_xx1 /. 
          vars\  // N, \*"\"\< \!\(mm\^4\)\>\""]\)}], "Input"],

Cell[BoxData[
    \(\(-\(h\/2\)\) + \(5\ h\^2\)\/\(2\ \((6\ b + 5\ h)\)\)\)], "Output"],

Cell[BoxData[
    \(\(5\ h\^3\ t\_2\)\/18 + \(125\ h\^5\ t\_2\)\/\(6\ \((12\ b + 10\ \
h)\)\^2\) - \(25\ h\^4\ t\_2\)\/\(6\ \((12\ b + 10\ h)\)\)\)], "Output"],

Cell[BoxData[
    \("For thin-walls neglect all higher order terms of thickness"\)], "Print"],

Cell[BoxData[
    \(\((\(5\ h\^3\)\/18 + \(125\ h\^5\)\/\(6\ \((12\ b + 10\ h)\)\^2\) - \
\(25\ h\^4\)\/\(6\ \((12\ b + 10\ h)\)\))\)\ t\_2\)], "Output"],

Cell[BoxData[
    \("The second moment of area about the neutral axis for the web is: "\)], \
"Print"],

Cell[BoxData[
    InterpretationBox[\("\!\(\[CapitalIota]\_xx1\) \
="\[InvisibleSpace]1.2617079889807163`*^7\[InvisibleSpace]" \!\(mm\^4\)"\),
      SequenceForm[ 
      "\!\(\[CapitalIota]\_xx1\) =", 12617079.889807163, " \!\(mm\^4\)"],
      Editable->False]], "Print"]
}, Open  ]],

Cell["Section 2: Flange", "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(dy\_2 = y\&_\_G;\)\), "\[IndentingNewLine]", 
    \(\(\[CapitalIota]\_xc2 = \(\(\ \)\(\(t\_2\^3\) b\)\)\/12;\)\), "\
\[IndentingNewLine]", 
    \(\[CapitalIota]\_xx2p = \ 
      ExpandAll[\ \[CapitalIota]\_xc2\  + \ \(A\_2\) 
            dy\_2\^2]\ \), "\[IndentingNewLine]", 
    \(Print["\<For thin-walls neglect all higher order terms of thickness\>"]\
\), "\[IndentingNewLine]", 
    \(\[CapitalIota]\_xx2 = 
      Normal[Series[\[CapitalIota]\_xx2p, {t\_2, 0, 
            1}]]\), "\[IndentingNewLine]", 
    \(Print["\<The second moment of area about the neutral axis for the \
flange is: \>"]\), "\[IndentingNewLine]", 
    \(Print[\*"\"\<\!\(\[CapitalIota]\_xx2\) =\>\"", \[CapitalIota]\_xx2 /. 
          vars\  // N, \*"\"\< \!\(mm\^4\)\>\""]\)}], "Input"],

Cell[BoxData[
    \(\(25\ b\ h\^4\ t\_2\)\/\(144\ b\^2 + 240\ b\ h + 100\ h\^2\) + \(b\ \
t\_2\%3\)\/12\)], "Output"],

Cell[BoxData[
    \("For thin-walls neglect all higher order terms of thickness"\)], "Print"],

Cell[BoxData[
    \(\(25\ b\ h\^4\ t\_2\)\/\(144\ b\^2 + 240\ b\ h + 100\ h\^2\)\)], \
"Output"],

Cell[BoxData[
    \("The second moment of area about the neutral axis for the flange is: \
"\)], "Print"],

Cell[BoxData[
    InterpretationBox[\("\!\(\[CapitalIota]\_xx2\) \
="\[InvisibleSpace]4.958677685950413`*^6\[InvisibleSpace]" \!\(mm\^4\)"\),
      SequenceForm[ 
      "\!\(\[CapitalIota]\_xx2\) =", 4958677.6859504133, " \!\(mm\^4\)"],
      Editable->False]], "Print"]
}, Open  ]],

Cell["Total", "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\[CapitalIota]\_xxp = 
      ExpandAll[\ \[CapitalIota]\_xx1p + \[CapitalIota]\_xx2p]\), "\
\[IndentingNewLine]", 
    \(\[CapitalIota]\_xx = 
      ExpandAll[\ \[CapitalIota]\_xx1 + \[CapitalIota]\_xx2]\), "\
\[IndentingNewLine]", 
    \(Print[\*"\"\<For thin-walls neglect all higher order terms of \
thickness: \!\(\[CapitalIota]\_xx\) =\>\"", \[CapitalIota]\_xx]\), "\
\[IndentingNewLine]", 
    \(Print["\<The second moment of area about the neutral axis is: \>"]\), "\
\[IndentingNewLine]", 
    \(Print[\*"\"\<\!\(\[CapitalIota]\_xx\) =\>\"", \[CapitalIota]\_xx /. 
          vars\  // N, \*"\"\< \!\(mm\^4\)\>\""]\)}], "Input"],

Cell[BoxData[
    \(\(5\ h\^3\ t\_2\)\/18 - \(25\ h\^4\ t\_2\)\/\(72\ b + 60\ h\) + \(25\ b\
\ h\^4\ t\_2\)\/\(144\ b\^2 + 240\ b\ h + 100\ h\^2\) + \(125\ h\^5\ t\_2\)\/\
\(864\ b\^2 + 1440\ b\ h + 600\ h\^2\) + \(b\ t\_2\%3\)\/12\)], "Output"],

Cell[BoxData[
    \(\(5\ h\^3\ t\_2\)\/18 - \(25\ h\^4\ t\_2\)\/\(72\ b + 60\ h\) + \(25\ b\
\ h\^4\ t\_2\)\/\(144\ b\^2 + 240\ b\ h + 100\ h\^2\) + \(125\ h\^5\ t\_2\)\/\
\(864\ b\^2 + 1440\ b\ h + 600\ h\^2\)\)], "Output"],

Cell[BoxData[
    InterpretationBox[\("For thin-walls neglect all higher order terms of \
thickness: \!\(\[CapitalIota]\_xx\) ="\[InvisibleSpace]\(\(5\ h\^3\ \
t\_2\)\/18 - \(25\ h\^4\ t\_2\)\/\(72\ b + 60\ h\) + \(25\ b\ h\^4\ \
t\_2\)\/\(144\ b\^2 + 240\ b\ h + 100\ h\^2\) + \(125\ h\^5\ t\_2\)\/\(864\ b\
\^2 + 1440\ b\ h + 600\ h\^2\)\)\),
      SequenceForm[ 
      "For thin-walls neglect all higher order terms of thickness: \!\(\
\[CapitalIota]\_xx\) =", 
        Plus[ 
          Times[ 
            Rational[ 5, 18], 
            Power[ h, 3], 
            Subscript[ t, 2]], 
          Times[ -25, 
            Power[ h, 4], 
            Power[ 
              Plus[ 
                Times[ 72, b], 
                Times[ 60, h]], -1], 
            Subscript[ t, 2]], 
          Times[ 25, b, 
            Power[ h, 4], 
            Power[ 
              Plus[ 
                Times[ 144, 
                  Power[ b, 2]], 
                Times[ 240, b, h], 
                Times[ 100, 
                  Power[ h, 2]]], -1], 
            Subscript[ t, 2]], 
          Times[ 125, 
            Power[ h, 5], 
            Power[ 
              Plus[ 
                Times[ 864, 
                  Power[ b, 2]], 
                Times[ 1440, b, h], 
                Times[ 600, 
                  Power[ h, 2]]], -1], 
            Subscript[ t, 2]]]],
      Editable->False]], "Print"],

Cell[BoxData[
    \("The second moment of area about the neutral axis is: "\)], "Print"],

Cell[BoxData[
    InterpretationBox[\("\!\(\[CapitalIota]\_xx\) \
="\[InvisibleSpace]1.7575757575757574`*^7\[InvisibleSpace]" \!\(mm\^4\)"\),
      SequenceForm[ 
      "\!\(\[CapitalIota]\_xx\) =", 17575757.575757574, " \!\(mm\^4\)"],
      Editable->False]], "Print"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
  "Second Moment of Area ",
  Cell[BoxData[
      \(\[CapitalIota]\_yy\)]],
  " about the neutral axis ( which is located at centroid C)"
}], "Section",
  FontFamily->"Courier"],

Cell["Section 1: Web", "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(dx\_1 = \(-\((x\&_\_G\ )\)\);\)\), "\[IndentingNewLine]", 
    \(\(\[CapitalIota]\_yc1 = \(\(t\_1\^3\) h\)\/12;\)\), \
"\[IndentingNewLine]", 
    \(\[CapitalIota]\_yy1p = \ 
      ExpandAll[\ \[CapitalIota]\_yc1\  + \ \(A\_1\) 
            dx\_1\^2]\), "\[IndentingNewLine]", 
    \(Print["\<For thin-walls neglect all higher order terms of thickness\>"]\
\), "\[IndentingNewLine]", 
    \(\[CapitalIota]\_yy1 = 
      Normal[Series[\[CapitalIota]\_yy1p, {t\_2, 0, 
            1}]]\), "\[IndentingNewLine]", 
    \(Print["\<The second moment of area about the neutral axis for the web \
is: \>"]\), "\[IndentingNewLine]", 
    \(Print[\*"\"\<\!\(\[CapitalIota]\_yy1\) =\>\"", \[CapitalIota]\_yy1 /. 
          vars\  // N, \*"\"\< \!\(mm\^4\)\>\""]\)}], "Input"],

Cell[BoxData[
    \(\(15\ b\^4\ h\ t\_2\)\/\(72\ b\^2 + 120\ b\ h + 50\ h\^2\) + \(125\ h\ \
t\_2\%3\)\/2592\)], "Output"],

Cell[BoxData[
    \("For thin-walls neglect all higher order terms of thickness"\)], "Print"],

Cell[BoxData[
    \(\(15\ b\^4\ h\ t\_2\)\/\(72\ b\^2 + 120\ b\ h + 50\ h\^2\)\)], "Output"],

Cell[BoxData[
    \("The second moment of area about the neutral axis for the web is: "\)], \
"Print"],

Cell[BoxData[
    InterpretationBox[\("\!\(\[CapitalIota]\_yy1\) \
="\[InvisibleSpace]5.950413223140496`*^6\[InvisibleSpace]" \!\(mm\^4\)"\),
      SequenceForm[ 
      "\!\(\[CapitalIota]\_yy1\) =", 5950413.2231404958, " \!\(mm\^4\)"],
      Editable->False]], "Print"]
}, Open  ]],

Cell["Section 2: Flange", "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(dx\_2 = b\/2 - x\&_\_G;\)\), "\[IndentingNewLine]", 
    \(\(\[CapitalIota]\_yc2 = \(\(\ \)\(\(t\_2\) b\^3\)\)\/12;\)\), "\
\[IndentingNewLine]", 
    \(\[CapitalIota]\_yy2p = \ 
      ExpandAll[\ \[CapitalIota]\_yc2\  + \ \(A\_2\) 
            dx\_2\^2]\), "\[IndentingNewLine]", 
    \(Print["\<For thin-walls neglect all higher order terms of thickness\>"]\
\), "\[IndentingNewLine]", 
    \(\[CapitalIota]\_yy2 = 
      Normal[Series[\[CapitalIota]\_yy2p, {t\_2, 0, 
            1}]]\), "\[IndentingNewLine]", 
    \(Print["\<The second moment of area about the neutral axis for the \
flange is: \>"]\), "\[IndentingNewLine]", 
    \(Print[\*"\"\<\!\(\[CapitalIota]\_yy2\) =\>\"", \[CapitalIota]\_yy2 /. 
          vars\  // N, \*"\"\< \!\(mm\^4\)\>\""]\)}], "Input"],

Cell[BoxData[
    \(\(b\^3\ t\_2\)\/3 + \(9\ b\^5\ t\_2\)\/\((6\ b + 5\ h)\)\^2 - \(3\ b\^4\
\ t\_2\)\/\(6\ b + 5\ h\)\)], "Output"],

Cell[BoxData[
    \("For thin-walls neglect all higher order terms of thickness"\)], "Print"],

Cell[BoxData[
    \(\((b\^3\/3 + \(9\ b\^5\)\/\((6\ b + 5\ h)\)\^2 - \(3\ b\^4\)\/\(6\ b + \
5\ h\))\)\ t\_2\)], "Output"],

Cell[BoxData[
    \("The second moment of area about the neutral axis for the flange is: \
"\)], "Print"],

Cell[BoxData[
    InterpretationBox[\("\!\(\[CapitalIota]\_yy2\) \
="\[InvisibleSpace]1.2958677685950413`*^7\[InvisibleSpace]" \!\(mm\^4\)"\),
      SequenceForm[ 
      "\!\(\[CapitalIota]\_yy2\) =", 12958677.685950413, " \!\(mm\^4\)"],
      Editable->False]], "Print"]
}, Open  ]],

Cell["Total", "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\[CapitalIota]\_yyp = 
      ExpandAll[\ \[CapitalIota]\_yy1p + \[CapitalIota]\_yy2p]\), "\
\[IndentingNewLine]", 
    \(\[CapitalIota]\_yy = 
      ExpandAll[\ \[CapitalIota]\_yy1 + \[CapitalIota]\_yy2]\), "\
\[IndentingNewLine]", 
    \(Print[\*"\"\<For thin-walls neglect all higher order terms of \
thickness: \!\(\[CapitalIota]\_yy\) =\>\"", \[CapitalIota]\_yy]\), "\
\[IndentingNewLine]", 
    \(Print["\<The second moment of area about the neutral axis is: \>"]\), "\
\[IndentingNewLine]", 
    \(Print[\*"\"\<\!\(\[CapitalIota]\_yy\) =\>\"", \[CapitalIota]\_yy /. 
          vars\  // N, \*"\"\< \!\(mm\^4\)\>\""]\)}], "Input"],

Cell[BoxData[
    \(\(b\^3\ t\_2\)\/3 - \(3\ b\^4\ t\_2\)\/\(6\ b + 5\ h\) + \(9\ b\^5\ \
t\_2\)\/\(36\ b\^2 + 60\ b\ h + 25\ h\^2\) + \(15\ b\^4\ h\ t\_2\)\/\(72\ \
b\^2 + 120\ b\ h + 50\ h\^2\) + \(125\ h\ t\_2\%3\)\/2592\)], "Output"],

Cell[BoxData[
    \(\(b\^3\ t\_2\)\/3 - \(3\ b\^4\ t\_2\)\/\(6\ b + 5\ h\) + \(9\ b\^5\ \
t\_2\)\/\(36\ b\^2 + 60\ b\ h + 25\ h\^2\) + \(15\ b\^4\ h\ t\_2\)\/\(72\ \
b\^2 + 120\ b\ h + 50\ h\^2\)\)], "Output"],

Cell[BoxData[
    InterpretationBox[\("For thin-walls neglect all higher order terms of \
thickness: \!\(\[CapitalIota]\_yy\) ="\[InvisibleSpace]\(\(b\^3\ t\_2\)\/3 - \
\(3\ b\^4\ t\_2\)\/\(6\ b + 5\ h\) + \(9\ b\^5\ t\_2\)\/\(36\ b\^2 + 60\ b\ h \
+ 25\ h\^2\) + \(15\ b\^4\ h\ t\_2\)\/\(72\ b\^2 + 120\ b\ h + 50\ h\^2\)\)\),
      
      SequenceForm[ 
      "For thin-walls neglect all higher order terms of thickness: \!\(\
\[CapitalIota]\_yy\) =", 
        Plus[ 
          Times[ 
            Rational[ 1, 3], 
            Power[ b, 3], 
            Subscript[ t, 2]], 
          Times[ -3, 
            Power[ b, 4], 
            Power[ 
              Plus[ 
                Times[ 6, b], 
                Times[ 5, h]], -1], 
            Subscript[ t, 2]], 
          Times[ 9, 
            Power[ b, 5], 
            Power[ 
              Plus[ 
                Times[ 36, 
                  Power[ b, 2]], 
                Times[ 60, b, h], 
                Times[ 25, 
                  Power[ h, 2]]], -1], 
            Subscript[ t, 2]], 
          Times[ 15, 
            Power[ b, 4], h, 
            Power[ 
              Plus[ 
                Times[ 72, 
                  Power[ b, 2]], 
                Times[ 120, b, h], 
                Times[ 50, 
                  Power[ h, 2]]], -1], 
            Subscript[ t, 2]]]],
      Editable->False]], "Print"],

Cell[BoxData[
    \("The second moment of area about the neutral axis is: "\)], "Print"],

Cell[BoxData[
    InterpretationBox[\("\!\(\[CapitalIota]\_yy\) \
="\[InvisibleSpace]1.890909090909091`*^7\[InvisibleSpace]" \!\(mm\^4\)"\),
      SequenceForm[ 
      "\!\(\[CapitalIota]\_yy\) =", 18909090.90909091, " \!\(mm\^4\)"],
      Editable->False]], "Print"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
  "Second Moment of Area ",
  Cell[BoxData[
      \(\[CapitalIota]\_xy\)]],
  " about the neutral axis ( which is located at centroid C)"
}], "Section",
  FontFamily->"Courier"],

Cell["Section 1: Web", "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(\[CapitalIota]\_xcyc1 = 0;\)\), "\[IndentingNewLine]", 
    \(\[CapitalIota]\_xy1p = \ 
      ExpandAll[\ \[CapitalIota]\_xcyc1\  + \ \(A\_1\) \((dx\_1)\) \
\((dy\_1)\)]\), "\[IndentingNewLine]", 
    \(Print["\<For thin-walls neglect all higher order terms of thickness\>"]\
\), "\[IndentingNewLine]", 
    \(\[CapitalIota]\_xy1 = 
      Normal[Series[\[CapitalIota]\_xy1p, {t\_2, 0, 
            1}]]\), "\[IndentingNewLine]", 
    \(Print["\<The second moment of area about the neutral axis for the web \
is: \>"]\), "\[IndentingNewLine]", 
    \(Print[\*"\"\<\!\(\[CapitalIota]\_xy1\) =\>\"", \[CapitalIota]\_xy1 /. 
          vars\  // N, \*"\"\< \!\(mm\^4\)\>\""]\)}], "Input"],

Cell[BoxData[
    \(\(-\(\(25\ b\^2\ h\^3\ t\_2\)\/\((12\ b + 10\ h)\)\^2\)\) + \(5\ b\^2\ \
h\^2\ t\_2\)\/\(2\ \((12\ b + 10\ h)\)\)\)], "Output"],

Cell[BoxData[
    \("For thin-walls neglect all higher order terms of thickness"\)], "Print"],

Cell[BoxData[
    \(\((\(-\(\(25\ b\^2\ h\^3\)\/\((12\ b + 10\ h)\)\^2\)\) + \(5\ b\^2\ \
h\^2\)\/\(2\ \((12\ b + 10\ h)\)\))\)\ t\_2\)], "Output"],

Cell[BoxData[
    \("The second moment of area about the neutral axis for the web is: "\)], \
"Print"],

Cell[BoxData[
    InterpretationBox[\("\!\(\[CapitalIota]\_xy1\) \
="\[InvisibleSpace]5.950413223140496`*^6\[InvisibleSpace]" \!\(mm\^4\)"\),
      SequenceForm[ 
      "\!\(\[CapitalIota]\_xy1\) =", 5950413.2231404958, " \!\(mm\^4\)"],
      Editable->False]], "Print"]
}, Open  ]],

Cell["Section 2: Flange", "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(\[CapitalIota]\_xcyc2 = 0;\)\), "\[IndentingNewLine]", 
    \(\[CapitalIota]\_xy2p = \ 
      ExpandAll[\ \[CapitalIota]\_xcyc2\  + \ \(A\_2\) 
            dx\_2\ dy\_2]\), "\[IndentingNewLine]", 
    \(Print["\<For thin-walls neglect all higher order terms of thickness\>"]\
\), "\[IndentingNewLine]", 
    \(\[CapitalIota]\_xy2 = 
      Normal[Series[\[CapitalIota]\_xy2p, {t\_2, 0, 
            1}]]\), "\[IndentingNewLine]", 
    \(Print["\<The second moment of area about the neutral axis for the \
flange is: \>"]\), "\[IndentingNewLine]", 
    \(Print[\*"\"\<\!\(\[CapitalIota]\_xy2\) =\>\"", \[CapitalIota]\_xy2 /. 
          vars\  // N, \*"\"\< \!\(mm\^4\)\>\""]\)}], "Input"],

Cell[BoxData[
    \(\(5\ b\^2\ h\^2\ t\_2\)\/\(2\ \((12\ b + 10\ h)\)\) - \(15\ b\^3\ h\^2\ \
t\_2\)\/\(\((6\ b + 5\ h)\)\ \((12\ b + 10\ h)\)\)\)], "Output"],

Cell[BoxData[
    \("For thin-walls neglect all higher order terms of thickness"\)], "Print"],

Cell[BoxData[
    \(\((\(5\ b\^2\ h\^2\)\/\(2\ \((12\ b + 10\ h)\)\) - \(15\ b\^3\ h\^2\)\/\
\(\((6\ b + 5\ h)\)\ \((12\ b + 10\ h)\)\))\)\ t\_2\)], "Output"],

Cell[BoxData[
    \("The second moment of area about the neutral axis for the flange is: \
"\)], "Print"],

Cell[BoxData[
    InterpretationBox[\("\!\(\[CapitalIota]\_xy2\) \
="\[InvisibleSpace]4.958677685950413`*^6\[InvisibleSpace]" \!\(mm\^4\)"\),
      SequenceForm[ 
      "\!\(\[CapitalIota]\_xy2\) =", 4958677.6859504133, " \!\(mm\^4\)"],
      Editable->False]], "Print"]
}, Open  ]],

Cell["Total", "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\[CapitalIota]\_xy = 
      ExpandAll[\ \[CapitalIota]\_xy1 + \[CapitalIota]\_xy2]\), "\
\[IndentingNewLine]", 
    \(Print[\*"\"\<For thin-walls neglect all higher order terms of \
thickness: \!\(\[CapitalIota]\_xy\) =\>\"", \[CapitalIota]\_xy]\), "\
\[IndentingNewLine]", 
    \(Print["\<The second moment of area about the neutral axis is: \>"]\), "\
\[IndentingNewLine]", 
    \(Print[\*"\"\<\!\(\[CapitalIota]\_xy\) =\>\"", \[CapitalIota]\_xy /. 
          vars\  // N, \*"\"\< \!\(mm\^4\)\>\""]\)}], "Input"],

Cell[BoxData[
    \(\(10\ b\^2\ h\^2\ t\_2\)\/\(24\ b + 20\ h\) - \(15\ b\^3\ h\^2\ \
t\_2\)\/\(72\ b\^2 + 120\ b\ h + 50\ h\^2\) - \(25\ b\^2\ h\^3\ t\_2\)\/\(144\
\ b\^2 + 240\ b\ h + 100\ h\^2\)\)], "Output"],

Cell[BoxData[
    InterpretationBox[\("For thin-walls neglect all higher order terms of \
thickness: \!\(\[CapitalIota]\_xy\) ="\[InvisibleSpace]\(\(10\ b\^2\ h\^2\ \
t\_2\)\/\(24\ b + 20\ h\) - \(15\ b\^3\ h\^2\ t\_2\)\/\(72\ b\^2 + 120\ b\ h \
+ 50\ h\^2\) - \(25\ b\^2\ h\^3\ t\_2\)\/\(144\ b\^2 + 240\ b\ h + 100\ \
h\^2\)\)\),
      SequenceForm[ 
      "For thin-walls neglect all higher order terms of thickness: \!\(\
\[CapitalIota]\_xy\) =", 
        Plus[ 
          Times[ 10, 
            Power[ b, 2], 
            Power[ h, 2], 
            Power[ 
              Plus[ 
                Times[ 24, b], 
                Times[ 20, h]], -1], 
            Subscript[ t, 2]], 
          Times[ -15, 
            Power[ b, 3], 
            Power[ h, 2], 
            Power[ 
              Plus[ 
                Times[ 72, 
                  Power[ b, 2]], 
                Times[ 120, b, h], 
                Times[ 50, 
                  Power[ h, 2]]], -1], 
            Subscript[ t, 2]], 
          Times[ -25, 
            Power[ b, 2], 
            Power[ h, 3], 
            Power[ 
              Plus[ 
                Times[ 144, 
                  Power[ b, 2]], 
                Times[ 240, b, h], 
                Times[ 100, 
                  Power[ h, 2]]], -1], 
            Subscript[ t, 2]]]],
      Editable->False]], "Print"],

Cell[BoxData[
    \("The second moment of area about the neutral axis is: "\)], "Print"],

Cell[BoxData[
    InterpretationBox[\("\!\(\[CapitalIota]\_xy\) \
="\[InvisibleSpace]1.0909090909090908`*^7\[InvisibleSpace]" \!\(mm\^4\)"\),
      SequenceForm[ 
      "\!\(\[CapitalIota]\_xy\) =", 10909090.909090908, " \!\(mm\^4\)"],
      Editable->False]], "Print"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Bending stress ", "Section",
  FontFamily->"Courier"],

Cell["Moment resultants", "Text",
  FontFamily->"Times New Roman"],

Cell[BoxData[{
    \(\(Mx[z_] = \(-P\)\ \((L - z)\);\)\), "\[IndentingNewLine]", 
    \(\(My[z_] = \(-\(P\/2\)\) \((L - z)\);\)\)}], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["PART A.  Deflection components", "Section",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell["Horizontal deflection", "Subsection",
  FontFamily->"Times New Roman"],

Cell[" (equation 9.17; page 285, Megson's text)", "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(upp = 
      ExpandAll[
          FullSimplify[\ \(Mx[z]\ \[CapitalIota]\_xy - My[z]\ \
\[CapitalIota]\_xx\)\/\(\[CapitalEpsilon] \((\(\[CapitalIota]\_xx\) \
\[CapitalIota]\_yy - \[CapitalIota]\_xy\^2)\)\)]]\ \  /. 
        vars\), "\[IndentingNewLine]", 
    \(up = 
      ExpandAll[FullSimplify[\ \[Integral]upp \[DifferentialD]z]]\  + 
        A1\), "\[IndentingNewLine]", 
    \(u = 
      ExpandAll[FullSimplify[\ \[Integral]up \[DifferentialD]z]] + 
        A2\), "\[IndentingNewLine]", 
    \(\(BC1 = \((up\  /. \ z \[Rule] \ 0)\) \[Equal] 
          0;\)\), "\[IndentingNewLine]", 
    \(\(BC2 = \((u\  /. \ z \[Rule] \ 0)\)\  \[Equal] 
          0;\)\), "\[IndentingNewLine]", 
    \(Cons1 = Flatten[Solve[{BC1, BC2}, {A1, A2}]]\), "\[IndentingNewLine]", 
    \(Print["\<u(z)=\>", \ 
      ExpandAll[u\  /. vars\ \  // N]]\), "\[IndentingNewLine]", 
    \(u\  = \ ExpandAll[\(u\  /. vars\)\ \  /. Cons1\  // N]\)}], "Input"],

Cell[BoxData[
    \(\(-\(\(21\ P\)\/\(704000\ \[CapitalEpsilon]\)\)\) + \(7\ P\ \
z\)\/\(704000000\ \[CapitalEpsilon]\)\)], "Output"],

Cell[BoxData[
    \(A1 - \(21\ P\ z\)\/\(704000\ \[CapitalEpsilon]\) + \(7\ P\ \
z\^2\)\/\(1408000000\ \[CapitalEpsilon]\)\)], "Output"],

Cell[BoxData[
    \(A2 + 
      A1\ z - \(21\ P\ z\^2\)\/\(1408000\ \[CapitalEpsilon]\) + \(7\ P\ \
z\^3\)\/\(4224000000\ \[CapitalEpsilon]\)\)], "Output"],

Cell[BoxData[
    \({A1 \[Rule] 0, A2 \[Rule] 0}\)], "Output"],

Cell[BoxData[
    InterpretationBox[\("u(z)="\[InvisibleSpace]\(A2 + 
          A1\ z - \(0.000014914772727272727`\ P\ z\^2\)\/\[CapitalEpsilon] + \
\(1.6571969696969698`*^-9\ P\ z\^3\)\/\[CapitalEpsilon]\)\),
      SequenceForm[ "u(z)=", 
        Plus[ A2, 
          Times[ A1, z], 
          Times[ -.000014914772727272727, P, 
            Power[ z, 2], 
            Power[ \[CapitalEpsilon], -1]], 
          Times[ .16571969696969698*^-8, P, 
            Power[ z, 3], 
            Power[ \[CapitalEpsilon], -1]]]],
      Editable->False]], "Print"],

Cell[BoxData[
    \(\(-\(\(0.000014914772727272727`\ P\ z\^2\)\/\[CapitalEpsilon]\)\) + \
\(1.6571969696969698`*^-9\ P\ z\^3\)\/\[CapitalEpsilon]\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(TeXForm[Expand[v]]\)], "Input"],

Cell["\<\
\\frac{0.0000946023\\,P\\,z^2}{\\CapitalEpsilon } - 
  \\frac{1.05114\\,{10}^{-8}\\,P\\,z^3}{\\CapitalEpsilon }\
\>", "Output"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Vertical deflection", "Subsection",
  FontFamily->"Times New Roman"],

Cell[" (equation 9.17; page 285, Megson's text)", "Text",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(vpp = \ 
      FullSimplify[\ \(My[z]\ \[CapitalIota]\_xy - Mx[z]\ \
\[CapitalIota]\_yy\)\/\(\[CapitalEpsilon] \((\(\[CapitalIota]\_xx\) \
\[CapitalIota]\_yy - \[CapitalIota]\_xy\^2)\)\)]\  /. 
        vars\), "\[IndentingNewLine]", 
    \(vp = 
      ExpandAll[FullSimplify[\ \[Integral]vpp \[DifferentialD]z]]\  + 
        B1\), "\[IndentingNewLine]", 
    \(v = 
      ExpandAll[FullSimplify[\ \[Integral]vp \[DifferentialD]z]] + 
        B2\), "\[IndentingNewLine]", 
    \(\(BC3 = \((vp\  /. \ z \[Rule] \ 0)\) \[Equal] 
          0;\)\), "\[IndentingNewLine]", 
    \(\(BC4 = \((v\  /. \ z \[Rule] \ 0)\)\  \[Equal] 
          0;\)\), "\[IndentingNewLine]", 
    \(Cons2 = Flatten[Solve[{BC3, BC4}, {B1, B2}]]\), "\[IndentingNewLine]", 
    \(Print["\<v(z)=\>", \ 
      ExpandAll[v\  /. vars\ \  // N]]\), "\[IndentingNewLine]", 
    \(v\  = \ ExpandAll[\(v\  /. vars\)\ \  /. Cons2\  // N]\)}], "Input"],

Cell[BoxData[
    \(\(-\(\(111\ P\ \((\(-3000\) + 
                z)\)\)\/\(1760000000\ \[CapitalEpsilon]\)\)\)\)], "Output"],

Cell[BoxData[
    \(B1 + \(333\ P\ z\)\/\(1760000\ \[CapitalEpsilon]\) - \(111\ P\ z\^2\)\/\
\(3520000000\ \[CapitalEpsilon]\)\)], "Output"],

Cell[BoxData[
    \(B2 + 
      B1\ z + \(333\ P\ z\^2\)\/\(3520000\ \[CapitalEpsilon]\) - \(37\ P\ \
z\^3\)\/\(3520000000\ \[CapitalEpsilon]\)\)], "Output"],

Cell[BoxData[
    \({B1 \[Rule] 0, B2 \[Rule] 0}\)], "Output"],

Cell[BoxData[
    InterpretationBox[\("v(z)="\[InvisibleSpace]\(B2 + 
          B1\ z + \(0.00009460227272727273`\ P\ z\^2\)\/\[CapitalEpsilon] - \
\(1.0511363636363636`*^-8\ P\ z\^3\)\/\[CapitalEpsilon]\)\),
      SequenceForm[ "v(z)=", 
        Plus[ B2, 
          Times[ B1, z], 
          Times[ .000094602272727272731, P, 
            Power[ z, 2], 
            Power[ \[CapitalEpsilon], -1]], 
          Times[ -.10511363636363636*^-7, P, 
            Power[ z, 3], 
            Power[ \[CapitalEpsilon], -1]]]],
      Editable->False]], "Print"],

Cell[BoxData[
    \(\(0.00009460227272727273`\ P\ z\^2\)\/\[CapitalEpsilon] - \
\(1.0511363636363636`*^-8\ P\ z\^3\)\/\[CapitalEpsilon]\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(v\)], "Input"],

Cell[BoxData[
    \(\(0.00009460227272727273`\ P\ z\^2\)\/\[CapitalEpsilon] - \
\(1.0511363636363636`*^-8\ P\ z\^3\)\/\[CapitalEpsilon]\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(u\)], "Input"],

Cell[BoxData[
    \(\(-\(\(0.000014914772727272727`\ P\ z\^2\)\/\[CapitalEpsilon]\)\) + \
\(1.6571969696969698`*^-9\ P\ z\^3\)\/\[CapitalEpsilon]\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(u/v\)], "Input"],

Cell[BoxData[
    \(\(\(-\(\(0.000014914772727272727`\ P\ z\^2\)\/\[CapitalEpsilon]\)\) + \
\(1.6571969696969698`*^-9\ P\ \
z\^3\)\/\[CapitalEpsilon]\)\/\(\(0.00009460227272727273`\ P\ z\^2\)\/\
\[CapitalEpsilon] - \(1.0511363636363636`*^-8\ P\ z\^3\)\/\[CapitalEpsilon]\)\
\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(FullSimplify[u/v]\)], "Input"],

Cell[BoxData[
    \(\(-0.15765765765765766`\)\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(v\  /. \ z \[Rule] \ L\)], "Input"],

Cell[BoxData[
    \(\(567.6136363636365`\ P\)\/\[CapitalEpsilon]\)], "Output"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Total deflection", "Subsection",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(w[z_] = \ 
        FullSimplify[\ \@ExpandAll[v\^2 + u\^2]];\)\), "\[IndentingNewLine]", \

    \(Print["\<w(z)=\>", \ 
      FullSimplify[w[z] /. vars\ \ ] // 
        N, \*"\"\< N/\!\(m\^2\)\>\""]\)}], "Input"],

Cell[BoxData[
    InterpretationBox[\("w(z)="\[InvisibleSpace]\(1.0641196703953959`*^-8\ \@\
\(\(P\^2\ \((\(-9000.00013356889`\) + z)\)\ \((\(-8999.999866431113`\) + z)\)\
\ z\^4\)\/\[CapitalEpsilon]\^2\)\)\[InvisibleSpace]" N/\!\(m\^2\)"\),
      SequenceForm[ "w(z)=", 
        Times[ .10641196703953959*^-7, 
          Power[ 
            Times[ 
              Power[ P, 2], 
              Plus[ -9000.0001335688903, z], 
              Plus[ -8999.9998664311133, z], 
              Power[ z, 4], 
              Power[ \[CapitalEpsilon], -2]], 
            Rational[ 1, 2]]], " N/\!\(m\^2\)"],
      Editable->False]], "Print"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(Angle\  = 
      ArcTan[FullSimplify[ExpandAll[v\/u]\ \  /. vars]] 
        180. \/\[Pi]\)], "Input"],

Cell[BoxData[
    \(\(-81.04062709282768`\)\)], "Output"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["AT THE TIP:  z=L ", "Subsection",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(wx\_tip = 
      ExpandAll[\(u\  /. \ z \[Rule] \ L\)\  /. \ 
          vars]\), "\[IndentingNewLine]", 
    \(wy\_tip = ExpandAll[\(v\  /. \ z \[Rule] \ L\)\  /. \ vars]\)}], "Input"],

Cell[BoxData[
    \(\(-\(\(89.48863636363637`\ P\)\/\[CapitalEpsilon]\)\)\)], "Output"],

Cell[BoxData[
    \(\(567.6136363636365`\ P\)\/\[CapitalEpsilon]\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(wL = 
      FullSimplify[\((\@ExpandAll[\ \((w[z]\ )\)\^2]\ \  /. 
                z \[Rule] \ L\ \ )\) /. vars // N\ ]\)], "Input"],

Cell[BoxData[
    \(574.6246220135139`\ \@\(P\^2\/\[CapitalEpsilon]\^2\)\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(AngleL\  = 
      ArcTan[FullSimplify[ExpandAll[v\/u /. \ z \[Rule] \ L]\ \  /. vars]] 
        180. \/\[Pi]\)], "Input"],

Cell[BoxData[
    \(\(-81.04062709282768`\)\)], "Output"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["AT THE MID-POINT:  z=L/2 ", "Subsection",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(wx\_mid = 
      ExpandAll[\(u\  /. \ z \[Rule] \ L/2\)\  /. \ 
          vars]\), "\[IndentingNewLine]", 
    \(wy\_mid = 
      ExpandAll[\(v\  /. \ z \[Rule] \ L/2\)\  /. \ vars]\)}], "Input"],

Cell[BoxData[
    \(\(-\(\(27.965198863636367`\ P\)\/\[CapitalEpsilon]\)\)\)], "Output"],

Cell[BoxData[
    \(\(177.37926136363637`\ P\)\/\[CapitalEpsilon]\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(wmid = 
      FullSimplify[\((\@ExpandAll[\ \((w[z]\ )\)\^2]\ \  /. \ {z \[Rule] \ 
                    L/2}\ )\) /. vars // N\ ]\)], "Input"],

Cell[BoxData[
    \(179.5701943792231`\ \@\(P\^2\/\[CapitalEpsilon]\^2\)\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(Anglemid\  = 
      ArcTan[FullSimplify[ExpandAll[v\/u]\ \  /. vars]] 
        180. \/\[Pi]\)], "Input"],

Cell[BoxData[
    \(\(-81.04062709282768`\)\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Neutral Axis", "Subsubsection",
  FontFamily->"Times New Roman"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(\(\[Sigma]\_zz = \ 
        ExpandAll[\(My[z]\ \[CapitalIota]\_xx - Mx[z]\ \[CapitalIota]\_xy\)\/\
\(\(\[CapitalIota]\_xx\) \[CapitalIota]\_yy - \[CapitalIota]\_xy\^2\)] x + \ 
          ExpandAll[\(Mx[z]\ \[CapitalIota]\_yy - My[z]\ \[CapitalIota]\_xy\)\
\/\(\(\[CapitalIota]\_xx\) \[CapitalIota]\_yy - \[CapitalIota]\_xy\^2\)] 
            y;\)\), "\[IndentingNewLine]", 
    \(Print[\*"\"\<\!\(\[Sigma]\_zz\)=\>\"", \ \[Sigma]\_zz\  /. vars\ \  // 
        N]\), "\[IndentingNewLine]", 
    \(Print[\*"\"\<At the tip:  \!\(\[Sigma]\_zz0\)=\>\""]\), "\
\[IndentingNewLine]", 
    \(\[Sigma]\_zz0 = \ \(\[Sigma]\_zz\  /. vars\)\  /. z \[Rule] \ 0\  // 
        N\)}], "Input"],

Cell[BoxData[
    InterpretationBox[\("\!\(\[Sigma]\_zz\)="\[InvisibleSpace]\(x\ \
\((0.000029829545454545455`\ P - 9.943181818181819`*^-9\ P\ z)\) + 
          y\ \((\(-0.00018920454545454546`\)\ P + 
                6.306818181818181`*^-8\ P\ z)\)\)\),
      SequenceForm[ "\!\(\[Sigma]\_zz\)=", 
        Plus[ 
          Times[ x, 
            Plus[ 
              Times[ .000029829545454545455, P], 
              Times[ -.99431818181818188*^-8, P, z]]], 
          Times[ y, 
            Plus[ 
              Times[ -.00018920454545454546, P], 
              Times[ .63068181818181815*^-7, P, z]]]]],
      Editable->False]], "Print"],

Cell[BoxData[
    \("At the tip:  \!\(\[Sigma]\_zz0\)="\)], "Print"],

Cell[BoxData[
    \(0.000029829545454545455`\ P\ x - 
      0.00018920454545454546`\ P\ y\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[{
    \(Ans = 
      ExpandAll[Flatten[Solve[\[Sigma]\_zz0 \[Equal] 0\  /. vars\ , \ y]]] // 
        N\), "\[IndentingNewLine]", 
    \(tan\[Alpha] = \(Ans[\([1]\)]\)[\([2]\)]/x\)}], "Input"],

Cell[BoxData[
    \({y \[Rule] 0.15765765765765766`\ x}\)], "Output"],

Cell[BoxData[
    \(0.15765765765765766`\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(\[Alpha] = ArcTan[tan\[Alpha]] 180\/\[Pi]\)], "Input"],

Cell[BoxData[
    \(8.95937290717233`\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
    \(Anglemid - \[Alpha]\)], "Input"],

Cell[BoxData[
    \(\(-90.`\)\)], "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
FrontEndVersion->"4.1 for Microsoft Windows",
ScreenRectangle->{{0, 1400}, {0, 1005}},
WindowToolbars->{"RulerBar", "EditBar"},
WindowSize->{1133, 729},
WindowMargins->{{0, Automatic}, {Automatic, 54}},
PrintingCopies->1,
PrintingPageRange->{Automatic, Automatic},
PageHeaders->{{Inherited, Inherited, Inherited}, {Cell[ 
        TextData[ {"Soln to Homework 7"}], "Header"], Inherited, Inherited}},
ShowCellLabel->False,
PrivateFontOptions->{"FontType"->"Outline"},
Magnification->1.5,
StyleDefinitions -> "Classroom.nb"
]

(*******************************************************************
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[1705, 50, 167, 3, 101, "Title"],
Cell[1875, 55, 88, 1, 71, "Input"],

Cell[CellGroupData[{
Cell[1988, 60, 88, 1, 90, "Section"],
Cell[2079, 63, 362, 7, 197, "Input"]
}, Open  ]],

Cell[CellGroupData[{
Cell[2478, 75, 90, 1, 90, "Section"],
Cell[2571, 78, 64, 1, 39, "Text"],

Cell[CellGroupData[{
Cell[2660, 83, 284, 5, 122, "Input"],
Cell[2947, 90, 401, 8, 47, "Print"]
}, Open  ]],
Cell[3363, 101, 66, 1, 39, "Text"],

Cell[CellGroupData[{
Cell[3454, 106, 287, 5, 122, "Input"],
Cell[3744, 113, 374, 7, 35, "Print"]
}, Open  ]],
Cell[4133, 123, 59, 1, 39, "Text"],

Cell[CellGroupData[{
Cell[4217, 128, 569, 11, 188, "Input"],
Cell[4789, 141, 61, 1, 83, "Output"],
Cell[4853, 144, 93, 1, 34, "Print"],
Cell[4949, 147, 505, 11, 47, "Print"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{
Cell[5503, 164, 259, 9, 90, "Section"],
Cell[5765, 175, 201, 6, 44, "Text"],

Cell[CellGroupData[{
Cell[5991, 185, 889, 17, 395, "Input"],
Cell[6883, 204, 78, 1, 34, "Print"],
Cell[6964, 207, 51, 1, 87, "Output"],
Cell[7018, 210, 62, 1, 34, "Print"],
Cell[7083, 213, 62, 1, 87, "Output"],
Cell[7148, 216, 93, 1, 34, "Print"],
Cell[7244, 219, 440, 11, 51, "Print"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{
Cell[7733, 236, 259, 9, 90, "Section"],
Cell[7995, 247, 201, 6, 44, "Text"],

Cell[CellGroupData[{
Cell[8221, 257, 983, 18, 420, "Input"],
Cell[9207, 277, 78, 1, 34, "Print"],
Cell[9288, 280, 55, 1, 87, "Output"],
Cell[9346, 283, 62, 1, 34, "Print"],
Cell[9411, 286, 71, 1, 90, "Output"],
Cell[9485, 289, 93, 1, 34, "Print"],
Cell[9581, 292, 62, 1, 34, "Print"],
Cell[9646, 295, 474, 12, 54, "Print"]
}, Open  ]],

Cell[CellGroupData[{
Cell[10157, 312, 74, 1, 71, "Input"],
Cell[10234, 315, 53, 1, 70, "Output"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{
Cell[10336, 322, 201, 6, 90, "Section"],
Cell[10540, 330, 63, 1, 39, "Text"],

Cell[CellGroupData[{
Cell[10628, 335, 799, 15, 273, "Input"],
Cell[11430, 352, 87, 1, 90, "Output"],
Cell[11520, 355, 159, 2, 90, "Output"],
Cell[11682, 359, 93, 1, 34, "Print"],
Cell[11778, 362, 153, 2, 90, "Output"],
Cell[11934, 366, 102, 2, 34, "Print"],
Cell[12039, 370, 271, 5, 35, "Print"]
}, Open  ]],
Cell[12325, 378, 66, 1, 39, "Text"],

Cell[CellGroupData[{
Cell[12416, 383, 789, 15, 252, "Input"],
Cell[13208, 400, 117, 2, 87, "Output"],
Cell[13328, 404, 93, 1, 34, "Print"],
Cell[13424, 407, 96, 2, 87, "Output"],
Cell[13523, 411, 105, 2, 34, "Print"],
Cell[13631, 415, 270, 5, 35, "Print"]
}, Open  ]],
Cell[13916, 423, 54, 1, 39, "Text"],

Cell[CellGroupData[{
Cell[13995, 428, 656, 13, 173, "Input"],
Cell[14654, 443, 245, 3, 87, "Output"],
Cell[14902, 448, 224, 3, 87, "Output"],
Cell[15129, 453, 1400, 41, 77, "Print"],
Cell[16532, 496, 88, 1, 34, "Print"],
Cell[16623, 499, 269, 5, 35, "Print"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{
Cell[16941, 510, 193, 6, 128, "Section"],
Cell[17137, 518, 63, 1, 39, "Text"],

Cell[CellGroupData[{
Cell[17225, 523, 787, 15, 252, "Input"],
Cell[18015, 540, 122, 2, 87, "Output"],
Cell[18140, 544, 93, 1, 34, "Print"],
Cell[18236, 547, 92, 1, 87, "Output"],
Cell[18331, 550, 102, 2, 34, "Print"],
Cell[18436, 554, 270, 5, 35, "Print"]
}, Open  ]],
Cell[18721, 562, 66, 1, 39, "Text"],

Cell[CellGroupData[{
Cell[18812, 567, 794, 15, 274, "Input"],
Cell[19609, 584, 132, 2, 90, "Output"],
Cell[19744, 588, 93, 1, 34, "Print"],
Cell[19840, 591, 122, 2, 90, "Output"],
Cell[19965, 595, 105, 2, 34, "Print"],
Cell[20073, 599, 271, 5, 35, "Print"]
}, Open  ]],
Cell[20359, 607, 54, 1, 39, "Text"],

Cell[CellGroupData[{
Cell[20438, 612, 656, 13, 173, "Input"],
Cell[21097, 627, 237, 3, 87, "Output"],
Cell[21337, 632, 209, 3, 87, "Output"],
Cell[21549, 637, 1378, 41, 77, "Print"],
Cell[22930, 680, 88, 1, 34, "Print"],
Cell[23021, 683, 267, 5, 35, "Print"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{
Cell[23337, 694, 193, 6, 128, "Section"],
Cell[23533, 702, 63, 1, 39, "Text"],

Cell[CellGroupData[{
Cell[23621, 707, 706, 13, 199, "Input"],
Cell[24330, 722, 147, 2, 90, "Output"],
Cell[24480, 726, 93, 1, 34, "Print"],
Cell[24576, 729, 147, 2, 90, "Output"],
Cell[24726, 733, 102, 2, 34, "Print"],
Cell[24831, 737, 270, 5, 35, "Print"]
}, Open  ]],
Cell[25116, 745, 66, 1, 39, "Text"],

Cell[CellGroupData[{
Cell[25207, 750, 709, 13, 199, "Input"],
Cell[25919, 765, 158, 2, 90, "Output"],
Cell[26080, 769, 93, 1, 34, "Print"],
Cell[26176, 772, 158, 2, 90, "Output"],
Cell[26337, 776, 105, 2, 34, "Print"],
Cell[26445, 780, 270, 5, 35, "Print"]
}, Open  ]],
Cell[26730, 788, 54, 1, 39, "Text"],

Cell[CellGroupData[{
Cell[26809, 793, 535, 10, 148, "Input"],
Cell[27347, 805, 211, 3, 87, "Output"],
Cell[27561, 810, 1354, 40, 77, "Print"],
Cell[28918, 852, 88, 1, 34, "Print"],
Cell[29009, 855, 269, 5, 35, "Print"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{
Cell[29327, 866, 59, 1, 90, "Section"],
Cell[29389, 869, 66, 1, 39, "Text"],
Cell[29458, 872, 141, 2, 117, "Input"]
}, Open  ]],

Cell[CellGroupData[{
Cell[29636, 879, 82, 1, 90, "Section"],

Cell[CellGroupData[{
Cell[29743, 884, 76, 1, 74, "Subsection"],
Cell[29822, 887, 90, 1, 39, "Text"],

Cell[CellGroupData[{
Cell[29937, 892, 954, 20, 326, "Input"],
Cell[30894, 914, 133, 2, 83, "Output"],
Cell[31030, 918, 136, 2, 87, "Output"],
Cell[31169, 922, 155, 3, 87, "Output"],
Cell[31327, 927, 62, 1, 70, "Output"],
Cell[31392, 930, 554, 13, 51, "Print"],
Cell[31949, 945, 159, 2, 87, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[32145, 952, 51, 1, 71, "Input"],
Cell[32199, 955, 137, 3, 89, "Output"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{
Cell[32385, 964, 74, 1, 74, "Subsection"],
Cell[32462, 967, 90, 1, 39, "Text"],

Cell[CellGroupData[{
Cell[32577, 972, 932, 19, 326, "Input"],
Cell[33512, 993, 126, 2, 83, "Output"],
Cell[33641, 997, 140, 2, 87, "Output"],
Cell[33784, 1001, 157, 3, 87, "Output"],
Cell[33944, 1006, 62, 1, 70, "Output"],
Cell[34009, 1009, 553, 13, 51, "Print"],
Cell[34565, 1024, 149, 2, 87, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[34751, 1031, 34, 1, 71, "Input"],
Cell[34788, 1034, 149, 2, 87, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[34974, 1041, 34, 1, 71, "Input"],
Cell[35011, 1044, 159, 2, 87, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[35207, 1051, 36, 1, 71, "Input"],
Cell[35246, 1054, 288, 5, 111, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[35571, 1064, 50, 1, 71, "Input"],
Cell[35624, 1067, 59, 1, 70, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[35720, 1073, 55, 1, 71, "Input"],
Cell[35778, 1076, 78, 1, 83, "Output"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{
Cell[35905, 1083, 71, 1, 74, "Subsection"],

Cell[CellGroupData[{
Cell[36001, 1088, 235, 6, 105, "Input"],
Cell[36239, 1096, 629, 14, 67, "Print"]
}, Open  ]],

Cell[CellGroupData[{
Cell[36905, 1115, 121, 3, 90, "Input"],
Cell[37029, 1120, 57, 1, 70, "Output"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{
Cell[37135, 1127, 72, 1, 74, "Subsection"],

Cell[CellGroupData[{
Cell[37232, 1132, 205, 4, 98, "Input"],
Cell[37440, 1138, 87, 1, 83, "Output"],
Cell[37530, 1141, 78, 1, 83, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[37645, 1147, 152, 3, 78, "Input"],
Cell[37800, 1152, 86, 1, 103, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[37923, 1158, 141, 3, 90, "Input"],
Cell[38067, 1163, 57, 1, 70, "Output"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{
Cell[38173, 1170, 80, 1, 74, "Subsection"],

Cell[CellGroupData[{
Cell[38278, 1175, 216, 5, 96, "Input"],
Cell[38497, 1182, 88, 1, 83, "Output"],
Cell[38588, 1185, 79, 1, 83, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[38704, 1191, 162, 3, 78, "Input"],
Cell[38869, 1196, 86, 1, 103, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[38992, 1202, 124, 3, 90, "Input"],
Cell[39119, 1207, 57, 1, 70, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[39213, 1213, 70, 1, 65, "Subsubsection"],

Cell[CellGroupData[{
Cell[39308, 1218, 698, 12, 176, "Input"],
Cell[40009, 1232, 639, 15, 35, "Print"],
Cell[40651, 1249, 68, 1, 34, "Print"],
Cell[40722, 1252, 103, 2, 70, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[40862, 1259, 205, 4, 96, "Input"],
Cell[41070, 1265, 69, 1, 70, "Output"],
Cell[41142, 1268, 54, 1, 70, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[41233, 1274, 74, 1, 90, "Input"],
Cell[41310, 1277, 51, 1, 70, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[41398, 1283, 52, 1, 71, "Input"],
Cell[41453, 1286, 43, 1, 70, "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)



(*******************************************************************
End of Mathematica Notebook file.
*******************************************************************)