(* Content-type: application/mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 7.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 145, 7] NotebookDataLength[ 116643, 2563] NotebookOptionsPosition[ 113558, 2498] NotebookOutlinePosition[ 114103, 2517] CellTagsIndexPosition[ 114060, 2514] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell["Kinetics defs 2022", "Section", CellChangeTimes->{{3.542948801634627*^9, 3.5429488340054784`*^9}, { 3.54294888700751*^9, 3.5429488892326374`*^9}, {3.603154693094841*^9, 3.603154707836867*^9}, 3.6044766362316465`*^9, {3.666884612660963*^9, 3.666884616329775*^9}, {3.6991607763526573`*^9, 3.6991607820826654`*^9}, { 3.730426901167432*^9, 3.730426903433354*^9}, 3.761887783998542*^9, { 3.792556151930075*^9, 3.7925561617596765`*^9}, {3.7926459798323584`*^9, 3.792645984158855*^9}, 3.8245242759641743`*^9, {3.8560529418040543`*^9, 3.856052943338414*^9}},ExpressionUUID->"f80f1a14-5e19-4d5b-b3db-\ 8ffc5f31ec0d"], Cell[CellGroupData[{ Cell["Constitutive equation defs", "Subsection", CellChangeTimes->{ 3.4710857577642*^9, {3.4710858112098*^9, 3.4710858239082003`*^9}, { 3.4710859049969997`*^9, 3.4710859061046*^9}, {3.4710892678536*^9, 3.4710892701*^9}, 3.4710893817648*^9, {3.4710894793115997`*^9, 3.4710895039596*^9}, {3.4710905118050003`*^9, 3.471090512741*^9}, { 3.4725252603414*^9, 3.472525264101*^9}, {3.4725253064862003`*^9, 3.4725253370934*^9}, {3.4747629998582*^9, 3.4747630055678*^9}, { 3.4752273878873997`*^9, 3.4752273904926*^9}, {3.4753663741484003`*^9, 3.4753664116352*^9}, {3.4753773567911997`*^9, 3.475377356916*^9}, 3.4753775870472*^9, {3.4754620951394*^9, 3.4754621031266003`*^9}, { 3.4779795207265997`*^9, 3.4779795250576*^9}, 3.537621939388029*^9, 3.537622934627954*^9, {3.5376230222229643`*^9, 3.537623029743394*^9}, { 3.537623460058007*^9, 3.537623470980632*^9}, {3.537624955779557*^9, 3.5376249594157653`*^9}, {3.6000028730424147`*^9, 3.6000028738224154`*^9}, {3.600051011531698*^9, 3.6000510174285083`*^9}, { 3.600051611598405*^9, 3.600051612175606*^9}, {3.601345436806925*^9, 3.6013454383201275`*^9}, {3.601345528035885*^9, 3.6013455313118906`*^9}, { 3.601345570093559*^9, 3.601345572137162*^9}, {3.601345611714432*^9, 3.6013456172368417`*^9}, 3.6016459822022843`*^9, {3.6031079936638145`*^9, 3.603107998499823*^9}, {3.603155726549856*^9, 3.6031557303874626`*^9}, { 3.604483923857771*^9, 3.6044839281321783`*^9}, {3.6044841812582226`*^9, 3.6044841822098246`*^9}, {3.6668845995307226`*^9, 3.666884605295741*^9}, { 3.666884878599985*^9, 3.666884887130411*^9}, {3.667525927678743*^9, 3.6675259326551514`*^9}, {3.667526020327306*^9, 3.6675260219809084`*^9}, { 3.7928123912503653`*^9, 3.7928124033550367`*^9}},ExpressionUUID->"16ef144c-4d13-402e-96b4-\ c4e44eec02f8"], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"elastic", " ", "solid"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"ESOLID", "[", RowBox[{"EE_", ",", "\[Nu]\[Nu]_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"e", ",", "ee", ",", "eee", ",", "EM", ",", "GM"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"e", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"ee", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Outer", "[", RowBox[{"Times", ",", RowBox[{"e", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"e", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"eee", "=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"ee", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}], "+", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{"ee", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}], "+", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{"ee", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}], "+", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"1", ",", "3"}], "]"}], "]"}]}]}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"EM", "=", RowBox[{ RowBox[{"Inverse", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", RowBox[{"-", "\[Nu]\[Nu]"}], ",", RowBox[{"-", "\[Nu]\[Nu]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Nu]\[Nu]"}], ",", "1", ",", RowBox[{"-", "\[Nu]\[Nu]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Nu]\[Nu]"}], ",", RowBox[{"-", "\[Nu]\[Nu]"}], ",", "1"}], "}"}]}], "}"}], "]"}], "*", "EE"}]}], ";", "\[IndentingNewLine]", RowBox[{"GM", "=", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}], "*", RowBox[{"EE", "/", RowBox[{"(", RowBox[{"2", "+", RowBox[{"2", "*", "\[Nu]\[Nu]"}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"EM", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "*", RowBox[{"Outer", "[", RowBox[{"Times", ",", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"i", ",", "i"}], "]"}], "]"}], ",", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"j", ",", "j"}], "]"}], "]"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], "+", "\[IndentingNewLine]", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"GM", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "*", RowBox[{"Outer", "[", RowBox[{"Times", ",", RowBox[{"eee", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"eee", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}]}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"elastic", " ", "plate"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"EPLATE", "[", RowBox[{"EE_", ",", "\[Nu]\[Nu]_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"e", ",", "ee", ",", "eee", ",", "EM", ",", "GM"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"e", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"ee", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Outer", "[", RowBox[{"Times", ",", RowBox[{"e", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"e", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"eee", "=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"ee", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}], "+", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{"ee", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}], "+", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{"ee", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}], "+", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"1", ",", "3"}], "]"}], "]"}]}]}], "}"}]}], ";", "\[IndentingNewLine]", RowBox[{"EM", "=", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "\[Nu]\[Nu]", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Nu]\[Nu]", ",", "1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], "}"}], "*", RowBox[{"EE", "/", RowBox[{"(", RowBox[{"1", "-", RowBox[{"\[Nu]\[Nu]", "^", "2"}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{"GM", "=", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}], "*", RowBox[{"EE", "/", RowBox[{"(", RowBox[{"2", "+", RowBox[{"2", "*", "\[Nu]\[Nu]"}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"EM", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "*", RowBox[{"Outer", "[", RowBox[{"Times", ",", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"i", ",", "i"}], "]"}], "]"}], ",", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"j", ",", "j"}], "]"}], "]"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], "+", "\[IndentingNewLine]", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"GM", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "*", RowBox[{"Outer", "[", RowBox[{"Times", ",", RowBox[{"eee", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"eee", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}]}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{ RowBox[{"\[IndentingNewLine]", RowBox[{"(*", RowBox[{"elastic", " ", "beam"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"EBEAM", "[", RowBox[{"EE_", ",", "\[Nu]\[Nu]_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"e", ",", "ee", ",", "EM", ",", "GM"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"e", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}]}], ";", " ", "\[IndentingNewLine]", RowBox[{"ee", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Outer", "[", RowBox[{"Times", ",", "e1", ",", "e2"}], "]"}], ",", RowBox[{"{", RowBox[{"e1", ",", "e"}], "}"}], ",", RowBox[{"{", RowBox[{"e2", ",", "e"}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"EM", "=", RowBox[{"EE", "*", RowBox[{"Inverse", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", RowBox[{"-", "\[Nu]\[Nu]"}], ",", RowBox[{"-", "\[Nu]\[Nu]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Nu]\[Nu]"}], ",", "1", ",", RowBox[{"-", "\[Nu]\[Nu]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", "\[Nu]\[Nu]"}], ",", RowBox[{"-", "\[Nu]\[Nu]"}], ",", "1"}], "}"}]}], "}"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"GM", "=", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "1"}], "}"}]}], "}"}], "*", RowBox[{"EE", "/", RowBox[{"(", RowBox[{"2", "+", RowBox[{"2", "*", "\[Nu]\[Nu]"}]}], ")"}]}]}]}], ";", "\[IndentingNewLine]", RowBox[{ RowBox[{"Simplify", "[", "\[IndentingNewLine]", RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"EM", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}], "*", RowBox[{"Outer", "[", RowBox[{"Times", ",", RowBox[{"e", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"j", ",", "j"}], "]"}], "]"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], "+", "\[IndentingNewLine]", RowBox[{ RowBox[{"GM", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], "*", RowBox[{"(", RowBox[{ RowBox[{"Outer", "[", RowBox[{"Times", ",", RowBox[{"ee", "[", RowBox[{"[", "1", "]"}], "]"}], ",", RowBox[{"ee", "[", RowBox[{"[", "2", "]"}], "]"}]}], "]"}], "+", RowBox[{"Outer", "[", RowBox[{"Times", ",", RowBox[{"ee", "[", RowBox[{"[", "2", "]"}], "]"}], ",", RowBox[{"ee", "[", RowBox[{"[", "1", "]"}], "]"}]}], "]"}]}]}]}]}], "]"}], " ", "+", RowBox[{ RowBox[{"GM", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}], "*", RowBox[{"Outer", "[", RowBox[{"Times", ",", RowBox[{ RowBox[{"ee", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}], "+", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{"ee", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}], "+", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}]}], "]"}]}], " ", "+", "\[IndentingNewLine]", RowBox[{ RowBox[{"GM", "[", RowBox[{"[", RowBox[{"3", ",", "3"}], "]"}], "]"}], "*", RowBox[{"Outer", "[", RowBox[{"Times", ",", RowBox[{ RowBox[{"ee", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}], "+", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"1", ",", "3"}], "]"}], "]"}]}], ",", RowBox[{ RowBox[{"ee", "[", RowBox[{"[", RowBox[{"3", ",", "1"}], "]"}], "]"}], "+", RowBox[{"ee", "[", RowBox[{"[", RowBox[{"1", ",", "3"}], "]"}], "]"}]}]}], "]"}]}]}]}]}], "]"}]}], "]"}], ";"}], " ", "*)"}]}]}]], "Input", InitializationCell->True, CellChangeTimes->{ 3.667525913825918*^9, 3.667525951874385*^9, {3.667526014180895*^9, 3.667526133989105*^9}, {3.6675278740286736`*^9, 3.6675279303758965`*^9}, { 3.6675452156286526`*^9, 3.6675452299144697`*^9}, 3.667559087584382*^9, { 3.7306191461426644`*^9, 3.7306192201817408`*^9}, {3.73061932111433*^9, 3.730619366837903*^9}, 3.730619435892133*^9, 3.7306194715677967`*^9, { 3.792556190513768*^9, 3.7925563019282928`*^9}, {3.7926460008279743`*^9, 3.792646002988407*^9}, {3.792646715621252*^9, 3.792646793692767*^9}, { 3.7926468336260786`*^9, 3.792646863732798*^9}, {3.792646919221965*^9, 3.7926469217462187`*^9}, {3.792646981602483*^9, 3.792647023293395*^9}, { 3.7926471113659925`*^9, 3.792647216354287*^9}, {3.7926472509556*^9, 3.792647329350418*^9}, 3.792647381625019*^9, {3.7926476107221904`*^9, 3.792647672247836*^9}, {3.7926478553866844`*^9, 3.79264814545584*^9}, { 3.7926482115981035`*^9, 3.7926482623345594`*^9}, {3.792648323009469*^9, 3.792648500433292*^9}, {3.792648589495519*^9, 3.7926486573317585`*^9}, { 3.792648690372342*^9, 3.792648694702925*^9}, {3.7926493608296776`*^9, 3.7926493631972895`*^9}, 3.792649923700011*^9, 3.7926500450992155`*^9, { 3.7928124160498257`*^9, 3.7928124337486877`*^9}, {3.84243210085151*^9, 3.842432101128209*^9}, {3.8424322071124315`*^9, 3.842432352128103*^9}, { 3.842432415616482*^9, 3.8424324409940157`*^9}, {3.84243253225478*^9, 3.8424326191061816`*^9}, {3.8424326787318554`*^9, 3.842432694798433*^9}}, CellLabel->"In[60]:=",ExpressionUUID->"718a3088-f79e-4e0f-855c-1eb4a3539f04"] }, Closed]], Cell[CellGroupData[{ Cell["Kinematics defs", "Subsection", CellChangeTimes->{ 3.4710857577642*^9, {3.4710858112098*^9, 3.4710858239082003`*^9}, { 3.4710859049969997`*^9, 3.4710859061046*^9}, {3.4710892678536*^9, 3.4710892701*^9}, 3.4710893817648*^9, {3.4710894793115997`*^9, 3.4710895039596*^9}, {3.4710905118050003`*^9, 3.471090512741*^9}, { 3.4725252603414*^9, 3.472525264101*^9}, {3.4725253064862003`*^9, 3.4725253370934*^9}, {3.4747629998582*^9, 3.4747630055678*^9}, { 3.4752273878873997`*^9, 3.4752273904926*^9}, {3.4753663741484003`*^9, 3.4753664116352*^9}, {3.4753773567911997`*^9, 3.475377356916*^9}, 3.4753775870472*^9, {3.4754620951394*^9, 3.4754621031266003`*^9}, { 3.4779795207265997`*^9, 3.4779795250576*^9}, 3.537621939388029*^9, 3.537622934627954*^9, {3.5376230222229643`*^9, 3.537623029743394*^9}, { 3.537623460058007*^9, 3.537623470980632*^9}, {3.537624955779557*^9, 3.5376249594157653`*^9}, {3.6000028730424147`*^9, 3.6000028738224154`*^9}, {3.600051011531698*^9, 3.6000510174285083`*^9}, { 3.600051611598405*^9, 3.600051612175606*^9}, {3.601345436806925*^9, 3.6013454383201275`*^9}, {3.601345528035885*^9, 3.6013455313118906`*^9}, { 3.601345570093559*^9, 3.601345572137162*^9}, {3.601345611714432*^9, 3.6013456172368417`*^9}, 3.6016459822022843`*^9, {3.6031079936638145`*^9, 3.603107998499823*^9}, {3.603155726549856*^9, 3.6031557303874626`*^9}, { 3.604483923857771*^9, 3.6044839281321783`*^9}, {3.6044841812582226`*^9, 3.6044841822098246`*^9}, {3.6668845995307226`*^9, 3.666884605295741*^9}, { 3.666884878599985*^9, 3.666884887130411*^9}, {3.667525927678743*^9, 3.6675259326551514`*^9}, {3.667526020327306*^9, 3.6675260219809084`*^9}, { 3.7928123912503653`*^9, 3.7928124033550367`*^9}, {3.7943964597734575`*^9, 3.794396462038081*^9}},ExpressionUUID->"a7d5c6f1-3a2f-43a7-a35a-\ 4129cabb329b"], Cell[BoxData[{ RowBox[{ RowBox[{ RowBox[{"Clear", "[", RowBox[{ "F", ",", "FF", ",", "HH", ",", "DF", ",", "GG", ",", "GRAD", ",", "DIV"}], "]"}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"HH", "[", RowBox[{"{", RowBox[{"rr_", ",", "var_", ",", "cnd_"}], "}"}], "]"}], ":=", RowBox[{"Thread", "[", RowBox[{ RowBox[{"Flatten", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"H", "[", RowBox[{"i", ",", "j"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}]}], "]"}], "]"}], " ", "\[Rule]", RowBox[{"Flatten", "[", RowBox[{"Map", "[", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[PartialD]", "#"], "rr"}], "&"}], ",", "var"}], "]"}], "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"hh", "[", RowBox[{"{", RowBox[{"rr_", ",", "var_", ",", "cnd_"}], "}"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", "HHH", "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"HHH", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"H", "[", RowBox[{"i", ",", "j"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}]}], "]"}], " ", "/.", RowBox[{"HH", "[", RowBox[{"{", RowBox[{"rr", ",", "var", ",", "cnd"}], "}"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Thread", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"h", "[", "i", "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}]}], "]"}], "->", RowBox[{"Map", "[", RowBox[{ RowBox[{ RowBox[{"Simplify", "[", RowBox[{ RowBox[{"Sqrt", "[", RowBox[{"#", ".", "#"}], "]"}], ",", "cnd"}], "]"}], "&"}], ",", "HHH"}], "]"}]}], "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"FF", "[", RowBox[{"{", RowBox[{"rr_", ",", "var_", ",", "cnd_"}], "}"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"HHH", ",", "hhh"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"HHH", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"H", "[", RowBox[{"i", ",", "j"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}]}], "]"}], " ", "/.", RowBox[{"HH", "[", RowBox[{"{", RowBox[{"rr", ",", "var", ",", "cnd"}], "}"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"hhh", " ", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"h", "[", "i", "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}]}], "]"}], "/.", RowBox[{"hh", "[", RowBox[{"{", RowBox[{"rr", ",", "var", ",", "cnd"}], "}"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Thread", "[", RowBox[{ RowBox[{"Flatten", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"F", "[", RowBox[{"i", ",", "j"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}]}], "]"}], "]"}], "\[Rule]", " ", RowBox[{"Flatten", "[", RowBox[{"Simplify", "[", RowBox[{ RowBox[{ RowBox[{"DiagonalMatrix", "[", RowBox[{"1", "/", "hhh"}], "]"}], ".", "HHH"}], ",", "cnd"}], "]"}], "]"}]}], "]"}]}]}], "]"}]}], ";"}], " ", "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"DF", "[", RowBox[{"{", RowBox[{"rr_", ",", "var_", ",", "cnd_"}], "}"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"FFF", ",", "FFFINV"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"FFF", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"F", "[", RowBox[{"i", ",", "j"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}]}], "]"}], " ", "/.", RowBox[{"FF", "[", RowBox[{"{", RowBox[{"rr", ",", "var", ",", "cnd"}], "}"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"FFFINV", "=", RowBox[{"Inverse", "[", "FFF", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Thread", "[", RowBox[{ RowBox[{"Flatten", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"dF", "[", RowBox[{"i", ",", "j", ",", "k"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"k", ",", "var"}], "}"}]}], "]"}], "]"}], "\[Rule]", RowBox[{"Flatten", "[", RowBox[{"Simplify", "[", RowBox[{ RowBox[{"Map", "[", RowBox[{ RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ SubscriptBox["\[PartialD]", "#"], "FFF"}], ")"}], ".", "FFFINV"}], "&"}], ",", "var"}], "]"}], ",", "cnd"}], "]"}], "]"}]}], "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"GG", "[", RowBox[{"{", RowBox[{"rr_", ",", "var_", ",", "cnd_"}], "}"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"FFF", ",", "HHH"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"FFF", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"F", "[", RowBox[{"i", ",", "j"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}]}], "]"}], " ", "/.", RowBox[{"FF", "[", RowBox[{"{", RowBox[{"rr", ",", "var", ",", "cnd"}], "}"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"HHH", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"H", "[", RowBox[{"i", ",", "j"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}]}], "]"}], " ", "/.", RowBox[{"HH", "[", RowBox[{"{", RowBox[{"rr", ",", "var", ",", "cnd"}], "}"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Thread", "[", RowBox[{ RowBox[{"Flatten", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"G", "[", RowBox[{"i", ",", "j"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}]}], "]"}], "]"}], "->", RowBox[{"Flatten", "[", RowBox[{"Simplify", "[", RowBox[{ RowBox[{"Inverse", "[", RowBox[{"HHH", ".", " ", RowBox[{"Transpose", "[", "FFF", "]"}]}], "]"}], ",", "cnd"}], "]"}], "]"}]}], "]"}]}]}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"\[Kappa]\[Kappa]", "[", RowBox[{"rr_", ",", "var_", ",", "cnd_"}], "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"GGT", ",", "DFT", ",", "GN", ",", "en"}], "}"}], ",", " ", "\[IndentingNewLine]", RowBox[{ RowBox[{"DFT", "=", RowBox[{"DF", "[", RowBox[{"rr", ",", "var", ",", "cnd"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"en", " ", "=", RowBox[{"D", "[", RowBox[{"rr", ",", RowBox[{"var", "[", RowBox[{"[", "2", "]"}], "]"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"GN", "=", RowBox[{"If", "[", RowBox[{ RowBox[{"AllTrue", "[", RowBox[{ RowBox[{"PossibleZeroQ", "[", "e\[Beta]", "]"}], ",", RowBox[{"#", "&"}]}], "]"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"DFT", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}], ",", RowBox[{"DFT", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}], ",", RowBox[{"DFT", "[", RowBox[{"[", RowBox[{"3", ",", "2"}], "]"}], "]"}]}], "}"}], ",", "\[IndentingNewLine]", RowBox[{"{", RowBox[{ RowBox[{"DFT", "[", RowBox[{"[", RowBox[{"1", ",", "3"}], "]"}], "]"}], ",", RowBox[{"DFT", "[", RowBox[{"[", RowBox[{"2", ",", "3"}], "]"}], "]"}], ",", RowBox[{"DFT", "[", RowBox[{"[", RowBox[{"3", ",", "3"}], "]"}], "]"}]}], "}"}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"GGT", "=", RowBox[{"GG", "[", RowBox[{"r", ",", "var", ",", "cnd"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"Return", "[", RowBox[{"Simplify", "[", RowBox[{"Transpose", "[", RowBox[{"GGT", ".", "GN"}], "]"}], "]"}], "]"}]}]}], "]"}]}], ";"}]}], "Input", InitializationCell->True, CellChangeTimes->{{3.6044841869366326`*^9, 3.6044842193846893`*^9}, { 3.60448429961563*^9, 3.604484425679452*^9}, 3.604484459219511*^9, { 3.6046620267627306`*^9, 3.604662028993534*^9}, {3.6048381720233164`*^9, 3.6048381756737227`*^9}, {3.666868290197764*^9, 3.6668683122058334`*^9}, { 3.6668683597467775`*^9, 3.6668683818964453`*^9}, {3.66686845639587*^9, 3.666868470308913*^9}, {3.666877098337613*^9, 3.666877181422058*^9}, { 3.6675282960618124`*^9, 3.6675283148698883`*^9}, {3.6675444752493052`*^9, 3.6675444758283386`*^9}, {3.667555712966066*^9, 3.6675557974839*^9}, { 3.6675561975497828`*^9, 3.667556197880801*^9}, {3.668252663901308*^9, 3.6682526700009193`*^9}, {3.6683165093185673`*^9, 3.668316509802168*^9}, 3.668316703897709*^9, {3.7618881447895975`*^9, 3.7618881491945267`*^9}, 3.7618882178660393`*^9, {3.7618882741028605`*^9, 3.761888277789795*^9}, { 3.7618916470099993`*^9, 3.7618916834392595`*^9}, 3.7618917639517007`*^9, { 3.7618919844468417`*^9, 3.761891984649885*^9}, {3.761892018345189*^9, 3.76189208954007*^9}, {3.761901785232708*^9, 3.7619018131345997`*^9}, { 3.761901904392881*^9, 3.761901904785979*^9}, {3.761902018863076*^9, 3.7619020338496556`*^9}, {3.761902399886441*^9, 3.761902402396004*^9}, { 3.7619026504706573`*^9, 3.7619026636835194`*^9}, {3.7619059915496955`*^9, 3.761905994645768*^9}, {3.761906262916795*^9, 3.76190627890534*^9}, { 3.7619064043186855`*^9, 3.761906453396896*^9}, {3.792556354018784*^9, 3.792556466846212*^9}, 3.7925565150676746`*^9, {3.7925565486936865`*^9, 3.7925567497115335`*^9}, {3.7925567825923357`*^9, 3.792556783363819*^9}, { 3.792556837086212*^9, 3.79255699870599*^9}, 3.7925570855341873`*^9, 3.7925571536529093`*^9, {3.7925571847033415`*^9, 3.792557576394307*^9}, { 3.792557614184314*^9, 3.792557618568714*^9}, {3.7925576591535816`*^9, 3.7925577052436485`*^9}, {3.7925578691672688`*^9, 3.792558448892994*^9}, { 3.792558581176793*^9, 3.7925586490883245`*^9}, {3.792558780174613*^9, 3.7925588005876756`*^9}, {3.7925588416690073`*^9, 3.7925588576355443`*^9}, {3.7925592749058166`*^9, 3.7925593555929585`*^9}, {3.792560913853984*^9, 3.7925610021344237`*^9}, { 3.7925610326253443`*^9, 3.7925610575251074`*^9}, {3.792561113591499*^9, 3.792561155917766*^9}, {3.792561845232957*^9, 3.7925619041510344`*^9}, { 3.7925638292356277`*^9, 3.792563985063381*^9}, {3.792564082099863*^9, 3.7925640979076595`*^9}, {3.792564137326541*^9, 3.7925642253326325`*^9}, { 3.792564264685545*^9, 3.7925643012659616`*^9}, {3.792564342645465*^9, 3.7925643627598143`*^9}, {3.7925644271683135`*^9, 3.792564429236191*^9}, { 3.7925645301916933`*^9, 3.7925645402830744`*^9}, {3.7925646092162127`*^9, 3.792564618505151*^9}, {3.792564652915941*^9, 3.7925647334391403`*^9}, 3.792564765383103*^9, {3.792564801734785*^9, 3.792564802731366*^9}, { 3.7925651879037533`*^9, 3.792565199057434*^9}, {3.7925670381837444`*^9, 3.792567050647067*^9}, {3.792567154812933*^9, 3.7925672131939507`*^9}, { 3.7925672602503347`*^9, 3.7925673758752804`*^9}, {3.792567496290472*^9, 3.792567573807705*^9}, {3.7925676048022065`*^9, 3.7925676480552025`*^9}, { 3.7925677589260807`*^9, 3.7925677901653624`*^9}, {3.792567876047678*^9, 3.7925679051049666`*^9}, 3.792567965501421*^9, {3.7925680094757185`*^9, 3.7925680224511566`*^9}, {3.792568053498392*^9, 3.7925681027827396`*^9}, { 3.7925681436539965`*^9, 3.792568179477368*^9}, {3.7925682425561633`*^9, 3.7925682542817435`*^9}, {3.7925683015096903`*^9, 3.792568337794861*^9}, { 3.792568410590354*^9, 3.7925684585433517`*^9}, {3.792568500069618*^9, 3.792568544260894*^9}, {3.7925685830060225`*^9, 3.792568606765068*^9}, { 3.79256866734807*^9, 3.7925687120921717`*^9}, {3.7925687483801003`*^9, 3.792568758891363*^9}, {3.7925688054327574`*^9, 3.792568806433818*^9}, { 3.7925688434385233`*^9, 3.7925688509576387`*^9}, {3.7925688816057577`*^9, 3.79256888601202*^9}, {3.7925689708945107`*^9, 3.792569003050635*^9}, { 3.792569092363428*^9, 3.792569137513811*^9}, {3.792569191561468*^9, 3.792569334535511*^9}, {3.7925694282303543`*^9, 3.792569482118329*^9}, { 3.79256953209522*^9, 3.7925697481922584`*^9}, 3.792569808336525*^9, { 3.7925698526772623`*^9, 3.792569876772537*^9}, {3.7925699286646357`*^9, 3.7925699317716227`*^9}, {3.7925699673075094`*^9, 3.7925703455642757`*^9}, {3.7925704404381704`*^9, 3.7925704577130213`*^9}, 3.7925705245001535`*^9, {3.7925706575427184`*^9, 3.7925706739489913`*^9}, { 3.792570713466522*^9, 3.792570742092165*^9}, 3.79257087493649*^9, 3.792570948822445*^9, {3.7925709833973784`*^9, 3.792570983627946*^9}, { 3.792571383695715*^9, 3.792571417201741*^9}, {3.7925715287632074`*^9, 3.7925716890877604`*^9}, {3.7925718082846317`*^9, 3.7925718896185102`*^9}, {3.792590997795863*^9, 3.79259099855244*^9}, { 3.7925910400445538`*^9, 3.7925910817225924`*^9}, 3.792591130587617*^9, { 3.7925911975707383`*^9, 3.7925912333964806`*^9}, {3.7925949482023687`*^9, 3.792594957289979*^9}, {3.792595025280615*^9, 3.792595030882008*^9}, { 3.7925950712628264`*^9, 3.7925950800117626`*^9}, {3.7925951872348475`*^9, 3.7925951915721235`*^9}, {3.792635717158695*^9, 3.79263576744442*^9}, { 3.792635829712722*^9, 3.7926359217639484`*^9}, {3.7926359620309753`*^9, 3.792636140458949*^9}, {3.7926362081677284`*^9, 3.7926362587135468`*^9}, { 3.79263636347906*^9, 3.792636382892848*^9}, {3.792636842824526*^9, 3.7926371232319517`*^9}, {3.792637219820684*^9, 3.7926372259977136`*^9}, { 3.792637514141102*^9, 3.7926375176812363`*^9}, {3.7926378322869935`*^9, 3.7926378362954063`*^9}, {3.7926379061322355`*^9, 3.7926379723075356`*^9}, {3.792638249529179*^9, 3.7926382737812505`*^9}, { 3.7926383101353693`*^9, 3.7926383440075893`*^9}, {3.792638402502085*^9, 3.7926384078678255`*^9}, 3.7926384548231983`*^9, 3.7926384993447423`*^9, { 3.792638534264426*^9, 3.792638535105315*^9}, {3.7926386277631364`*^9, 3.792638628039978*^9}, {3.792639096007407*^9, 3.792639169501439*^9}, { 3.792639280134539*^9, 3.792639311964221*^9}, {3.7926394283350706`*^9, 3.7926394483726883`*^9}, {3.792639703218527*^9, 3.792639723938527*^9}, { 3.792639767781361*^9, 3.792639768071415*^9}, {3.7926422179783792`*^9, 3.792642302493026*^9}, {3.7926425281846733`*^9, 3.7926425386729965`*^9}, { 3.792643430070241*^9, 3.792643458696355*^9}, {3.792643493295534*^9, 3.792643496908493*^9}, {3.792643624400555*^9, 3.792643777959143*^9}, { 3.792643827290126*^9, 3.792643860120751*^9}, {3.792644641764208*^9, 3.7926446853323097`*^9}, {3.7926452058923297`*^9, 3.792645375276899*^9}, { 3.7926454233150563`*^9, 3.7926454244698086`*^9}, {3.7926456987270765`*^9, 3.792645715973406*^9}, {3.7926460230589604`*^9, 3.7926460552369814`*^9}, { 3.7926560559283648`*^9, 3.7926560942680063`*^9}, {3.7926565894328995`*^9, 3.7926565928040247`*^9}, {3.792657705491242*^9, 3.7926577384722776`*^9}, { 3.7926766587432857`*^9, 3.7926766722659645`*^9}, {3.792676746206684*^9, 3.792676751416193*^9}, {3.7928126377010927`*^9, 3.792812638325923*^9}, { 3.7943973400145607`*^9, 3.7943973656987944`*^9}, 3.794397434029624*^9, { 3.7943974991666927`*^9, 3.7943975061932535`*^9}, {3.7943978360324554`*^9, 3.794397839768346*^9}, {3.794397891267982*^9, 3.794397892284954*^9}}, CellLabel->"In[62]:=",ExpressionUUID->"58c30491-20f6-423f-9696-44262d8befe7"], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"gradient", " ", "and", " ", "divergence", " ", "operators"}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{ RowBox[{"GRAD", "[", RowBox[{"{", RowBox[{"rr_", ",", "var_", ",", "cnd_"}], "}"}], "]"}], "[", "sca_", "]"}], ":=", RowBox[{ RowBox[{"Module", "[", RowBox[{ RowBox[{"{", "GGG", "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"GGG", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"G", "[", RowBox[{"i", ",", "j"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}]}], "]"}], " ", "/.", RowBox[{"GG", "[", RowBox[{"{", RowBox[{"rr", ",", "var", ",", "cnd"}], "}"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"GGG", ".", RowBox[{"Map", "[", RowBox[{ RowBox[{ RowBox[{ SubscriptBox["\[PartialD]", "#"], "sca"}], "&"}], " ", ",", "var"}], "]"}]}]}]}], "]"}], "/;", RowBox[{"Not", "[", RowBox[{"ArrayQ", "[", "sca", "]"}], "]"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"GRAD", "[", RowBox[{"{", RowBox[{"rr_", ",", "var_", ",", "cnd_"}], "}"}], "]"}], "[", "vec_", "]"}], ":=", RowBox[{ RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"DFF", ",", "GGG"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"DFF", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"dF", "[", RowBox[{"i", ",", "j", ",", "k"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"k", ",", "var"}], "}"}]}], "]"}], "/.", RowBox[{"DF", "[", RowBox[{"{", RowBox[{"rr", ",", "var", ",", "cnd"}], "}"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"GGG", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"G", "[", RowBox[{"i", ",", "j"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}]}], "]"}], " ", "/.", RowBox[{"GG", "[", RowBox[{"{", RowBox[{"rr", ",", "var", ",", "cnd"}], "}"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"GGG", "[", RowBox[{"[", RowBox[{"i", ",", "k"}], "]"}], "]"}], "*", RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"var", "[", RowBox[{"[", "k", "]"}], "]"}]], RowBox[{"vec", "[", RowBox[{"[", "j", "]"}], "]"}]}]}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "3"}], "}"}]}], "]"}], "+", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"GGG", "[", RowBox[{"[", RowBox[{"i", ",", "l"}], "]"}], "]"}], "*", RowBox[{"vec", "[", RowBox[{"[", "k", "]"}], "]"}], "*", RowBox[{"DFF", "[", RowBox[{"[", RowBox[{"l", ",", "k", ",", "j"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"l", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "3"}], "}"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}]}]}], "]"}], "/;", " ", RowBox[{"VectorQ", "[", "vec", "]"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"GRAD0", "[", RowBox[{"{", RowBox[{"rr_", ",", "var_", ",", "cnd_"}], "}"}], "]"}], "[", "vec_", "]"}], ":=", RowBox[{ RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"DFF", ",", "GGG"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"DFF", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"dF", "[", RowBox[{"i", ",", "j", ",", "k"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"k", ",", "var"}], "}"}]}], "]"}], "/.", RowBox[{"DF", "[", RowBox[{"{", RowBox[{"rr", ",", "var", ",", "cnd"}], "}"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"GGG", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"G", "[", RowBox[{"i", ",", "j"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}]}], "]"}], " ", "/.", RowBox[{"GG", "[", RowBox[{"{", RowBox[{"rr", ",", "var", ",", "cnd"}], "}"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"GGG", "[", RowBox[{"[", RowBox[{"i", ",", "k"}], "]"}], "]"}], "/.", RowBox[{"n", "\[Rule]", "0"}]}], ")"}], "*", RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"var", "[", RowBox[{"[", "k", "]"}], "]"}]], RowBox[{"vec", "[", RowBox[{"[", "j", "]"}], "]"}]}]}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "3"}], "}"}]}], "]"}], "+", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"GGG", "[", RowBox[{"[", RowBox[{"i", ",", "l"}], "]"}], "]"}], "/.", RowBox[{"n", "\[Rule]", "0"}]}], ")"}], "*", RowBox[{"vec", "[", RowBox[{"[", "k", "]"}], "]"}], "*", RowBox[{"DFF", "[", RowBox[{"[", RowBox[{"l", ",", "k", ",", "j"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"l", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "3"}], "}"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}]}]}], "]"}], "/;", " ", RowBox[{"VectorQ", "[", "vec", "]"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"DIV", "[", RowBox[{"{", RowBox[{"rr_", ",", "var_", ",", "cnd_"}], "}"}], "]"}], "[", "mat_", "]"}], ":=", RowBox[{ RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"DFF", ",", "GGG"}], "}"}], ",", "\[IndentingNewLine]", RowBox[{ RowBox[{"DFF", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"dF", "[", RowBox[{"i", ",", "j", ",", "k"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"k", ",", "var"}], "}"}]}], "]"}], "/.", RowBox[{"DF", "[", RowBox[{"{", RowBox[{"rr", ",", "var", ",", "cnd"}], "}"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"GGG", "=", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{"G", "[", RowBox[{"i", ",", "j"}], "]"}], ",", RowBox[{"{", RowBox[{"i", ",", "var"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "var"}], "}"}]}], "]"}], " ", "/.", RowBox[{"GG", "[", RowBox[{"{", RowBox[{"rr", ",", "var", ",", "cnd"}], "}"}], "]"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"GGG", "[", RowBox[{"[", RowBox[{"j", ",", "k"}], "]"}], "]"}], "*", RowBox[{"(", RowBox[{ SubscriptBox["\[PartialD]", RowBox[{"var", "[", RowBox[{"[", "k", "]"}], "]"}]], RowBox[{"mat", "[", RowBox[{"[", RowBox[{"j", ",", "i"}], "]"}], "]"}]}], ")"}]}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "3"}], "}"}]}], "]"}], "+", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"GGG", "[", RowBox[{"[", RowBox[{"k", ",", "l"}], "]"}], "]"}], "*", RowBox[{"mat", "[", RowBox[{"[", RowBox[{"j", ",", "i"}], "]"}], "]"}], "*", RowBox[{"DFF", "[", RowBox[{"[", RowBox[{"l", ",", "j", ",", "k"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"l", ",", "1", ",", "3"}], "}"}]}], "]"}], "\[IndentingNewLine]", "+", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"GGG", "[", RowBox[{"[", RowBox[{"k", ",", "l"}], "]"}], "]"}], "*", RowBox[{"mat", "[", RowBox[{"[", RowBox[{"k", ",", "j"}], "]"}], "]"}], "*", RowBox[{"DFF", "[", RowBox[{"[", RowBox[{"l", ",", "j", ",", "i"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"k", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"l", ",", "1", ",", "3"}], "}"}]}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}]}]}], "]"}], " ", "/;", RowBox[{"MatrixQ", "[", "mat", "]"}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"double", " ", "inner", " ", "product"}], "*)"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"DDOT", "[", RowBox[{"a_", ",", "b_"}], "]"}], ":=", "\[IndentingNewLine]", RowBox[{ RowBox[{"TensorContract", "[", RowBox[{ RowBox[{"a", ".", "b"}], ",", RowBox[{"{", RowBox[{"3", ",", "4"}], "}"}]}], "]"}], "/;", RowBox[{ RowBox[{ RowBox[{"TensorRank", "[", "a", "]"}], "\[Equal]", "4"}], "&&", RowBox[{ RowBox[{"TensorRank", "[", "b", "]"}], "\[Equal]", "2"}]}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"DDOT", "[", RowBox[{"a_", ",", "b_"}], "]"}], ":=", "\[IndentingNewLine]", RowBox[{ RowBox[{"TensorContract", "[", RowBox[{ RowBox[{"a", ".", "b"}], ",", RowBox[{"{", RowBox[{"1", ",", "2"}], "}"}]}], "]"}], "/;", RowBox[{ RowBox[{ RowBox[{"TensorRank", "[", "a", "]"}], "\[Equal]", "2"}], "&&", RowBox[{ RowBox[{"TensorRank", "[", "b", "]"}], "\[Equal]", "4"}]}]}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"DDOT", "[", RowBox[{"a_", ",", "b_"}], "]"}], ":=", "\[IndentingNewLine]", RowBox[{ RowBox[{"TensorContract", "[", RowBox[{ RowBox[{"a", ".", "b"}], ",", RowBox[{"{", RowBox[{"1", ",", "2"}], "}"}]}], "]"}], "/;", RowBox[{ RowBox[{ RowBox[{"TensorRank", "[", "a", "]"}], "\[Equal]", "2"}], "&&", RowBox[{ RowBox[{"TensorRank", "[", "b", "]"}], "\[Equal]", "2"}]}]}]}], ";"}]}]}]], "Input", InitializationCell->True, CellChangeTimes->{{3.6044841869366326`*^9, 3.6044842193846893`*^9}, { 3.60448429961563*^9, 3.604484425679452*^9}, 3.604484459219511*^9, { 3.6046620267627306`*^9, 3.604662028993534*^9}, {3.6048381720233164`*^9, 3.6048381756737227`*^9}, {3.666868290197764*^9, 3.6668683122058334`*^9}, { 3.6668683597467775`*^9, 3.6668683818964453`*^9}, {3.66686845639587*^9, 3.666868470308913*^9}, {3.666877098337613*^9, 3.666877181422058*^9}, { 3.6675282960618124`*^9, 3.6675283148698883`*^9}, {3.6675444752493052`*^9, 3.6675444758283386`*^9}, {3.667555712966066*^9, 3.6675557974839*^9}, { 3.6675561975497828`*^9, 3.667556197880801*^9}, {3.668252663901308*^9, 3.6682526700009193`*^9}, {3.6683165093185673`*^9, 3.668316509802168*^9}, 3.668316703897709*^9, {3.7618881447895975`*^9, 3.7618881491945267`*^9}, 3.7618882178660393`*^9, {3.7618882741028605`*^9, 3.761888277789795*^9}, { 3.7618916470099993`*^9, 3.7618916834392595`*^9}, 3.7618917639517007`*^9, { 3.7618919844468417`*^9, 3.761891984649885*^9}, {3.761892018345189*^9, 3.76189208954007*^9}, {3.761901785232708*^9, 3.7619018131345997`*^9}, { 3.761901904392881*^9, 3.761901904785979*^9}, {3.761902018863076*^9, 3.7619020338496556`*^9}, {3.761902399886441*^9, 3.761902402396004*^9}, { 3.7619026504706573`*^9, 3.7619026636835194`*^9}, {3.7619059915496955`*^9, 3.761905994645768*^9}, {3.761906262916795*^9, 3.76190627890534*^9}, { 3.7619064043186855`*^9, 3.761906453396896*^9}, {3.792556354018784*^9, 3.792556466846212*^9}, 3.7925565150676746`*^9, {3.7925565486936865`*^9, 3.7925567497115335`*^9}, {3.7925567825923357`*^9, 3.792556783363819*^9}, { 3.792556837086212*^9, 3.79255699870599*^9}, 3.7925570855341873`*^9, 3.7925571536529093`*^9, {3.7925571847033415`*^9, 3.792557576394307*^9}, { 3.792557614184314*^9, 3.792557618568714*^9}, {3.7925576591535816`*^9, 3.7925577052436485`*^9}, {3.7925578691672688`*^9, 3.792558448892994*^9}, { 3.792558581176793*^9, 3.7925586490883245`*^9}, {3.792558780174613*^9, 3.7925588005876756`*^9}, {3.7925588416690073`*^9, 3.7925588576355443`*^9}, {3.7925592749058166`*^9, 3.7925593555929585`*^9}, {3.792560913853984*^9, 3.7925610021344237`*^9}, { 3.7925610326253443`*^9, 3.7925610575251074`*^9}, {3.792561113591499*^9, 3.792561155917766*^9}, {3.792561845232957*^9, 3.7925619041510344`*^9}, { 3.7925638292356277`*^9, 3.792563985063381*^9}, {3.792564082099863*^9, 3.7925640979076595`*^9}, {3.792564137326541*^9, 3.7925642253326325`*^9}, { 3.792564264685545*^9, 3.7925643012659616`*^9}, {3.792564342645465*^9, 3.7925643627598143`*^9}, {3.7925644271683135`*^9, 3.792564429236191*^9}, { 3.7925645301916933`*^9, 3.7925645402830744`*^9}, {3.7925646092162127`*^9, 3.792564618505151*^9}, {3.792564652915941*^9, 3.7925647334391403`*^9}, 3.792564765383103*^9, {3.792564801734785*^9, 3.792564802731366*^9}, { 3.7925651879037533`*^9, 3.792565199057434*^9}, {3.7925670381837444`*^9, 3.792567050647067*^9}, {3.792567154812933*^9, 3.7925672131939507`*^9}, { 3.7925672602503347`*^9, 3.7925673758752804`*^9}, {3.792567496290472*^9, 3.792567573807705*^9}, {3.7925676048022065`*^9, 3.7925676480552025`*^9}, { 3.7925677589260807`*^9, 3.7925677901653624`*^9}, {3.792567876047678*^9, 3.7925679051049666`*^9}, 3.792567965501421*^9, {3.7925680094757185`*^9, 3.7925680224511566`*^9}, {3.792568053498392*^9, 3.7925681027827396`*^9}, { 3.7925681436539965`*^9, 3.792568179477368*^9}, {3.7925682425561633`*^9, 3.7925682542817435`*^9}, {3.7925683015096903`*^9, 3.792568337794861*^9}, { 3.792568410590354*^9, 3.7925684585433517`*^9}, {3.792568500069618*^9, 3.792568544260894*^9}, {3.7925685830060225`*^9, 3.792568606765068*^9}, { 3.79256866734807*^9, 3.7925687120921717`*^9}, {3.7925687483801003`*^9, 3.792568758891363*^9}, {3.7925688054327574`*^9, 3.792568806433818*^9}, { 3.7925688434385233`*^9, 3.7925688509576387`*^9}, {3.7925688816057577`*^9, 3.79256888601202*^9}, {3.7925689708945107`*^9, 3.792569003050635*^9}, { 3.792569092363428*^9, 3.792569137513811*^9}, {3.792569191561468*^9, 3.792569334535511*^9}, {3.7925694282303543`*^9, 3.792569482118329*^9}, { 3.79256953209522*^9, 3.7925697481922584`*^9}, 3.792569808336525*^9, { 3.7925698526772623`*^9, 3.792569876772537*^9}, {3.7925699286646357`*^9, 3.7925699317716227`*^9}, {3.7925699673075094`*^9, 3.7925703455642757`*^9}, {3.7925704404381704`*^9, 3.7925704577130213`*^9}, 3.7925705245001535`*^9, {3.7925706575427184`*^9, 3.7925706739489913`*^9}, { 3.792570713466522*^9, 3.792570742092165*^9}, 3.79257087493649*^9, 3.792570948822445*^9, {3.7925709833973784`*^9, 3.792570983627946*^9}, { 3.792571383695715*^9, 3.792571417201741*^9}, {3.7925715287632074`*^9, 3.7925716890877604`*^9}, {3.7925718082846317`*^9, 3.7925718896185102`*^9}, {3.792590997795863*^9, 3.79259099855244*^9}, { 3.7925910400445538`*^9, 3.7925910817225924`*^9}, 3.792591130587617*^9, { 3.7925911975707383`*^9, 3.7925912333964806`*^9}, {3.7925949482023687`*^9, 3.792594957289979*^9}, {3.792595025280615*^9, 3.792595030882008*^9}, { 3.7925950712628264`*^9, 3.7925950800117626`*^9}, {3.7925951872348475`*^9, 3.7925951915721235`*^9}, {3.792635717158695*^9, 3.79263576744442*^9}, { 3.792635829712722*^9, 3.7926359217639484`*^9}, {3.7926359620309753`*^9, 3.792636140458949*^9}, {3.7926362081677284`*^9, 3.7926362587135468`*^9}, { 3.79263636347906*^9, 3.792636382892848*^9}, {3.792636842824526*^9, 3.7926371232319517`*^9}, {3.792637219820684*^9, 3.7926372259977136`*^9}, { 3.792637514141102*^9, 3.7926375176812363`*^9}, {3.7926378322869935`*^9, 3.7926378362954063`*^9}, {3.7926379061322355`*^9, 3.7926379723075356`*^9}, {3.792638249529179*^9, 3.7926382737812505`*^9}, { 3.7926383101353693`*^9, 3.7926383440075893`*^9}, {3.792638402502085*^9, 3.7926384078678255`*^9}, 3.7926384548231983`*^9, 3.7926384993447423`*^9, { 3.792638534264426*^9, 3.792638535105315*^9}, {3.7926386277631364`*^9, 3.792638628039978*^9}, {3.792639096007407*^9, 3.792639169501439*^9}, { 3.792639280134539*^9, 3.792639311964221*^9}, {3.7926394283350706`*^9, 3.7926394483726883`*^9}, {3.792639703218527*^9, 3.792639723938527*^9}, { 3.792639767781361*^9, 3.792639768071415*^9}, {3.7926422179783792`*^9, 3.792642302493026*^9}, {3.7926425281846733`*^9, 3.7926425386729965`*^9}, { 3.792643430070241*^9, 3.792643458696355*^9}, {3.792643493295534*^9, 3.792643496908493*^9}, {3.792643624400555*^9, 3.792643777959143*^9}, { 3.792643827290126*^9, 3.792643860120751*^9}, {3.792644641764208*^9, 3.7926446853323097`*^9}, {3.7926452058923297`*^9, 3.792645375276899*^9}, { 3.7926454233150563`*^9, 3.7926454244698086`*^9}, {3.7926456987270765`*^9, 3.792645715973406*^9}, {3.7926460230589604`*^9, 3.7926460552369814`*^9}, { 3.7926560559283648`*^9, 3.7926560942680063`*^9}, {3.7926565894328995`*^9, 3.7926565928040247`*^9}, {3.792657705491242*^9, 3.7926577384722776`*^9}, { 3.7926766587432857`*^9, 3.7926766722659645`*^9}, {3.792676746206684*^9, 3.792676751416193*^9}, {3.7928126377010927`*^9, 3.792812638325923*^9}, { 3.7943973400145607`*^9, 3.7943973656987944`*^9}, 3.794397434029624*^9, { 3.7944488232294135`*^9, 3.794448874045909*^9}, {3.7944489122712097`*^9, 3.7944489178283005`*^9}, {3.794450507644116*^9, 3.7944505223033066`*^9}}, CellLabel->"In[69]:=",ExpressionUUID->"ab668980-2ce3-4ad7-95ec-4e45a51ac766"] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell["Example problems", "Section", CellChangeTimes->{ 3.4710857577642*^9, {3.4710858112098*^9, 3.4710858239082003`*^9}, { 3.4710859049969997`*^9, 3.4710859061046*^9}, {3.4710892678536*^9, 3.4710892701*^9}, 3.4710893817648*^9, {3.4710894793115997`*^9, 3.4710895039596*^9}, {3.4710905118050003`*^9, 3.471090512741*^9}, { 3.4725252603414*^9, 3.472525264101*^9}, {3.4725253064862003`*^9, 3.4725253370934*^9}, {3.4747629998582*^9, 3.4747630055678*^9}, { 3.4752273878873997`*^9, 3.4752273904926*^9}, {3.4753663741484003`*^9, 3.4753664116352*^9}, {3.4753773567911997`*^9, 3.475377356916*^9}, 3.4753775870472*^9, {3.4754620951394*^9, 3.4754621031266003`*^9}, { 3.4779795207265997`*^9, 3.4779795250576*^9}, 3.537621939388029*^9, 3.537622997047524*^9, {3.5376238658302155`*^9, 3.5376238672312956`*^9}, { 3.600002366135124*^9, 3.6000023692707295`*^9}, {3.6000479572325335`*^9, 3.6000479595569377`*^9}, {3.6004096645198383`*^9, 3.600409717357131*^9}, { 3.600571492488657*^9, 3.600571494735061*^9}, {3.6005936284438515`*^9, 3.6005936292394524`*^9}, {3.6006741159692173`*^9, 3.6006741341492577`*^9}, {3.6009244661154003`*^9, 3.6009244668018017`*^9}, {3.6009251420333877`*^9, 3.60092514354659*^9}, 3.600925365035779*^9, 3.6043821695544815`*^9, {3.666877582611085*^9, 3.6668775848264914`*^9}, {3.6682521885528736`*^9, 3.668252217319324*^9}, 3.668252975775611*^9, {3.6682530237788477`*^9, 3.668253025596376*^9}, { 3.7304275902874823`*^9, 3.7304275953971214`*^9}, {3.792636476950823*^9, 3.792636484957548*^9}},ExpressionUUID->"40a00b8b-dcfa-4191-b3c9-\ 72923821af9d"], Cell[CellGroupData[{ Cell["Exercise problem 1", "Subsection", CellChangeTimes->{ 3.4710857577642*^9, {3.4710858112098*^9, 3.4710858239082003`*^9}, { 3.4710859049969997`*^9, 3.4710859061046*^9}, {3.4710892678536*^9, 3.4710892701*^9}, 3.4710893817648*^9, {3.4710894793115997`*^9, 3.4710895039596*^9}, {3.4710905118050003`*^9, 3.471090512741*^9}, { 3.4725252603414*^9, 3.472525264101*^9}, {3.4725253064862003`*^9, 3.4725253370934*^9}, {3.4747629998582*^9, 3.4747630055678*^9}, { 3.4752273878873997`*^9, 3.4752273904926*^9}, {3.4753663741484003`*^9, 3.4753664116352*^9}, {3.4753773567911997`*^9, 3.475377356916*^9}, 3.4753775870472*^9, {3.4754620951394*^9, 3.4754621031266003`*^9}, { 3.4779795207265997`*^9, 3.4779795250576*^9}, {3.4781723763912*^9, 3.4781723803691998`*^9}, {3.4781776626736*^9, 3.4781776635628*^9}, { 3.4781800950726*^9, 3.4781800976154003`*^9}, 3.4785220883716*^9, { 3.4802390003438*^9, 3.4802390181902*^9}, {3.60315471527808*^9, 3.6031547167444825`*^9}, 3.603155521690296*^9, {3.6660843947173185`*^9, 3.6660843958873205`*^9}, {3.6660852732892294`*^9, 3.6660852822904453`*^9}, {3.6964934969791903`*^9, 3.696493529669236*^9}, 3.696493658969417*^9, {3.696500606350078*^9, 3.696500624720104*^9}, { 3.6965019445919514`*^9, 3.6965019665919824`*^9}, {3.696502040442086*^9, 3.6965020414520874`*^9}, {3.791952497157584*^9, 3.791952504410621*^9}, { 3.791953177295227*^9, 3.791953184330231*^9}, {3.792636513237234*^9, 3.79263653883904*^9}, 3.792640960748376*^9, 3.7926423623390675`*^9},ExpressionUUID->"058f235a-37f5-47dc-83c4-\ eccf8470e38e"], Cell[BoxData[{ RowBox[{ RowBox[{"E\[Sigma]", "=", RowBox[{ FractionBox["EE", RowBox[{"1", "-", SuperscriptBox["\[Nu]", "2"]}]], RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "\[Nu]", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"\[Nu]", ",", "1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{ RowBox[{"(", RowBox[{"1", "-", "\[Nu]"}], ")"}], "/", "2"}]}], "}"}]}], "}"}]}]}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"ex", ",", "ey"}], "}"}], "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"1", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "1"}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{ RowBox[{"{", RowBox[{"er", ",", "e\[Phi]"}], "}"}], "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Cos", "[", "\[Phi]", "]"}], ",", RowBox[{"Sin", "[", "\[Phi]", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"-", RowBox[{"Sin", "[", "\[Phi]", "]"}]}], ",", RowBox[{"Cos", "[", "\[Phi]", "]"}]}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"eexy", "=", RowBox[{"{", RowBox[{ RowBox[{"Outer", "[", RowBox[{"Times", ",", "ex", ",", "ex"}], "]"}], ",", RowBox[{"Outer", "[", RowBox[{"Times", ",", "ey", ",", "ey"}], "]"}], "\[IndentingNewLine]", ",", RowBox[{ RowBox[{"Outer", "[", RowBox[{"Times", ",", "ex", ",", "ey"}], "]"}], "+", RowBox[{"Outer", "[", RowBox[{"Times", ",", "ey", ",", "ex"}], "]"}]}]}], "}"}]}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"eer\[Phi]", "=", RowBox[{"{", RowBox[{ RowBox[{"Outer", "[", RowBox[{"Times", ",", "er", ",", "er"}], "]"}], ",", RowBox[{"Outer", "[", RowBox[{"Times", ",", "e\[Phi]", ",", "e\[Phi]"}], "]"}], "\[IndentingNewLine]", ",", RowBox[{ RowBox[{"Outer", "[", RowBox[{"Times", ",", "er", ",", "e\[Phi]"}], "]"}], "+", RowBox[{"Outer", "[", RowBox[{"Times", ",", "e\[Phi]", ",", "er"}], "]"}]}]}], "}"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"A", "=", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"Outer", "[", RowBox[{"Times", ",", RowBox[{"eexy", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"eexy", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}], "*", RowBox[{"E\[Sigma]", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}]}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"B", "=", RowBox[{"Simplify", "[", RowBox[{"Sum", "[", RowBox[{ RowBox[{ RowBox[{"Outer", "[", RowBox[{"Times", ",", RowBox[{"eer\[Phi]", "[", RowBox[{"[", "i", "]"}], "]"}], ",", RowBox[{"eer\[Phi]", "[", RowBox[{"[", "j", "]"}], "]"}]}], "]"}], "*", RowBox[{"E\[Sigma]", "[", RowBox[{"[", RowBox[{"i", ",", "j"}], "]"}], "]"}]}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "3"}], "}"}], ",", RowBox[{"{", RowBox[{"j", ",", "1", ",", "3"}], "}"}]}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]"}], "\[IndentingNewLine]", RowBox[{"Simplify", "[", RowBox[{"A", "-", "B"}], "]"}], "\[IndentingNewLine]"}], "Input", CellChangeTimes->{{3.6660967834123464`*^9, 3.666096788482355*^9}, { 3.695949560175663*^9, 3.695949567405673*^9}, {3.695949627465757*^9, 3.695949628675759*^9}, {3.6964935354592447`*^9, 3.696493668959431*^9}, 3.696501778351719*^9, {3.6965019414619474`*^9, 3.6965019419319477`*^9}, { 3.6965019731219916`*^9, 3.6965020166720524`*^9}, {3.696502051582102*^9, 3.6965020719521303`*^9}, {3.6965021140421886`*^9, 3.6965021543522453`*^9}, {3.6965772421616664`*^9, 3.6965772681417027`*^9}, {3.7925584845952363`*^9, 3.7925584911826506`*^9}, {3.7925585217498474`*^9, 3.7925585218910875`*^9}, {3.792558554075327*^9, 3.792558557784437*^9}, { 3.792558671407984*^9, 3.7925586746519556`*^9}, {3.792564245668003*^9, 3.7925642462950478`*^9}, 3.792564281534704*^9, {3.7925647475013676`*^9, 3.792564751126097*^9}, {3.7925648201521907`*^9, 3.7925648225030465`*^9}, { 3.792567063790918*^9, 3.7925671329203243`*^9}, {3.792567224777611*^9, 3.79256723192507*^9}, 3.7925674025156918`*^9, {3.7925678042836637`*^9, 3.7925678140810623`*^9}, {3.79256794434678*^9, 3.7925679480860615`*^9}, { 3.792568117048511*^9, 3.792568122685511*^9}, {3.7925682814597254`*^9, 3.792568283122607*^9}, {3.792568624525566*^9, 3.792568627034419*^9}, { 3.7925703544308853`*^9, 3.7925703814536676`*^9}, {3.792571007487669*^9, 3.7925710209880543`*^9}, {3.7925710584190063`*^9, 3.79257106562472*^9}, { 3.7925712518126144`*^9, 3.7925712519613543`*^9}, {3.792571324042125*^9, 3.792571327361492*^9}, {3.7925714525204506`*^9, 3.7925715067808056`*^9}, { 3.792571718427247*^9, 3.792571742304722*^9}, {3.7925912657573986`*^9, 3.7925912703227863`*^9}, {3.79263627957973*^9, 3.7926363311980457`*^9}, { 3.7926365498859687`*^9, 3.7926366205451093`*^9}, 3.792636657586319*^9, { 3.7926366954346275`*^9, 3.792636781020611*^9}, {3.792637252272473*^9, 3.7926372533447304`*^9}, {3.7926373604491634`*^9, 3.792637386305845*^9}, { 3.7926374451634507`*^9, 3.792637461335726*^9}, {3.7926400740533457`*^9, 3.7926401339694033`*^9}, {3.7926405567165103`*^9, 3.792640587997743*^9}, { 3.7926406659703007`*^9, 3.792640747055319*^9}, {3.792640879065442*^9, 3.7926409040375195`*^9}, {3.79264547630871*^9, 3.7926454769081364`*^9}, { 3.7926455166247683`*^9, 3.7926455334238358`*^9}, 3.792645578236719*^9, { 3.792646082252921*^9, 3.792646110693123*^9}, {3.7928126729140964`*^9, 3.792812713483203*^9}, {3.792855629731773*^9, 3.7928556354480305`*^9}, { 3.7929783299879456`*^9, 3.7929783502009363`*^9}, {3.8563178875738525`*^9, 3.856317916228176*^9}, {3.8563180064783087`*^9, 3.856318732086397*^9}, { 3.856318922639045*^9, 3.8563189660646563`*^9}, {3.8563190192659626`*^9, 3.85631903145503*^9}},ExpressionUUID->"e77e5056-5e15-4c76-8c25-\ 5a6bd2b9002b"], Cell[BoxData["\[IndentingNewLine]"], "Input", CellChangeTimes->{ 3.8563189399804244`*^9},ExpressionUUID->"f676107b-b7b3-4d05-85ac-\ 634b3fece0e9"] }, Closed]], Cell[CellGroupData[{ Cell["Exercise problem 2", "Subsection", CellChangeTimes->{ 3.4710857577642*^9, {3.4710858112098*^9, 3.4710858239082003`*^9}, { 3.4710859049969997`*^9, 3.4710859061046*^9}, {3.4710892678536*^9, 3.4710892701*^9}, 3.4710893817648*^9, {3.4710894793115997`*^9, 3.4710895039596*^9}, {3.4710905118050003`*^9, 3.471090512741*^9}, { 3.4725252603414*^9, 3.472525264101*^9}, {3.4725253064862003`*^9, 3.4725253370934*^9}, {3.4747629998582*^9, 3.4747630055678*^9}, { 3.4752273878873997`*^9, 3.4752273904926*^9}, {3.4753663741484003`*^9, 3.4753664116352*^9}, {3.4753773567911997`*^9, 3.475377356916*^9}, 3.4753775870472*^9, {3.4754620951394*^9, 3.4754621031266003`*^9}, { 3.4779795207265997`*^9, 3.4779795250576*^9}, {3.4781723763912*^9, 3.4781723803691998`*^9}, {3.4781776626736*^9, 3.4781776635628*^9}, { 3.4781800950726*^9, 3.4781800976154003`*^9}, 3.4785220883716*^9, { 3.4802390003438*^9, 3.4802390181902*^9}, {3.60315471527808*^9, 3.6031547167444825`*^9}, 3.603155521690296*^9, {3.6660843947173185`*^9, 3.6660843958873205`*^9}, {3.6660852732892294`*^9, 3.6660852822904453`*^9}, {3.6964934969791903`*^9, 3.696493529669236*^9}, 3.696493658969417*^9, {3.696500606350078*^9, 3.696500624720104*^9}, { 3.6965019445919514`*^9, 3.6965019665919824`*^9}, {3.696502040442086*^9, 3.6965020414520874`*^9}, {3.791952497157584*^9, 3.791952504410621*^9}, { 3.791953177295227*^9, 3.791953184330231*^9}, {3.792636513237234*^9, 3.79263653883904*^9}, {3.7926423723440523`*^9, 3.792642372540087*^9}, 3.7929771449783726`*^9},ExpressionUUID->"96611b31-3a2d-4d88-ad4d-\ b57525805571"], Cell[BoxData[ RowBox[{"\[IndentingNewLine]", RowBox[{"(*", RowBox[{"cylindrical", " ", "coordinate", " ", "system", " ", "def"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"xyz", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", RowBox[{"{", RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", " ", RowBox[{"r", ">", "0"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"displacement", " ", "and", " ", "volume", " ", "force"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"bb", "=", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"uu", "=", RowBox[{"{", RowBox[{ RowBox[{"u", "[", "y", "]"}], ",", RowBox[{"v", "[", "y", "]"}], ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"thin", " ", "slab", " ", "model"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"\[Sigma]\[Sigma]", "=", RowBox[{"DDOT", "[", RowBox[{ RowBox[{"EPLATE", "[", RowBox[{"\[CapitalEpsilon]", ",", "\[Nu]"}], "]"}], ",", RowBox[{ RowBox[{"GRAD", "[", "xyz", "]"}], "[", "uu", "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"eqs", "=", RowBox[{ RowBox[{ RowBox[{"DIV", "[", "xyz", "]"}], "[", "\[Sigma]\[Sigma]", "]"}], "+", "bb"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"boundary", " ", "value", " ", "problem"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"eq", "[", "1", "]"}], "=", RowBox[{"eqs", "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"eq", "[", "2", "]"}], "=", RowBox[{"eqs", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"eq", "[", "3", "]"}], "=", RowBox[{"u", "[", "0", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"eq", "[", "4", "]"}], "=", RowBox[{"v", "[", "0", "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"eq", "[", "5", "]"}], "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"\[Sigma]\[Sigma]", "[", RowBox[{"[", RowBox[{"2", ",", "1"}], "]"}], "]"}], "/.", RowBox[{"y", "\[Rule]", "h"}]}], ")"}], "-", "\[Tau]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"eq", "[", "6", "]"}], "=", RowBox[{"(", RowBox[{ RowBox[{"\[Sigma]\[Sigma]", "[", RowBox[{"[", RowBox[{"2", ",", "2"}], "]"}], "]"}], "/.", RowBox[{"y", "\[Rule]", "h"}]}], ")"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"sol", "=", RowBox[{ RowBox[{"DSolve", "[", RowBox[{ RowBox[{"Table", "[", RowBox[{ RowBox[{ RowBox[{"eq", "[", "i", "]"}], "\[Equal]", "0"}], ",", RowBox[{"{", RowBox[{"i", ",", "1", ",", "6"}], "}"}]}], "]"}], ",", RowBox[{"{", RowBox[{"u", ",", "v"}], "}"}], ",", "y"}], "]"}], "[", RowBox[{"[", "1", "]"}], "]"}]}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"Simplify", "[", RowBox[{"\[Sigma]\[Sigma]", " ", "/.", "sol"}], "]"}], "\[IndentingNewLine]"}]}]], "Input", CellChangeTimes->{{3.6660967834123464`*^9, 3.666096788482355*^9}, { 3.695949560175663*^9, 3.695949567405673*^9}, {3.695949627465757*^9, 3.695949628675759*^9}, {3.6964935354592447`*^9, 3.696493668959431*^9}, 3.696501778351719*^9, {3.6965019414619474`*^9, 3.6965019419319477`*^9}, { 3.6965019731219916`*^9, 3.6965020166720524`*^9}, {3.696502051582102*^9, 3.6965020719521303`*^9}, {3.6965021140421886`*^9, 3.6965021543522453`*^9}, {3.6965772421616664`*^9, 3.6965772681417027`*^9}, {3.7925584845952363`*^9, 3.7925584911826506`*^9}, {3.7925585217498474`*^9, 3.7925585218910875`*^9}, {3.792558554075327*^9, 3.792558557784437*^9}, { 3.792558671407984*^9, 3.7925586746519556`*^9}, {3.792564245668003*^9, 3.7925642462950478`*^9}, 3.792564281534704*^9, {3.7925647475013676`*^9, 3.792564751126097*^9}, {3.7925648201521907`*^9, 3.7925648225030465`*^9}, { 3.792567063790918*^9, 3.7925671329203243`*^9}, {3.792567224777611*^9, 3.79256723192507*^9}, 3.7925674025156918`*^9, {3.7925678042836637`*^9, 3.7925678140810623`*^9}, {3.79256794434678*^9, 3.7925679480860615`*^9}, { 3.792568117048511*^9, 3.792568122685511*^9}, {3.7925682814597254`*^9, 3.792568283122607*^9}, {3.792568624525566*^9, 3.792568627034419*^9}, { 3.7925703544308853`*^9, 3.7925703814536676`*^9}, {3.792571007487669*^9, 3.7925710209880543`*^9}, {3.7925710584190063`*^9, 3.79257106562472*^9}, { 3.7925712518126144`*^9, 3.7925712519613543`*^9}, {3.792571324042125*^9, 3.792571327361492*^9}, {3.7925714525204506`*^9, 3.7925715067808056`*^9}, { 3.792571718427247*^9, 3.792571742304722*^9}, {3.7925912657573986`*^9, 3.7925912703227863`*^9}, {3.79263627957973*^9, 3.7926363311980457`*^9}, { 3.7926365498859687`*^9, 3.7926366205451093`*^9}, 3.792636657586319*^9, { 3.7926366954346275`*^9, 3.792636781020611*^9}, {3.792637252272473*^9, 3.7926372533447304`*^9}, {3.7926373604491634`*^9, 3.792637386305845*^9}, { 3.7926374451634507`*^9, 3.792637461335726*^9}, {3.7926400740533457`*^9, 3.7926401339694033`*^9}, {3.7926405567165103`*^9, 3.792640587997743*^9}, { 3.7926406659703007`*^9, 3.792640747055319*^9}, {3.792640879065442*^9, 3.7926409040375195`*^9}, {3.7926409891805267`*^9, 3.792641237656807*^9}, { 3.792641354447501*^9, 3.7926413563803883`*^9}, {3.7926442358456583`*^9, 3.7926443740368595`*^9}, {3.7926456139386654`*^9, 3.792645616404057*^9}, { 3.792645729643696*^9, 3.792645736788617*^9}, {3.792645773768033*^9, 3.7926458162367473`*^9}, {3.792645915457271*^9, 3.7926459227492557`*^9}, { 3.792656605770455*^9, 3.792656626882081*^9}, {3.792854686676328*^9, 3.7928547156452847`*^9}, 3.792854758388832*^9, 3.7928556213082223`*^9, { 3.7929771568643274`*^9, 3.7929772763831253`*^9}, {3.7929773416939363`*^9, 3.7929774820406485`*^9}, {3.792977523908347*^9, 3.7929778639627104`*^9}, { 3.7929778961215143`*^9, 3.792977905920352*^9}, {3.7929779520035353`*^9, 3.7929781479014225`*^9}, {3.7929782105074205`*^9, 3.792978275879098*^9}, 3.79297831350446*^9, {3.8563194862995367`*^9, 3.856319492280815*^9}}, CellLabel-> "In[103]:=",ExpressionUUID->"00b16eb9-0392-48cc-81a8-2c8f0f6c7186"] }, Closed]], Cell[CellGroupData[{ Cell["Exercise problem 5", "Subsection", CellChangeTimes->{ 3.4710857577642*^9, {3.4710858112098*^9, 3.4710858239082003`*^9}, { 3.4710859049969997`*^9, 3.4710859061046*^9}, {3.4710892678536*^9, 3.4710892701*^9}, 3.4710893817648*^9, {3.4710894793115997`*^9, 3.4710895039596*^9}, {3.4710905118050003`*^9, 3.471090512741*^9}, { 3.4725252603414*^9, 3.472525264101*^9}, {3.4725253064862003`*^9, 3.4725253370934*^9}, {3.4747629998582*^9, 3.4747630055678*^9}, { 3.4752273878873997`*^9, 3.4752273904926*^9}, {3.4753663741484003`*^9, 3.4753664116352*^9}, {3.4753773567911997`*^9, 3.475377356916*^9}, 3.4753775870472*^9, {3.4754620951394*^9, 3.4754621031266003`*^9}, { 3.4779795207265997`*^9, 3.4779795250576*^9}, {3.4781723763912*^9, 3.4781723803691998`*^9}, {3.4781776626736*^9, 3.4781776635628*^9}, { 3.4781800950726*^9, 3.4781800976154003`*^9}, 3.4785220883716*^9, { 3.4802390003438*^9, 3.4802390181902*^9}, {3.60315471527808*^9, 3.6031547167444825`*^9}, 3.603155521690296*^9, {3.6660843947173185`*^9, 3.6660843958873205`*^9}, {3.6660852732892294`*^9, 3.6660852822904453`*^9}, {3.6964934969791903`*^9, 3.696493529669236*^9}, 3.696493658969417*^9, {3.696500606350078*^9, 3.696500624720104*^9}, { 3.6965019445919514`*^9, 3.6965019665919824`*^9}, {3.696502040442086*^9, 3.6965020414520874`*^9}, {3.791952497157584*^9, 3.791952504410621*^9}, { 3.791953177295227*^9, 3.791953184330231*^9}, {3.792636513237234*^9, 3.79263653883904*^9}, {3.7926423723440523`*^9, 3.792642372540087*^9}, 3.7926465436707745`*^9, 3.7926466656202493`*^9},ExpressionUUID->"1d2c0176-782d-4b0e-86cc-\ 6f9e6aaf1797"], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"cylindrical", " ", "coordinate", " ", "system", " ", "def"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"r\[Phi]z", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"r", " ", RowBox[{"Cos", "[", "\[Phi]", "]"}]}], ",", RowBox[{"r", " ", RowBox[{"Sin", "[", "\[Phi]", "]"}]}], ",", "z"}], "}"}], ",", RowBox[{"{", RowBox[{"r", ",", "\[Phi]", ",", "z"}], "}"}], ",", " ", RowBox[{"r", ">", "0"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{ RowBox[{"uu", "=", RowBox[{"{", RowBox[{ RowBox[{"u", "[", RowBox[{"r", ",", "\[Phi]"}], "]"}], ",", RowBox[{"v", "[", RowBox[{"r", ",", "\[Phi]"}], "]"}], ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{"Simplify", "[", RowBox[{"t", "*", RowBox[{"DDOT", "[", RowBox[{ RowBox[{"EPLATE", "[", RowBox[{"\[CapitalEpsilon]", ",", "\[Nu]"}], "]"}], ",", RowBox[{ RowBox[{"GRAD", "[", "r\[Phi]z", "]"}], "[", "uu", "]"}]}], "]"}]}], "]"}], "]"}], "\[IndentingNewLine]"}]}]], "Input", CellChangeTimes->{{3.6660967834123464`*^9, 3.666096788482355*^9}, { 3.695949560175663*^9, 3.695949567405673*^9}, {3.695949627465757*^9, 3.695949628675759*^9}, {3.6964935354592447`*^9, 3.696493668959431*^9}, 3.696501778351719*^9, {3.6965019414619474`*^9, 3.6965019419319477`*^9}, { 3.6965019731219916`*^9, 3.6965020166720524`*^9}, {3.696502051582102*^9, 3.6965020719521303`*^9}, {3.6965021140421886`*^9, 3.6965021543522453`*^9}, {3.6965772421616664`*^9, 3.6965772681417027`*^9}, {3.7925584845952363`*^9, 3.7925584911826506`*^9}, {3.7925585217498474`*^9, 3.7925585218910875`*^9}, {3.792558554075327*^9, 3.792558557784437*^9}, { 3.792558671407984*^9, 3.7925586746519556`*^9}, {3.792564245668003*^9, 3.7925642462950478`*^9}, 3.792564281534704*^9, {3.7925647475013676`*^9, 3.792564751126097*^9}, {3.7925648201521907`*^9, 3.7925648225030465`*^9}, { 3.792567063790918*^9, 3.7925671329203243`*^9}, {3.792567224777611*^9, 3.79256723192507*^9}, 3.7925674025156918`*^9, {3.7925678042836637`*^9, 3.7925678140810623`*^9}, {3.79256794434678*^9, 3.7925679480860615`*^9}, { 3.792568117048511*^9, 3.792568122685511*^9}, {3.7925682814597254`*^9, 3.792568283122607*^9}, {3.792568624525566*^9, 3.792568627034419*^9}, { 3.7925703544308853`*^9, 3.7925703814536676`*^9}, {3.792571007487669*^9, 3.7925710209880543`*^9}, {3.7925710584190063`*^9, 3.79257106562472*^9}, { 3.7925712518126144`*^9, 3.7925712519613543`*^9}, {3.792571324042125*^9, 3.792571327361492*^9}, {3.7925714525204506`*^9, 3.7925715067808056`*^9}, { 3.792571718427247*^9, 3.792571742304722*^9}, {3.7925912657573986`*^9, 3.7925912703227863`*^9}, {3.79263627957973*^9, 3.7926363311980457`*^9}, { 3.7926365498859687`*^9, 3.7926366205451093`*^9}, 3.792636657586319*^9, { 3.7926366954346275`*^9, 3.792636781020611*^9}, {3.792637252272473*^9, 3.7926372533447304`*^9}, {3.7926373604491634`*^9, 3.792637386305845*^9}, { 3.7926374451634507`*^9, 3.792637461335726*^9}, {3.7926400740533457`*^9, 3.7926401339694033`*^9}, {3.7926405567165103`*^9, 3.792640587997743*^9}, { 3.7926406659703007`*^9, 3.792640747055319*^9}, {3.792640879065442*^9, 3.7926409040375195`*^9}, {3.7926409891805267`*^9, 3.792641237656807*^9}, { 3.792641354447501*^9, 3.7926413563803883`*^9}, {3.7926442358456583`*^9, 3.7926443740368595`*^9}, {3.7926456139386654`*^9, 3.792645616404057*^9}, { 3.792645729643696*^9, 3.792645736788617*^9}, {3.792645773768033*^9, 3.7926458162367473`*^9}, {3.792645915457271*^9, 3.7926459227492557`*^9}, { 3.7926465462747645`*^9, 3.79264661200599*^9}, {3.792648710710764*^9, 3.7926487393932557`*^9}, {3.7926487839391785`*^9, 3.7926488421205935`*^9}, {3.792648880612414*^9, 3.7926489676569724`*^9}, { 3.7926490182167296`*^9, 3.7926490604132338`*^9}, {3.792649100534929*^9, 3.7926491155468817`*^9}, {3.7926491950873065`*^9, 3.7926493286080837`*^9}, {3.792649456834825*^9, 3.7926495394590435`*^9}, { 3.792649574061387*^9, 3.79264970713194*^9}, {3.7926499408120832`*^9, 3.7926499822071095`*^9}, {3.7926500321391106`*^9, 3.7926500341628647`*^9}, {3.7926501699380007`*^9, 3.792650184540125*^9}, { 3.792656674198734*^9, 3.7926566883533077`*^9}, {3.7926567279317493`*^9, 3.79265681422954*^9}, {3.7926568546305704`*^9, 3.7926568591991196`*^9}, { 3.7926571186526995`*^9, 3.7926571435643706`*^9}, {3.792657538909108*^9, 3.792657547639739*^9}, {3.7926575815540295`*^9, 3.7926576218596487`*^9}, { 3.7926576698131356`*^9, 3.7926576985868692`*^9}, {3.792676791307905*^9, 3.792676824346179*^9}, {3.7926768584199357`*^9, 3.7926768652744293`*^9}, { 3.7928127965970325`*^9, 3.7928128247524414`*^9}, {3.7928138970496793`*^9, 3.792813897721137*^9}, 3.792813930657861*^9, 3.792853558460622*^9, { 3.7928548524002066`*^9, 3.7928548571900845`*^9}, {3.7928548902208395`*^9, 3.7928548915953207`*^9}, {3.7928549691341887`*^9, 3.7928549800426197`*^9}, {3.792978449181644*^9, 3.792978457120002*^9}}, CellLabel->"In[76]:=",ExpressionUUID->"28aa3560-a750-4049-bc05-7d02a9fa356c"] }, Closed]], Cell[CellGroupData[{ Cell["Exercise problem 6", "Subsection", CellChangeTimes->{ 3.4710857577642*^9, {3.4710858112098*^9, 3.4710858239082003`*^9}, { 3.4710859049969997`*^9, 3.4710859061046*^9}, {3.4710892678536*^9, 3.4710892701*^9}, 3.4710893817648*^9, {3.4710894793115997`*^9, 3.4710895039596*^9}, {3.4710905118050003`*^9, 3.471090512741*^9}, { 3.4725252603414*^9, 3.472525264101*^9}, {3.4725253064862003`*^9, 3.4725253370934*^9}, {3.4747629998582*^9, 3.4747630055678*^9}, { 3.4752273878873997`*^9, 3.4752273904926*^9}, {3.4753663741484003`*^9, 3.4753664116352*^9}, {3.4753773567911997`*^9, 3.475377356916*^9}, 3.4753775870472*^9, {3.4754620951394*^9, 3.4754621031266003`*^9}, { 3.4779795207265997`*^9, 3.4779795250576*^9}, {3.4781723763912*^9, 3.4781723803691998`*^9}, {3.4781776626736*^9, 3.4781776635628*^9}, { 3.4781800950726*^9, 3.4781800976154003`*^9}, 3.4785220883716*^9, { 3.4802390003438*^9, 3.4802390181902*^9}, {3.60315471527808*^9, 3.6031547167444825`*^9}, 3.603155521690296*^9, {3.6660843947173185`*^9, 3.6660843958873205`*^9}, {3.6660852732892294`*^9, 3.6660852822904453`*^9}, {3.6964934969791903`*^9, 3.696493529669236*^9}, 3.696493658969417*^9, {3.696500606350078*^9, 3.696500624720104*^9}, { 3.6965019445919514`*^9, 3.6965019665919824`*^9}, {3.696502040442086*^9, 3.6965020414520874`*^9}, {3.791952497157584*^9, 3.791952504410621*^9}, { 3.791953177295227*^9, 3.791953184330231*^9}, {3.792636513237234*^9, 3.79263653883904*^9}, {3.7926423723440523`*^9, 3.792642372540087*^9}, 3.7926465436707745`*^9, 3.7926466656202493`*^9, 3.7928550062227163`*^9},ExpressionUUID->"6c57edae-2113-4493-9e69-\ 23824389683a"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{"cylindrical", " ", "coordinate", " ", "system", " ", "def"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"r\[Phi]z", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"r", " ", RowBox[{"Cos", "[", "\[Phi]", "]"}]}], ",", RowBox[{"r", " ", RowBox[{"Sin", "[", "\[Phi]", "]"}]}], ",", "z"}], "}"}], ",", RowBox[{"{", RowBox[{"r", ",", "\[Phi]", ",", "z"}], "}"}], ",", " ", RowBox[{"r", ">", "0"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"thin", " ", "slab", " ", "model"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"uu", "=", RowBox[{"{", RowBox[{ RowBox[{"u", "[", "r", "]"}], ",", RowBox[{"v", "[", "r", "]"}], ",", "0"}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"NN", "=", RowBox[{"Simplify", "[", RowBox[{"t", "*", RowBox[{"DDOT", "[", RowBox[{ RowBox[{"EPLATE", "[", RowBox[{"\[CapitalEpsilon]", ",", "\[Nu]"}], "]"}], ",", RowBox[{ RowBox[{"GRAD", "[", "r\[Phi]z", "]"}], "[", "uu", "]"}]}], "]"}]}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"eqs", "=", RowBox[{"Simplify", "[", RowBox[{ RowBox[{"DIV", "[", "r\[Phi]z", "]"}], "[", "NN", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"boundary", " ", "value", " ", "problem"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"eq1", "=", RowBox[{"eqs", "[", RowBox[{"[", "1", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"eq2", "=", RowBox[{"eqs", "[", RowBox[{"[", "2", "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"eq3", "=", RowBox[{ RowBox[{"(", RowBox[{ RowBox[{"NN", "[", RowBox[{"[", RowBox[{"1", ",", "2"}], "]"}], "]"}], "-", "\[Tau]"}], ")"}], "/.", RowBox[{"r", "\[Rule]", "R"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"eq4", "=", RowBox[{ RowBox[{"(", RowBox[{"NN", "[", RowBox[{"[", RowBox[{"1", ",", "1"}], "]"}], "]"}], ")"}], "/.", RowBox[{"r", "\[Rule]", "R"}]}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"eq5", "=", RowBox[{"v", "[", RowBox[{"\[Epsilon]", " ", "R"}], "]"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"eq6", "=", RowBox[{"u", "[", RowBox[{"\[Epsilon]", " ", "R"}], "]"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"DSolve", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"eq1", "\[Equal]", "0"}], ",", RowBox[{"eq2", "\[Equal]", "0"}], ",", RowBox[{"eq3", "\[Equal]", "0"}], ",", RowBox[{"eq4", "\[Equal]", "0"}], ",", RowBox[{"eq5", "\[Equal]", "0"}], ",", RowBox[{"eq6", "\[Equal]", "0"}]}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"u", "[", "r", "]"}], ",", RowBox[{"v", "[", "r", "]"}]}], "}"}], ",", "r"}], "]"}], "\[IndentingNewLine]"}]}]], "Input", CellChangeTimes->{{3.6660967834123464`*^9, 3.666096788482355*^9}, { 3.695949560175663*^9, 3.695949567405673*^9}, {3.695949627465757*^9, 3.695949628675759*^9}, {3.6964935354592447`*^9, 3.696493668959431*^9}, 3.696501778351719*^9, {3.6965019414619474`*^9, 3.6965019419319477`*^9}, { 3.6965019731219916`*^9, 3.6965020166720524`*^9}, {3.696502051582102*^9, 3.6965020719521303`*^9}, {3.6965021140421886`*^9, 3.6965021543522453`*^9}, {3.6965772421616664`*^9, 3.6965772681417027`*^9}, {3.7925584845952363`*^9, 3.7925584911826506`*^9}, {3.7925585217498474`*^9, 3.7925585218910875`*^9}, {3.792558554075327*^9, 3.792558557784437*^9}, { 3.792558671407984*^9, 3.7925586746519556`*^9}, {3.792564245668003*^9, 3.7925642462950478`*^9}, 3.792564281534704*^9, {3.7925647475013676`*^9, 3.792564751126097*^9}, {3.7925648201521907`*^9, 3.7925648225030465`*^9}, { 3.792567063790918*^9, 3.7925671329203243`*^9}, {3.792567224777611*^9, 3.79256723192507*^9}, 3.7925674025156918`*^9, {3.7925678042836637`*^9, 3.7925678140810623`*^9}, {3.79256794434678*^9, 3.7925679480860615`*^9}, { 3.792568117048511*^9, 3.792568122685511*^9}, {3.7925682814597254`*^9, 3.792568283122607*^9}, {3.792568624525566*^9, 3.792568627034419*^9}, { 3.7925703544308853`*^9, 3.7925703814536676`*^9}, {3.792571007487669*^9, 3.7925710209880543`*^9}, {3.7925710584190063`*^9, 3.79257106562472*^9}, { 3.7925712518126144`*^9, 3.7925712519613543`*^9}, {3.792571324042125*^9, 3.792571327361492*^9}, {3.7925714525204506`*^9, 3.7925715067808056`*^9}, { 3.792571718427247*^9, 3.792571742304722*^9}, {3.7925912657573986`*^9, 3.7925912703227863`*^9}, {3.79263627957973*^9, 3.7926363311980457`*^9}, { 3.7926365498859687`*^9, 3.7926366205451093`*^9}, 3.792636657586319*^9, { 3.7926366954346275`*^9, 3.792636781020611*^9}, {3.792637252272473*^9, 3.7926372533447304`*^9}, {3.7926373604491634`*^9, 3.792637386305845*^9}, { 3.7926374451634507`*^9, 3.792637461335726*^9}, {3.7926400740533457`*^9, 3.7926401339694033`*^9}, {3.7926405567165103`*^9, 3.792640587997743*^9}, { 3.7926406659703007`*^9, 3.792640747055319*^9}, {3.792640879065442*^9, 3.7926409040375195`*^9}, {3.7926409891805267`*^9, 3.792641237656807*^9}, { 3.792641354447501*^9, 3.7926413563803883`*^9}, {3.7926442358456583`*^9, 3.7926443740368595`*^9}, {3.7926456139386654`*^9, 3.792645616404057*^9}, { 3.792645729643696*^9, 3.792645736788617*^9}, {3.792645773768033*^9, 3.7926458162367473`*^9}, {3.792645915457271*^9, 3.7926459227492557`*^9}, { 3.7926465462747645`*^9, 3.79264661200599*^9}, {3.792648710710764*^9, 3.7926487393932557`*^9}, {3.7926487839391785`*^9, 3.7926488421205935`*^9}, {3.792648880612414*^9, 3.7926489676569724`*^9}, { 3.7926490182167296`*^9, 3.7926490604132338`*^9}, {3.792649100534929*^9, 3.7926491155468817`*^9}, {3.7926491950873065`*^9, 3.7926493286080837`*^9}, {3.792649456834825*^9, 3.7926495394590435`*^9}, { 3.792649574061387*^9, 3.79264970713194*^9}, {3.7926499408120832`*^9, 3.7926499822071095`*^9}, {3.7926500321391106`*^9, 3.7926500341628647`*^9}, {3.7926501699380007`*^9, 3.792650184540125*^9}, { 3.7926502271748133`*^9, 3.7926503633961535`*^9}, {3.7926508310875545`*^9, 3.7926508377005816`*^9}, {3.792650944957585*^9, 3.7926510091838975`*^9}, { 3.7926510427120757`*^9, 3.7926510600627584`*^9}, {3.7926556398687887`*^9, 3.792655851964652*^9}, {3.7926559284174895`*^9, 3.7926560034677305`*^9}, { 3.7926577637695856`*^9, 3.792657845404994*^9}, {3.792657880388002*^9, 3.7926578985820055`*^9}, {3.792659349700811*^9, 3.7926593509119062`*^9}, { 3.792676918830967*^9, 3.7926770375327706`*^9}, {3.792677115843172*^9, 3.7926771428524466`*^9}, {3.7926771824245033`*^9, 3.792677243012619*^9}, { 3.7928128469893346`*^9, 3.7928129438050303`*^9}, {3.7928130059865217`*^9, 3.792813050454248*^9}, {3.7928131170215435`*^9, 3.7928131232015967`*^9}, { 3.7928131673140187`*^9, 3.792813227358411*^9}, {3.7928132976091924`*^9, 3.7928132980153475`*^9}, {3.7928135331501904`*^9, 3.7928138782638683`*^9}, {3.7928139358555284`*^9, 3.7928139868561187`*^9}, {3.7928140638390837`*^9, 3.792814406652067*^9}, { 3.792814447287648*^9, 3.7928146356452684`*^9}, {3.7928147506552734`*^9, 3.792814754436676*^9}, {3.792814820025283*^9, 3.7928148268464622`*^9}, { 3.7928148580874405`*^9, 3.7928148741312075`*^9}, {3.7928150834579782`*^9, 3.7928150878199387`*^9}, {3.792855039108314*^9, 3.792855321260748*^9}, { 3.7928553524079556`*^9, 3.792855409872696*^9}, {3.7928554724369173`*^9, 3.7928555666095023`*^9}, {3.792978481561473*^9, 3.792978507684194*^9}}, CellLabel->"In[79]:=",ExpressionUUID->"71aadbba-f62a-42af-ba97-97231434b6e0"], Cell[BoxData[ RowBox[{"{", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"u", "[", "r", "]"}], "\[Rule]", "0"}], ",", RowBox[{ RowBox[{"v", "[", "r", "]"}], "\[Rule]", FractionBox[ RowBox[{ RowBox[{"(", RowBox[{ SuperscriptBox["r", "2"], "-", RowBox[{ SuperscriptBox["R", "2"], " ", SuperscriptBox["\[Epsilon]", "2"]}]}], ")"}], " ", RowBox[{"(", RowBox[{"1", "+", "\[Nu]"}], ")"}], " ", "\[Tau]"}], RowBox[{"r", " ", "t", " ", SuperscriptBox["\[Epsilon]", "2"], " ", "\[CapitalEpsilon]"}]]}]}], "}"}], "}"}]], "Output", CellChangeTimes->{3.856319326579201*^9}, CellLabel->"Out[89]=",ExpressionUUID->"875db314-8128-47d8-8e68-c5180c0767fd"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell["Home assignment 2", "Subsection", CellChangeTimes->{ 3.4710857577642*^9, {3.4710858112098*^9, 3.4710858239082003`*^9}, { 3.4710859049969997`*^9, 3.4710859061046*^9}, {3.4710892678536*^9, 3.4710892701*^9}, 3.4710893817648*^9, {3.4710894793115997`*^9, 3.4710895039596*^9}, {3.4710905118050003`*^9, 3.471090512741*^9}, { 3.4725252603414*^9, 3.472525264101*^9}, {3.4725253064862003`*^9, 3.4725253370934*^9}, {3.4747629998582*^9, 3.4747630055678*^9}, { 3.4752273878873997`*^9, 3.4752273904926*^9}, {3.4753663741484003`*^9, 3.4753664116352*^9}, {3.4753773567911997`*^9, 3.475377356916*^9}, 3.4753775870472*^9, {3.4754620951394*^9, 3.4754621031266003`*^9}, { 3.4779795207265997`*^9, 3.4779795250576*^9}, {3.4781723763912*^9, 3.4781723803691998`*^9}, {3.4781776626736*^9, 3.4781776635628*^9}, { 3.4781800950726*^9, 3.4781800976154003`*^9}, 3.4785220883716*^9, { 3.4802390003438*^9, 3.4802390181902*^9}, {3.60315471527808*^9, 3.6031547167444825`*^9}, 3.603155521690296*^9, {3.6660843947173185`*^9, 3.6660843958873205`*^9}, {3.6660852732892294`*^9, 3.6660852822904453`*^9}, {3.6964934969791903`*^9, 3.696493529669236*^9}, 3.696493658969417*^9, {3.696500606350078*^9, 3.696500624720104*^9}, { 3.6965019445919514`*^9, 3.6965019665919824`*^9}, {3.696502040442086*^9, 3.6965020414520874`*^9}, {3.791952497157584*^9, 3.791952504410621*^9}, { 3.791953177295227*^9, 3.791953184330231*^9}, {3.792636513237234*^9, 3.79263653883904*^9}, 3.792640960748376*^9, 3.7926423623390675`*^9, { 3.792851420801773*^9, 3.792851434697543*^9}},ExpressionUUID->"dd42abb7-371e-484b-a521-\ 795b673612a4"], Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "spherical", " ", "shell", " ", "coordinate", " ", "system", " ", "def"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"spe", "=", RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"R", "-", "n"}], ")"}], RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"Sin", "[", "\[Theta]", "]"}], " ", RowBox[{"Cos", "[", "\[Phi]", "]"}]}], ",", RowBox[{ RowBox[{"Sin", "[", "\[Theta]", "]"}], " ", RowBox[{"Sin", "[", "\[Phi]", "]"}]}], ",", RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "}"}]}], ",", RowBox[{"{", RowBox[{"\[Phi]", ",", "\[Theta]", ",", "n"}], "}"}], ",", " ", RowBox[{ RowBox[{ RowBox[{"R", "-", "n"}], ">", "0"}], "&&", " ", RowBox[{"\[Theta]", ">", "0"}], " ", "&&", " ", RowBox[{"\[Theta]", "<", "Pi"}]}]}], "}"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"stress", " ", "and", " ", "volume", " ", "force"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"NN", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"N\[Phi]\[Phi]", ",", "0", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "N\[Theta]\[Theta]", ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b", "=", RowBox[{"{", RowBox[{"0", ",", "0", ",", RowBox[{"-", "\[CapitalDelta]p"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{ RowBox[{"FullSimplify", "[", RowBox[{ RowBox[{ RowBox[{"DIV", "[", "spe", "]"}], "[", "NN", "]"}], "+", "b"}], "]"}], " ", "/.", RowBox[{"n", "\[Rule]", "0"}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.6660967834123464`*^9, 3.666096788482355*^9}, { 3.695949560175663*^9, 3.695949567405673*^9}, {3.695949627465757*^9, 3.695949628675759*^9}, {3.6964935354592447`*^9, 3.696493668959431*^9}, 3.696501778351719*^9, {3.6965019414619474`*^9, 3.6965019419319477`*^9}, { 3.6965019731219916`*^9, 3.6965020166720524`*^9}, {3.696502051582102*^9, 3.6965020719521303`*^9}, {3.6965021140421886`*^9, 3.6965021543522453`*^9}, {3.6965772421616664`*^9, 3.6965772681417027`*^9}, {3.7925584845952363`*^9, 3.7925584911826506`*^9}, {3.7925585217498474`*^9, 3.7925585218910875`*^9}, {3.792558554075327*^9, 3.792558557784437*^9}, { 3.792558671407984*^9, 3.7925586746519556`*^9}, {3.792564245668003*^9, 3.7925642462950478`*^9}, 3.792564281534704*^9, {3.7925647475013676`*^9, 3.792564751126097*^9}, {3.7925648201521907`*^9, 3.7925648225030465`*^9}, { 3.792567063790918*^9, 3.7925671329203243`*^9}, {3.792567224777611*^9, 3.79256723192507*^9}, 3.7925674025156918`*^9, {3.7925678042836637`*^9, 3.7925678140810623`*^9}, {3.79256794434678*^9, 3.7925679480860615`*^9}, { 3.792568117048511*^9, 3.792568122685511*^9}, {3.7925682814597254`*^9, 3.792568283122607*^9}, {3.792568624525566*^9, 3.792568627034419*^9}, { 3.7925703544308853`*^9, 3.7925703814536676`*^9}, {3.792571007487669*^9, 3.7925710209880543`*^9}, {3.7925710584190063`*^9, 3.79257106562472*^9}, { 3.7925712518126144`*^9, 3.7925712519613543`*^9}, {3.792571324042125*^9, 3.792571327361492*^9}, {3.7925714525204506`*^9, 3.7925715067808056`*^9}, { 3.792571718427247*^9, 3.792571742304722*^9}, {3.7925912657573986`*^9, 3.7925912703227863`*^9}, {3.79263627957973*^9, 3.7926363311980457`*^9}, { 3.7926365498859687`*^9, 3.7926366205451093`*^9}, 3.792636657586319*^9, { 3.7926366954346275`*^9, 3.792636781020611*^9}, {3.792637252272473*^9, 3.7926372533447304`*^9}, {3.7926373604491634`*^9, 3.792637386305845*^9}, { 3.7926374451634507`*^9, 3.792637461335726*^9}, {3.7926400740533457`*^9, 3.7926401339694033`*^9}, {3.7926405567165103`*^9, 3.792640587997743*^9}, { 3.7926406659703007`*^9, 3.792640747055319*^9}, {3.792640879065442*^9, 3.7926409040375195`*^9}, {3.79264547630871*^9, 3.7926454769081364`*^9}, { 3.7926455166247683`*^9, 3.7926455334238358`*^9}, 3.792645578236719*^9, { 3.792646082252921*^9, 3.792646110693123*^9}, {3.7928126729140964`*^9, 3.792812713483203*^9}, {3.792851484852304*^9, 3.7928515343430758`*^9}, { 3.7928516322342525`*^9, 3.792851653605426*^9}, {3.792851684772189*^9, 3.79285169134925*^9}, {3.792851766532589*^9, 3.7928518018168983`*^9}, { 3.7928520571796103`*^9, 3.7928521060526395`*^9}, {3.7928521361087284`*^9, 3.792852138746013*^9}, {3.7928521873947916`*^9, 3.7928521903234606`*^9}, { 3.792852247517789*^9, 3.7928522514516816`*^9}, {3.792852318646655*^9, 3.7928523228525033`*^9}, {3.7928523631618133`*^9, 3.79285239772397*^9}, { 3.7928525201503925`*^9, 3.792852547727299*^9}, 3.7928534959804125`*^9, { 3.7928537828547506`*^9, 3.7928537833100705`*^9}}, CellLabel-> "In[116]:=",ExpressionUUID->"302181a3-46d9-4c5a-8bba-029d22831c9a"] }, Closed]], Cell[CellGroupData[{ Cell["Home assignment 3", "Subsection", CellChangeTimes->{ 3.4710857577642*^9, {3.4710858112098*^9, 3.4710858239082003`*^9}, { 3.4710859049969997`*^9, 3.4710859061046*^9}, {3.4710892678536*^9, 3.4710892701*^9}, 3.4710893817648*^9, {3.4710894793115997`*^9, 3.4710895039596*^9}, {3.4710905118050003`*^9, 3.471090512741*^9}, { 3.4725252603414*^9, 3.472525264101*^9}, {3.4725253064862003`*^9, 3.4725253370934*^9}, {3.4747629998582*^9, 3.4747630055678*^9}, { 3.4752273878873997`*^9, 3.4752273904926*^9}, {3.4753663741484003`*^9, 3.4753664116352*^9}, {3.4753773567911997`*^9, 3.475377356916*^9}, 3.4753775870472*^9, {3.4754620951394*^9, 3.4754621031266003`*^9}, { 3.4779795207265997`*^9, 3.4779795250576*^9}, {3.4781723763912*^9, 3.4781723803691998`*^9}, {3.4781776626736*^9, 3.4781776635628*^9}, { 3.4781800950726*^9, 3.4781800976154003`*^9}, 3.4785220883716*^9, { 3.4802390003438*^9, 3.4802390181902*^9}, {3.60315471527808*^9, 3.6031547167444825`*^9}, 3.603155521690296*^9, {3.6660843947173185`*^9, 3.6660843958873205`*^9}, {3.6660852732892294`*^9, 3.6660852822904453`*^9}, {3.6964934969791903`*^9, 3.696493529669236*^9}, 3.696493658969417*^9, {3.696500606350078*^9, 3.696500624720104*^9}, { 3.6965019445919514`*^9, 3.6965019665919824`*^9}, {3.696502040442086*^9, 3.6965020414520874`*^9}, {3.791952497157584*^9, 3.791952504410621*^9}, { 3.791953177295227*^9, 3.791953184330231*^9}, {3.792636513237234*^9, 3.79263653883904*^9}, 3.792640960748376*^9, 3.7926423623390675`*^9, { 3.792851420801773*^9, 3.792851434697543*^9}, 3.7928526533298545`*^9},ExpressionUUID->"7ee9c706-1c89-4e5e-8fe6-\ e8a1c6a1785e"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"(*", RowBox[{ "cylindrical", " ", "shell", " ", "coordinate", " ", "system", " ", "def"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{ RowBox[{"cyl", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{ RowBox[{"(", RowBox[{"R", "-", "n"}], ")"}], " ", RowBox[{"Cos", "[", "\[Phi]", "]"}]}], ",", RowBox[{ RowBox[{"(", RowBox[{"R", "-", "n"}], ")"}], RowBox[{"Sin", "[", "\[Phi]", "]"}]}], ",", "z"}], "}"}], ",", RowBox[{"{", RowBox[{"z", ",", "\[Phi]", ",", "n"}], "}"}], ",", " ", RowBox[{ RowBox[{"R", ">", "0"}], "&&", RowBox[{ RowBox[{"R", "-", "n"}], ">", "0"}]}]}], "}"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"(*", RowBox[{"stress", " ", "and", " ", "volume", " ", "force"}], "*)"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"NN", "=", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Nzz", "[", RowBox[{"z", ",", "\[Phi]"}], "]"}], ",", RowBox[{"Nz\[Phi]", "[", RowBox[{"z", ",", "\[Phi]"}], "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"Nz\[Phi]", "[", RowBox[{"z", ",", "\[Phi]"}], "]"}], ",", RowBox[{"N\[Phi]\[Phi]", "[", RowBox[{"z", ",", "\[Phi]"}], "]"}], ",", "0"}], "}"}], ",", RowBox[{"{", RowBox[{"0", ",", "0", ",", "0"}], "}"}]}], "}"}]}], ";"}], "\[IndentingNewLine]", RowBox[{ RowBox[{"b", "=", RowBox[{"{", RowBox[{"bz", ",", "b\[Phi]", ",", "bn"}], "}"}]}], ";"}], "\[IndentingNewLine]", "\[IndentingNewLine]", RowBox[{"MatrixForm", "[", RowBox[{ RowBox[{"FullSimplify", "[", RowBox[{ RowBox[{ RowBox[{"DIV", "[", "cyl", "]"}], "[", "NN", "]"}], "+", "b"}], "]"}], " ", "/.", RowBox[{"n", "\[Rule]", "0"}]}], "]"}]}]}]], "Input", CellChangeTimes->{{3.6660967834123464`*^9, 3.666096788482355*^9}, { 3.695949560175663*^9, 3.695949567405673*^9}, {3.695949627465757*^9, 3.695949628675759*^9}, {3.6964935354592447`*^9, 3.696493668959431*^9}, 3.696501778351719*^9, {3.6965019414619474`*^9, 3.6965019419319477`*^9}, { 3.6965019731219916`*^9, 3.6965020166720524`*^9}, {3.696502051582102*^9, 3.6965020719521303`*^9}, {3.6965021140421886`*^9, 3.6965021543522453`*^9}, {3.6965772421616664`*^9, 3.6965772681417027`*^9}, {3.7925584845952363`*^9, 3.7925584911826506`*^9}, {3.7925585217498474`*^9, 3.7925585218910875`*^9}, {3.792558554075327*^9, 3.792558557784437*^9}, { 3.792558671407984*^9, 3.7925586746519556`*^9}, {3.792564245668003*^9, 3.7925642462950478`*^9}, 3.792564281534704*^9, {3.7925647475013676`*^9, 3.792564751126097*^9}, {3.7925648201521907`*^9, 3.7925648225030465`*^9}, { 3.792567063790918*^9, 3.7925671329203243`*^9}, {3.792567224777611*^9, 3.79256723192507*^9}, 3.7925674025156918`*^9, {3.7925678042836637`*^9, 3.7925678140810623`*^9}, {3.79256794434678*^9, 3.7925679480860615`*^9}, { 3.792568117048511*^9, 3.792568122685511*^9}, {3.7925682814597254`*^9, 3.792568283122607*^9}, {3.792568624525566*^9, 3.792568627034419*^9}, { 3.7925703544308853`*^9, 3.7925703814536676`*^9}, {3.792571007487669*^9, 3.7925710209880543`*^9}, {3.7925710584190063`*^9, 3.79257106562472*^9}, { 3.7925712518126144`*^9, 3.7925712519613543`*^9}, {3.792571324042125*^9, 3.792571327361492*^9}, {3.7925714525204506`*^9, 3.7925715067808056`*^9}, { 3.792571718427247*^9, 3.792571742304722*^9}, {3.7925912657573986`*^9, 3.7925912703227863`*^9}, {3.79263627957973*^9, 3.7926363311980457`*^9}, { 3.7926365498859687`*^9, 3.7926366205451093`*^9}, 3.792636657586319*^9, { 3.7926366954346275`*^9, 3.792636781020611*^9}, {3.792637252272473*^9, 3.7926372533447304`*^9}, {3.7926373604491634`*^9, 3.792637386305845*^9}, { 3.7926374451634507`*^9, 3.792637461335726*^9}, {3.7926400740533457`*^9, 3.7926401339694033`*^9}, {3.7926405567165103`*^9, 3.792640587997743*^9}, { 3.7926406659703007`*^9, 3.792640747055319*^9}, {3.792640879065442*^9, 3.7926409040375195`*^9}, {3.79264547630871*^9, 3.7926454769081364`*^9}, { 3.7926455166247683`*^9, 3.7926455334238358`*^9}, 3.792645578236719*^9, { 3.792646082252921*^9, 3.792646110693123*^9}, {3.7928126729140964`*^9, 3.792812713483203*^9}, {3.792851484852304*^9, 3.7928515343430758`*^9}, { 3.7928516322342525`*^9, 3.792851653605426*^9}, {3.792851684772189*^9, 3.79285169134925*^9}, {3.792851766532589*^9, 3.7928518018168983`*^9}, { 3.7928520571796103`*^9, 3.7928521060526395`*^9}, {3.7928521361087284`*^9, 3.792852138746013*^9}, {3.7928521873947916`*^9, 3.7928521903234606`*^9}, { 3.792852247517789*^9, 3.7928522514516816`*^9}, {3.792852318646655*^9, 3.7928523228525033`*^9}, {3.7928523631618133`*^9, 3.79285239772397*^9}, { 3.7928525201503925`*^9, 3.792852547727299*^9}, {3.792852703673498*^9, 3.7928527069711623`*^9}, {3.792852818435492*^9, 3.7928528779791985`*^9}, { 3.7928529379096475`*^9, 3.792852994298066*^9}, {3.7928530946194897`*^9, 3.792853169663046*^9}, {3.7928533129828734`*^9, 3.7928533407625628`*^9}, { 3.792853424532394*^9, 3.7928534252051926`*^9}, 3.7928535024469404`*^9}, CellLabel-> "In[120]:=",ExpressionUUID->"dd9669d4-eac6-41a5-8386-e8be94bc2d62"], Cell[BoxData[ TagBox[ RowBox[{"(", "\[NoBreak]", TagBox[GridBox[{ { RowBox[{"bz", "+", FractionBox[ RowBox[{ SuperscriptBox["Nz\[Phi]", TagBox[ RowBox[{"(", RowBox[{"0", ",", "1"}], ")"}], Derivative], MultilineFunction->None], "[", RowBox[{"z", ",", "\[Phi]"}], "]"}], "R"], "+", RowBox[{ SuperscriptBox["Nzz", TagBox[ RowBox[{"(", RowBox[{"1", ",", "0"}], ")"}], Derivative], MultilineFunction->None], "[", RowBox[{"z", ",", "\[Phi]"}], "]"}]}]}, { RowBox[{"b\[Phi]", "+", FractionBox[ RowBox[{ SuperscriptBox["N\[Phi]\[Phi]", TagBox[ RowBox[{"(", RowBox[{"0", ",", "1"}], ")"}], Derivative], MultilineFunction->None], "[", RowBox[{"z", ",", "\[Phi]"}], "]"}], "R"], "+", RowBox[{ SuperscriptBox["Nz\[Phi]", TagBox[ RowBox[{"(", RowBox[{"1", ",", "0"}], ")"}], Derivative], MultilineFunction->None], "[", RowBox[{"z", ",", "\[Phi]"}], "]"}]}]}, { RowBox[{"bn", "+", FractionBox[ RowBox[{"N\[Phi]\[Phi]", "[", RowBox[{"z", ",", "\[Phi]"}], "]"}], "R"]}]} }, GridBoxAlignment->{"Columns" -> {{Center}}, "Rows" -> {{Baseline}}}, GridBoxSpacings->{"Columns" -> { Offset[0.27999999999999997`], { Offset[0.5599999999999999]}, Offset[0.27999999999999997`]}, "Rows" -> { Offset[0.2], { Offset[0.4]}, Offset[0.2]}}], Column], "\[NoBreak]", ")"}], Function[BoxForm`e$, MatrixForm[BoxForm`e$]]]], "Output", CellChangeTimes->{3.8563198122445683`*^9}, CellLabel-> "Out[123]//MatrixForm=",ExpressionUUID->"9e6a9c2e-d400-40ed-914f-\ 11370466d016"] }, Open ]] }, Closed]] }, Open ]] }, AutoGeneratedPackage->Automatic, WindowSize->{791., 463.75}, WindowMargins->{{70.5, Automatic}, {Automatic, 121}}, PrivateNotebookOptions->{"VersionedStylesheet"->{"Default.nb"[8.] -> False}}, Magnification:>1.5 Inherited, FrontEndVersion->"12.3 for Microsoft Windows (64-bit) (July 9, 2021)", StyleDefinitions->"Default.nb", ExpressionUUID->"6381a8e6-3515-4192-acd6-2147746b0a6a" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[567, 22, 639, 9, 100, "Section",ExpressionUUID->"f80f1a14-5e19-4d5b-b3db-8ffc5f31ec0d"], Cell[CellGroupData[{ Cell[1231, 35, 1842, 27, 81, "Subsection",ExpressionUUID->"16ef144c-4d13-402e-96b4-c4e44eec02f8"], Cell[3076, 64, 17119, 433, 1437, "Input",ExpressionUUID->"718a3088-f79e-4e0f-855c-1eb4a3539f04", InitializationCell->True] }, Closed]], Cell[CellGroupData[{ Cell[20232, 502, 1879, 27, 55, "Subsection",ExpressionUUID->"a7d5c6f1-3a2f-43a7-a35a-4129cabb329b"], Cell[22114, 531, 17855, 398, 1380, "Input",ExpressionUUID->"58c30491-20f6-423f-9696-44262d8befe7", InitializationCell->True], Cell[39972, 931, 21194, 469, 1242, "Input",ExpressionUUID->"ab668980-2ce3-4ad7-95ec-4e45a51ac766", InitializationCell->True] }, Closed]] }, Closed]], Cell[CellGroupData[{ Cell[61215, 1406, 1630, 24, 77, "Section",ExpressionUUID->"40a00b8b-dcfa-4191-b3c9-72923821af9d"], Cell[CellGroupData[{ Cell[62870, 1434, 1606, 24, 81, "Subsection",ExpressionUUID->"058f235a-37f5-47dc-83c4-eccf8470e38e"], Cell[64479, 1460, 6682, 157, 584, "Input",ExpressionUUID->"e77e5056-5e15-4c76-8c25-5a6bd2b9002b"], Cell[71164, 1619, 149, 3, 71, "Input",ExpressionUUID->"f676107b-b7b3-4d05-85ac-634b3fece0e9"] }, Closed]], Cell[CellGroupData[{ Cell[71350, 1627, 1632, 24, 55, "Subsection",ExpressionUUID->"96611b31-3a2d-4d88-ad4d-b57525805571"], Cell[72985, 1653, 6792, 151, 699, "Input",ExpressionUUID->"00b16eb9-0392-48cc-81a8-2c8f0f6c7186"] }, Closed]], Cell[CellGroupData[{ Cell[79814, 1809, 1660, 25, 55, "Subsection",ExpressionUUID->"1d2c0176-782d-4b0e-86cc-6f9e6aaf1797"], Cell[81477, 1836, 5399, 94, 214, "Input",ExpressionUUID->"28aa3560-a750-4049-bc05-7d02a9fa356c"] }, Closed]], Cell[CellGroupData[{ Cell[86913, 1935, 1684, 25, 55, "Subsection",ExpressionUUID->"6c57edae-2113-4493-9e69-23824389683a"], Cell[CellGroupData[{ Cell[88622, 1964, 8081, 163, 557, "Input",ExpressionUUID->"71aadbba-f62a-42af-ba97-97231434b6e0"], Cell[96706, 2129, 765, 22, 119, "Output",ExpressionUUID->"875db314-8128-47d8-8e68-c5180c0767fd"] }, Open ]] }, Closed]], Cell[CellGroupData[{ Cell[97520, 2157, 1655, 25, 55, "Subsection",ExpressionUUID->"dd42abb7-371e-484b-a521-795b673612a4"], Cell[99178, 2184, 5152, 102, 271, "Input",ExpressionUUID->"302181a3-46d9-4c5a-8bba-029d22831c9a"] }, Closed]], Cell[CellGroupData[{ Cell[104367, 2291, 1679, 25, 55, "Subsection",ExpressionUUID->"7ee9c706-1c89-4e5e-8fe6-e8a1c6a1785e"], Cell[CellGroupData[{ Cell[106071, 2320, 5480, 109, 300, "Input",ExpressionUUID->"dd9669d4-eac6-41a5-8386-e8be94bc2d62"], Cell[111554, 2431, 1964, 62, 158, "Output",ExpressionUUID->"9e6a9c2e-d400-40ed-914f-11370466d016"] }, Open ]] }, Closed]] }, Open ]] } ] *)