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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 99行代码中的有限元分析独立能运行吗？
121/1返回列表
查看: 1289  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

最后的火花

木虫 (正式写手)

[交流] 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行代码中的有限元分析独立能运行吗？我只想看有限元分析的结果，暂时不进行拓扑，可以做到这样吗？

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

» 猜你喜欢

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

查看全部散金贴
1楼 2018-05-31 15:32:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huvhfcu75

金虫 (小有名气)
★
最后的火花(金币+1): 谢谢参与
::::
回复此楼
4楼2018-05-31 16:08:23
已阅   回复此楼   关注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的回帖
简单回复
krxepv202楼
2018-05-31 15:59   回复  
最后的火花(金币+1): 谢谢参与
mgshr533楼
2018-05-31 16:00   回复  
最后的火花(金币+1): 谢谢参与
xuxtex5205楼
2018-05-31 16:10   回复  
最后的火花(金币+1): 谢谢参与
eeumuas1676楼
2018-05-31 16:11   回复  
最后的火花(金币+1): 谢谢参与
wehege9827楼
2018-05-31 16:13   回复  
最后的火花(金币+1): 谢谢参与
wezfr158楼
2018-05-31 16:39   回复  
最后的火花(金币+1): 谢谢参与
xmeufny0419楼
2018-05-31 16:41   回复  
最后的火花(金币+1): 谢谢参与
顶
qkvjub72710楼
2018-05-31 16:43   回复  
最后的火花(金币+1): 谢谢参与
sfvmv6411楼
2018-05-31 16:49   回复  
最后的火花(金币+1): 谢谢参与
相关版块跳转 我要订阅楼主 最后的火花 的主题更新
121/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:49 by eObJPa4gWmp6
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 j5mmowWXNWxv 2026-09-05 3/150 2026-09-06 12:49 by eObJPa4gWmp6
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 j5mmowWXNWxv 2026-09-05 4/200 2026-09-06 12:29 by eObJPa4gWmp6
[找工作] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:40 by eObJPa4gWmp6
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 7/350 2026-09-06 11:20 by eObJPa4gWmp6
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 8/400 2026-09-06 11:09 by eObJPa4gWmp6
[考研] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-04 7/350 2026-09-06 11:00 by eObJPa4gWmp6
[论文投稿] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-04 8/400 2026-09-06 11:00 by eObJPa4gWmp6
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 4/200 2026-09-06 07:20 by E1iiBLMU2XnD
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6 7ZYIEW2YWn72 2026-09-04 9/450 2026-09-06 06:00 by E1iiBLMU2XnD
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +8 CRJ0Aq2FtVH0 2026-09-03 15/750 2026-09-06 05:48 by E1iiBLMU2XnD
[硕博家园] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:39 by E1iiBLMU2XnD
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 jPDp0sc2B7zQ 2026-09-05 5/250 2026-09-06 02:19 by E1iiBLMU2XnD
[基金申请] 面上没中，邀请各位路过的虫友分析一下分数 +12 jackleilei 2026-09-03 13/650 2026-09-05 18:10 by fireguard
[文学芳草园] 两块石头 +6 阿美_Lml888 2026-09-03 8/400 2026-09-05 13:52 by 阿美_Lml888
[文学芳草园] 初秋的晨风 +4 阿美_Lml888 2026-09-04 6/300 2026-09-05 13:51 by 阿美_Lml888
[基金申请] 学科评审组评审是指会评吗？ +6 瞬息宇宙 2026-08-31 9/450 2026-09-05 13:11 by haizai99
[博后之家] 广西大学-广州大学招聘博士后 欢迎广大优秀人才！！！ +4 SCSIOyxguo 2026-09-03 7/350 2026-09-05 10:39 by Honglianlian
[硕博家园] Ei源刊怎么投 +4 Fengshun9711 2026-09-04 4/200 2026-09-04 19:30 by 研途知予
[基金申请] 要骂人了，新模版改版就是要淡化问题凝练这种虚的东西，结果有个评委还在说凝练得不够 +9 瞬息宇宙 2026-08-31 17/850 2026-09-02 14:04 by anjeeshine
信息提示
关闭
请填处理意见
关闭