首页 > 其他分享 >【Abaqus】材料行为的非均匀空间分布

【Abaqus】材料行为的非均匀空间分布

时间：2023-08-02 13:45:48浏览次数：38  
标签：10 材料 Abaqus Material Elset FourElem 12.5 均匀 MODULUS

设想一种情况:在有限元分析中，一个区域或者整个网格中，每个单元的材料行为都是单独的。这时在ABAQUS中应该如何设置？

两种办法：

  • 给每个单元创建一个集合，然后一一赋予SECTION.
  • 使用*Distribution关键字，实现空间分布的材料行为，再将SECTION属性赋予给单元。这一种方法好处是，减少后处理的压力，没有那么多的SET和SECTION

两种方法对比：

  • 第一种方法可以在CAE界面中手动设置，也可以用Python脚本自动设置，但是一旦单元数量过多，就会导致CAE界面卡顿，且Python速度也不是很快。电脑性能不行的话，后处理压力比较大，因为有很多的Set,Section.
  • 第二种方法，不支持CAE，只能修改INP文件实现。优点是后处理软件压力小点，没有那么多set

下面以下面的四个实体单元为例子，分别使用两种方法，直接编写INP文件，命令行中提交计算，进行静态分析。

img

方法一

*Heading
** JOB:OneByOneSet_AnisotropicMaterial_FourELEM
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=FourElem
*Node
      1,         -25.,       -11.25,          10.
      2,         -25.,       -0.625,          10.
      3,         -25.,          10.,          10.
      4,         -25.,       -11.25,           0.
      5,         -25.,       -0.625,           0.
      6,         -25.,          10.,           0.
      7,        -12.5,       -11.25,          10.
      8,        -12.5,       -0.625,          10.
      9,        -12.5,          10.,          10.
     10,        -12.5,       -11.25,           0.
     11,        -12.5,       -0.625,           0.
     12,        -12.5,          10.,           0.
     13,           0.,       -11.25,          10.
     14,           0.,       -0.625,          10.
     15,           0.,          10.,          10.
     16,           0.,       -11.25,           0.
     17,           0.,       -0.625,           0.
     18,           0.,          10.,           0.
*Element, type=C3D8R
1,  7,  8, 11, 10,  1,  2,  5,  4
2,  8,  9, 12, 11,  2,  3,  6,  5
3, 13, 14, 17, 16,  7,  8, 11, 10
4, 14, 15, 18, 17,  8,  9, 12, 11
*Elset, elset=E1
1
*Elset,elset=E2
2
*Elset,elset=E3
3
*Elset,elset=E4
4
**Gobal orien
*Orientation, name=Ori-1
1., 0., 0., 0., 1., 0.
1, 0.
**Section Assign
*Solid Section,Elset=E1,Material=M1,Orientation=Ori-1
*Solid Section,Elset=E2,Material=M2,Orientation=Ori-1
*Solid Section,Elset=E3,Material=M3,Orientation=Ori-1
*Solid Section,Elset=E4,Material=M4,Orientation=Ori-1
*End Part
** ASSEMBLY
**
*Assembly, name=Assembly
**  
*Instance, name=FourElem-1, part=FourElem
*End Instance
**  
*Nset, nset=Set-1, instance=FourElem-1, generate
 13,  18,   1
*Elset, elset=Set-1, instance=FourElem-1
 3, 4
*Elset, elset=_Surf-1_S2, internal, instance=FourElem-1
 1, 2
