#include "stdafx.h"
#include<stdio.h>
#include<math.h>
#include <conio.h>
#define NE 5
#define E 210000000000
#define u 0.3
#define Fp -1000
#define h 0.1
#define A 1
int l[4]={0},m[4]={0},n[4]={0};
float Fd[4]={0};
float uv[4]={0};
float q3[3]={0};
static float uva[10]={0};
int jm[NE+1][3]={{0,0,0},{0,0,0},{0,0,1},{0,2,1},{0,0,2},{0,2,2}};//坐标
static float ke[4][4]={{0}},/*基本方程,最后一列为约束条件*/krs[3][3]={{0}},k[10][10]={{0}},x[4]={0};/*总纲*/
int i,j,p,q,a1,b1,c1,z,a,b,c;
void jdugd(int a,int b,int c);//刚度矩阵
void bjtj(int a1,int b1);//边界条件
void fun(float ax,float ay,float bx,float by);//计算位移
void SDB(int a,int b,int c);//计算应力
void main()
{
jdugd(1,3,2);//一单元
jdugd(2,3,5);//二单元
jdugd(2,5,4);//三单元
printf("=====================总刚度矩阵=====================\n");
for(i=0;i<=9;i++)
{if (j=9)printf("\n");
for(j=0;j<=9;j++)
printf("%g ",k[i][j]);//总刚度矩阵
}
printf("\n");
bjtj(3,5);//引入边界条件//三点总刚度矩阵生成
fun(0,0,0,-1000);//给出位移
printf("=====================基本方程=====================\n");
for( i=0;i<=3;i++)//基本方程
{if (j=9)printf("\n");
for( j=0;j<=3;j++)
printf("%g ",ke[i][j]);}
printf("\n ");
printf("=====================位移=====================\n");
printf("3X= 3Y= 5X= 5Y= \n");
for (i=0;i<4;i++)
{ printf(" %g",uv[i]); }
printf("\n ");
printf("=====================单元一应力=====================\n");
printf("SX= SY= SXY=\n");
SDB(1,3,2);//单元一
for (i=0;i<3;i++)
{ printf(" %g",q3[i]); }
printf("\n ");
printf("=====================单元二应力=====================\n");
- 1
- 2
前往页