【Abaqus】材料行为的非均匀空间分布
设想一种情况:在有限元分析中，一个区域或者整个网格中，每个单元的材料行为都是单独的。这时在ABAQUS中应该如何设置？
两种办法：
- 给每个单元创建一个集合，然后一一赋予SECTION.
- 使用*Distribution关键字，实现空间分布的材料行为，再将SECTION属性赋予给单元。这一种方法好处是，减少后处理的压力，没有那么多的SET和SECTION
两种方法对比：
- 第一种方法可以在CAE界面中手动设置，也可以用Python脚本自动设置，但是一旦单元数量过多，就会导致CAE界面卡顿，且Python速度也不是很快。电脑性能不行的话，后处理压力比较大，因为有很多的Set,Section.
- 第二种方法，不支持CAE，只能修改INP文件实现。优点是后处理软件压力小点，没有那么多set
下面以下面的四个实体单元为例子，分别使用两种方法，直接编写INP文件，命令行中提交计算，进行静态分析。
方法一
*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
结果对比
结果一致
资料
2.8.1 Distribution definition
*ELASTIC
*SOLID SECTION
*SOLID SECTION
*DISTRIBUTION TABLE
1.2.1 Input syntax rules
【Abaqus】材料行为的非均匀空间分布的更多相关文章
- 非均匀B样条拟合MATLAB程序
直接上代码,多的不再说了. %------------------非均匀B样条拟合MATLAB程序----------------- clear k=; x=load('data.txt'); [n, ...
- C语言数组篇（三）字符空间 和 非字符空间
一维数组和字符串 首先是字符数组(区别字符串) ] = {'a','b','c'}; //这只是单纯的字符数组,不是字符串 字符串最重要的标志就是结尾有一个'\0' ...
- 非memory空间有地址分配
对于非memory空间有地址分配,是由于有寄存器配置,比如AHB.APB.一些外设.
- OpenFOAM——设置非均匀边界方法总结
在使用OpenFOAM求解的时候我们经常需要设置非均匀的边界,比如我们在计算层流的时候,很多时候需要入口为充分发展的入口边界,下面我们就以入口处为充分发展的层流速度分布为总结OpenFOAM当中设定不 ...
- OpenFOAM——设置自定义非均匀场区域
在使用OpenFOAM进行计算的时候,我们需要对计算域设置非均匀场,比如最典型的溃坝算例,在开始计算以前,我们需要首先设定某一区域的水的体积分数为1,就是下面这样的: 有可能我们在计算传热问题的时候, ...
- 非均匀B样条离散点的加密与平滑
非均匀B样条离散点的加密与平滑 离散点的预处理是点云网格化很关键的一步,主要就是离散点的平滑.孔洞修补:本文是基于非均匀B样条基函数进行离散点云的加密和平滑的,一下为初步实现结果. 算法步骤: 1.数 ...
- 全国土壤阳离子交换量CEC空间分布数据
数据下载链接:百度云下载链接 土壤阳离子交换量,简称CEC,是指土壤胶体所能吸附各种阳离子的总量.土壤阳离子交换量 cation exchange capacity 即CEC 是指土壤胶体所能吸附各 ...
- 【NLP 系列】Bert 词向量的空间分布
作者:京东零售 彭馨 1. 背景 我们知道Bert 预训练模型针对分词.ner.文本分类等下游任务取得了很好的效果,但在语义相似度任务上,表现相较于 Word2Vec.Glove 等并没有明显的提升. ...
- Linux进程空间分布 & 进程控制块 PCB
p.p1 { margin: 0.0px 0.0px 0.0px 0.0px; font: 14.0px Verdana; color: #555555 } span.s1 { } Linux使用两级 ...
- Linux进程空间分布 & 上下文
Linux使用两级保护机制:0级供内核使用,3级供用户程序使用.从图中可以看出,每个进程有各自的私有用户空间(0~3G),这个空间对系统中的其他进程是不可见的.最高的1GB字节虚拟内核空间则为所有进程 ...
随机推荐
- 【PHP】读取本地文件夹中所有图片并显示
<? //获取文件夹下的所有文件 $dir_str = ''; $imgType = array('gif','png','jpg','jpeg','bmp'); $handle = opend ...
- 全网最详细的 tcpdump 使用指南
https://www.cnblogs.com/wongbingming/p/13212306.html 今天要给大家介绍的一个 Unix 下的一个 网络数据采集分析工具,也就是我们常说的抓包工具. ...
- Qt/C++地图高级绘图/指定唯一标识添加删除修改/动态显示和隐藏/支持天地图高德地图百度地图
一.前言说明 已经有了最基础的接口用来添加覆盖物,而且还有通过进入覆盖物模式动态添加覆盖物的功能,为什么还要来个高级绘图?因为又有新的需求,给钱就搞,一点底线都没有.无论哪个地图厂家,提供的接口都是没 ...
- [转]WorldWind开发中WorldWindowGLCanvas .setPreferredSize()函数找不到
值高温假期,无意翻到了csdn中三维GIS开发的专栏,讲的是worldwind Java三维GIS系统开发的东西,十分感兴趣.恰巧要求的环境已经存在,直接耍起来.将最新的Worldwind和JOGL下 ...
- 陌陌技术分享：陌陌IM在后端KV缓存架构上的技术实践
本文由冀浩东分享,原题"单核QPS近6000S,陌陌基于OceanBase的持久化缓存探索与实践",为了阅读便利,本文进行了排版和内容优化等. 1.引言 挚文集团于 2011 年 ...
- 跟着源码学IM(十)：基于Netty，搭建高性能IM集群（含技术思路+源码）
本文原题"搭建高性能的IM系统",作者"刘莅",内容有修订和改动.为了尊重原创,如需转载,请联系作者获得授权. 1.引言 相信很多朋友对微信.QQ等聊天软件的实 ...
- AI实战营第二期 | 环境配置及安装
AI实战营第二期 | 环境配置及安装 强烈推荐使用 mim 来管理 OpenMMLab repo OpenMMLab repo 不要 mim 和 pip 混用 OpenMMLab 1.0 和 Open ...
- Python学习（一）——配套《PyTorch深度学习实战》
记录一下Python学习过程中的一些问题: 1. 在JupyterLab中查询当前文件的地址 import os print(os.getcwd()) #查询该文件的地址 2. 新建cell 在 Ju ...
- 今天记录一下小程序使用微信客服api，而不是小程序客服
小程序客服缺少很多东西,并且只能使用button的开放能力,所以尝试使用一下微信客服,自己开发客服又比较麻烦,秉着能免费绝不花钱的想法,接下来就直接写代码,也就是api,记录下来方便使用 wx.ope ...
- SqlServer数据插入性能小记
对于SQL性能,已经很久没关注了.由于近期项目既没有几千万的海量数据也没有过多的性能要求,所以在性能上也就没下太多功夫.然而,前不久和朋友闲谈间话题一转就说到程序上了,他说他用Delphi做了个数据导 ...