有限元程序SOLID2D分析报告Word下载.docx
《有限元程序SOLID2D分析报告Word下载.docx》由会员分享,可在线阅读,更多相关《有限元程序SOLID2D分析报告Word下载.docx(15页珍藏版)》请在冰点文库上搜索。
B、课程设计思路的描述;
C、程序的理论基础;
D、程序设计思路、框图、数据结构、数据流;
E讨论(含存在的问题、希望进一步开展的工作,需要完善之处,心得体会等)。
时间安排:
第1周。
指导教师签名:
2012年6月12日
系主任(或责任教师)签名:
有限元程序SOLID2D分析
1课程设计分析
1.1课程设计的目的
课程设计是在学过相关知识的情况下为了对该知识有更好的理解和掌握所进行的一种设计。
它的目的在于使学生灵活运用所学的知识解决特定的实际问题,学会综合运用各种知识来理解和掌握该科目的实际运用方法。
此次《计算力学》课程设计的目的是:
●培养学生的自我独立思考能力,进一步加深学生对《计算力学》知识的理解;
●提高学生综合运用知识的能力,学会用多方面的知识解决问题;
●学会对实际问题有较正确的理解,然后从更好的方向解决问题;
●对学生程序的读写能力的培养;
●培养学生查阅资料的好习惯,使课程设计更加完善。
1.2SOLID2D课程设计的分析
1.2.1SOLID2D单元介绍
SOLID2D是一个从三维单元退化而成的二维实体单元(类似板、壳单元),可以解决线性平面问题及轴对称问题。
它是具有四个或八个结点的分析线性问题的单元,在热、电、力研究方面都具有良好的运用能力。
1.2.2程序目的
该程序分为
(1),
(2)两个部分,在这里仅对第
(2)部分的目的进行描述。
计算刚度比例瑞利阻尼余量;
计算平面与轴对称单元余量;
计算增强位移模式;
计算平面与轴对称弹性单元的应力应变,瑞利阻尼等;
计算平面/轴对称FEAPpv单元增强应力。
1.2.3程序中所运用的相关知识
概念性知识:
应力、应变、瑞利阻尼、加速度、矩阵、惯性、位移模式、刚度、高斯积分、形函数、应力发散、权、雅克比行列式、切线模量、收敛性、耦合、余量等等。
计算性知识:
矩阵的行循环,列循环、矩阵的转置、应力,应变的计算、矩阵的乘积、计算高斯积分点和权、局部迭代、验证收敛性、消除增强模式等等。
2程序流程图
2.1刚度比例瑞利阻尼余量的流程图
结束
由于阻尼原因修正应力
3个应变量(厚度/圈)
赋初值
定义数据
类型
调用函数
开始
2.1.1由于阻尼原因修正应力流程图
1=>
i
j
sig(i)=sig(i)+dd(i,j)*eps(j)
1+i=>
i>
4?
N
Y
1+j=>
j>
N
end
2.2平面与轴对称余量程序流程图
定义数据类型
计算应力矢量发散性(P)
为元素存储时间历史点数据
计算加速度
瑞利质量阻尼:
计算速度
行循环
计算重力,热,惯性和应力贡献
2.2.1计算应力矢量发散性(P)流程图
int(d(74))>
0?
YN
b2=d(12)*dm
int(d(75))>
b1=d(11)*dm
b1=d(11)+prldv(int(d(74)))*d(71)
YN
b2=d(12)+prldv(int(d(75)))*d(72)
rr=d(4)
2.2.2瑞利质量阻尼:
计算速度流程图
d(77)=0.0d0?
vc
(1)=0.0d0
vc
(2)=0.0d0
vc
(1)=vc
(1)+shp(3,j)*vl(1,j)
vc
(2)=vc
(2)+shp(3,j)*vl(2,j)
J>
nel?
N
vc
(1)=rr*vc
(1)*cfac*d(77)
vc
(2)=rr*vc
(2)*cfac*d(77)
Y
aj0=lfac*d(77)*rr
p(1,j)=p(1,j)-(vc
(1)+aj0*vl(1,j))*shp(3,j)*xsj
p(2,j)=p(2,j)-(vc
(2)+aj0*vl(2,j))*shp(3,j)*xsj
2.3增强位移模式流程图
计算增强应变'
不符'
形函数
2.4平面与直线轴弹性单元程序流程图
计算B的转置乘以D乘以j乘以w
2.4.1计算高斯积分点和权流程图
nel=6ornel=7?
l=d(5)
quad=.true.
l=7
quad=.false.
calltint2d(l,lint,el)
下一步
callint2d(l,lint,sg)
l*l=lint?
2.4.2将对称半径乘以雅克比行列式
xsj0=0.0d0
Stype=3?
xsj0=xsj(l)
xsj(l)=xsj(l)*xx
shpr(j)=shp(3,j,1)/xx
1+j=j
2.4.3切模量中形成质量矩阵流程图
k
mass(k,j)=mass(j,k)
1+k=>
k>
N
j1=0
j1=j1+1
k1=0
k1=k1+1
s(j,k)=s(j,k)+mass(j1,k1)
s(j+1,k+1)=s(j+1,k+1)+mass(j1,k1)
ndf+k=>
ndf*nel?
ndf+j=>
ndf*nel?
2.5计算平面/轴对称FEAPpv单元增强应力流程图
输入数据
收回增强模式(节省上次迭代)
积分计算和权
仅初始化历史变量
计算形函数
轴对称平面
用局部迭代计算强化模式
初始化历史应力应变
计算增强应变形状函数
计算点的应变
贮存时间积历史数据元素
规模模量和应力
计算应力和切线模
增强余量计算
刚度计算
计算d*b矩阵=a
计算增量位移增强增强模式
验证收敛性
时间积分,以确定静态或动态
瑞利阻尼效应
集中质量的影响
增强耦合阵列
消除增强模式
构建静态凝结
计算减少残余
计算并输出单元变量
计算主应力值
积应力值
输出格式
2.5.1计算增强余量流程图
bb(2*j-1)=bb(2*j-1)-sigl
(1)*shpi(1,j)-sigl(4)*shpi(2,j)
bb(2*j)=bb(2*j)-sigl
(2)*shpi(2,j)-sigl(4)*shpi(1,j)
3课程设计心得
通过这次《计算力学》课程设计我所做的题目是SOLID2DNFEAPPY20程序分析。
该程序是SOLID2D在许多方面上的应用。
例如:
热、电、受力等等。
并且给出了在各个情况下该实体单元内部的力学性能解答方法与过程。
SOLID2D是一个二维的实体单元,它是由三维实体单元退化而来。
但是它也同样具有一般实体单元所具有的性能。
它是一个具有4个结点或8个结点的单元,解决平面与轴对称单元的线性问题。
但该单元也有一定的局限性。
它不解决非线性问题,也不再具有三维实体所特有的性质等等。
在分析过程中,我也遇到了许多困难。
首先是程序的理解。
程序是在FORTRAN环境下编写的,对我们的FORTRAN运用能力有较高的要求,而且许多注解是英文,又对我们的英语特别是专业英语要求较高。
在这个方面我的能力是比较弱的,但是通过老师对我们的帮助,再加上对相关资料的查阅,我也大致的了解了全部的程序过程,并对其目的、步骤有了较好的理解。
其次是程序中所运用的相关数学知识和力学知识。
它要求较高的数学能力,例如高斯积分,权,雅克比行列式,迭代,收敛性的验算等等。
还有瑞利阻尼,切模量,耦合等等的物理概念也比较不常见。
在上个学期学习的《计算力学》课程中,我们学到的相关知识都在课程设计中得到很好的运用,让我们不仅仅是在理性上了解这些知识,更是从课程设计中直观的了解到了这些知识的用途,可以解决哪些实际的问题。
让我们对该知识的掌握提供了更好的基础条件。
虽然这次的课程设计由于时间紧迫,我们没有对相关的知识和问题进行深入的探讨与研究,但却已经给我们带来了很多收获。
我的许多能力得到了提高,特别是分析问题和综合运用的能力。
我想今后每次有这样学习的机会都要好好珍惜,尽自己最大的努力做到最好,我们自己才能提高,获得更多的知识来充实自己。
4课程设计参考资料
王勖成编著,《有限单元法》清华大学出版社2001年;
T.R.Chandrupatla,A.D.Belegundu,《IntroductiontoFiniteElementsinEngineering》机械工业出版社2005年1月;
O.C.Zienkiewicz,R.L.Taylor,《TheFiniteElementMethod》ButterworthHeinemann,2000
潘金贵等,《现代计算机常用数据结构和算法》南京大学出版社,1994;
Diomidis.Spinellis著,赵学良译,《代码阅读方法与实践》,北京,清华大学出版社,2004年3月。