24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1086  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

最后的火花

木虫 (正式写手)


[交流] 99行代码中的有限元分析独立能运行吗?

有限元求解子程序
1 function [U]=FE(nelx,nely,x,penal)
2 [KE] = lk;                %单元刚度矩阵
3 K = sparse(2*(nelx+1)*(nely+1), 2*(nelx+1)*(nely+1));   %总体刚度矩阵的稀疏矩阵
4 F = sparse(2*(nely+1)*(nelx+1),1); U = zeros(2*(nely+1)*(nelx+1),1); %力矩阵的稀疏矩阵
5 for elx = 1:nelx
6   for ely = 1:nely
7     n1 = (nely+1)*(elx-1)+ely;
8     n2 = (nely+1)* elx   +ely;              
9     edof = [2*n1-1; 2*n1; 2*n2-1; 2*n2; 2*n2+1; 2*n2+2; 2*n1+1; 2*n1+2];   %这里的Y轴是反向的,但是不影响最后的结果,详情请见二楼TYNGOD这位高手的解释,感谢TYNGOD。
10    K(edof,edof) = K(edof,edof) + x(ely,elx)^penal*KE;   %将单元刚度矩阵组装成总的刚度矩阵
11  end
12end

问题1:99行代码中有限元求解子程序最终求得的结果是第10行的K,总刚度矩阵吗?

问题2:99行代码中的有限元分析独立能运行吗?我只想看有限元分析的结果,暂时不进行拓扑,可以做到这样吗?

还有朋友在研究这个程序吗?欢迎私信我再继续交流
回复此楼

» 猜你喜欢

» 抢金币啦!回帖就可以得到:

查看全部散金贴

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ShiningEast

铁虫 (初入文坛)


★ ★
小木虫: 金币+0.5, 给个红包,谢谢回帖
jjdg: 金币+1, 感谢参与 2018-06-13 17:14:42
99行中的K是经过拓扑处理后的总刚度矩阵,在程序第10行中,把x中所有的值均设为1,得到的就是正常的总体刚度矩阵了。
99行中的有限元程序是可以独立运行的。
12楼2018-06-13 16:50:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答
普通表情 高级回复(可上传附件)
信息提示
请填处理意见