Orthotropy/Anisotropy as material choice

[considinefarm]

Is it possible to add orthotropy/anisotropy stiffness tensor? I'm not sure how to include them, but would contribute with some guidance. #enhancement

[nicoguaro]

@considinefarm, the assembler of SolidsPy has an optional parameter (uel) that allows you to pass a function (callable) that defines the local stiffness matrix. You could use this option to assemble elements with anisotropic materials.

In that case, you would a script like the following

import matplotlib.pyplot as plt
import numpy as np
import solidspy.preprocesor as pre
import solidspy.postprocesor as pos
import solidspy.assemutil as ass
import solidspy.solutil as sol

def aniso_el():
	"""Anisotropic element"""
	pass

nodes, mats, elements, loads = pre.readin(folder=folder)
DME, IBC , neq = ass.DME(nodes, elements)
KG = ass.assembler(elements, mats, nodes, neq, DME, uel=aniso_el)
RHSG = ass.loadasem(loads, IBC, neq)
UG = sol.static_sol(KG, RHSG)
UC = pos.complete_disp(IBC, nodes, UG)
E_nodes, S_nodes = pos.strain_nodes(nodes , elements, mats, UC)
pos.fields_plot(elements, nodes, UC, E_nodes=E_nodes, S_nodes=S_nodes)

[considinefarm]

Thanks – I’ll give it a try. From: Nicolás Guarín-Zapata <notifications@github.com> Sent: Monday, May 18, 2020 4:49 PM To: AppliedMechanics-EAFIT/SolidsPy <SolidsPy@noreply.github.com> Cc: JOHN M CONSIDINE <john.considine@wisc.edu>; Mention <mention@noreply.github.com> Subject: Re: [AppliedMechanics-EAFIT/SolidsPy] Orthotropy/Anisotropy as material choice (#21) @considinefarm<https://github.com/considinefarm>, the assembler of SolidsPy has an optional parameter (uel) that allows you to pass a function (callable) that defines the local stiffness matrix. You could use this option to assemble elements with anisotropic materials. In that case, you would a script like the following import matplotlib.pyplot as plt import numpy as np import solidspy.preprocesor as pre import solidspy.postprocesor as pos import solidspy.assemutil as ass import solidspy.solutil as sol def aniso_el(): """Anisotropic element""" pass nodes, mats, elements, loads = pre.readin(folder=folder) DME, IBC , neq = ass.DME(nodes, elements) KG = ass.assembler(elements, mats, nodes, neq, DME, uel=aniso_el) RHSG = ass.loadasem(loads, IBC, neq) UG = sol.static_sol(KG, RHSG) UC = pos.complete_disp(IBC, nodes, UG) E_nodes, S_nodes = pos.strain_nodes(nodes , elements, mats, UC) pos.fields_plot(elements, nodes, UC, E_nodes=E_nodes, S_nodes=S_nodes) — You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub<#21 (comment)>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/ABKFVLWAR6GKCMYSV6O7DS3RSGUNLANCNFSM4M725LAA>.