【讨论】处女之作，出了问题！

hxzdavid

学习abaqus两月左右，终于写了一个自己的.inp文件，导入cae中分析却出了问题。error是：406 elements have missing property definitions
下面是源代码：
请大牛指出哪里出了问题，小生这相有礼了！

chenhu

修改了一下，你问题出在你定义的是平面应力单元，应该使用solid section去定义单元属性，而你却使用壳单元的shell section去定义了，另外，你nfill中15应为14，这表示两线之间共分为多少份 ！
*HEADING  
THE PLANAR(TWO DIMENSIONAL PROBLEM)  
UNITSENGTH-MM FORCE-N STRESS-N/MM2  
**  
**THE DEFINITION OF NODE  
**  
**LEFT EDGE  
*NODE  
1,0,0  
30,30,0  
*NGEN,NSET=BOTTOMEDGE  
1,30,1  
**RIGHT EDGE  
*NODE  
421,0,30  
450,30,30  
*NGEN,NSET=TOPEDGE  
421,450,1  
*NFILL  
BOTTOMEDGE,TOPEDGE,14,30  
**  
*NSET,NSET=LEFTEDGE,GENERATE  
1,421,30  
*NSET,NSET=RIGHTEDGE,GENERATE  
30,450,30  
**  
**THE DEFINITION OF THE ELEMENT  
**  
*ELEMENT,TYPE=CPS4I  
1,1,2,32,31  
*ELGEN,ELSET=ALL  
1,29,1,1,14,30,29  
*ELSET,ELSET=RIGHTELEMENT,GENERATE  
29,406,29  
**  
*SOLID SECTION,ELSET=ALL,MATERIAL=M1
4,  
**  
**THE DEFINITION OF THE MATERIAL  
*MATERIAL,NAME=M1  
*ELASTIC  
2.2e4,0.3  
**  
*BOUNDARY  
LEFTEDGE,ENCASTRE  
*BOUNDARY,TYPE=DISPLACEMENT  
BOTTOMEDGE,2  
**  
**STEP  
*STEP,NAME=STRETCHING  
*STATIC  
0.01,1,1.E-5  
*DLOAD  
RIGHTELEMENT,P2,-1000000  
**  
**OUTPUT  
*NODE PRINT,FREQ=1,NSET=TOPEDGE  
U,  
V,  
*EL PRINT,FREQ=1  
S,  
E,  
*END STEP  
  
下图是计算后的mises应力分布！

hxzdavid

问题解决了，多谢多谢！

ilovenili

hxzdavid 兄弟，告诉我你是怎么学写input文件的
我也是学了两个月，还是不会写，可能是没有掌握
方法，拜托了！！

hxzdavid

我记得论坛上曾经贴过输入文件的模板，找出来看看，很有用处！

hxzdavid

我也是新手，很多问题同样不懂！

milo

  
个人观点：
        关于网格生成部分就不要花太多的时间去手工写
输入文件了，采用合适的前处理器可以省很多时间。重点
应该化在如何定义history data上。
        

chenhu

但有时手写网格对于你了解模型及方便的进行单元选择以便加载、约束等也有一定作用，对于规则、尤其是有规律的模型网格分布个人感觉还是手写的好！
当然，重中之重也是history data，这不仅对控制加载、分析步能较好控制，而且还可以更方便的控制变量输出以控制结果文件的大小！
找几个例子好好琢磨一下，并不只是看，而且要自己试着写，从头至尾的写，多看看keyword帮助，首先要保证你对一般的inp能看懂，不懂的keyword要查，搞清楚它的意思及用法，这样才能对你自己写打下基础！
关于模型方面，我的作法就是结点及单元写完一点就倒入CAe下看看对不对，不对的尽快改，这样效率会高一点！结点单元也不是一个一个的输，而要掌握分布规律，多用*ncopy、*ngen、*nfill、*elgen来完成结点、单元的批量生成！
关于history data，建议在分析之前针对分析的类型看看analysis manual相关内容，这其中要用到什么*keyword它一般都会告诉你！这样，作到心中有数，再试着写写，看看相关例子对照，有可能的话，还可以“抄”一下例子的写法，为已作用！
另外，如果你还不太明白inp的布局的话，多看看abaqus/standard get started with keyword与with interactive两个手册对照着看，收获会很大（庄老师的那本中文翻译就是老版本abaqus的帮助手册，有些还能用).
熟悉CAE操作的话，可以先在CAE下完整的作几个分析，将它的write input文件打开看看，除了模型的不看外，其它的加载等还是可以参考的，不过困难的就是它的各种对part 、set等取的名字有点让人头晕，当然仔细看还是可以看明白 的！
  
我想这样学下去，简单一点还是可以写出来的！我也是9月初才接触abaqus的，所以也只用了点皮毛，大家相互交流！也请大牛讲一讲！:）

seansheu

milo wrote:
   
个人观点：
        关于网格生成部分就不要花太多的时间去手工写
