map solution（mesh-to-mesh solution mapping）中遇到的问题

drunkalcohol

背景：使用standard求解器，模型分析中，出现网格扭曲过大的提示，计算也不能收敛。

思路：网格扭曲过大是由于模型大变形引起的，所以打算利用map solution（mesh-to-mesh solution mapping）在step间网格重画，克服扭曲过大的问题。

难题：使用map solution（mesh-to-mesh solution mapping），第一步是提取出变形后的几何体。虽然可以通过先import ODB得到几何体，但是所得到的是网格几何体，而Abaqus/CAE并不支持网格几何体再划分网格，除非通过第三方的例如HM，但是我恰恰不会HM。也在网上看到通过自己写脚本，利用PartFromOdb得到几何体的方法，但是我在Abaqus/CAE中用命令行搞了半天，也搞不定，而且网上的示例是二维的，我要得到三维零件的外形。
    有懂得在Abaqus/CAE中提取变形后的零件几何体的（非网格几何体）帮帮忙，给点指导。如果能给出step by step的操作指南，那就更美妙了。感激不尽。

1147haoran

whskycom 发表于 2012-4-5 07:28
现举一个简单例子，如下：
例子：有一个初始为1*1*1的立方体，step1结束后此立方体变为1*2*0.5，step2结 ...

1) 不必相等
2）我想应该是可以的，毕竟最终ABAQUS 所需要的只是inp文件，之前的过程类似“black box”吧，
3）参考用户手册：
Solution mapping operates by interpolating results from nodes in the old mesh to points (either nodes or integration points) in the new mesh. The first step, therefore, involves associating solution variables with nodes in the old mesh. For nodal solution variables, such as nodal temperature or pore pressure, the association is already made. For integration point variables Abaqus obtains the solution variables at the nodes of the old mesh by extrapolating values from the integration points to the nodes of each element and then averaging these values over all similar elements abutting each node.

Next, the location of each point in the new mesh is obtained with respect to the old mesh. The new mesh points include integration points in all cases and nodes in procedures that record nodal state in addition to displacements (for example, nodal temperatures in coupled temperature-displacement procedures).
The element (in the old mesh) in which the point lies is found, and the point's location in that element is obtained. (This procedure assumes that all points in the new mesh lie within the bounds of the old mesh: warning messages are issued if this is not so, and the values of the variables are set to zero.)
The variables are then interpolated from the nodes of the old element to the points in the new model.
All necessary variables are interpolated automatically in this way so that the solution can proceed with the new mesh.

drunkalcohol

wcjhit110 发表于 2011-9-1 11:04
我是搞大变形的，也想用一下Map solution。麻烦楼主能否把网上关于 “ 在网上看到通过自己写脚本，利用Part ...

估计对你没有太多帮助

drunkalcohol

没人帮忙么？帖子都沉了

敦诚

你要是不会hm也没有关系，hm中就有个功能是把网格模型转化为几何模型的，应用也比较简单，你可以试着学一学，然后再对于新生成的几何模型在abaqus里面划分来完成map solution。

JingheSu

回复 1# drunkalcohol
如果是3维的,还是想想其他办法吧.
好像那个python命令只能得到网格几何体(PartFromOdb:This method creates an orphan mesh Part object by reading an output database)
而且似乎只有这一个命令可以得到非网格几何partart2DGeomFrom2DMesh,但仅限于2Dmesh.

drunkalcohol

回复 3# 敦诚
我还是去搞个HM研究下吧。
不过兄弟有HM10.0（*破*解*版*）的下载链接不？网上关于这个的好像被清理过了。

drunkalcohol

回复 4# JingheSu
好像还真是，用户手册里的例子全是二维的

wcjhit110

我是搞大变形的，也想用一下Map solution。麻烦楼主能否把网上关于 “ 在网上看到通过自己写脚本，利用PartFromOdb得到几何体的方法，网上的示例是二维的”发给连接给我或者脚本及用法！万分感谢！wcjhit@163.com

wcjhit110

drunkalcohol 发表于 2011-9-19 20:16
估计对你没有太多帮助

非常感谢！我研究一下！

aqshouhe

好像可以先编辑成一个带网格变形3D模型，然后将这个模型倒入三维画图软件中，保留轮廓线，重新生成几何图形。至于怎么把前面的计算结果导进来继续后续分析，就不会了

whskycom

敦诚 发表于 2011-8-8 21:33
你要是不会hm也没有关系，hm中就有个功能是把网格模型转化为几何模型的，应用也比较简单，你可以试着学一学 ...

现举一个简单例子，如下：
例子：有一个初始为1*1*1的立方体，step1结束后此立方体变为1*2*0.5，step2结束后变为1*5*0.2。假设由于变形过大，step1结束后我将变形后的网格提取出来并重新划分，然后再重新导入abaqus中计算。几个问题如下：
1）重新划分的网格数和原网格数是否必须相等？
2）如果step1结束后，立方体还是立方体（只是长宽高有所变化）即：变形比较规则，那么有可否不进行网格提取，而直接新建一个立方体并划分网格？
3）最重要的问题：重新划分后的网格（或者inp文件）与step1产生的.odb文件映射（或关联）的机制（或原理，不是流程）是什么，怎么就能将原来的.odb文件重新加入到.inp中呢？是通过节点编号吗？还是通过instance名称？还是通过什么其他方式？
非常感谢