。。。

刚刚给罗马尼亚小妞讲模拟，居然给她讲困了。

mein scheisse Deutsch!

炮！

Z： 我看到了你的结婚照。

M： 嘿嘿，美女吧！

Z： 哈哈。。。

M： 都不像了！

Z： 那必须的啊！（我回话之前，没看到炮写的那句“都不像了！”）

Z： 你是我见过得最漂亮的炮！

M： 那是啊，咱底板儿好！

M： 哈哈。。。

星期五的晚上，心情大好，拍夜景！

这两张照片是好看啊？是好看啊？还是好看啊？

巴黎印象

今天心情不错，找到了曾经的照片，传上来些！

今天两件事儿！

上午接到旅行社的电话，告诉我回国的航班要推迟一天！

火了！一共只有那么两天假期，还要在德国干靠三天。国内的机票我都订好了啊，还要改日期。索性给汉莎打了电话，原来原来的航班取消了，我可以选择提前一天或者推迟一天的航班。嘿嘿，坏事儿变好事儿，我可以早回家一天了，周末告诉老妈。不过国内的航班还是要改的，不爽。。。需要跟moment取得联系！

下午接到研究所的通知，过两天我们部门要来一个罗马尼亚女学生，我跟小L要带着她做论文。

论文是关于有限元模拟的，分明就是我自己一个人带嘛。小L却非要过来分一杯羹，还是他妈的一大杯。

无论如何，我也要带学生了，罗马尼亚小妞儿！

^_^

卓越拓扑之路----锥体的结构化网格 II

昨天晚上下班，一辆奥迪的Q7强行左转，结果划到了刚刚驶过的U-Bahn身上，我坐在U-Bahn上。Q7身上被狠狠的画出了一道深沟，真是暴殄天物！不过这道深沟倒是给了我一点启发。在家里，我同样把三棱锥肢解了。

拦腰切一刀，三棱锥变成了两个三棱柱。到此为止，拓扑也不再有难度了。
重新建立一个BLOCK，也拦腰切一刀！

把棱柱相接面上的节点关联，两边的节点合并，再次关联！现在BLOCK跟三棱锥一模一样了！

用两个y-block，再次把棱柱劈开！

拓扑关系至此完成！剩下的，就是定义每条edge上面的单元数量了。农民都会干！

清晰的拓扑概念是完美离散品质保证！


神奇的金字塔！

感谢那部无私的Q7！

卓越拓扑之路----锥体的结构化网格 I

说起来搞有限元已经有段时间了，深知空间离散对于仿真的重要性，在这方面下了不少功夫，不断的追求更高的网格品质。有些几何体看起来是那么简单，但是画起网格来却总是让我崩溃。锥体就是个典型！当然，对于简单的集合体，非结构化的网格不在考虑范围之内。

今天有些时间，干脆就研究研究三棱锥。













常规做法，先做CAD，然后建立一个BLOCK包围离散集合体。



















对于三棱锥而言，底面只有三个节点，所以需要对BLOCK其中一面上的节点进行合并。为了保证网格的连贯性，索性我把相对应的BLOCK上表面的节点合并了。




















把相应的节点跟几何体相关联以后，再用一个简单得Y-BLOCK，差不多就可以设定网格的大小了。

记得我再做渗透分析的时候，模型的每条边上我都定义了很多的网格，以便计算收敛的好些，模拟的更贴近实际，结果看到了千丝万缕的空间，如同万箭穿心一般。为了不至于眼花缭乱，对于这个集合体，我只定义了有限几个单元。

现在问题已经很简单了，锥体上端只有一个节点，只有BLOCK上表面的所有节点都合并起来，问题似乎就应该解决了。

上面的网格，看似漂亮，结构清晰明了，质量很高。但是没用，最后的一步是一直以来困扰我的难题。虽然，我可以把节点合并，却最终改变了最顶端网格的性质，违背了结构化网格的基本定义。冥思苦想一上午，还是没有想到好地解决办法。

难道锥体很难建立高质量的结构化空间离散吗？

未完，待续。。。