西北工业大学学报
西北工業大學學報
서북공업대학학보
JOURNAL OF NORTHWESTERN POLYTECHNICAL UNIVERSITY
2015年
2期
197-203
,共7页
束一秀%李亚智%姜薇%贾雨轩
束一秀%李亞智%薑薇%賈雨軒
속일수%리아지%강미%가우헌
ABAQUS%裂尖%疲劳裂纹扩展%有限元%流程图%矩阵代数%刚度矩阵%应力强度因子%交互积分%水平集法%扩展有限元法
ABAQUS%裂尖%疲勞裂紋擴展%有限元%流程圖%矩陣代數%剛度矩陣%應力彊度因子%交互積分%水平集法%擴展有限元法
ABAQUS%렬첨%피로렬문확전%유한원%류정도%구진대수%강도구진%응력강도인자%교호적분%수평집법%확전유한원법
ABAQUS,crack tips,fatigue crack propagation,finite element method,flowcharting,matrix algebra,stiffness matrix,stress intensity factors%interaction integral,level set method,X-FEM ( Extended Fi-nite Element Method)
研究基于扩展有限元的多裂纹扩展分析方法以及ABAQUS环境下的程序实现。在位移函数中增加扩充项以描述裂纹周围的不连续位移场。初始裂纹以点集或方程的形式给出,并使用水平集函数将多裂纹信息离散到单元节点上,水平集函数还用来追踪裂纹扩展路径。使用局部水平集更新方法减小了计算规模,改进了裂尖单元的判断准则,借助商业软件Tecplot软件实现了裂纹扩展的动态显示功能。使用交互积分法计算混合模式下的应力强度因子,用最大周向拉应力准则判断裂纹扩展方向。探讨了网格密度和积分域尺寸对方法的影响,数值算例表明扩展有限元方法能够模拟任意形状的裂纹,并且能够反映多裂纹扩展的规律。
研究基于擴展有限元的多裂紋擴展分析方法以及ABAQUS環境下的程序實現。在位移函數中增加擴充項以描述裂紋週圍的不連續位移場。初始裂紋以點集或方程的形式給齣,併使用水平集函數將多裂紋信息離散到單元節點上,水平集函數還用來追蹤裂紋擴展路徑。使用跼部水平集更新方法減小瞭計算規模,改進瞭裂尖單元的判斷準則,藉助商業軟件Tecplot軟件實現瞭裂紋擴展的動態顯示功能。使用交互積分法計算混閤模式下的應力彊度因子,用最大週嚮拉應力準則判斷裂紋擴展方嚮。探討瞭網格密度和積分域呎吋對方法的影響,數值算例錶明擴展有限元方法能夠模擬任意形狀的裂紋,併且能夠反映多裂紋擴展的規律。
연구기우확전유한원적다렬문확전분석방법이급ABAQUS배경하적정서실현。재위이함수중증가확충항이묘술렬문주위적불련속위이장。초시렬문이점집혹방정적형식급출,병사용수평집함수장다렬문신식리산도단원절점상,수평집함수환용래추종렬문확전로경。사용국부수평집경신방법감소료계산규모,개진료렬첨단원적판단준칙,차조상업연건Tecplot연건실현료렬문확전적동태현시공능。사용교호적분법계산혼합모식하적응력강도인자,용최대주향랍응력준칙판단렬문확전방향。탐토료망격밀도화적분역척촌대방법적영향,수치산례표명확전유한원방법능구모의임의형상적렬문,병차능구반영다렬문확전적규률。
A numerical method and its implementation in ABAQUS using X?FEM to analyze multiple crack propa?gation was studied. The displacement discontinuity was approximated by appending enriched items to standard field. We give the initial cracks in the form of point?set or curve functions and dispersed them to the element nodes using Level Set Method ( LSM) . LSM was also used to track crack propagation paths. The local level set update technique was used to reduce the computing scale and the Judgmental Approach has been improved. The display of the process of crack propagation was realized using the commercial software Tecplot. Interaction integral technique was used to calculate mixed mode stress intensity factors. The maximal circumferential stress criterion was used to calculate the kinking angles of the propagating cracks. Influence of mesh and integral zone on simulation results were taken into consideration. Numerical examples were presented to demonstrate the benefits of the proposed implementation.