*Surface, type=ELEMENT, name=Surf-1
_Surf-1_S2, S2
*End Assembly
**Material:M1
*Material,Name=M1
*Elastic,TYPE=ANISOTROPIC
8929.4608,1526.5216,8923.6134,3635.3632,3476.3026,120679.2644,8.3944,7.7689,
439.6758,3701.6801,-145.4016,-135.4002,-7614.9276,-29.5989,4212.6626,-122.1556,
-112.4105,-6361.0837,-24.5995,433.279,4055.9163
**Material:M2
*Material,Name=M2
*Elastic,TYPE=ANISOTROPIC
17550.0133,13611.2277,25852.2543,13021.7401,17628.6822,24242.5311,9788.8481,13699.1544,
13033.5732,14762.1216,-9538.4395,-13349.2807,-12699.5968,-10779.3782,14203.3879,-11372.6959,
-15914.7528,-15140.8525,-12851.6274,12522.8552,18629.8078
**Material:M3
*Material,Name=M3
*Elastic,TYPE=ANISOTROPIC
8929.8392,1529.3591,8928.7933,3642.9557,3613.164,120481.7162,9.95,9.8574,
517.3742,3702.3839,-147.4303,-146.0569,-7669.087,-35.2767,4222.4903,-143.1982,
-141.2357,-7414.3551,-34.4302,509.844,4193.1319
**Material:M4
*Material,Name=M4
*Elastic,TYPE=ANISOTROPIC
9058.1198,1940.0409,10089.9165,5393.1754,11874.9043,98667.1446,203.4878,541.258,
4678.4804,3943.7999,639.8814,1702.0259,14711.8024,766.6464,6110.7723,1221.6736,
3249.5219,28087.9378,1463.6933,4602.6886,12487.
** 
** STEP: Step-1
** 
*Step, name=Step-1, nlgeom=YES
*Static
0.1, 1., 1e-05, 0.2
** 
** BOUNDARY CONDITIONS
** 
** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary
Set-1, ENCASTRE
** 
** LOADS
** 
** Name: Load-1   Type: Pressure
*Dsload
Surf-1, P, 1.
** 
** OUTPUT REQUESTS
** 
*Restart, write, frequency=0
** 
** FIELD OUTPUT: F-Output-1
** 
*Output, field
*Node Output
U, 
*Element Output, directions=YES
S, 
** 
** HISTORY OUTPUT: H-Output-1
** 
*Output, history, variable=PRESELECT
*End Step

方法二

*Heading
**JOB:DistributionAnisotropicMaterial_FourELEM
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=FourElem
*Node
      1,         -25.,       -11.25,          10.
      2,         -25.,       -0.625,          10.
      3,         -25.,          10.,          10.
      4,         -25.,       -11.25,           0.
      5,         -25.,       -0.625,           0.
      6,         -25.,          10.,           0.
      7,        -12.5,       -11.25,          10.
      8,        -12.5,       -0.625,          10.
      9,        -12.5,          10.,          10.
     10,        -12.5,       -11.25,           0.
     11,        -12.5,       -0.625,           0.
     12,        -12.5,          10.,           0.
     13,           0.,       -11.25,          10.
     14,           0.,       -0.625,          10.
     15,           0.,          10.,          10.
     16,           0.,       -11.25,           0.
     17,           0.,       -0.625,           0.
     18,           0.,          10.,           0.
**定义单元和单元类型
*Element, type=C3D8R
1,  7,  8, 11, 10,  1,  2,  5,  4
2,  8,  9, 12, 11,  2,  3,  6,  5
3, 13, 14, 17, 16,  7,  8, 11, 10
4, 14, 15, 18, 17,  8,  9, 12, 11
**创建单元集合E1
*Elset, elset=E1
1,2,3,4
**定义全局坐标系
*Orientation, name=Ori-1
1., 0., 0., 0., 1., 0.
1, 0.
**赋予截面给实体单元集E1，指定材料和材料局部坐标
*Solid Section,Elset=E1,Material=M1,Orientation=Ori-1
*End Part
*Assembly, name=Assembly
**  
*Instance, name=FourElem-1, part=FourElem
**在实例下定义材料空间分布
*Distribution,Name=dist1,LOCATION=ELEMENT,TABLE=tab1
 ,8929.4608,1526.5216,8923.6134,3635.3632,3476.3026,120679.2644,8.3944
