BilinearTraction

Source: models/solid_mechanics/traction_separation_law/BilinearTraction.py

Assemble bilinear cohesive zone traction from an externally updated damage: \(T_n = K(1-d)\delta_n^+ + K\delta_n^-,\;T_{si} = K(1-d)\delta_{si}\).

Inputs

damage — input · Scalar · required

Current damage

normal_separation — input · Scalar · required

Normal separation (typically the Macaulay-positive part of the normal jump)

normal_penetration — input · Scalar

Optional normal penetration (typically the Macaulay-negative part of the normal jump). When supplied, K times this is added to the normal traction as a penalty term resisting interpenetration.

tangential_separation_1 — input · Scalar · required

First tangential separation

tangential_separation_2 — input · Scalar · required

Second tangential separation

Outputs

traction — output · Vec · required

Traction vector

Parameters

penalty_stiffness — parameter · Scalar · required

Penalty stiffness