Learning Forum

[lschuler]

Hello,

I am trying to define my own finite element type using the subroutine USERELEM. I can see in the example written in the subroutine that the function ElemShpFn is called to obtain the shape functions, their derivatives and other information. What are the input arguments of this function? Can I use it to obtain shape functions of other elements that the ones used in the example?

Thank you

[Surya Prakash]

Can you share the example that you are referring to as an image?

Regards Surya

[lschuler]

Here is the part of the code corresponding to the call to the shape functions. I am guessing that the first input is the order of the element (1 linear, 2 quadratic) because KeyOpt(1)=0 corresponds to linear quads elements and KeyOpt(1)=1 to quadratic cubes elements. The two following inputs are the integration point number, the desired output (1 for the shape functions, 2 for their derivatives, 3 for the integration weight). Is that right?
Thank you Louis