ansys 中画弯矩图（ansys 中画弯矩图）.docVIP

ansys 中画弯矩图（ansys 中画弯矩图） then itfinum = 6 itfjnum = 12 * * * * * * you * endif ! 对不需平均的节点进行处理 if, unassembly, none, then ! 选出不进行处理的节点包并获取不进行处理节点的数目 cmsel, unassembly * get, unnodnum, node 0, count ! 定义长度为unnodnum的数组 (unnod), 以存放选中单元的单元编号 * unnod, array, unnodnum ! 将选中单元的编号按顺序存入数组unnod do i, 0, unnodnum - 1.1 unnod (i + 1) = ndnext (i) * enddo * * * * * * you unnodnum = 0 * endif ! 选出所需的单元和节点包 cmsel, eassembly cmsel, nassembly ! 获得当前选中单元总数 (存入变量selelenum) * get, selelenum, elem, 0, count ! 定义长度为selelenum的数组 (elenum), 以存放选中单元的单元编号 * elenum, array, selelenum ! 将选中单元的编号按顺序存入数组elenum do i, 0, selelenum - 1.1 elenum (i + 1) = elnext (i) * enddo ! 获得当前选中节点总数 (存入变量selnodnum) * get, selnodnum, node 0, count ! 定义长度为selnodnum的数组 (nodnum), 以存放选中单元的单元编号 * nodnum, array, selnodnum ! 将选中单元的编号按顺序存入数组nodnum do i, 0, selnodnum - 1.1 nodnum (i + 1) = ndnext (i) * enddo ! 定义所需的线性单元内力etable, 节点i的内力存入数组itnfi. ! 节点j的内力存入数组itnfj itnfi, smisc, itfinum barn. itnfj, smisc, itfjnum barn. ! 定义所需的结果数组, 并将其置零 itnfineo barn, smisc, 5. itnfineo, sadd, itnfi, 1 itnfjneo barn, smisc, 11. itnfjneo, sadd, itnfj, 1 * * * * * * ce, k, 1, selnodnum, 1 ! 处理不需平均的节点 index = 0. if, unnodnum, ge, 1,. do i, 1, unnodnum if, nodnum eq (k), (j), then, unnod index = 1. * * * * * * you * endif * enddo * * * * * * you * endif do i selelenum, 1, 1, ! 选出和节点k相连的线性单元中, i节点 为节点k的单元编号 (对线性单元而言) if, nelem (elenum (j), (k), eq, nodnum). (j) elenum elei. * i * * * * * * you * endif * enddo do i selelenum, 1, 1, ! 选出和节点k相连的线性单元中, j节点 为节点k的单元编号 (对线性单元而言) if, nelem (elenum (j), (k), eq, which nodnum elej = elenum (j). * i * * * * * * you * endif * enddo if index, eq, 0,. if, elej,, 0,, then. 有可能出现elej为0的情况 ! 取出i节点为节点k的单元的i节点端的内力放入参数etelei * get. Etelei, elem, election, etab, itnfi ! 取出j节点为节点k的单元的j节点端的内力放入参数etelej * get, etelej, elem, elej, etab, itnfj ! 平均节点k的单元的i节点端的内力和节点k的单元的j节点端的内力 Etave = (etelei + etelej) / 2 ! 将平均后的内力存入结果数组中 Detab, election, itnfineo, etave Detab, elej, itnfjneo, etave * else * endif * else * endif * enddo