elasticity (beam).py
4c9bfd2ccc95master

elasticity (beam).py
View Options

			from rrompy.hfengines.vector_linear_problem import \
			LinearElasticityBeamPoissonRatio as LEBPR

			# Scaled variables
			scaling = 10

			L = 5.
			g = 9.8 / scaling
			E = 1e6
			nu = .47
			rho = 1.5e3
			lambda_ = E * nu / (1. + nu) / (1. - 2. * nu)
			mu_ = E / (1. + nu)
			n = 50

			# Create mesh and define function space
			solver = LEBPR(n = n, rho_ = rho, g = g, E = E, nu0 = nu, length = L)

			uh = solver.solve(nu)
			solver.plot(uh, what = ["REAL"], figsize = (6, 3))