输入文件了，采用合适的前处理器可以省很多时间。重点
应该化在如何定义history data上。
        

  
milo兄說的好，的確網格前處理可交由專業之前處理軟件來做，
心力專注在history data上就好了...
  
不過，我也提供一個小心得，供大家參考...
在劃分網格時，可就分析的需求，或是後處理欲觀察某些網格的需求，
能多設SET就多設一些，無論是NODE亦或是ELEMENT，
此外大SET還可許多包含小SET，彈性變化大，這可是ABAQUS的一大特色，端視你個人自己的變化智慧，
如此一來，於history data的定義時你會發現方便許多，
後處理時，要找出local的一些data亦會快速方便

seansheu

不知chenhu兄專業是在那個領域？
其實手寫網格對一些簡單的例子，或是學理驗證簡單模型可能遠較前處理軟件快速，
  
不過就以小弟經驗來說，因多接觸電子業或機械業的委託顧問案，
手上拿到的多是複雜的CAD檔案，單一零件有時還好，若是還包含組裝那可就頭痛了，
再加上模型通常有加上許多奇怪的特徵，例如：肋、導圓角、曲面...
這時可就非得求助專業的前處理軟體，要不然無法處理。
  
通常，網格劃分可分為Top-Down和Bottom-Up兩種方式，
以Geometry-base的方式去切割體積，然後做網格，這屬Top-Down，ABAQUS/CAE屬此類，
若是先將Geometry各項特徵投影至一個面，然後做2d mesh，再抽拉成3d mesh，這則屬Bottom-Up
ABAQUS/CAE很可惜無法用Bottom-UP做網格，但ANSYS就兩種都行，
若是專業之前處理軟件，如：HyperMesh、I-Deas 或是 Patran 就都可以兩種方式處理網格，
重點是，端視你欲處理的模型，再決定該採用Top-Down或是Bottom-Up何者較有效率。
  
個人一點小心得，供大家參考...

chenhu

老大高见！足有受益！
我所说的也仅限于初学者（因为我就是初学者），所以简单的还能对付，在老大面前就不足提了！
我本人是化工机械出身的，化工设备不像动力机械那样的精密或复杂，所以手写也还不算慢，这样倒也省了学其它前处理软件的时间了！：）
老大说的Top-Down和Bottom-Up很重要，我这阵子就对CAE唯一感到不爽的也就是它由面网生成体网这一功能缺陷，对于版里的那个三维裂纹网格的划分要是能用这个功能就会很好的解决！

seansheu

呵~~~我那經驗也只是一些前人教授於我，現在再轉述罷了，不足掛齒矣...
不過，不管方式為何，只要能快速解決問題（當然答案要對），這才是重點。
正所謂：黑貓白貓，能捉耗子就是好貓
所以大家也就不必太拘泥於前處理的形式了!
  
BTW,我覺得ABAQUS若要好好推動它CAE的marketing，Bottom-Up的處理方式應該是未來版本必要的功能

chenhu

seansheu wrote:
呵~~~我那經驗也只是一些前人教授於我，現在再轉述罷了，不足掛齒矣...
不過，不管方式為何，只要能快速解決問題（當然答案要對），這才是重點。
正所謂：黑貓白貓，能捉耗子就是好貓
所以大家也就不必太拘泥於前處理的形式了!
  
BTW,我覺得ABAQUS若要好好推動它CAE的marketing，Bottom-Up的處理方式應該是未來版本必要的功能

是呀，我也这么想的，这么好的一种分网思路它就没有想到？它和Altair公司的合作应该可以借鉴一下吧！它的CAE要是再强点，那现在ansys在中国所占有的市场可能就有点危险了！：（（开个玩笑，各有各有用法）
但自下而上的分网功能的实现我想涉及到的问题就是当up网格生成后bottom的那些临时网格怎么处理到它？要知道，abaqus对单元、结点的要求是较高的！

seansheu

chenhu wrote:
[quote]seansheu wrote:
呵~~~我那經驗也只是一些前人教授於我，現在再轉述罷了，不足掛齒矣...
不過，不管方式為何，只要能快速解決問題（當然答案要對），這才是重點。
正所謂：黑貓白貓，能捉耗子就是好貓
所以大家也就不必太拘泥於前處理的形式了!
  
BTW,我覺得ABAQUS若要好好推動它CAE的marketing，Bottom-Up的處理方式應該是未來版本必要的功能

是呀，我也这么想的，这么好的一种分网思路它就没有想到？它和Altair公司的合作应该可以借鉴一下吧！它的CAE要是再强点，那现在ansys在中国所占有的市场可能就有点危险了！：（（开个玩笑，各有各有用法）
但自下而上的分网功能的实现我想涉及到的问题就是当up网格生成后bottom的那些临时网格怎么处理到它？要知道，abaqus对单元、结点的要求是较高的！ [/quote]
  
