Values at an arbitrary point

jam3pa

Could someone help me calculate the von Mises stress value at an arbitrary point in my analysis, which does not necessarily have to be at a node of the mesh? Therefore, interpolation would be necessary. I need to do this for a list of points.
Any assistance would be greatly appreciated.

mf

jam3pa
Hello,

in ParaVis under 'Filters' there are tools like 'Plot Over Line' and 'Point Line Interpolator'. Both give results on lines rather than single points, though. I don't know how this is done with single points or lists of points, but I am sure it can be done somehow,

Mario.

jacob

ProbeLocation in paraview interpolates value

jam3pa

Hi,
@mf and @jacob thank you very much for your help.
However, I was actually looking for an alternative that could be implemented using a function in Code_Aster, so I can program it directly in the .comm file. The reason is that I need to repeat the process multiple times.

jonasloenartz

Using MARCR_LIGN_COUPE you can get interpolated values on a line. This macro creates nodes along a line and projects a grid on these nodes using shape-function interpolation (the part about shape function interpolation is from memory, so take it with a grain of salt). Maybe this could be a solution for you as well. Using CREA_MAILLAGE with the keyword CREA_POI1 to create nodes in a second mesh, and then use PROJ_CHAMP to project the result-filed onto those nodes.

jam3pa

Hello,

@jonasloenartz thank you very much for your response. It seems to be exactly what I need. However, after reviewing the Code_Aster manuals, I'm still a bit confused. Could you please provide me with an example to better understand how to apply these functions? Specifically, I have a non-linear path of points where I want to evaluate the von Mises stresses.

Thank you for your help!

jonasloenartz

jam3pa I have a command-file somewhere, where I use this command (but on a straight line). I can share it if I find it, I will look into it later

jonasloenartz

I had a look in older models, but I couldn't find it, sorry. I hope you managed somehow.

Jonas

jam3pa

Don't worry @jonasloenartz , your help has been very useful. In the end, I performed the calculations using the PROJ_CHAMP. Since I needed the stress along a non-linear path, I calculated the stress along several straight lines defined by three points and took the stress at the last point of each line, where the last point of each line represents my desired path. It's not efficient at all, but I managed to solve it this way. Here’s the code I used:

for i in range(n):
table = MACR_LIGN_COUPE(LIGN_COUPE=_F(COOR_EXTR=(Z,x2,y2),
COOR_ORIG=(Z,X ,Y),
NB_POINTS=3,
OPERATION='EXTRACTION',
TYPE='SEGMENT'),
NOM_CHAM='SIEQ_ELNO',
RESULTAT=RE_VISU)
Esf = table.EXTR_TABLE()
P_V_Ks[i,:] = [TablaK.ABSC_CURV, Esf.VMIS[-1], TablaK.K1, TablaK.K2, TablaK.K3]