I am looking in the development of the Finite Elasticity model and the applied linearization (for the Newton's method), and I wish to ask if someone can confirm or correct my understanding of the code.
In the definition of the LocalMatrix at file ElasticSolve.src, the force vestor is defined in the following way:
May I assume that the term DOT_PRODUCT(dBasisdx(p,:),Stress1(i,:)) presents the internal force vector? And the term DOT_PRODUCT(dBasisdx(p,:),dStress1U(i,:)) presents the vector which is result of multiplication of the jacobian and the dispalcement vector?ForceVector(dim*(p-1)+i) = ForceVector(dim*(p-1)+i) &
+(Basis(p)*Force(i)*DetDefG &
+Basis(p)*InertialForce(i)*Density &
-DOT_PRODUCT(dBasisdx(p,:),Stress1(i,:)) &
+DOT_PRODUCT(dBasisdx(p,:),dStress1U(i,:)))*s
The stiffness matrix is defined in this way:
Does it mean that DOT_PRODUCT(dBasisdx(q,:),dStress1(j,:)) is the jacobian of the internal force vector?StiffMatrix(dim*(p-1)+i,dim*(q-1)+j) &
= StiffMatrix(dim*(p-1)+i,dim*(q-1)+j) &
+ DOT_PRODUCT(dBasisdx(q,:),dStress1(j,:))*s
Regards,
Stan