Improve efficiency of AssembleFEMatrices?

7 次查看（过去 30 天）
Elliot Bontoft
Elliot Bontoft 2021-3-30
评论： Ravi Kumar 2021-3-31
Dear all,
I am using the PDE Toolbox to perform a FEM simulation. I am using the following toolbox function to aquire the matrices I am after:
FEM = assembleFEMatrices(model,'nullspace')
However, the only matrices I require are the global stiffness matrix (FEM.Kc) and nodal force vector (FEM.Fc).
The issue is that the assembleFEMatrices function is quite costly. I was wondering if there is a way improve the efficiency of this? Could the stiffness matrix and force vector be calculated another way, i.e. could they be calculated from the results of the 'solve' function that is less time consuming:
results = solve(model);
Or if it is possible to assemble ONLY the matrices you desire to reduce the computational expense?
Thank you for any advices/solutions that can be offered

请先登录，再回答此问题。

采纳的回答

Ravi Kumar
Ravi Kumar 2021-3-30
If you are using R2020b or newer version of MATLAB, you should be able to specify required matrices as input:
https://www.mathworks.com/help/pde/ug/assemblefematrices.html
Regards,
Ravi
  3 个评论
显示 1更早的评论
Elliot Bontoft
Elliot Bontoft 2021-3-31
Hi all,
Ravi was right and this was possible through using specified requested matrices in the input. You just have to manually impose the boundary conditions to build the stiffness matrix and force vector, as follows:
FEM = assembleFEMatrices(model,'HKAQFG');
[FEM.B,Or] = pdenullorth(FEM.H); % compute the nullspace of columns of H
FEM_Fc = FEM.B'*(FEM.F + FEM.G);
FEM_Kc = FEM.B'*(FEM.K + FEM.A + FEM.Q)*FEM.B;
This is notably faster than using
FEM = assembleFEMatrices(model,'nullspace')
to access the FEM.Fc and FEM.Kc matrices
Ravi Kumar
Ravi Kumar 2021-3-31
Eaxctly! Selecting a subset of matrices or specifying BC appliction option are sort of orthogonal. You need to do the way you are doing. I would also add, you can assemble some of these matrices once and some that have some variation over time or solution repeatedly.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 General PDEs 的更多信息

标签

产品


版本

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by