1、加权残值法程序设计大作业重力作用下的等截面圆柱两端固定 ,请采用加权残值配点法通过编程计算给出 :1.位移沿高度分布曲线2.应力沿高度分布曲线3.比较加权残值解与精确解的误差精确解：任务一：给出精确解位移 /应力沿高度分布曲线步骤 1： c+编程形成计算结果数据文件0 0.000.1 0.022 0.0x vFILE *fp;if(fp=fopen(”results.txt”,”w”)=NULL) printf(“cannot open file).n);exit(1);fprintf(fp,”%f/n”,v);fclose(fp);步骤 2：将数据文件导入 EXCEL画曲线图加权残值配点法：
2、取试函数：配点：代定系数方程：#include “stdafx.h“#define E 1.0e6#define RO 2500#define g 9.8FILE *fp;void main()int i;int n=100;double L=2.0;double x;double dx;double v;fp=fopen(“displace.txt“,“w“);dx=L/n;for(i=0;in+1;i+)x=dx*i;v=0.5*RO*g*x*(x-L)/E;fprintf(fp,“%0.2ft %en“,x,v);任务二：用高斯消元法程序计算代定系数任务三：给出加权残值法位移沿高度分布曲线任务四：给出加权残值法应力沿高度分布曲线任务五：将加权残值法和精确解位移沿高度分布曲线画到同一张图上，进行比较。
