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

文档仅供参考

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

栏目 节点位移 单元应力的三个分量 节点应力的三个分量 实际输出的数据 I RESULT_N(2*I_ 1) RESULT_N(2*I) 节点号 x方向位移 y方向位移 IE STE_ELE(IE,1) STE_ELE(IE,2) STE_ELE(IE,3) 单元号 x方向应力 y方向应力 剪切应力 I STS-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的数据为

1 0.0 2.0 2 0.0 1.0 3 1.0 1.0 4 0. 0. 5 1.0 0. 6 2.0 0.

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

9

2020年4月19日

文档仅供参考

3 5 6 6 0 0 0 0 1 1 1 1 0 4 （4）源程序文件chengxu.for为：

PROGRAM FEM2D DIMENSION

IJK_ELE(500,3),X(500),Y(500),IJK_U(50,3),P_IJK(50,3), &RESULT_N(500),AK(500,100) 20

DIMENSION STS_ELE(500,3),STS_ND(500,3) OPEN(4,FILE='BASIC.IN') OPEN(5,FILE='NODE.IN') OPEN(6,FILE='ELEMENT.IN') OPEN(8,FILE='DATA.OUT') OPEN(9,FILE='FOR_POST.DAT')

READ(4,*)ID,N_ELE,N_NODE,N_BC,N_LOAD IF(ID.EQ.1)WRITE(8,20) IF(ID.EQ.2)WRITE(8,25)

FORMAT(/5X,'=========PLANE

STRESS

PROBLEM========') 25

FORMAT(/5X,'=========PLANE

STRAIN

PROBLEM========')

CALL

READ_IN(ID,N_ELE,N_NODE,N_BC,N_BAND,N_LOAD,PE,PR,PT, & IJK_ELE,X,Y,IJK_U,P_IJK)

10

2020年4月19日

文档仅供参考

CALL BAND_K(N_DOF,N_BAND,N_ELE,IE,N_NODE,

& IJK_ELE,X,Y,PE,PR,PT,AK)

CALL

FORM_P(N_ELE,N_NODE,N_LOAD,N_DOF,IJK_ELE,X,Y,P_IJK, & RESULT_N)

CALL DO_BC(N_BC,N_BAND,N_DOF,IJK_U,AK,RESULT_N) CALL SOLVE(N_NODE,N_DOF,N_BAND,AK,RESULT_N) CALL

CAL_STS(N_ELE,N_NODE,N_DOF,PE,PR,IJK_ELE,X,Y,RESULT_N, & STS_ELE,STS_ND) c to putout a data file

WRITE(9,70)REAL(N_NODE),REAL(N_ELE)

70 FORMAT(2f9.4)

WRITE(9,71)(X(I),Y(I),RESULT_N(2*I-1),RESULT_N(2*I), & STS_ND(I,1),STS_ND(I,2),STS_ND(I,3),I=1,N_NODE) 71 FORMAT(7F9.4)

WRITE(9,72)(REAL(IJK_ELE(I,1)),REAL(IJK_ELE(I,2)), &REAL(IJK_ELE(I,3)),REAL(IJK_ELE(I,3)),

&STS_ELE(I,1),STS_ELE(I,2),STS_ELE(I,3),I=1, N_ELE) 72 c

CLOSE(4)

11

2020年4月19日

FORMAT(7f9.4)

文档仅供参考

CLOSE(5) CLOSE(6) CLOSE(8)

CLOSE(9) END c

c to get the original data in order to model the problem

SUBROUTINE

READ_IN(ID,N_ELE,N_NODE,N_BC,N_BAND,N_LOAD,PE,PR, &PT,IJK_ELE,X,Y,IJK_U,P_IJK) DIMENSION

IJK_ELE(500,3),X(N_NODE),Y(N_NODE),IJK_U(N_BC,3), & P_IJK(N_LOAD,3),NE_ANSYS(N_ELE,14) REAL ND_ANSYS(N_NODE,3) READ(4,*)PE,PR,PT

READ(4,*)((IJK_U(I,J),J=1,3),I=1,N_BC) READ(4,*)((P_IJK(I,J),J=1,3),I=1,N_LOAD) READ(5,*)((ND_ANSYS(I,J),J=1,3),I=1,N_NODE) READ(6,*)((NE_ANSYS(I,J),J=1,14),I=1,N_ELE)

DO 10 I=1,N_NODE

2020年4月19日

12