關於你這問題，很簡單...
首先，2d mesh品質一定要達要求，如此3d mesh才會抽拉成好的質，
然後，在HM、I-DEAS或是ANSYS，當抽拉完3d mesh後，就把2d mesh全砍掉就好啦
砍掉那些2d mesh是不會影響到那最底層3d mesh的網格節點的
這在軟體內他會自動判斷，若是節點屬於3d element，那麼它不會讓節點被砍掉
以I-deas為例，它若是先砍掉 3D element，那麼node還會留下來，所以還要再砍一次節點
  
通常，抽拉3D網格完後，最好再做一下共點的檢查，讓重複的點做merge動作，
並且check一下是否有重複的element，然後把重複的給清掉...
  
我以前在HM 2.0b版做網格時都需要做這些動作，現在的HM 6.0還沒深入摸透，我就不知是否還需要如此麻煩了...

chenhu

seansheu wrote:
[quote]chenhu wrote:
[quote]seansheu wrote:
呵~~~我那經驗也只是一些前人教授於我，現在再轉述罷了，不足掛齒矣...
不過，不管方式為何，只要能快速解決問題（當然答案要對），這才是重點。
正所謂：黑貓白貓，能捉耗子就是好貓
所以大家也就不必太拘泥於前處理的形式了!
  
BTW,我覺得ABAQUS若要好好推動它CAE的marketing，Bottom-Up的處理方式應該是未來版本必要的功能

是呀，我也这么想的，这么好的一种分网思路它就没有想到？它和Altair公司的合作应该可以借鉴一下吧！它的CAE要是再强点，那现在ansys在中国所占有的市场可能就有点危险了！：（（开个玩笑，各有各有用法）
但自下而上的分网功能的实现我想涉及到的问题就是当up网格生成后bottom的那些临时网格怎么处理到它？要知道，abaqus对单元、结点的要求是较高的！ [/quote]
  
關於你這問題，很簡單...
首先，2d mesh品質一定要達要求，如此3d mesh才會抽拉成好的質，
然後，在HM、I-DEAS或是ANSYS，當抽拉完3d mesh後，就把2d mesh全砍掉就好啦
砍掉那些2d mesh是不會影響到那最底層3d mesh的網格節點的
這在軟體內他會自動判斷，若是節點屬於3d element，那麼它不會讓節點被砍掉
以I-deas為例，它若是先砍掉 3D element，那麼node還會留下來，所以還要再砍一次節點
  
通常，抽拉3D網格完後，最好再做一下共點的檢查，讓重複的點做merge動作，
並且check一下是否有重複的element，然後把重複的給清掉...
  
我以前在HM 2.0b版做網格時都需要做這些動作，現在的HM 6.0還沒深入摸透，我就不知是否還需要如此麻煩了... [/quote]
  
那我们就看CAE未来能不能实现这样的功能吧！感觉seansheu大侠真的是大牛呀！不知自己啥时候才能有这水平呀！
另外，CAE下merge结点的功能 在哪，我好像没找到？
另外，CAE下对网格质量的评价好像也没有吧，其只能靠自己来判定了，这好像就不如hm那样有一套完整全面的判定标准了，不过，话也说过了，软件的侧重点不一样呀！

seansheu

chenhu wrote:
那我们就看CAE未来能不能实现这样的功能吧！感觉seansheu大侠真的是大牛呀！不知自己啥时候才能有这水平呀！
另外，CAE下merge结点的功能 在哪，我好像没找到？
另外，CAE下对网格质量的评价好像也没有吧，其只能靠自己来判定了，这好像就不如hm那样有一套完整全面的判定标准了，不过，话也说过了，软件的侧重点不一样呀！

  
ABAQUS/CAE下是沒有merge的功能，我指的這是對專業前處理軟件而言，
另外，CAE也是有對網格質量評價的功能，你可在mesh那模組中找到一個打勾的按鈕，
那就是評價網格質量的工具。
  
其實，我只能算是軟件奴隸，軟件功能懂些皮毛，熟用罷了，
FEM/CAE工作最重要的還是學理的了解,要不然可是會GIGO...
學理上的應用我可是還差wjytj,cafe0519遠的多呢
我學習的路還很長,大家共勉吧

ghostfql

几位大牛真是厉害！差的很远呀，需要努力了！对了，自己手写inp文件是不是一定是必须的呢？我觉得让cae自己写inp，然后再改动一下，加入一些cae无法直接嵌入的东西，不是很好吗？为什莫要花那末多时间学习写inp呢？请执教，谢谢

hxzdavid

都说得挺有道理，那要看你达到什么目的！
如果自己面对的工作就是与abaqus打交道，为了精通软件，当然什么都得做；如果是做毕业课题，也许以后根本用不上，化那么多时间在input上，得不尝失，也没有这个必要。

yzf1111

建议模型用尽可能简单的方式如HM（偶很喜欢的说），或ghostfql说的方法，不要在这上面上浪费太多时间，History　data 手写，有助于加深对问题的理解，一家之言呵呵

ghostfql

是亚，我自己觉得应该吧更多时间花在user subroutine上，如果要做研究的话，不是吗？可能大牛有好的意见，望批评指正，洗耳恭听呀！