用伪代码方式写出有限元求解步骤（Write the finite element solution steps in pseudo code）.docVIP

用伪代码方式写出有限元求解步骤（Write the finite element solution steps in pseudo code） ==================================================== * * 有限元程序* * * * ==================================================== 主程序 昏暗的NJ，东北、NZ、苗、NJ2；*定义变量 昏暗的JZ（NJ，2）、ZC（新西兰，3）、普（非常任法官，3），JM（NE，3）； * ZC（I，1）、ZC（I，2）、ZC（I，3）分别表示第我单元的编号，水平约束，竖向约束；* JM（I，1），JM（I，2），JM（I，3）分别表示地我单元的编号，水平载荷，竖向载荷； 读NJ、NE、NZ、苗； 对于E = 1，NE； 电话输入（E、JZ、ZC、NJ、JM，EO，穆山，娄，jxm、NP、NL）； 电话dygd（E，问）； 下一个E 电话ztgd（NJ2、NE、JM）； 叫救命（NJ2，非常任法官）； 叫ISC（新西兰，NJ2，KZ，P）； 打电话给高斯（KZ，P）； 呼叫压力（P，pp，s，NE，e，JM）； 结束 ---------------------------------------------------------------------------- 子程序输入（E、JZ、ZC、NJ、JM，EO，穆山，娄，jxm、NP、NL）*输入模块 定义x（NJ），Y（NJ），JM（NE，3）、ZC（新西兰，3）、普（非常任法官，3），PL（非常任法官，5）*定义数组 *输入数据 读E0，穆山，T，柔*输入弹性模量，泊松比，厚度，容重 读X（我），Y（我），我= 1。新泽西*输入节点坐标 读JM（I，J），I = 1。NE，j = 1，3×输入节点号码数组 ZC（i，j）读，我= 1，NR，j = 1，3×输入节点约束 如果NP＞0×输入节点荷载 读PJ（I，J），i = 1，NP，J = 1 3。 最后如果 如果NL 0×输入线荷载 读取PL（i，j），i = 1…NL，j = 1….. 5 最后如果 结束 ----------------------------------------------------------------------------- 子程序dygd（E，问单元刚度模块）* 二= JM（E，1），JJ = JM（E，2），MM = JM（E，3）；*节点坐标 双= JZ（JJ，2）- JZ（毫米，2）*求BI、CI（I，J，M） 北京= JZ（毫米，2）- JZ（II，2） BM = JZ（II，2）- JZ（JJ，2） CI = - JZ（JJ，1）+ JZ（毫米，1） CJ = - JZ（毫米1）+ JZ（II，1） CM = - JZ（II，1）+ JZ（JJ，1） = 0.5 *（双* * *求面积bj-bj CI） 如果问 1然后 我= 3；J = 1.6 B（i，j）＝0，s（i，j）＝0； 下一个我； 最后如果 B（1,1）=双，B（1,3）= BJ，B（1,5）= BM； B（2,2）= CI，B（2,4）= CJ，B（2,5）=厘米； B（3,1）= CI，B（3,3）= CJ，B（3,5）=厘米； B（3,2）=双，B（3,4）= BJ，B（3,6）= BM； 如果jxm = 1 *判断问题类型 E0 = E0 /（1×亩），亩=亩/（1）； 最后如果； 我= 3；J = 1 d（i，j）＝0； 下一个 D（1,1）= E0 /（1×亩）； D（2,2）= D（1,1），D（2,1）= D（1,1）*亩 D（1,2）= D（2,1），D（3,3）= D（1,1）*（1）/ 2； 我= 3；J = 1.6 K = 1 S（i）， (j) = s (i, j) + d (i, k) * b (k, j); next k next i, j if 2, then 求单元刚度矩阵 * box for i = 1, j = 1.6 ke (i, j) = 0; next i, j; for i = 1, j = 1.6 for k = 1.3; ke (i, j) = m (i, j) + b (k, i) * s (k, j); next k next i, j; end if end ---------------------------------------------------------------------- 子程序 ztgd (nj2, ne, jm) * 总刚度模块 for i = 1, nj2; j = 1, nj2