-
Notifications
You must be signed in to change notification settings - Fork 1
/
matlab00.jl
43 lines (36 loc) · 912 Bytes
/
matlab00.jl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
using Compat, ClassicalLaminateTheory
forcesAndMoments =
@compat Dict(
:Nx => 0,
:Ny => 0,
:Nxy => 0,
:Mx => 0,
:My => 0,
:Mxy => 0)
materialProperties =@compat Dict(
:E1 => 20.0*10^6,
:E2 => 1.5*10^6,
:G12 => 1.0*10^6,
:ν12 => .34,
:alpha1 => 0.2e-6,
:alpha2 => 20e-6,
:beta1 => 0.2e-4,
:beta2 => 20e-4,
:thickness => 0.0075)
laminateProperties =@compat Dict(
:nplies => 10,
:repeats => 3,
:symmetric => true,
:orientation => [90, 45, -45, 90, 45, 45, 90, -45, 45, 90],
:deltaTemp => -280.0,
:deltaMoisture => 0
)
println()
#createLaminate!(laminateProperties, materialProperties)
createLaminate!(laminateProperties, materialProperties, forcesAndMoments)
show(laminateProperties)
println("ABD matrix:")
round(laminateProperties[:ABD], 3) |> display
println()
println("Total forces and moments:")
round(laminateProperties[:fTotal], 7) |> display