neuper@37906
|
1 |
(* chapter 'Biegelinie' from the textbook:
|
neuper@37906
|
2 |
Timischl, Kaiser. Ingenieur-Mathematik 3. Wien 1999. p.268-271.
|
neuper@37999
|
3 |
author: Walther Neuper 050826,
|
neuper@37906
|
4 |
(c) due to copyright terms
|
neuper@37906
|
5 |
*)
|
neuper@37906
|
6 |
|
bonzai@41919
|
7 |
theory Biegelinie imports Integrate Equation EqSystem Atools begin
|
neuper@37906
|
8 |
|
neuper@37906
|
9 |
consts
|
neuper@37906
|
10 |
|
neuper@42385
|
11 |
qq :: "real => real" (* Streckenlast *)
|
neuper@37999
|
12 |
Q :: "real => real" (* Querkraft *)
|
neuper@37999
|
13 |
Q' :: "real => real" (* Ableitung der Querkraft *)
|
neuper@37999
|
14 |
M'_b :: "real => real" ("M'_b") (* Biegemoment *)
|
neuper@37999
|
15 |
M'_b' :: "real => real" ("M'_b'") (* Ableitung des Biegemoments *)
|
neuper@37999
|
16 |
y'' :: "real => real" (* 2.Ableitung der Biegeline *)
|
neuper@37999
|
17 |
y' :: "real => real" (* Neigung der Biegeline *)
|
neuper@37999
|
18 |
(*y :: "real => real" (* Biegeline *)*)
|
neuper@37999
|
19 |
EI :: "real" (* Biegesteifigkeit *)
|
neuper@37906
|
20 |
|
neuper@37906
|
21 |
(*new Descriptions in the related problems*)
|
neuper@37999
|
22 |
Traegerlaenge :: "real => una"
|
neuper@37999
|
23 |
Streckenlast :: "real => una"
|
neuper@37999
|
24 |
BiegemomentVerlauf :: "bool => una"
|
neuper@37999
|
25 |
Biegelinie :: "(real => real) => una"
|
neuper@37999
|
26 |
Randbedingungen :: "bool list => una"
|
neuper@37999
|
27 |
RandbedingungenBiegung :: "bool list => una"
|
neuper@37999
|
28 |
RandbedingungenNeigung :: "bool list => una"
|
neuper@37999
|
29 |
RandbedingungenMoment :: "bool list => una"
|
neuper@37999
|
30 |
RandbedingungenQuerkraft :: "bool list => una"
|
neuper@37999
|
31 |
FunktionsVariable :: "real => una"
|
neuper@37999
|
32 |
Funktionen :: "bool list => una"
|
neuper@37999
|
33 |
Gleichungen :: "bool list => una"
|
neuper@37906
|
34 |
|
neuper@37906
|
35 |
(*Script-names*)
|
neuper@37906
|
36 |
Biegelinie2Script :: "[real,real,real,real=>real,bool list,
|
neuper@37906
|
37 |
bool] => bool"
|
neuper@37906
|
38 |
("((Script Biegelinie2Script (_ _ _ _ _ =))// (_))" 9)
|
neuper@37906
|
39 |
BiegelinieScript :: "[real,real,real,real=>real,bool list,bool list,
|
neuper@37906
|
40 |
bool] => bool"
|
neuper@37906
|
41 |
("((Script BiegelinieScript (_ _ _ _ _ _ =))// (_))" 9)
|
neuper@37906
|
42 |
Biege2xIntegrierenScript :: "[real,real,real,bool,real=>real,bool list,
|
neuper@37906
|
43 |
bool] => bool"
|
neuper@37906
|
44 |
("((Script Biege2xIntegrierenScript (_ _ _ _ _ _ =))// (_))" 9)
|
neuper@37906
|
45 |
Biege4x4SystemScript :: "[real,real,real,real=>real,bool list,
|
neuper@37906
|
46 |
bool] => bool"
|
neuper@37906
|
47 |
("((Script Biege4x4SystemScript (_ _ _ _ _ =))// (_))" 9)
|
neuper@37906
|
48 |
Biege1xIntegrierenScript ::
|
neuper@37906
|
49 |
"[real,real,real,real=>real,bool list,bool list,bool list,
|
neuper@37906
|
50 |
bool] => bool"
|
neuper@37906
|
51 |
("((Script Biege1xIntegrierenScript (_ _ _ _ _ _ _ =))// (_))" 9)
|
neuper@37906
|
52 |
Belastung2BiegelScript :: "[real,real,
|
neuper@37906
|
53 |
bool list] => bool list"
|
neuper@37906
|
54 |
("((Script Belastung2BiegelScript (_ _ =))// (_))" 9)
|
neuper@37906
|
55 |
SetzeRandbedScript :: "[bool list,bool list,
|
neuper@37906
|
56 |
bool list] => bool list"
|
neuper@37906
|
57 |
("((Script SetzeRandbedScript (_ _ =))// (_))" 9)
|
neuper@37906
|
58 |
|
neuper@52148
|
59 |
axiomatization where
|
neuper@37906
|
60 |
|
neuper@52148
|
61 |
Querkraft_Belastung: "Q' x = -qq x" and
|
neuper@52148
|
62 |
Belastung_Querkraft: "-qq x = Q' x" and
|
neuper@37906
|
63 |
|
neuper@52148
|
64 |
Moment_Querkraft: "M_b' x = Q x" and
|
neuper@52148
|
65 |
Querkraft_Moment: "Q x = M_b' x" and
|
neuper@37906
|
66 |
|
neuper@52148
|
67 |
Neigung_Moment: "y'' x = -M_b x/ EI" and
|
neuper@52148
|
68 |
Moment_Neigung: "M_b x = -EI * y'' x" and
|
neuper@37906
|
69 |
|
neuper@37906
|
70 |
(*according to rls 'simplify_Integral': .. = 1/a * .. instead .. = ../ a*)
|
neuper@37983
|
71 |
make_fun_explicit: "Not (a =!= 0) ==> (a * (f x) = b) = (f x = 1/a * b)"
|
neuper@37906
|
72 |
|
neuper@37954
|
73 |
ML {*
|
neuper@37972
|
74 |
val thy = @{theory};
|
neuper@37972
|
75 |
|
neuper@37954
|
76 |
(** theory elements for transfer into html **)
|
neuper@55420
|
77 |
store_thes
|
neuper@55420
|
78 |
[make_thy @{theory}
|
neuper@55420
|
79 |
["Walther Neuper 2005 supported by a grant from NMI Austria"],
|
neuper@55420
|
80 |
make_isa @{theory} ("IsacKnowledge", "Theorems")
|
neuper@55420
|
81 |
["Walther Neuper 2005 supported by a grant from NMI Austria"],
|
neuper@55420
|
82 |
make_thm thy "IsacKnowledge" ("Belastung_Querkraft", prop_of @{thm Belastung_Querkraft})
|
neuper@55420
|
83 |
["Walther Neuper 2005 supported by a grant from NMI Austria"],
|
neuper@55420
|
84 |
make_thm thy "IsacKnowledge" ("Moment_Neigung", prop_of @{thm Moment_Neigung})
|
neuper@55420
|
85 |
["Walther Neuper 2005 supported by a grant from NMI Austria"],
|
neuper@55420
|
86 |
make_thm thy "IsacKnowledge" ("Moment_Querkraft", prop_of @{thm Moment_Querkraft})
|
neuper@55420
|
87 |
["Walther Neuper 2005 supported by a grant from NMI Austria"],
|
neuper@55420
|
88 |
make_thm thy "IsacKnowledge" ("Neigung_Moment", prop_of @{thm Neigung_Moment})
|
neuper@55420
|
89 |
["Walther Neuper 2005 supported by a grant from NMI Austria"],
|
neuper@55420
|
90 |
make_thm thy "IsacKnowledge" ("Querkraft_Belastung", prop_of @{thm Querkraft_Belastung})
|
neuper@55420
|
91 |
["Walther Neuper 2005 supported by a grant from NMI Austria"],
|
neuper@55420
|
92 |
make_thm thy "IsacKnowledge" ("Querkraft_Moment", prop_of @{thm Querkraft_Moment})
|
neuper@55420
|
93 |
["Walther Neuper 2005 supported by a grant from NMI Austria"],
|
neuper@55420
|
94 |
make_thm thy "IsacKnowledge" ("make_fun_explicit", prop_of @{thm make_fun_explicit})
|
neuper@55420
|
95 |
["Walther Neuper 2005 supported by a grant from NMI Austria"]
|
neuper@55420
|
96 |
]
|
neuper@55420
|
97 |
*}
|
neuper@37954
|
98 |
|
neuper@37954
|
99 |
(** problems **)
|
s1210629013@55359
|
100 |
setup {* KEStore_Elems.add_pbts
|
s1210629013@55339
|
101 |
[(prep_pbt thy "pbl_bieg" [] e_pblID
|
s1210629013@55339
|
102 |
(["Biegelinien"],
|
s1210629013@55339
|
103 |
[("#Given" ,["Traegerlaenge l_l", "Streckenlast q_q"]),
|
s1210629013@55339
|
104 |
(*("#Where",["0 < l_l"]), ...wait for < and handling Arbfix*)
|
s1210629013@55339
|
105 |
("#Find" ,["Biegelinie b_b"]),
|
s1210629013@55339
|
106 |
("#Relate",["Randbedingungen r_b"])],
|
s1210629013@55339
|
107 |
append_rls "e_rls" e_rls [], NONE, [["IntegrierenUndKonstanteBestimmen2"]])),
|
s1210629013@55339
|
108 |
(prep_pbt thy "pbl_bieg_mom" [] e_pblID
|
s1210629013@55339
|
109 |
(["MomentBestimmte","Biegelinien"],
|
s1210629013@55339
|
110 |
[("#Given" ,["Traegerlaenge l_l", "Streckenlast q_q"]),
|
s1210629013@55339
|
111 |
(*("#Where",["0 < l_l"]), ...wait for < and handling Arbfix*)
|
s1210629013@55339
|
112 |
("#Find" ,["Biegelinie b_b"]),
|
s1210629013@55339
|
113 |
("#Relate",["RandbedingungenBiegung r_b","RandbedingungenMoment r_m"])
|
s1210629013@55339
|
114 |
],
|
s1210629013@55339
|
115 |
append_rls "e_rls" e_rls [], NONE, [["IntegrierenUndKonstanteBestimmen"]])),
|
s1210629013@55339
|
116 |
(prep_pbt thy "pbl_bieg_momg" [] e_pblID
|
s1210629013@55339
|
117 |
(["MomentGegebene","Biegelinien"], [], append_rls "e_rls" e_rls [], NONE,
|
s1210629013@55339
|
118 |
[["IntegrierenUndKonstanteBestimmen","2xIntegrieren"]])),
|
s1210629013@55339
|
119 |
(prep_pbt thy "pbl_bieg_einf" [] e_pblID
|
s1210629013@55339
|
120 |
(["einfache","Biegelinien"], [], append_rls "e_rls" e_rls [], NONE,
|
s1210629013@55339
|
121 |
[["IntegrierenUndKonstanteBestimmen","4x4System"]])),
|
s1210629013@55339
|
122 |
(prep_pbt thy "pbl_bieg_momquer" [] e_pblID
|
s1210629013@55339
|
123 |
(["QuerkraftUndMomentBestimmte","Biegelinien"], [], append_rls "e_rls" e_rls [], NONE,
|
s1210629013@55339
|
124 |
[["IntegrierenUndKonstanteBestimmen","1xIntegrieren"]])),
|
s1210629013@55339
|
125 |
(prep_pbt thy "pbl_bieg_vonq" [] e_pblID
|
s1210629013@55339
|
126 |
(["vonBelastungZu","Biegelinien"],
|
s1210629013@55339
|
127 |
[("#Given" ,["Streckenlast q_q","FunktionsVariable v_v"]),
|
s1210629013@55339
|
128 |
("#Find" ,["Funktionen funs'''"])],
|
s1210629013@55339
|
129 |
append_rls "e_rls" e_rls [], NONE, [["Biegelinien","ausBelastung"]])),
|
s1210629013@55339
|
130 |
(prep_pbt thy "pbl_bieg_randbed" [] e_pblID
|
s1210629013@55339
|
131 |
(["setzeRandbedingungen","Biegelinien"],
|
s1210629013@55339
|
132 |
[("#Given" ,["Funktionen fun_s","Randbedingungen r_b"]),
|
s1210629013@55339
|
133 |
("#Find" ,["Gleichungen equs'''"])],
|
s1210629013@55339
|
134 |
append_rls "e_rls" e_rls [], NONE, [["Biegelinien","setzeRandbedingungenEin"]])),
|
s1210629013@55339
|
135 |
(prep_pbt thy "pbl_equ_fromfun" [] e_pblID
|
s1210629013@55339
|
136 |
(["makeFunctionTo","equation"],
|
s1210629013@55339
|
137 |
[("#Given" ,["functionEq fu_n","substitution su_b"]),
|
s1210629013@55339
|
138 |
("#Find" ,["equality equ'''"])],
|
s1210629013@55339
|
139 |
append_rls "e_rls" e_rls [], NONE, [["Equation","fromFunction"]]))] *}
|
neuper@37999
|
140 |
ML {*
|
neuper@37954
|
141 |
(** methods **)
|
neuper@37954
|
142 |
|
neuper@37954
|
143 |
val srls = Rls {id="srls_IntegrierenUnd..",
|
neuper@37954
|
144 |
preconds = [],
|
neuper@37954
|
145 |
rew_ord = ("termlessI",termlessI),
|
neuper@37954
|
146 |
erls = append_rls "erls_in_srls_IntegrierenUnd.." e_rls
|
neuper@37999
|
147 |
[(*for asm in NTH_CONS ...*)
|
neuper@38045
|
148 |
Calc ("Orderings.ord_class.less",eval_equ "#less_"),
|
neuper@37999
|
149 |
(*2nd NTH_CONS pushes n+-1 into asms*)
|
neuper@38014
|
150 |
Calc("Groups.plus_class.plus", eval_binop "#add_")
|
neuper@37954
|
151 |
],
|
neuper@42451
|
152 |
srls = Erls, calc = [], errpatts = [],
|
neuper@37999
|
153 |
rules = [Thm ("NTH_CONS",num_str @{thm NTH_CONS}),
|
neuper@38014
|
154 |
Calc("Groups.plus_class.plus", eval_binop "#add_"),
|
neuper@37999
|
155 |
Thm ("NTH_NIL",num_str @{thm NTH_NIL}),
|
neuper@37954
|
156 |
Calc("Tools.lhs", eval_lhs"eval_lhs_"),
|
neuper@37954
|
157 |
Calc("Tools.rhs", eval_rhs"eval_rhs_"),
|
neuper@37954
|
158 |
Calc("Atools.argument'_in",
|
neuper@37954
|
159 |
eval_argument_in "Atools.argument'_in")
|
neuper@37954
|
160 |
],
|
neuper@37954
|
161 |
scr = EmptyScr};
|
neuper@37954
|
162 |
|
neuper@37954
|
163 |
val srls2 =
|
neuper@37954
|
164 |
Rls {id="srls_IntegrierenUnd..",
|
neuper@37954
|
165 |
preconds = [],
|
neuper@37954
|
166 |
rew_ord = ("termlessI",termlessI),
|
neuper@37954
|
167 |
erls = append_rls "erls_in_srls_IntegrierenUnd.." e_rls
|
neuper@37999
|
168 |
[(*for asm in NTH_CONS ...*)
|
neuper@38045
|
169 |
Calc ("Orderings.ord_class.less",eval_equ "#less_"),
|
neuper@37999
|
170 |
(*2nd NTH_CONS pushes n+-1 into asms*)
|
neuper@38014
|
171 |
Calc("Groups.plus_class.plus", eval_binop "#add_")
|
neuper@37954
|
172 |
],
|
neuper@42451
|
173 |
srls = Erls, calc = [], errpatts = [],
|
neuper@37999
|
174 |
rules = [Thm ("NTH_CONS",num_str @{thm NTH_CONS}),
|
neuper@38014
|
175 |
Calc("Groups.plus_class.plus", eval_binop "#add_"),
|
neuper@37999
|
176 |
Thm ("NTH_NIL", num_str @{thm NTH_NIL}),
|
neuper@37954
|
177 |
Calc("Tools.lhs", eval_lhs "eval_lhs_"),
|
neuper@37954
|
178 |
Calc("Atools.filter'_sameFunId",
|
neuper@37954
|
179 |
eval_filter_sameFunId "Atools.filter'_sameFunId"),
|
neuper@37954
|
180 |
(*WN070514 just for smltest/../biegelinie.sml ...*)
|
neuper@37954
|
181 |
Calc("Atools.sameFunId", eval_sameFunId "Atools.sameFunId"),
|
neuper@37969
|
182 |
Thm ("filter_Cons", num_str @{thm filter_Cons}),
|
neuper@37969
|
183 |
Thm ("filter_Nil", num_str @{thm filter_Nil}),
|
neuper@37969
|
184 |
Thm ("if_True", num_str @{thm if_True}),
|
neuper@37969
|
185 |
Thm ("if_False", num_str @{thm if_False}),
|
neuper@37969
|
186 |
Thm ("hd_thm", num_str @{thm hd_thm})
|
neuper@37954
|
187 |
],
|
neuper@37954
|
188 |
scr = EmptyScr};
|
neuper@37999
|
189 |
*}
|
neuper@37999
|
190 |
|
s1210629013@55373
|
191 |
setup {* KEStore_Elems.add_mets
|
s1210629013@55373
|
192 |
[prep_met thy "met_biege" [] e_metID
|
s1210629013@55373
|
193 |
(["IntegrierenUndKonstanteBestimmen"],
|
s1210629013@55373
|
194 |
[("#Given" ,["Traegerlaenge l_l", "Streckenlast q__q", "FunktionsVariable v_v"]),
|
s1210629013@55373
|
195 |
(*("#Where",["0 < l_l"]), ...wait for < and handling Arbfix*)
|
s1210629013@55373
|
196 |
("#Find" ,["Biegelinie b_b"]),
|
s1210629013@55373
|
197 |
("#Relate",["RandbedingungenBiegung r_b", "RandbedingungenMoment r_m"])],
|
s1210629013@55373
|
198 |
{rew_ord'="tless_true",
|
s1210629013@55373
|
199 |
rls' = append_rls "erls_IntegrierenUndK.." e_rls
|
s1210629013@55373
|
200 |
[Calc ("Atools.ident",eval_ident "#ident_"),
|
s1210629013@55373
|
201 |
Thm ("not_true",num_str @{thm not_true}),
|
s1210629013@55373
|
202 |
Thm ("not_false",num_str @{thm not_false})],
|
s1210629013@55373
|
203 |
calc = [], srls = srls, prls = Erls, crls = Atools_erls, errpats = [], nrls = Erls},
|
s1210629013@55373
|
204 |
"Script BiegelinieScript " ^
|
s1210629013@55373
|
205 |
"(l_l::real) (q__q::real) (v_v::real) (b_b::real=>real) " ^
|
s1210629013@55373
|
206 |
"(r_b::bool list) (r_m::bool list) = " ^
|
s1210629013@55373
|
207 |
" (let q___q = Take (qq v_v = q__q); " ^
|
s1210629013@55373
|
208 |
" q___q = ((Rewrite sym_neg_equal_iff_equal True) @@ " ^
|
s1210629013@55373
|
209 |
" (Rewrite Belastung_Querkraft True)) q___q; " ^
|
s1210629013@55373
|
210 |
" (Q__Q:: bool) = " ^
|
s1210629013@55373
|
211 |
" (SubProblem (Biegelinie',[named,integrate,function], " ^
|
s1210629013@55373
|
212 |
" [diff,integration,named]) " ^
|
s1210629013@55373
|
213 |
" [REAL (rhs q___q), REAL v_v, REAL_REAL Q]); " ^
|
s1210629013@55373
|
214 |
" Q__Q = Rewrite Querkraft_Moment True Q__Q; " ^
|
s1210629013@55373
|
215 |
" (M__M::bool) = " ^
|
s1210629013@55373
|
216 |
" (SubProblem (Biegelinie',[named,integrate,function], " ^
|
s1210629013@55373
|
217 |
" [diff,integration,named]) " ^
|
s1210629013@55373
|
218 |
" [REAL (rhs Q__Q), REAL v_v, REAL_REAL M_b]); " ^
|
s1210629013@55373
|
219 |
(*([5], Res), M_b x = c_2 + c * x + -1 * q_0 / 2 * x ^^^ 2*)
|
s1210629013@55373
|
220 |
" e__1 = NTH 1 r_m; " ^
|
s1210629013@55373
|
221 |
" (x__1::real) = argument_in (lhs e__1); " ^
|
s1210629013@55373
|
222 |
" (M__1::bool) = (Substitute [v_v = x__1]) M__M; " ^
|
s1210629013@55373
|
223 |
(*([6], Res), M_b 0 = c_2 + c * 0 + -1 * q_0 / 2 * 0 ^^^ 2*)
|
s1210629013@55373
|
224 |
" M__1 = (Substitute [e__1]) M__1; " ^
|
s1210629013@55373
|
225 |
(*([7], Res), 0 = c_2 + c * 0 + -1 * q_0 / 2 * 0 ^^^ 2*)
|
s1210629013@55373
|
226 |
" M__2 = Take M__M; " ^
|
s1210629013@55373
|
227 |
(*([8], Frm), M_b x = c_2 + c * x + -1 * q_0 / 2 * x ^^^ 2*)
|
s1210629013@55373
|
228 |
(*without above Take 'Substitute [v_v = x__2]' takes _last formula from ctree_*)
|
s1210629013@55373
|
229 |
" e__2 = NTH 2 r_m; " ^
|
s1210629013@55373
|
230 |
" (x__2::real) = argument_in (lhs e__2); " ^
|
s1210629013@55373
|
231 |
" (M__2::bool) = (Substitute [v_v = x__2]) M__M; " ^
|
s1210629013@55373
|
232 |
(*([8], Res), M_b L = c_2 + c * L + -1 * q_0 / 2 * L ^^^ 2*)
|
s1210629013@55373
|
233 |
" M__2 = (Substitute [e__2]) M__2; " ^
|
s1210629013@55373
|
234 |
" (c_1_2::bool list) = " ^
|
s1210629013@55373
|
235 |
" (SubProblem (Biegelinie',[LINEAR,system],[no_met]) " ^
|
s1210629013@55373
|
236 |
" [BOOL_LIST [M__1, M__2], REAL_LIST [c,c_2]]); " ^
|
s1210629013@55373
|
237 |
" M__M = Take M__M; " ^
|
s1210629013@55373
|
238 |
" M__M = ((Substitute c_1_2) @@ " ^
|
s1210629013@55373
|
239 |
" (Try (Rewrite_Set_Inst [(bdv_1, c),(bdv_2, c_2)] " ^
|
s1210629013@55373
|
240 |
" simplify_System False)) @@ " ^
|
s1210629013@55373
|
241 |
" (Rewrite Moment_Neigung False) @@ " ^
|
s1210629013@55373
|
242 |
" (Rewrite make_fun_explicit False)) M__M; " ^
|
s1210629013@55373
|
243 |
(*----------------------- and the same once more ------------------------*)
|
s1210629013@55373
|
244 |
" (N__N:: bool) = " ^
|
s1210629013@55373
|
245 |
" (SubProblem (Biegelinie',[named,integrate,function], " ^
|
s1210629013@55373
|
246 |
" [diff,integration,named]) " ^
|
s1210629013@55373
|
247 |
" [REAL (rhs M__M), REAL v_v, REAL_REAL y']); " ^
|
s1210629013@55373
|
248 |
" (B__B:: bool) = " ^
|
s1210629013@55373
|
249 |
" (SubProblem (Biegelinie',[named,integrate,function], " ^
|
s1210629013@55373
|
250 |
" [diff,integration,named]) " ^
|
s1210629013@55373
|
251 |
" [REAL (rhs N__N), REAL v_v, REAL_REAL y]); " ^
|
s1210629013@55373
|
252 |
" e__1 = NTH 1 r_b; " ^
|
s1210629013@55373
|
253 |
" (x__1::real) = argument_in (lhs e__1); " ^
|
s1210629013@55373
|
254 |
" (B__1::bool) = (Substitute [v_v = x__1]) B__B; " ^
|
s1210629013@55373
|
255 |
" B__1 = (Substitute [e__1]) B__1 ; " ^
|
s1210629013@55373
|
256 |
" B__2 = Take B__B; " ^
|
s1210629013@55373
|
257 |
" e__2 = NTH 2 r_b; " ^
|
s1210629013@55373
|
258 |
" (x__2::real) = argument_in (lhs e__2); " ^
|
s1210629013@55373
|
259 |
" (B__2::bool) = (Substitute [v_v = x__2]) B__B; " ^
|
s1210629013@55373
|
260 |
" B__2 = (Substitute [e__2]) B__2 ; " ^
|
s1210629013@55373
|
261 |
" (c_1_2::bool list) = " ^
|
s1210629013@55373
|
262 |
" (SubProblem (Biegelinie',[LINEAR,system],[no_met]) " ^
|
s1210629013@55373
|
263 |
" [BOOL_LIST [B__1, B__2], REAL_LIST [c,c_2]]); " ^
|
s1210629013@55373
|
264 |
" B__B = Take B__B; " ^
|
s1210629013@55373
|
265 |
" B__B = ((Substitute c_1_2) @@ " ^
|
s1210629013@55373
|
266 |
" (Rewrite_Set_Inst [(bdv, x)] make_ratpoly_in False)) B__B " ^
|
s1210629013@55373
|
267 |
" in B__B)"),
|
s1210629013@55373
|
268 |
prep_met thy "met_biege_2" [] e_metID
|
s1210629013@55373
|
269 |
(["IntegrierenUndKonstanteBestimmen2"],
|
s1210629013@55373
|
270 |
[("#Given" ,["Traegerlaenge l_l", "Streckenlast q__q", "FunktionsVariable v_v"]),
|
s1210629013@55373
|
271 |
(*("#Where",["0 < l_l"]), ...wait for < and handling Arbfix*)
|
s1210629013@55373
|
272 |
("#Find" ,["Biegelinie b_b"]),
|
s1210629013@55373
|
273 |
("#Relate",["Randbedingungen r_b"])],
|
s1210629013@55373
|
274 |
{rew_ord'="tless_true",
|
s1210629013@55373
|
275 |
rls' = append_rls "erls_IntegrierenUndK.." e_rls
|
s1210629013@55373
|
276 |
[Calc ("Atools.ident",eval_ident "#ident_"),
|
s1210629013@55373
|
277 |
Thm ("not_true",num_str @{thm not_true}),
|
s1210629013@55373
|
278 |
Thm ("not_false",num_str @{thm not_false})],
|
s1210629013@55373
|
279 |
calc = [],
|
s1210629013@55373
|
280 |
srls = append_rls "erls_IntegrierenUndK.." e_rls
|
s1210629013@55373
|
281 |
[Calc("Tools.rhs", eval_rhs"eval_rhs_"),
|
s1210629013@55373
|
282 |
Calc ("Atools.ident",eval_ident "#ident_"),
|
s1210629013@55373
|
283 |
Thm ("last_thmI",num_str @{thm last_thmI}),
|
s1210629013@55373
|
284 |
Thm ("if_True",num_str @{thm if_True}),
|
s1210629013@55373
|
285 |
Thm ("if_False",num_str @{thm if_False})],
|
s1210629013@55373
|
286 |
prls = Erls, crls = Atools_erls, errpats = [], nrls = Erls},
|
s1210629013@55373
|
287 |
"Script Biegelinie2Script " ^
|
s1210629013@55373
|
288 |
"(l_l::real) (q__q::real) (v_v::real) (b_b::real=>real) (r_b::bool list) = " ^
|
s1210629013@55373
|
289 |
" (let " ^
|
s1210629013@55373
|
290 |
" (fun_s:: bool list) = " ^
|
s1210629013@55373
|
291 |
" (SubProblem (Biegelinie',[vonBelastungZu,Biegelinien], " ^
|
s1210629013@55373
|
292 |
" [Biegelinien,ausBelastung]) " ^
|
s1210629013@55373
|
293 |
" [REAL q__q, REAL v_v]); " ^
|
s1210629013@55373
|
294 |
" (equ_s::bool list) = " ^
|
s1210629013@55373
|
295 |
" (SubProblem (Biegelinie',[setzeRandbedingungen,Biegelinien], " ^
|
s1210629013@55373
|
296 |
" [Biegelinien,setzeRandbedingungenEin]) " ^
|
s1210629013@55373
|
297 |
" [BOOL_LIST fun_s, BOOL_LIST r_b]); " ^
|
s1210629013@55373
|
298 |
" (con_s::bool list) = " ^
|
s1210629013@55373
|
299 |
" (SubProblem (Biegelinie',[LINEAR,system],[no_met]) " ^
|
s1210629013@55373
|
300 |
" [BOOL_LIST equ_s, REAL_LIST [c,c_2,c_3,c_4]]); " ^
|
s1210629013@55373
|
301 |
" B_B = Take (lastI fun_s); " ^
|
s1210629013@55373
|
302 |
" B_B = ((Substitute con_s) @@ " ^
|
s1210629013@55373
|
303 |
" (Rewrite_Set_Inst [(bdv, v_v)] make_ratpoly_in False)) B_B " ^
|
s1210629013@55373
|
304 |
" in B_B)"),
|
s1210629013@55373
|
305 |
prep_met thy "met_biege_intconst_2" [] e_metID
|
s1210629013@55373
|
306 |
(["IntegrierenUndKonstanteBestimmen","2xIntegrieren"], [],
|
s1210629013@55373
|
307 |
{rew_ord'="tless_true", rls'=Erls, calc = [], srls = e_rls, prls=e_rls, crls = Atools_erls,
|
s1210629013@55373
|
308 |
errpats = [], nrls = e_rls},
|
s1210629013@55373
|
309 |
"empty_script"),
|
s1210629013@55373
|
310 |
prep_met thy "met_biege_intconst_4" [] e_metID
|
s1210629013@55373
|
311 |
(["IntegrierenUndKonstanteBestimmen","4x4System"], [],
|
s1210629013@55373
|
312 |
{rew_ord'="tless_true", rls'=Erls, calc = [], srls = e_rls, prls=e_rls, crls = Atools_erls,
|
s1210629013@55373
|
313 |
errpats = [], nrls = e_rls},
|
s1210629013@55373
|
314 |
"empty_script"),
|
s1210629013@55373
|
315 |
prep_met thy "met_biege_intconst_1" [] e_metID
|
s1210629013@55373
|
316 |
(["IntegrierenUndKonstanteBestimmen","1xIntegrieren"], [],
|
s1210629013@55373
|
317 |
{rew_ord'="tless_true", rls'=Erls, calc = [], srls = e_rls, prls=e_rls, crls = Atools_erls,
|
s1210629013@55373
|
318 |
errpats = [], nrls = e_rls},
|
s1210629013@55373
|
319 |
"empty_script"),
|
s1210629013@55373
|
320 |
prep_met thy "met_biege2" [] e_metID
|
s1210629013@55373
|
321 |
(["Biegelinien"], [],
|
s1210629013@55373
|
322 |
{rew_ord'="tless_true", rls'=Erls, calc = [], srls = e_rls, prls=e_rls, crls = Atools_erls,
|
s1210629013@55373
|
323 |
errpats = [], nrls = e_rls},
|
s1210629013@55373
|
324 |
"empty_script"),
|
s1210629013@55373
|
325 |
prep_met thy "met_biege_ausbelast" [] e_metID
|
s1210629013@55373
|
326 |
(["Biegelinien", "ausBelastung"],
|
s1210629013@55373
|
327 |
[("#Given" ,["Streckenlast q__q", "FunktionsVariable v_v"]),
|
s1210629013@55373
|
328 |
("#Find" ,["Funktionen fun_s"])],
|
s1210629013@55373
|
329 |
{rew_ord'="tless_true",
|
s1210629013@55373
|
330 |
rls' = append_rls "erls_ausBelastung" e_rls
|
s1210629013@55373
|
331 |
[Calc ("Atools.ident", eval_ident "#ident_"),
|
s1210629013@55373
|
332 |
Thm ("not_true", num_str @{thm not_true}),
|
s1210629013@55373
|
333 |
Thm ("not_false", num_str @{thm not_false})],
|
s1210629013@55373
|
334 |
calc = [],
|
s1210629013@55373
|
335 |
srls = append_rls "srls_ausBelastung" e_rls
|
s1210629013@55373
|
336 |
[Calc ("Tools.rhs", eval_rhs "eval_rhs_")],
|
s1210629013@55373
|
337 |
prls = e_rls, crls = Atools_erls, errpats = [], nrls = e_rls},
|
s1210629013@55373
|
338 |
"Script Belastung2BiegelScript (q__q::real) (v_v::real) = " ^
|
s1210629013@55373
|
339 |
" (let q___q = Take (qq v_v = q__q); " ^
|
s1210629013@55373
|
340 |
" q___q = ((Rewrite sym_neg_equal_iff_equal True) @@ " ^
|
s1210629013@55373
|
341 |
" (Rewrite Belastung_Querkraft True)) q___q; " ^
|
s1210629013@55373
|
342 |
" (Q__Q:: bool) = " ^
|
s1210629013@55373
|
343 |
" (SubProblem (Biegelinie',[named,integrate,function], " ^
|
s1210629013@55373
|
344 |
" [diff,integration,named]) " ^
|
s1210629013@55373
|
345 |
" [REAL (rhs q___q), REAL v_v, REAL_REAL Q]); " ^
|
s1210629013@55373
|
346 |
" M__M = Rewrite Querkraft_Moment True Q__Q; " ^
|
s1210629013@55373
|
347 |
" (M__M::bool) = " ^
|
s1210629013@55373
|
348 |
" (SubProblem (Biegelinie',[named,integrate,function], " ^
|
s1210629013@55373
|
349 |
" [diff,integration,named]) " ^
|
s1210629013@55373
|
350 |
" [REAL (rhs M__M), REAL v_v, REAL_REAL M_b]); " ^
|
s1210629013@55373
|
351 |
" N__N = ((Rewrite Moment_Neigung False) @@ " ^
|
s1210629013@55373
|
352 |
" (Rewrite make_fun_explicit False)) M__M; " ^
|
s1210629013@55373
|
353 |
" (N__N:: bool) = " ^
|
s1210629013@55373
|
354 |
" (SubProblem (Biegelinie',[named,integrate,function], " ^
|
s1210629013@55373
|
355 |
" [diff,integration,named]) " ^
|
s1210629013@55373
|
356 |
" [REAL (rhs N__N), REAL v_v, REAL_REAL y']); " ^
|
s1210629013@55373
|
357 |
" (B__B:: bool) = " ^
|
s1210629013@55373
|
358 |
" (SubProblem (Biegelinie',[named,integrate,function], " ^
|
s1210629013@55373
|
359 |
" [diff,integration,named]) " ^
|
s1210629013@55373
|
360 |
" [REAL (rhs N__N), REAL v_v, REAL_REAL y]) " ^
|
s1210629013@55373
|
361 |
" in [Q__Q, M__M, N__N, B__B])"),
|
s1210629013@55373
|
362 |
prep_met thy "met_biege_setzrand" [] e_metID
|
s1210629013@55373
|
363 |
(["Biegelinien", "setzeRandbedingungenEin"],
|
s1210629013@55373
|
364 |
[("#Given" , ["Funktionen fun_s", "Randbedingungen r_b"]),
|
s1210629013@55373
|
365 |
("#Find" , ["Gleichungen equs'''"])],
|
s1210629013@55373
|
366 |
{rew_ord'="tless_true", rls'=Erls, calc = [], srls = srls2, prls=e_rls, crls = Atools_erls,
|
s1210629013@55373
|
367 |
errpats = [], nrls = e_rls},
|
s1210629013@55373
|
368 |
"Script SetzeRandbedScript (fun_s::bool list) (r_b::bool list) = " ^
|
s1210629013@55373
|
369 |
" (let b_1 = NTH 1 r_b; " ^
|
s1210629013@55373
|
370 |
" f_s = filter_sameFunId (lhs b_1) fun_s; " ^
|
s1210629013@55373
|
371 |
" (e_1::bool) = " ^
|
s1210629013@55373
|
372 |
" (SubProblem (Biegelinie',[makeFunctionTo,equation]," ^
|
s1210629013@55373
|
373 |
" [Equation,fromFunction]) " ^
|
s1210629013@55373
|
374 |
" [BOOL (hd f_s), BOOL b_1]); " ^
|
s1210629013@55373
|
375 |
" b_2 = NTH 2 r_b; " ^
|
s1210629013@55373
|
376 |
" f_s = filter_sameFunId (lhs b_2) fun_s; " ^
|
s1210629013@55373
|
377 |
" (e_2::bool) = " ^
|
s1210629013@55373
|
378 |
" (SubProblem (Biegelinie',[makeFunctionTo,equation]," ^
|
s1210629013@55373
|
379 |
" [Equation,fromFunction]) " ^
|
s1210629013@55373
|
380 |
" [BOOL (hd f_s), BOOL b_2]); " ^
|
s1210629013@55373
|
381 |
" b_3 = NTH 3 r_b; " ^
|
s1210629013@55373
|
382 |
" f_s = filter_sameFunId (lhs b_3) fun_s; " ^
|
s1210629013@55373
|
383 |
" (e_3::bool) = " ^
|
s1210629013@55373
|
384 |
" (SubProblem (Biegelinie',[makeFunctionTo,equation]," ^
|
s1210629013@55373
|
385 |
" [Equation,fromFunction]) " ^
|
s1210629013@55373
|
386 |
" [BOOL (hd f_s), BOOL b_3]); " ^
|
s1210629013@55373
|
387 |
" b_4 = NTH 4 r_b; " ^
|
s1210629013@55373
|
388 |
" f_s = filter_sameFunId (lhs b_4) fun_s; " ^
|
s1210629013@55373
|
389 |
" (e_4::bool) = " ^
|
s1210629013@55373
|
390 |
" (SubProblem (Biegelinie',[makeFunctionTo,equation]," ^
|
s1210629013@55373
|
391 |
" [Equation,fromFunction]) " ^
|
s1210629013@55373
|
392 |
" [BOOL (hd f_s), BOOL b_4]) " ^
|
s1210629013@55373
|
393 |
" in [e_1, e_2, e_3, e_4])"
|
s1210629013@55373
|
394 |
(* filter requires more than 1 sec !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
|
s1210629013@55373
|
395 |
"Script SetzeRandbedScript (fun_s::bool list) (r_b::bool list) = " ^
|
s1210629013@55373
|
396 |
" (let b_1 = NTH 1 r_b; " ^
|
s1210629013@55373
|
397 |
" f_s = filter (sameFunId (lhs b_1)) fun_s; " ^
|
s1210629013@55373
|
398 |
" (e_1::bool) = " ^
|
s1210629013@55373
|
399 |
" (SubProblem (Biegelinie',[makeFunctionTo,equation]," ^
|
s1210629013@55373
|
400 |
" [Equation,fromFunction]) " ^
|
s1210629013@55373
|
401 |
" [BOOL (hd f_s), BOOL b_1]); " ^
|
s1210629013@55373
|
402 |
" b_2 = NTH 2 r_b; " ^
|
s1210629013@55373
|
403 |
" f_s = filter (sameFunId (lhs b_2)) fun_s; " ^
|
s1210629013@55373
|
404 |
" (e_2::bool) = " ^
|
s1210629013@55373
|
405 |
" (SubProblem (Biegelinie',[makeFunctionTo,equation]," ^
|
s1210629013@55373
|
406 |
" [Equation,fromFunction]) " ^
|
s1210629013@55373
|
407 |
" [BOOL (hd f_s), BOOL b_2]); " ^
|
s1210629013@55373
|
408 |
" b_3 = NTH 3 r_b; " ^
|
s1210629013@55373
|
409 |
" f_s = filter (sameFunId (lhs b_3)) fun_s; " ^
|
s1210629013@55373
|
410 |
" (e_3::bool) = " ^
|
s1210629013@55373
|
411 |
" (SubProblem (Biegelinie',[makeFunctionTo,equation]," ^
|
s1210629013@55373
|
412 |
" [Equation,fromFunction]) " ^
|
s1210629013@55373
|
413 |
" [BOOL (hd f_s), BOOL b_3]); " ^
|
s1210629013@55373
|
414 |
" b_4 = NTH 4 r_b; " ^
|
s1210629013@55373
|
415 |
" f_s = filter (sameFunId (lhs b_4)) fun_s; " ^
|
s1210629013@55373
|
416 |
" (e_4::bool) = " ^
|
s1210629013@55373
|
417 |
" (SubProblem (Biegelinie',[makeFunctionTo,equation]," ^
|
s1210629013@55373
|
418 |
" [Equation,fromFunction]) " ^
|
s1210629013@55373
|
419 |
" [BOOL (hd f_s), BOOL b_4]) " ^
|
s1210629013@55373
|
420 |
" in [e_1,e_2,e_3,e_4])"*)),
|
s1210629013@55373
|
421 |
prep_met thy "met_equ_fromfun" [] e_metID
|
s1210629013@55373
|
422 |
(["Equation","fromFunction"],
|
s1210629013@55373
|
423 |
[("#Given" ,["functionEq fu_n","substitution su_b"]),
|
s1210629013@55373
|
424 |
("#Find" ,["equality equ'''"])],
|
s1210629013@55373
|
425 |
{rew_ord'="tless_true", rls'=Erls, calc = [],
|
s1210629013@55373
|
426 |
srls = append_rls "srls_in_EquationfromFunc" e_rls
|
s1210629013@55373
|
427 |
[Calc("Tools.lhs", eval_lhs"eval_lhs_"),
|
s1210629013@55373
|
428 |
Calc("Atools.argument'_in", eval_argument_in "Atools.argument'_in")],
|
s1210629013@55373
|
429 |
prls=e_rls, crls = Atools_erls, errpats = [], nrls = e_rls},
|
s1210629013@55373
|
430 |
(*(M_b x = c_2 + c * x + -1 * q_0 / 2 * x ^^^ 2) (M_b L = 0) -->
|
s1210629013@55373
|
431 |
0 = c_2 + c * L + -1 * q_0 / 2 * L ^^^ 2*)
|
s1210629013@55373
|
432 |
"Script Function2Equality (fu_n::bool) (su_b::bool) =" ^
|
s1210629013@55373
|
433 |
" (let fu_n = Take fu_n; " ^
|
s1210629013@55373
|
434 |
" bd_v = argument_in (lhs fu_n); " ^
|
s1210629013@55373
|
435 |
" va_l = argument_in (lhs su_b); " ^
|
s1210629013@55373
|
436 |
" eq_u = (Substitute [bd_v = va_l]) fu_n; " ^
|
s1210629013@55373
|
437 |
(*([1], Res), M_b L = c_2 + c * L + -1 * q_0 / 2 * L ^^^ 2*)
|
s1210629013@55373
|
438 |
" eq_u = (Substitute [su_b]) eq_u " ^
|
s1210629013@55373
|
439 |
" in (Rewrite_Set norm_Rational False) eq_u) ")]
|
s1210629013@55373
|
440 |
*}
|
s1210629013@55373
|
441 |
|
neuper@37906
|
442 |
end
|
neuper@37906
|
443 |
|