首页 >
论文范文 >
自然科学论文 >
力学论文 >
固体力学论文 >
探讨复合材料层合板在三点弯曲工况下的应力反映
摘要：为研究三点弯曲工况下复合材料层合板的应力问题,基于连续介质损伤力学以及Hashin失效准则,建立复合材料层合板三点弯曲有限元模型。模拟三点弯曲成形过程,讨论成形过程中应力分布及变化规律。结果表明:0°铺层应力值始终高于90°铺层,下表面的应力值小于上表面,但下表面的应力值达到最大值的时间要早于上表面;沿X轴方向,越接近压头位置,上下表面的应力差越小,在板料端部上下表面应力值相等;沿Y轴方向,上下表面的应力值变化幅度较小,上下表面的应力差基本保持不变。
加入收藏
纤维增强复合材料比强度、比刚度高,抗疲劳性能好[1,2,3],目前已在航空航天、汽车、建筑、机械等领域得到应用[4,5]。然而,复合材料在外力作用下容易发生纤维损伤、基体损伤、层间损伤等失效模式[6,7,8],导致复合材料零件的力学性能严重退化,使用寿命降低,限制了纤维增强复合材料的进一步推广应用[9,10]。因此,研究纤维增强复合材料在外力作用下的应力状态,获取失效发生的力学机制,对于提高复合材料的成形性能具有重要意义。
三点弯曲试验在机械、材料等领域被广泛采用,是材料力学性能试验的基本方法之一[11]。纤维增强复合材料在弯曲成形过程中受力情况较为复杂,与压头接触一侧受压,对侧受拉,板料内部则是由受压逐渐过渡为受拉,往往发生复合材料失效模式中的一种或几种[12]。然而,复合材料三点弯曲成形损伤机理非常复杂,理论分析由于做了较多的假设而与实际情况误差较大,试验研究又不可能考虑全部因素对成形过程的影响[13]。因此,对于复合材料三点弯曲的数值模拟得到了广泛重视[14,15,16]。
杨慧[17]建立了一种基于实体单元的损伤计算模型,对复合材料层合板三点弯曲失效过程中的基体失效、纤维断裂和层间分层3种基本损伤失效过程进行了准确预测。施建伟等[18]将双线性粘结域单元层应用于复合材料的分层损伤建模,模拟分析了分层的起始以及演化,得到了不同位置分层损伤起始及损伤扩展情况。王伟等[19]利用黏接元法与虚拟裂纹闭合技术对离位增韧复合材料帽型长桁加筋板三点弯曲进行了数值模拟,认为离位增韧法可以改善层合板的抗分层性能。Zhou等[20]通过扫描电镜研究了玻璃纤维增强塑料在三点弯曲条件下的损伤演化,获取了力学性能参数与纤维铺层角度的关系。Santiuste等[21]在三点低速碰撞工况下对比分析了2种不同失效准则对复合材料梁的失效预测结果。然而,已有研究主要针对复合材料层合板的损伤演化,关于损伤力学机制方面的研究目前还少有报道。
本文利用Abaqus软件平台,采用连续介质损伤力学模型建立复合材料层合板三点弯曲有限元模型,通过Hashin失效准则模拟三点弯曲过程中的拉伸以及压缩失效,分析讨论应力分布情况及应力变化规律。
1、有限元模型的建立
在Abaqus中建立复合材料层合板三点弯曲的几何模型,压头和底座为刚体,层合板为变形体,并将层合板与压头、底座等部件装配,压头半径5mm,底座半径2mm,层合板尺寸为80mm×15mm×2mm。复合材料层合板的铺层信息为[90/0/90/0/90/0/90/0],共8个铺层,具体铺层方向如图1所示,铺层从最底层开始依次向上铺设。文中所用复合材料为玻璃纤维增强环氧树脂,其单向复合材料层合板的性能参数如表1所示。
图1复合材料层合板三维模型铺层示意图
采用动态显示算法进行求解,分析时间10s。压头、底座与层合板之间的接触为面面硬接触,摩擦系数0.3。30mm位移载荷施加在压头上,底座的所有自由度均被约束。层合板的单元类型为SC8R,采用扫略方式划分网格。层合板中间区域的网格进行了加密处理,以便更加精确地分析压头附近的应力分布。有限元模型如图2所示。
2、结果与讨论
2.1层合板纵向应力变化规律
在板料上沿X轴方向(板料长度方向)选取由中心至端部包含承载与支撑位置的8个单元,每个单元之间间隔5mm,层合板上下表面对应的单元号从左至右分别为1185(承载位置),855,525,195,3076,3256,3436(支撑位置),3586。
图2复合材料层合板三点弯曲有限元模型
图3给出了层合板上表面8个单元的应力-时间曲线。从图3可看出,层合板上表面不同位置处应力变化趋势基本一致,其中单元1185的应力在6s左右首先达到最大值,此后略有下降,但仍保持在较高的应力水平,单元855的应力值在7s时达到最大值,而单元525在9s左右达到最大值,这3个单元的最大应力值非常接近,先后发生纤维压缩失效。随着单元离板料中心距离的增大,应力值逐渐降低,此时仅发生纤维和基体压缩损伤,并未达到完全失效。而位于板料远端部的单元3436和单元3586,应力变化很小,并且一直保持接近于0,此时基本不发生层合板的损伤。从以上分析可知,由承载位置的板料中心沿X轴方向至板料端部,各个位置的应力值逐渐降低。
图3层合板上表面沿板料长度方向上的应力分布
图4为层合板下表面8个单元的应力-时间曲线。从图4可看出,层合板下表面不同位置的应力变化趋势波动较大,但总的变化趋势与上表面相同,即距离板料中心越远,应力值越小。有3个位置较为特殊,即单元1185、单元855和单元3256。层合板下表面对应着90°铺层,4s时压头与板料接触后,各单元均开始发生基体拉伸损伤,应力值产生波动,其中对应于压头加载位置的单元1185的应力值最先达到最大,随后急剧下降,此时在该单元位置发生基体拉伸失效,随后应力值又有所上升。单元3256的情况与单元1185比较相似,也发生了基体拉伸失效。而单元855发生了较为严重的损伤,虽然没有完全失效,但应力值产生了一定幅度的降低。层合板下表面的应力值波动较大,原因在于与下表面90°铺层1相邻近的为0°铺层2,铺层2在三点弯曲成形过程中会发生纤维拉伸损伤的扩展,与铺层1的基体拉伸损伤扩展一起,导致了下表面应力较大幅度的波动。
图4层合板下表面沿板料长度方向上的应力分布
图5为加载结束时(10s)层合板顶层和底层8个单元的应力分布情况。从图5可看出,在承载位置处层合板上表面的应力远大于下表面,沿X轴方向随着与承载位置间距离l增加,上下表面的应力差减小,在距承载位置30mm处上下表面的应力值相等,此后距离继续增加,应力差值始终为零。
图5沿板料长度方向上上下表面应力对比
由图5结合图3和图4可以看出,虽然上表面的应力值比下表面更大,但下表面的应力值达到最大值的时间更早。此时对应的情况是下表面为90°铺层(铺层1),承载后在4.5s时便发生了基体拉伸失效,此时应力值达到峰值,随后应力值下降;上表面为0°铺层(铺层8),加载后在6.5s时发生纤维压缩失效,此时的应力值达到最大值,随后略有下降。显然,下表面完全失效的时间更早,即达到应力最大值的时间也更早。
2.2层合板横向应力变化规律
在板料上沿Y轴方向(板料宽度方向)选取由中心至端部的8个单元,每个单元之间间隔1.5mm,层合板上下表面对应的单元号依次为1185,1183,1181,1179,1177,1175,1173,1171。
图6为板料上表面不同单元位置处应力的变化情况。从图6可看出,层合板沿着Y轴方向上各位置的单元应力保持相同的变化趋势,4～6s时,应力急剧增加至最大值,随后有所降低,此时对应着层合板上表面(0°铺层的铺层8)在加载位置出现了纤维压缩失效。由图6还可以看出,由板料中心沿着Y轴方向至顶端,应力值略有增加,但是总体来看各单元应力值之间的差异非常小,这也意味着沿板料宽度方向上应力分布非常均匀,从而导致纤维压缩损伤以及失效状态沿板料宽度方向上也非常接近,即损伤以及失效均是沿着整个宽度方向几乎同时发生,这与前期的研究结论一致。
图6层合板上表面沿板料宽度方向上的应力分布
图7为层合板下表面不同单元位置处应力随加载时间的变化。可以看出,加载位置下表层各位置的应力值及其变化趋势存在较大差异。应力沿着Y方向由中心向端部逐渐增大,在板料顶端达到最大值,这主要是由于板料边缘效应的影响。在所选择的8个单元中,靠近板料中心的4个单元的应力变化趋势为先增大,后减小,再增加;而靠近顶端的4个单元的变化趋势为先增大而后减小。各单元的应力差异较大的主要原因在于基体拉伸失效最先开始于板料中心位置,该位置单元应力增加至极限值后迅速降低,随后失效扩展至板料边缘,此时板料中心位置应力又有所上升,此时临近0°铺层(铺层2)中的纤维开始发挥作用承受载荷;而板料边缘的应力则是随着加载过程的进行,应力值逐渐增大,随着失效扩展至板料边缘,应力值逐渐降低,此时附近的纵向纤维还没有进入承载状态。
图7层合板下表面沿板料宽度方向上的应力分布
图8为板料成形结束时上下表面应力分布状态。从图8可看出,板料上下表面的应力在受载位置基本保持稳定,且上表面的应力值远大于下表面。上表面应力值维持在200MPa左右且变化较小,下表面的应力值虽然相对较小,但各位置的应力值波动略大一些。此时,上表面对应着0°铺层发生纤维压缩失效,而下表面90°铺层发生基体拉伸失效。
图8沿板料宽度方向上上下表面应力对比
根据表1可知,基体的拉伸强度与纤维的压缩强度存在着数量级的差别,三点弯曲成形后的应力分布情况也正是基体和纤维强度差异的真实反映。
图9为压头位置下不同铺层的应力变化曲线。从图9可看出,0°铺层的铺层2的应力最大,其次是铺层8和铺层6,而90°铺层的铺层7、铺层5、铺层3、铺层1的应力依次减小。此外,由图9还可看出,0°铺层的应力始终大于90°铺层,层合板上下表面的应力大于内部。因此,各种失效损伤均先发生于层合板的一个表面,随后向内部扩展,最终到达另一个表面。本文研究的复合材料层合板三点弯曲失效模式包括了纤维的拉伸和压缩失效以及基体的拉伸和压缩失效,0°铺层的主要失效模式为纤维压缩和纤维拉伸失效,90°铺层的主要失效模式为基体压缩和基体拉伸失效。各种失效模式发生的先后顺序依次为基体拉伸,纤维压缩,基体压缩,纤维拉伸。原因在于,0°铺层中的纤维能够更好地强化铺层的承载能力,应力值增至较高水平也不会发生铺层破坏;而90°铺层中的纤维对基体的强化作用较弱,铺层在较低的应力水平下便发生失效。基体拉伸失效发生最早,铺层1在4.5s拉应力仅为17MPa时便发生基体拉伸失效,5s时铺层3在13MPa应力值作用下也发生基体拉伸失效,6.5s时铺层2在应力值为246MPa情况下发生基体拉伸失效,0°铺层的铺层2中纤维对基体的强化效果非常显著,应力达到较高水平;此后纤维压缩失效出现,铺层8和铺层6分别在6.5s和8.5s时完全失效,失效时的应力均为232MPa;紧接着发生基体压缩失效,铺层7和铺层5分别在7s和10s时发生完全失效,此时的应力值分别为68MPa和70MPa;完全纤维拉伸失效在文中并未发生,最接近失效状态的铺层2在成形结束时的损伤达到最大值95%,此时的应力值为530MPa。
图9压头位置处各铺层应力随时间
的变化规律
3、结论
1)对于层合板上表面,由加载位置沿X轴方向向板料端部,应力值逐渐降低,不同单元位置的应力值随加载时间的变化趋势也存在较大差异,而沿Y轴方向上各单元的应力值非常接近,变化趋势基本一致。
2)对于层合板下表面,沿X轴方向距离板料中心越远,应力值越小,而沿Y轴方向则是距离板料中心越远(越接近板料端部)应力值越大。
3)沿Y轴方向,层合板上表层的应力大于下表层,而沿着X轴方向,虽然整体趋势也是上表层应力大于下表层,但是随着远离板料中心,上下表层的应力差逐渐减小,直至两者应力值相等。
4)无论沿着Y轴方向还是X轴方向,虽然下表层的应力值小于上表层,但是下表层的应力值达到最大值的时间要早于上表层。
参考文献：
[2]杨燕宁,孟家光,程燕婷,等.三种纬编针织增强体复合材料的力学性能[J].西安工程大学学报,2018,32(2):127-131.
[3]张照.碳纤维织物增强聚醚醚酮基(CFF/PEEK)航空复合材料的制备及其界面改性[D].上海:东华大学,2017:1-15.
[5]丁露露,成小乐,祁志旭,等.工艺参数对B4C颗粒增强铝基复合材料板材挤压过程的影响[J].纺织高校基础科学学报,2018,31(4):430-437.
[6]许玮,李炳宏,贾妍,等.碳纳米管增强铜基复合材料的拉伸性能[J].西安工程大学学报,2017,31(4):533-536.
[7]张琦,高强,赵升吨.碳纤维复合材料板热冲压成形试验研究[J].机械工程学报,2012,48(18):72-77.
[9]张衡,严飙,龚友坤,等.碳纤维机织物增强热塑性树脂复合材料热冲压叠层模型[J].复合材料学报,2017,34(12):2741-2746.
[11]陈同海,周正亮,张守玉,等.夹层复合材料三点弯曲试验和有限元模拟分析对比[J].工程塑料应用,2017,45(7):111-114.
[12]牛智林.三维正交机织复合材料三点弯曲疲劳行为实验研究与有限元计算[D].上海:东华大学,2012:28-29.
[13]堵同亮,彭雄奇,郭早阳,等.碳纤维编织复合材料冲压成形实验与仿真分析[J].功能材料,2013,44(16):2401-2405.
[17]杨慧.复合材料层合板三点弯曲损伤失效数值研究[J].燃气涡轮试验与研究,2006,19(4):38-42.
[18]施建伟,张彦飞,杜瑞奎,等.复合材料层合板三点弯曲分层损伤有限元模拟[J].工程塑料应用,2015,43(2):60-63.
[19]王伟,陈普会,李念.离位增韧复合材料加筋板三点弯曲试验研究与数值分析[J].复合材料学报,2016,33(11):2500-2509.
张华伟,吴佳璐.三点弯曲工况下复合材料层合板的应力变化[J].纺织高校基础科学学报,2020,33(01):32-37.
基金:国家自然科学基金(51475086);中央高校基本科研业务费专项资金(N172304036).
分享：
固体废物是指在生产建设、日常生活和其他活动中产生的污染环境的固态、半固态废弃物质，包括工业固体废物、危险废物、城市生活垃圾、建筑垃圾和农业固体废物等类型。固体废物量多面广、成分复杂，不当处置将使毒害成分以直接或间接途径进入环境，导致长期性、叠加性、复杂性污染，严重威胁生态环境安全和公众健康。
2023-08-28
钢渣是炼钢工程中产生的一种工业固体废渣，经过化学处理和提纯后的钢渣颗粒是一种良好的胶凝材料。据统计，我国钢渣废料的利用率仅为30%左右，远低于西方发达国家，产生了诸多负面效应。钢渣过多堆积会占用大量土地资源，其内部含有硫、磷等有害物质会污染土壤和地下水源。因此，提高钢渣的资源化利用率不仅使其变废为宝，还对生态环境产生积极作用。
2023-08-25
材料力学实验对材料力学教学至关重要。随着时代发展，目前的材料力学实验教学逐渐无法满足本科实验教学需要，无法调动学生的积极性和主观能动性，难以使学生的动手能力和创新能力得到真正的培养和提升[1]。本文首先分析当前材料力学实验教学中存在的问题，然后提出一些改革措施，希望能够进一步完善实验教学，提升学生的整体实力。
2020-09-09
已有研究多关注于有限元软件ANSYS，其建模计算过程相对复杂，本文以材料力学课程梁的弯曲内力图为例，引入有限元分析软件ADINA，该软件建模过程更加简单，使用方便。将之用于梁的剪力图和弯矩图的绘制教学，可快捷、形象地展示梁的弯曲变形动态过程，提高课堂效率。
2020-09-09
材料的力学性能实验是一门材料学科的基础实验课程。通过材料的力学性能实验能够促进理论知识理解，在教学中，不仅要教育学生学习、掌握其基本概念和基本理论，而且还要培养学生灵活运用知识的能力和创新的能力[3]；同时，需要了解和掌握有关材料的力学性能的基本理论知识，能对不同失效产生的原因进行理论分析，从而针对的采取改进手段[4]。
2020-07-14
为了构建以工程为主的特色教育体系，需要融合主体学科各专业特色，重新梳理知识体系，以问题为导向，注重培养学生的核心能力。以“塑性加工力学”课程为例，基于冶金特色的工程教育教学体系，教改工作需要从课程属性与课程目标定位、教学内容设置与安排、教学模式与手段、学生实践能力培养、考核方式等方面进行探索。
2020-06-13
本文依据弹性力学的特点，提出了可视化和参与式教学的相结合的教学方法，以孔口应力集中为案例进行了相关教学探讨，利用有限元软件对多组开孔构件进行有限元分析，采用图形清晰反映外部荷载作用下孔口应力变化规律和特点。同时，让学生参与模型试验测试，更直观的反映孔口应力集中的特点，易于学生接受和理解。
2020-06-13
实验教学是材料力学课程教学中的一个重要环节，但随着时代的不断发展，现行教学方法、手段、模式等已不能满足本科实验教学需要，学生的主观能动性得不到发挥，不利于学生动手能力、创新能力和探索精神的培养，文中根据当前材料力学实验教学中存在的问题进行了分析，并提出了相应的解决对策，旨在进一步完善实验教学，从而促进学生整体实力的提升。
2020-06-08
材料力学课程是机械类及近机类专业的一门重要专业基础课程，该课程开设的实验项目主要有拉伸实验（铸铁和低碳钢）、压缩实验（铸铁和低碳钢）、扭转实验和纯弯曲梁的正应力实验等4个实验，实验学时一般是4～6个学时，主要以验证性实验为主，一般采用在理论课程里面包含实验课时，或者单独与其它课程一起成为单独的一门实验课程，我校采用后一种方法。
2020-06-08
《材料力学》主要是研究杆件的强度、刚度和稳定性问题。材料力学是一门重要的专业基础课，包括机械、土木、水利和交通等专业的学生都会学习到这门课程。由于材料力学本身内容繁杂，概念抽象，对数学工具的应用要求较高，再加上力学课程的课时安排在很多高校中并不充裕，在实际的教学实践中发现相当多的学生对这门课的掌握并不理想。
2020-06-08
人气：6743
人气：5448
人气：4006
人气：3681
人气：2633
我要评论
期刊名称：固体力学学报
期刊人气：1149
主管单位：中国科学技术协会
主办单位：中国力学学会
出版地方：湖北
专业分类：科学
国际刊号： 0254-7805
国内刊号：42-1250/O3
邮发代号：38-44
创刊时间：1975年
发行周期：双月刊
期刊开本：16开
见刊时间：一年半以上
影响因子：0.963
影响因子：1.730
影响因子：0.542
影响因子：0.744
影响因子：1.028
学术期刊咨询、科普杂志阅读、专利申请方案、课题研究设计、专著挂名咨询、论文下载、文献查询、原创文章检测、SCI实验设计方案等多项服务。
服务热线
您的论文已提交，我们会尽快联系您，请耐心等待！
你的密码已发送到您的邮箱，请查看！
2020-06-08
272
上传者：管理员
