def get_pars(dim, lam, mu): c = stiffness_from_lame(3, lam, mu) if dim == 2: tr = TransformToPlane() try: c = tr.tensor_plane_stress(c3=c) except: sym = (dim + 1) * dim // 2 c = nm.zeros((sym, sym), dtype=nm.float64) return c
def get_pars(dim, lam, mu): c = stiffness_from_lame(3, lam, mu) if dim == 2: tr = TransformToPlane() try: c = tr.tensor_plane_stress(c3=c) except: sym = (dim + 1) * dim / 2 c = nm.zeros((sym, sym), dtype=nm.float64) return c
def get_pars( lam, mu, dim, full = False ): from sfepy.mechanics.matcoefs import stiffness_tensor_lame, TransformToPlane if full: c = stiffness_tensor_lame( 3, lam, mu ) if dim == 2: tr = TransformToPlane() try: c = tr.tensor_plane_stress( c3 = c ) except: sym = (dim + 1) * dim / 2 c = nm.zeros( (sym, sym), dtype = nm.float64 ) return c else: return lam, mu