7.7689,439.6758,3701.6801,-145.4016,-135.4002,-7614.9276,-29.5989,4212.6626
-122.1556,-112.4105,-6361.0837,-24.5995,433.279,4055.9163
1,8929.4608,1526.5216,8923.6134,3635.3632,3476.3026,120679.2644,8.3944
7.7689,439.6758,3701.6801,-145.4016,-135.4002,-7614.9276,-29.5989,4212.6626
-122.1556,-112.4105,-6361.0837,-24.5995,433.279,4055.9163
2,17550.0133,13611.2277,25852.2543,13021.7401,17628.6822,24242.5311,9788.8481
13699.1544,13033.5732,14762.1216,-9538.4395,-13349.2807,-12699.5968,-10779.3782,14203.3879
-11372.6959,-15914.7528,-15140.8525,-12851.6274,12522.8552,18629.8078
3,8929.8392,1529.3591,8928.7933,3642.9557,3613.164,120481.7162,9.95,
9.8574,517.3742,3702.3839,-147.4303,-146.0569,-7669.087,-35.2767,4222.4903
-143.1982,-141.2357,-7414.3551,-34.4302,509.844,4193.1319
4,9058.1198,1940.0409,10089.9165,5393.1754,11874.9043,98667.1446,203.4878,
541.258,4678.4804,3943.7999,639.8814,1702.0259,14711.8024,766.6464,6110.7723,
1221.6736,3249.5219,28087.9378,1463.6933,4602.6886,12487.
*End Instance
**  
*Nset, nset=Set-1, instance=FourElem-1, generate
 13,  18,   1
*Elset, elset=Set-1, instance=FourElem-1
 3, 4
*Elset, elset=_Surf-1_S2, internal, instance=FourElem-1
 1, 2
*Surface, type=ELEMENT, name=Surf-1
_Surf-1_S2, S2
*End Assembly
**确定分布表格中的数据格式
*Distribution Table,Name=tab1
MODULUS , MODULUS , MODULUS , MODULUS , MODULUS , MODULUS , MODULUS
MODULUS , MODULUS , MODULUS , MODULUS , MODULUS , MODULUS , MODULUS, MODULUS
MODULUS , MODULUS , MODULUS , MODULUS , MODULUS , MODULUS 
**Material:M1
*Material,Name=M1
**各向异性材料，空间分布按照dist1中给定
*Elastic,TYPE=ANISOTROPIC
FourElem-1.dist1
** ----------------------------------------------------------------
** 
** STEP: Step-1
** 
*Step, name=Step-1, nlgeom=YES
*Static
0.1, 1., 1e-05, 0.2
** 
** BOUNDARY CONDITIONS
** 
** Name: BC-1 Type: Symmetry/Antisymmetry/Encastre
*Boundary
Set-1, ENCASTRE
** 
** LOADS
** 
** Name: Load-1   Type: Pressure
*Dsload
Surf-1, P, 1.
** 
** OUTPUT REQUESTS
** 
*Restart, write, frequency=0
** 
** FIELD OUTPUT: F-Output-1
** 
*Output, field
*Node Output
U, 
*Element Output, directions=YES
S, 
** 
** HISTORY OUTPUT: H-Output-1
** 
*Output, history, variable=PRESELECT
*End Step

结果对比

结果一致

img

资料

2.8.1 Distribution definition
*ELASTIC
*SOLID SECTION
*SOLID SECTION
*DISTRIBUTION TABLE
1.2.1 Input syntax rules

标签：10,材料,Abaqus,Material,Elset,FourElem,12.5,均匀,MODULUS
From： https://www.cnblogs.com/aksoam/p/17600445.html

