文档介绍：万方数据
{o{)+)+(1))+(3)≤,ANSYS因子中的应用1科学技术与工程周伟江王生楠分析问题。由于增加了一个时间变量,不仅在数学处理上困难得多,在物理上也复杂得多。由于对于其是表面裂纹前缘与自由表面相交的角点,目前尚无精确的解析解。而三维问题的数值分析,由于裂使之研究受到一定的限制。本文主要介绍有限元软ANSYS{u}{}{i}{F}Newmark衶五,,馨铩{}(1){0)={F())c{){)(4)AICAI三维裂纹问题在动态断裂力学中很少受到注意,因为运动载荷下的三维裂纹提出了一个困难的三维问题应力强度因子沿裂纹前缘线是变换的,尤纹尖端奇异性的复杂性,以及计算工作量大等因素,对所得结果的正确性做了一定的论证。运动方程与逐步积分法将所研究的物体用有限个单元进行离散,其运o【縶五【縶五【縶(1)[][c][K]移、速度、加速度和等效结点力矢量。在时刻出由运动方程得(3)f{)={)+f(1d){)+Ata{+)式衋与口为控制算法精度和稳定性的两个参数。t=0{}再根据式皇,求出下一时刻缸的瑊,礁讅,由此即可得到所有时间离散点上的位移、速度与加速度,进而可以求得各个时刻的应力、应变和应力强度因子。C裂纹顶端渐近位移场与应力强度因子对于线性弹性均匀介质,当载荷随时间变化,稳定裂纹顶端的渐近位移场与静态情况完全类似。对于平面臀侍狻】7920075摘要三维裂纹在动态断裂力学中由于其数学和物理上的复杂性,求解其动态应力强度因子受到一定的限制。文中主要ANSYSANSYA03461B{iI+)+C1t+}+[K1t"I+)={())(ltT710072)年“’.
万方数据
标准分享网载
=姒牡半∽印号(1+y)in芒,平面应力矗气一.,ANSYS3—,平面应变,通过裂纹面上各点处垂直于裂纹面的位移,在平面应变情况下,由式可得r=0K(t)(7)()由于裂纹顶端区域应变场具有庀阶的奇异性。为此在裂纹顶端采用四分之一中点奇异等参元‘。r的奇异性。其他部位采用标准八节点等参元。J模态叠加法。直接积分法又可以分为完全法和缩减法。这里主要讨论直接积分法中的完全法,它是三ANSYSNewmark全分析法步骤计算所有时间离散点上的位移、应力等信息的基础上,再利用软件求解静态应力强度因最后,将所有取出的时刻和对应的应力强度因子的和时间的关系图,即动态应力强度因子。考虑到对称结构,建立八分之一立方体的中心裂纹模型。有限立方的八分之一裂纹模型的尺寸脉冲扰动ⅰ跋煜拢缤所示,其中选择单元偷ピ猄,在单元属性中将纳柚梦狵。。(6)(7)由上述的有限元法求得。单元形式的选取本文求解三维动态应力强度的思想是在按照完子的方法,计算出某个时间点上的应力强度因子。结果导人作图软件,我们就可以得到应力强度因子为了模型的普遍性,采用断裂力学中常见的有1maawI=05lw2=1hwl=2P=1kgm3E=106,泊松比3=100最后进入通用后处理器求解得各个时间点和相应的应力强度因子导人作图软件,其中纵坐标为将9ANSYS1221=。7
万方数据
p足。耐/,其中%是相应载荷下的静态rii45动态应力强度因子转化为无量纲的系数,即CthC=本文的求解思想是基于完全分析法。理论上是以逐步积分法为指导,在