def test_pm45(): # {{{1
A = lamina(hf, hr, 100, 45, 0.5)
B = lamina(hf, hr, 100, -45, 0.5)
pm45 = laminate('pm45', [A, B, B, A])
assert math.isclose(pm45.thickness, 0.4545, rel_tol=0.01)
assert math.isclose(pm45.ρ, 1.43, rel_tol=0.01)
assert pm45.vf == 0.5
assert math.isclose(pm45.wf, 0.615, rel_tol=0.01)
assert math.isclose(pm45.Ex, 17899.5, rel_tol=0.01)
assert math.isclose(pm45.Ey, 17899.5, rel_tol=0.01)
assert math.isclose(pm45.Ez, 16979.8, rel_tol=0.01)
assert math.isclose(pm45.Gxy, 31928.0, rel_tol=0.01)
assert math.isclose(pm45.Gxz, 4608.06, rel_tol=0.01)
assert math.isclose(pm45.Gyz, 4608.06, rel_tol=0.01)
assert math.isclose(pm45.νxy, 0.756, rel_tol=0.01)
assert math.isclose(pm45.νyx, 0.756, rel_tol=0.01)
assert math.isclose(pm45.αx, 6.4269e-06, rel_tol=0.01)
assert math.isclose(pm45.αy, 6.4269e-06, rel_tol=0.01)
assert math.isclose(pm45.tEx, 17899.5, rel_tol=0.01)
assert math.isclose(pm45.tEy, 17899.5, rel_tol=0.01)
assert math.isclose(pm45.tEz, 19301.8, rel_tol=0.01)
assert math.isclose(pm45.tGxy, 31935.1, rel_tol=0.01)
assert math.isclose(pm45.tGyz, 5529.6, rel_tol=0.01)
assert math.isclose(pm45.tGxz, 5529.6, rel_tol=0.01)
assert math.isclose(pm45.tνxy, 0.75, rel_tol=0.01)
assert math.isclose(pm45.tνxz, 0.102, rel_tol=0.01)
assert math.isclose(pm45.tνyz, 0.102, rel_tol=0.01)
def test_ud(): # {{{1
la = lamina(hf, hr, 100, 0, 0.5)
ud = laminate('ud', [la, la, la, la])
assert math.isclose(ud.thickness, 0.4545, rel_tol=0.01)
assert math.isclose(ud.ρ, 1.43, rel_tol=0.01)
assert ud.vf == 0.5
assert math.isclose(ud.wf, 0.615, rel_tol=0.01)
assert math.isclose(ud.Ex, 118810, rel_tol=0.01)
assert math.isclose(ud.Ey, 16979.8, rel_tol=0.01)
assert math.isclose(ud.Ez, 16979.8, rel_tol=0.01)
assert math.isclose(ud.Gxy, 5095.588, rel_tol=0.01)
assert math.isclose(ud.Gxz, 4246.32, rel_tol=0.01)
assert math.isclose(ud.Gyz, 4969.80, rel_tol=0.01)
assert math.isclose(ud.νxy, 0.28, rel_tol=0.01)
assert math.isclose(ud.νyx, 0.04, rel_tol=0.01)
assert math.isclose(ud.αx, 2.75e-07, rel_tol=0.01)
assert math.isclose(ud.αy, 4.14e-5, rel_tol=0.01)
assert math.isclose(ud.tEx, 118810.0, rel_tol=0.01)
assert math.isclose(ud.tEy, 16979.80, rel_tol=0.01)
assert math.isclose(ud.tEz, 16979.80, rel_tol=0.01)
assert math.isclose(ud.tGxy, 5095.58, rel_tol=0.01)
assert math.isclose(ud.tGyz, 5963.76, rel_tol=0.01)
assert math.isclose(ud.tGxz, 5095.58, rel_tol=0.01)
assert math.isclose(ud.tνxy, 0.2799, rel_tol=0.01)
assert math.isclose(ud.tνxz, 0.2799, rel_tol=0.01)
assert math.isclose(ud.tνyz, 0.4235, rel_tol=0.01)
def test_plain_weave(): # {{{1
A = lamina(hf, hr, 100, 0, 0.5)
B = lamina(hf, hr, 100, 90, 0.5)
pw = laminate('pw', [A, B, B, A])
assert math.isclose(pw.thickness, 0.4545, rel_tol=0.01)
assert math.isclose(pw.ρ, 1.43, rel_tol=0.01)
assert pw.vf == 0.5
assert math.isclose(pw.wf, 0.615, rel_tol=0.01)
assert math.isclose(pw.Ex, 68327.5, rel_tol=0.01)
assert math.isclose(pw.Ey, 68327.5, rel_tol=0.01)
assert math.isclose(pw.Ez, 16979.8, rel_tol=0.01)
assert math.isclose(pw.Gxy, 5095.588, rel_tol=0.01)
assert math.isclose(pw.Gxz, 4743.715, rel_tol=0.01)
assert math.isclose(pw.Gyz, 4472.41, rel_tol=0.01)
assert math.isclose(pw.νxy, 0.07, rel_tol=0.01)
assert math.isclose(pw.νyx, 0.07, rel_tol=0.01)
assert math.isclose(pw.αx, 6.427e-06, rel_tol=0.01)
assert math.isclose(pw.αy, 6.427e-06, rel_tol=0.01)
assert math.isclose(pw.tEx, 68335.7, rel_tol=0.01)
assert math.isclose(pw.tEy, 68335.7, rel_tol=0.01)
assert math.isclose(pw.tEz, 19301.8, rel_tol=0.01)
assert math.isclose(pw.tGxy, 5095.588, rel_tol=0.01)
assert math.isclose(pw.tGyz, 5529.675, rel_tol=0.01)
assert math.isclose(pw.tGxz, 5529.675, rel_tol=0.01)
assert math.isclose(pw.tνxy, 0.07, rel_tol=0.01)
assert math.isclose(pw.tνxz, 0.39, rel_tol=0.01)
assert math.isclose(pw.tνyz, 0.39, rel_tol=0.01)
def test_qi(): # {{{1
A = lamina(hf, hr, 100, 0, 0.5)
B = lamina(hf, hr, 100, 90, 0.5)
C = lamina(hf, hr, 100, 45, 0.5)
D = lamina(hf, hr, 100, -45, 0.5)
qi = laminate('qi', [A, B, C, D, D, C, B, A])
assert math.isclose(qi.thickness, 0.9090, rel_tol=0.01)
assert math.isclose(qi.ρ, 1.43, rel_tol=0.01)
assert qi.vf == 0.5
assert math.isclose(qi.wf, 0.615, rel_tol=0.01)
assert math.isclose(qi.Ex, 49236.42, rel_tol=0.01)
assert math.isclose(qi.Ey, 49236.42, rel_tol=0.01)
assert math.isclose(qi.Ez, 16979.8, rel_tol=0.01)
assert math.isclose(qi.Gxy, 18511.8, rel_tol=0.01)
assert math.isclose(qi.Gxz, 4658.93, rel_tol=0.01)
assert math.isclose(qi.Gyz, 4557.19, rel_tol=0.01)
assert math.isclose(qi.νxy, 0.3298, rel_tol=0.01)
assert math.isclose(qi.νyx, 0.3298, rel_tol=0.01)
assert math.isclose(qi.αx, 6.42695e-06, rel_tol=0.01)
assert math.isclose(qi.αy, 6.42695e-06, rel_tol=0.01)
assert math.isclose(qi.tEx, 49242.7, rel_tol=0.01)
assert math.isclose(qi.tEy, 49242.7, rel_tol=0.01)
assert math.isclose(qi.tEz, 19301.8, rel_tol=0.01)
assert math.isclose(qi.tGxy, 18515.3, rel_tol=0.01)
assert math.isclose(qi.tGyz, 5529.6, rel_tol=0.01)
assert math.isclose(qi.tGxz, 5529.6, rel_tol=0.01)
assert math.isclose(qi.tνxy, 0.32977, rel_tol=0.01)
assert math.isclose(qi.tνxz, 0.28244, rel_tol=0.01)
assert math.isclose(qi.tνyz, 0.28244, rel_tol=0.01)