有上下封底的圆柱壳如何且要求上下封底是蝴蝶网格，该如何画？

icola

看了timothy1985兄的代码，可以画出上下底面为蝴蝶网格的实心圆柱，要是画上下底面为蝴蝶网格的圆柱壳该如何画呢？以下为timothy1985兄的代码：

1. 
   
2. block 1 11 21;1 11 21;1 11; 0 1 1; 0 1 1; 0 1;
   
3. dei 2 3;2 3;;
   
4. pb 2 2 1 2 2 2 xy .9 .9;
   
5. sfi 1 -3;1 -3;;zcy 2;
   
6. lct 3 rz 90; repe 3;
   
7. lrep 0:3;
   
8. endpart
   

复制代码

freebeita

block 1 0 21;1 0 21;.....
中间的0表示建立壳体

nickcoom

楼主也不是壳，画一个，我怎么不喜欢用pb呢，不熟悉，lz可以讲讲不。有啥好处……

nickcoom

材料不一样的，想知道pb的作用啊……

icola

nickcoom，我就是不知道怎么画壳才问的啊，你倒是画了一个，把代码贴上来让兄弟我学习学习呗！

icola

请教一下，你是怎么画出的壳？ 3# nickcoom

icola

pb的作用是调节节点坐标
pb 2 2 2 2 2 1 xy .9 .9表示将 i 索引范围2:2, j 索引范围2:2, k 索引范围1:2 的点的x和y坐标调节为0.9和0.9，而z坐标不变。 4# nickcoom

icola

不对啊，这样做会出现错误。 2# freebeita

nickcoom

7# icola 我的意思，它调节那个有什么左右，可以直接投影也是这种效果，我试过的。一直不知道他的作用，命令流在笔记本里面，晚上贴吧……

icola

不要pb也确实可产生蝴蝶网格，此处的作用是调节与中间矩形相接网格的内角，使之更接近90度，以提高网格质量。 9# nickcoom

freebeita

只能说pb可以更好的调整中间部分的网格节点，使图形更接近正交，但最后都要光滑处理，tme或者esm都行

icola

TG 2.2 Tutorial 的第60页讲了利用pb、tricent和relax进行网格优化的问题。可进行参考。第一章图为优化前，第二张图为优化后。

nickcoom

楼上说的没有错，如果你本身就定义那个block 1 11 21;1 11 21;1 11; 0 1 1; 0 1 1; 0 1;我觉得没有必要再用pb了，这个里面本来就是一样的坐标，你pb一下，有啥用呢，所以我才这么问的……不明白，那个杯子是有gap才，才pb在一起的。没有回去，我现在直接在你的那个命令上改下吧。这么发楼主那样的附近不会啊

nickcoom

楼上说的没有错，如果你本身就定义那个block 1 11 21;1 11 21;1 11; 0 1 1; 0 1 1; 0 1;我觉得没有必要再用pb了，这个里面本来就是一样的坐标，你pb一下，有啥用呢，所以我才这么问的……不明白，那个杯子是有gap才，才pb在一起的。没有回去，我现在直接在你的那个命令上改下吧。这么发楼主那样的附近不会啊

nickcoom

如何上传lz那样的命令流啊……求解……

icola

谢谢，命令流是按照图中的方法粘贴的。 14# nickcoom

nickcoom

谢谢了，真是细心啊……很高兴，tg学的人越来越多啊

mikebaby

1. cylinder -1;1 17;1 10;1;0 360;0 10;
   
2. block 1 3 7 9;1 3 7 9;-1;-0.5 -0.5 0.5 0.5;-0.5 -0.5 0.5 0.5;0;
   
3. dei 1 2 0 3 4;1 2 0 3 4;;
   
4. sd 1 cy 0 0 0 0 0 1 1;
   
5. sfi -1 0 -4;;-1;sd 1;
   
6. sfi ;-1 0 -4;-1;sd 1;
   
7. lct 1 mz 10;;
   
8. lrep 0 1;
   
9. endpart
   
10. merge

复制代码

我也发一个，初学乍练，请多指教！

nickcoom

楼上用了俩个块啊，哈哈，本来我也是这个思路的……

icola

用你的命令流确实可以生成柱壳，想用nseti将该柱壳上的所有节点加入一个名为aaa的节点集，用eseti命令将柱壳上的所有单元加入一个名为bbb的单元集，然后输出为NE/NASTRAN格式的网格文件，但网格文件中aaa节点集中的节点编号不是从1到300，而是从171到300；bbb单元集中的节点编号也不是从1到240，而是从145到240，这是为何？
定义节点和单元集并输出为NE/NASTRAN格式网格文件的命令流如下：

1. c ****输出为NE/NASTRAN格式****
   
2. ne/nastran
   
3. cylinder -1;1 17;1 10;1;0
   
4. 360;0 10;
   
5. block 1 3 7 9;1 3 7 9;-1;-0.5 -0.5 0.5 0.5;-0.5 -0.5 0.5
   
6. 0.5;0;
   
7. dei 1 2 0 3 4;1 2 0 3 4;;
   
8. sd 1 cy 0 0 0 0 0 1 1;
   
9. sfi -1 0
   
10. -4;;-1;sd 1;
    
11. sfi ;-1 0 -4;-1;sd 1;
    
12. lct 1 mz 10;;
    
13. lrep 0 1;
    
14. c
    
15. ***定义名为aaa的节点集合名为bbb的单元集***
    
16. nseti ;;; = aaa ;
    
17. eseti ;;; = bbb ;
    
18. endpart
    
19. merge
    
20. write

复制代码

输出的网格文件中的节点集和单元集的内容为：
$
$----------     NODE SETS     ----------
$
$ Node set aaa     
SET 1 = 171 THRU 300
$ Element set bbb     
SET 2 = 145 THRU 240

18# mikebaby