层板胶合木层间断裂韧度研究

层板胶合木层间断裂韧度研究

ID:35022216

大小:16.20 MB

页数:77页

时间:2019-03-16

上传者:U-24835
层板胶合木层间断裂韧度研究_第1页
层板胶合木层间断裂韧度研究_第2页
层板胶合木层间断裂韧度研究_第3页
层板胶合木层间断裂韧度研究_第4页
层板胶合木层间断裂韧度研究_第5页
资源描述:

《层板胶合木层间断裂韧度研究》由会员上传分享,免费在线阅读,更多相关内容在学术论文-天天文库

犬拔?^大fDALIANUNIVERSITYOFTECHNOLOGY硕_士害位1£文MASTERALDISSERTATION_层板胶合木层间断裂初度研究学科专业结构工程作者姓名王亚敢指导教师徐博答辨日期2015jg_6R 硕士学位论文层板胶合木层间断裂_度研究StudyonInterlaminarFractureTouhnessg-inGluedLaminatedTimber??.作者姓名:王亚勒学科、专业:结构工程学号:21206042指导教师:徐博输完成日期:2015年5月大理工大旁DalianUniversityofTechnology 大连理工大学学位论文独创性声明作者郑重声明:所呈交的学位论文,是本人在导师的指导下进行研究工作所取得的成果。尽我所知,除文中巳经注明引用内容和致谢的地方外,本论文不包含其他个人或集体已经发表的研究成果,也不包含其他已申请学位或其他用途使用过的成果一。与我同工作的同志对本研究所做的贡献均巳在论文中做了明确的说明并表示了谢意。若有不实之处,本人愿意承担相关法律责任。学位论文题目:层板胶合水层间断裂初度研究作者签名:日期:^月1日 大连理工大学硕士学位论文.摘要层板胶合木凭借其良好的力学性能而被广泛应用于现代木结构建筑中。在许多情况下层板胶合木由于出现断裂而发生破坏,现阶段已有的相关研究大多集中在其层内断裂一行为方面,仅有少量关于其层间断裂行为的研宄。木材是种正交各向异性材料,不同年轮角度的层板组合可能会对层板胶合木层间断裂初度产生影响,因此对该问题进行断裂力学分析是必要的一,进而对层板胶合木的生产工艺及性能研究也会提供定的指导。基于层板胶合木层间断裂韧度问题研究的必要性,本文开展了如下的研宄工作:(1)层板胶合木层间断裂朝度试验测定分别对不同年轮层板组合的层板胶合木进行了I型层间断裂的双悬臂梁试验和11型层间断裂的端部切口三点弯曲试验,在得到临界荷载Pic和Piic后基于经典梁理论计算了其层间断裂勒度双悬臂梁试验结果表明,该种试验所得的层板胶合木层间断裂初度G一ic随着组合层板的年轮角度差值的增大有定的增大趋势;端部切口三点弯曲试验结果表明,该种试验所得的层板胶合木层间断裂铜度C??c随着组合层板的年轮角度差值的増大有轻微增大的趋势。对不同年轮层板组合的层板胶合木进行了I型层间断裂的单边切口三点弯曲试验,试验结果表明,该种试验测得的层板胶合木层间断裂临界荷载Pic随着组合层板的年轮角度差值的增大有轻微的减小趋势。(2)层板胶合木层间断裂韧度GG有限元分析ic和iic的分别对不同年轮层板组合的层板胶合木进行了I型层间断裂的双悬臂梁试验和II“型层间断裂的端部切口三点弯曲试验的有限元模拟,釆用ABAQUS程序利用回路积”计算得到J积分值,分法,在线弹性情况下其与断裂軔度Gic和Gnc是等价的有限元计算的结果与试验测定结果吻合良好。对不同年轮层板组合的层板胶合木进行了I型层间断裂的单边切口三点弯曲试验“”的有限元模拟,采用ABAQUS程序利用回路积分法计算得到?/积分值,在线弹性情况下其与断裂朝度G,,ic是等价的结果表明该种试验不同年轮角度差值的层板组合对于层板胶合木层间断裂钥度G,c无显著影响。双悬臂梁试验的层间断裂軔度G,c有限元计算结果与单边切口三点弯曲试验的层间断裂朝度Gic有限元计算结果相比偏小,这可能是由于木材性质的离散性、试件的形式与尺寸差别、裂纹制作误差、提取临界荷载的误差误差等因素造成的。(3)层板胶合木层间断裂钥度iiTic和的有限元分析-I- 层板胶合木层间断裂擁研宄对不同年轮层板组合的层板胶合木进行了I型层间断裂的双悬臂梁试验和II型层间断裂的端部切口三点弯曲试验的有限元模拟,采用应力外推法计算得到层间断裂軔度尤A:。,IC和?C双悬臂梁有限元计算结果表明该种试验所得的层板胶合木层间断裂初度iiTic随着组合层板的年轮角度差值的增大有轻微的增大趋势;端部切口三点弯曲有限元计算结果表明,该种试验所得的层板胶合木层间断裂韧度K?c随着组合层板的年轮角度差值的增大有轻微的增大趋势。对不同年轮层板组合的层板胶合木进行了I型层间断裂的单边切口三点弯曲试验的有限元模拟,,,采用应力外推法计算得到层间断裂軔度尺IC结果表明该种试验中截面方向上不同年轮角度差值的层板组合对于层板胶合木层间断裂临界荷载无显著影响。双悬臂梁试验的层间断裂朝度Kic有限元计算结果与单边切口三点弯曲试验的层间断裂朝度if元计算结果相比偏小,这可能是由于木材性质的离散性、试件的形ic有限式与尺寸差别、裂纹制作误差、提取临界荷载的误差以及采用有限元应力外推法进行计算应力强度因子的误差等因素造成的。关键词:层板胶合木;断裂初度;层间断裂;年轮方向有限元;--II 大连理工大学硕士学位论文n-StudyonIterlaminarFractureTouhnessofGluedLaminatedTimbergAbstractBecauseofitsenvironmentalandeconomicbenefits,glulamiswidelyusedinconstructionsectors.Howeverthefracturebehavioroflulamismorecomlexthanthatof,gpwood.Therearetwotypesoffracturethatcanoccuringlulammembers:intralaminarfractureandinterlaminarfracture.Intralaminarfractureisinessencethefractureinwood,whileinterlaminarfractureoccursattheinterfaceoftwoadacentlaminae.Woodisanoithotroicjpmaterial,thefracturetoughnesspropertiesdependondifferentgrowthringorientations,thustecestttltlncetceofiisnsaryodoheineraminarfracureanalsisofluamcorninghedifferenyggrowthringorientationsandthenrovidessomeuidancetoroductionrocessand,pgppperformancestudyoflulam.Therelativestudieswerecarriedoutasfollows:g(1)ExperimentalinvestigationontheinterlaminarfracturetoughnessoflulamgThedoubDCB-lecantileverbeamtestsforModeIfractureandthreeointbendin()pgend-tcexure3tModIInohedflENFestsforefracturewereerformedtoobtainthe()pinterlaminarfracturetoughnessGicandGncoflulamwithfocusinontheinfluenceofgggrowthrinorientationstheG\candGmcwerecalculatedbelementarbeamtheorafterg,yyya*ccordintothecriticaloadsPand/obtainedintheests.testsresultssglicuctDCBhowedGincreaseatacertaindereeastheincreaseofthedifferencebetweentwoadacenticgjlaminae.3ENFtestsresultsshowedGcincreaseslihtlastheincreaseofthedifferencengybetweentwoadacentlaminae.j-Thethreeo--intbendintestsofasinleedenotchedsecimen3SENBforModeIpgggp()fracturewereerformedwithfocusinontheinfluenceofrowthrinorientations,thetestspgggresultsshowedPdecreaselihtlastheincreaseofthedifferencebetweentwoadacenticsgyjlaminae.(2)DeterminationoftheinterlaminarfracturetoughnessGicandGncofglulaminfiniteelementanalysisThefiniteelementanalysisofDCBtests3ENFtestswereperformedwithfocusingon“"theinfluenceofrowthrinorientationstheABsoftwareusesthecontourinteralg,AQUSggm??whichcouethodtocalculate/interalldbeinterretedasbeineuivalenttotheenerg,pgqgyreleaserateGunderlinearelasticmaterialassumtion.Accordintothecriticalloads尸candpgi尸ic,theGicandGwccalculatedbyfiniteelementmethod(FEM)werefoundbegoodareementwiththetestsresults,gIII 层板胶合木层间断裂郁度研究The-mentsamefiniteelementanalysisaboveionedwerealsoconductedon3SENBteststoobtaintheGicthenumericalresultsshowedthattheeffectofrowthrinorientationsof,ggtwoadacentlaminaeonGcwassliht.jigThenumericalresultsofGicobtainedfromDCBweresmallerthanthatobtainedfrom3SENBtests,whichmightbecausedbythewoodpropertiesdiscreteness,thediferenceof'specimenssizeandformthecrackroducinerrortheerrorofcriticalloadsdetermination,,pg,andsoon.(3)DeterminationoftheinterlaminarfracturetoughnessKicandKncofglulaminfiniteelementanalysisThefiniteelementanalysisofDCBtests3ENFtestswereerformedwithfocusingonptheinfluenceofrowthrinorientationstheKand尺werecalculatedusinthestressgg\cnc,gextrapolation.DCBtestsnumericalresultsshowedKicincreaseatacertaindegreeastheincreaseoftdifferenceetweentwoadacentamnae.3ENFtestsnumericalresultsshowedhebjliA^increaseslihtlastheincreaseofthedifferencebetweentwoadacentlaminae.icgyjTh-esamefiniteelementanalysisabovementionedwerealsoconductedon3SENBteststoobtaintheKic,thenumericalresultsshowedthattheeffectofgrowthringorientationsoftwoadacentlaminaeonwasslight.jThenumericalresultsofKicobtainedfromDCBweresmallerthanthatobtainedfrom3SENBtestswhichmihtbecausedbythewoodropertiesdiscretenessthedifferenceof,,,gpsecimenssizeandformthecrackroducinerrortheerrorofcriticalloadsdeterminationp,pg,,theerrorofthestressextrapolationandsoon.,KeyWords:Glulam;FractureToughness;InterlaminarFracture;GrowthRingOrientations;FiniteElementMethodIV 大连理工大学硕士学位论文目录mmIAbstractIll--1舰1-11概况-1.--12.断裂力学2--11.2.断裂力学的起源与发展22--1.2.能量释放率与脆性断裂的能量准则4-123-..应力强度因子与断裂准则613-6-.木材顺纹断裂--1.3.1木材顺纹断裂行为的研宄92--1.3.木材胶结面断裂行为的研究15--1317.3.层板胶合木断裂行为的研宄-14本文的-17.研究内容--1.4.1层板17胶合木层间断裂勃度Gic和Guc的研究--1.4.2层板胶合木层间断裂朝度和iiCnc的研宄18-2-1层板胶合木层间断裂试验9--2.1木材弹性常数测定19--2.2层板年轮方向角度代表值22-2-.3层板胶合木层间断裂DCB试验27-231DCB-..试验概况27--2.3.2DCB试验结果30--2.4层板胶合木层间断裂3ENF试验32-213EN-.4.F试验概况322423ENF--..试验结果34--2.5层板胶合木层间断裂3SENB试验35--2.5.13SENB试验概况35--2.5.23SENB试验结果372--.6本章小结38--3:.有限元分析40-31DCB试验的有限元分析-42.层板胶合木层间断裂--3.1.1DCB试验能量释放率G计算42-V- 层板胶合木层间断裂韧度研究.--3.1.2DCB试验应力强度因子计算44---32ENF46.层板胶合木层间断裂3试验的有限元分析.--3213E46..NF试验能量释放率G计算:...3--.2.23ii:计算48ENF试验应力强度因子-33SENB-51.3层板胶合木层间断裂试验的有限元分析-331351-..SENB试验能量释放率G计算…B--3.3.23SEN试验应力强度因子ii:计算53--3.4试验结果与有限元结果对比553-55-.4.1DCB试验结果与有限元结果对比-3ENF55-.4.23试验结果与有限元结果对比343DCB3SENB-56-..有限元与有限元结果对比3--.5本章小结57--4结论与展望59--4159.誠--4.2展望60--^^Xm61-66-攻读硕士期间发表学术论文情况--致谢67-68-大连理工大学学位论文版权使用授权书VI 大连理工大学硕士学位论文1绪论1.1概况木材是人类从大自然获得的一种天然材料,被广泛应用于建筑和室内装修等领域,在欧美等国家,大量的住宅、学校、办公楼等建筑主要采用现代技术生产的工程木来建⑴造。与时下被普遍应用的砖混结构、钢筋混凝土结构以及钢结构相比,木结构有很多2[]优点:节能环保、木材可再生、轻质高强等。虽然由于历史原因目前我国林木资源相对匮乏,但是经过国家相关部门多年的经营,林木资源日益丰富,尤其是速生林的种植规模不断壮大,这些都为利用木质建筑材料发展绿色建筑提供了契机。一一结构用木材般分为两类:针叶材和阔叶材。针叶材质地般较软,又被称为软木;3一][而阔叶材质地般较硬,所以又被称为硬木。结构中的承重构件大多采用针叶材。目前已经列入我国《木结构设计规范》的结构用木材,有落叶松、红松、云杉等国产针叶材和桦木和水曲柳等国产阔叶材,此外还有北美花旗松、欧洲云杉和俄罗斯红松等进口木材。随着胶合技术的发展,出现了各种以不同大小天然木材单元体为原料进行重组胶一[4]工程木有层板胶合木合而成的木质材料。目、,般被称为工程木前应用比较广泛的旋切板胶合木、平行木片胶合木和层叠木片胶合木等。20mn ̄45n层板胶合木,是以厚度为iim的板材,沿顺纹方向叠层胶合而成的工程木5[]制品。由于在层板胶合木的制作过程中,人工剔除了木材自然生长带有或者后天因素造成的各种缺陷,因此层板胶合木的抗拉和抗压强度均优于天然木材,并且在材料质量6[,”的均匀化方面也优于天然木材。在西方的发达国家和地区,,木结构的应用已经有很长的历史并且胶合木的应用也十分普遍。由联合国欧洲经济委员会(UNECE)和联合国粮农组织(FAO)共同发布—8[]的《20072008林产品年度市场报告》显示,在北美地区,用于民用建筑或者学校、商店等其他非民用建筑的工程木占其全部产量的60%以上;欧洲地区2007年的工程木3制产品的总产量大约为250万m,这些产品也主要被用于桥梁和体育馆等大型建筑的建造。随着新型木结构用胶粘剂和连接件的出现与发展,各种新建桥梁、车站、商场以及体育馆等公共建筑和商业建筑也大多采用胶合木结构,展现了胶合木结构的经济性、‘实用性以及艺术价值。随着我国经济的快速发展,人民物质生活条件的不断提高,低碳环保的意识也逐渐在人民的心中生根发芽苗壮成长,因此,越来越多的人开始对既经济环保又美观舒适的木结构建筑产生了偏爱。目前在苏州、成都、上海等许多城市,许多--1 ‘--层板胶合木层间断裂朝度研宄.层板胶合木建筑相继建起,如广西柳州开元寺、苏州胥河上即将建成的仿古桥、成都欢9口门[】乐谷入厅、上海欢乐谷过山车结构体系等。在按照传统的材料力学强度理论进行材料强度校核时,木材被视为均质、连续的理想固体,忽略了其内部存在的各种初始裂纹,因此在使用中有可能发生初始裂纹扩展引起的意外断裂行为。天然木材是自带缺陷(裂纹和木节等)的,并且木材在加工和水分循环变化过程中也会产生裂纹,在这些缺陷附近或者裂纹尖端都可能出现应力的奇异’1G]证实了木材内部微小裂紋的起裂与扩展是其失效的主要原因【性。实践。由于传统旳材料力学强度准则并不能有效的预测含有初始裂纹木材的破坏,通常需要采用断裂力学准["]则评估含裂纹构件在低应力下的强度和安全。层板胶合木与传统原木、方木相比最大的优势在于层板胶合木构件的尺寸可以不受。限制,而这种优势得益于其制造工艺过程中的接长和层积手段鉴于这两项技术在整个一,,以获得科学的接长和层积技术参数制造阶段的重要性国内外学者直在研究,进而12[]制造出理想的结构用木材产品。由于层积制作的要求,层板胶合木中存在大量的胶结一一面,而胶结面上的边界条件实际上代表了侧的材料对另侧的变形约束,这种约束可13能会使胶结面端部[】、折点、裂纹等处也产生应力奇异性,引起层板胶合木层间断裂行为的出现,进而导致层板胶合木构件发生分层破坏。一R木材作为种正交各向异性材料,,,其三个弹性主轴(纵向L轴径向轴弦向T)轴方向上弹性常数以及强度差异较大,所以在进。由于木材弦向和径向的胀缩率不同行层板胶合木层板拼合时需要依次改变年轮布置方向。如果年轮方向垂直于板面,可能会使木材湿度发生变化时所引起的翘曲度减小。层板胶合木胶结面两侧很可能出现不同年轮角度的层板组合,这代表了不同的边界条件,也可能会引起层板胶合木层间断裂韧度的变化。因此,本文基于断裂力学理论,釆用试验与有限元模拟相结合的方法,研宄层板胶合木胶结面两侧不同年轮角度层板组合对其层间断裂軔度的影响,本研宄对层板胶合木的生产工艺及性能研究具有一定的理论价值和指导意义。’1.2断裂力学1.2.1断裂力学的起源与发展一断裂力学出现之前,人们般釆用传统的强度准则进行构件或者结构的设计,通过保证构件不产生强度失效从而保证整体结构可以安全工作。传统强度准则即经典强度理15论的基本思想是保证构件的工作应力不超过材料的许用应力[],可表示为:(1.1)-2- 大连理工大学硕士学位论文经典强度理论作为结构设计的基本理论被沿用至今,并且在进行设计时考虑了由于荷载计算、材料特性和分析方法等诸多不确定因素引起的误差,通过釆用量值大于1的安全系数来保证构件和结构的安全但是,在生产实践中不断出现了构件的工作应力低于材料的许用应力情况下的脆断事故,尤其是高强度和超高强度材料构件、低温或腐、。烛环境中的结构爆接结构等部分上述情况的典型脆断事故如下:(1)1885年长岛的铆接水塔在静水压力试验中脆断。“”(2)第二次世界大战期间,美国全岸接自由号货轮连续发生多起脆断事故,到1948年,其中的238艘因损伤严重而完全报废。(3)1944年,美国克里夫兰发生了液化气罐大破裂的灾难性事故。“”(4)1954年,英国的deHaviland公司研制的旅客机彗星I号在印度上空失事。“”(5)1958年,美国北极星导弹的发动机壳体在进行试车时突然发生爆炸。一系列的脆断事故产生的原因,主要是基于经典强度准则的基本设计理论将材料作为无缺陷的均勾连续体进行研究,忽略了材料和构件以及结构内部的缺陷、煌接处存在的爆缝以及疲劳荷载等因素的作用,使得微观裂纹发生扩展造成的低应力状态下的脆断事故。为了保证结构和构件的安全工作,以裂纹体为研宄对象的断裂力学应运而生。i7i8,[1921年和1924年,Grifith%现由于内部细小的裂纹缺陷导致玻璃的实际强度远低于采用分子结构理论所预期的理论强度的现象,并基于能量原理提出了裂纹的失稳扩展条件:当裂纹扩展所释放出来的应变能(裂纹扩展力)大于形成新裂纹所需的表面,裂纹就会失稳扩展能时,这项研宄奠定了现代断裂力学发展的基础。严格地说一,Grifith理论只适用于理想脆性材料,对于金属材料等般的工程材料一i9[2G[]应用误差较大。为了将Griffith理论转化为门工程科学,Irwin^Orowax各自独立的考虑了裂纹尖端区域塑性耗散功的存在,然后对裂纹扩展表面能进行了补充与修G一正rifith理论适用于,使得修正的般工程材料的断裂问题分析,并引入了能量释放2i[]率G,为线弹性断裂力学的快速发展奠定了理论基础。1957年in基于裂纹尖端,Irw附近应力场奇异性的分析,提出了应力强度因子尤的概念,并将能量释放率G和裂纹*一i:尖端应力强度因子i通过等效模量五联系起来,证明了上述两个参数的等效性,进步推动了线弹性断裂力学的发展与成熟。因此一,现代线弹性断裂力学中用来度量材料中裂纹推动力的参数般有两个,即能量释放率G和应力强度因子A:,在线弹性材料I型,但是这两个参数并不是相互独立的和11型断裂类型中存在以下的对应关系:r2jG=2)(1.^E-3- 层板胶合木层向断裂賴3度研究*[22]式中:五为等效弹性模量。对于木材的顺纹断裂则有式(1.3):士?(1.3)提谓女*式中:X和y分别代表木材试件长度方向和高度方向的弹性主轴,fx、分别为两个方向的弹性模量,GW为剪切模量,#x为泊松比。y一3在断裂力学的理论分析中,根据裂纹的受力和位移特点,般把裂纹抽象为种基本类型(参见图1.1)(1)张开型裂纹(I型):外荷载是垂直于裂纹平面的正应力,裂纹面相对位移,也是垂直于裂纹平面的则其裂纹面上表面和下表面沿着y方向的位移分量不连续;(2)滑幵型裂纹(II型):外荷载是平行与裂纹平面的剪力,裂纹面在其自身平面内做垂直于裂纹前缘的剪切滑动,则其裂纹面上表面和下表面沿着X方向的位移分量不连续;(3),撕开型裂纹(m型):外荷载是离面的剪力裂纹面在其自身平面内作平行Z于裂纹前缘的错动,则其裂纹面上表面和下表面沿着方向的位移分量不连续。M■_A"斷(J丨(a)张开型裂纹(b)滑开塑裂纹(C)撕裂型裂纹图1.1裂纹的3种基本类型F.1.1Teebacofighrsitypescrack1.2.2能量释放率与脆性断裂的能量准则能量释放率G是线弹性断裂力学中裂纹扩展推动力的重要参数之一,其代表着裂纹扩展时释放出的能量提供了形成新裂纹表面的表面能和裂纹尖端附近的塑性形变功。由于能■释放率G是通过能量分析所得出的概念,所以具有普遍适用的特点。假定一裂纹体的初始裂纹面积为儿当其裂纹面积扩展了&4时,在裂纹扩展的过程中外荷载做功为裂纹表面能增量为3r,塑性功增量为3A,整个系统的弹性应变能增量为5U一。根据能量第定律,当不考虑能量损失时,整个体统的能量是守恒的,即外荷载所做的功等于整个系统能量的增加,那么:-4- 大连理工大学硕士学位论文5W=5U+ST+SK(1.4)-/:況<5^,记『,为裂纹面积扩展时整个弹性系统所释放的能量则上式看写为-5U=-=SWdU5T+SK(1.5)将裂纹扩展单位面积时所释放的全部能量定义为其能量释放率G’则有:見a(1.6)dAdAdA能量释放率G代表了在裂纹扩展单位面积时向外提供的推动裂纹扩展的能力,因此能量释放率G也被称为裂纹t展力。在具体试验中可以通过测定试件柔度去确定能量释放率G。在外荷载位移A恒定的A==情况下0,0。,在裂纹扩展的过程中,施力点的位移无变化,即i/故外力功增量^^^^(1则由式.6)可得:G一担(1.7)dAdA在线弹性情况下,弹性体的变形能等于外力所做的功,并且加载点位移A与外荷载P成正比,即:=-UFA(1.8)2A=)CAP(1.9()式中,,为裂纹体的柔度它只是裂纹面积乂的函数。=(1.9)可得由£M0和式:===-dAACdP+PdCA0PdC(1.10){A)即CA)dP){)((=再由式(1,.9)和式(1.10)<M0,可得:.8)并考虑式(1以及^=-----^^^^--MP=dUPdh^=CAPPdCA(111)MP+i).i))(2222C(A)2111)17将式(.再代入式(.)可得:G==型(112)dAdA2dA则基于能量释放率的断裂准则可表示为:G>G(1.13)c其中用Gc为裂纹矿展单位面积时所消耗的能量,当能量释放率G超过材料对裂纹扩展抵抗能力Gc时即认为发生断裂。G和Gc的单位相同,其国际单位为N/m;工程上2常用单位为J/m。5 层板胶合木层间断裂韧度研宄一定材料而言,对于其内部裂纹扩展过程所消耗的总能量被认为是材料常数,与外荷载情况、试件形状与尺寸以及裂纹的几何形状与尺寸等因素无关,因此Gc反映了材料抵抗断裂破环的能力,,通常被统称为材料旳断裂軔度需要通过试验确定。对于I型和II型裂纹,能量释放率分别记为和Gn,相应的断裂初度分别记为Gic和Guc。1.23.应力强度因子与断裂准则基于裂纹尖端附近应力场的分析表明,裂纹尖端应力场具有奇异性,此时采用裂纹尖端附近某一点处应力大小来衡量材料在裂纹尖端是否安全是没有意义的,因此引入了应力强度因子尤来衡量整个裂纹尖端附近应力集中的强弱程度一。对于定的裂纹尖一端附近的应力场i:,其应力强度因子i是个有限量,它的大小与坐标选择无关,与裂纹一和构件的几何形状、作用力的大小和形式相关。般情况下的应力强度因子可由式(1.14)来定义?,即.='=Krlimlnra「j1^0^°(1.14)\_I=L'=limJlTcrr00“\r-^OI应力强度因子作为线弹性断裂力学的一个重要的物理参量,除了有明确的定义外,23[必须能够由计算和试验来测定]。常用的应力强度因子的计算方法有:复变函数法、边2425627283G][>2,]29,]l,32,33,34[][[P]界配置法、权函数法、积分变换法、奇异积分方程以及有限元法等。一些含裂纹的特殊构型对于,学者们己经求得了应力强度因子闭合形式的解,这些应力强度因子的解已经汇编成了几本应力强度因子手册供工程师参考使用基于力学原理的断裂准则可表示为:K>K(1.15)c其中用Xc为应力强度因子的临界值,当应力强度因子X超过材料对裂纹扩展的抵抗能力尺C时即认为发生断裂。尺和iTc的单位相同,其国际单位为工程上常用单位为。一一对于定材料而言,应力强度因子临界值尺C也是个常数,其反映了材料抵抗断裂破环的能力,通常也被统称为材料的断裂韧度,需要通过试验确定。对于I型和II型i^ii:裂纹,应力强度因子分别记为尺,Xn。I和n相应的断裂靭度分别记为,c和c1.3木材顺纹断裂木材是由中空的细胞组成的蜂窝状结构,细胞壁的主体是厚度最大的次生壁中层,°°?其微细纤维紧密靠拢,与纵轴有约呈1030的交角,这是木材呈各向异性的基本原6 大连理工大学硕士学位论文⑴一。因木材在般情况下表现出线弹性的本构关系,在进行断裂力学研究时,通常将木材当作正交各向异性材料,并忽略早晚材性质差异引起的非均质性,从而采用线弹性断裂力学分析。木材的3个弹性主轴分别为纵向L轴,径向R轴和弦向T轴(参见图1.2),3个弹性主轴相互垂直,每个主轴方向的弹性常数以及强度互不相同。木材清材试件的顺纹??抗拉强度可达120MPa150MPa,为顺纹抗压强度的23倍;木材清材试件的横纹抗一38?[]3040。拉强度比顺纹抗拉强度低得多,般只有顺纹抗拉强度的1/1/T/【.Z■乂图1.2木材的3个弹性主轴Fi.1.2Threerincialaxisofwoodgpp木材中裂纹的扩展受到材料弹性主轴方向函数的影响一,通常并不定沿着其初始方39[]向扩展,而是沿着纤维方向扩展。Wu最早将木材作为正交各向异性材料开展了断裂力学研宄,并应用线弹性断裂力学预测了含裂纹木材的强度根据Wu的理论,将应力强度因子概念用于各向异性材料的断裂预测需基于以下假定:(1)裂纹相对于对称主轴的方向是固定的;(2)各向异性情况下应力强度因子的确定方法与各向同性材料相同;一41重合[](3)临界裂纹扩展方向与弹性主轴方向之。Wu认为当裂纹与弹性主轴重合且其沿着初始裂缝方向扩展时,各向同性材料的断41[]裂力学可用于各向异性材料。一6木材断裂模式般由两个弹性主轴字母组合代表种断裂模式,、RL、,共计即TLLR、LT、RT、TR,其中第1个字母代表裂纹平面的法线方向,第2个字母代表预期的裂纹扩展方向(参见图1.3)。顺纹断裂为TL模式和RL模式。,7 层板胶合木层间断裂10度研究:^TLRLLR?1?TRRTLT图1.3木材的6种断裂模式Fi.13Sixrincialcrackroaaionssemsofwoodg.ppppgtyt断裂朝度是表征材料对裂纹扩展抵抗能力的重要力学参数in对辖射松。Kg等人(Montereyine)的断裂朝度进行了试验测定,发现其横纹断裂軔度比顺纹断裂軔度大p42[]10倍。邵卓平等人也对杉木进行了横纹断裂的三点弯曲试验,发现其具有很高的抵[“]抗横纹断裂的能力。木材顺纹长木质纤维的构造特性决定了其具有很高的抵抗横纹断裂的能力,含横纹裂纹的木材构件不会因为裂纹尖端应力奇异性而发生低应力破坏,可以采用传统的材料力学强度准则对其进行评价。由于木材在树木生长过程中和加工过程中形成的裂纹和缺陷大都在纤维方向上,而木材又在沿纤维方向上对裂纹扩展的抵抗能43量是最小的[]因此,木材断裂研宄主要集中在顺纹断裂方面。,一一在其他材料的研究中断裂初度般被认为是个材料常数,但是随着材料试件的几一441何形式或尺寸的改变],其某些力学性能参数可能会发生定的变化,为了确定断裂朝一一度是否是木材的个材料常数,些学者开展了木材断裂靭度尺寸效应的研宄。Boatriht等人在南美松inusradiata)顺口四g(P纹断裂的非对称加载的单边切点剪切(4SENB)试验中,设计了厚度在Imm到50mm之间的9种试件,发现当试件厚度大451]0mm-于,ST。StalTsch1时试件厚度对其■ic值无显著影响nzegg等人对不同朝带长度和试件厚度的云杉(Spmce)试件进行了顺纹断裂的楔入勞拉(WS)试验,发现当軔带长度和试件厚度分别大于70mm和30nim时,试件尺寸满足平面应变和断裂过程区的要_求,其顺纹断裂值无明显尺寸效应。Moura等人采用了6种几何相似但尺寸不同的试件对挪威云杉CPiceaabiesL.)顺纹断裂的单边切口三点弯曲(3SENB)试验进行了/:数值模拟,得到了当试件尺寸满足断裂过程区要求时,该种木材顺纹断裂显ic值无明47的尺寸效应的结论[]。邵卓平采用6种厚度在10mm到45mm之间的杉木和马尾松试件进行了TL断裂模式的紧凑拉伸试验(CT),试验结果也表明当试件厚度大于10mm时3943这两种木材值受试件厚度变化的影响不大[,]。王坚伟对3种不同厚度的杉木和-8- 大连理工大学硕士学位论文--S.rinr),.PF.材(SpucePeFi试件进行了顺纹断裂的紧凑拉伸试验(CT〉也发现了两种48[]木材的iTic值随着厚度的变化无明显变化的现象。Yoshihara等人采用5种不同高度的试件,对西部铁杉Westernhemlock)开展顺纹断裂的双悬臂梁(DCB)试验和端部切(口三点弯曲(3ENF)试验,发现了Gc值随着试件高度的变化没有显著变化的现象,并49[]?证明了经裂纹长度修正后的断裂朝度Gc值也无明显尺寸效应。上述研究表明,当满足平面应变状态与断裂过程区要求时,木材顺纹断裂初度无明显的尺寸效应。1.31.木材顺纹断裂行为的研究G常用的测定木材顺纹断裂軔度ic值的试验方法是双悬臂梁(DCB)试验。邵卓平对杉木、云杉和杨木进行了顺纹断裂的DCB试验,发现木材顺纹断裂是否出现纤维桥接现象与木材种类相关,杉木没有桥接现象而云杉和杨木有不同程度的桥接现象釆用经典梁理论计算断裂韧度时忽略了裂纹尖端变形和横向剪力的影响,误差较大。为此Yoshihara在釆用DCB试验测定西部铁杉(Westernhemlock)TL断裂模式的G一,c值时,提出了通过附加应变片测定试件某点应变进而计算G,c值的柔度组合法。51利用该方法可以不需要测定木材的弹性常数即可简便地得到其顺纹断裂?[](,(:值,而且通过相互对比发现柔度组合法的计算结果与附加裂纹长度修正的Timoshenko梁理论的49[]计算结果以及采用修正的虚拟裂纹闭合有限元法所得的计算结果接近。,^玲[?图1.4特殊试件截面Fig.1.4Specialcrosssectionofspecimens^ ̄curv)为了得到中密度纤维板TL断裂模式的裂纹扩展阻力曲线Ue,Yoshihara还采用了特殊截面(参见图1.4)的试件进行DCB试验,该形式的试件截面保证了裂纹-9- _.层板胶合木层间断裂軔度研究.的稳定扩展,试验结果展现了裂纹两侧的木质纤维桥接作用引起的扩展断裂朝度随着裂纹扩展而增大,最后趋于稳定值的现象(参见图1.5).??>.6「'■S/3I1^2iTL断裂模式■”L_1__._iI.__I._I ̄ ̄“?o‘“020406080100120140Aa/mm52[]图1.5I型断裂DCB试验裂纹扩展阻力曲线-52[Fi.1.5Ticalcurweobi]gypmodeIRtanedbytheDCBtestMoura等人将海岸松(Maritimepine)RL断裂模式的DCB试验与粘聚力损伤模型(Cohesivedamagemodel)相结合,把C?ic分为裂纹扩展和纤维桥接(Gib)两部分(参见图1.6),利用该方法可不需要监测裂纹的扩展情况,只通过分析试验的荷-53[]载位移曲线即可简便而准确得到0。10值‘./ ̄—+,GicGiGibftj“裂纹扩展)l\>..?s.......V.ftb(纤维桥接)\..:x.乙—Zi80515c裂纹张开位移1’?53图[】1.6粘聚力损伤模型53F[].mlig.16Cohesivedamageode在DCB试验基础上Morel等人开展了锥形双悬臂梁(TDCB)试验,特殊形状的锥一形悬臂试件使得裂纹的扩展速度与加载点速度致,保证了I型断裂裂纹的稳定扩展,--10 大连理工大学硕士学位论文[54进而可以得到 ̄]i?curve,并且当试件设计合理时可以使得试件柔度与裂纹长度呈线性[55]关系,那么值就可以不需测定裂纹的长度而直接由荷载值计算得到。(?三点弯曲(ENF。测定木材顺纹断裂,ic值的试验方法主要是端部缺口3)试验为研宄试件跨高比对3ENF试验测定G,ic值的影响,Yoshihara采用了8种不同跨高比的西部铁杉(Westernhemlock)试件进行试验,并应用式(1.16)和(1.17)计算Ghc值,发现含有弹性模量的式(1.16)的计算结果稳定,而含有柔度值的式(1.17)在跨高比6[5]较小时计算结果偏大(参见图.7)。因此,3ENF,(1.17)计1在试验中当采用式算Giic值时,通常需要限定试件跨高比,最小跨高比约为12。-0.6r0.6二9CP-0.5王「_y0.5G--与0.4土£0.4■?STTt3--一0.3II一0.3TTIT“"-0-20.2?:o1^i:;1IJJ::--0.10.1III■II■I‘I■丨III■??■?_I■■QQQQ5101520253051015202530试件跨高比试件跨高比65图【]1.7断裂朝度Giic与试件跨高比关系F.1.7eaoneweenheacureouhnessGlculaedbeacheuaonandesaneigRltibttfrttgnccatyqtithp/dpth56[]ratio9户C(1.16)lie?3229CpG=(117)C.?3'3'252Z+3a()柔度组合法同样可用于计算木材3ENF试验的顺纹断裂軔度。Yoshihara对西部铁杉(Westernhemlock)TL断裂模式进行了3ENF试验,并采用柔度组合法、附加裂纹长度修正的Timoshenko梁理论和修正的虚拟裂纹闭合有限元法计算所得G?c值,发现3一49,57[]种方法计算结果致,验证了柔度组合法计算G。ic值的准确性Silva在海岸松(Maritimeine)的3ENF试验中分别采用柔度组合法和其他2种修p正的梁理论以及虚拟裂纹闭合有限元法计算了Gi,c值,发现柔度组合法与有限元法所得-?11 层板胶合木层间断裂搬研究.结果吻合,验证了柔度组合法计算G?c值的准确性,同时发现裂纹表面的摩擦力对断裂一58[]朝度G?C值计算结果影响很小。Yoshihara采用特殊试件截面(参见图1.4)的试件对中密度纤维板(MediumDensityFiberboard)进行了TL断裂模式旳3ENF试验,得到了其II型断裂的curve,发现了由于裂纹两侧的木质纤维桥接作用引起的扩展断裂靭度随着裂纹扩展而持续增大的现52[]象(参见图1.8)。10「Z8一.、[yA6-一'ZZ、4-HZO?初试非线性2TL断裂模式\IIIII.IJ0I01020304050Aa/mm52[]图1.8II型断裂3ENF试验裂纹扩展阻力曲线52-[F].1,8Tlcuobiedbhe3ENFesigypicamodeIIRr\etanyttt端部荷载勞裂(ELS)试验可以保证木材顺纹II型断裂试验中裂纹的稳定扩展进而_"Curve。Siii)得到i?lva等人对海岸松(Pnuspnaster进行了RL断裂模式的ELS试验,一/?"K,,并将试验得到的:urve与虚拟裂纹闭合有限元法所得的对比发现两者致59[验证了ELS]试验测定木材断裂軔度的准确性。常用的测定I型断裂尤,C值的试验方法有单边切口试件的三点弯曲(3SENB)试验和楔入勞拉(WS)试验。Watanabe等人米用sandwich型试件对日本柳杉Crtomeria.(ypjaponica)进行了3SENB试验,通过改变中间块木材的预制裂纹形式得到不同的木材断裂9)-模式(参见图1,TLTR.测定了该种木材从断裂模式变化到断裂模式下的荷载裂纹张幵位移曲线(参见图1.10),发现i:TL断裂模式变化到TR断裂模式逐渐,c值从111)TL降低(参见图.。这表明了该种木材断裂模式起裂勒度高、脆性强的特点;同时也解释了该种木材在自然干燥时容易出现与TR断裂模式类似的裂纹但并不持续扩展_的现象。--12 大连理工大学硕士学位论文.B_°°0TL30()■__。。6090TR()6[印图1.9sandwich型试件断裂模式的变化F'60.1.9Thiddli[]igemeecewoodsanlevariationofsandwichsecmensipgp■— ̄。0TL()「力■:歡_II‘10IIII0.00.51.0..1520COD/mm1-图.10不同断裂模式的荷载裂纹张开位移曲线60F-[]ig.1.10LoadCODcurveofdifferentfracturesstemsy0.30「'...T...-0.25M_2-B0.20-II0.15I“-^0.10-0.05■■Li.1‘1‘L0.000306090裂纹面与。LT平面夹角/图1.11不同断裂模式的断裂軔度值6o]Fi[g1.11Fracturetounessofdifferentfracturesstemsghy--13 层板胶合木层间断裂初度研究Ardalany也米用sandwich型试件对福射松(Montereine)锯材和由福射松制成的yp单板层积材(LVL)进行了顺纹断裂3SENB试验,发现单板层积材的尤显高ic值要明于福射松(Montereine)锯材A^,ypic值并在试验实测值的基础上采用有限元法对一该试验进行了数值模拟,结果显示数值模拟预测的临界荷载与试验实测值致,表明线61[]弹性断裂力学方法可以有效的预测带裂纹木梁的断裂破坏。王坚伟采用sandwich型试件的3SENB试验分别测定了杉木和S.P.F材S-P-(ruceineFir)TLRLiT断裂模式和断裂模式的,pic值并结合声发射技术确定临界荷载尸C以保证测试结果的精确性,然后与同种材料的紧凑拉伸)(CT试验结果进行了48[]对比,发现两者吻合较好。Reiterer等人采用WS试验研究了含水率变化对云杉(spruce)TL断裂模式仏c值的影响,发现当木材含水率从7%增大到12%(标准含水率)时其l显著减小,木ie值材含水率从12%(标准含水率)增大到55%(含水率高于纤维饱和点)时其K,c值保621]持稳定(参见图.12)。106.「:0.5^-土T1T-£0-^“I民.]2丄1丄蓉0-::.3截-由^j-0.2I—I_I_I_I_I_I_I_I__IJ_I_I_I_I_1_I_I_I_I_IQ*51015202530354045505560木材含水率/%图1.12不同含水率木材断裂軔度^F]i.1.12FracturetouhnessATofdiffetitt#ggicrenmosureconenSchachner等人在应用有限元法对云杉(spruce)3SENB试验进行数值模拟时,在裂纹尖端布置退化的1/4节点的奇异单元,通过节点位移值求解值的方法。他们发现木材纵向长纤维的构造特性使得从RL断裂模式到LR断裂模式其断裂朝度不断增大,一63[]R这现象证实了木材突出的各向异性性质eiterer。等人将上述有限元方法用于像木(oak)和水曲柳(ash)顺纹断裂WS试验的数值模拟中,发现两种木材的RL断裂模式的尺IC值均高于TL断裂模式的Kic,这表明木材径向木质纤维对顺纹断裂朝度具有64[]强化作用。RetStanzl-Tierer和schegg采用同样的方法对比了软木云杉(spruce)与硬-4-1 大连理工大学硕士学位论文木桤木(alder)、橡木(oak)和水曲柳(ash)顺纹断裂籾度的试验结果,发现硬木的木质纤维长度较软木的短,,,因此纤维桥接作用较弱导致裂纹起裂后不稳定扩展而软木云杉(spruce)有稳定的裂纹扩展;同时还发现同种木材顺纹断裂值随着其密度增大而增大;并再次验证了木材径向木质纤维强化作用引起的RL断裂模式i^于,c值高一6566TL断裂模式的值这[,】结论(参见图1.13)。■_TL断裂模式1.2[Y■□RL^断裂模式.V.。:-1:_0.40.2450500550600650700750木材密度'/(kg/m)?■[]图1.13不同种类木材的断裂軔尺C(从左到右度:云、桤木、1杉橡树、水曲柳)F65i.1.13Frachnfdiffesifri[]turetouessATorenwoodecesomleftorh:srucealderoakashggictp(tgtp,,,)测定II型断裂軔度J^uC值的试验方法主要是非对称加载的单边切口试件四点弯曲(4SENB)试验。Yoshihara采用非对称加载的4SENB试验与虚拟裂纹闭合有限元法相结合的方法计算了北美云杉(Sitkaspruce)TL断裂模式的G?c值与/^iic值,并与3ENF-试验结果进行了对比,发现非对称加载的4SENB试验中裂纹长度a的变化对荷载位移 ̄曲线的影响很小,因此不适宜采用柔度标定法计算其G,但是当(a/H)值在07085,ic值..之间时两种方法测得的值是近似的,验证了非对称加载的4SENB试验在测定木材67[】顺纹断裂&,c值的可行性与准确性。Susanti等人则对贝壳杉(Agathissp.)短裂纹锥形截面的试件进行了TL断裂模式非对称加载的4SENB试验,发现含短裂纹试件的剪切控制应力明显低于无裂纹试件,并且在进行裂纹长度修正后非对称加载的4SENB试验所得/Gic值与3ENF试验所得值吻合,再次证实了非对称加载的4SENB试验测定K,iic值的可行性与准确性也发现在进行裂纹长度进行修正后即使无裂纹的试件也可应68用断裂力学计算其破坏的控制剪应力[】。1.3.2木材胶结面断裂行为的研究工程木中存在大量的胶结面,其中气泡、残澄和收缩过程产生的裂纹处都可能出现应力奇异性。众多学者采用能量原理对木材胶结面的断裂问题进行了研宄。--15 层板胶合木层间断裂韧度研宄CCBD-lfi对于木材胶结面的I型断裂,onrad等人釆用D试验测定了花旗松(ougasr)-C?Flexohiiti酷酸树脂胶合剂胶结面的,racrnnic值并结合了柔性版印刷技术(gppgtechnique)分析胶粘剂微滴直径和间距对该种胶结面的G,c值影响,发现了其随着微滴直径和粒径与间距的比值增大而增大,随着胶合剂微滴间距增大而减小的现象(参见图_dl1.14和图1.15)。Vge等人也采用DCB试验方法,研究了纤维素添加剂对于云杉(spruce)的脲酸胶合剂胶结面Gic值的影响,发现纤维素的加入可以增强该种胶结面7?[](?0。之Veie【uce)胶结面,值后g等人采用同样的方法研究了不同胶合剂对云杉(spr’断裂靭度的影响,发现聚氨酷胶合剂的云杉(spruce)胶结面Gic值大约是脲酸胶合剂的1.79倍,并得出了断裂试验能较好的区分胶合剂之间的差异,断裂軔度适合用来表征[71]胶合剂胶结性能的结论。300■■微粒间距l.Omm[■。O250/微粒间距.Smm'.s1/y■參150/。、。。.A■0:色_I.11I02004006008001000微粒直径.图114不同微滴直径的断裂軔度G,69F[]i.1.14FracturetounessGofdifferentdroledametergghicpti‘?400r、-l350■微粒间距.Omm□0.5mm.微粒间距2QQ250.言<200-X^.i5oy;i^"-10050、■oL_._-II0.00.00.20.40.6.81微粒直径/微粒间距比值图1.15不同微粒直径/微粒间距的断裂朝度69F[]i.1.15FracturetouhnessGofdifferentdroletdiameterdivdedbrolesacinggicpiydptpg--16 大连理工大学硕士学位论文对于云杉(spruce)胶结面的.II型断裂,Qiao等人则采用锥形端部缺口三点弯曲(3TENF) ̄curveTENF试验克试验测定了其Giic值和i?,3服了裂纹初始长度难以精确测定以及当初始裂纹长度小于0.35倍支座跨度时裂纹不稳定扩展的缺点。与TDCB试验类似,当试件形状被合理设计时可以使试件柔度与裂纹长度呈线性关系,从而简化了一值的计算,,,并将计算结果值并与数值模拟结果值相对比发现两者致验证了该试7273,[]验方法测定木材胶结面G?c值的可行性和准确性。1.33.层板胶合木断裂行为的研究,对于层板胶合木,其断裂行为可能是发生在层板内部的层内断裂也有可能是发生在层间胶结面处的层间断裂。层板胶合木的层内断裂实际上就是木材的顺纹断裂,相关的研究己经比较充分。但是,对于层板胶合木层间断裂行为的研究十分缺乏。层板胶合木制作过程中木质层板的胶结存在一定的随机性,那么层间胶结面可能出现两侧不同年轮角度的木质层板组合,Wang等人采用图像分析技术结合DCB试验对层板胶合木进行了层间胶结面断裂韧度测定,发现胶结面两侧不同年轮角度的层板组合会引起断裂軔度74一G值的变化【],但是由,ic于试验数据较少并未得到定性的结论。木材是种正交各向,异性材料,不同年轮角度的层板组合可能会对层板胶合木层间断裂韧度产生影响因此对该问题进行断裂力学分析是必要的。1.4本文的研究内容本文将层板胶合木作为研宄对象,通过对全部试件截面的数码照片进行图像分析确定层板年轮角度的代表值,并对年度角度进行了统计分析。然后,分别对不同年轮角度层板组合的层板胶合木进行了层板胶合木层间断裂軔度的双悬臂梁(DCB)试验、端部切口三点弯曲(3ENF)试验和单边切口三点弯曲(3SENB)试验,在试验的基础上对3种层间断裂试验进行有限元模拟,分析不同年轮角度层板组合对层板胶合木层间断裂朝度的影响一,为层板胶合木的断裂理论分析和工业制作工艺提供定的指导。1.4.1层板胶合木层间断裂朝度仏C和汉CI的研究对不同年轮角度层板组合的层板胶合木分别进行了I型层间断裂的双悬臂梁(DCB)试验和n型层间断裂的端部切口三点弯曲(3ENF)试验,以层板胶合木胶结面两侧层板年轮角度的差值为参量,研究了角度差值对两种断裂试验的临界荷载的影响,在得到两种断裂试验的临界荷载后基于经典梁理论分别计算了层板胶合木层间断裂朝度G,,c和G,ic进而研究角度差值对两种断裂试验的层间断裂勒度的影响。17 层板胶合木层间断裂朝度研究;应用ABAQUS程序分别对上述两种断裂试验进行了不同年轮角度层板组合的三维有限元模拟,采用ABAQUS自带模块计算了两种层间断裂试验在临界荷载下的/积分,在线弹性情况下的J积分即为层间断裂勃度(Jic和G?c,研宄不同年轮角度差值对于上述两种试验所得层板胶合木的层间断裂軔度G,ic和Giic的影响。最后将有限元所得结果与经典梁理论所得结果进行对比分析。此外,对不同年轮角度层板组合的层板胶合木进行了I型层间断裂的sandwich型试3SENB一件单边切口三点弯曲()试验,.基于试验所得到的临界荷载对该试验进行了系+列不同年轮角度层板组合的三维有限元模拟并通过回路积分法计算得到层间断裂韧度(7,,C研究不同年轮角度差值对于该种试验所得层板胶合木的层间断裂韧度Gic的影响。最后,将双悬臂梁(DCB)试验和单边切口三点弯曲(3SENB)试验所得层间断裂铜度G.c的有限元计算结果进行对比分析。1.4.2层板胶合木层间断裂钢度K、和/r?。的研究。应用ABAQUS程序对不同年轮角度层板组合的层板胶合木分别进行了I型层间断裂的双悬臂梁(DCB)试验和II型层间断裂的端部切口三点弯曲(3ENF)试验的三维有限元模拟,基于在临界荷载下裂纹尖端节点的应力场的分析釆用应力外推法分别计算iiCZiTH了其层间断裂韧度ic和C,研宄不同年轮角度差值对于上述两种试验所得层板胶合木的层间断裂韧度尺iT。IC和iiic的影响此外应用ABAQUS程序对不同年轮角度层板组合的层板胶合木分别进行了I型层间断裂的sandwich型试件单边切口三点弯曲(3SENB)试验的三维有限元模拟,基于对临界荷载下裂纹尖端的应力场的分析采用应力外推法计算了其层间断裂軔度Kic’研究不同年轮角度差值对于上述试验所得层板胶合木的层间断裂初度的影响。最后,将双悬臂梁(DCB)试验和单边切口三点弯曲(3SENB)试验所得层间断裂軔度K、c的有限元计算结果进行对比分析。--18 大连理工大学硕士学位论文2层板胶合木层间断裂试验本文对不同年轮角度层板组合分别进行了I型层间断裂的DCB试验和3SENB试一验,II型层间断裂的3ENF试验。本文试验所用的材料为工厂批量生产的同批层板胶一合木,材料均性较好,因此忽略了胶粘剂性能、胶结面状况、粘结工艺和胶层厚度等其他因素的差异对于层板胶合木层间断裂初度的影响,重点关注胶结面两侧年轮角度对于层板胶合木层间断裂勒度的影响。己有的层板胶合木层间断裂研宄并未充分考虑木材的正交各向异性性质,而是忽略的木材径向与切向弹性性质的差别,因此本研宄将年轮角度作为影响参数,主要考虑层板年轮角度差值(即充分考虑木材径向与切向弹性性质的差别)对层板胶合木层间断裂初度的影响。2.1木材弹性常数测定层板胶合木层间断裂试验使用的试验材料是由大连双华木业有限公司生产的同一><><批层板胶合木00mm100rnm60mni。,层板基材为樟子松,尺寸全部为10该批层板胶30合木的层板组成形式是同等组合,其在厚度方向上由两层mtn厚的木质层板胶结而成,并且在长度方向上无指接,在宽度方向上无拼接。层板胶合木制作采用的胶粘剂是--:KR134AX200)。21水性高分子乙炼基聚氨酯胶粘剂(主剂,固化剂:图.为本文所用的层板胶合木以及其截面示意图。层板胶合木平均含水率约为9.6%,平均密度约^58358mo为.kg/I:图2.1本文采用的层板胶合木lFi.21TheGluamusedinthetestsg.--19 层板胶合木层间断裂铜度研究;一本研宄忽略了木材早晚材之间的性质差异,将其作为种线弹性的正交各向异性材料进行研究。正交各向异性材料具有3个相互正交的弹性对称面,其独立的弹性常数有759[】个,在进行层板胶合木层间断裂试验的有限元模拟时需要输入木材特定的9个弹性常数。本文采用了单向轴压试验和变跨距三点弯曲试验测定了层板木材樟子松的12个弹性常数,其中包含3个杨氏模量及3,6、Er.El,个剪切模量GLR、GLT.GRT个泊松比/化R、/^LT、/<RT、//RL、/<TL、"TR。;本文测定木材弹性常数所用的为清材试件,试件不包含裂纹和木结等缺陷,也不包一含胶结面,木材性质均性相对较好3。采用单向轴压试验可以测定层板木材个方向上£EE22,的杨氏模量1、R、JT和泊松比《lr、ult^/RT以及剪切模量GRT。如图.所不单;/j20><2())><40mm()向轴压试验试件的尺寸全部为rnrn(长)rnm(宽高,其中(a)型、(b)型和(C)型试件高度方向分别为木材的T向、R向和L向,(d)型试件高度方°向为年轮45方向,每种类型试件各3个。荷载沿着试件长度方向施加,加载速度为Imm/min,并保证加载时试件处于线弹性范围内。隠於_隱晒___abc(d()()())图2.2单向轴压试验试件F.ig2.2Samplesofthecompressiontests)2对相互垂直的应变片-(a)型、(b)型和(c型试件表面分别贴有,通过应力一应变曲线可以得到杨氏模量,通过对应变片应变的比值可以得到泊松比。从应变片的布置可以知道,(a)型试件可以测得Et’;/lt和/^T,(b)型试件可以测得Er,fiiR一和<TR,(c)型试件可以测得五,T和Z。(d)型/!^/TL/RL试件上贴有对相互垂直的应变°-片,根据其应力应变曲线和应变的比值可以分别得到45年轮方向上的杨氏模量五45和泊松比//45,那么剪切模量GRT可以根据式(2.1)计算:-__^Q一RT0/,\211+(")45(2.1)-20- 大连理工大学硕士学位论文采用变跨距三点弯曲试验可以测定木材的剪切模量Gui、Glt。试件尺寸为300mm><20m><20mm(高(长)in(宽))L向,宽,其中长度方向为木材的度方向和高度方向分别为木材的R向和T向。荷载沿着试件高度方向施加,加载速度为Imm/min,并保证加载时试件处于线弹性范围内一。当荷载施加方向与木材R向致时可以测得剪切模一量(?LR,当荷载施加方向与木材T向致时可以测得剪切模量GLT。变跨距三点弯曲试验要求试件的跨度与高度的比值应在5-20之间.5,在本文中分别选取了240mni,2l0nim>180mm.150mm,120mm这5种跨距进行三点弯曲试驗,根6[7]据得到的不同跨距下木材的抗弯弹性模量MOE,采用式(2.2)求出对应的剪切模量:21.2A(今>G=—^(2.2)a()moeAPfi=MOE(2.3),通过单向轴压试验所测得的杨氏模量、泊松比和剪切模量以及通过变跨距三点弯曲试验所测得的剪切模量见表2.1和表2.2。试验测得木材弹性常数结果的变异系数不大,因此将上述两种试验所测定的木材弹性常数用于第3章的有限元模拟是合适的。表2.1试验木材层板弹性常数Tab.2.1Elasticconstantsofwoodlaminaeinthetests-弹性常数序a?Et/MpjEg/Mpafit/MpaGRi/MpaGbR/MpaGtr/Mpa110779.681390.86634.1352.81734.35632.29212038581164188...735.1563.760351731.45311572310053720484432717477.5.2...7621.4////774.65705.38平均值11463.601186.76696.59slw757.56672.72标准差636.48193.805458^30553.99 ̄变异系数sJsiMiil?5J3-2-1 层板胶合木层间断裂軔度研宄表2.2试验木材层板泊松比T'ab.2.2Poissonsratioofwoodlaminaeinthetests泊松比序号ULRi^"TL/fVTMRl/RTMTRy10.3970.62600380.8330.0250.406,20.3660.6690.0270.6840.0330.422303460.6190.0450.77800320.333.. ̄平均值0^0.6380^0?^0.0300.387标准差0^0.0270.0090.0750.0040.047变异系数^^24?M^14^5312^22.层板年轮方向角度代表值在试验中对3种断裂试验的所有试件截面进行了拍照,然后通过对图片的处理和分析确定胶结面两侧层板木材年轮角度的代表值。如图2.2所示,将试件截面图片导入AutoCAD2013程序中,然后在胶结面中点处分别做上下层板年轮的切线,取切线与胶,结面的夹角作为年轮角度代表值。为便于统计与分析并根据对称性,试验中将较小年轮角度的层板定义为胶结面的下层板,较大年轮角度的层板为胶结面的上层板。im7图2.2试件年轮方向角度定义Fig.2.2Examplesofthephotographsusedtoidentifythegrowthringorientationsofspecimens一由于试验所使用的层板胶合木为工厂批量生产的,其层板组合般遵循了破心下料的原则,因此层板年轮角度组合以及其差值分布存在随机性。DCB试件、3ENF试件和3SENB试件的总数目分别分76个、90个和70个,其层板角度分布统计分别见图2.3、°°图2.4和图2.5。DCB试件的《1和6(2平均值分别为20.21和38.71。3ENF试件的cti和°°°°ct2平均值分别为25.49和40.03。3SENB试件的ai和平均值分别为16.29和26.27。-22- 大连理工大学硕士学位论文20r60?■rrs■m■?狗■15■.40、?:■■1。.■顏?、i?<f批?—V ̄ ̄:B20门JIIj「■■._5??.■■门H■■?.■■■1■- ̄■“■”0Un、义i'p'--— ̄‘ ̄■ ̄■■■”■■_—J ̄■ ̄ ̄I—I_.I.111li1.1U1.11.1.1I11Irriq1f..010203040805060700102030405060。试件编号丨(a)下层板角度cti10012「「■_80-I■9-n-60-卿:PDbn.:■—_睡>40■!^^?、■*?、rr__?■?■:■:」^20■■3门,二■‘■_■?■■0,1■■」_贝n0L_01020304050607080900102030405060708090‘。"丨试件编号1(b)上层板角度a2图2.3DCB试件层扳年轮角度a)(下层板角度…(b)上层板角度a2'F.2.Rig3eresentationofDCBecpsimensrowthrnorientationsanlpgigesatheowerlaminaeag()l\(b)theuerlaminaeappi20r6■_0p!■--15“-■\V!?"4■0■?”读■■■■?■、gio-h■■:_--握;■■--rannn20.-■;!.■■"■■■5-ram■■‘-n ̄.f.I■!】Hili1.111IIr^.lim.qI01020304050708090100060102030405060。a试件编号j(a)下层板角度ct丨-23- 层板胶合木层间断裂勒度研究21001「「80-I‘■9rn_■-■60■:n?■■■■■「■■■、■■、--6■‘|门“?--■頃40—??nnS■■■-■‘■■、■■-20-“3门''0nn■卜II11IIII1IIIIIIII11,■■■■q01020304050607080900102030405060708090。试件编号aj(b)上层板角度0C2图2.43ENF试件层板年轮角度(a)下层板角度(b)上层板角度a2F'2i..4Representation3ENFsecnrowrnrenangofpimesgthigoittiosanglesathelowerlaminaeab\()()heuerlamnaeatppii60r25r■50■■20■目40-■■15■■?30i■■栽9^■-_暖=■20.",,10一"■?h5■[(j-10■■■O门—"门nrr-n■IJ■??隱■?-■(J""“:、1■■■■■■lil{.fIIU11illIIfL賴I■tiAl....l..0102030405060700102030405060。Ol试件编号L\(a)下层板角度ai7012■[60-■1。-■■?■-■.50.?.?_■广■■‘-nn■Rrrr,400i__■??tT?载6iir_网:pj3Q■““—=^■>‘H产暖3??■20?4一一■I1II■“■■■、■■‘_nn:10?^门^■■■r■?睡?\?■■门1■__II11ilIIiII1.I.1__1._II_.Q1-.1l-■104050607000203010203040506070试件编号(b)上层板角度《2图2.43SENB试件层板年轮角度(a)下层板角度a()上ib层板角度F'i.2.4Reresentationof3SENBsecngimensrowthrorentaionsanlesatheloppgigitgwerlaminaeaib()()eulinaethpperamai-24- 大连理工大学硕士学位论文DCB试件、3ENF试件和3SENB试件的层板角度《丨与a2的关系见图2.5。,80r80r■70-70-■■6060■‘??-50-50.:■"。=^40V?"““_■.-■-30.30‘■■.?2020."■“■?_S■■-1010-■■0■-^00102030405060708000304070102506080a,.(a)DCB3ENF试件(b)试件70r60-■■■■■50-?■_■■■_■■■<*■■、40-■■?cT譯30-.■.??■■_20??■■10-:050102030405060oa/,(c)3SENB试件图2.5层板年轮角度a与《2(a)DCB试件()3ENFc)3SENBib试件(试件'-F..ig25RepresentationofspecimensrowthrinorientationsanlesaaaDCBsecimensb3ENFggg\2()p()specimensc3SENBsecimens()p为了便于分析和处理数据-,本文将角度《2《《差值,那么丨定义为角度2与角度1的角度《2-可以用来代表胶结面两侧层板年轮角度的差异性。DCB试件、3ENF试件和1°°°ENB263S试件的角度分布见图.,其平均值分别为17.29、14.63和9.99。-25- 层板胶合木层间断裂朝度研究.6020「「■■■50?.W■■■■15■40-|:?^■■親30inn10-■與於」■_?Q■?1F■1■■20.■: ̄n■■■■'■■■.■■5■■■■■■10■■■_‘■‘‘网…i...nn.II11II11lnn』n,p010203040506070800102030405060试件编号?尸2-1(a)DCB试件50r30「-40■?25T■“■‘20:、30..“■■■■:■15■■續n20■■■■.^-__!!J__10■,n"—n—-"I-■■■fr1■10门门■■5'■■■■_■■0'■‘■“「■_■JI.II^0020406080100010203040。《试件编号/2■丨(b)3ENF试件50r30r—■2540-:■20、30-■■p3■錢.15。■■、20..展■?"?110门?■‘■■■—^ ̄■—^10■,■=““■■■■c■,''?■??ok::::IIIII.Iin.rin.n.q010203040506070010203040试件编号oaj(c)3SENB试件图2.6《2-(a)DCB3ENFc)3SENB层板年轮角度1试件(b)试件(试件F.2.6Reresentaionofanesa-aDCBib3ENFs3SENBsisigptglsecmensecimenscecmenix()p()p()p-26- 大连理工大学硕士学位论文23DCB.层板胶合木层间断裂试验2.3.1DCB试验概况本文采用DCB试验测定层板胶合木层间I型断裂軔度G,c,DCB试验试件的几何形状如图2.7所示,试件长度I为300mm,高度2//为60mm,厚度5为30mm,预制裂纹总长(《+如)为120mrn,其中flo代表加载点距试件边缘的距离,为20mm。预制裂,2nim118,05mm缝位置处于层板胶结面处先使用厚度为的带据预制mm再使用厚度为.2mn,120,(L)。的刀片切割i得到总长为mm的尖锐裂纹裂纹的初始方向为顺纹方向nP应变片2^/间胶结面一―-"士/〒口:丨IaaT?)WL_啦1.HI舰2-图.7双悬臂梁(DCB)试验模型Fi.2.7SamleofthedoublecantileverbeamDCBtestgp()DCB试验采用的试验机为济南时代试金试验机有限公司制造的WAW-300型微机控制电液伺服万能试验机(参见图2.8),并为DCB试验订做了专用的加载设备(参见2图.9)。试验采用位移控制的加载方式,加载速度为2mm/min。为提升测量数据的精lN的±P和确度,分别附加了Ok拉压力传感器和量程为50mm的位移计进行荷载加载点-27- 层板胶合木层间断裂朝度研究;IMC。DCB线位移(5的测量,并采用数据采集器进行数据的釆集试验中裂纹随着荷载P下0?的增加是持续扩展的,试验在荷载值降至时停止。2-图.8WAW300型微机控制电液伺服万能试验机F--ig.2.8WAW300typecomputercontrolelectrohydraulicservouniversaltestingmachineHgTvfi.^w.脅^”M-_:f丨,__图2.9DCB试验加载设备Fig.2.9LoadingequipmentsofDCBtests-28- 大连理工大学硕士学位论文2-.10DCB图为试验典型的荷载加载点线位移⑷曲线,从该曲线可以看出,一P一裂纹在起裂之后存在定程度的扩展,荷载值在裂纹起裂后段时间内继续增加,达 ̄到峰值荷载后缓慢下降,直至荷载值P降为0。本文在DCB试验中PS曲线线性段的终点所对应的荷载定义为临界荷载Pc’进而计算层间I型断裂初度值。由于试验所得一试验曲线的形式是致的,并且由于需要在试验曲线上指明临界荷载的定义,为了便于一阅读,因此试验给出了条具有代表性的试验曲线211CB。图.为D试验曲线线性-纵向应变-。段对应的荷载、PeO曲线600r5。。…4。。.RS300-广/1200/100I*0—‘‘‘‘’0510152025S/mm2-图10DCB.试验荷载加载点线位移(尸4)曲线F-i.2.10Examleoflliliiligptheoadoadngnedspacement(PS)curvesobtanedfromtheDCBtest600r500.400.300■I200-?100*^.000.050.100150.200.25,tJW图2--.11DCB试验荷载纵向应变(PSL)曲线F--i.2.11ExamleofloadlonitudinalstrainPeicurvesobaedomegptheg()tinfrthDCBtests图2.12为DCB试验试件破坏形式示意图,在荷载作用下初始裂纹基本沿着胶结面处做自相似的扩展,为层间I型断裂破坏。由于层板胶合木制作时使用的为水性胶粘剂,01mm左右)并且胶层厚度较薄(.,胶层破坏没有明显现象。-29- 层板胶合木层间断裂初度研究--—mlI1Ii11Mr1TfIiriHfii1nffTfri1fTTini;av??—am—!^■‘: ̄ ̄ ̄^.图2.12DCB试件破坏形式Fi.2.12DamaeformsofDCBsecimensggp2.3.2DCB试验结果2-。图.13DCB试验P2.13可以看出由于木为的临界荷载ic与角度《21的关系从图材性能自身离散性较大使得户,(:的分布较为分散,/^(:的平均值和标准差分别为474.51^和92.05N。在不考虑《2-随机性引起不同角度差的试件数目不平均的情况下对1分布的一-Pa.,可以发现P-,当c-添加直线型趋势线《t\ci\IC随着2】的增大有定的增大趋势2i°。从0增大至57时趋势线上所对应的Pic值从448.99N增大至539.33N,约增大了20%。800r700■■‘■600■??■■?■■■■“■一5500?_■':-■^00?4■.■■;:■-■_—■■.?300■200■,■100?1—————————■—0?‘I‘I‘”■—‘——■‘‘0102030405060。<hJ-图2.13DCB试验《2.1与临界荷载Pic关系FifrCBti.2.13RelationshibetweentheanelaiAandthecriticalloadPobtanedomtheDestsgpg\c77[]G-K,可采用Irwin,根据能量原理断裂试验中材料的c值ies关系式计算计算式如下:-30- 大连理工大学硕士学位论文G=…(2.1)c()告盖丨C:荷载式中Pc为临界,5为试件厚度,C为试件的柔度,^;为裂纹长度,flc为临界裂纹长度。ST的标准试验方法之一美国材料与试验协会(AM)将柔度标定法作为计算Gc值,因此众多学者采用该方法计算木材的顺纹断裂Gc值。试件的柔度值可以基于经典梁理7879一[,论4十算,也有些文献对经典梁理论进行了裂纹尖端变形修正或者考虑84了横向剪力的修正[:"等。本文中忽略了裂纹尖端变形和横向剪力对计算结果的影响。基于经典梁理论,IG可由式(22)层板胶合木型层间断裂钥度.,c12jGy(2.2)图2.14为基于经典梁理论计算的DCB试验Gic与0(2-1关系,从图2.14可以看出由,G0061于木材自身离散性较大引起的Gic分布较为分散ic的平均值和标准差分别为1.22J/m和36.81J/m。在不考虑《2-i分布的随机性引起不同角度差的试件数目不平均的情一,可以发现(?的增大有况下对添加直线型趋势线,,C随着定的增大趋势°°22当《2-0578901J2753J/m1从增大至时趋势线上所对应的Gic值从./m增大至1.,约增大了43%。并且从公式(2.2)可以确定,0应/^,(:与所对0的平方值是呈线性关系的。200?「■■■■■■-■150*■■霧■■■■I■5‘100■*O??■“”"■50■Q1‘‘‘■0102030405060。?丨2-1图2.14DCB试验与断裂韧度Gic关系Fi.2.14Relationshibetweentheanel?2-andfracturetouhnessGobtainedfromthegpiDCBtestsg1gc--31 层板胶合木层间断裂初度研宄■243ENF.层板胶合木层间断裂试验2.4.13ENF试验概况本文采用3ENF试验测定层板胶间II型断裂靭度C?,3ENF佘木层,ic试验试件的几2.15Z240mm30何形状如图所示,试件总长300mm,跨度2为,高度为rmTi,厚度B%30mm,预制裂纹总长(a+oo)为90mm,其中oo代表支座距试件边缘的距离,为30mm。预制裂缝位置处于层板胶结面处,先使用厚度为2mm的带锅预制88mm,再使用厚度为0.5mtn的刀片切割2nim,得到总长为90mrn的尖锐裂纹,裂纹的方向为顺纹(L)方向。P:zxpr_图2.15端部切口三点弯曲试验(3ENF)试验模型F-i.2.15Samleofthitbfl3ENFtgpreeponendendnotchedexure()tes3ENF试验采用的试验机仍然为济南时代试金试验机有限公司制造的WAW-300型微机控制电液伺服万能试验机(参见图2.8)。试验采用位移控制的加载方式,加载速度为2mm/min。为提升测量数据的精确度,分别附加了30kN的拉压力传感器和量程为±50mm的位移计进行荷载P和加载点线位移^的测量,并采用IMC数据采集器进行数-32- 大连理工大学硕士学位论文据的采集。在3ENF试验中,裂缝不会随着荷载的增加而持续扩展,其裂纹最多扩展至跨中位置,如果继续加载试件最终会发生受弯破坏3ENF。试验在目测试件跨中下部由于受弯而出现破坏时停止加载。1163ENF- ̄d图为试验典型的荷载加载点线位移、P、曲线,从该曲线可以看出,该3ENF试验中裂纹是不存在扩展的脆性断裂,荷载值尸在裂纹起裂后先骤降继而增大,并且在裂纹起裂后由于裂缝长度的改变,荷载值可能会超过起裂荷载,然后可能发生裂 ̄纹再次扩展的现象,如果继续加载试件最终会发生受弯破坏。本文在3ENF试验中将Pd曲线初始线性段的终点所对应的荷载定义为临界荷载Pc’进而计算层间I型断裂初度值一。由于试验所得试验曲线的形式是致的,并且由于需要在试验曲线上指明临界荷载一的定义,2173ENF,为了便于阅读因此试验给出了条具有代表性的试验曲线。图.为--eL试验曲线对应的荷载纵向应变)曲线。4500「4000-一-3500-U,--30005y71,2500-nc/!2000-/1500-/-/1000500-/0246810121416<mm5/一2_图.163ENF试验荷载加载点线位移(Pj)曲线-Fi.2.6Examllg1eoftheoadloadinlinedislacement{PScurvesobtanedome3ENFesspgp)ifrthtt4000r3500-3000.A52500-\/*2000-/-/1500-1000/5007L_—I_?__-Q■I012345678910-3VlO图2173ENF一-.试验荷载纵向应变(P曲线--Fi.2.17ExamleofeloadlonitudinalstrainPcurvesobtadth3ENFtgpthg(si)inefromeests-33- 层板胶合木层间断裂初度研究图2.18为3ENF试验试件破坏形式示意图,由于是剪切型破坏,裂纹上下表面仍然是闭合的,因此用黑笔标识出了裂纹扩展的长度,在荷载作用下初始裂纹基本沿着胶结,I型断裂破坏面处做自相似的扩展为层间I。由于层板胶合木制作时使用的为水性胶粘剂,并且胶层厚度较薄(O.hTim左右),胶层破坏没有明显现象。::。…十…w-V.t??-w-?-‘:fv■?<?(.:K‘圓■■隱一?,..、IJ?‘‘‘21ENF试件破坏形式图.83Fig.2.丨8Damageformof3ENFspecimen24.23ENF.试验结果_21ENF-图.9为3试验的临界荷载P?c与角度《2关系。从图2.19可以看出由于木1的材自身离散性较大引起^^,(,3495.01N:的分布较为分散的平均值和标准差分别为和60902N。-.在不考虑a2布的随机性引起不同角度差的试件数目不平均的情况下对i分尸^直线型趋势线PC-^-t-,《(^2^添加,可以发现iiC随着21的增大有轻微增大的趋势当21°°从0增大至44时趋势线上所对应的Piic值从3407.87N增大至3724.72N,约增大了9.30%。6000「_-_5000?■■■4。。。■"5一:■-■?■3000■■吕...:!..:^■??■.2000■-1000°01020304050?尸-2丨21NF试验尸图.93E《2.丨与临界荷载lie关系Fi.21Rltonshtweenenea-anderiiPobainedfromtheE1S1Fttg.9eaiipbethagliithctcalloadnct3ess852[基于经典粱理论层板胶合木II型层间断裂朝度G?c可由式(.3)H十算:-34- 大连理工大学硕士学位论文■'9七=(23).Iic巧图2.20为基于经典梁理论计算的3ENF试验G与《2.,G?平均值和标,ic1的关系c:的22准差分别为731.70J/m和264.40J/m。在不考虑《2-引起不同角度差的试1分布的随机性件数目不平均的情况下对G?C>42-1添加直线型趋势线,可以发现G?C随着《2-1的增大有°°一2定的增大趋势,0-044G.当C268461J/m1从增大至时趋势线上所对应的iic值从283224J,2156%23)?增大至./m约增大了.。并且从公式(.可以确定,C?c与所对应的.7^。1(:的平方值是呈线性关系的.1800r丨1500■■■■■一-?1200■_…?I■■■■-_■5900■■■-r0.丄rr+n\■■-_.0600■‘?■?■‘300■-°01020304050?尸2-丨图2.203ENF试验《2-1与断裂籾度Giic关系F.2.2Renan-hneig0latioshipbetweentheelajandfracturetoussGobtainedfromheFessgignct3ENtt2.5层板胶合木层间断裂3SENB试验2.5.13SENB试验概况..本文同时也采用了3SENB试验进行层板胶合木层间I型断裂行为研宄。为克服树一86[]木直径的限制,在木材的3SENB试验中单边切口试件般釆用由Wang提出的sandw,ich型试件本文也对工厂出产的层板胶合木进行了再加工和粘合,得到如图221.所示的sandwich型单边切口试件。3SENB试验试件的几何形状如图2.22所示,试件总长300mm,跨度2Z为240mm,高度2//为60inm,厚度5为30nim,预制裂纹总长a为27nim,裂缝高度比例U/2//)为0.45。预制裂缝位置处于层板胶结面处,先使用厚度为2mm旳带锯预制25mm,再使用厚度为0.5mm的刀片切割2mm,得到总长为27mm的尖锐裂纹,裂纹的方向为顺纹(L)方向。-35- .?层板胶合木层间断裂郁度研究;?.\、、^-?-r^T^Ai,-h.5r^.一,^4‘.‘、、L___?____i"慕I.」图2.21sandwich型单边切口试件--Fi.2.21Thesandwichteofaileedenotchedsecimengypsnggp层间胶结面j^ ̄_^^.LLIB\\i|小瀰mIt'■\LJ.图2.22单边切口三点弯曲试验(3SENB)试验模型F222Same--e-secmenig..plofthreeointbendintestofasinleednotchedi3SENBpgggp()3SENB试验釆用的试验机仍然为济南时代试金试验机有限公司制造的WAW-300?型微机控制电液伺服万能试验机(参见图2.8)。试验采用位移控制的加载方式,加载■?速度为Imm/m,OlN的拉压力传感器和量程in。为提升测量数据的精确度分别附加了lc为±50mm的位移计进行荷载P和加载点线位移S的测量,并采用IMC数据采集器进行数据的采集。3SENB试验中裂纹随着荷载的增加是持续扩展的,试验在荷载值P下降至0时停止。一2233SENB ̄)曲试验典型的荷载加载点线位移CPS线,在达到峰值荷载图.为 ̄^曲线近似于直线3SENB前P,当达到峰值荷载后起裂荷载P骤降,因此本文取试验一中的峰值荷载作为临界荷载P,,c。由于试验所得试验曲线的形式是致的并且由于需一,要在试验曲线上指明临界荷载的定义,为了便于阅读因此试验给出了条具有代表性--36 大连理工大学硕士学位论文的试验曲线.243SENB,。图2为试验试件破坏形式示意图.在荷载作用下初始裂纹基本.?沿着胶结面处做自相似的扩展,为层间I型断裂破坏。由于层板胶合木制作时使用的为0?水性胶粘剂.1mm,,并且胶层厚度较薄(左右)胶层破坏没有明显现象。.210r-……0:>18二八?150/\?./\(.1::/v初始非线性30-I—1——f——‘—.—■■0..00.51.01.520SImm2一图.233SENB试验荷载加载点线位移曲线F-mentcurvesi.2.23Examleofeloadloadinlinilainedfromthe3SENBtegpthgedspceobtasts‘"■‘4‘M#‘:-聽廳越蘇纖A海!儀:!纖編着,1___麵、图2,243SENB试件破坏形式‘Fi.2.24Damaeformof3SsecimenggENBp-2.5.23SENB试验结果2253SENBP-225图.为试验的临界荷载0(2。.ic与角度1的关系从图可以看出由于木材自身离散性较大引起/^〔的分布较为分散,Pic的平均值和标准差分别为351.28N和79.50N。在不考虑a2-i分布的随机性引起不同角度差的试件数目不平均的情况下对-/^c^2^添加直线型趋势线,可以发现P《-的趋势,当a2.ic随着21的增大有轻微减小i°°从0增大至42时趋势线上所对应的Pic值从357.08N减小至333.48N,约减小了6.61%。-37- .层板胶合木层‘间断裂籾度研宄;?-■.‘..,,700■「.、-,■.60--0■■500-...5,■■^■.°■.300.■■■:■■■■“■■200■■■100■01020304050图2.253SENBa2-P试验丨与临界荷载c关系Fig.2.25Relationshipbetweentheangela2-iandthecriticalload尸cobtainedfromthe3SENBtests26.本章小结对层板胶合木进行了层间I型断裂的双悬臂梁(DCB)试验和单边切口三点弯曲(3SENB)试验、层间II型断裂的端部切口三点弯曲(3ENF)试验,通过对全部试件截面的图片分析确定层板胶合木胶结面两侧年轮角度值,研究了不同年轮角度层板组合对于层板胶合木层间断裂临界荷载和断裂朝度的影响,得到如下结果和结论:(1)由于试验所使用的层板胶合木全部出自于工厂制作,胶结面两侧层板的年轮角度是随机出现的,通过对全部试件的截面图片进行分析,得到了不同试验方案试件的层板年轮角度a.。i和CC2以及其差值6(21的分布图(2)层板胶合木层间I型断裂的DCB试验所得临界荷载Pic随着层板年轮角度差一值的增大有定的增大趋势,基于经典梁理论计算所得的层板胶合木层间I型断裂朝度一G也有-引起不同角度差,C定的增大趋势。在不考虑层板年轮角度差C(2l分布的随机性°的试件数目不平均的情况下对添加直线型趋势线,发现随着《2-01从°447G340787N372472N增大至时趋势线上所对应的^0和ic值分别从.增大至.、22U.684.2420.6/m增大至832J/m,分别约增大了%、43%。(3)层板胶合木层间II型断裂的3ENF试验所得临界荷载P?c随着层板年轮角度差值《2-趋势,基于经典梁理论计算所得层板胶合木层间II型断裂1的增大有轻微的增大初度G?c也有轻微的增大趋势。在不考虑层板年轮角度差a2_的随机性引起不同角i分布-度差的试件数目不平均的情况下对G?Ca2-直线型趋势线,l添加发现随着-38- 大连理工大学硕士学位论文°°N增从0增大至57时趋势线上所对应的Phc和Gic值分别从448.99大至539.33N、2289.01J/m增大至127.53J/m,分别约增大了9.30%、21.56%。(4)层板胶合木层间I型断裂的3SENB试验所得临界荷载i^c随着层板年轮角度差值的增大有轻微的减小趋势。在不考虑层板年轮角度差《2-分布的随机性引起不同角1。度差的试件数目不平均的情况下对i^(^^2-^添加直线型趋势线,发现随着《2-1从0增°大至42时趋势线上所对应的Pic值从357.08N减小至333.48N,约减小了6.61%。-39- 层板胶合木层间断裂初度研究3.有限元分析一有限元法是种先进并且高效的数值分析方法,在当今科学技术发展和工程分析中获得了广泛的应用。并且,,伴随着电子计算机科学与技术的迅猛发展有限元法已经成87[]为计算机辅助设计(CAD)和计算机辅助制造(CAM)的重要组成部分。有限元法在力学分析研宄方面的功能十分强大,可以应用于弹性力学以及塑性力学。、静力学、动力学问题和固体力学、流体力学、热力学和磁力学等问题的研究中随着计算机技术的不断发展和提升,,多种多样的有限元法分析软件也不断涌现从早期的NASTRAN,AKSK,SAP到当前的ANSYS,MARC和ABAQUS都拥有众多的用户。上述有限元软件不仅功能强大,而且使用方便,因此适用范围广本文将采用ABAQUS6.12有限元分析软件对层板胶合木的上述层间断裂试验进行三维建模和分析与研宄。一一在有限元分析中,木材般被作为种正交各向异性的复合材料进行研宄。木材具一有粘弹性性质,在承受短期载荷作用时其材料特性般可以被当作线弹性的进行研宄,但实际上其力学性能随着加载速度、加载持续时间和木材含水率变化等因素是可能会变化的。并且,木材的线弹性假设只适合于小变形情况下的力学分析,甚至该假设在木材三个弹性主轴方向上也有一定的差异因此,对层板胶合木进行断裂力学研宄是相对困难的。目前将断裂力学应用于层板胶合木层间断裂行为的研宄较少,应用有限元法进行层板胶合木层间断裂軔度计算尚未出现,对其建立三维有限元模型难度较大,因此必须对木材的性质进行一定的简化。一1mm在对工程木材料进行有限元分析的过程中,由于其含有的胶层厚度般很小(0.左右),所以在大多数情况下可以不考虑胶层对其整体力学性能的影响。特别是脆性胶,89[]由于胶粘剂的性能与木材类似,所以大多可以忽略胶层的影响。本文研究的重点内容为层板胶合木胶结面商侧不同年轮角度层板组合对其层间断一,裂韧度的影响问题,因此忽略了木材的非均质性、粘弹性和非线性将其作为种均质并连续的正交各向异性材料进行弹性分析,同时也忽略了胶层的影响,重点关注胶结面两侧不同层板年轮角度的差值与层板胶合木层间断裂初度的关系。PG一R]ice基于全量塑性理论得到的个与裂纹尖端路径无关的积分值,即积分值。■/积分不仅具有守恒性,对于线弹性材料而言,■/积分和能量释放率G是等价的,是计算能量释放率的有效工具。/积分的计算通常采用等效积分区域法,,即通过散度定理一91[]将J积分中的回路积分转变为个包含裂纹尖端的区域积分。40 大连理工大学硕士学位论文=应力强度因子if是对应于裂纹尖端r0时的值,然而在采用数值方法计算Z时,=由于裂纹尖端应力的奇异性使得无法取得r0时的应力值,因此而无法完成K的计算。(3因此,参考公式.1),本文采用应力外推法来计算/i:,以1型断裂/^,的计算为例,其基本思路如下。=-Wmsa=尺jlnr100「,^A_(3.1)==LK,,Wm-Jl^rrr00“巧|I--|>0在裂纹尖端前缘直线上=0一>0一,即0时,对应于每个ri,有个非奇异的应力iT值以及对应的:ii(3.2)然后,构造数据对(n,Kn),并采用最小二乘法对已经构造成的数据对进行拟合。假定n和尺:li之间可以用线性关系来近似,则有K=Ar+B(\3.3)===r ̄当0时,K、ii:ir0;B。()因此,需在有限元程序中提取裂纹尖端前缘直线上若干单元积分点的主应力和坐,标值构造若干组(n,K“)的数据对,然后进行线性拟合并得到拟合直线的截距,即为所需的应力强度因子K”那么对应于临界荷载Pc时的应力强度因子即为断裂初度X,Kic。需要注意的是,在接近裂纹尖端时,(nu)的数据对可能会有较大的振9[1荡],应当通过剔除靠近裂纹尖端的数据来提高计算精度。本文釆用ABAQUS自带积分模块计算能量释放率G时,在裂纹尖端设置了退化的1/4节点奇异单元,其能使常规的20节点二次等参元(C3D20)的应力呈现的奇异性进而提高计算精度。本文在应用应力外推法计算应力强度因子ii:时采用ABAQUS的八节点常规单元(C3D8),由于裂纹尖端奇异性的存在,为得到较为可靠的结果,对所有有限元裂纹尖端进行了不同程度的网格细化。°°0-由于实测的层板年轮角度《,a丨和在50之间,因此在三维有限元模型中丨°°°。。。0120040>和a2取为、0、、3、、50并满足a2a〗,那么得到的角度差值41 层板胶合木层间断裂軔度研究°°°°°0-也为(T、10、20、30、40、50。-2根据角度差值c(12i与试验实测临界荷,平均值作为限元分析中的施加荷载P载Pc的关系取对应于附近Pc的。建立有限元模型时输入的材料参数取自本文第2章第1节试验测定所得的木材弹性常数。本文有限元分析部分主要包含以下内容:首先,采用ABAQUS有限元分析软件建立不同年轮角度层板组合的层板胶合木层间断裂试验有限元模型?,,利用ABAQUS自带模块计算对应于临界荷载Pc时的/积分进而得到层间断裂朝度Go分析不同年轮角度层板组合对其层间断裂顿度的影响Gc;然后,基于对不同年轮角度层板组合的层板胶合木层间断裂试验有限元模型对应于临界荷载Pc时裂纹尖端的应力场进行分析,提取裂纹尖端前部附近的若干单元积分点的主应力和坐标值,进而采用应力外推法计算得到相应的层板胶合木层间断裂初度Kc;最后,将基于基本梁理论的Gc计算结果与有限元的Gc计算结果进行对比分析,并将不同试验方法通过有限元所得同种断裂类型的断裂朝度Gc或iCc进行对比分析。3.1层板胶合木层间断裂DCB试验的有限元分析3.1DCBG.1试验能量释放率计算31ABADCB图.为本文采用QUS有限元程序进行试验能量释放率G计算时建立的三维模型。该模型在裂纹尖端设置奇异单元部分网格尺寸大约为L67mm,其余部分网"”格尺寸为4mm。ABAQUS有限元程序采用了回路积分法进行J积分的计算,而对于线弹性材料J积分与能量释放率G是相等的。(a)试件整体模型(b)裂纹尖端图3.1DCB试件有限元模型FTi!ii.3.1hefinteeletdeofheDCBtstedecmegmenmotespn°上述三维有限元模型共计有16548个单元。图3.2为和a30i同时为时的有限元模型变形图。-42- 大连理工大学硕士学位论文■(a)试件整体模型(b)裂纹尖端图°°、3.2年轮角度组合为30和30时的DCB试验有限元模型变形Fig.3.2ThemodeldeformationofDCBtestsintheFEcalculationsinthestateofanglescombinationis°°30and30采用有限元法计算不同年轮角度层板组合的DCB试验层间断裂初度Gic的施加荷312-33。载尸见表.,Gc.,Ga2见图.i计算结果见表3ic与i关系表3.1不同年轮角度组合时DCB试验有限元施加荷载PTab.3.1AppliedloadPinFEMwithdifferentangelcombinationsofDCBtestsP/N。ail0fo20^40^I0466447.67.80429.70571.96538.92.9150910466.91447.67509.80429.70571.962066.944.67509.8.70417042930466.9144767509.80.40466.91447.6750466.91表3.2不同年轮角度组合时DCB试验有限元计算所得断裂初度Tab.3.2G]ccalculatedbyFEMwithdifferentangelcombinationsofDCBtests^G2/J/mic()oar010W^4050I084.3079.06107.912..137^14286991087,7284.34113,9682.59146.412096.1592.22122.2186.8230104.8898.60127.7140109.73100.7850109.53-4-3 层板胶合木层间断裂軔度研究200「-150a"Sg一一^^^-__1001___MSgI50-00102030405060。?2-"3.3DCBa2_图有限元试验i与断裂朝度Gic关系Fi.3.3Relationshibetween?-andGiDCBttbFEM2obtanedfromtheessgp1icy直线型趋势线G《--对DCB试验有限元结果的G,C^21添加,可以发现ic随着21的°°一G93402.J/m增大有定的增大趋势,当《2050.1从增大至时趋势线上所对应的ic从225增大至1.14J/m,约增大了34%。3.1.2DCB试验应力强度因子/r计算图3.4为本文釆用ABAQUS有限元程序进行DCB试验应力强度因子计算时建立的三维模型。该模型在裂纹尖端进行了网格细化,裂纹前缘直线方向上网格细化部分长度为5mm,其网格尺寸为0.05mm,其余部分网格尺寸为2mm。.hiA孤ib(a)试件整体模型()裂纹尖端图3.4DCB试件有限元模型Fii.3.4ThefiniteelemenmodelofheDCBestedsecmengtttp°上述三维有限元模型共计有77870个单元。图3.5为cti和ai同时为30时的有限元模型变形图。-44- 大连理工大学硕士学位论文(a)试件整体模型(b)裂纹尖端°°图3.5年轮角度组合为30和30时的DCB试验有限元模型变形F..hldraiCBessillaoninheeofanlescoinaioiig35TemodeefomtonofDttntheFEcacutiststatgmbtns。°30and30采用有限元法计算不同年轮角度层板组合的DCB试验层间断裂朝度的施加荷°载尸见表3.1,当ai和同时为30时应力外推法的方案见图3.6,尺ic计算结果见3.3iST37。表,c.i与丨关系见图0.4r=.+.y02001x02122?0=.3Rf0.9867i-f0.2^?0.1°1.00.20.40.60.8r/mm°。图3,6年轮角度组合为30和30时的DCB试验应力外推法Fi.3TlBteEationg.6hestressextrapoationofDCstsintheFcalculationsinthestateofanglescombinis。。30and30表3.3不同年轮角度组合时DCB试验有限元计算所得断裂朝度Tab.3.3KccalculatedbyFEMwithdifferentanelcombinationsofDCBtests\gKlMa^]c(Jmp)air0io20304050I ̄“ ̄00.2140.2060.2370.2000,2490.226100.2140.2060.2360.1980.246-45- 层板胶合木层间断裂軔度研宄 ̄"““ar0W:MM4050tI ̄200.2140.205033,20^300.228.2120.2010400.2000.192500.1850.30「?-‘.?B^0.25^”?=?n,n■D0?r.201。。.^0.15^-40.10■0.05‘ ̄‘ ̄‘―‘ ̄‘ ̄‘ ̄‘ ̄ ̄■_I■ ̄I0000102030405060。?2-"37DCB试验a2-/T图.有限兀与断裂韧度cii关系F..li-ig37Reatonshibetweena2andKobainedfromtheDCBtestbFEMpi\ctsy"C-《-对DCB试验有限元结果的尺iCt2l添加直线型趋势线,可以发现ric随着21的°°一增大有定的增大趋势.050,当《21从增大至时趋势线上所对应的尺IC值从"0.205Ma>0237M,51p/iir增大至.pa>/iii约增大了1.6%。3.2层板胶合木层间断裂3ENF试验的有限元分析3.213ENF<?.试验能暈释放率计算.图3.8为本文采用ABAQUS有限元程序进行3ENF试验能量释放率G计算时建立的三维模型。在建立3ENF试验模型过程中,在裂纹面设置了接触对以避免裂纹表面相互侵入,,法向的接触属性硬接触切向的接触属性为无摩擦。该模型在裂纹尖端设置奇异单元部分网格尺寸大约为1,.25mm其余部分网格尺寸为4tnm。46 -大连理工大学硕士学位^论文-圖.(a)试件整体模型(b)裂纹尖’'ENF试件有限元模型?图3.83Fi.3.8Thefiniteelementmodelofthe3ENFtestedsecimengp°上述三维有限元模型共计有9048个单元。图3.9为和同时为30时的有限元模型变形图。隱.(a)试件整体模型(b)裂纹尖端°°图3,9年轮角度组合为30和30时的3ENF试验有限元模型变形Fi.3.9ThemodeldeformationofNFesitFEcalculaonseaeanlescombinaonsg3Ettsnhetiinthsttofgtii°°30and30采用有限元法计算不同年轮角度层板组合的3ENF试验层间断裂初度G?c的施加荷3-.4,G?计算结果见表3.5,G?c023.10。载见表c与1关系见图表3.4不同年轮角度组合时3ENF试验有限元施加荷载尸Tab.3.4AliedloadPinFEMwithdifferentanelcombinationsof3ENFtestsppg01020^40^^^03411..393400.093552573675.63317636.103411.393400.093552.573675.63317636203411.393400.093552.573675.6330411.393400.093552573.403411.393400.09503411.39-47- 层板胶合木层间断裂軔度研究■^°°°a3ENF试验中实测所得最大值为44,因此,对于角度组合为0和50的情况,由于缺少临界荷载A:而并未进行有限元计算表3.5不同年轮角度组合时3ENF试验有限元计算所得断裂朝度Tab3.5GiiccalculatedbyFEMwithdifferentangelcombinationsof3ENFtestsG2/Jmic(/) ̄"“。\2/0fO20^4^50^^0820—.50814.98890.3955.54716.82110820.88816.60893.66959.87716.82.89520820.88817.64.95960.10.30820.88818.69894.9140、^823.18819.45..5082702..1200r■0^900一a0今°、600■。aO300!—,—.I—.—I—.—I—‘—I——■>Q—01020304050。?丨2-丨.’“‘图3.10有限元3ENF试验a2-与断裂朝度Gic关系iiFiionsh-i.3.10Relatibetween0:2andGobtanedfromthe3ENFtestsbFEMgp1hcy对3ENF试验有限元结果的添加直线型趋势线,可以发现Gc随,着的°°2增大有轻微的增大趋势,当从0增大至40时趋势线上所对应C?nc从837.30J/m28573J增大至.5/m,约增大了2.41%。3.2.23ENF试验应力强度因子/T计算图3.11为本文采用ABAQUS有限元程序进行3ENF试验应力强度因子计算时建立的三维模型。同样的,在建立3ENF试验模型过程中,在裂纹面设置了接触对以避免裂纹表面相互侵入,法向的接触属性为硬接触,切向的接触属性为无摩擦。该模型在裂纹-48- 大连理工大学硕士学位论文尖端进行了网格细化,裂纹前缘直线方向上网格细化部分长度为4mm,其网格尺寸为0.02mm,其余部分网格尺寸为2mm。.":^■5變853:1(a)试件整体模型(b)裂纹尖端图3.113ENF试件有限元模型Fig.3.11Thefiniteelementmodelofthe3ENFtestedspecimen上述三维有限元模型共计有125640个单元。图3.12为和ai同时为3{r时的有限元模型变形图。(a)试件整体模型(b)裂纹尖端°°图3.12年轮角度组合为30和30时的3ENF试验有限元模型变形Fig.3.12Themodeldeformationof3ENFtestsintheFEcalculationsinthestateofanglescombinationis。。30and30采用有限元法计算不同年轮角度层板组合的3ENF试验层间断裂初度的施加荷°载P见表3.4,当ai和a_i同时为30时应力外推法的方案见图3.13,/Gic计算结果见3-表.6,iC《3.14。丨丨C与2丨关系见图"对3ENF试验有限元结果的;^?-加直线型趋势线,可以发现兄着《2-,,£^21添,c随1的°°一增大有定的增大趋势,当从0增大至50时趋势线上所对应的ii:?c值从0585.66316.Mpa^/^^增大至02MpaN/i^,约增大了1.%。-49- 层板胶合木层间断裂軔度研究.1.5「=4--y2.0776x0.515_12-IR=0.980^09【0-:.6■^0-.3'‘‘‘‘o.o0.00.10.20.30.4r/mm。。图3.13年轮角度组合为30和30时的3ENF试验应力外推法Fi.314Tsressexraolaionof3ENFessinheFEcaculationsnhesaeofanlescombinaiong.hettpttttlitttgt。。is30and30表3.6不同年轮角度组合时3ENF试验有限元计算所得断裂韧度Tab.3.6ATcalculaedbFEMwihdifferenanelcombinationsof3ENFtessictyttgt^IIc/(Mpa^/m)0air10^304050I00—.6480.6770.7000.603100.6420.6400.6930,6890.593200.6180.6150.6400.659300.5750.5710.593400.5140.509500.4351.0「一0-奶.8。曰80^“"一0.6一=a!!。S。“0-%.40-.2——————‘‘■—‘■—‘■—‘—■‘—■‘0.001020304050。?丨2-1图3.14有限元3ENF试验《2-1与断裂朝度/^uc关系Fi.3.14Relatonhieweena-andKENFispbtzicobtainedfromthe3testsbFEMgiiy-50- 大连理工大学硕士学位论文33SENB.3层板胶合木层间断裂试验的有限元分析3.3.13SENB试验能量释放率G计算图3.15为本文采用ABAQUS有限元程序进行3SENB试验能量释放率G计算时建立的三维模型。该模型在裂纹尖端设置奇异单元部分网格尺寸大约为2nim,其余部分网格尺寸为4mm。舰七瓶瓜(a)试件整体模型(b)裂纹尖图3.153SENB试件有限元模型F.3.15fiitltl3SENBedsecenigTheneeemenmodeofthetstepim°上述三维有限元模型共计有18228个单元。图3.16为0C1和《2同时为30时的有限元模型变形图。(a)试件整体模型(b)裂纹尖端°。图3.16年轮角度组合为30和30时的3SENB试验有限元模型变形F.defoltiithtofliig.316Themoldermationof3SENBtestsintheFEcalcuaonsnestaeangescombnation°°s30an30id采用有限元法计算不同年轮角度层板组合的3SENB试验层间断裂朝度Gic的施加荷载P见表3.7,Gic计算结果见表3.8,Gic与a2.i关系见图3.17。-5-1 层板胶合木层间断裂韧度研宄表3.7不同年轮角度组合时3SENB试验有限元施加荷载PTab.3.7AliedPFhdiffelcipploadinEMwitrentaneombinatonsof3SENBessgtt。ail01020304050I036146.二.33623327.96360.19399.661036146336.23327.963609399.66..12036.46336.123327.96360.19306.3146336.23327.9640361.46336.2350361.46°°a°3SENB试验中实测所得a2-i最大值为42,因此,对于角度组合为0和50的情况,由于缺少临界荷载Pc而并未进行有限元计算?表3.8不同年轮角度组合时3SENB试验有限元计算所得断裂軔度Tab3.8GcalculaedbFEMwithdiffenelcombinaonsofictyrentagti3SENBtestsG2/J/mic() ̄0air10^304^50^=0113.73111.16123.59166.05209.6910157.84153.09150.47.1827223.3220189.18.7172.177308205.2630224.80202.9491.51140248.08214.9750253.48300r250-B"°°200-?i?j。。.°。i150°°O-10050■QI■01020304050。?2-"图3.17有限元3SENB试验a2-i与断裂軔度G系ic关F.3.17iiigRelatonshbetweenaz-andGobainedfromthe3SENBtessFEMpjicttby-52- 大连理工大学硕士学位论文-a-,可以发现G对DCB试验有限元结果的G,c2i添加直线型趋势线,c随着的增大是近似保持恒定的。3.3.23SENB试验应力强度因子/T计算图3.18为本文采用ABAQUS有限元程序进行3SENB试验应力强度因子计算时建立的三维模型。该模型在裂纹尖端进行了网格细化,裂纹前缘直线方向上网格细化部分长度为3mm,其网格尺寸为0.05rnrn,在试件高度方向上网格细化部分以上部分的网格尺寸为,其余部分网格尺寸为2rnm。Imm(a)试件整体模型(b)裂纹尖端图3.183SENB试件有限元模型Fi3.hefinielemodelfheedecimen.18Ttementot3SENBessgttp°上述三维有限元模型共计有102960个单元。图3.19为ai和ai同时为30时的有限元模型变形图。(a)试件整体模型(b)裂纹尖端°°图3.19年轮角度组合为30和30时的3SENB试验有限元模型变形Fi.3.19Themodeldeformationof3SENBteststheEcitettefancotginFalculatonsinhsaoglesmbinaion。°is30and30--53 层板胶合木层间断裂钥度研究采用有限元法计算不同年轮角度层板组合的3SENB试验层间断裂初度Z,c的施加°P见3.7.荷载表,当和ai同时为30时应力外推法的方案见图320,ATic计算结果3.9,iCa-321见表.。,c与2i关系见图1.5r=1.7745X+0.3657■y_1.222=R0.962.。..9#?6■0-.31‘———?——‘‘——‘—0■‘‘‘.00.00.10.20.30.40.5r/mm°°图3.20年轮角度组合为30和30时的3SENB试验应力外推法Fig.3.20Thestressextrapolationof3SENBtestsintheFEcalculationsinthestateofanlescombinationg°=^and30is30表3.9不同年轮角度组合时3SENB试验有限元计算所得断裂朝度Tab.3.9K\ccalculatedbyFEMwithdifferentangelcombinationsof3SENBtestsAr/Maiic(p\/n)ar010^304050tI°a/i=00.2470.2360.2480.2940.346100.2670.2670.2800.3230.367200.3160.3140.3210.358300.3660.3560.353400.4050.384500425.3S对ENB试验有限,可以发a2-元结果的添加直线型趋势线现随着i的增大是近似保持恒定。-54- 大连理工大学硕士学位论文■0.6「广■0.50-°0.40。。。。令三。“0-^.3I°0aa0°°00■.2?0.10'—‘—‘—‘—‘—‘—■—‘—■—‘—‘—‘001020304050。?丨2.1图3.21有限元3SENB试验ct2-i与断裂初度/Tic关系32libifrbFi..1Reatonshietweena2-and尺obtanedomthe3SENBtestsFEMgpiicy.3.4试验结果与有限元结果对比341DCB..试验结果与有限兀结果对比图3.22为DCB试验基于经典梁理论(EBT)的Gic计算结果与有限元法(FEM)的G,c计算结果的对比。基于经典梁理论和采用有限元法的G,c计算结果平均值分别为22一10435J033J./m和11./m,并且两种方法所得结果的直线型趋势线十分致。300「-Td■EBTlifEBTrenneo-250□---FEMreneofFTdlinEM200-C■“^■?i■^150-■■B“■■■-50-广-■■■0I_1___I_I__I_I__I_1_I_1__I_I_I0102030405060。?2-"图3.22DCB试验经典梁理论Gic结果与有限元Gic结果对比Fig.3.22ComparisonofthefracturetoughnessGicobtainedfromtheDCBtestsbyEBTandFEM3423ENF..试验结果与有限元结果对比图3.23为3ENF试验基于经典梁理论(EBT)的计算结果与有限元法(FEM)的Gi,c计算结果的对比。基于经典梁理论和采用有限元法的G?c计算结果平均值分别为-5-5 层板胶合木层间断裂靭度研宄22731.70J/m和845.25J/m由于裂缝制作的误差很能引起I/II复合型断裂的出现,因此导致了试验基于经典梁理论的G?c计算结果略微偏小,但是两种方法所得结果的趋势线近似一致。2400「?TrendiineofEBT■丁20。0-口FEM---TrendlFEMineof-_^1600_*■■^1200-■■:■■■■U“,?口“■D■。8。。?"■■■:■■400■■■_■IIIII1II0II01020304050?尸2-丨图3.233ENF试验经典梁理论Gnc结果与有限元Gic结果对比Fi.3.23ComarisonofthefracureouhnessGobainedfromthe3ENFtestsbEBTandFEMgpttgncty3.4.3DCB有限元与3SENB有限元结果对比图3.24为DCB试验和3SENB试验采用有限元法计算所得Gic结果的对比。从图3.24可以发现,3SENB试验与DCB试验相比有限元法计算Gic的结果离散性较大。并,3SENB试验结果的直线型趋势线随着角度差《2-DCB试验结果且1的增大保持稳定,一的直线型趋势线随着随着角度差a2-i的增大有定的增大趋势。350“.—Trend3SENBiineof3SENB[-300□DCB—TrendiineofDCB250-■■5:“200-._o__.__:;■“O50-■-1a°-S一_100--§十850■°0102030405060。?I2-l图3.24DCB试验与3SENB试验有限元Gic结果对比F.3.24ComarisonofheactureouhnessGobtainedfromtheDCBand3SENBtessFEMigptfrtgictby-56- 大连理工大学硕士学位论文...^3.25DCB3SENB尺。图为试验和试验采用有限元法计算所得ic结果的对比从图3.25可以发现,3SENB试验与DCB试验相比有限元法计算尺,c的结果离散性较大。并且,3SENB试验结果的直线型趋势线随着角度差的增大保持稳定,DCB试验结果一的直线型趋势线随着随着角度差a2-i的增大有定的增大趋势。0.6r.—Trend3SENBlineof3SENB-—Td0linefDCB.5。DCBreno.-4::.V.:0.3:::;身:■&5.n-O-0u.2&0-lt^?0.1^—————————‘‘‘■I_‘II—■‘—■‘000261003040500。?2-"图3.25DCB试验与3SENB试验有限元Kic结果对比Fig.3.25ComparisonofthefracturetoughnessK丨cobtainedfromtheDCBand3SENBtestsbyFEMD2CB试验和3SENB试验采用有限元法的Gic计算结果平均值分别为110.33_l/m和"2174/m..84J,^^:果平均值分别为0215Ma>fi0.3243SENB试验结10计算结p/ii^果较DCB试验结果偏大,这可能是由于木材性质的离散性、试件的形式与尺寸差别、裂纹制作误差以及采用有限元应力外推法进行计算应力强度因子的误差等原因造成的。3.5本章小结采用ABAQUS程序对层板胶合木层间I型断裂的双悬臂梁(DCB)试验和单边切口三点弯曲(3SENB)试验、层间II型断裂的端部切口三点弯曲(3ENF)试验进行了三维数值模拟,通过定义不同材料的方向得到不同年轮角度层板组合情况的三维有限元模型,结合临界荷载值分别采用回路积分法和应力外推法计算得到层间断裂朝度Gc和Kc.研究了不同年轮角度层板组合对于层板胶合木层间断裂行为断裂初度的影响,并将有限元计算结果与基于经典梁理论计算结果进行了对比,以及将不同试验方式G,c和c,可以得到如下结果和结论:,的有限元计算结果进行对比(1)层板胶合木层间I型断裂的DCB试验有限元结果表明,层间断裂初度随°°着年轮角度差值的增大有轻微的增大趋势,当从0增大至50时趋势线上所对应M:M^的冗0.205a^i^i0.237a>,约增大了15.61%。1(:值从p/增大至py-57- 层板胶合木层间断裂初度研宄(2)层板胶合木层间II型断裂的3ENF试验有限结果表明,层间断裂軔度尺随°°着年轮角度差值的增大有轻微的增大趋势,当从0增大至50时趋势线上所对应的i^0.58aiii0.662约3.16%。iic值从5Mp>/增大至增大了1(3)层板胶合木层间I型断裂的3SENB试验有限结果表明,层间ff裂勒度Gic和尺IC随着年轮角度差值的增大近似保持恒定。(4)层板胶合木层间I型断裂的DCB试验基于经典梁理论与有限元法的Gic计算一致,基于经典梁理论和釆用有限元法的(7平均值结果的直线型趋势线十分,C计算结果22分别为104.35J/m和110.33J/m,两者相差很小,有限元计算结果验证了试验结果的准确性。(5)层板胶合木层间II型断裂的3ENF试验基于经典梁理论与有限元法的(?,,c计一算结果的直线型趋势线近似,有限元计算结果验证了试验结果的准确性致。但是,基22于经典梁理论和采用有限元法的CJiic计算结果平均值分别为731.70J/m和845.25J/m,试验基于经典梁理论的G?c计算结果略微偏小,这可能是由于裂缝制作的误差很能引起I/II复合型断裂的出现引起的。(5)层板胶合木层间I型断裂的DCB试验和3SENB试验采用有限元法计算所得的Gic和r,c结果进行对比,发现3SENB试验与DCB试验相比Gic和Kic的结果离散性较大。并且,3SENB试验有限元结果的直线型趋势线随着角度差的增大保持稳一定.,DCB试验有限元结果的直线型趋势线随着随着角度差《21的增大有定的增大趋2。DC3SENBCJ11033J势B试验和试验采用有限元法的./mic计算结果平均值分别为和"2174.84J/m,^:1(:计算结果平均值分别为0.215Mpa^/iii和0.3243SENB试验结果较DCB试验结果偏大,这可能是由于木材性质的离散性、试件的形式与尺寸差别、裂纹制作误差、提取临界荷载的误差以及采用有限元应力外推法进行计算应力强度因子的误差等因素填成的。58 大连理工大学硕士学位论文4结论与展望41.结论本文主要进行了不同年轮角度层板组合层板胶合木的层间I型断裂的双悬臂梁(DCB)试验和单边切口三点弯曲(3SENB)试验,以及层间II型断裂的端部切口三点弯曲(3ENF)试验,并釆用有限元法对试验进行了有限元模拟,研究了层板胶合木的层间断裂軔度与胶结面两侧层板年轮角度的关系。本文的主要结论如下:(1)层板胶合木层间I型断裂的双悬臂梁(DCB)试验基于经典梁理论计算所得其层间断裂勒度G一,c随着层板年轮角度差值的增大有定的增大趋势。有限元计算结果‘与试验结果吻合较好。(2)层板胶合木层间II型断裂的端部切口三点弯曲(3ENF)试验基于经典梁理论计算所得其层间断裂朝度(??0随着层板年轮角度差值的増大有轻微的增大趋势。有限元计算结果与试验所得Guc结果相比略微偏小,可能是由裂纹制作误差引起的。但是,两种方法所得结果的趋势线近似一致。(3)层板胶合木层间I型断裂的双悬臂梁(DCB)试验的有限计算所得其层间断裂軔度ICiC随着年轮角度差值的增大有轻微的增大趋势。(4)层板胶合木层间II型断裂的端部切口三点弯曲(3ENF)试验的有限元计算所得其层间断裂朝度iCnc随着年轮角度差值的增大有轻微的增大趋势。(5)层板胶合木层间I型断裂的单边切口三点弯曲(3SENB)试验的有限元计算所得其层间断裂初度Gi^ic随着年轮角度差值的增大近似保持恒定。ic和(6)层板胶合木层间B)试验和单边切口三点弯曲I塑断裂的双悬臂梁(DC(3SENB)试验采用有限元法计算所得Gic和尺IC结果的对比发现,单边切口三点弯曲(3SENB)试验结果较双悬臂梁(DCB)试验结果偏大,这可能是由于木材性质的离散性、试件的形式与尺寸差别、裂纹制作误差、提取临界荷载的误差以及采用有限元应力外推法进行计算应力强度因子的误差等因素造成的。一批次由工厂生产的用于实际工程的层板胶合木本文的试验材料全部采用同,因此年轮角度的分布具有一定的局限性,并没有得到覆盖全部年轮角度差值的试件。但是,本文对于工程实际使用的层板胶合木中出现概率较大的年轮角度情况的试件进行了层间断裂韧度测定试验,得到的试验结果离散性较大,只进行了相关定性的趋势分析,但是本文为层板胶合木的断裂研宄提供了一个新的方向一,通过试验和有限元法得到了定一数目的断裂韧度值,为进步开展该问题的研究奠定了基础。-59- 层板胶合木层间断裂初度研宄4.2展望应用线弹性断裂力学研宄木材顺纹断裂行为的成果相对较多,这在传统的材料力学之外为评价木材的强度问题提供了一种新的并且更加有效的方法一。层板胶合木作为种集成材,其中存在大量的胶结面,充分考虑胶结面两侧木材年轮角度不同对其层间断裂軔度影响的相关研宄较为缺乏。本文对层板胶合木的I型层间断裂和n型层间断裂分别进行了比较全面的试验研宄和有限元模拟,填补了该方面研宄的空白。虽然本文取得了一—一定的成果,些问题与不足:但仍然存在,需要进步的进行探索与研究(1)层板胶合木层间断裂朝度不仅与胶结面两侧层板的性质有关,也与层板胶合木制作时所采用的胶合剂种类、胶合剂用量、胶结加压压强大小与持续时间和胶合剂添一加剂(如纤维素)等原因相关,可以进行进步的研究;(2)由于本文采用的层板胶合木出自工厂制作,胶结面两侧层板年轮角度是随机出现的,没有得到角度组合分布均句且广泛的的目标试件,并且裂纹的预制是通过切割而成误差较大,后续研宄可对工厂制作层板胶合木的层板组合情况和预留裂纹等方面提一出定的要求;,便于对层板胶合木进行较为准确的层间断裂行为研究(3)本文在进行两种层板胶合木层间I型断裂试验的断裂窃度有限元计算结果对一一比时,发现存在定的差异,对于产生该差异的具体原因需要进步研究。60 大连理工大学硕士学位论文参考文献1]潘景龙祝恩淳.木结构设计原理[M].北京:中国建筑工业出版社2009.[,,一2-周海宾.中国木结构产业的现状及发展J.木材工业2012261:710.,[]()[],()[3]樊承谋,王永维,潘景龙.木结构[M].北京:高等教育出版社,2009.[4]谢力生.木结构材料与设计基础[M].北京:科学出版社,2013.5-[]中华人民共和国国家标准.GB/T507082012胶合木结构技术规范[S].北京:中国建筑工业出版社,2012.6-陆继圣.日本胶合木的制造与应用[J].木材工业9966:3031.[],1,()7樊承谋.北:[],张盛东陈松来等.木结构基本原理[M]京中国建筑工业出版社2008.,,—New8UNECE/FAO.Forestroductsannualmarketreview20072008R.YorkandGeneva:[],[]pGenevaTimberandForestStudyPaper,23,2008.9J-李登华杨学兵:2011:6266.[]王永维.胶合木结构特点评述[.36,,]四川建筑科学研究,(5)[10]SmithI,LandisE,GongM.Fractureandfetigueinwood[M].Hoboken,NewJersey:JohnWileyandSons,2003.-111.邵卓平..200339119125[],任海青,江泽慧木材横纹理断裂及强度准则[J]林业科学,,():12J2014周海宾..3:[],孙华林,张俊珍等结构用胶合木性能若干影响因素[]建筑材料学报,,()553-558.13.:.许金泉界面力学[M].北京科学出版社,2006[] ̄141郭营洁.J,世20124:4348.[],任海青结构用胶合木生产工艺研究进展[]界林业研究,,-15薛世峰侯密山.工程断裂力学[M].东营:中国石油大学出版社2012:11.[],,-16志东.北京:201211.郦正能,张纪奎等.应用断裂力学M北京航空航天大学出版社:[],关[],h.henomenadflowinlidhilohicalf17GrifitAATheofrutureansosJ.Psotransactionsothe[】pp[]proalsocietyofLondon.SeriesAcontaininaersofamathematicalorhysicalcharacterA1921y,g,,ppp-221:163198.18GriffithAA.ThetheorofrutureJ.ProceedinsoftheFirstConressofAliedMechanics[y]g,]p[gpp4-192:5563.-19Irw.inGRFracturednamicsJ.Fracturinofmetals1948:147196.[]y[]g,[20]OrowaxE.Fractureandstrengthofsolids[J].ReportsonProgressinPhysics,1949,12;185.21IrwinGR.AnalsisofstressandstrainsneartheendofacracktraversinlateJ.Alied[]ygap[]pp-Mechanics195724:109114.,,22GustafssonP.fracturemechanicsmodelsforstrenthanalsisoftimberbeamswithaholeoranotch:[]gya-reortofRILEMTC133trrandiniilcrrreort:meansessappoachtaackappoachaperinJ.Volp(pp)[]Reo-rtTVS7134LundUniveritn1pMsy993.,Swede,Muskhe.23lishviliNL数学弹性力学的几个基本问题赵惠元.M.:1958.[],译[]北京科学出版社,[24]WilliamsML.Onthestressdistributionatthebaseofastationarycrack[J],ASME,JournalofA-ppliedMechanics1957:109114.,-6-1 层板胶合木层间断裂初度研究[25】BuecknerH,BuecknerH.AnovelprincipleforthecomputationofstressintensityfactorsJ.Z.[]n-AewMa.Mec.1:52546.g.thh9709,[26]WuXR,CarissonJ.Weightfunctionsandstressintensityfactorsolutions[M].Pergamon,1991.一-27赵伟学仁.J空材料学报122.航990:9.[],吴种高效准确的三维权函数法及其应用[],[28]SneddonIN,LowengrubM.Crackproblemsintheclassicaltheoryofelasticity[J].JohnWiley8cs--SonInc.NewYorkLondonSdney1969.,,y,ErF-29doanGutaG.ThetressanalifmultilaeredcomositeswithaflawJ.International,ssso[]gpyyp[]Journa-lofSoandSrres117:3961lidstuctu97.,,[30]ErdoganF,GuptaGD,CookTS.Numericalsolutionofsingularintegralequations[M]//Methodsofanandsounckner-lsNeerands1973:36425.yisaltiosofcraproblems.Sprigthl8,[31]WilsonWK.Finiteelementmethodsforelasticbodiescontainingcracks[M]//Methodsofanalysis-andouonofcrackr.ringetiernd1973:484515sltispoblemsSprNeilas.,32HenshellRDShawKG.CracktifiniteelementsareunnecessarJIttillfo[.nernaonaJournar],py[]Me-umercatEi1975345507Nilhodsinngneering9;9.,,()33BarsoumRSOntheiticfiniteelementsinlinearfractremechaicsJ.[.useofsoaramerun]p[]Me-InternationalJournalforNumericalthodsinEngineering1976101:2537.,,()[34]王克仁,徐纪林,高桦.CalculationofStressIntensityFactorsForCombinedModeBendec-SpimensJ.FractureandSociety19785:123133.[],,()35Hres-nSihGC.andbookofstsitensityfactorsM.LehighUniversityInstituteofFractureandSolid[][],Mechanics1973.,36RookeDP,CartwrightDJ.Compendiumofstressintensityfactors[J.ProcurementExecutive,[]]MniDe.H.MS.1B.1istryfence.O.976330ook976.of,p(),37TadaH,PraisPC,IrwinGT.TheStressIntensityFactorhandbookM.Hellertwon:DelReasearch[][]Corporation,1985.38鲍甫成.中国主要人工林树种木材性质北京:中国林业出版社1998.],[J-39邵卓平.木材断裂初性测试研宄.力学与实践2002244:4951.[][],,()[40]WuEM.ApplicationoffractureMechanicstoanisotropicplates[J].JournalofApplied344-Mechanics.1967967974.,():-41孙艳玲振友.木材断裂力学的研究J.1997193:8592.[],鹿[】北京林业大学学报,,()42K-PingMJSutherlandIJLenocL.FracturetouhnessofwetanddrinusradiateJ.Holzals[],,ggy[]-rkof-RohundWest1999,574:235240.,()43邵卓平.:09木材与竹材的断裂与损伤[D].安徽安徽农业大学.20.[]44P-lanasJGuineaViMeneraledszeeffecteuationforuaibritlialMaial[GElces.GiisteMaters.ters],,qqes-andSructur1997205:671.,,()687[45]BoatrightSWJ,GarrrenttGG.Theeffectofmicrostructureandstressstateonfracturebehaviourof-wood..JournalofMaterialsSciencd198318;21812199,,46S-TsceSETanDMTeEKNmewooduretanzlhsch.ewsplittinthodforfract[]gg,,gggcharace-trizationTec1929:50..WoodScihnol9531;62 大连理工大学硕士学位论文[47]DEMouraM,DouradoN,MoraisJJL.Crackequivalentbasedmethodappliedtowoodfracture-characterizationusingthesingleedgenotchedthreepointbendingtestJ.EnineerinFracture[]gg-Mechanics201077310520:5.,,()48王坚伟.D.:.2008木材断裂行为及其声发射特性研究[江苏南京林业大学.[]]49YoshiharaHSatohA.Shearandcratideformationcorrectionforthedoubletilbeam[],ckpcanever--andthreelIIIfrpointendnotchedfexurespecimensformodeandmodeacturetoughnessmeasuremenwoodJEnneernFr-tMeccs200976:335346of,giigacturehani3.[],,()50邵卓平任海青J.200372:[],江泽慧.柔度法标定木材断裂軔性的研究[林业科学1,],,()-1.11216[51]YoshiharaH,KawamuraT.ModeIfracturetoughnessestimationofwoodbyDCBtest[J].Cosarndnufacn-omA:ednceur200611:212113itesPtaliscieamati3705.pppg,,()52YoshiharaH.ModeIandmodeIIinitiationfracturetoughnessandresistancecurveofmedium[]dens--ityfiberboardmeasuredboublecantileverbeamandthreeintbendendnotchedflexuretestsydpo-JEnnFechan..gieeringractureMics20107713;25372549[],,()53DEMouraMMoraisJJLDouradoN.AnewdatareductionschemeformodeIwoodfracture[],,characterizationusinthedoublecantileverbeamtestJ.EnineerinFractureMechanics2008g[]gg,,7513-:38523865.()54MeBEJi-lmenorlSouchaudSmitbuhletal.?curvebehairandrouhnessdeveotoffracture[],,ch,vogpsursJnernaonaJournaofFure2002-face,Ittillract1144:307325.[],,()[55]MorelS,MourotG,SchmitbuhlJ.Influenceofthespecimengeometryon及-curvebehaviorandacceIJourna21--42rougheningoffrturesurfesJ.nternationallofFracture,2003,112:23.[]()56YoshiharaH.Influenceofsan/dethratioonthemeasurementofmodeIIfracturetouhnessof[]ppg-W-woodlJ47bendnotchedfexuretest.JournalofoodScience11:812.y,200,[]()57Yosh.MIIhnbinedbENFiharaHodeinitiationfracturetougessanalsisforwoodota3test[J,[]yy]Cs-omoitesscienceandechnolo20056514:212207tgy98.p,,()[58]SilvaMAL,DEMouraM,MoraisJJL.NumericalanalysisoftheENFtestformodeIIwoodfrac-tureJ.CompositesPartA:AliedScienceandManufacturin,2006,379:13341344.[]ppg()[59]SilvaMAL,MoraisJJL,DEMouraM,etal.ModeHwoodfracturecharacterizationusinthegELSesJnM-ttirinFracture200774141247.,Egneegechanics:2331[],,()nabnof-60WataeKShidaSOhtaM.Evaluatio?ndcheckoationbasedonmodeIfracture[】,,prpgaJlSc2013-376touhnessof.sui(Crtomeriaaonica.JournaofWoodience1575:71ggyp,,jp)[]()61ArdalanMDeamBFraiacomoM.Exerimentalresultsoffractureenerandfracture[]y,,gpgytouhnessofRadiataPinelaminatedveneerlumber(LVL)inmodeI(oenmg)[J].Materialsandgprucure201245-Stts8:11891205.,,()[62]ReitererA,TscheggS.TheinfluenceofmoisturecontentonthemodeIfracturebehaviourofrucewJ-spood.JournalofMaterialsScience20023720:44874491.[],,()rerer-woodna63SchachneAchertrocouss.HRiteStanzlTsSE.OhoticfraturetofJJourlof[],,ggpghne[]e-rialsScienceLeters20001920Mat:17831785.,,()63 层板胶合木层间断裂軔度研究64ReererABureIeaTrercemeititSinnGtl.heradialinfontofthewoodstructureandits[],g,,me" ̄imlicationonmechancaandfracicririsreepilturechanalpopertesacompaonbetweentwot-seciesJ.JournalofMtl2002535940.p[]aeriasScience,,37():965R-eitererASinnGStanzlTscheSEFracurecharacrisicdiff.ttetsoferentwoodspeciesunder[],,ggerend.A2002modeIloadiniculartotherainJMaterialsScienceandEnineerin:21:gppg[]gg,33,()29-36.-66StanzlTscheSE..MicrostructureandfracturemechanicalresonseofwoodJInternational[]ggp[]JournaFrcure20061393--:495508lofat4.,,()-67YoshiharawtH,ModeIIfracturemechanicsroertiesofoodmeasuredbyheasymmetricfouroint[]ppp--bendesunneedechedcmen.EnneerinFracureMechancs2008751:inttsigsilnotsiJiti6gaggpegg,,[]()4727-4739.[68SusantiCME,NakaoT,YoshiharaH.ExaminationoftheMode11fracturebehaviourofwoodwith]mme-ineerashortcrackinanasytricfourpointbendingtestJ.EninFractureMechanics2011[]gg,,7816-:27752788.()69ConradMPCSmF-ithGDFlundG.racturedisctiswoodadheivJ,[],,emofonnuousebonds[]Inernaona-ttilournaladesonandadhesves2003231:3ofhii947.,,j()[70]VeigelS,MiillerU,KeckesJ,etal.Cellulosenanofibrilsasfillerforadhesives:effectonspecific-fractureenerf-e2271gyosolidwoodadhesivebondsJ.Cellulos2011185:1237.[],,()7nd-munmean1VeigelSFollrichJGilAlttterWetal.Comarisonoffractureenertestibsof,,,pgygy[]-ecanerbeam-doubltilevDCBspecimensandlapointtestingmethodforthecharacterizationof()jWooW--adhesivelbondedoodJEuroeanJournaProduct2012701:310.yw.lofdandoods3[]p,,()[72iaoPizhon,WanJiali,DavalosJF.AnalsisoftaeredENFsecimenandcharacterizationof]QggypprfacreunrMoe-UoadnJbondedinteefractudedlig.InternationalJournalofSolidsandStructures[],20408-03:18651884.,() ̄73WangJaaoPzhon.FrandFRPndedceiliiiacturetoughnessofwoodwoodwoodbointerfas[],Qgr--IIJJC200337175undemodeloadin.ournalofomositematerials0;8897.g[]p,,()[74WangVZ,GingerJD,NarayanK.Intralaminarandinterlaminarfracturecharacterizationin]-gluedlaminatedtimbermembersusingimageanalsisJ.EnineerinFractureMechanics201282:y[]gg,,73-84.75陈明祥.弹塑性力学[M].北京:,2010.[]科学出版社76DivosFTanakaTNaaoHetal.Determiionofshearduluoniiij.natmosconstructonszetmber[],,g,[]Wscien-oodceandtech1998:402.nology,,326393()[77IrwinGR,KiesJA.CriticalenergyrateanalysisoffracturestrengthJ.SieMilestoneSeriesMs,][]p1137-997:1311,64.78BenbowJJRoeslerFC.Erim^tsoncontrolledfracturesJ.ProceedinsthePhical,xpeofs[][]gySociety.SectionB1957702:201.,,()[79]WhitneyJM,BrowningCE,MairA.AnalysisoftheflexuretestforlaminatedcompositematerialsJosen-MaaTesandDen145.,Compitterilstigsi974:30[]g,64 大连理工大学硕士学位论文80WilliamsJG,OnthecalculationofenergyreleaseratesforcrackedlaminatesJ.International][][-Journa:11.lofFracture19883620119,,()[81]KanninenMF.Anaugmenteddoublecantileverbeammodelforstudyingcrackpropagationand-arrestJ.Internat1ionalJournaloffracture97391;8392.[],,()82KanninenMF.AdnamicanalsisofunstablecrackroaationandairestintheDCBtest[]yyppgsecimenJlF741045-430.InternationalJournaofracture193:1.p[],,()83OzdCaonLamanae-amnaesecmensJComsilFrlssA.BelsisofanlllitDCBi.ites[],ygpyp[]poScence-andTech159230inology999:5315.,,()[84]WhitneyJM,BrowningCE,MairA.Analysisoftheflexure:testforlaminatedcompositematerials-J.ComositeMaterialsTestinandDesign1974:3045.[]pg,[85]RusselA,StreetKN.Factorsaffectingtheinterlaminarfracturetoughnessofgraphite/epoxy-.inscienceandengesInProrsofICCMIVlaminatesC//Proressinerinofcompoite.:gesdins,[]gggT1-okyo982:279286.,-86nZndw-modeWa.AsaichthreeointbendspecimenfortestinIinterlaminarfracture[]gQpgouneorr-renforcedcomostghssffibeipitematerialsJ.Internationalournaloffracture1997853):[],,j(23-2401.王勖成.200387.有限单元法M北京:清华大学出版社有限公司.[][],88-李重伟.有限元分析方法综述m.天津建设科技2006.增刊:14.[],[89]何盛.樟子松指接材有限元模型分析研宄[D].南京林业大学硕士论文,2011.90JRximateAnlinConceniRice.APathIndeendentInteralandtheAroasisofStratratonb[]pgppyyNCkJA-otchesandracs.JournalofppliedMechanics196835:379386.,,[][91]解德,钱勤,李长安.断裂力学中的数值计算方法及工程应用[M].北京:科学出版社,2009.65 层板胶合木层间断裂初度研究攻读硕士期间发表学术论文情况1徐博激王亚励赵艳华.木材顺纹断裂韧度的研究进展.(己录用).,,力学与实践主办单位一:中国力学学会与中国科学院力学研宄所。本硕士学位论文第章相关内容)(-2Bo-HanXuYa-XunWanYanHuaZhaoInWeiDon.fluenceofrowthrin,,,ggggor-ientationontheinterlaminarfracturetoughnessingluedlaminatedtimber.JournalofMaterialsinCivilEnineerinUnderthesecondreview.TheoranizerThAmigg(:eercan)gSocietyofCivilEnineers.(本硕士学位论文第二、三章相关内容)g-66- 大连理工大学硕士学位论文致谢.时光飞逝,白夠过隙,转眼之间三年的硕士研究生生活就要敲响结束的钟声一,本以为我在写下这段致谢时定是内心充满着雏鸟破壳而出的欣喜和对将要自由飞翔的憧憬一丝感激与不舍的情愫眼睛里浮起了。但是,现在我心中多了,一片雾气。为什么我眼里常含着泪水,因为我对这校园爱得深沉。三年的硕士研,究生生活教会我了太多的东西不只是在学术学业上让我得,更.到了长足的进步是在我正式步入社会进入工作之前给我一段很长的时间去让自己的心智更加成熟。感谢我的导师徐博激副教授三年的言传身教,他以自己突出的学术能力指导着我在科研的道路上不断攀登髙峰,更以自己独特的文化素养和人格魅力感染着.我,让我认识到需要提升自己的综合素质。在向老师表达最真擎的感谢的同时,祝他身体健康一,工作顺利,家庭美满,切顺利。感谢以赵艳华副教授和董伟副教授为代表的结构工程研宄所全体老师对我在.科研瓶颈时给与的指点,帮我解决了诸多难题,使得我可以顺利完成硕士研究生阶段的科研任务。感谢结构试验大厅刘毅老师,王志刚老师和曲秀华老师在我进行力学试验期间给与的帮助与支持,让我能够快速并完美的完成自己的试验任务一。祝全体结构工程专业的老师切顺利。感谢师兄蔡竞,,师弟袁东伟和姚学峰,室友王俊、江培情和陈智以及全体结构工程专业硕博研究生好友对我生活和学业上的关心与照顾,祝你们学业有成,前程似锦。感谢父母二十七年来对我的养育和支持,感谢全部家人对我的关心与照顾。一我参加工作后,定努力工作,尽快成长为家庭里最有力的顶梁柱。最后,感谢百忙之中抽空评阅本文的各位老师。王亚励2015年5月6日于大连理工大学67 层板胶合木层间断裂钥度研究大连理工大学学位论文版权使用授权书本人完全了解学校有关学位论文知识产权的规定,在校攻读学位期间论文工作的知识产权属于大连理工大学,允许论文被查阅和借阅。学校有权保留论文并向国家有关部门或机构送交论文的复印件和电子版,可以将本学位论文的全部或部分内容编入有关数据库进行检索,可以釆用影印、缩印、或扫描等复制手段保存和C编本学位论文。:层板胶合太层间断裂初度研究学位论文题目''^:日期:年作者签名^3一^月」^日导师签名:/辨織日期:讀年月丨0日68

当前文档最多预览五页,下载文档查看全文

此文档下载收益归作者所有

当前文档最多预览五页,下载文档查看全文
温馨提示:
1. 部分包含数学公式或PPT动画的文件,查看预览时可能会显示错乱或异常,文件下载后无此问题,请放心下载。
2. 本文档由用户上传,版权归属用户,天天文库负责整理代发布。如果您对本文档版权有争议请及时联系客服。
3. 下载前请仔细阅读文档内容,确认文档内容符合您的需求后进行下载,若出现内容与标题不符可向本站投诉处理。
4. 下载文档时可能由于网络波动等原因无法下载或下载错误,付费完成后未能成功下载的用户请联系客服处理。
大家都在看
近期热门
关闭