Nonlinear Imperfect Interface - Revision history

Nairnj: /* Interface Law Properties */

2026-03-28T00:51:58Z

Interface Law Properties

← Older revision		Revision as of 00:51, 28 March 2026
Line 108:		Line 108:
	: Return derivative of the force due to the traction law at the provided displacement discontinuity <tt>delta</tt>. The derivative must be in [[ConsistentUnits Command#Legacy and Consistent Units\|Force units/Length units]] for given displacement in [[ConsistentUnits Command#Legacy and Consistent Units\|Length units]]. Imperfect interface laws give a traction or stress derviative per unit discontinuity. To convert such a law to force derivative units, multiply the traction derivative by the contact area (provided in the <tt>area</tt>).		: Return derivative of the force due to the traction law at the provided displacement discontinuity <tt>delta</tt>. The derivative must be in [[ConsistentUnits Command#Legacy and Consistent Units\|Force units/Length units]] for given displacement in [[ConsistentUnits Command#Legacy and Consistent Units\|Length units]]. Imperfect interface laws give a traction or stress derviative per unit discontinuity. To convert such a law to force derivative units, multiply the traction derivative by the contact area (provided in the <tt>area</tt>).

	; <tt>double GetFnEnergy(double delta,double area) const</tt><br><tt>double GetFtPrime(double delta,double tr) const</tt>		; <tt>double GetFnEnergy(double delta,double area) const</tt><br><tt>double GetFtPrime(double delta,double area) const</tt>
	: Return area under the traction law up to given displacement discontinuity in <tt>delta</tt> (in [[ConsistentUnits Command#Legacy and Consistent Units\|Length units]]). The units of the return value should be [[ConsistentUnits Command#Legacy and Consistent Units\|Energy units]]. ~~Many~~ imperfect interface laws give a traction or stress ~~which have derivative for stress per unit length~~. To convert such a law ~~proper force derivative units~~, multiply ~~the traction derivative~~ by the contact area (provided in the <tt>area</tt>. Note interfacial energy can be archived in the [[MPM Global Archiving Options\|global archiving interface energy option]], but has no other affect on calculations. It does not contribute to dissipated energy and does not cause any temperature changes.		: Return area under the traction law up to given displacement discontinuity in <tt>delta</tt> (in [[ConsistentUnits Command#Legacy and Consistent Units\|Length units]]). The units of the return value should be [[ConsistentUnits Command#Legacy and Consistent Units\|Energy units]]. Most imperfect interface laws give a traction or stress. To convert integration of such a law to energy, multiply by the contact area (provided in the <tt>area</tt>). Note: interfacial energy can be archived in the [[MPM Global Archiving Options\|global archiving interface energy option]], but has no other affect on calculations. It does not contribute to dissipated energy and does not cause any temperature changes.

	== Examples ==		== Examples ==

Nairnj: /* Interface Law Properties */

2026-03-28T00:49:38Z

Interface Law Properties

← Older revision		Revision as of 00:49, 28 March 2026
Line 106:		Line 106:

	; <tt>double GetFnPrime(double delta,double area) const</tt><br><tt>double GetFtPrime(double delta,double area) const</tt>		; <tt>double GetFnPrime(double delta,double area) const</tt><br><tt>double GetFtPrime(double delta,double area) const</tt>
	: Return derivative of the force due to the traction law at the provided displacement discontinuity <tt>delta</tt>. The derivative must be in [[ConsistentUnits Command#Legacy and Consistent Units\|Force units/Length units]] for given displacement in [[ConsistentUnits Command#Legacy and Consistent Units\|Length units]]. Imperfect interface laws give a traction or stress per unit discontinuity. To convert such a law to force derivative units, multiply the traction derivative by the contact area (provided in the <tt>area</tt>).		: Return derivative of the force due to the traction law at the provided displacement discontinuity <tt>delta</tt>. The derivative must be in [[ConsistentUnits Command#Legacy and Consistent Units\|Force units/Length units]] for given displacement in [[ConsistentUnits Command#Legacy and Consistent Units\|Length units]]. Imperfect interface laws give a traction or stress derviative per unit discontinuity. To convert such a law to force derivative units, multiply the traction derivative by the contact area (provided in the <tt>area</tt>).

	; <tt>double GetFnEnergy(double delta,double area) const</tt><br><tt>double GetFtPrime(double delta,double tr) const</tt>		; <tt>double GetFnEnergy(double delta,double area) const</tt><br><tt>double GetFtPrime(double delta,double tr) const</tt>

Nairnj: /* Interface Law Properties */

2026-03-28T00:48:44Z

Interface Law Properties

← Older revision		Revision as of 00:48, 28 March 2026
Line 106:		Line 106:

	; <tt>double GetFnPrime(double delta,double area) const</tt><br><tt>double GetFtPrime(double delta,double area) const</tt>		; <tt>double GetFnPrime(double delta,double area) const</tt><br><tt>double GetFtPrime(double delta,double area) const</tt>
	: Return derivative of the force due to the traction law at the provided displacement discontinuity <tt>delta</tt>. The derivative must be in [[ConsistentUnits Command#Legacy and Consistent Units\|Force units/Length units]] for given displacement in [[ConsistentUnits Command#Legacy and Consistent Units\|Length units]]. ~~Many imperfect~~ interface laws give a traction or ~~stress which have derivative for~~ stress per unit ~~length~~. To convert such a law ~~proper~~ force derivative units, multiply the traction derivative by the contact area (provided in the <tt>area</tt>).		: Return derivative of the force due to the traction law at the provided displacement discontinuity <tt>delta</tt>. The derivative must be in [[ConsistentUnits Command#Legacy and Consistent Units\|Force units/Length units]] for given displacement in [[ConsistentUnits Command#Legacy and Consistent Units\|Length units]]. Imperfect interface laws give a traction or stress per unit discontinuity. To convert such a law to force derivative units, multiply the traction derivative by the contact area (provided in the <tt>area</tt>).

	; <tt>double GetFnEnergy(double delta,double area) const</tt><br><tt>double GetFtPrime(double delta,double tr) const</tt>		; <tt>double GetFnEnergy(double delta,double area) const</tt><br><tt>double GetFtPrime(double delta,double tr) const</tt>

Nairnj: /* Interface Law Properties */

2026-03-28T00:45:51Z

Interface Law Properties

← Older revision		Revision as of 00:45, 28 March 2026
Line 100:		Line 100:
	==== Interface Law Properties ====		==== Interface Law Properties ====

	Interface law calculations are done in the following methods. Each one has a <tt>switch(style)</tt> block to add any implemented law. When adding a new law, insert a <tt>case</tt> block (with code defined [[#Define Shape Code\|above]]) to return to following values. If you are using [[ConsistentUnits Command\|legacy units, be sure to return values using the '''mm-g-s''' units.		Interface law calculations are done in the following methods. Each one has a <tt>switch(style)</tt> block to add any implemented law. When adding a new law, insert a <tt>case</tt> block (with code defined [[#Define Shape Code\|above]]) to return to following values. If you are using [[ConsistentUnits Command\|legacy units]], be sure to return values using the '''mm-g-s''' units.

	; <tt>double GetFn(double delta,double area) const</tt><br><tt>double GetFt(double delta,double area) const</tt>		; <tt>double GetFn(double delta,double area) const</tt><br><tt>double GetFt(double delta,double area) const</tt>

Nairnj: /* Add Interface Law Parameters */

2026-03-28T00:44:18Z

Add Interface Law Parameters

← Older revision		Revision as of 00:44, 28 March 2026
Line 83:		Line 83:
	to check entered values. If they are not allowed, return a string that explains the error.		to check entered values. If they are not allowed, return a string that explains the error.

	Finally, to record use of the new properties in output ~~file~~, add code to		Finally, to record use of the new properties in simulation output, add code to

	void NonlinearInterface::PrintContactLaw(void) const		void NonlinearInterface::PrintContactLaw(void) const

Nairnj: /* Add Interface Law Parameters */

2026-03-28T00:43:55Z

Add Interface Law Parameters

← Older revision		Revision as of 00:43, 28 March 2026
Line 83:		Line 83:
	to check entered values. If they are not allowed, return a string that explains the error.		to check entered values. If they are not allowed, return a string that explains the error.

	Finally, to record use of the new properties, add code to		Finally, to record use of the new properties in output file, add code to

	void NonlinearInterface::PrintContactLaw(void) const		void NonlinearInterface::PrintContactLaw(void) const

Nairnj: /* Add Interface Law Parameters */

2026-03-28T00:43:17Z

Add Interface Law Parameters

← Older revision		Revision as of 00:43, 28 March 2026
Line 81:		Line 81:
	const char *NonlinearInterface::VerifyAndLoadProperties(int np)		const char *NonlinearInterface::VerifyAndLoadProperties(int np)

	to ~~checked~~ entered values. If ~~enter values~~ are not ~~allows~~ return a string that explains the error.		to check entered values. If they are not allowed, return a string that explains the error.

	Finally, to record use of the new properties, add code to		Finally, to record use of the new properties, add code to

Nairnj: /* Add Interface Law Parameters */

2026-03-28T00:42:47Z

Add Interface Law Parameters

← Older revision		Revision as of 00:42, 28 March 2026
Line 77:		Line 77:
	If you create new material properties, they should be added to the <tt>NairnMPM.dtd</tt> file. The best approach is to search for an existing property (like <tt>Npeak</tt>) and insert new property with new name (usually needed two place in the DTD file).		If you create new material properties, they should be added to the <tt>NairnMPM.dtd</tt> file. The best approach is to search for an existing property (like <tt>Npeak</tt>) and insert new property with new name (usually needed two place in the DTD file).

	If any new properties ~~are restrictions~~ on allowed values, add code to		If any new properties haverestrictions on allowed values, add code to

	const char *NonlinearInterface::VerifyAndLoadProperties(int np)		const char *NonlinearInterface::VerifyAndLoadProperties(int np)

Nairnj: /* Define Shape Code */

2026-03-28T00:41:38Z

Define Shape Code

← Older revision		Revision as of 00:41, 28 March 2026
Line 61:		Line 61:
	=== Define Shape Code ===		=== Define Shape Code ===

	In <tt>NonlinearInterface.hpp</tt> define a new code for the new shape being implemented~~. In other word, add~~ a new constant to the enumerated list		In <tt>NonlinearInterface.hpp</tt> define a new code for the new shape being implemented by adding a new constant to the enumerated list

	enum { NL_LINEAR_INTERFACE=0,MORSE_POTENTIAL,MAX_STYLES };		enum { NL_LINEAR_INTERFACE=0,MORSE_POTENTIAL,MAX_STYLES };

Nairnj: /* Define Shape Code */

2026-03-28T00:31:53Z

Define Shape Code

← Older revision		Revision as of 00:31, 28 March 2026
Line 61:		Line 61:
	=== Define Shape Code ===		=== Define Shape Code ===

	In <tt>NonlinearInterface.hpp</tt> define a new code for the new shape being ~~implement~~. In other word, add a new constant to the enumerated list		In <tt>NonlinearInterface.hpp</tt> define a new code for the new shape being implemented. In other word, add a new constant to the enumerated list

	enum { NL_LINEAR_INTERFACE=0,MORSE_POTENTIAL,MAX_STYLES };		enum { NL_LINEAR_INTERFACE=0,MORSE_POTENTIAL,MAX_STYLES };