CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > General Forums > Main CFD Forum

How to use shape function to get the displacement of area from vertex

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree1Likes
  • 1 Post By FMDenaro

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 26, 2022, 13:01
Default How to use shape function to get the displacement of area from vertex
  #1
New Member
 
J Powell
Join Date: Feb 2022
Posts: 12
Rep Power: 4
miraboreasu is on a distinguished road
Hello,

Sorry for another thread.
There is a triangle in the space, I have the displacement (vector) of each vertex, like

d1: -0.000229128 0.000131475 -0.000706995
d2: -0.000252541 6.02E-05 -0.000662299
d3: -0.000270627 0.000147876 -0.000661116

1. How to calculate the displacement in the normal direction of the whole triangle?
2. Is something called shape function related to this, and how?
miraboreasu is offline   Reply With Quote

Old   September 26, 2022, 21:32
Default
  #2
Senior Member
 
Sayan Bhattacharjee
Join Date: Mar 2020
Posts: 495
Rep Power: 8
aerosayan is on a distinguished road
By displacement, do you mean, you want to move/translate the triangle in the direction of its normal?

If yes, then convert the normal into a unit vector, scale it up by the distance you want to move everything, and apply the translation to each point of the triangle.

https://www.gamedev.net/forums/topic...ector/4734660/

I'm new to FEM, but shape functions are used to approximate a solution in a cell. They're polynomial and denote the variation of a property inside each cell's domain.

I'm not an expert, by any means, so I don't know if they're applicable here.
aerosayan is offline   Reply With Quote

Old   September 27, 2022, 05:31
Default
  #3
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,887
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Shape functions are used to define a functional variation in the triangle (or other figures). If you have a triangle (1,2,3) then you can define a linear polynomial:


f(x,y) = c0+c1*x+c2*y


that can be rearranged in terms of the values in the vertices


f(x,y)=N1(x,y)*f1+N2(x,y)*f2+N3(x,y)*f3


where the Nj are the linear shape functions.


As you can see, you can get both derivatives and integrals by simple analytical integration.
aerosayan likes this.
FMDenaro is offline   Reply With Quote

Old   September 28, 2022, 11:21
Default
  #4
New Member
 
J Powell
Join Date: Feb 2022
Posts: 12
Rep Power: 4
miraboreasu is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Shape functions are used to define a functional variation in the triangle (or other figures). If you have a triangle (1,2,3) then you can define a linear polynomial:


f(x,y) = c0+c1*x+c2*y


that can be rearranged in terms of the values in the vertices


f(x,y)=N1(x,y)*f1+N2(x,y)*f2+N3(x,y)*f3


where the Nj are the linear shape functions.


As you can see, you can get both derivatives and integrals by simple analytical integration.


So here is my vertices displacement,

d1: -0.00023 0.000131 -0.00071
d2: -0.00025 6.02E-05 -0.00066
d3: -0.00027 0.000148 -0.00066

I already have the pressure,p, and area, a, of this area.

is d1 the f1? Sorry I didn't understand.
miraboreasu is offline   Reply With Quote

Old   September 28, 2022, 11:39
Default
  #5
Senior Member
 
Lucky
Join Date: Apr 2011
Location: Orlando, FL USA
Posts: 5,756
Rep Power: 66
LuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura aboutLuckyTran has a spectacular aura about
f1 f2 and f3 are the values of f at the vertices 1, 2, and 3.

When you do FEM, at the very beginning, before you do anything else, you do two things:
1) you create a computational grid
2) you (not me, you!) decide a basis for the shape functions, e.g. piecewise linear

It's very weird to say you are doing FEM and then ask somehow else how to evaluate the integral because you must have defined the integral before you can say you are doing FEM. All we can do is state the general principle. Or you must tell us exactly what is your FEM approach and we can help show some of the mechanical steps.

Work comes from integrating the mechanical pressure over the displacement, the integral being determined by your FEM approach. So how did you define the FEM integral? I hope you see the issue.


Unless.. you are not doing FEM at all but something else entirely...
LuckyTran is offline   Reply With Quote

Old   September 28, 2022, 14:04
Default
  #6
New Member
 
J Powell
Join Date: Feb 2022
Posts: 12
Rep Power: 4
miraboreasu is on a distinguished road
Quote:
Originally Posted by LuckyTran View Post
f1 f2 and f3 are the values of f at the vertices 1, 2, and 3.

When you do FEM, at the very beginning, before you do anything else, you do two things:
1) you create a computational grid
2) you (not me, you!) decide a basis for the shape functions, e.g. piecewise linear

It's very weird to say you are doing FEM and then ask somehow else how to evaluate the integral because you must have defined the integral before you can say you are doing FEM. All we can do is state the general principle. Or you must tell us exactly what is your FEM approach and we can help show some of the mechanical steps.

Work comes from integrating the mechanical pressure over the displacement, the integral being determined by your FEM approach. So how did you define the FEM integral? I hope you see the issue.


Unless.. you are not doing FEM at all but something else entirely...

I am post-processing results from a combined finite–discrete element method software. So I don't touch on the detail of the method. This software output the displacement of vertices as
d1: -0.00023 0.000131 -0.00071
d2: -0.00025 6.02E-05 -0.00066
d3: -0.00027 0.000148 -0.00066

I know the pressure I applied to the surface meshed by many triangles. And I can calculate the area of each triangle from the nodal coordinates, but I don't know what to do next.

I want to use pressure*area*displacement, but I don't know what displacement I should use.

Can you give me some advice?
miraboreasu is offline   Reply With Quote

Reply

Tags
fem

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[mesh manipulation] refineMesh Error mohsen.boojari OpenFOAM Meshing & Mesh Conversion 3 March 1, 2018 23:07
ParaView for OF-1.6-ext Chrisi1984 OpenFOAM Installation 0 December 31, 2010 07:42
Version 15 on Mac OS X gschaider OpenFOAM Installation 113 December 2, 2009 11:23
CFX Solver Memory Error mike CFX 1 March 19, 2008 08:22
Droplet Evaporation Christian Main CFD Forum 2 February 27, 2007 07:27


All times are GMT -4. The time now is 11:32.