耦合通用小模块

kawich

下面是一个节点耦合的小模块
具体用法：将要进行耦合的两类节点分别放在cm_1集合和cm_2集合中，cm_1为主体节点集合，cm_2为载体节点集合。然后就可以直接调用这个模块了。一般来说cm_2节点密度大于等于cm_1节点密度时耦合效果较好。


----------------------------------------------------------------------
cmsel,s,CM_1
*get,count1_nd,node,0,count  !201
!--------存放集合1的节点号 cm1_nd
*del,cm1_ndc
*dim,cm1_ndc,array,count1_nd ! 6
*do,i,1,count1_nd
    *if,i,EQ,1,then
        *get,nnum,node,0,num,min
    *else
        *get,nnum,node,nnum,nxth
    *endif
    cm1_ndc(i)=nnum
*enddo
!--------存放所有集合2的节点号 cm2_nd
allsel,all
cmsel,s,CM_2

*get,count2_nd,node,0,count
*del,cm2_ndc
*dim,cm2_ndc,array,count2_nd
*do,i,1,count2_nd
    *if,i,EQ,1,then
        *get,nnum,node,0,num,min
    *else
        *get,nnum,node,nnum,nxth
    *endif
    cm2_ndc(i)=nnum
*enddo
!*********************去重
allsel,all
cmsel,s,CM_2
cm,fdzctnc,node
*del,cm2cp_ndc
*dim,cm2cp_ndc,array,count1_nd
*do,i,1,count1_nd
    cm2cp_ndc(i)=0
*enddo

     allsel,all
     cmsel,s,fdzctnc  

*do,i,1,count1_nd
    k=nnear(cm1_ndc(i))
    cm2cp_ndc(i)=k
    flag=0
                                !--------防止集合2的节点被重复耦合，
                                !---即要求每个集合1的节点与不同的集合2节点耦合在一起  
  *do,j,1,count1_nd
       *if,i,EQ,j,then         
           *cycle
       *endif
        
       *if,cm2cp_ndc(j),EQ,0,then
         *exit
       *endif

        *if,cm2cp_ndc(i),EQ,cm2cp_ndc(j),then   !在除去重合的节点以外的集合重

新选取
             nsel,u,node,,cm2cp_ndc(j)
             flag=1  
             *exit         
        *endif
    *enddo   
!*******************如果原节点集合发生改变则重选，反之则否  
    *if,flag,eq,1,then   
    k=nnear(cm1_ndc(i))
    cm2cp_ndc(i)=k
    *endif
        
    nsel,u,node,,k           
!*******************        
*enddo
!!----------开始耦合
allsel,all
/prep7
*do,i,1,count1_nd
    cp,next,UX,cm1_ndc(i),cm2cp_ndc(i)
    cp,next,UY,cm1_ndc(i),cm2cp_ndc(i)
    cp,next,UZ,cm1_ndc(i),cm2cp_ndc(i)
*enddo

kawich

自己顶一下 ~~  

pfsr

感谢分享，如果是原创可以加技术分。

kawich

也算不上原创，主体思想框架来源于一个工程项目的命令流文件，我只是做了些修改增强了它的可移植性，使它能够通用

woodylacurs

好东西，先收藏了～

sky25

不错，收藏了。。。谢谢诶

wangzxiaoh

好 我也收藏了

尔冭美

感谢分享