工学南京林业大学本科生毕业论文.docx
- 文档编号:12739342
- 上传时间:2023-06-07
- 格式:DOCX
- 页数:36
- 大小:297.23KB
工学南京林业大学本科生毕业论文.docx
《工学南京林业大学本科生毕业论文.docx》由会员分享,可在线阅读,更多相关《工学南京林业大学本科生毕业论文.docx(36页珍藏版)》请在冰点文库上搜索。
工学南京林业大学本科生毕业论文
南京林业大学
本科毕业设计(论文)
题目:
木构件钉结合持久
蠕变研究
学院:
木材工业学院
专业:
木材科学与工程
学号:
070409126
学生姓名:
曾元
指导教师:
杨小军
职称:
讲师
二O一一年六月
木构件钉结合持久蠕变研究
摘要
本文介绍了木构件钉结合持久蠕变在工程中的重要性;阐述了粘弹性本构理论的相关知识;分析了木构件钉结合持久蠕变的影响因素;采用握钉力测试方法对钉结合持久蠕变性能进行了研究,结果表明:
应力水平和温湿度是影响蠕变的决定因素;应力越大,温湿度越高,蠕变破坏越快,蠕变变形量越大;Burgers四元件模型可以用来模拟木构件与钉结合在不同环境、不同应力条件下的蠕变特性;为安全考虑,建议长期使用中应使螺钉应力水平设计为握钉力40%以内。
关键词:
胶合板SPF钉结合持久蠕变
StudyonCreepofNailConnectioninWoodStructureMembers
Abstract
ThisarticledescribestheimportanceofcreepofNailConnectioninWoodStructureMembersinengineering;describestheknowledgeofviscoelasticconstitutivetheory;Italsoanalyzesthefactorsofcreep,usingtheMethodoftestingnailholdingpowerofwoodtostudyoncreep.Theresultsshowedthat:
thestresslevelandtemperatureandhumidityarethedecisivefactorofcreep.Thestressisgreaterandthetemperatureandhumidityishigherthenthefasterthecreepdamaged,creepdeformationislarger;Burgersfour-componentmodelcanbeusedtosimulatethewoodcomponentsandScrewindifferentenvironments,differentcreepcharacteristicsunderstressconditions;Tosecurityconsiderations,theproposedlong-termusestresslevelsshouldbedesignedtoscrewnailholdingpowerof40%orless.
Keywords:
PlywoodSPFNail-connectionCreep
1前言
1.1蠕变研究状况
所谓“蠕变”即在恒定应力作用下,材料的应变随时间的延长而逐渐增大的特性[1]。
随着金属在高温的、电磁的和化学动力条件下的广泛使用。
聚合物、生物体、地质材料、复合材料和智能材料的迅建发展。
材料和结构与时间相关的力学性能、破坏机理及失效准则的研究日益追切。
并取得了很多的新进展。
这也为力学和材料科学提出了许多的重要研究内容和新的学科分支。
蠕变力学即为其中之一蠕变力学考虑材料时间相关的力学行为。
表征材料的蠕变强度,研究一定载荷下物体的变形、应力和位移随时问变化的现象与过程,估算材料和结构的蠕变寿命,进行合理设计。
蠕变力学在航空航天、动力、化工机械、能源地下工程中有广泛的应用,我国蠕变力学的研究,在金属材料、聚合物及其复合材料、岩土、生物、地学等方面都有相当多的研究成果[2]。
1.2课题研究背景
现代木结构住宅建筑在世界上许多国家已很普遍,在美国和北欧,木结构住宅处于市场主导地位;在日本的新建住宅房屋中,有半数以上也是木结构。
20世纪80年代以后,由于我国可用于木结构建筑的木材逐渐减少,用纯木结构建造的建筑被迫停建。
近年来,随着国外木结构住宅进入中国市场,以及国内木结构建筑相关规范的相继出台,在我国停滞了20多年的木结构建筑开始复苏,现代木结构建筑逐渐引起了人们的关注,得到了认可[3]。
由于2×4规格的木材在木结构框架中被大量使用,因此木结构框架俗称为“2×4”结构。
它为北美洲的人们提供了全世界最高的住房标准。
“2×4”房屋之所以倍受青睐,是由于它有如下诸多优点:
a、冬暖夏凉。
这是由于木材具有绝缘、不导热的天然特征。
这一特点带来的好处是大大节省了用于房屋供暖和降温的能源支出。
b、优越的抗震、抗风性能。
c、卓越的透气性。
房屋内的空气得以保持清新洁净。
d、卓越的吸音性能。
一个典型的例子是:
木制住房内的回音很小。
e、利于环保。
生产木制建材所消耗的能源远低于生产水泥和钢材,所产生的固体废弃物空气污染也最少。
同时,加拿大所产木材均来自可再生的森林,并能够被有效地循环使用。
f、易于扩建、改造和修缮。
这是由于木材易于割锯、搬运、操作和安装。
g、易于满足对外观式样与实际功用性的多样需求[4]。
对于木结构建筑而言,我们往往注重研究其防火、防虫、防腐以及防水性能,希望延长其使用寿命。
然而我们却忽略了木材本身的蠕变性能以及用于连接的钉的持久蠕变研究。
对于木材蠕变性能的研究已经相对很成熟,影响木质材料蠕变性能的因素很多,大致可分以下几种:
材料种类、应力水平、加载方式(即载荷特性、加载方向)和环境条件(温度和相对湿度)[5]。
测试木质材料的蠕变特性,可采用压缩、拉伸、弯曲等加载方式,但考虑到试验的可行性,许多学者往往采用三点或四点抗弯试验方式,采取的应力水平范围一般为30%~70%。
环境的相对湿度与木材的含水率有很大关系,而含水率不同,木材的蠕变性能也不相同。
水分赋予木材极大的塑性,尤其是对填充物质的增塑作用。
由于水分的介入,能破坏木材组织内的氢键结构,使木材内部纤维素分子链间距增大,屏蔽大分子间的相互作用力,从而增加了大分子的柔度,使链段易于运动,这给木材组分的运动提供了足够的空间,使得木质材料在长期外载作用下更易于变形。
环境温度影响木质材料的温度。
木质材料内部温度升高,会有两方面的变化:
(1)内含能量水平的瞬间变化;
(2)永久性结构重组。
这是因为,木质材料是由部分结晶的纤维素和部分非晶态纤维素、木质素聚合物组成,当加热到玻璃态转化点温度时,木质材料的非结晶区部分会发生由玻璃质状态向橡胶状态转化,此时的木质素聚合物也失去粘结硬固作用,材料的物理力学性能发生巨大变化,力学强度急剧下降,木质材料的蠕变,尤其是变形中的粘弹性变形和粘性变形部分会显著增加。
在木质材料的玻璃态转化点,它的每个聚合单体都获得足够的能量来减弱相互间的吸引力并激发内部分子运动。
据研究,木材半纤维素和木质素在160~200℃范围软化,一般木材使用环境的温度远远低于此温度,但是往往环境温度和湿度是同时影响材料蠕变特性的含水率增加,木材主要组分的玻璃化转变温度会相应降低,如湿材的半纤维素在
54~56℃、木质素在72~128℃即可发生玻璃化转变现象,这就给木质材料的安全正常使用带来了潜在危险。
J.M.Dinwoodie等学者[6]的研究表明,木质刨花板受弯时,环境温度与相对湿度极大地影响了试材的蠕变形变,且当温度从20℃到30℃时,蠕变变形比从10℃到20℃更大;相对湿度从65%~90%的蠕变比从30%~65%时蠕变更大,并且环境湿度比温度对蠕变的影响更大。
同时针对木质材料的这一特性也研究出了一些关于降低木质材料蠕变变形的措施,如降低木材的吸湿性、加强木材内部组分间连接强度等[7]。
连接是木结构重要的组成部分,而钉连接具有充分的紧密性和韧性,制作简单,安全可靠,是木结构中常用的一种连接。
如屋面板、楼面板采用钉与梁相连,轻木结构中墙骨采用斜钉与底板相连、剪力墙的蒙皮结构以及临时支撑等[8]。
但往往钉的蠕变现象却被我们大多数人所忽视。
既然钉存在此现象我想就应该被研究与探讨,以此促进木结构建筑的进一步发展与完善。
1.3国内外研究现状
国内外的专家一直致力于对金属的持久蠕变性能的研究。
尤其是合金与金属基符合材料的研究。
当然对于蠕变变形机理、蠕变断裂机理以及如何提高蠕变性能的方法已研究的相对透彻与成熟[9]。
由于木结构建筑在我国仍处于初步发展阶段,无论是设计还是施工技术很大程度上都需要提高和完善,在阅读大量文献资料之后发现对于木构件钉合持久蠕变的研究也处于空白阶段,即使是北美、日本、加拿大等木结构建筑体系成熟的国家亦是如此。
1.4研究目的、意义
钉结合是一种操作简便的连接方式,被广泛应用于我国传统木工和现代木结构建筑中,过去使用的钉主要主要是木钉和竹钉,现在主要采用金属钉,包括只直钉、螺钉等。
钉结合会破坏木材纤维且连接的强度相对较低,故主要用在木制品抗剪部位。
有时钉结合仅起辅助作用。
钉结合的牢靠程度常用握钉力来检验,然而握钉力的试验方法周期很短,钉结合在木构件中使用时处于长期受力状态,依据国标测试的结果并不能代表长期受力行为。
木材及木质材料都是粘弹性材料,钉结合也将表现粘弹性
[10]。
为此本课题将针对螺钉研究在木构件中持久受力蠕变特性,以便于钉结合更好应用于建筑木构件的节点处。
同时对于延长木结构建筑的使用年限以及提高其安全性也有一定的借鉴和帮助。
当然这在一定程度上也保证了住户的人生以及财产安全。
1.5课题研究创新点
a、首次对木构件的钉结合持久蠕变进行研究;
b、首次研究不同环境条件下钉结合的蠕变特性规律。
2粘弹性理论
为了模拟粘弹性材料的形变规律,可以将材料视作非均质的不同的质点组成的。
其中一部分质点是纯粘性的,另一部分质点是纯弹性的。
这些质点的不同组配,就构成模拟不同粘弹性特性的各种模型。
作为纯粘性质点的元件一般是阻尼,其特性为形变速率与施力的大小成正比(符合
关系)。
作为纯弹性质点的典型元件是线弹簧(spring),符合
关系。
这里所说的粘弹性模型都是一维的。
即全部元件都在同一个方向上联系,而且力的作用也是沿着这个方向【11】。
2.1Maxwell、Kelvin、Burgers模型
最简单的组合是一个弹性元件和一个粘性元件串联或并联起来。
如图1的(a)和(b)所示。
图2.1蠕变模型
两元件串联的结构成为Maxwell模型,并联称为Kelvin模型(或Voigt模型)。
假设元件组成的模型都是具有单位横截面积和长度的,那么载荷也就是应力σ,位移也就是应变ε。
Maxwell模型(a)受拉伸应力σ时,弹簧和阻尼器中产生的应力是相等的。
整个系统的总伸长等于二元件各自的伸长之和,即:
…………
(1)
其中,
,
,将
带入式
(1),则得到Maxwell模型所描述的应力-应变-时间关系:
上式也可写成:
其中
n/K,称为松弛时间。
Kelvin模型(b)受拉伸时,两元件的伸长是一致的,并且就等于系统的总伸长。
(4)
系统的总应力(载荷)由两元件共同承担。
所以
,或写为一般形式:
…(6)
Maxwell模型由一个弹簧和一个阻尼器串联而成,所以蠕变特性实际是上符合粘性元件特性的,即具有长时间的粘性流动。
即所施加的应力直接由粘性元件(阻尼器)连续不断地响应所传递,而弹簧则始终给出某个固定的弹性应变,其蠕变量随时间线性地增加,在这一点上没有充分反映出粘弹性特性,因为粘弹性特性反映的应变率应该是随时间而减弱的。
这就是Maxwell模型的局限性所在。
再看Maxwell模型的松弛特性:
当应变保持为常值,随着时间的增加,粘性元件起作用。
粘性元件逐渐伸长同时负荷逐渐降低,而弹簧则逐渐收缩,直至趋于完全恢复。
应力以逐渐减小的比率降低而趋于零。
由于Maxwell模型能描述松弛特性而不能确切的描述蠕变特性,所以Maxwell模型又叫做松弛模型。
Kelvin模型在施加恒定载荷时,因为起始加载时弹簧受到阻尼器限制,而后随阻尼器逐渐伸长,同时弹簧也逐渐拉伸,直到全部弹性应变完全发挥时停止。
Kelvin基本体模型在恒定应变的载荷条件下,粘性元件作用不再存在,形变为单纯弹簧状态下的应力响应,故不会有松弛。
因此Kelvin模型又称蠕变模型。
Burgers模型(c),这种模型的连续变形是由瞬间弹性变形、延迟弹性变形和黏性流动变形组成[7,8],即:
(7)
式中:
为瞬时弹性模量,
为延迟弹性模量,
、
分别为黏性系数和延迟黏性系数,τ为延迟时间。
2.2Kelvin、Maxwell标准线性固体模型
图2.2Kelvin、Maxwell标准线性固体模型
图1.4(a)为Kelvin标准线性固体模型,表示弹簧和Kelvin模型串联的情况,力σ在串联的两个环节中是一样的。
伸长ε则为二者之和
,
弹簧中:
Kelvin模型中:
对上式微分后带入各关系式:
(8)
图1.4(b)为Maxwell标准线性固体,表示弹簧和Maxwell基本体模型并联的情况。
两环节的ε相同,而
。
Maxwell模型的伸长由弹簧K1和阻尼器η1共同构成。
通乘以K1/n1后,整理得:
(9)
2.3蠕变曲线特性
采用理论模型是为了把复合材界面蠕变用较直观的方法表示出来。
由于界面材料包括胶粘剂及木材浅表层和碳纤维层,其中木材和胶粘剂均为典型的粘弹性材料,因此界面具有典型的高分子所具备的粘弹性。
在一定的应力水平下,这两种基体物质都会有蠕变的特性,但其复合材料界面是否有这样的特性,在确定蠕变理论模型前,需要对试验曲线作简要分析。
通过对复合材界面蠕变试验的时间曲线进行分析,得出界面蠕变曲线有以下4个特性:
(1)当施加应力后,立刻产生瞬时的弹性应变,所以流变模型中应包括弹性元件。
材料的弹性可用虎克体,即用一弹性元件(弹簧)来表示,其力学响应瞬间完成,本构关系即为虎克定律的表达式:
。
(2)试验中的蠕变曲线都反映出应变随时间增大而增大的趋势,因此流变模型中应包括粘性元件。
材料的线性粘性可用牛顿液体,即一个假想的粘性元件(粘壶)来表示。
其本构关系可写为:
。
式中
为剪切粘滞系数,一般简称为“粘度”。
(3)当施加应力水平后,应变随时间增大,存在着塑性变形,所以流变模型应包括塑性元件。
材料的塑性可用塑性体,即两个相互接触的摩擦滑块来表示。
它们之间在受力时产生摩擦力,最大值为f,当作用应力小于时,应变为零;当达到或大于时,变形为任意值。
(4)在低应力水平作用下的蠕变曲线,随时间的增长,变形将趋于某一定值,与广义Kelvin模型曲线相似。
其蠕变变形由三个部分组成:
第一部分是弹性变形,第二部分是延迟弹性变形,这种变形是随时间而变化的,可以用弹簧和阻尼器并联起来去模拟,第三部分是高分子的相互滑移引起的粘性流动,可以用一个阻尼器模拟。
总变形等于三部分变形的总和。
2.4流变模型选择原则
材料流变模型的选择与确定,应该遵循两个原则:
(1)模型能够很好地反映材料的力学特性;
(2)模型应尽可能简单、直观,便于工程应用。
根据界面材料时间挠度曲线,可用上述三种模型:
Maxwell标准线性固体模型、Kelvin标准线性固体模型和Burgers四元件模型来描述界面材料在短时应力作用下的粘弹性行为。
2.5burgers模型选用
图2.3木材蠕变曲线3阶段图
一般而言,木材在恒定应力作用下典型蠕变曲线可以分为如图2.3所示的3个阶段。
图中第一阶段为初期蠕变(暂时蠕变)阶段,ε-t曲线向上弯曲,其应变速率逐渐递减,因此亦称衰减蠕变阶段。
第二阶段为稳态蠕变(等速蠕变),该阶段内,曲线斜率保持不变,应变速率呈定值稳定状态。
这个区段历时的长短主要取决于应力水平。
第三阶段为试件达到破坏前应变速率呈加速增长的第三期蠕变,这种蠕变将导致试件迅速破坏。
木构件钉结合持久程度主要是钉与木材的摩擦阻力,它们之间相互影响,但通过握钉力试验观察,发现螺钉的钉齿未发生变形,也没有任何破损的情况,木材与钉结合部分有少量木屑产生。
推断其持久蠕变主要取决于木材本身的持久蠕变。
然而描述木材的蠕变特性时,常常选用Burgers蠕变模型【12】,蠕变曲线开始时一般存在一定的瞬时变形,然后剪应变指数递减的速率增长,最后应变速率逐渐趋近于稳定,其蠕变特性与典型的Burgers模型蠕变特性曲线(图2.5)较接近。
Burgers蠕变模型由Kelvin模型与Maxwell模型串联而成。
图2.4Burgers模型图2.5Burgers模型曲线
因此本试验拟采用Burgers模型描述曲线的一、二阶段,根据具体的试验数据检验是否吻合Burgers模型曲线。
3试验材料与方法
3.1试验材料
SPF——云杉(Spruce)松木(Pine)冷杉(Fir),规格:
40x90x4000mm,密度:
4000g/m³,一等材,木材外观明亮、洁净,颜色由白色到浅黄色不等,纹路细微笔直、质地光滑。
胶合板,7层,密度0.648g/cm3。
螺纹钉:
长度40mm,直径2.5mm;长度30mm,直径2.5mm。
3.2试验用仪器与设备
高低温湿热试验箱,型号:
GDS-100;电源电压:
AC380V;工作室尺寸:
270X400X680;电源频率:
50HZ;温度范围:
0~100℃;湿度范围:
15~98%RH;温度波动度:
±0.5℃;温度均匀度:
±2℃。
微机控制电子万能试验机:
型号UTM4304,最大试验力30KN。
其它加工工具:
电动切割机、直尺、螺旋测微仪、螺丝刀、转孔机、铅笔等。
3.3试验方法
3.3.1试验设计
本试验欲通过不同材料在相同应力水平不同温湿度、相同温湿度不同应力水平的蠕变曲线测试,采用相关蠕变模型分析钉在持久作用力下的蠕变特型。
表3.1试验设计数据表
编号
材料
钉长
应力
干球温度
相对湿度
状态表示
SPF-1
SPF
40
70%
25
91%
状态一
SPF-2
SPF
40
85%
25
91%
状态一
SPF-3
SPF
40
40%
60
90%
状态四
SPF-4
SPF
40
60%
60
90%
状态四
SPF-5
SPF
40
70%
40
40%
状态二
SPF-6
SPF
40
40%
50
60%
状态三
SPF-7
SPF
40
50%
60
90%
状态四
SPF-8
SPF
40
60%
60
90%
状态四
J-1
胶合板
40
60%
25
91%
状态一
J-2
胶合板
40
40%
60
90%
状态四
J-3
胶合板
40
60%
60
90%
状态四
J-4
胶合板
40
60%
50
60%
状态三
J-5
胶合板
40
40%
50
60%
状态三
3.3.2试验方法
基材加工:
对SPF、胶合板原料进行锯切粗加工,刨切机精加工,使试件平整光滑,SPF尺寸为90x80x38mm,胶合板尺寸为90x80x18mm。
钻孔:
用钻孔机在试件中心钻口。
SPF孔径2.5mm,孔深25mm;胶合板孔径2.5mm,孔深15mm。
螺钉固定:
用螺丝刀将螺钉小心拧如孔内,保证钉的笔直,SPF拧入深度为25mm,胶合板为18mm。
加应力、改变温湿度:
根据试验设定相应应力与温湿度。
在微机控制电子万能试验机上,根据《GB/T14018-1992木材握钉力试验方法》分别测出SPF、胶合板的握钉力。
图3.1握钉力测试装置
试验机蠕变程序控制设计:
a 始控方式-位移控制:
5mm/min,终控方式-力控制:
xxxN;
b 始控方式-采样间隔:
30;
c 始控方式-力保载:
32400s;
d 始控方式-位移控制:
-50mm/min,终控方式-力控制:
0N。
4试验结果与分析
4.1试验现象
a、试件在锯切时沿顺纹方向的木材较沿横纹方向的难锯切,横纹方向毛刺多。
b、钻孔位置在节子处时较在无节子处加工难度加大,节子处拧螺钉时较困难。
c、蠕变试验开始阶段随着应力的增加,试验机发出“咯吱咯吱”的响声,当达到设定应力时声响逐渐消失。
当听到试验机发出连续吱吱声,从电脑上显示蠕变曲线呈现快速上升趋势,应力逐渐降低时,表明实验结束。
d、取出试件,观察到试件表面与钉结合处有少量木屑。
在高温高湿的环境下试件呈现暗灰色,由于含水率增加,表面变得松软。
e、SPF及胶合板握钉力试验曲线呈现同样规律,随钉拔出位移的增加,握钉力呈线性递增,达极限值后逐渐下降。
胶合板握钉力达极限值后下降较为缓慢,胶合板厚度方向材性差异较小,曲线图上表现明显。
4.2试验结果
4.2.1握钉力试验结果
图4.1SPF握钉力试验曲线图
图4.2胶合板握钉力试验曲线图
表4.1握钉力数据表
编号
材料
握钉力(N)
平均值(N)
变异系数
SPF-a
SPF
2907.16
2892.5
0.052
SPF-b
2646.80
SPF-c
2724.27
SPF-d
2920.99
SPF-e
3048.38
SPF-f
2729.4
J-a
胶合板
1816.12
1815.5
0.044
J-b
1702.04
J-c
1867.38
J-d
1954.26
J-e
1782.87
J-f
1770.33
4.2.2蠕变试验结果
表4.2试验结果表
编号
应力
环境状态
初始变形(mm)
总变形量(mm)
SPF-1
70%
状态一
0.53183
0.75964
SPF-2
85%
状态一
0.94181
2.97802
SPF-3
40%
状态四
1.13557
2.42837
SPF-4
60%
状态四
0.64565
2.32195
SPF-5
70%
状态二
0.74402
2.13181
SPF-6
40%
状态三
0.72119
2.22709
SPF-7
50%
状态四
0.52033
0.26947
SPF-8
60%
状态四
1.2308
2.40719
J-1
60%
状态一
0.71297
0.15172
J-2
40%
状态四
0.65291
2.90589
J-3
60%
状态四
0.38194
3.09898
J-4
60%
状态三
0.82069
2.31984
J-5
40%
状态三
1.88837
0.14867
4.3结果分析
4.3.1蠕变试验曲线分析
图4.3J-4蠕变曲线
图4.4SPF-3蠕变曲线
由图可看出胶合板和SPF握钉蠕变曲线分布规律几乎一致,大致可分为3个阶段:
第一阶段,蠕变是非稳定的蠕变阶段,它的特点是开始蠕变速度较大,但随着时间的推移,蠕变速度逐步减小,第二阶段是稳定阶段的蠕变,它的特点是蠕变以固定的但是对于该应力和温度下是最小的蠕变速度进行,这就在蠕变曲线上表现为一具有一定倾斜角度的直线段。
第三阶段蠕变以迅速增大的速度进行,这是一种失稳状态。
直到发生断裂,至此整个蠕变过程结束
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 工学 南京 林业大学 本科生 毕业论文
文档标签
- 北京大学本科生毕业论文格式
- 中南林业科技大学毕业论文
- 年级北京林业大学本科毕业论文
- 江西理工大学本科生毕业论文
- 新疆农业大学本科生毕业论文
- 南京工业大学生制药
- Bzishtq河南大学本科毕业论文
- 南京林业大学论文工作方案
- 南京工业大学第二十学生
- 中南民族大学本科生毕业论文
- 西南大学毕业生
- 大学本科毕业生就业工作
- 南京工业大学本科生毕业论文
- 南昌大学本科生论文模版
- Hsrwvu河南大学本科毕业论文
- 中南民族大学本科毕业论文
- 桂林理工大学毕业论文正文
- 南京林业大学
- 中南林业科技大学本科生
- 南京农业大学本科生毕业论文
- 燕山大学本科毕业论文格式
- 南京财经大学本科毕业论文格式
- 南京财经大学本科毕业论文格式
- 河北大学本科生毕业论文设计
- 东北农业大学大学生毕业论文
- 河北师范大学本科毕业论文格式