满应力法设计超静定.docxVIP

河南科技大学学 院 土木工程学院班 级 ＿＿＿工程力学112＿＿学生姓名 刘冬日 期 2014.11.13 题目：应用满应力法求解三杆超静定问题三杆超静定结构桁架，其几何尺寸如图1-1所示。其中，要考虑材料密度的影响，材料常数为：图1-1解：题目中所定义A1为根据荷载对称，工况1和2反对称，工况3对称下的杆件1和杆件2的截面面积，A2为杆件3的截面面积。约束条件（st）：工况2条件下可以不必重复计算，和工况1结果相似。工况3条件下杆1、2的应力相同。程序设计思路:计算每种工况下各个杆件的应力比，对每个杆件的应力比取最大值作为该杆件的应力比r1,r2,r3.并且取r1、r2的最大值作为1号杆件和2号杆件的应力比，与r3一起进入到下一次迭代过程。精度设置为1e-6,比较小。观察下面两表可以得出结论在精度较小如1e-6的情况下，密度不同也会取得相同的截面面积；但是在精度为0.1时，迭代次数明显减少，杆件最后优化面积也不同。r1=max2(r1,r2);A1=A1*r1;A2=A2*r3;i++;W=100*(g1*2*sqrt(2)*A1+g2*A2);fprintf(fp,%d\t\t%f\t%f\t%f\t%f\t%f\t\n,i,A1,A2,r1,r2,W);}while (fabs(W-W1)Q);//fabs为取绝对值return 0;}//取三个数最大值double max3(double a,double b,double c){double max2(double a,double b);return max2(max2(a,b),c);}//取两个数最大值double max2(double x,double y){if (x=y)return x;else return y;}# include stdio.h# include math.hint main (){ FILE *fp; fp=fopen(answer.txt,w+);//结果输出到answer.txt文档中double max3(double x,double y,double z);double max2(double x,double y);int C1=2000,C2=1500,P=2000,i=0;double A1=1.0,A2=1.0,A,r11,r12,r13,r21,r22,r23,r31,r32,r33,r1,r2,r3,W,W1,g1,g2,Q=1e-6;//上面各变量含义为A杆件界面面积，r为应力比，W为重量，Q为精度printf(请输入密度g1,g2,用空格隔开\n);scanf(%lf%lf,g1,g2);printf(g1=%f,g2=%f\n,g1,g2);fprintf(fp,g1=%f,g2=%f\n,g1,g2);fprintf(fp,循环次数\tA1\t\tA2\t\t应力比1\t\t应力比2\t\t重量\t\n);do{A=2*A2*A1+sqrt(2)*A1*A1;W1=100*(g1*2*sqrt(2)*A1+g2*A2);r11=P*(A2+sqrt(2)*A1)/(A*C1);r12=P*sqrt(2)*A1/(A*C1);r13=P*A2/(A*C2);r21=sqrt(2)*P/(A*C1);r22=P*(A2+A1*sqrt(2))/(A*C1);r23=P*A2/(A*C2);r31=P*A1/(A*C1);r32=r31;r33=2*P*A1/(A*C1);r1=max3(r11,r21,r31);r2=max3(r12,r22,r32);r3=max3(r13,r23,r33);A.杆件1、2号的密度为0.1，杆件3的密度为为0.1表1-1迭代次数K1、2杆件截面面积A13杆截面面积A21、2杆应力比r1、r23杆应力比r3重量W0110.707107028427110.7071070.5857861.032726085786420.7302480.5395041.019545004957130.7445210.5109581.012375016781340.7537350.4925311.008125024414250.7598590.4802831.0054610294875............……………………360.7734590.4530821.0000001407543A.杆件1、2号的密度为0.1，杆件3的密度为为0.3表1-2迭代次数K1、2杆件截面面积A13杆截面面积A21、2杆应力比r1