stiffness_matrix.py
No OneTemporary
Actions

Subscribers

None

File Metadata

	#!/usr/bin/env python3
	""" stiffness_matrix.py: Example showing how to access the stifness matrix in
	python"""

	__author__ = "Guillaume Anciaux"
	__credits__ = [
	"Guillaume Anciaux <guillaume.anciaux@epfl.ch>",
	]
	__copyright__ = "Copyright (©) 2018-2021 EPFL (Ecole Polytechnique Fédérale" \
	" de Lausanne) Laboratory (LSMS - Laboratoire de Simulation" \
	" en Mécanique des Solides)"
	__license__ = "LGPLv3"

	import akantu as aka


	def getStiffnessMatrix(material_file, mesh_file, traction):
	aka.parseInput(material_file)
	spatial_dimension = 2

	# --------------------------------------------------------------------------
	# Initialization
	# --------------------------------------------------------------------------
	mesh = aka.Mesh(spatial_dimension)
	mesh.read(mesh_file)

	model = aka.SolidMechanicsModel(mesh)
	model.initFull(_analysis_method=aka._static)

	model.assembleStiffnessMatrix()
	K = model.getDOFManager().getMatrix('K')
	stiff = aka.AkantuSparseMatrix(K).toarray()

	return stiff


	# --------------------------------------------------------------------------
	# main
	# --------------------------------------------------------------------------
	def main():
	mesh_file = 'plate.msh'
	material_file = 'material.dat'

	traction = 1.
	mat = getStiffnessMatrix(material_file, mesh_file, traction)
	print(mat)


	# --------------------------------------------------------------------------
	if __name__ == "__main__":
	main()