有限元大作业matlab课程设计例子范文
相关主题
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
有限元大作业matlab课程设计例
子
有限元大作业程序设计
学校:天津大学
院系:建筑工程与力学学院
专业:01级工程力学
姓名:刘秀
学号:\\\\\\\\\\\
指导老师:
连续体平面问题的有限元程序分析
[题目]:
如图所示的正方形薄板四周受均匀载荷的作用,该结构在边界上受正向分布压力,m kN p 1=,同时在沿对角线y 轴上受一对集中压力,载荷为2KN ,若取板厚1=t ,泊松比0=v 。
[分析过程]:
由于连续平板的对称性,只需要取其在第一象限的四分之一部分参加分析,然后人为作出一些辅助线将平板“分割”成若干部分,再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性,四分之一部分的边界约束,载荷可等效如图所示。
[程序原理及实现]:
用FORTRAN程序的实现。由节点信息文件NODE.IN和单元信息文件ELEMENT.IN,经过计算分析后输出一个一般性的文件DATA.OUT。模型基本信息由文件为BASIC.IN生成。
该程序的特点如下:
问题类型:可用于计算弹性力学平面问题和平面应变问题
单元类型:采用常应变三角形单元
位移模式:用用线性位移模式
载荷类型:节点载荷,非节点载荷应先换算为等效节点载荷
材料性质:弹性体由单一的均匀材料组成
约束方式:为“0”位移固定约束,为保证无刚体位移,弹性体至少应有对三个自由度的独立约束
方程求解:针对半带宽刚度方程的Gauss消元法
输入文件:由手工生成节点信息文件NODE.IN,和单元信息文件ELEMENT.IN
结果文件:输出一般的结果文件DATA.OUT
程序的原理如框图: