Diffusion akantu (bugfixes/all_bugs_found_by_working_on_tutorials) Browse

all_bugs_found…

ntn_base_contact.hh
View Options

			/**
			* @file ntn_base_contact.hh
			*
			* @author David Simon Kammer <david.kammer@epfl.ch>
			*
			* @date creation: Fri Mar 16 2018
			* @date last modification: Tue Sep 29 2020
			*
			* @brief base contact for ntn and ntrf contact
			*
			*
			* @section LICENSE
			*
			* Copyright (©) 2015-2021 EPFL (Ecole Polytechnique Fédérale de Lausanne)
			* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
			*
			* Akantu is free software: you can redistribute it and/or modify it under the
			* terms of the GNU Lesser General Public License as published by the Free
			* Software Foundation, either version 3 of the License, or (at your option) any
			* later version.
			*
			* Akantu is distributed in the hope that it will be useful, but WITHOUT ANY
			* WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
			* A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more
			* details.
			*
			* You should have received a copy of the GNU Lesser General Public License
			* along with Akantu. If not, see <http://www.gnu.org/licenses/>.
			*
			*/

			/* -------------------------------------------------------------------------- */
			#ifndef AST_NTN_BASE_CONTACT_HH_
			#define AST_NTN_BASE_CONTACT_HH_

			/* -------------------------------------------------------------------------- */
			// akantu
			#include "aka_csr.hh"
			#include "solid_mechanics_model.hh"

			// simtools
			#include "synchronized_array.hh"

			namespace akantu {

			class NTNBaseContact;

			/* -------------------------------------------------------------------------- */
			// class NTNContactSynchElementFilter : public SynchElementFilter {
			// public:
			// // constructor
			// NTNContactSynchElementFilter(NTNBaseContact & contact);

			// // answer to: do we need this element ?
			// virtual bool operator()(const Element & e);

			// private:
			// const NTNBaseContact & contact;
			// const ElementTypeMapArray<UInt> & connectivity;
			// };

			/* -------------------------------------------------------------------------- */
			class NTNBaseContact : public DataAccessor<Element>, public Dumpable {
			/* ------------------------------------------------------------------------ */
			/* Constructors/Destructors */
			/* ------------------------------------------------------------------------ */
			public:
			NTNBaseContact(SolidMechanicsModel & model, const ID & id = "contact");
			~NTNBaseContact() override;

			/* ------------------------------------------------------------------------ */
			/* Methods */
			/* ------------------------------------------------------------------------ */
			public:
			/// initializes ntn contact parallel
			virtual void initParallel();

			/// add split node
			virtual void addSplitNode(UInt node, UInt /unused/ = 0);

			/// update normals, lumped boundary, and impedance
			virtual void updateInternalData();

			/// update (compute the normals)
			virtual void updateNormals() = 0;

			/// update the lumped boundary B matrix
			virtual void updateLumpedBoundary();

			/// update the impedance matrix
			virtual void updateImpedance() = 0;

			/// compute the normal contact force
			virtual void computeContactPressure();

			/// impose the normal contact force
			virtual void applyContactPressure();

			/// register synchronizedarrays for sync
			virtual void registerSynchronizedArray(SynchronizedArrayBase & array);

			/// dump restart file
			virtual void dumpRestart(const std::string & file_name) const;

			/// read restart file
			virtual void readRestart(const std::string & file_name);

			/// compute the normal gap
			virtual void computeNormalGap(Array<Real> & gap) const = 0;

			/// compute relative normal field (only value that has to be multiplied with
			/// the normal)
			/// relative to master nodes
			virtual void
			computeRelativeNormalField(const Array<Real> & field,
			Array<Real> & rel_normal_field) const = 0;

			/// compute relative tangential field (complet array)
			/// relative to master nodes
			virtual void
			computeRelativeTangentialField(const Array<Real> & field,
			Array<Real> & rel_tang_field) const = 0;

			/// function to print the contain of the class
			virtual void printself(std::ostream & stream, int indent = 0) const;

			/// computes the acceleration
			void computeAcceleration(Array<Real> & acceleration) const;

			protected:
			/// updateLumpedBoundary
			virtual void
			internalUpdateLumpedBoundary(const Array<UInt> & nodes,
			const ElementTypeMapArray<UInt> & elements,
			SynchronizedArray<Real> & boundary);

			// to find the slave_elements or master_elements
			virtual void findBoundaryElements(const Array<UInt> & interface_nodes,
			ElementTypeMapArray<UInt> & elements);

			/// synchronize arrays
			virtual void syncArrays(SyncChoice sync_choice);

			/* ------------------------------------------------------------------------ */
			/* Constructors/Destructors */
			/* ------------------------------------------------------------------------ */
			public:
			inline UInt getNbData(const Array<Element> & elements,
			const SynchronizationTag & tag) const override;

			inline void packData(CommunicationBuffer & buffer,
			const Array<Element> & elements,
			const SynchronizationTag & tag) const override;

			inline void unpackData(CommunicationBuffer & buffer,
			const Array<Element> & elements,
			const SynchronizationTag & tag) override;

			/* ------------------------------------------------------------------------ */
			/* Dumpable */
			/* ------------------------------------------------------------------------ */
			public:
			void addDumpFieldToDumper(const std::string & dumper_name,
			const std::string & field_id) override;

			/* ------------------------------------------------------------------------ */
			/* Accessors */
			/* ------------------------------------------------------------------------ */
			public:
			AKANTU_GET_MACRO(Model, model, SolidMechanicsModel &)

			AKANTU_GET_MACRO(Slaves, slaves, const SynchronizedArray<UInt> &)
			AKANTU_GET_MACRO(Normals, normals, const SynchronizedArray<Real> &)
			AKANTU_GET_MACRO(ContactPressure, contact_pressure,
			const SynchronizedArray<Real> &)
			AKANTU_GET_MACRO(LumpedBoundarySlaves, lumped_boundary_slaves,
			const SynchronizedArray<Real> &)
			AKANTU_GET_MACRO(Impedance, impedance, const SynchronizedArray<Real> &)
			AKANTU_GET_MACRO(IsInContact, is_in_contact, const SynchronizedArray<bool> &)

			AKANTU_GET_MACRO(SlaveElements, slave_elements,
			const ElementTypeMapArray<UInt> &)

			AKANTU_GET_MACRO(SynchronizerRegistry, *synch_registry,
			SynchronizerRegistry &)

			/// get number of nodes that are in contact (globally, on all procs together)
			/// is_in_contact = true
			virtual UInt getNbNodesInContact() const;

			/// get index of node in either slaves or masters array
			/// if node is in neither of them, return -1
			virtual Int getNodeIndex(UInt node) const;

			/// get number of contact nodes: nodes in the system locally (on this proc)
			/// is_in_contact = true and false, because just in the system
			virtual UInt getNbContactNodes() const { return this->slaves.size(); }

			bool isNTNContact() const { return this->is_ntn_contact; }

			/* ------------------------------------------------------------------------ */
			/* Class Members */
			/* ------------------------------------------------------------------------ */
			protected:
			using SurfacePtrSet = std::set<const ElementGroup *>;

			ID id;

			SolidMechanicsModel & model;

			/// array of slave nodes
			SynchronizedArray<UInt> slaves;
			/// array of normals
			SynchronizedArray<Real> normals;
			/// array indicating if nodes are in contact
			SynchronizedArray<Real> contact_pressure;
			/// array indicating if nodes are in contact
			SynchronizedArray<bool> is_in_contact;
			/// boundary matrix for slave nodes
			SynchronizedArray<Real> lumped_boundary_slaves;
			/// impedance matrix
			SynchronizedArray<Real> impedance;

			/// contact surface
			SurfacePtrSet contact_surfaces;

			/// element list for dump and lumped_boundary
			ElementTypeMapArray<UInt> slave_elements;
			CSR<Element> node_to_elements;

			/// parallelisation
			std::unique_ptr<SynchronizerRegistry> synch_registry;
			std::unique_ptr<ElementSynchronizer> synchronizer;

			bool is_ntn_contact{true};
			};

			/// standard output stream operator
			inline std::ostream & operator<<(std::ostream & stream,
			const NTNBaseContact & _this) {
			_this.printself(stream);
			return stream;
			}

			} // namespace akantu

			/* -------------------------------------------------------------------------- */
			/* inline functions */
			/* -------------------------------------------------------------------------- */
			#include "ntn_base_contact_inline_impl.hh"

			#endif /* AST_NTN_BASE_CONTACT_HH_ */

rAKA/extra_packages/traction-at-split-node-contact/src/ntn_contact/ntn_base_contact.hhbf256d3a02cdbugfixes/all_bugs_found…

ntn_base_contact.hhView Options

rAKA/extra_packages/traction-at-split-node-contact/src/ntn_contact/ntn_base_contact.hh
bf256d3a02cdbugfixes/all_bugs_found…

ntn_base_contact.hh
View Options