interpolation strain tensor field
显示 更早的评论
I have a 3D+t strain tensor field (3D over time) delivered by FAM software (Abaqus). What I need is to interpolate the strain tensor on a fixed plane to gest a 2D tensor maps over time. Can I use the intrep3 function for that interpolating each component of a tensor? Is there any better solution for this problem. I have read that it is better to interpolate eigenvalues and eigenvectors but I did not found details about how to do this
采纳的回答
Nicola Bombace
2018-10-11
1 个投票
You can interpolate the strain components easily with the interp3 function, and the method linear.
Is it an elastic material? If so you could reconstruct the strain tensor by interpolating only the displacement and using the shape functions...
Which kind of elements are you using?
8 个评论
BLP
2018-10-11
Nicola Bombace
2018-10-12
编辑：Nicola Bombace
2018-10-12
Hi. In this case, there is a stream of literature that considers the interpolation of second order tensors with eigenvalues and eigenvectors (for example here and here ). The functions you need to add to your workflow for such methodology is eig. To use it, you will have to construct the strain tensor in Matlab as a matrix. The interp3 function with linear option eig uses functions with the same order of continuity of your elements. Personally I have always used interpolation directly from the components, since is a very general approach. Hope this helps.
Stephen23
2018-10-12
Hi Nicola. Thank you for your answer. I have found the first paper, you suggested, but I did not find the second one. I am not very fluent in tensor calculus, and the second paper semms as a remedy for this? However, I will be grateful for explaining one thing to me. I am not sure what do you mean by constructing tensor as a matrix. Do you mean a 3D table?
Nicola Bombace
2018-10-12
Hi Beata,
the second paper is
"Kratz, Andrea, et al. "Visualization and analysis of second‐order tensors: Moving beyond the symmetric positive‐definite case." Computer Graphics Forum. Vol. 32. No. 1. Oxford, UK: Blackwell Publishing Ltd, 2013."
For the strain tensor what I meant is that Abaqus gives you the stress in Voigt notation while you want to transfer them in matrix notation as explained in this Wikipedia page. Let's suppose that your strains, at one moment in time, are saved as
[eps1, eps2, eps3 ,eps4, eps5, eps6]; %Voigt notation
E = [ eps1 eps6*0.5 eps5*0.5;
eps6*0.5 eps2 eps4*0.5;
eps5*0.5 eps4*0.5 eps3]; % Matrix notation
eigenVals = eig(E); % Vector with eigenvalues
From here on you can continue with any algorithm that uses eigenvalues of the strain tensor.
BLP
2018-10-12
Nicola Bombace
2018-10-12
编辑：Nicola Bombace
2018-10-12
If I understand correctly you are using Abaqus Explicit. In this case, you should really use C3D8R, for possible volumetric locking of the elements. Using this element would give you three advantages:
- Robust simulation, because of no locking.
- Smaller Computational time, because of the reduced integration
- You can get the strain directly from the integration point of the element (is an output of Abaqus), and this integration point is in the centre of the hexahedron.
Ps You are really welcome, I hope all this information helps :)
BLP
2018-10-12
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Numerical Integration and Differential Equations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
