1 |
REQUIRE "atoms.a4l"; |
2 |
|
3 |
(* a parameterised simply-supported beam model *) |
4 |
|
5 |
MODEL beam_parameterised( |
6 |
n WILL_BE integer_constant; |
7 |
E WILL_BE youngs_modulus; |
8 |
I WILL_BE second_moment_of_area; |
9 |
L WILL_BE distance; |
10 |
); |
11 |
v[1..n] IS_A distance; |
12 |
x[1..n] IS_A distance; |
13 |
|
14 |
P IS_A force; |
15 |
a,b IS_A distance; |
16 |
a + b = L; |
17 |
|
18 |
isrightp[1..n] IS_A boolean_var; |
19 |
|
20 |
FOR i IN [1..n] CREATE |
21 |
v_left[i]: v[i] = P*b/(6*E*I*L)*((L^2 - b^2)*x[i] - x[i]^3); |
22 |
v_right[i]: v[i] = P*b/(6*E*I*L)*((L^2 - b^2)*x[i] - x[i]^3 + (L/b)*(x[i]-a)^3 ); |
23 |
END FOR; |
24 |
|
25 |
CONDITIONAL |
26 |
FOR i IN [1..n] CREATE |
27 |
isright[i]: x[i] > a; |
28 |
END FOR; |
29 |
END CONDITIONAL; |
30 |
|
31 |
FOR i IN [1..n] CREATE |
32 |
isrightp[i] == SATISFIED(isright[i]); |
33 |
|
34 |
WHEN (isrightp[i]) |
35 |
CASE TRUE: |
36 |
USE v_right[i]; |
37 |
CASE FALSE: |
38 |
USE v_left[i]; |
39 |
END WHEN; |
40 |
END FOR; |
41 |
METHODS |
42 |
METHOD specify; |
43 |
FIX P, a, x[1..n]; |
44 |
END specify; |
45 |
METHOD values; |
46 |
a.upper_bound := 100 {m}; |
47 |
b.upper_bound := 100 {m}; |
48 |
L.upper_bound := 100 {m}; |
49 |
END values; |
50 |
|
51 |
END beam_parameterised; |
52 |
|
53 |
(* |
54 |
Model of a simply-supported beam of length L |
55 |
with a single vertical point load P at 0 < a < L |
56 |
*) |
57 |
|
58 |
MODEL beam; |
59 |
n IS_A integer_constant; |
60 |
n :== 1; |
61 |
E IS_A youngs_modulus; |
62 |
I IS_A second_moment_of_area; |
63 |
L IS_A distance; |
64 |
B IS_A beam_parameterised(n,E,I,L); |
65 |
|
66 |
P ALIASES B.P; |
67 |
a ALIASES B.a; |
68 |
b ALIASES B.b; |
69 |
|
70 |
x ALIASES B.x[1]; |
71 |
v ALIASES B.v[1]; |
72 |
|
73 |
METHODS |
74 |
METHOD specify; |
75 |
FIX E, I, L; |
76 |
FIX P, a; |
77 |
FIX x; |
78 |
END specify; |
79 |
|
80 |
METHOD values; |
81 |
E := 200 {GPa}; (* approx, for steel *) |
82 |
I := 142e6 {mm^4}; (* AISC Universal Beam 360UB50.7 *) |
83 |
L := 3.5 {m}; |
84 |
P := 14.0 {kN}; |
85 |
a := 1.75 {m}; |
86 |
x := 2.0 {m}; |
87 |
END values; |
88 |
|
89 |
METHOD on_load; |
90 |
RUN reset; |
91 |
RUN values; |
92 |
END on_load; |
93 |
|
94 |
END beam; |