有限元大作业matlab---课程设计例子.doc

有限元大作业matlab---课程设计例子

有限元大作业程序设计

学校：天津大学

院系：建筑工程与力学学院

专业：01级工程力学

姓名：刘秀

学号：\\\\\\\\\\\

指导老师：

连续体平面问题的有限元程序分析

[题目]：

如图所示的正方形薄板四周受均匀载荷的作用，该结构在边界上受正向分布压力，，同时在沿对角线y轴上受一对集中压力，载荷为2KN，若取板厚，泊松比。

2kN

2kN

2kN

1kN/m

[分析过程]：

由于连续平板的对称性，只需要取其在第一象限的四分之一部分参加分析，然后人为作出一些辅助线将平板“分割”成若干部分，再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性，四分之一部分的边界约束，载荷可等效如图所示。

源程序文件：

chengxu.for

程序需读入的数据文件：

BASIC.IN，NODE.IN，ELEMENT.IN（需要手工生成）

程序输出的数据文件：

DATA.OUT

（4）数据文件格式：

需读入的模型基本信息文件BASIC.IN的格式如下表

栏目

格式说明

实际需输入的数据

基本模型数据

第1行，每两个数之间用“，”号隔开

问题类型，单元个数，节点个数，有约束的节点数，有载何的节点数

材料性质

第2行，每两个数之间用“，”号隔开

弹性模量，泊松比，单元厚度

节点约束信息

在材料性质输入行之后另起行，每两个数之间用“，”号隔开

LJK_U(N_BC,3)

位移约束的节点编号，该节点x方向约束代码，该节点y方向代码，

节点荷载信息

在节点约束信息输入行之后另起行，每两个数之间用“，”号隔开

P_IJK(N_LOAD,3)

载荷作用的节点编号，该节点x主向载荷，该节点y方向载荷，……

需读入的节点信息文件NODE.IN的格式如下表

栏目

格式说明

实际需输入的数据

节点信息

每行为一个节点的信息（每行三个数，每两个数之间用空格或“，”分开）

ND_ANSYS(N_NIDE)

节点号，该节点的x坐标，该节点y方向坐标

需读入的单元信息文件ELEMENT.IN的格式如下表

栏目

格式说明

实际需输入的数据

单元信息

每行为一个单元的信息（每行有14个整型数，前4个为单元节点编号，对于3节点编号，第4个节点编号与第3个节点编号相同，后10个数无用，可输入“0”，每两个整型数之间用至少一个空格分开）

NE_ANSYS(N_ELE,14)

单元的节点号1（空格）单元的节点号2（空格）单元的节点号3（空格）单元的节点号4（空格）

0（空格）0（空格）0（空格）0（空格）0（空格）0（空格）0（空格）0（空格）0（空格）0

输出结果文件DATA.OUT格式如下表

栏目

实际输出的数据

节点位移

IRESULT_N(2*I_1)RESULT_N(2*I)

节点号x方向位移y方向位移

单元应力的三个分量

IESTE_ELE(IE,1)STE_ELE(IE,2)STE_ELE(IE,3)

单元号x方向应力y方向应力剪切应力

节点应力的三个分量

ISTS-ND(I,1)STS-ND(I,2)STS-ND(I,3)

节点号x方向应力y方向应力剪切应力

[算例原始数据和程序分析]：

（1）模型基本信息文件BASIC.IN的数据为

1,4,6,5,3

1.,0.,1.

1,1,0,2,1,0,4,1,1,5,0,1,6,0,1

1,-0.5,-1.5,3.,-1.,-1,6,-0.5,-0.5

（2）手工准备的节点信息文件NODE.IN的数据为

10.02.0

20.01.0

31.01.0

40.0.

51.00.

62.00.

（3）手工准备的单元信息文件ELEMENT.IN的数据为

12330000111101

24550000111102

53220000111103

3