指定材料的平面应变和平面应力有限元fortran程序
第一次在这里发新贴，请多支持。该程序使用三节点单元(triangular element)mesh一个左端与墙固接的梁，上表面均布载荷为QP。
求解了弹性矩阵D，local stiffness matrix (SLM)，global stiffness matrix (GSM)等
该程序在linux下用fortran77编译通过，运行正常，给刚入门的朋友们借鉴一下。
这个程序不是我写的，是我老板以前写的，不过我正在学习写4节点和8节点的fortran程序，但愿能成功。
4楼至6楼是4节点单元的fortran程序
如果有问题或建议，欢迎跟贴讨论。
[ 本帖最后由 eldfxm 于 2008-5-14 10:47 编辑 ]
回复 5# 的帖子
续楼上。附件是从第１页到第１５页的全部内容。
其中第２页到第１３页是推导。
第１４页是节点与单元关联的子程序。
第１５页是代码流程图。
[ 本帖最后由 eldfxm 于 2008-5-21 08:44 编辑 ] 为什么是linux的环境呢，之前有试着用linux，结果连个无线上网的驱动都没装好。。。 原帖由 luzheng111 于 2008-5-10 13:35 发表 http://www.simwe.com/forum/images/common/back.gif
为什么是linux的环境呢，之前有试着用linux，结果连个无线上网的驱动都没装好。。。
因为我这里Windows只能买正版的，所以就用免费的linux了。linux下的fortran77和windows下的fortran90的compiler不同，但是代码的基本思路还是相同的。这个代码只作为参考，要在windows下编译的话需要用fortran90改写。
我本人也很希望能在windows下使用fortran，我在等能同时支持Windows Vista和Visual Studio 2008的Inter Fortran Compiler，等出来后再改写一下。 4节点单元的fortran程序写完了。仍然是指定材料的平面应力和平面应变的情况。
先上张图作为题目描述。前面3节点单元求解的问题和这个一样。
[ 本帖最后由 eldfxm 于 2008-5-14 07:01 编辑 ]
回复 4# 的帖子
续楼上。接下来描述的是４节点单元的思路，包括有限元公式推导
节点坐标矩阵
节点与单元关联矩阵
弹性矩阵
单元刚度矩阵
高斯数值积分
总体刚度矩阵
边值问题：载荷和边界条件
求解子程序
显示节点位移和节点应力
先上第２页作为样图，共有１５页。
回复 6# 的帖子
然后上代码。这个代码还是在linux下用fortran77编译通过运行正常，仅供参考。求解用的是带矩阵的方法，比较古老的方法了，但对于这种小问题还游刃有余。
欢迎回帖讨论。
[ 本帖最后由 eldfxm 于 2008-5-15 11:19 编辑 ] 哦。。。。还是谢了。。 原帖由 xinyi826 于 2008-5-20 17:47 发表 http://www.simwe.com/forum/images/common/back.gif
哦。。。。还是谢了。。
谢谢支持:handshake 看看，谢谢Lz 看看我有没有贡献积分 看看能下不 linux下编译环境用起来顺手吗？ 原帖由 eldfxm 于 2008-5-10 04:15 发表 http://forum.simwe.com/images/common/back.gif
第一次在这里发新贴，请多支持。
该程序使用三节点单元(triangular element)mesh一个左端与墙固接的梁，上表面均布载荷为QP。
求解了弹性矩阵D，local stiffness matrix (SLM)，global stiffness matrix (GSM)等 ...
很有F77风格的代码，Q4和Q8平面膜元写起来并不复杂的，Q8膜元的精度已经很高了，不过建议用F90/95风格写，这样更容易看懂
F77还是有很多现在看来不合理的地方。
其实，纯粹搞研究而不是做项目的话，用matlab写是最简单快捷方便的 一般了，不大实用了，不过还是谢谢楼主 非常感谢楼主分享 仔细看了下楼主的手稿,个人觉得(24)公式有点问题,T应该是T,可能是笔误. 楼主,真乃牛人也 非常感谢！ 看看，谢谢Lz
页:
[1]
2