[分享]例子－－解超静定梁，多质点问题
<P>！解超静定梁，多质点<br><br>/PREP7!进入前处理模块PREP7</P><P>N, 1, 0, 0!定义各个结点</P>
<P>N, 2, 3, 0</P>
<P>N, 3, 6, 0</P>
<P>N, 4, 9, 0</P>
<P>N, 5, 12, 0</P>
<P>N, 6, 15, 0 </P>
<P>N, 7, 18, 0</P>
<P>N, 8, 21, 0 </P>
<P>N, 9, 24, 0</P>
<P>N, 10, 27, 0 </P>
<P>N, 11, 30, 0 </P>
<P>NPLOT ! 图形中不显示结点号码</P>
<P>NPLOT, 1 ! 图形中显示结点号码</P>
<P>NLIST ! 在直角坐标系统下,列出结点资料</P>
<P>DSYS, 1 ! 列出资料时,转换至圆柱坐标</P>
<P>NLIST ! 在圆柱坐标系统下,列出结点资料</P>
<P>ET, 1, BEAM3 !定义第一类单元为平面梁单元BEAM3</P>
<P>MP, EX, 1, 207e5 !定义第一类材料的弹性模量EX</P>
<P>R, 1, 1, 0.020833, 0.5!平面梁的实常数:截面积,抗弯惯性矩和截面高度</P>
<P>E, 1, 2 !定义各个单元两端结点的连接关系</P>
<P>E, 2, 3</P>
<P>E, 3, 4 </P>
<P>E, 4, 5</P>
<P>E, 5, 6 </P>
<P>E, 6, 7</P>
<P>E, 7, 8 </P>
<P>E, 8, 9</P>
<P>E, 9,10 </P>
<P>E, 10, 11</P>
<P>EPLOT !在图形窗口显示单元</P>
<P>ELIST !列出单元信息</P>
<P>FINISH!退出前处理模块</P>
<P>/SOLU !进入求解模块SOLUTION</P>
<P>ANTYPE, STATIC !申明求解类型是静力分析</P>
<P>OUTPR, BASIC, ALL ! 在输出结果中，列出元素的结果</P>
<P>D, 1, UX, 0, , , , , UY, ROTZ !定义1号结点约束所有自由度UX,UY和ROTZ </P>
<P>D, 11, UX, 0, , , , ,UY !约束11号结点的UX和UY自由度</P>
<P>SFBEAM, 1, 1, PRES, 0, 60 !定义各个单元上的线性分布荷载</P>
<P>SFBEAM, 2, 1, PRES, 60, 120</P>
<P>SFBEAM, 3, 1, PRES, 120, 180</P>
<P>SFBEAM, 4, 1, PRES, 180, 240</P>
<P>SFBEAM, 5, 1, PRES, 240, 300</P>
<P>SFBEAM, 6, 1, PRES, 300, 360</P>
<P>SFBEAM, 7, 1, PRES, 360, 420</P>
<P>SFBEAM, 8, 1, PRES, 420, 480</P>
<P>SFBEAM, 9, 1, PRES, 480, 540</P>
<P>SFBEAM, 10, 1, PRES, 540, 600</P>
<P>F, 6, FY, -1000 !定义在6号结点上的沿-Y方向大小为1000的集中力</P>
<P>SOLVE !开始求解</P>
<P>FINISH!退出求解模块</P>
<P>/POST1!进入后处理模块POST1</P>
<P>PRDISP ! 列出结点位移数据</P>
<P>ETABLE,IMOMENT, SMISC, 6 ! 建立单元结果表，元素I点力矩</P>
<P>ETABLE,JMOMENT, SMISC, 12 ! 建立单元结果表，元素J点力矩</P>
<P>ETABLE, ISHEAR, SMISC, 2 ! 建立单元结果表，元素I点剪力</P>
<P>ETABLE, JSHEAR, SMISC, 8 ! 建立单元结果表，元素J点剪力</P>
<P>PRETAB ! 列出单元表计算结果</P>
<P>/TITLE, SHEAR FORCE DIAGRAM !定义图形窗口标题</P>
<P>PLLS, ISHEAR, JSHEAR ! 绘制结构剪力分布图</P>
<P>/TITLE, BENDING MOMENT DIAGRAM!定义图形窗口标题</P>
<P>PLLS, IMOMENT, JMOMENT ! 结构弯矩分布图</P>
<P>FINISH !退出后处理模块POST1</P>
[此贴子已经被作者于2006-4-25 13:57:05编辑过]
页:
[1]