| 查看: 1288 | 回复: 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行代码中的有限元分析独立能运行吗？我只想看有限元分析的结果，暂时不进行拓扑，可以做到这样吗？ 还有朋友在研究这个程序吗？欢迎私信我再继续交流 |
回复此楼» 猜你喜欢
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
科研人应该花精力去思考如何解决问题，而不是去凝练问题
已经有18人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有15人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
面上没中，邀请各位路过的虫友分析一下分数
已经有13人回复
-
两块石头
已经有8人回复
-
初秋的晨风
已经有6人回复
-
学科评审组评审是指会评吗？
已经有9人回复
:: 4楼2018-05-31 16:08:23
12楼2018-06-13 16:50:45
简单回复
krxepv202楼
2018-05-31 15:59
回复
最后的火花(金币+1): 谢谢参与
mgshr533楼
2018-05-31 16:00
回复
最后的火花(金币+1): 谢谢参与
2018-05-31 16:10
回复
最后的火花(金币+1): 谢谢参与
2018-05-31 16:11
回复
最后的火花(金币+1): 谢谢参与
2018-05-31 16:13
回复
最后的火花(金币+1): 谢谢参与
wezfr158楼
2018-05-31 16:39
回复
最后的火花(金币+1): 谢谢参与
2018-05-31 16:41
回复
最后的火花(金币+1): 谢谢参与
顶
qkvjub72710楼
2018-05-31 16:43
回复
最后的火花(金币+1): 谢谢参与
sfvmv6411楼
2018-05-31 16:49
回复
最后的火花(金币+1): 谢谢参与