Did anyone developed a code/function to create Stiffness matrix of 8-Node Brick (Hexahedral) in Finite Element Analysis ?

Question

Asif Arshid 2017-10-22

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/362654-did-anyone-developed-a-code-function-to-create-stiffness-matrix-of-8-node-brick-hexahedral-in-fini

编辑： Precise Simulation 2024-11-14，3:33

I am developing a model, which will use 8-nodded brick elements. I am looking for an already built code/function, which can create its stiffness matrix.

Thanks in advance,

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Precise Simulation 2017-10-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/362654-did-anyone-developed-a-code-function-to-create-stiffness-matrix-of-8-node-brick-hexahedral-in-fini#answer_287176

编辑：Precise Simulation 2024-11-14，3:33

在 MATLAB Online 中打开

Q1 hexahedral elements are available and pre-defined with the FEATool toolbox. Once installed, the following code assembles a bilinear diffusion matrix with 8 node brick elements on a 5x5x5 unit cube:

    grid = blockgrid( 5, 5, 5 );
    form = [2 3 4; 2 3 4];   % Bilinear form specification.
    sfun = { 'sf_hex_Q1'; 'sf_hex_Q1' };   % Q1 elements.
    coef = 1;   % Bilinear form coefficient. 
    icub = 2;   % Quadrature rule/order.
    [vrowinds,vrowcols,vavals,n_rows,n_cols] = ...
        assemblea( form, sfun, coef, icub, grid.p, grid.c, grid.a);
    help assemblea   % See help text.
    A = sparse(vrowinds,vrowcols,vavals,n_rows,n_cols);   % Create sparse matrix.
    spy( A )

To see the element/basis function definition simply open the featool/ellib/sf_hex_Q1.m file in a text editor.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Asif Arshid 2017-10-25

Thanks for sharing this resource. I downloaded and installed FEATool but I am unable to view/print the stiffness matrix of (any) single brick element, which was my objective.

请先登录，再进行评论。