Traction load on the 2-dimensional element face

[SothunN]

Hello everyone!

This is my outward normal traction vector expression: y2-y1/((x2-x1)^2+(y2-y1)^2)^2. And I want to impose this expression in ABAQUS but I have no idea to deal with it. So is there any way that can help me to solve the problem? Since I want to impose the traction load on the element face as in the sketch below.

 

[FEA way]

Like I said in previous threads, expression field editor likely won’t be sufficient for this and you may have to use scripting, subroutine or mapped field instead.

 

[SothunN]

Thank you!