连续体平面问题有限元分析.docVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
连续体平面问题的有限元分析 【题目】:正方形薄板四周受均匀载荷的作用,该结构在边界上受正向分布压力,P=1kN/m,同时在沿对角线y轴上受一对集中压力,载荷为2kN,板厚t=1,泊松比ν=0,见下图: Y 2kN 1kN/m x o 2kN 摘要:有限单元法作为一门课程在现实工程中的应用已经十分的广泛,利本文用计算机,结合FORTRAN语言和有限单元法的理论课知识,对平面应力问题的薄板进行内力值的计算和相应的结构分析,并可以将由计算机计算的实际结果与由用弹性力学理论知识计算出来的解答相比较。 关键词:连续体;平面应力;应力;节点;单元 理论依据与分析:此问题,为弹性力学里的平面应力问题,在板的内部,到处都有【1】:σz=0,τyz=0,τxz=0;σx=f1(x,y),σy=f2(x,y),τxy=f3(x,y),应力具有这种性质的问题,称为平面应力问题。弹性薄板在工程中应用很广泛,对于一些简单的情况,如等厚、单跨、无大孔口,外形规则(如矩形,圆形等)的薄板,已有一些解答和表格可资利用盘【2】。由于连续平板的连续性,仅需要取其在第一象限的四分之一部分研究计算,然后做出一些辅助线将平板分成若干部分,在为每个部分选择分子单元,采用此模型化为4个全等的三角形单元,利用其对称性,四分之一的边界约束,荷载可等效如图所示。 1kN/m 程序原理及实现: 用FORTRAN程序的实现。有节点信息文件NODE.IN和单元信息文件ELEMENT.IN,经过计算分析后输出一个一般性的文件DATA.OUT。模型基本信息由文件BASIC.IN生成。该程序的特点如下:问题类型:可用于计算弹性力学平面应力问题和平面应变问题。单元类型:采用常应变三角形单元。位移模式:用线性位移模式。载荷类型:节点载荷,非节点载荷应先转换为等效节点载荷。材料性质:弹性体由单一的均匀材料组成。约束方式:为“0”位移固定约束,为保证无刚体位移,弹性体至少应有三个自由度的独立约束。方程求解:针对半带宽刚度方程的GUASS消去法。输出文件:由手工生成节点信息文件NODE.IN,和单元信息文件ELEMENT,IN。结果文件:输出一般的结果文件DATA.OUT。 程序的原理如框图: 开始 输入数据(子程序READ_IN) BASIC.IN(基本信息文件) NODE.IN(节点信息文件) ELEMEENT.IN(单元信息文件) ↓ 形成单元刚度矩阵(子程序FORM_KE) ↓ 以半带存储方式形成整体刚度矩阵 (BAND K) ↓ 形成节点载荷向量(子程序FORM_P) ↓ 处理边界条件(子程序DO_BC) ↓ 求解方程获得节点位移(子程序SOLVE) ↓ 计算单元及节点应力(子程序) ↓ 输出文件 DATA.OUT ↓ 结束 说明:(1)主要变量: ID: 问题类型码:,ID=1时的平面应力问题,ID=2时的平面应变问题 N_NODE: 节点个数 N_LOAD: 节点荷载个数 N_DOF: 自由度,N_DOF=N_NODE*2(平面问题) N_ELE: 单元个数 N_BAND: 矩阵半带宽 N_BC: 有约束的节点个数 PE: 弹性模量 PR: 泊松比 PT: 厚度 LJK_ELE(I,3):单元节点编号数组,LJK_ELE(I,1),LJK_ELE(I,2),LJK_ELE(I,3)分别放单元I的三个节点的整体编号。X(N_NODE),Y(N_NODE):节点坐标数组,X(I),Y(I)分别存放节点I的

文档评论(0)

nnh91 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档