abaqus如何设定当Miseabaqus应力云图大于某个值时停止计算?

ABAQUS中应力值显示的太小,附加图片了,怎么调大点,请高手指教_百度知道
ABAQUS中应力值显示的太小,附加图片了,怎么调大点,请高手指教
我有更好的答案
viewport -& viewport annotation options -& legend -& set fonts
采纳率:38%
为您推荐:
其他类似问题
abaqus的相关知识
换一换
回答问题,赢新手礼包
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。查看: 19540|回复: 27
后处理中avg75%的真正含义,及理论基础!
以前后处理的时候看应力不连续就设置avg为100%,但不能真正理解其含义,现在我们可以通过下面的解释,明白这是怎么回事了,希望对这个有疑惑的朋友们,赶紧学习吧!
一. 应力解的性质
应用位移元进行有限元分析时,未知的场函数是位移,即求解的未知量是节点位移。我们知道
应变 strain = B * u, 应力 stress = D * strain = D*B*u
其中u是节点位移向量,B是应变矩阵,D是弹性矩阵
而B是通过单元的插值函数(形函数)N对坐标求导得到的,每求导一次,插值多项式的次数就降低一次,而求导运算使得应变和应力的精度较位移解的精度降低。
因此求解得到的应力与真实解可能有较大误差,其近似性表现在:
1. 单元内部一般不满足平衡方程;
2. 单元之间的界面上应力一般不连续 ;
3. 在力的边界上一般不满足边界条件.
近似应力解的性质是:在精确解上下振荡,并在某些点上近似解等于精确解,即在单元内存在最佳应力点。可以利用近似应力解的这个性质来提高其精度。
最佳应力点的位置:位移近似解是p次多项式,求导的微分算子是m阶,那么应力近似解显然n=p-m次多项式。若应力精确解是n+1次多项式,那么在n+1阶的高斯积分点上应力近似解的精度高于其他位置(有些情况下与精确解相等)。因此单元的最佳应力点就是高斯积分点,这也是为什么应力结果存储于高斯积分点上的原因。
由位移元得到的位移解是在全域连续,应变和应力解在单元内部连续而在单元间一般不连续,即前面所说的应力解第2条近似性:在单元边界上发生突跳。因此会出现的一个现象是:同一个节点由于围绕它的不同单元计算得到的应变和应力通常是不同的。此外,在力的边界条件上也不完全符合。由于工程中我们关心单元边缘和节点上的应力,所以需要对计算得到的应力进行处理。
二. 处理方法
最简单的处理方法就是进行单元平均或节点平均,一般采用算术平均或加权平均的方法。
此外一般采用的方法是应力磨平,包括总体应力磨平、单元应力磨平以及分片应力磨平。
以总体应力磨平为例来介绍应力磨平的基本思路:
构造一个改进的应力解S*,它是在全域连续的。它与有限元求解的近似应力解之间满足加权最小二乘原则(其数学形式是使一个泛函取驻值),通过这个原则进行计算可以解出改进的应力解S*.
这个方法的缺点是计算量太大,以至于比有限元求解位移解的计算量还要大。因此又分别提出了计算量较小的单元应力磨平和分片应力磨平方法。
具体到abaqus的实现,对于存储于积分点的结果(Element-Based Filed Output Results, such as stress, strain,etc.)它默认所采用的是一种条件平均方法(conditional averaging method):
1. 首先将应力张量的分量从高斯积分点上外插到单元的节点上;
2. 默认条件下将在节点上计算不变量(标量),最后再对外插结果进行平均。
‘75%’就是进行平均时的一个参数(叫averaging threshold),理解它需要先介绍另一个节点量即 relative nodal variation 节点相对变化量,在指定区域内进行平均时其定义为
relative nodal variation = (max@node - min@node) / (max over active regions - min over active regions)
当 relative nodal variation & averaging threshold 时,由单元贡献的这个值将参与平均,否则不参与平均。
引用一个例子
&首先某节点是n个单元的公共节点,每个单元在该节点都有一个应力值,比如节点1为周围四个四边形单元的公共节点,那么每个四边形单元在节点1都有一个应力的外推值,假定为(1e6,1.5e6,1.7e6,2e6),最大与最小应力差值为1e6;假定在整个模型显示区域的最大应力与最小应力差值为4e6,那么节点应力相对变化率为1e6/4e6=0.25。在节点1的4个应力值是否会进行平均计算呢?由设定的avg值决定,节点应力相对变化率低于设定avg值的将会被平均。对于节点1,如果设定avg=0.3,那么节点1 输出值为平均应力(1e6+1.5e6+1.7e6+2e6)/4=1.55e6,如果设定avg=0.2,那么节点1处的应力值将不参与平均运算。avg设定为1 ,则意味着所有显示区域所有节点的应力计算都参与平均,所得云图越光滑&
楼主太厉害了,以前看到那个数字都没有多想,学习了;
有个疑问请教下楼主,最后一段话中,如果节点1处的应 ...
应该取的是四个值中最大的一个,因为这样得到的结果会更加保守。
另外自己也亲自测试了一下首先,把Avg设为0%,这样在每个节点处的计算值都不会被平均,如图所示节点被1,2,3,4单元共用,应该有四个应力值,而显示应力计算值为1.112e3。
然后利用查询工具得到该节点上的四个应力,分别为724.557,933.288,550.704,1112.即结果提取的是最大值。
本帖子中包含更多资源
才可以下载或查看,没有帐号?
这个帖子必须支持,lz强大了。
深有收获,感谢
这个帖子,用abaqus的人都应该细细品味。
强烈支持lz。
我调成100%怎么还不连续,求解...
本帖子中包含更多资源
才可以下载或查看,没有帐号?
我调成100%怎么还不连续,求解...
调到100%还不连续,就说明这个情况不是因为差值引起的云图不连续。初步判断,可能是你这里的不是共节点的,也就是这个地方有两个节点,没有paste到一块。如方便的话,可以把文件传上来,更好的交流。
真是有深度的解析,学习到了,多谢分享
理解的这么深入,楼主强大
楼主无敌啊
楼主费心了,细微的问题做详尽的解释。下面的举例说明,更加深了对此概念的解释。利用Python对Abaqus进行后处理结果输出
利用Python对Abaqus进行后处理结果输出
在Abaqus的二次开发过程中,通常需要采用Python脚本语言将Abaqus的计算结果进行输出,然后再进行处理。Python使Abaqus的内核语言,使用较为方便,Abaqus运行Python语言的方式有多种,可以直接命令窗口,也可以读入脚本,还可以采用类似批处理的方式。
本次以一个例子细说Python语言在Abaqus后处理中的应用,模型的计算结果云图如图1所示。
输出所有节点的Mises应力
直接上Python代码:
myodb=openOdb(path='Job-1.odb')
cpFile=open('artlcF1.txt','w')
RF=myodb.steps['Step-1'].frames[1].fieldOutputs['S'].values
for i in range(len(RF))
cpFile.write('.3F\n' % (RF[i].mises))
cpFile.close()
#引入模块,因为需要打开结果文件
#打开结果文件,并复制给变量myodb
#打开一个txt文件
#将输出场赋值给RF
#循环语句,向txt文件逐行写入mises应力
Abaqus的结构层次分的很细,比如结果文件下分如下:
使用过Abaqus的都知道step表示载荷步,frame表示载荷子步,因而在读取Mises应力时需要详细地指定输出哪一步的应力,而应力结果是输出场数据(fieldOutput)的中一种,需要指定是何种应力,程序才知道怎么读取并写入。
由于Abaqus里面涉及的变量特别多,通常很难记清楚那一项下面都有哪些量可以调用,此时比较好的方式是采用print
函数查看,例如查看myodb.steps['Step-1'].frames[1].fieldOutputs下面有哪些变量可以调用,在窗口输入:
myodb.steps['Step-1'].frames[1].fieldOutputs
{'CF': 'FieldOutput object',
'E': 'FieldOutput object', 'RF': 'FieldOutput object', 'S':
'FieldOutput object', 'U': 'FieldOutput object'}
各种不同的结果,包括位移、应力和支反力等等,因此可以知道通过如下的方式读取应力:
myodb.steps['Step-1'].frames[1].fieldOutputs ['S']
此时读取的信息特别多,我们想要的是其中的数值信息,因此可以:
myodb.steps['Step-1'].frames[1].fieldOutputs
['S'].values
通过此句能够读取所有节点的应力数据,输出其中一个:
myodb.steps['Step-1'].frames[1].fieldOutputs
['S'].values[1]
({'baseElementType': 'CPS3',
'conjugateData': None, 'conjugateDataDouble': 'unknown',
array([-4.4, -0.307, 0.0,
1.97], 'f'), 'dataDouble': 'unknown', 'elementLabel':
2, 'face': None, 'instance': 'OdbInstance object',
'integrationPoint': 1, 'inv3': -5.23,
'localCoordSystem': None, 'localCoordSystemDouble': 'unknown',
'magnitude': None, 'maxInPlanePrincipal': -0.041,
'maxPrincipal': 0.0, 'midPrincipal': -0.041,
'minInPlanePrincipal': -5.7, 'minPrincipal':
-5.7, 'mises':
5.11, 'nodeLabel': None, 'outOfPlanePrincipal':
0.0, 'position': INTEGRATION_POINT, 'precision': SINGLE_PRECISION,
'press': 1.46, 'sectionPoint': None, 'tresca':
5.7, 'type': TENSOR_2D_PLANAR})
输出的信息特别多,但是可以看到有mises这一项。
最终得到的txt文件如下:
输出所有节点的位移
RS=myodb.steps['Step-1'].frames[1].fieldOutputs['U'].values
cpFile=open('distance.txt','w')
for i in RS :
cpFile.write('%d %7.4f %7.4f\n' %
(i.nodeLabel,i.data[0],i.data[1]))
cpFile.close()
同样的,输出位移场U,txt写入时第一列是节点编号,第二列和第三列分别是x方向位移和y方向位移。
如上面代码中的红色标记,有nodeLable和data两项,nadeLable使节点编号,data则是位移数据。
在窗口输入代码:
myodb.steps['Step-1'].frames[1].fieldOutputs['U'].values[1].data
分别是两个方向的位移。
再通过data[0]和data[1]读取这两个位移值,最终得到的txt文件如图4。
图4位移输出txt文件
除了读取结果数据,还可以再次基础上做一些运算,同样的输出结果最大值则只需要写一个遍历,比较前后数据大小即可找到最大值,这和其他任何语言都是一样。
文中对输出格式没有解释,可以通过其他书籍和资料了解。
已投稿到:
以上网友发言只代表其个人观点,不代表新浪网的观点或立场。abaqus中应力最大值怎么查看_中华文本库
在ABAQUS 中对应力的部分理解关于 abaqus 中 mises, s11 s22 s33 , s12, tresca pressure, max principal, mid principal,min principal。简单地理解, 在 ABAQUS...
时,应尽可能让其中最大的真实应力和塑性应变大于模型中可能出现的应力和应变值...能用它去核算一些应力内力位移等,abaqus 提供了交互良好的后处理,可以查看各种...
在ABAQUS中实现大体积混凝土徐变温度应力计算_含作者姓名_修改后 071219 - 第XX 卷第X 期 2007 年12月 河海大学学报( 自然科学版) Journal of Hohai ...
实验二 平面问题有限元解的收敛性一、实验目的和要求: 在 ABAQUS 软件中用有限...模型σ x 应力云图: (2)模型σ y 应力云图: (3)底边中点σ x 最大值:...
所 以换算后实际工程总的混凝土在龄期为 3 天的时候的压应力为 σ c = Ecσ s =4MPa,这与图 5 Es 中 ABAQUS 计算结果中的最大压应力的值 3MPa 也...
ABAQUS 需要这些值并相 应地在输入文件中解释这些数据。 大多数实验数据常常是用名义应力和名义应变的值给出的。这时,必须应用 下列公式把塑性材料的数据从名义应力...
在Abaqus 中使用梁单元进行计算在 Abaqus 中使用梁...这可从 mises 应力与 S11 的比较看出来;而为什么...通过放大系数的选择使稳定时间增量最大 化;梁横截面...
基于ABAQUS的热应力分析_机械/仪表_工程科技_专业...完成后打开工作目录下对应的 odb 文件查看分析结果。...节点坐标值的文本文件和对应节点最大 位移值,并以此...
ABAQUS加重力,地应力,初应力,输出图片格式_建筑/土木_工程科技_专业资料。ABAQUS...A:使用 rebar element Q:请问各位高手,如何在 viewer 里将 contour 最大值与...
关于ABAQUS 历史输出变量 在帮助文档 Analysis User’s Manual 中 4.2.2 节 ...(用 *Node Output 关键字表示)可以包括的变量,从中可以看到应力(S)是不能...关于ABAQUS的运用体会及一些经验总结_百度文库
您的浏览器Javascript被禁用,需开启后体验完整功能,
享专业文档下载特权
&赠共享文档下载特权
&10W篇文档免费专享
&每天抽奖多种福利
两大类热门资源免费畅读
续费一年阅读会员,立省24元!
关于ABAQUS的运用体会及一些经验总结
阅读已结束,下载本文需要
想免费下载本文?
定制HR最喜欢的简历
你可能喜欢}

我要回帖

更多关于 abaqus输出应力云图 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信