def elem_dep_data_qt10ms(na, nb, nt):
fens, fes = t4_block(a, b, t, na, nb, nt, orientation='a')
fens, fes = t4_to_t10(fens, fes)
femm = FEMMDeforLinearQT10MS(material=m, material_csys=mcsys, fes=fes)
bfes = mesh_boundary(femm.fes)
fesel = fe_select(fens, bfes, box=[a, a, 0, b, 0, t], inflate=htol)
tsfes = bfes.subset(fesel)
sfemm = FEMMDeforLinear(fes=tsfes, integration_rule=TriRule(npts=3))
return fens, femm, sfemm
E = 1000
nu = 0.4999
W = 2.5
H = 5
L = 50
# nW, nL, nH = 20, 20, 20
nW, nL, nH = 4, 20, 4
htol = min(L, H, W) / 1000
magn = -0.2 * 12.2334 / 4
Force = magn * W * H * 2
Force * L**3 / (3 * E * W * H**3 * 2 / 12)
uzex = -12.0935378981478
m = MatDeforTriaxLinearIso(e=E, nu=nu)
start = time.time()
fens, fes = t4_block(W, L, H, nW, nL, nH, orientation='ca')
fens, fes = t4_to_t10(fens, fes)
print('Mesh generation', time.time() - start)
model_data = {}
model_data['fens'] = fens
model_data['regions'] = [{'femm': FEMMDeforLinearQT10MS(material=m, fes=fes)}]
algo_common.plot_mesh(model_data)
model_data['boundary_conditions'] = {}
# Clamped face
cn = fenode_select(fens, box=numpy.array([0, W, 0, 0, 0, H]), inflate=htol)
essential = [
{
'node_list': cn,
