×
INTELLIGENT WORK FORUMS
FOR ENGINEERING PROFESSIONALS

Log In

Come Join Us!

Are you an
Engineering professional?
Join Eng-Tips Forums!
  • Talk With Other Members
  • Be Notified Of Responses
    To Your Posts
  • Keyword Search
  • One-Click Access To Your
    Favorite Forums
  • Automated Signatures
    On Your Posts
  • Best Of All, It's Free!
  • Students Click Here

*Eng-Tips's functionality depends on members receiving e-mail. By joining you are opting in to receive e-mail.

Posting Guidelines

Promoting, selling, recruiting, coursework and thesis posting is forbidden.

Students Click Here

Shape function gradients with 3D boudary elements

Shape function gradients with 3D boudary elements

Shape function gradients with 3D boudary elements

(OP)
Hello all,

I am developing a 3D boundary element program where I have only surface elements (4-node quads and 3-node tris), but that are oriented in 3D space. Likewise, I want to compute x,y,z gradients on each of these faces.

I have a nodal solution (phi), and I want to compute dphi/dx, dphi/dy, dphi/dz.

I have successfully done this for quads, but cannot seem to figure it out for triangles.

For the quad-4 with linear shape functions:
N = [1/4*(1+xi)*(1+eta) 1/4(1+xi)*(1-eta) 1/4(1-xi)*(1-eta) 1/4(1-xi)*(1-eta)]
I can compute the differential operators in x,y,z (dNdx, dNdy, and dNdz) by simply computing 2 tangential vectors (dxdxi, dxdeta, dydxi, dydeta, dzdxi, dzdeta):
e1 = dNdxi*[ex,ey,ez]
and
e2 = dNdeta*[ex,ey,ez]

Taking the cross product of these gives me my surface gradient:
e3 = cross(e1,e2)

Then using these 3 vectors as my jacobian:
Jac = [e1; e2; e3];

I can invert my jacobian to compute dxidx,dxidy,dxidz, detadx,detady,detadz, and finally post multiply with dNdxi and dNdeta again to compute dNdx, dNdy, dNdz:

dNdx = Jac(1,1)^-1*dNdxi + Jac(1,2)^-1*dNdeta
dNdy = Jac(2,1)^-1*dNdxi + Jac(2,2)^-1*dNdeta
dNdz = Jac(3,1)^-1*dNdxi + Jac(3,2)^-1*dNdeta

These work well, and I get the correct gradients.

But now for the 3-node triangle, it's not so straight forward. I have formulated the TRI-3 with 3 local coordinates (xi, eta, and zeta), where my shape functions are just:
N = [xi eta zeta],
and so:
dN = I (identity).

How can I compute dNdx, dNdy, and dNdz for the linear 3-node triangle?

Any help is appreciated, and thank you in advance.

RE: Shape function gradients with 3D boudary elements

Being not so familiar with Gauss Quadrature in the last 10 years, I found this which also has some directions for a triangle. Not sure if it is explained there but it feels like you are entering a world of pain :) Good luck!
http://math2.uncc.edu/~shaodeng/TEACHING/math5172/...

Please let me know if you find the solution later, I would be very interested as well..

Spaceship!!


Aerospace Engineer, M.Sc. / Aircraft Stress Engineer with 7 years of experience
(United States)

RE: Shape function gradients with 3D boudary elements

(OP)
hi aerostress82, the issue I am talking about is not Guass Quadrature. Although, if I were to integrate my gradient (which I am talking about) over the element, I would use Gauss Quadrature.

The issue I am referring to is computing the discrete exterior gradient on the element face (in 3D, 3 components). Since it's oriented in 3D space, standard 2D formulations do not apply, and since it's a boundary element, standard 3D volume formulations also do not apply.

Sort of tricky, and cannot seem to get the correct derivation, and cannot find anything online about it; though I know somebody has at least attempted it before.

RE: Shape function gradients with 3D boudary elements

Thanks for the clarification jfieldstone. I'm not familiar with boundary element method but the quad-4 equation looked very similar to finite element - that's why I sent that link.
Maybe if you happen to have some research assistant friend(s) in university, you could send him/her the sciencedirect links related to your case? (after checking them from sciencedirect)

Other than that, hope someone can help you on this. Good luck with it!

Spaceship!!
Aerospace Engineer, M.Sc. / Aircraft Stress Engineer

Red Flag This Post

Please let us know here why this post is inappropriate. Reasons such as off-topic, duplicates, flames, illegal, vulgar, or students posting their homework.

Red Flag Submitted

Thank you for helping keep Eng-Tips Forums free from inappropriate posts.
The Eng-Tips staff will check this out and take appropriate action.

Reply To This Thread

Posting in the Eng-Tips forums is a member-only feature.

Click Here to join Eng-Tips and talk with other members! Already a Member? Login


Resources

White Paper - Implementing a Multi-Domain System
IoT systems are multi-domain designs that often require AMS, Digital, RF, photonics and MEMS elements within the system. Tanner EDA provides an integrated, top-down design flow for IoT design that supports all these design domains. Learn more about key solutions that the Tanner design flow offers for successful IoT system design and verification. Download Now

Close Box

Join Eng-Tips® Today!

Join your peers on the Internet's largest technical engineering professional community.
It's easy to join and it's free.

Here's Why Members Love Eng-Tips Forums:

Register now while it's still free!

Already a member? Close this window and log in.

Join Us             Close