相关文章

  • ABAQUS 在长城汽车发动机开发中的应用
    作为FEA软件行业的领导者，ABAQUS在长城汽车股份有限公司动力研究院的结构强度、刚度等分析中起到了举足轻重的作用，为发动机的研发贡献了不可磨灭的力量。本文通过实例介绍Abaqus软件在发动机开发过程中的应用。 长城汽车股份有限公司动力研究院高举自主研发的旗帜，推出多......
  • 【专题】2023新能源电池材料发展概览报告PDF合集分享（附原数据表）
    全文链接：http://tecdat.cn/?p=33303原文出处：拓端数据部落公众号《新能源电池材料发展概览报告合集》综合分析了电池下游应用场景以及当前材料体系所面临的问题，重点关注了光伏、储能和氢能领域的电池产品，预测了相关电池材料体系的发展趋势和投资价值。报告合集涵盖了光伏HJT电池......
  • ABAQUS-循环对称条件的详解
    概括anlysisofmodelthatexhibitcyclicsymmetry循环对称分析技术用于Standard求解器。makesitpossibletoanalyzethebehaviorofa360°structurewith cyclicsymmetrybasedonamodelofarepetitivesector;可以在静态、准静态、热传递分析中确定循环......
  • Abaqus 中的步进、增量、迭代和尝试概念 硕迪科技
    Abaqus中的步进、增量、迭代和尝试等可能会在概念上让Abaqus初学者感到困惑。清楚地了解分析步骤、荷载增量和迭代之间的区别非常重要。在这篇文章中快速了解Abaqus步骤和增量迭代。在ABAQUS中，步进增量迭代是解决非线性问题的一种数值计算方法。这种方法通常用于模拟材料的非......
  • Abaqus添加初始缺陷
    主要介绍通过施加节点位移的方法步骤一：复制model，新建Step，staticlinearperturbationTools->AnalyticalField定义场函数,例如：Asin(pi*X/Lx)*sin(pi*Y/Ly)一个四边形区域，A是位移最大值Load->Uz=1,选择场函数，即在Z方向施加位移Model->EditKeywords****Outputrequests......
  • 第二届先进与新兴材料国际学术会议（AEM2023）
    第二届先进与新兴材料国际会议（AEM2023）将于2023年9月16-17日在中国湖北武汉举行。该会议每年由湖北省众科地质与环境技术服务中心承办，我们诚挚地邀请您将论文全文投稿至EI和SCI期刊。截稿日期：2023年9月1日接受/拒稿通知：投稿后1-2周收录检索：EI/Scopus在线投稿★主讲嘉宾0......
  • 镭拓科普塑料激光焊接机可否焊接两种不同材料
    编辑：镭拓激光在上一篇中，我们给大家分享了关于塑料激光焊接机能够焊接哪些材料的小知识，本篇我们接着探讨有关塑激光料焊接机的内容。塑料激光焊接机是一种引人关注的新技术，它可以实现两种不同材料的无缝连接。然而，很多人好奇塑料激光焊接机能否将两种不同的塑料材焊接在一起？塑料激光......
  • 半导体晶圆衬底抛光研磨材料哪个好用？
    半导体晶圆衬底抛光研磨是半导体工艺中至关重要的一环，对晶圆衬底的加工质量和性能起着决定性的影响。在半导体制造过程中，抛光研磨工序是必不可少的，它能够平整晶圆衬底表面、去除缺陷并提供所需的光洁度和精度。在进行半导体晶圆衬底抛光研磨时，有几个关键要点需要注意。首先......
  • 塑料材料性能
    一、ABS塑料 英文名称:AcrylonitrileButadieneStyrene(丙烯腈-丁二烯-苯乙烯) 比重:1.05克/立方厘米        成型收缩率:0.4-0.7%        成型温度：200-240℃         干燥条件：80-90℃ 2小时物料性能1、综合性能较好,冲击强度......
  • ChatGPT炒股：监控水泥等原材料价格每日变动情况
    水泥、钢铁等大宗商品原材料价格对于很多企业的业绩影响较大，这些产品的价格变动情况需要高频监控，来评估企业未来业绩变动趋势。下面以水泥网的水泥指数，来说明如何通过ChatGPT来编写Python程序，实现监控价格变动的功能。打开水泥网，右键检查：查看水泥指数的对应源代码。找到后，就可以通......

赞助商

阅读排行

网站主页关于我们联系我们网站地图

本网站内容转载自其他媒体，侵权联系[admin##ips99.com]。

Copyright © 2020-2023 IPS99 版权所有 IPS99