Hypermesh二次开发系列003

最新推荐文章于 2024-08-10 20:30:00 发布

这个就很牛啊

最新推荐文章于 2024-08-10 20:30:00 发布

阅读量3.9k

点赞数 4

分类专栏： HyperMesh二次开发

本文链接：https://blog.csdn.net/gb11235/article/details/115498330

版权

HyperMesh二次开发专栏收录该内容

4 篇文章 24 订阅

订阅专栏

1 *createmark 应用案例

1.1 需求描述

一个面被划分过网格了，需要在某一个指定位置，扩个孔出来，但现在面的几何信息已经被删除了。如图1所示
在这里插入图片描述

图1 功能描述

1.2 问题分析

这类问题暂且线定义为网格编辑。没有几何信息，只能直接操作和编辑单元与节点。使用界面操作的思路：

在指定位置创建临时节点：Geom>>nodes
据此为圆心做圆（创建几何元素：线）：Geom>>lines
删除完全在圆内的单元，可能会留下锯齿状单元。
使用project功能把上一步生成的锯齿状单元的节点投影到圆周上：Tool>>Project>>to line

上述是操作步骤，过程效果如图2。
在这里插入图片描述

图2 过程效果示意图是的，此过程GUI操作过程繁琐，而且容易出错。最好的办法就是不要遇到这样的需求，如果不可以的话，退而求次按照上述步骤细心操作，实在没办法就学习下hypermesh的二次开发吧，至少能做出如图3的效果，只要源码跑起来，指哪里就可以扩哪里，想扩多大就扩多大。

在这里插入图片描述

图3 TCL批量操作效果图

1.3代码剖析

指定圆心位置和半径的功能块。

*clearmark elem 1
*clearmark node 1
#x0 y0 z0 为圆心坐标，rad为孔半径
set y0 0.0
set z0 0.0
set rad 20.0
set h 50
set elemSize 0
set r1 [expr $rad+$elemSize]
*nodemodify 1 0 0 0

选出半径范围内的所有节点，然后同节点选出单元，并获取单元编号。

*createmark node 1 "by cylinder" $x0 $y0 $z0 0 0 1 [expr $rad+$elemSize] $h "inside" 1 1 0
*isolateentitybymark 1
*findmark nodes 1 257 1 elements 0 1
set elemList [hm_getvalue elems mark=1 dataname=id]

然后是一顿猛如虎的操作：删除单元，投影节点。代码如下：

foreach elemId $elemList {
    set count 0
    set ndId [hm_getvalue elems id=$elemId dataname=nodes]
    foreach i $ndId {
        foreach { dis x y z } [hm_getdistance nodes 1 $i 0] {}
        if {[expr $dis-$rad] <= 0.000001} {
            incr count 
        } 
    }    
    if {$count ==4} {
        *createmark elements 1 $elemId
        *deletemark elements 1
    } elseif {$count ==3}  {
        foreach i $ndId {
            foreach { dis x y z } [hm_getdistance nodes 1 $i 0] {}
            if {[expr $dis-$rad] >= 0.000001} {
                *nodemodify $i [expr $rad/$dis*($x)+$x0] [expr $rad/$dis*($y)+$y0] [expr $rad/$dis*($z)+$z0] 
            } 
        }
        *createmark elements 1 $elemId
        *deletemark elements 1
    } elseif {$count ==1} {
        foreach i $ndId {
            foreach { dis x y z } [hm_getdistance nodes 1 $i 0] {}
            if {[expr $dis-$rad] <= 0.000001} {
                *nodemodify $i [expr $rad/$dis*($x)+$x0] [expr $rad/$dis*($y)+$y0] [expr $rad/$dis*($z)+$z0] 
            }
        }
    } else {        
        set distList []
        foreach i $ndId {
            foreach { dis x y z } [hm_getdistance nodes 1 $i 0] {}
            set t1 [expr abs($dis-$rad)]
            lappend distList [list $i $t1]
        }
        set distListSorted [lsort -real -index 1 $distList]
        set nd1 [lindex $distListSorted 0 0] 
        set nd2 [lindex $distListSorted 1 0] 
        foreach { dis x y z } [hm_getdistance nodes 1 $nd1 0] {}
        *nodemodify $nd1 [expr $rad/$dis*($x)+$x0] [expr $rad/$dis*($y)+$y0] [expr $rad/$dis*($z)+$z0] 
        foreach { dis x y z } [hm_getdistance nodes 1 $nd2 0] {}
        *nodemodify $nd2 [expr $rad/$dis*($x)+$x0] [expr $rad/$dis*($y)+$y0] [expr $rad/$dis*($z)+$z0]
    }
}