第一篇:Abaqus裂纹模拟心得(Contour Integral不是XFEM)
Abaqus裂纹模拟心得(Contour Integral不是XFEM)
最近由于项目需要,做了一些裂纹相关的模拟,在此把一些心得体会贴到论坛上与大家分享,如有不当之处,欢迎大家指正!
本帖主要侧重于介绍裂纹定义过程中各个选项的意义,具体的操作过程论坛里已经有高手做了很好的教程,至于断裂力学理论推荐大家看一下沈成康写的《断裂力学》一书。裂纹的定义和输出需要用到interaction模块和step模块:
一、Interaction模块
1.1 预制裂纹(步骤:菜单/special/crack/assign seam)
注意:并不是作裂纹分析都要定义seam,如果你的裂纹不是一条缝,而是一个缺口,则不需要assign seam,直接走下一步(定义裂纹)就行。
1.2 创建裂纹(步骤:菜单/special/crack/create,type:contour integral)
—crack front:crack front是用来定义第一围线积分的区域,2D下我们可以选择包围裂尖点的面,3D则选择包围裂尖线的面;另外还有一种定义crack front的方法,就是直接选择裂尖点(2D)或裂尖线3D),用这个方法定义crack front不需要再定义下一步的crack tip/line,比较简便,两种方法算出的结果没有明显的差别,其实只是影响积分路线的问题,但是J积分值是路径无关的,看个人喜好吧
—crack tip/line:这个比较好理解就是裂尖点(2D)或线(3D),如果我们在上一步中用方法二定义crack front,这一步就直接跳过了
—crack extension direction(定义裂纹扩展方向):这里定义的其实是一个虚拟的裂纹扩展方向,定义了这个参考方向后,我们才能通过输出的角度判断裂纹扩展方向,可以通过两种方法:
o q vector:输入一个方向,用来作为计算裂纹的扩展方向的参考方向;
o normal to crack plane:crack plane表示裂纹的对称面(当裂纹在一个平面内时,可能需要分开定义多个裂纹),这种方法下我们只需定义裂纹面的法线方向,通过(t表示裂纹尖端的切线), 会在每个节点得出一个q方向(如下图);
o 注意:q的方向对输出的应力强度因子,J积分等都会有影响,一般情况下,q最好在裂纹平面内,且垂直于裂尖线的切线,否则算出的应力强度因子,J积分值等等在不同围线积分中会差别较大。
二、step模块
定义好了裂纹相关参数后,我们需要返回step模块定义输出变量: 步骤:菜单/output/history output requests/create,domain:crack,可以输出的值包括:J-integral,Ct-integral,stress intensity factor,T-stress
—J-integral :用于应变率无关材料的准静态分析过程,包括线弹性,非线性弹性,弹塑性材料(单调加载工况)的静态分析。J-integral的优点是和积分路径无关,从而可以避开尖端塑性区的影响。
—Ct-integral:用于蠕变分析(一般较少用到)—应力强度因子:
o 只能用于分析线弹性材料,表示裂纹尖端的应力场强度;
o 有三个应力强度因子K1,K2,K3,分别对应于张开型,滑开型和撕开型裂纹的应力强度因子
o 在输出应力强度因子时也会输出一个J-integral值,因为算法不同,这个值和直接输出的J-integral会略有差异;
o 方向判断准则:Maximum tangential stress(在dat文件中输出的MTS值就是通过这个准则算出的裂纹扩展方向),Maximum energy release rate(dat中用MERR表示),K2=0(dat中的K20)
—T-stress(表示裂纹尖端平行于裂纹面方向的应力)
第二篇:裂纹模拟心得
Abaqus裂纹模拟心得(Contour Integral不是XFEM)
最近由于项目需要,做了一些裂纹相关的模拟,在此把一些心得体会贴到论坛上与大家分享,如有不当之处,欢迎大家指正!
本帖主要侧重于介绍裂纹定义过程中各个选项的意义,具体的操作过程论坛里已经有高手做了很好的教程,至于断裂力学理论推荐大家看一下沈成康写的《断裂力学》一书。裂纹的定义和输出需要用到interaction模块和step模块:
一、Interaction模块
1.1预制裂纹(步骤:菜单/special/crack/assign seam)
注意:并不是作裂纹分析都要定义seam,如果你的裂纹不是一条缝,而是一个缺口,则不需要assign seam,直接走下一步(定义裂纹)就行。
1.2创建裂纹(步骤:菜单/special/crack/create,type:contour integral)
—crack front:crack front是用来定义第一围线积分的区域,2D下我们可以选择包围裂尖点的面,3D则选择包围裂尖线的面;另外还有一种定义crack front的方法,就是直接选择裂尖点(2D)或裂尖线3D),用这个方法定义crack front不需要再定义下一步的crack tip/line,比较简便,两种方法算出的结果没有明显的差别,其实只是影响积分路线的问题,但是J积分值是路径无关的,看个人喜好吧
—crack tip/line:这个比较好理解就是裂尖点(2D)或线(3D),如果我们在上一步中用方法二定义crack front,这一步就直接跳过了
—crack extension direction(定义裂纹扩展方向):这里定义的其实是一个虚拟的裂纹扩展方向,定义了这个参考方向后,我们才能通过输出的角度判断裂纹扩展方向,可以通过两种方法:
o q vector:输入一个方向,用来作为计算裂纹的扩展方向的参考方向;
o normal to crack plane:crack plane表示裂纹的对称面(当裂纹在一个平面内时,可能需要分开定义多个裂纹),这种方法下我们只需定义裂纹面的法线方向,通过(t表示裂纹尖端的切线),会在每个节点得出一个q方向(如下图);
o注意:q的方向对输出的应力强度因子,J积分等都会有影响,一般情况下,q最好在裂纹平面内,且垂直于裂尖线的切线,否则算出的应力强度因子,J积分值等等在不同围线积分中会差别较大。
二、step模块
定义好了裂纹相关参数后,我们需要返回step模块定义输出变量:
步骤:菜单/output/history output requests/create,domain:crack,可以输出的值包括:J-integral,Ct-integral,stress intensity factor,T-stress
—J-integral:用于应变率无关材料的准静态分析过程,包括线弹性,非线性弹性,弹塑性材料(单调加载工况)的静态分析。J-integral的优点是和积分路径无关,从而可以避开尖端塑性区的影响。
—Ct-integral:用于蠕变分析(一般较少用到)—应力强度因子:
o 只能用于分析线弹性材料,表示裂纹尖端的应力场强度;
o 有三个应力强度因子K1,K2,K3,分别对应于张开型,滑开型和撕开型裂纹的应力强度因子
o 在输出应力强度因子时也会输出一个J-integral值,因为算法不同,这个值和直接输出的J-integral会略有差异;
o 方向判断准则:Maximum tangential stress(在dat文件中输出的MTS值就是通过这个准则算出的裂纹扩展方向),Maximum energy release rate(dat中用MERR表示),K2=0(dat中的K20)
—T-stress(表示裂纹尖端平行于裂纹面方向的应力)
xfem表面上看属于断裂力学的裂纹扩展问题,实际上它却是用材料弱化的机制研究裂纹的扩展,这点首先放着不提,我来说说断裂与损伤。断裂力学早于损伤力学,格里菲斯最早研究的,断裂力学是为了研究纯在初始裂纹时材料的反映,由于会有应力集中现象出现,应力表达式具有奇异性与裂纹尖端离某一点半径的平方根倒数有关。为了解决这一问题才有了裂纹强度因子,但是无法使用这种比较简单的假设分析整个断裂的过程。以后断裂力学的发展虽然得到了扩充,比如引入了塑性区域和paris模型研究疲劳,但是对于整个断裂尖端材料的性质还是进行了不切合实际的简化,比如理想弹塑性,线弹性等(但是对于金属比较合理)。损伤力学晚于断裂力学卡切诺夫和他的合作者(名字忘了)研究金属蠕变的时候进入了这一研究方法,他可以充分反映材料的弱化,同时从大的范围上来看,损伤力学是涵盖在连续介质力学下面的一个力学直系里面,这一点与断裂力学截然不同,人们无法将断裂力学合理的归纳为连续介质的范畴中,而往往单独叫断裂力学力学。损伤力学出现以后大大抢了断裂力学的风头,因为断裂力学更多偏重静态问题(这种静态是指无法确切的研究整个断裂过程),而损伤力学可以做到,只要定义损伤形式和演化规律,一条或者几条损伤代就可能出现并实现材料或者结构性能的评估。但是损伤力学所能反映的特有性质:材料弱化又成为这种力学体系的一个致命弱点,网格敏感称为损伤力学无法摆脱和倍受攻击的“伤口”,为了完善这一问题,aifantis和他的合作者提出了划时代的非局部化模型理论,而deborst,peerlings,geerling在公元2000年以前将这种非局部化模型体系完全的建立起来,人们发现了使用损伤力学可以真正的实现材料断裂问题的解决方案,随着计算机和有限元力学的发展,现在利用损伤力学的体系计算断裂力学问题已经可以做到了。所以在abaqus软件中或有那么多的损伤力学模型,因为每一种对应一种损伤机理,可以很好的反映材料的弱化以致断裂的模拟(比如使用单元删除技术)。
但是断裂力学与损伤力学两者的结合问题却一直没有得到很好的解决,究其原因在于断裂力学存在裂纹尖端的应力奇异现象,而用损伤力学的方法研究裂纹尖端的时候,会使得裂纹尖端既要保持应力奇异又要保持材料的弱化,这样的耦合问题很难得到合理且唯一的解(这一部分可以参考余寿文老师和冯西桥老师的那本损伤力学)。但是有限元方式使得这种耦合机制得到改善,其实Xfem本身就是这两种力学机制相耦合的产物,xfem中可以预置crack来模拟裂纹尖端的应力奇异,使用材料弱化定义方法实现损伤机理,所以说xfem具有划时代意义也不为过。
ls所说的损伤力学可以实现裂纹发展,但是从机理上存在问题,第一点损伤代具有宽度,而断裂带没有。第二点损伤存在网格敏感问题,网格越细,损伤带越小,与实际不符。第三点,虽然可以利用非局部化模型消除网格敏感,但极其复杂,一般很少采用。第四点,损伤无法真正实现裂纹尖端的应力奇异效应
第三篇:ABAQUS中的断裂力学及裂纹分析总结
ABAQUS中的断裂力学及裂纹分析总结(转自simwe)
(1)
做裂纹ABAQUS有几种常见方法。最简单的是用debond命令, 定义 *FRACTURE CRITERION, TYPE=XXX, 参数。。** *DEBOND, SLAVE=XXX, MASTER=XXX, time increment=XX 0,1, ……......time,0
要想看到开裂特别注意需要在指定的开裂路径上定义一个*Nset,然后在 *INITIAL CONDITIONS, TYPE=CONTACT中定义 master, slave, 及指定的Nset 这种方法用途其实较为有限。(2)
另一种方法,在interaction模块,special, 定义crack seam, 网格最好细化,用collapse element模拟singularity.这种方法可以计算J积分,应力强度因子等常用的断裂力学参数.裂尖及奇异性定义: 在interaction-special,先定义crack, 定义好裂尖及方向, 然后在singularity选择: midside node parameter: 输入0.25, 然后选Collapsed element side, duplicate nodes,8节点单元对应(1/r)+(1/r^1/2)奇异性。
这里midside node parameter选0.25对应裂尖collapse成1/4节点单元。如果midside nodes不移动到1/4处, 则对应(1/r)奇异性, 适合perfect plasticity的情况.网格划分: 裂尖网格划分有一些技巧需要注意,partition后先处理最外面的正方形,先在对角线和边上 布点,记住要点constraint, 然后选第三个选项do not allow the number of elements to change不准seed变化,密度可以自己调整.最里面靠近圆的正方形可以只在对角线上布点.也可以进一步分割内圆及在圆周上布点.里面裂尖周围的内圆选free mesh, element type选cps6或者cpe6,外面四边形选sweep mesh, element type选cps8或者cpe8, 记住把quad下那个缩减积分的勾去掉。
这种方法的几个值得注意的问题,见不少朋友问过。主要是对断裂力学的理解问题。
1.为什么设置理想弹塑性(epp)分析的时候得到的xx,yy方向或者最大应力值Sxx, Syy会超过材料的屈服强度Sy呢, 这分析结果可能吗?
这是因为在ABAQUS中对应等于材料的屈服强度的是von Mises等效应力Se=Sy,因此在平面应变的条件下,xx方向的应力Sxx=Sy*pi/SRQT(3)>Sy, 而Syy=Sy*(2+pi)/SRQT(3), 大概是3倍的屈服应力。所以得到大于材料的屈服强度的xx及yy方向应力是正常的。
2.为什么设置collapse element的时候对弹性分析在中间就一个点而要把单元边上的中点移到1/4处,但弹塑性分析却要在中间设置一圈点并且保持单元边上的中点位置不变呢? 这个其实不是随便定的,在有限元中分析裂纹时,对弹性分析需要模拟裂尖1/SQRT(r)的奇异性,这样在把单元边上的中点移到1/4处后计算出来的等参单元拉格郎日型函数对应的u field正好包含1/ SQRT(r)项,事实上这一方法在断裂力学的数值模拟发展史上是很巧妙的一个发现,至今仍然被广泛采用。至于理想弹塑性分析需要模拟裂尖1/r的奇异性, 这样大家都知道在把单元边上的点放在到1/2处后计算出来的正常的等参单元拉格郎日型函数对应的u field包含1/ r项, 可以模拟弹塑性分析需要的裂尖1/r的奇异性。所以在看似动手点几下就能实现的分析模式后面有很清楚漂亮的理论作支持。
还有就是比较新的cohesive element单元。我仔细读了ABAQUS cohesive element的理论帮助,个人意见ABAQUS的cohesive element采用的是广泛应用于混凝土的类似fictitious crack的方法。只考虑了Dugdale-Barenblatt energy mechanism。这其中softening law 的影响是非常重要的。但ABAUQS似乎只提供了linear 或者exponential 的softening law,复杂的本构关系还需要另想办法。至於基于Griffith-Irwin energy dissipation mechanism的J-integral值可以在LEFM分析中单独算。(ABAQUS用的是Suo Zhigang 和Hutchinson在1990一篇论文中提出的方法)目前cohesive fracture mechnics已经被应用于各种材料。不过在使用到纳米或者更小数量级的研究中碰到了不少问题,可能需要结合位错和分子动力学 的一些理论。现有的cohesive element单元需要定义damage initiation和evolution的准则,softening准则目前好像只有linear和exponential,但对一般材料也够用了。然后通过设置后处理display group可以看到裂纹扩展情况。裂纹扩展不是ABAQUS的强项,目前比较方便的只能用cohesive element,我做过几个模型效果还可以,但对应的参数需要一定的实验数据支持,否则做出来了也不知道对不对。要注意geometric thickness和constitutive thickness; material stiffness和interfacial stiffness的区别 以及厚度与精度的影响。Cohesive element的核心主要是TS-Law,无论里面的数据如何选取,厚度如何变化,cohesive element的表现取决于TS-Law的定义。具体dava的popo10已经给过详细的解释的讨论,可以搜索他们的帖子。如图是我做的3d cohesive element interfacial crack 的例子。
第四篇:模拟外汇心得
模拟外汇心得
实盘交易与虚盘交易的差别
这学期开设了模拟外汇交易实验课,我感到很兴奋,终于可以体会24小时都在交易的全球外汇市场了。实盘交易是一种现货交易。它要求客户手有多少外汇才能做相应金额的交易,客户在完成交易后即持有另一货币,因此它除了提供利用汇率波动赚取差价的机会外,还可满足客户的兑换需求。虚盘交易也称外汇保证金交易。在这种交易中,客户只用缴纳一定保证金后便可与交易商拟定交易合约,进行通常10-400倍的融资,即理论上可做手中外汇金额的10-400倍的交易。由于可将手中外汇放大上百倍进行交易,因而虚盘客户所承担的风险和可得的赢利都被相应放大。
外汇市场的每日交易量已是天文数宇,高达一点九兆美元,使得外汇市场成为全世界最大、资本流动性最高的金融市场,其他金融市场的大小与交易量和外汇市场比较起来,就显得逊色多了。另外,和中国股市、期市相比,外汇市场没有涨跌停板限制,无论何时都可以进行交易或是停损,更不会出现国内股市所谓的庄家和控盘现象。
首先,实盘交易者只能在一定的方向中获利,例如,如果我的本金是美圆,那么在买入其他货币后,只有在美圆下跌的时候,才能获利。如果美圆继续上升,那么我的头寸就被“套牢”,要么斩仓止损,要么白白错过其他交易机会,做实盘就永远要持有一种币种,也就永远不能休息。而保证金操作可以作到对市场方向的收放自如,在其它货币上升时,我可以卖出美圆,买入其它货币(即为作空美圆)获利,而在其他货币下跌时,我可以买入美圆,卖出其他货币(即为做多美圆),增加了买卖的机会。其次还有一个最大的优势就是保证金的杠杆作用,在保证金交易中一般交易商都为客户交易提供一定的融资比例,这个比例越高,我所需要的本金也越少,获利的机会也越多。做外汇保证金交易的投资者,可以通过境外机构提供的交易平台与国际外汇市场实现无缝对接。风险是可以控制的,可预设限价和停止损点。资金是灵活的,随时可抽取资金。外汇市场是全世界资金流动性最强的市场虚盘:外汇保证金交易是一种比股市投资更优越的金融投资工具。其投资方向是全球外汇市场。由于汇价每日波动幅度不像股市那么大,若以实盘换汇交易投资的话,回报率小。利用保证金交易来投资,投资人可以运用杠杆原理,以小搏大、双向交易、操做灵活,而且由于是全球性的市场,行情被个人或机构控制的机会很小,所以相对于股市来说更加透明和公平。
当然,外汇保证金交易属于高风险的投资,不一定适用于每一个投资者,使用高额的融资比例来进行外汇操作有利也有弊。因此在决定进行外汇投资前务必慎重考虑投资目标,熟练程度以及风险承受度
炒外汇的心得体会
通过这次炒外汇我的体会就是要分两部分去修炼自己。一方面是练习心态、另一方面是适应交易系统也就是适应技术分析,并在实战中不断完善自己。这次炒外汇我最大的感受是两点:一是止损,二是控制自己想赚钱的想法。首先止损就是设置好一定的损失范围,超出这个范围就亏损出局,等待下一次机会的出现,每个人都有自己的生活方式和习惯。炒外汇也是同样的也都有自己的分析方法和技巧,这个说要上冲,那个说肯定要破支撑,各有各的道理,各有各的看法。最终总有一些人是对的,但是人毕竟是人,判断和分析总有对和错的时候,天下没有常胜不败的将军,但是失败了要学会怎么把损失控制在最小的程度才是最关键的根本,因为人的本性是贪婪的,赚的时候高兴。亏的时候心情变坏,而且不承认自己是错的,又不舍得割肉出局,总觉得这一切都是暂时的,只要能抗住这段时期就能见光明,殊不知,自己已经越走越远。外汇行业里有一句“止损永远是对的,错了也对,死抗永远是错的,对了也错,止损大于一切”要学会认输,机会还是有的,市场永远是对的,只要能学会顺势而为的道理,那么离成功的机会就越近,这就是止损的重要性,有更多机会去赢取更多的银子,调整好心态,分析失误的原因,找出根本。那么又增加了一些操作经验,丰富了自己的实战水平,避免下次再犯同样的错误,为下次的胜利打下了更坚实的基础。
在自己炒外汇赚钱的时候,要适时的控制自己。一般人都有一种这样的心态,赚了,想赚的更多,亏了想怎么样才能把亏损祢补回来,这是每个人的本性,不懂得控制自己的情绪和欲望做投机这个行业是很危险的,想赚钱是很正常的要求,那么赚多少才能满足自己?我刚开始没有一个计划。因此经常出现赚少了,后悔;亏了,更后悔。我举个自己的实例来证明下止欲的重要性,我在某一个交易日里做的很顺手,哪天的行情把握的很准,简直是感觉货币就是按照自己的想法来走的,结果一天下来资金翻番,当时心里那个美。得意的自己都摸不着北了,要是此时能控制自己的欲望那么一切都很完美,但是偶感觉既然这么顺手在继续多做几笔吧,结果第一笔被套,那么自己的原则忘记了,不设置止损(因为当时太相信自己的判断能力了)反而加仓摊平成本,又犯了第二个错误重仓操作,因为不相信自己会判断错误,又犯了第三个错误,跟市场作对,结果是把当天的赢利全部吃掉,不算,还亏了自己的本金,当时哪个后悔啊,都怪自己没有控制好自己的欲望最终导致自己丢失自己的原则,把自己经常挂在嘴边的话,都当耳旁风抛到银河系了,还好,最后自己知道自己错了,总结了教训,希望自己不要再犯同样的错误。
通过这次炒外汇,我学到了很多知识,了解了炒外汇的基本流程,对自己在以后现实中炒外汇打下了一个坚实的基础。
最后谢谢老师的指导。
第五篇:ERP模拟心得
ERP沙盘模拟实训心得
最终,我们的先锋公司的所有者权益为60,我们的资金量还比较大。第一次上课,是在学校后山,曹老师组织我们玩了一个有趣的“白菜萝卜蹲”的游戏,ERP沙盘以这样一个很特别的方式与我们见面。在游戏中最深刻的体会就是合作与竞争:团队成员之间的合作,团队与团队间的合作、竞争。这一课,终身难忘!
初期:
ERP课程开始的几周,我们什么都不懂,但不了解问题的重要性,碍于提问。直到实训的开始,我们太多的疑问以致于我们的企业无法正常的运行,这个时候才发现懂得的东西太少,这时候我们便不懂就问了,组员解决不了的问题我们就只有求助于老师了,老师也很耐心,讲解的很清楚,很容易理解。在这实训的六年中,我扮演的是一个财务总监的角色,对于对财务毫无概念的我,一开始存在很多问题,让我难以下手,连财务报表都不会填写,对此感到从未有过的无助,但万能的CEO解决了我头痛的问题,我也逐渐步入正轨。
中期:
由于对财务这一块,我没有作系统的、长远的规划,造成在第二年的时候资金断链,这时候我们都不知道该怎么办才好,只好求助于老师,老师告诉我们只有一个办法,就是把厂房卖掉。被迫无奈把厂房出售,暂时解决企业的燃眉之急。从这以后便不敢对财务这一块玩忽职守,在每一年的初期我都会对企业的资金支出与收入作系统的预测,保证财务不对企业造成负面的影响。
后期:
随着对比赛规则越来越了解,我们公司的生产经营得心应手,我们会在新一年的开始就对财务支出与收入作系统的预算,预防资金用到不必要的地方,对公司财务这一块来讲有利于助公司健康茁壮地成长。
总结:
我学到的不从知识获取方面讲,但从企业的竞争与最基本的运作模式,我想我们是有了很深层次的理解,明白一个企业的发展,需要的是强有力的计划与必不可少的团队力量。总之,在这一学期的课程中,我们学会了很多传统课程中学不到的知识,我们更加主动地去思考问题并解决问题,并且考虑问题更加周全。这样的课程的开设对我们存在很重要的意义。
湖南名族学院
0901物流:王莉