管的不规则压缩分析实例（原创）

shifang157300

PE塑料管（平面应变模型），受压后圆形截面（外径32mm，壁厚2.9mm）变成了椭圆形（长轴28mm，短轴18mm），计算管截面上的塑性区。

本例的关键是位移边界条件的施加，我的思路如下：
选择管截面最外面一圈节点（共Mnode个节点），获得其节点号，存入数组NUMP(Mnode)。使用循环语句分别对以上节点施加x方向和y方向位移。关键命令流如下：
!选择外面一周节点
!
lsel,s,,,1,4
NSLL,S,1
!
!
!获得已选择节点的节点号,存入NUMP()
*DIM,NUMP,ARRAY,1000
j=0
*do,i,1,Mnode
    *GET,KEY,NODE,i,NSEL
    *if,KEY,eq,1,then
        j=j+1
       NUMP(j)=i
    *endif
*enddo
/com, %J%
!在已选择的节点间循环
!
*do,i,1,j
    *GET,LOCX,NODE,NUMP(i),LOC,X   !获得节点x坐标
    *GET,LOCY,NODE,NUMP(i),LOC,Y   !获得节点y坐标
   *if,LOCX,LT,0,then
        D,NUMP(i),UX,-((D-SP)/2)*LOCX/(D/2)
   *endif
    *if,LOCX,GE,0,then
        D,NUMP(i),UX,-(D-SP)/2*LOCX/(D/2)
    *endif
!
   *if,LOCY,LT,0,then
        D,NUMP(i),UY,-((D-SZ)/2)*LOCY/(D/2)
   *endif
    *if,LOCY,GE,0,then
        D,NUMP(i),UY,-(D-SZ)/2*LOCY/(D/2)
    *endif
*enddo
等效塑性应变云图：

本例子的命令流：
欢迎各位提出其他方法。

watterhit

取楼主的1/2模型，做接触分析....