【转载】基于ANSYS APDL的有裂纹平板问题的断裂力学仿真(PLANE183)
对于一般的强度问题,我们总是用应力来度量其强度的。但是对于有裂纹的,高强度的构件,使用应力来度量其强度就是错误的,此时需要使用新的准则来考察其强度问题。
《断裂力学》提供了对于这种问题的强度计算方法,并给出了诸如能量释放率,应力强度因子,J积分等概念来度量含有裂纹构件的强度,以考察一个带有裂纹的构件,在某种外力作用下,它的裂纹是否会进一步扩展;或者如果想要它的裂纹不进一步扩展的话,其裂纹的长度应该是多少,等等。
本篇给出一个最经典的例子,就是一块平板上有一个裂纹,在平板上施加拉力,考虑在该力作用下平板强度的问题。
【问题描述】
一长平板在中间有一水平裂纹,现在板的上下边沿施加均布拉力如下图,要求该裂纹的应力强度因子。
其中材料参数,图中个尺寸的大小以及分布力系的大小如下表。
【问题分析】
1. 该例子来源于ANSYS 15.0 APDL帮助中的一个例子VM256 CINT Command>,帮助中对该例子依次使用PLANE183,SOLID185,SOLID186进行建模,并考察应力强度因子。本文只使用了其中的PLANE183建模部分,并对其中命令的顺序进行了部分整理,并删除了部分笔者以为不必要的程序。
2. 对于2-D裂纹,使用ANSYS所推荐的PLANE183单元。
3. 因为是一个对称问题,只取四分之一建模,并把裂纹尖端点作为坐标原点。
4. 几何建模时对于裂纹用直线表示,而由于裂纹尖端存在着很高的应力梯度,需要对此处仔细划分网格。这里用KSCON指明裂纹尖端,并说明如何在其周围划分网格。
5. 设置对称边界条件,并用CINT定义计算裂纹的相关参数。
6. 后处理中提取出应力强度因子。
7. 本文使用命令流的方式进行求解。
【求解过程】
1. 建模
1.1 创建单元类型
在命令窗口中输入
/PREP7
ET,1,PLANE183,,,2
上述命令确定用PLANE183来建模平面应变问题。PLANE183是ANSYS推荐的建模带裂纹的平面问题的单元。而对于3D中的裂纹建模,ANSYS所推荐的是SOLID186单元。
1.2 输入材料属性
在命令窗口中输入
MP,EX,1,30E6
MP,NUXY,1,0.3
上述命令定义了材料的弹性模量和泊松比。
1.3 创建有裂纹的几何模型
在命令窗口中输入
K,1 ! DEFINE KEYPOINTS AND LINE SEGMENTS
K,2,4
K,3,4,5
K,4,-1,5
K,5,-1
上述命令用于定义四分之一模型的几个关键点。命令输入后结果如下图。
在命令窗口中输入
L,1,2
L,2,3
LESIZE,2,,,4
上述命令绘制了两条直线,并把第二条直线设定为划分4等份。
在命令窗口中输入
L,3,4
LESIZE,3,,,4
上述命令又绘制了一条直线如下图,并设定其网格划分份数为4.
在命令窗口中输入
L,4,5,
LESIZE,4,,,6,.2
L,5,1
ESIZE,,5
上述命令又绘制了两条直线,并设定了其网格划分尺寸的大小。
在命令窗口中输入
KSCON,1,.15,0,8
AL,1,2,3,4,5
上述命令
第1行用于定义应力集中关键点,并确定如何在它周围划分网格。它指定关键点1为应力集中关键点,而围绕该关键点的第一行单元的半径是0.15,在该关键点周围第一圈单元的圆周方向的单元数目是8.裂纹尖端奇异性关键字是0,显然,这是在对裂纹尖端进行网格划分做准备的非常重要的一个命令。
第2行则用于把前面几条直线包围的区域形成一个面。结果如下
1.4 施加边界条件
在命令窗口中输入
DL,1,1,SYMM
DL,4,1,SYMM
上述命令确定边1和4上满足对称边界条件。因为这里取的是四分之一模型。命令执行完毕后在这两条边上有“S”的字符出现,表面是对称(symmetry)
在命令窗口中输入
SFL,3,PRES,-.5641895
上述命令用于在3号边上施加向上的均布载荷。正是这个载荷使得裂纹继续张开。命令运行后结果如下图。
在命令窗口中输入
AMESH,1
用于对面进行网格划分,结果如下
进一步放大裂纹尖端处
可见,围绕裂纹尖端,第一圈单元有8个。
在命令窗口中输入
ALLSEL,ALL
NSEL,S,LOC,X,0
NSEL,R,LOC,Y,0
D,ALL,UX,0
上述命令用于约束住裂纹尖端,使得它不能发生X方向的位移。
1.5 定义裂纹参数
在命令窗口中输入
ALLSEL,ALL
NSEL,S,LOC,X,0
NSEL,R,LOC,Y,0
CM,CRACKTIP,NODE
上述命令用于找到裂纹尖端点,并把它定义成一个组件,名字叫CRACKTIP。定义完毕后,该组件中只有一个节点,就是裂纹尖端处的这个节点。
在命令窗口中输入
ALLSEL,ALL
CINT,NEW,1
CINT,TYPE,SIFS
CINT,CTNC,CRACKTIP
CINT,SYMM,ON
CINT,NCON,6
CINT,NORM,0,2
上述命令首先选择所有对象,然后定义了一系列参数,以便用于后面的裂纹参数的计算。
它首先定义一个新裂纹(NEW),
然后指明要计算的裂纹参数(TPYE)是应力强度因子(SIFS),
接着说明该裂纹的尖端节点(CTNC)是前面定义的组件CRACKTIP,
然后说明该裂纹是对称的(SYMM,ON )
接着指明在进行积分时所用到的等值点的数目是6.
接着指明裂纹的坐标系是全局坐标系,而裂纹面的法线方向是Y轴方向。
在命令窗口中输入
ALLSEL,ALL
FINISH
上述命令用于选择所有对象,然后退出前处理器。
2. 求解
在命令窗口中输入
/SOLU
ANTYPE,STATIC
SOLVE
FINI
上述命令含义简单,就是设置为静力学分析,计算完后就退出求解器。
3. 后处理
在命令窗口中输入
/POST1
*GET,K,CINT,1,CTIP,1,,5,,K1
*STATUS,K
上述命令进入了后处理,并提取裂纹尖端处的应力集中因子,然后显示该值的大小。结果如下图。
可见,该点处的应力集中因子是1.041.
在测出该材料的断裂韧度后,用断裂韧度除以应力集中因子,就得到安全系数,从而可以判断,对于该裂纹构件,在此工作应力的作用下,裂纹是否会进一步扩展。
【转载】基于ANSYS APDL的有裂纹平板问题的断裂力学仿真(PLANE183)的更多相关文章
- 【转载】ANSYS的APDL与C语言混合编程(实例)
原文地址:http://www.cnblogs.com/lyq105/archive/2010/05/04/1727557.html 本文讨论的不是利用C语言为ANSYS写扩展(或者说是用户子程序), ...
- 【转载】ANSYS 动力分析 (9) - 瞬态动力分析 (1)
原文地址:http://htbbzzg.blog.163.com/blog/static/69725206201081663611208/ 第四章 瞬态动力分析 第一节:瞬态动力分析的定义和目的 ...
- 转载 基于Selenium WebDriver的Web应用自动化测试
转载原地址: https://www.ibm.com/developerworks/cn/web/1306_chenlei_webdriver/ 对于 Web 应用,软件测试人员在日常的测试工作中, ...
- [转载] 基于zookeeper、连接池、Failover/LoadBalance等改造Thrift 服务化
转载自http://blog.csdn.net/zhu_tianwei/article/details/44115667 http://blog.csdn.net/column/details/sli ...
- [转载] 基于Redis实现分布式消息队列
转载自http://www.linuxidc.com/Linux/2015-05/117661.htm 1.为什么需要消息队列?当系统中出现“生产“和“消费“的速度或稳定性等因素不一致的时候,就需要消 ...
- [转载] 基于Dubbo的Hessian协议实现远程调用
转载自http://shiyanjun.cn/archives/349.html Dubbo基于Hessian实现了自己Hessian协议,可以直接通过配置的Dubbo内置的其他协议,在服务消费方进行 ...
- [转载] 基于Dubbo框架构建分布式服务
转载自http://shiyanjun.cn/archives/1075.html Dubbo是Alibaba开源的分布式服务框架,我们可以非常容易地通过Dubbo来构建分布式服务,并根据自己实际业务 ...
- 【转载】Ansys中的阻尼
原文地址:http://www.cnblogs.com/ylhome/archive/2009/08/26/1554195.html ANSYS动力学分析中提供了各种的阻尼形式,这些阻尼在分析中是如何 ...
- [转载]基于TFS实践敏捷-项目管理
最近,一直想让团队加入敏捷开发,但TFS2010用下来,并不是令人满意,还好,TFS2012横空出世了.相比TFS2010,TFS2012改进太多了,主要体现在以下方面: Team Web Acces ...
随机推荐
- 解决Cannot find MySQL header files under /usr/include/mysql的错误
按照下面的步骤能成功,亲测.转帖,做笔记 编译php-5.5-6的mysql支持,出现Cannot find MySQL header files under /usr/include/mysql. ...
- caffe中各层的作用:
关于caffe中的solver: cafffe中的sover的方法都有: Stochastic Gradient Descent (type: "SGD"), AdaDelta ( ...
- C++ do{...}while(0)的好处
在开源软件里面经常可以看到这样的写法. #define X(a) do { f1(a); f2(a); } while(0) 1. 主要作用是放在宏定义里面,避免宏带来的语法问题. 比如 #defin ...
- Spring <bean> 参数意义
1.id bean的唯一标识, 2.class 类的完全限定名, 3.parent 父类bean定义的名字. 如果没有任何声明,会使用父bean,但是也可以重写父类.重写父类时,子bean 必须与父b ...
- php短路与 短路或
关于php短路的问题特性,三种写法的区别 $a = 1;$b=0;第一种: $a && $b = 'cccccccc';第二种 $a || $b = 'cccccccc';第三种 if ...
- google protobuf ios开发使用
简介: protobuf 即 google protocol buffer 是一种数据封装格式协议: 比如其他经常用的xml,json等格式:protobuf的优势是效率高,同样的一份数据使用prot ...
- LCT专题练习
[bzoj2049]洞穴勘测 http://www.cnblogs.com/Sdchr/p/6188628.html 小结 (1)LCT可以方便维护树的连通性,但是图的连通性的维护貌似很麻烦. [bz ...
- Mybatis 开启事务@Transactional
- linux修改密码的几种方法
1. 启动电脑 ,进入grub模式. 也就是下面这个模式: 按下e键,进入下面这个画面.... 选第二个(kernel的那个): 然后按下e键之后进入 下面这个版面: 之后敲入 single ...
- 动态SQL基础
动态SQL就是我们可以动态构造SQL代码,然后再执行这个批处理.有两种执行动态SQL的方法,分别是EXEC命令和sp_executesql存储过程.下面是示例代码. USE TSQLFundament ...