User login

You are here

A question about RPIM shapefunction

I am solving a simple solid mechanics problem by meshless method.  use RPIM to  calculate shape function.

when i calculate the derivatives of shape function, i found the derivative is not close to zero at the compute point, when the point is on the boundry of problem domain. 

 

My question is, should the derivative be zero if the node is on the boundry?

derivatives of shape function respect to x, when the point is in domain

fig.1 derivatives of shape function respect to x when the point for computing is in domain.

 the derivative at the comput point is very close to zero.derivatives of shape function respect to y when the point for computing is on lower boundry of domain.

 Fig.2 derivatives of shape function respect to Y when the point for computing is on the lower boundry of domain.

   the derivative at the comput point is not close to zero.

derivatives of shape function respect to x when the point for computing is at lower left corner of domain.

 also, the derivative at the comput point is much bigger than zero.

You are right, if the node is inside the domain then the shape function
derivative is zero at the node but if the node is on the boundary then
the derivative of the shape function is not zero at the node. If you
plot the shape function and look at its slope you will also see this and
the reason is that the shape function is symmetric about a node when it
is inside the domain but this is not true for a node on the boundary.

 

 Thank you zahoorswati!

 

I see your point. The shape function derivative is discontinuous on the normal direction of boundary.

My compute result shows that, I can get the right displacement and stress inside the domain, but can't get them on the boundary. So I thought may be the shape function has some problem.

I will check on my code. 

why can't I submit reply.

Subscribe to Comments for "A question about RPIM shapefunction"

Recent comments

More comments

Syndicate

Subscribe to Syndicate