版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (4293)
>学术会议 (145)
>导师招生 (138)
>论文投稿 (121)
>休闲灌水 (110)
>文献求助 (99)
>虫友互识 (99)
>博后之家 (78)
>考博 (48)
>硕博家园 (47)
>基金申请 (46)
>找工作 (27)
>考研 (26)
>教师之家 (24)
>公派出国 (24)
>论文道贺祈福 (18)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

最后的火花

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 3341
帖子: 685
在线: 114.5小时
虫号: 2679796

[交流] 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文章T0P，我:8O.551.O54,科目全,可十急已经有8人回复
两块石头已经有9人回复
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有6人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有9人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有8人回复
广西大学-广州大学招聘博士后欢迎广大优秀人才！！！已经有8人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有3人回复
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有3人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有4人回复
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有3人回复

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

查看全部散金贴

1楼 2018-05-31 15:32:43

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

ShiningEast

铁虫 (初入文坛)

应助: 1 (幼儿园)
金币: 13.3
帖子: 4
在线: 2小时
虫号: 3557255

★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
jjdg: 金币+1, 感谢参与 2018-06-13 17:14:42

99行中的K是经过拓扑处理后的总刚度矩阵，在程序第10行中，把x中所有的值均设为1，得到的就是正常的总体刚度矩阵了。
99行中的有限元程序是可以独立运行的。

赞一下