板的稳定(结构数值原理)
弹性力学及有限元法1

Elae Element Method
机械工程与自动化学院
现代设计与分析研究所
张瑞金 Rjzhang@
弹 性 力 学 及 有 限 元 法
第一章 绪论
了解弹性力学的定义;
了解弹性力学研究方法 ; 掌握有限单元法的基本思想; 了解常用有限元计算程序; 课程计划。
绪 论
现有网格基础上,根据有限元计算结果估计计算误差、重新划分网格和 再计算的一个循环过程。 3、由求解线性问题发展到求解非线性问题 许多工程问题如材料的破坏与失效、裂纹扩展等仅靠线性理论根本不能 解决,必须进行非线性分析求解,例如薄板成形就要求同时考虑结构的 大位移、大应变(几何非线性)和塑性(材料非线性);而对塑料、橡 胶、陶瓷、混凝土及岩土等材料进行分析,则必须考虑材料非线性。 4、由单一结构场求解发展到耦合场问题的求解 求解线性结构问题,只要离散单元足够小,所得的解就可足够逼近于精 确值。现在发展方向是结构非线性、流体动力学和耦合场问题的求解。 例如由于摩擦接触而产生的热问题,金属成形时由于塑性功而产生的热 问题,需要结构场和温度场的有限元分析结果交叉迭代求解,即“热力耦 合”的问题。 5、程序面向用户的开放性 商业化的提高要求给用户一个开放的环境。
解析法:得出精确的函数解
数值法: 差分法:采用差商代替微商,将弹力中导 出的微分方程及其边界条件化为差分方程 (代数方程)进行求解。 变分法:根据变形体的能量极值原理,导 出弹性力学的变分方程,并进行求解。 有限单元法:离散模型的数值解
绪 论
弹 性 3. 有限元法基本思想 力 学 及 有 将求解区域划分为有限个互不重叠的单元,单元 之间仅依靠节点连接,单元内部点的待求量可由 限 元 单元节点量通过选定的函数关系插值求得,建立 法
结构设计原理复习资料

二.填空题:1.我国钢材按化学成分可以分为、普通低合金钢两大类。
2.在钢筋混凝土构件中钢筋的作用是替混凝土受拉和。
3.混凝土的强度指标有混凝土的立方体强度、和混凝土轴心抗压强度。
4.混凝土的变形可分为受力变形和。
5.钢筋被混凝土包住,可以保护钢筋免于生锈,保证结构的。
6.公路桥涵设计中所采用的荷载有永久荷载、可变荷载和。
7.当永久作用的效应对结构安全不利时,其作用分项系数取。
8.当结构的状态函数Z服从正态分布时,其可靠指标与Z的成正比。
9.容许应力是以平截面和的假定为基础。
10.近几十年来钢筋混凝土结构计算理论的发展,主要是由容许应力法向发展。
11.钢筋混凝土受弯构件常用的截面形式有矩形、和T形等。
12.钢筋混凝土板可分为整体现浇板和。
13.混凝土保护层是具有足够厚度的混凝土层,它是取钢筋边缘至构件截面表面之间的。
14.肋板式桥的桥面板可分为周边支承板和。
15.梁内的钢筋常常采用骨架形式,一般分为绑扎钢筋骨架和两种形式。
16.为了避免少筋梁破坏,必须确定钢筋混凝土受弯构件的。
17.受弯构件在荷载作用下,各截面上除产生弯矩外,一般同时还有。
18.把配有纵向受力钢筋和腹筋的梁称为。
19.在矩形截面梁中,主拉应力的数值是沿着某一条主拉应力轨迹线逐步增大的。
20.随着剪跨比的变化,无腹筋简支梁沿斜截面破坏的主要形态有斜拉破坏、斜压破坏和。
21.当主拉应力超过混凝土的抗拉强度时,构件便会。
22.钢筋混凝土构件抗扭性能有两个重要衡量指标,它们分别是构件的开裂扭矩和构件的。
23.根据抗扭配筋率的多少,钢筋混凝土矩形截面受扭构件的破坏形态一般可分为少筋破坏、、超筋破坏和部分超筋破坏。
24.在纯扭作用下,构件的裂缝总是与构件纵轴成方向发展。
25.扭矩和抗扭刚度的大小在很大程度上取决于的数量。
26.普通箍筋的作用是防止纵向钢筋,并与纵向钢筋形成钢筋骨架,便于施工。
27.轴压柱中,螺旋箍筋的作用是使截面中间部分混凝土成为,从而提高构件的承载力和延性。
喷丝板结构对熔融纺丝挤出过程聚合物熔体流动特性的影响

第31卷㊀第6期2023年11月现代纺织技术AdvancedTextileTechnologyVol.31ꎬNo.6Nov.2023DOI:10.19398∕j.att.202304032喷丝板结构对熔融纺丝挤出过程聚合物熔体流动特性的影响沈泽坤ꎬ王㊀会ꎬ应起繁(东华大学环境科学与工程学院ꎬ上海㊀201620)㊀㊀摘㊀要:熔体在喷丝板微孔内流动时的稳定性和流场分布的均匀性是后续决定纤维成型质量的关键ꎬ这对熔体挤出前后的速度差㊁流道内的剪切速率分布和口模段内径向速度均匀性提出了要求ꎮ利用计算流体力学技术对涤纶工业熔融纺丝中的聚合物微孔挤出过程进行了数值模拟ꎬ得到了聚酯熔体在微孔内流动过程的速度㊁压力和剪切速率分布ꎬ讨论了熔体挤出前后的速度差和剪切速率分布对熔体流动稳定性的影响ꎮ提出了评价口模段内熔体径向速度分布的流动非均匀系数ꎬ指出不同长径比和收敛角分别通过改变流动充分发展段长度和口模段入口处径向速度分量来影响流动非均匀性ꎮ研究发现非均匀系数随长径比的增大而减小ꎬ随收敛角的增大先减后增ꎮ综合分析结果表明ꎬ长径比为3㊁收敛角为74ʎ的喷丝板最佳ꎮ关键词:喷丝板结构ꎻ熔融纺丝ꎻ非牛顿流体ꎻ数值模拟ꎻ流变学ꎻ非均匀系数中图分类号:TQ342.21㊀㊀㊀文献标志码:A㊀㊀㊀文章编号:1009 ̄265X(2023)06 ̄0080 ̄12收稿日期:20230425㊀网络出版日期:20230809作者简介:沈泽坤(1999 )ꎬ男ꎬ江苏南通人ꎬ硕士研究生ꎬ主要从事多相流热物理方面的研究ꎮ通信作者:王会ꎬE ̄mail:huiwang@dhu.edu.cn㊀㊀2022年中国化学纤维产量为6698万吨ꎬ占全球化纤总产量的70%以上[1]ꎬ对中国国民经济的发展起着至关重要的支撑作用ꎮ以聚酯(PET)纤维为代表的合成纤维[2]具有高强度㊁高模量㊁优异的力学性能及高热稳定性[3]ꎬ结构分布均匀㊁取向程度高的高性能纤维可满足国防[4]㊁军工[5]㊁海洋工程[6]㊁航空航天[7]及5G通信[8]等领域的应用ꎬ已成为重要的战略物资ꎮ目前高均匀性聚酯纤维的制备仍然是一个巨大的挑战ꎬ纤维的轴向丝径不均匀ꎬ限制了其力学性能的提高[9]ꎮ聚酯纤维主要采用熔融纺丝法[10]制备ꎬPET固体切片熔融后在喷丝板微孔流道内流动并挤出冷却成型ꎮ高端熔融纺丝技术的核心之一是提高喷丝板内熔体流动的稳定性和均匀性ꎬ强化熔体内部取向程度ꎮ流动过程中对速度场和应力场等物理场的控制至关重要[11]ꎮ熔体流动过程中速度的非均匀分布ꎬ会引起应力㊁取向[12]㊁结晶[13]等一系列不均匀ꎬ最终在成品纤维中表现出力学性能下降ꎮ中国高端熔融纺丝装备长期依赖国外进口的行业现状下ꎬ有必要对熔体微孔流动过程优化控制理论进行研究ꎮ喷丝板是纺丝组件中熔体最后流经并定型的部件ꎬ因此喷丝板内熔体的流动特性对后续纤维的品质有着至关重要的影响ꎮ针对聚合物熔体在喷丝板微孔内的流动过程ꎬ国内外学者已经开展了一定的研究ꎮ由于喷丝板微孔流道大多不超过5cmꎬ熔体在喷丝板内流动时间很短ꎬ多数研究者将这一过程近似认为等温流动ꎮGan等[14]利用Fluent对熔体在微孔内的速度和压力进行了模拟ꎬ认为熔体速度㊁压力的分布与喷丝板流道结构有着密切联系ꎮSuresh等[15]采用CFD模拟的方法ꎬ研究了不同幂律指数值和不同熔体入口流速下喷丝板内部熔体的流动特性ꎬ研究表明75ʎ收敛角下纤维膜性能更佳ꎬ且优化效果与熔体流动过程中熔体的流速㊁剪切速率等特性有关ꎻ张伟等[16]采用数值模拟对单个喷丝孔内的流动进行了分析ꎬ研究指出喷嘴长径比太大会使得熔体受口模段剪切作用增强ꎬ不利于速度均匀分布ꎻ吴金亮等[17]设计了不同的喷丝板结构分别开展实际生产试验ꎬ得到了制备某种特定型号的超细旦多孔聚酯预取向丝的最佳喷丝板结构ꎻ孙华平等[18]利用POLYFLOW软件对具有圆形导孔和锥形导孔的不同喷丝板结构流道内的聚合物熔体流动过程进行了数值模拟ꎬ对比分析了不同结构下熔体的压力㊁速度和剪切速率分布ꎬ该研究结果表明ꎬ锥形导孔相比圆形导孔流道内压力降较为平缓ꎬ速度分布和剪切速率变化更均匀ꎬ有利于流体稳定和纤维成型ꎮ付丽等[19]采用Phan ̄thien ̄tanner(PTT)本构模型ꎬ运用有限元分析方法ꎬ对高密度聚乙烯改性超高分子量聚乙烯共混物熔融法挤出初生丝的过程进行了数值模拟ꎬ获得了速度场和剪切速率场的分布ꎬ模拟中将熔体挤出过程近似成等温流动ꎬ不考虑温度对熔体黏性的影响ꎮ近年来有学者指出ꎬ熔体流动的稳定性与喷丝板微孔流道内熔体的速度分布和剪切速率分布有关ꎮ付丽[20]指出熔体在流经喷丝板时产生的速度差和速度波动越大ꎬ熔体的流动过程就越不稳定ꎮ赵力宁等[21]研究了熔体中晶体的生长演化行为ꎬ发现当熔体剪切速率分别高于或低于某一阈值时晶体以球状和枝晶形态生长ꎬ高剪切速率下熔体取向程度和一致性更佳ꎮDeKort等[22]发现液晶聚合物加工过程中剪切速率的增加ꎬ对应取向参数的松弛速率增加ꎬ冷却后性能提高ꎮ喷丝板内部熔体流动的均匀性普遍认为与流道内压力有关ꎬ更大的压力有利于熔体径向速度分布的均匀性ꎮ顾家耀[23]最早发现如果喷丝板内熔体压力较低ꎬ流道内空隙体积变多ꎬ将使得熔体流动的雷诺数增大ꎬ产生更多紊流ꎬ并进一步恶化径向方向上压力分布的均匀性ꎬ使熔体径向速度分布差异明显ꎬ最终导致纤维成品的不均匀ꎮ邹爱国[24]在此基础上对纺丝组件进行了优化ꎬ更换调配了喷丝板上方的过滤网ꎬ改善了喷丝板内的熔体压力ꎬ实验结果表明改善工艺后的纤维成品结构更均匀ꎬ性能更好ꎮ然而目前缺少直接衡量口模段内熔体径向速度分布均匀程度的参数ꎬ无法从径向速度分布均匀性的角度对喷丝板的设计提供指导ꎮ本文数值模拟聚酯熔融纺丝过程中PET聚合物熔体在喷丝板微孔流道内的速度场㊁压力场和剪切速率场ꎬ在分析不同入口流速和不同喷丝板微孔结构下熔体挤出前后速度差㊁流道压降和自由段剪切速率分布特性对熔体流动稳定性影响的基础上ꎬ提出与口模段内熔体径向速度有关的非均匀系数ꎬ进一步阐明微孔流道结构参数对熔体流动均匀性的影响机制ꎬ指出对应最佳熔体流动非均匀系数的喷丝板结构ꎬ为喷丝板的优化设计提供参考ꎮ1㊀物理模型与理论模型1.1㊀喷丝孔流道物理模型采用PRB64 ̄24 ̄0.25∗0.75L型喷丝板ꎬ其中包含24个喷丝孔ꎬ单个喷丝孔物理模型如图1所示ꎮ整个喷丝孔流道包括入口段㊁输送段㊁收敛段和口模段四部分ꎬ孔外为自由段ꎮ其中外径d1为3.8mmꎬ内径d2为2.5mmꎬ微孔直径D为0.25mmꎬ输送段长度L为7mmꎬ口模段长度l1为0.75mmꎬ长径比l1∕D为3ꎬ收敛角为74ʎꎬ自由段长度为3mmꎮ图1㊀喷丝孔二维模型Fig.1㊀Two ̄dimensionalmodelofspinneret1.2㊀数学模型聚合物熔体流动的特点主要是高黏度和低雷诺数ꎬ满足工程条件的同时为简化计算需对熔体流动作出如下必要假设ꎮPET熔体不可压缩的假塑性流体ꎬ在流道内做稳态流动ꎬ壁面处速度为零ꎬ只考虑黏性力ꎮ根据上述假设ꎬ可得方程形式如下所示:连续性方程:Ñ V=0(1)式中:Ñ为微分算子ꎻV为速度矢量ꎮ动量方程:dVdt=-Ñpρ+Ñτρ(2)式中:ρ为密度ꎻτ为应力张量ꎻp为压力ꎮ聚合物熔体流动本构方程采用Carreau黏弹型本构方程:η=η0[1+(λγ)2]n-12(3)式中:η为黏度ꎻη0为零剪切黏度ꎻλ为松弛时间ꎻγ为剪切速率ꎻn为非牛顿指数ꎮ该模型适用于较大的剪切速率变化范围ꎬ能够准确反映黏度变化特点ꎬ因此选择Carreau模型作为本构方程ꎮ模拟涉及的材料物性参数如表1所示ꎮ18第6期沈泽坤等:喷丝板结构对熔融纺丝挤出过程聚合物熔体流动特性的影响表1㊀PET物性参数Tab.1㊀PhysicalparametersofPETPET物性参数参数值密度ρ∕(kg m-3)1197.44比热容CP∕(J (kg K)-1)2043.54导热系数k∕(W (m K)-1)0.22非牛顿指数n0.91松弛时间λ∕s0.01零剪切黏度η0∕(Pa s)162.541.3㊀网格划分和边界条件由于入口段变化率和长度相对完整流道占比很小ꎬ对熔体后续流动几乎不产生影响ꎬ因此数值模拟以输送段起点作为流动入口ꎮ采用ICEM软件对喷丝孔流道二维模型计算域进行网格划分ꎬ对收敛段㊁口模段以及自由段网格进行加密ꎬ得到的结构化网格如图2所示ꎮ采用计算流体力学ANSYSFLUENT软件模拟熔融纺丝过程中PET熔体的等温稳态层流过程ꎬ边界条件设置如下ꎮa)入口(Inlet):V=1.0m/sꎬT=290ħꎮb)无滑移壁面(Walldie):Vn=Vs=0ꎬT=290ħꎮc)无剪切壁面(Wallfree):τx=τy=0ꎮd)出口(Outlet):P=0Paꎮ图2㊀计算域网格示意Fig.2㊀Griddiagramofcomputingdomain1.4㊀网格无关性与模型验证设置5套不同数量网格进行网格无关性验证ꎬ分别为8096㊁19734㊁39974㊁50094㊁60214ꎮ如图3所示ꎬ选取计算域网格中心对称轴线上的熔体速度作为监测数据ꎬ网格数量大于39974时ꎬ熔体速度几乎不再发生变化ꎬ最大误差不超过3%ꎮ因此综合考虑ꎬ确定后续模拟采用50094的网格数量ꎮ在张伟等[15]所选用的工况下进行模型验证ꎬ将口模段㊁收敛段和部分输送段的中心对称轴线速度分布与文献数据进行对比如图4所示ꎬ可以发现最大误差约为9.7%ꎬ具有较好的一致性ꎮ图3㊀网格无关性检验Fig.3㊀Gridindependencetesting㊀㊀㊀图4㊀模型验证Fig.4㊀Modelverification28 现代纺织技术第31卷2㊀结果与分析分别探究不同入口流速㊁不同收敛角和不同长径比对熔体在喷丝孔流道内流动特性的影响ꎬ各工况如表2所示ꎮ表2㊀模拟工况汇总Tab.2㊀Simulationconditionsummary序号入口流速∕(m s-1)长径比收敛角∕(ʎ)10.537421.037432.037443.037451.027461.037471.047481.057491.0354101.0374111.0396121.031162.1㊀入口流速对喷丝孔流道内PET熔体稳定流动的影响㊀㊀图5为入口流速对速度场的影响及计算域中心对称轴线上的速度分布ꎮ可以明显发现流道内整体速度分布都随入口速度的增大而增大ꎬ且存在着相同的变化趋势ꎬ收敛段之前流速基本不变ꎬ在收敛段开始增大ꎬ口模段内达到最大值ꎬ此处径向存在速度梯度ꎬ速度大小沿径向从表面向中心递增ꎮ挤出口模后ꎬ由于失去壁面束缚ꎬ速度降低至一稳定值不再发生变化ꎮ随着速度增大ꎬ熔体挤出后的速度变化即自由段内平均速度相比口模段分别减小9.16%㊁8.91%㊁9.17%㊁9.30%ꎬ入口流速为1.0m∕s时ꎬ挤出前后速度变化最小ꎬ流动最稳定ꎮ图6为入口流速对压力场的影响及计算域中心对称轴线上的压力分布ꎮ可以发现熔体在喷丝孔通道内流经输送段时的压力基本不变ꎬ而在收敛段由于速度的急速升高压力开始下降ꎬ在口模段沿流动方向压力不断减小ꎬ挤出后压力降低为恒定值ꎮ高入口流速放大了收敛段和口模段的剪切束缚作用ꎬ流道内最大压力随入口流速的增大而显著增大ꎮ不同熔体入口流速下流道最大压力均满足挤出要求ꎬ此时流道内熔体均较为密实ꎬ均匀性好ꎬ但出入口压降的增大也会造成熔体内应力的恶化ꎬ压降越小ꎬ所得初生纤维制品内应力分布更佳ꎬ无规则变形越小ꎮ因此流道内压力并不是越大越好ꎬ入口流速在0.5~1.0m∕s范围内时最大压力变化较为稳定ꎬ流道内压降最大相差24.39%ꎬ处于合理范围ꎬ表明应选择0.5~1.0m∕s的入口流速范围ꎮ图7为入口流速对剪切速率场的影响及计算域中心对称轴线上的剪切速率分布ꎮ可以看出ꎬ受到速度梯度增大的影响ꎬ熔体在自由段内的中心平均剪切速率随入口流速增大而增大ꎬ从2594.97s-1增加到17238.19s-1ꎬ熔体取向程度高ꎬ有利于成丝质量ꎮ而当流速为2.0m∕s和3.0m∕s时ꎬ流道内熔体中心最大剪切速率均达到105数量级ꎬ更易导致熔体破裂ꎬ不利于稳定流动ꎮ综合模拟结果宜选择1.0m∕s的入口流速ꎮ㊀㊀㊀图5㊀入口流速对速度场的影响及对称轴上的速度分布Fig.5㊀Influenceoftheinletvelocityonthevelocityfieldandthevelocitydistributionalongtheaxisofsymmetry38 第6期沈泽坤等:喷丝板结构对熔融纺丝挤出过程聚合物熔体流动特性的影响㊀图6㊀入口流速对压力场的影响及对称轴线上的压力分布Fig.6㊀Influenceoftheinletvelocityonthepressurefieldandpressuredistributionalongtheaxisofsymmetry图7㊀入口流速对剪切速率场的影响及对称轴上的剪切速率分布Fig.7㊀Influenceoftheinletvelocityontheshearratefieldandshearratedistributionalongtheaxisofsymmetry 48 现代纺织技术第31卷2.2㊀长径比对喷丝孔流道内PET熔体稳定流动的影响㊀㊀图8为长径比对速度场的影响及计算域中心对称轴线上的速度分布ꎮ可以发现随着长径比的增大ꎬ自由段内的高速区域减少ꎬ这是因为熔体在口模段内受到壁面束缚作用的程度最强ꎮ随着长径比增大ꎬ挤出前后自由段内平均速度相比口模段分别减小9.27%㊁8.91%㊁9.65%㊁9.88%ꎬ长径比为2和3时流动更稳定ꎮ图9为长径比对压力场的影响及计算域中心对称轴线上的压力分布ꎮ可以看出随着长径比增大ꎬ口模微孔流道变长ꎬ由于口模段的剪切束缚作用最强ꎬ整个喷丝孔流道内产生的最大压力升高ꎬ压降也相应增大ꎮ不同长径比下流道内最大压力均达到挤出要求ꎬ长径比为4和5时ꎬ喷丝孔流道压降相比长径比为2和3时分别增大45.45%㊁42.86%ꎬ考虑内应力分布宜选择长径比为2或3的喷丝板ꎮ图10为长径比对剪切速率场的影响及计算域中心对称轴线上的剪切速率分布ꎮ可以发现此时流道内最大剪切速率相近ꎬ不同长径比下自由段内平均剪切速率分别为5639.25㊁5939.03㊁8049.92㊁3456.71s-1ꎬ长径比为2㊁3㊁4时熔体取向程度较高ꎬ但长径比为2和4时自由段内分别存在一较高的剪切速率极大值15402.39s-1和15754.31s-1ꎬ熔体更易破裂ꎬ不利于稳定流动ꎮ综合模拟结果宜选择长径比为3的喷丝板ꎮ㊀㊀㊀㊀图8㊀长径比对速度场的影响及对称轴上的速度分布Fig.8㊀Influenceofthelength ̄diameterratioonthevelocityfieldandvelocitydistributionalongtheaxisofsymmetry㊀㊀㊀㊀图9㊀长径比对压力场的影响及对称轴上的压力分布Fig.9㊀Influenceofthelength ̄diameterratioonthepressurefieldandpressuredistributionalongtheaxisofsymmetry58 第6期沈泽坤等:喷丝板结构对熔融纺丝挤出过程聚合物熔体流动特性的影响图10㊀长径比对剪切速率场的影响及对称轴上的剪切速率分布Fig.10㊀Influenceofthelength ̄diameterratioontheshearratefieldandshearratedistributionalongtheaxisofsymmetry2.3㊀收敛角对喷丝孔流道内PET熔体稳定流动的影响㊀㊀图11为收敛角对速度场的影响及计算域中心对称轴线上的速度分布ꎮ口模段内的速度总体上随着收敛角的增大而增大ꎬ但收敛角为54ʎ和74ʎ时速度差异不大ꎬ收敛角增加到96ʎ时口模段内速度开始上升ꎬ此后继续增大收敛角ꎬ口模段内速度将显著提高ꎮ不同收敛角下熔体挤出后速度下降至一相似值不再变化ꎬ挤出前后自由段内平均速度相比口模段分别降低9.11%㊁8.91%㊁11.31%㊁18.48%ꎬ考虑到流动稳定性宜选择收敛角为54ʎ和74ʎ的喷丝板ꎮ图12为收敛角对压力场的影响及计算域中心对称轴线上的压力分布ꎮ可以发现随着收敛角增大ꎬ喷丝孔流道内最大压力减小ꎮ收敛角在54ʎ~96ʎ范围内时流道最大压力变化较为稳定ꎬ此收敛角范围内流道内压降最大相差14.29%ꎬ处于合理范围ꎬ有利于纺丝过程的稳定性ꎬ因此宜选择54ʎ~96ʎ的收敛角ꎮ图13为收敛角对剪切速率场的影响及计算域中心对称轴线上的剪切速率分布ꎮ随着收敛角增大ꎬ流道内熔体中心最大剪切速率从39679.53s-1增加到70936.11s-1ꎬ熔体更易破裂ꎬ流动越不稳定ꎮ自由段内熔体中心平均剪切速率随收敛角增大从6261.71s-1下降到2738.05s-1ꎬ熔体取向程度降低ꎮ综合模拟结果ꎬ宜选择收敛角为54ʎ和74ʎ的喷丝板ꎮ68 现代纺织技术第31卷㊀㊀㊀㊀图11㊀收敛角对速度场的影响及对称轴上的速度分布Fig.11㊀Influenceoftheconvergenceangleonvelocityfieldandvelocitydistributionalongtheaxisofsymmetry㊀㊀㊀图12㊀收敛角对压力场的影响及对称轴上的压力分布Fig.12㊀Influenceoftheconvergenceangleonthepressurefieldandthepressuredistributionalongtheaxisofsymmetry78 第6期沈泽坤等:喷丝板结构对熔融纺丝挤出过程聚合物熔体流动特性的影响图13㊀收敛角对剪切速率场的影响及对称轴上的剪切速率分布Fig.13㊀Influenceoftheconvergenceangleontheshearratefieldandshearratedistributiontheaxisofsymmetry2.4㊀流动非均匀性系数口模段内流场的径向均匀性直接决定了后续初生纤维的结构均匀性ꎬ对于成丝质量的影响至关重要ꎮ定义二分之一口模段截面处熔体最大速度与截面平均速度之比u∗为熔体流动非均匀系数ꎮ图14为入口流速㊁长径比和收敛角对流动非均匀系数的影响ꎮ可以看出ꎬ入口流速对口模段内流动均匀性影响较小ꎬ不同流速下非均匀系数最大仅相差0.22%ꎮ这是因为在只改变入口流速的情况下ꎬ熔体流动的均匀性主要取决于喷丝孔流道几何结构ꎮ增大长径比时ꎬ流动非均匀系数不断下降ꎮ这是由于熔体高黏度低热导率特性下普朗特数很高ꎬ达到106数量级ꎬ充分发展段长达数千米ꎮ而口模段长度受制于装备尺寸一般不会超过2mmꎬ因此流动均未充分发展ꎮ随着长径比的增大ꎬ口模段二分之一处距离流动入口越远ꎬ边界层越厚ꎬ所以流动越均匀ꎬ表现为非均匀系数减小ꎮ而收敛角对流动均匀性的影响较为明显ꎬ不同收敛角下非均匀系数最大相差达10%ꎮ可以发现收敛角大于96ʎ后流动非均匀性骤增ꎮ这是因为大角度情况下口模段入口处熔体径向速度分量变大ꎬ进入口模段后轴向速度减小ꎬ熔体高速流动区域从而减少ꎮ㊀㊀图14㊀入口流速㊁长径比和收敛角对流动非均匀系数的影响Fig.14㊀Influenceoftheinletvelocityꎬlength ̄diameterratioandconvergenceangleonnon ̄uniformitycoefficientofflow88 现代纺织技术第31卷3㊀结㊀论本文对于PET聚合物熔体在微孔流道内的流动过程ꎬ在分析不同喷丝板微孔结构和不同入口流速下熔体流变过程中速度分布㊁压力分布和剪切速率分布特性对熔体流动稳定性影响的基础上ꎬ提出了与熔体速度有关的非均匀系数ꎬ定量对比了不同喷丝板微孔结构和不同入口流速下口模段内熔体流动的均匀性ꎬ进一步阐明了微孔流道结构参数对熔体流动均匀性的影响机制ꎬ指出了对应最佳熔体流动非均匀系数的喷丝板结构ꎬ为喷丝板的优化设计提供了参考ꎮ主要结论如下:a)随着入口流速的上升熔体挤出前后的平均速度变化程度先减后增ꎬ1.0m∕s时流动最稳定ꎻ入口流速在0.5~1.0m∕s范围内时压降变化更为稳定ꎬ且更有利于初生纤维内应力ꎻ流速的上升增大了自由段内平均剪切速率ꎬ强化了熔体取向程度ꎬ但入口流速为2.0m∕s和3.0m∕s时流道内最大剪切速率达105数量级ꎬ熔体破裂的可能性更高ꎬ不利于流动稳定性ꎮ喷丝板几何结构不变的情况下入口流速对熔体流动非均匀系数影响不大ꎮ研究表明应选择1.0m∕s的入口流速ꎮb)口模长径比为2和3时挤出前后平均速度变化程度小ꎬ流动更稳定ꎻ长径比为2和3时压降更有利于初生纤维内应力ꎻ长径比为2㊁3㊁4时自由段内平均剪切速率高ꎬ即熔体取向程度较高ꎬ但长径比为2和4时自由段内仍明显存在剪切速率极大值ꎬ易导致熔体破裂不利于流动稳定性ꎮ同时ꎬ增大长径比时由于熔体始终处于充分发展段ꎬ因此边界层越厚ꎬ流动越均匀ꎬ表现为非均匀系数减小ꎮ研究表明长径比为3喷丝板最佳ꎮc)收敛角为54ʎ和74ʎ时挤出前后平均速度变化程度小ꎬ流动更稳定ꎻ收敛角在54ʎ~96ʎ范围内时压降变化较为稳定ꎬ最大相差14.29%ꎬ处于合理范围ꎻ自由段内熔体平均剪切速率随收敛角增大而下降ꎬ熔体取向程度降低ꎬ但收敛角为54ʎ和74ʎ时流道内熔体中心最大剪切速率较小ꎬ熔体不易破裂ꎬ流动稳定性更好ꎮ同时由于收敛角改变了口模入口处的径向速度分量ꎬ对流动均匀性的影响较为明显ꎮ研究表明收敛角为74ʎ的喷丝板最佳ꎮ参考文献:[1]国家统计局.2022年中国化学纤维产量[EB∕OL].[2023 ̄06 ̄15].https:∕∕data.stats.gov.cn∕easyquery.htm?cn=C01.NationalBureauofStatistics.Chemicalfiberproductionin2022ofChina[EB∕OL].[2023 ̄06 ̄15].https:∕∕data.stats.gov.cn∕easyquery.htm?cn=C01.[2]IRMATOVAMBꎬSODIQALIYEVAFM.Studyingthepropertiesoflocalpolyesterfiberꎬformedfrompolyethyleneterephthalategranules[J].ACADEMICIA:AnInternationalMultidisciplinaryResearchJournalꎬ2022ꎬ12(6):16 ̄22. [3]杨童童ꎬ李世君ꎬ刘峰ꎬ等.阻燃用聚酯纤维的性能及其纺纱工艺探究[J].合成技术及应用ꎬ2020ꎬ35(3):39 ̄44.YANGTongtongꎬLIShijunꎬLIUFengꎬetal.Researchonthepropertiesofflameretardantpolyesterfiberandspinningprocess[J].SyntheticTechnology&Applicationꎬ2020ꎬ35(3):39 ̄44.[4]万雷ꎬ吴文静ꎬ吕佳滨ꎬ等.我国对位芳纶产业链发展现状及展望[J].高科技纤维与应用ꎬ2019ꎬ44(3):21 ̄26.WANLeiꎬWUWenjingꎬLÜJiabinꎬetal.Thecurrentsituationandprospectofpara ̄aramidindustrychaininChina[J].Hi ̄TechFiberandApplicationꎬ2019ꎬ44(3):21 ̄26.[5]李杜ꎬ陈清清ꎬ张玲丽ꎬ等.芳纶1414装甲材料的制备及拉伸性能探讨[J].棉纺织技术ꎬ2022ꎬ50(9):19 ̄23.LIDuꎬCHENQingqingꎬZHANGLilingꎬetal.Preparationandtensilepropertydiscussionofaramid1414armourmaterial[J].CottonTextileTechnologyꎬ2022ꎬ50(9):19 ̄23.[6]赵智垒ꎬ丁永春ꎬ杨中桂ꎬ等.高性能纤维复合材料在海上风电的应用[J].船舶工程ꎬ2022ꎬ44(S1):51 ̄56.ZHAOZhileiꎬDINGYongchunꎬYANGZhongguiꎬetal.Applicationofhighperformancefibercompositesinoffshorewindpower[J].ShipEngineeringꎬ2022ꎬ44(S1):51 ̄56. [7]GALIZIAPꎬSCITIDꎬSARAGAFꎬetal.Off ̄axisdamagetoleranceoffiber ̄reinforcedcompositesforaerospacesystems[J].JournaloftheEuropeanCeramicSocietyꎬ2020ꎬ40(7):2691 ̄2698.[8]任明当ꎬ苏荣ꎬ杨屹.芳纶纤维在光电线缆中的应用探赜[J].光纤与电缆及其应用技术ꎬ2023(1):1 ̄6.RENGMingdangꎬSURongꎬYANGYi.Applicationdiscussionofaramidfiberinopticalandelectricalcables[J].OpticalFiber&ElectricCableandTheirApplicationsꎬ2023(1):1 ̄6.[9]YOUNESB.Astatisticalinvestigationoftheinfluenceofthemulti ̄stagehot ̄drawingprocessonthemechanicalpropertiesofbiodegradablelinearaliphatic ̄aromaticco ̄polyesterfibers[J].AdvancesinMaterialsScienceandApplicationsꎬ2014ꎬ3(4):186 ̄202.98第6期沈泽坤等:喷丝板结构对熔融纺丝挤出过程聚合物熔体流动特性的影响[10]RAOPNꎬSABAVATHGKꎬPAULSN.ImpactofMTAblend%inmeltspinningprocessandpolyesterproperties[J].SNAppliedSciencesꎬ2021ꎬ3(2):1 ̄11. [11]郝浩然.熔融纺丝组件内流体力学性能的研究[D].北京:北京化工大学ꎬ2021.HAOHaoran.StudyonHydrodynamicPropertiesofMeltSpinningPack[D].Beijing:BeijingUniversityofChemicalTechnologyꎬ2021.[12]PAPKOVDꎬDELPOUVENꎬDELBREILHLꎬetal.Quantifyingpolymerchainorientationinstrongandtoughnanofiberswithlowcrystallinity:Towardnextgenerationnanostructuredsuperfibers[J].ACSnanoꎬ2019ꎬ13(5):4893 ̄4927.[13]MISZTAL ̄FARAJBꎬPᶒCHERSKIRBꎬDENISPꎬetal.Modelingoforientedcrystallizationkineticsofpolymersintheentirerangeofuniaxialmolecularorientation[J].Polymerꎬ2019ꎬ173:141 ̄157.[14]GANXHꎬLIUQꎬMAXJ.Thecharacteristicsofmeltflowincompositespinningmicropore[J].AdvancedMaterialsResearchꎬ2011ꎬ383:110 ̄115.[15]SURESHKꎬSELVAMKꎬKARUNANITHIB.CFDsimulationstudiesontheflowbehaviorofpower ̄lawfluidsusedtoextrudethepolymerichollowfibermembranethroughanangularspinneret[C]∕∕AIPConferenceProceedings.ChennaiꎬIndia.AIPPublishingꎬ2019:160 ̄160. [16]张伟ꎬ成文凯ꎬ张先明.涤纶工业丝熔融纺丝过程的数值模拟[J].现代纺织技术ꎬ2022ꎬ30(5):52 ̄59.ZHANGWeiꎬCHENGWenkaiꎬZHANGXianming.Numericalsimulationofmeltspinningprocessofpolyesterindustrialyarn[J].AdvancedTextileTechnologyꎬ2022ꎬ30(5):52 ̄59.[17]吴金亮ꎬ王铁军ꎬ孙福ꎬ等.细旦多孔聚酯预取向丝Barmag环吹纺丝组件的设计[J].现代纺织技术ꎬ2021ꎬ29(6):49 ̄54.WUJinliangꎬWANGTiejunꎬSUNFuꎬetal.DesignofBarmagringblowingspinningassemblyoffinedenierporouspolyesterpre ̄orientedyarn[J].AdvancedTextileTechnologyꎬ2021ꎬ29(6):49 ̄54.[18]孙华平ꎬ冯培ꎬ杨崇倡.基于Polyflow对扁平复合导电纤维喷丝板孔道优化设计[J].合成纤维工业ꎬ2019ꎬ42(3):69 ̄73.SUNHuapingꎬFENGPeiꎬYANGChongchang.Designofspinneretpilotholeforflatcompositeconductivefibersbasedonpolyflowsoftware[J].ChinaSyntheticFiberIndustryꎬ2019ꎬ42(3):69 ̄73.[19]付丽ꎬ薛平ꎬ刘丽超ꎬ等.PE ̄UHMW∕PE ̄HD熔融挤出初生丝挤出胀大现象模拟分析[J].工程塑料应用ꎬ2019ꎬ47(4):53 ̄58.FULiꎬXUEPingꎬLIULichaoꎬetal.Simulationanalysisofextrusionswellinginas ̄spunfilamentofthePE ̄UHMW∕PE ̄HDblendsmeltextrusionprocess[J].EngineeringPlasticsApplicationꎬ2019ꎬ47(4):53 ̄58.[20]付丽.超高分子量聚乙烯熔融挤出初生丝过程的模拟与实验研究[D].北京:北京化工大学ꎬ2019.FULi.SimulationandExperimentalStudyontheProcessofMeltExtrusionofUltraHighMolecularWeightPolyethylene[D].Beijing:BeijingUniversityofChemicalTechnologyꎬ2019.[21]赵力宁ꎬ林鑫ꎬ黄卫东.较低剪切速率下过冷熔体非枝晶组织的形成与演化[J].金属学报ꎬ2011ꎬ47(4):403 ̄409.ZHAOLiningꎬLINXinꎬHUANGWeidong.Formationandevolutionofthenon ̄den ̄driticmorphologyinundercoolingmeltwithlowershearingrate[J].ActaMetallurgicaSinicaꎬ2011ꎬ47(4):403 ̄409.[22]DEKORTGꎬLEONÉNꎬSTELLAMANNSEꎬetal.Effectofshearrateontheorientationandrelaxationofavanillicacidbasedliquidcrystallinepolymer[J].Polymersꎬ2018ꎬ10(9):75 ̄76.[23]顾家耀.POY条干不匀率影响因素的探讨[J].合成纤维工业ꎬ1993ꎬ16(1):53 ̄56.GUJiayao.DiscussiononinfluencingfactorsofunevennessofPOY[J].ChinaSyntheticFiberIndustryꎬ1993ꎬ16(1):53 ̄56.[24]邹爱国.黑色锦纶6POY纺丝工艺[J].广东化纤ꎬ1999(2):1 ̄3.ZOUAiguo.Blacknylon6POYspinningprocess[J].GuandongChemicalFiberꎬ1999(2):1 ̄3.09 现代纺织技术第31卷EffectofspinneretstructureonflowcharacteristicsofpolymermeltinmeltspinningextrusionprocessSHENZekunꎬWANGHuiꎬYINGQifan(CollegeofEnvironmentalScienceandEngineeringꎬDonghuaUniversityꎬShanghai201620ꎬChina)Abstract:Atpresent thepreparationofhighlyhomogeneouspolyesterfibersisstillagreatchallenge.Thestabilityofmeltflowinthemeltspinningprocessandthedegreeofnon ̄uniformityofradialvelocitydistributionbeforeextrusionaffectthequalityofsubsequentspinning.Inthispaper theflowprocessofpolyestermeltinmicroporousflowchannelwasstudiedbynumericalsimulation.Themeltflowcharacteristicsunderdifferentmicroporousstructuresanddifferentinletvelocitiesofthespinneretwereanalyzed.Thenon ̄uniformitycoefficientrelatedtothemeltvelocitywasproposed andtheinfluencemechanismofmicroporousflowchannelstructureparametersonthestabilityanduniformityofmeltflowwasfurtherclarified.Thespinneretstructurecorrespondingtotheoptimumnon ̄uniformitycoefficientofmeltflowwaspointedout.Withtheincreaseoftheinletvelocity thevariationdegreeoftheaveragemeltvelocitybeforeandafterextrusionfirstdecreasesandthenincreases andtheflowisthemoststableat1m∕s.Thetimepressuredropismorestablewhentheinletvelocityrangesfrom0.5m∕sto1m∕s andismoreconducivetotheinternalstressoftheprimaryfiber.Theincreaseoftheflowvelocityincreasestheaverageshearrateinthefreesectionandstrengthensthemeltorientation.However whentheinletflowvelocityis2m∕sand3m∕s themaximumshearrateintheflowpassagereaches105ordersofmagnitude andthepossibilityofmeltfractureishigher whichisnotconducivetotheflowstability.Theinletvelocityhaslittleeffectonthenon ̄uniformitycoefficientofmeltflowwhenthespinneretgeometryisunchanged.Thestudyshowsthattheinletvelocityof1m∕sshouldbeselected.Whenthedielength ̄diameterratiois2and3 theaveragevelocitybeforeandafterextrusionchangesless andtheflowismorestable.Whenthelength ̄diameterratiois2and3 thepressuredropismoreconducivetotheinternalstressoftheprimaryfiber.Whenthelength ̄diameterratiois2 3 4 theaverageshearrateinthefreesectionishigh thatis themeltorientationdegreeishigh.However whenthelength ̄diameterratiois2and4 thereisstillamaximumshearrateinthefreesection whichiseasytoleadtomeltfractureandisnotconducivetoflowstability.Atthesametime asthemeltisalwaysinthefullydevelopedsectionwhenthelength ̄diameterratioisincreased thethickertheboundarylayeris themoreuniformtheflowwillbe andthenon ̄uniformitycoefficientwilldecrease.Theresultsshowthatthespinneretwithalength ̄diameterratioof3isthebest.Whentheconvergenceangleis54ʎand74ʎ theaveragevelocitybeforeandafterextrusionchangesless andtheflowismorestable.Whentheconvergenceangleisfrom54ʎto96ʎ thepressuredropchangesstablywiththemaximumdifferenceof14.29% whichiswithinthereasonablerange.Theaveragemeltshearrateinthefreesectiondecreaseswiththeincreaseofconvergenceangle andthedegreeofmeltorientationdecreases.However whentheconvergenceangleis54ʎand74ʎ themaximummeltshearrateinthecenteroftheflowchannelissmall themeltisnoteasytobreak andtheflowstabilityisbetter.Atthesametime becausetheconvergenceanglechangestheradialvelocitycomponentattheentranceofthedie theinfluenceontheflowuniformityisobvious.Theresultsshowthatthespinneretwithaconvergenceangleof74ʎisthebest.Keywords:spinneretstructure meltspinning non ̄Newtonianfluid numericalsimulation rheology non ̄uniformcoefficient19 第6期沈泽坤等:喷丝板结构对熔融纺丝挤出过程聚合物熔体流动特性的影响。
楼板荷载转换系数计算公式

楼板荷载转换系数计算公式楼板荷载转换系数计算公式是指在建筑结构设计中用来计算楼板荷载的一个重要公式。
楼板荷载是指楼板所承受的荷载,包括活荷载和死荷载。
在建筑结构设计中,需要根据楼板所承受的荷载来确定楼板的尺寸和材料,以确保楼板的安全性和稳定性。
因此,楼板荷载转换系数计算公式对于建筑结构设计具有重要的意义。
楼板荷载转换系数计算公式的基本形式为:Q = α G + β Qk。
其中,Q表示楼板荷载,α和β为楼板荷载转换系数,G表示楼板的死荷载,Qk表示楼板的活荷载。
通过这个公式,可以将楼板的活荷载和死荷载转换成统一的荷载值,从而进行楼板的设计计算。
在实际的建筑结构设计中,楼板荷载转换系数的计算需要根据具体的情况进行。
一般来说,楼板荷载转换系数的计算需要考虑以下几个因素:1. 楼板的跨度,楼板的跨度是指楼板所跨越的空间的距离,跨度越大,楼板荷载转换系数通常越小。
2. 楼板的材料,不同材料的楼板对荷载的承受能力不同,因此需要根据楼板的材料来确定楼板荷载转换系数。
3. 楼板的荷载类型,楼板所承受的荷载可以分为活荷载和死荷载,不同类型的荷载对楼板荷载转换系数的影响也不同。
根据以上因素,可以得出楼板荷载转换系数的计算公式为:α = f1 f2 f3。
其中,f1、f2、f3分别表示楼板的跨度系数、材料系数和荷载类型系数。
这些系数需要根据具体的情况进行计算,以确定楼板荷载转换系数的具体数值。
在实际的建筑结构设计中,楼板荷载转换系数的计算需要结合建筑的具体情况进行。
一般来说,可以通过建筑结构设计规范中的相关公式和表格来确定楼板荷载转换系数的数值。
同时,也可以通过有限元分析等计算方法来得出楼板荷载转换系数的准确数值。
总之,楼板荷载转换系数计算公式是建筑结构设计中一个重要的公式,它可以将楼板的活荷载和死荷载转换成统一的荷载值,从而进行楼板的设计计算。
在实际的建筑结构设计中,需要根据楼板的具体情况来确定楼板荷载转换系数的数值,以确保楼板的安全性和稳定性。
有限元静力分析基本原理

此外,随着大数据和人工智能技术的快速发展,有限元分析可以与这 些技术相结合,实现更加智能化、自动化的工程设计和管理。
THANKS
感谢观看
离散化
将连续的物理系统划分为有限个离散的单元, 每个单元具有一定的形状和大小。
集成
将所有单元的数学方程集成为一个整体的有 限元方程组。
单元分析
对每个离散单元进行数学建模,建立单元的 数学方程。
求解
通过求解有限元方程组,得到物理系统的近 似解。
有限元的数学基础
线性代数
01
有限元方法涉及大量的线性代数运算,如矩阵运算、线性方程
定不变的载荷作用下的响应。
它主要关注的是结构的平衡状态 和位移,而不考虑时间因素和动
态效应。
静力分析广泛应用于工程领域, 如建筑、机械、航空航天等,用 于评估结构的强度、刚度和稳定
性。
静力分析的基本步骤
建立数学模型
首先需要建立结构的数学模型,包括对结构的离散化、选 择合适的单元类型和确定边界条件等。
该方法基于离散化的思想,将 复杂的结构分解为简单的、相 互连接的单元,通过求解每个 单元的平衡方程来获得结构的
整体响应。
有限元静力分析在工程领域中 广泛应用于结构强度、刚度、 稳定性等方面的分析,为结构 设计提供了重要的理论依据和 实践指导。
随着计算机技术的发展,有限 元分析软件不断涌现,为工程 师提供了更加高效、精确的数 值分析工具。
施加载荷
根据实际工况,在结构上施加相应的载荷,包括重力、外 部力、压力等。
求解平衡方程
通过有限元方法,将连续的结构离散为有限个单元,并建 立平衡方程组。然后使用数值方法求解这个方程组,得到 各节点的位移和应力等结果。
结构力学部分名词解释

1、结构按其几何形状可分为杆件结构、薄壁板壳结构和实体结构。
2、结构力学的研究对象是杆件结构。
它是一门研究杆件结构强度、刚度、稳定性和合理组成的科学。
3、杆件结构按其受力特性可分为梁、拱、刚架、桁架、组合结构。
4、结点分为铰结点和刚结点。
铰结点之产生杆端轴力和剪力,不引起杆端弯矩;刚结点除产生杆端轴力和剪力,还引起杆端弯矩,当结构发生变形时,汇交于刚结点各杆端的切线之间的夹角将保持不变。
5、支座的类型:可动铰支座、固定铰支座、固定支座、定向滑动支座。
6、本来是几何可变,经微小位移后又成几何不变的体系称为几何瞬变体系。
7、顺便体系能否应用于工程结构?P8可见,即使荷载不大,也会使杆件产生非常大的内力和变形。
因此,瞬变体系在工程中不能采用,对于接近瞬变的体系也应避免。
8、凡减少一个自由度装置,称一个约束。
一根链杆相当于一个约束;一个单铰相当于两个约束;一个刚性联结相当于三个约束;联结n个刚片的复铰相当于(n-1)个单铰(n为刚片数)9、以刚片作为组成体系的基本部件进行计算的方法称为刚片法。
10、计算自由度W W=3m-2h-r (m刚片数 h 联结刚片的单铰数目r 支座链杆数目)11、平面体系几何不变的必要条件:W>0,表明体系缺少足够的约束,因此是几何可变的;W=0,表明体系具有成为几何不变所必须的最少约束数目;W<0,表明体系具有多余的约束。
12、体系本身为几何不变时必须满足W≤3的条件。
必须指出,W≤0只是几何不变的必要条件,不是充分条件。
13、静定结构与超静定结构的区别:静定结构的几何组成特征是几何不变且无多余约束;超静定结构的几何组成特征是几何不变且有多余约束;仅用静力平衡条件就可以求解的结构称为静定结构;综合运用平衡条件与位移协调条件求解的结构,称为超静定结构。
14、内力图绘制:梁上无荷载(q=0)的区段,Q图为一水平线,M图为一斜直线;梁上有均布荷载(q=常数)的区段,Q图为一斜直线,M图为二次抛物线;集中力作用点的两侧,剪力有突变,其差值等于该集中力,在集中力作用点处,M图是连续的,但因集中力偶两侧的剪力值相同,所以两侧M图的切线应相互平行;集中力偶作用处,剪力无变化,但在集中力偶两侧弯矩有突变,其差值等于该集中力偶,在M图中形成台阶,又因集中力偶两侧的剪力值相同,所以两侧M图的切线应相互平行。
偶应力理论层合梁的稳定性及尺度效应
为 了解释 这一 现 象 , 之发 展 的细 观 理 论 主 为
要有 偶应 力理论 和应 变梯 度理论 。
收 稿 日期 : 0 2—0 21 6—2 0
基金项 目:国家 自然科学基 金资助项 目 ( 目编号 :17 16 项 10 2 5 ) 作者简介 : 陈万吉( 9 1 14 一), , 宁鞍 山人 , 男 辽 教授 , 博士生导师 , 主要研究方 向: 飞机轻 质复合材料结构 高阶强论有限元软体研 发 非线性力学中数学规划法 , 细观偶压力/ 应变梯 度理论 和有 限元 。
1
材 料细 观材料 参数 。偶应 力 理论是 在 细观尺 度上 增 加绕 夹杂 的转 动 平衡 条 件 , 即在 平 衡 方程 中添 加 转动平 衡条 件 ( 即绕 细 观夹 杂 的转 动 平 衡 ) 这 , 是偶应 力理 论 的基 础 , 1 示在 ( , 面 内绕纤 表 Z) 维转 动细观 参数 , 纤维 断面 与基体 相互 作用 , 纤维 被看成 是 影 响 转 动 平 衡 的 夹 杂 。1 示在 ( 2表 Y, Z) 内基 体 内绕 夹 杂 的转 动细 观 参 数 。这 与经 面 典 的偶 应力 理论 引入材 料细 观参数 的方法是 一致
明显增大, 这种现象被称为尺度效应。传统 的弹 性力 学理 论不 考 虑材 料 的 微 观结 构 , 以无 法 对 所
这种 尺度 效应 给 出解 释 。
长度 参数 )N ue ( 4个 细观 材 料 长度 参 , ebr H 含
数 ) 。近 期 有 FekH t isn 等 l — uc no 【 ( 1个 细 c h 等 含 观材 料 长 度 参 数 ) 出 的 偶 应 力 理 论 和 Y g 提 n a 等 _ 提 出的修 正 对称 偶应 力 理论 ( 1个 细观 材 8 含 料长 度参 数 ) 。
结构力学二期末考试试题及答案
结构力学二期末考试试题及答案一、选择题(每题5分,共20分)1. 下列哪项不是结构力学的研究对象?A. 梁B. 板C. 杆D. 流体答案:D2. 静定结构是指:A. 内力和位移都可确定的结构B. 内力可确定,位移不可确定的结构C. 内力不可确定,位移可确定的结构D. 内力和位移都不可确定的结构答案:A3. 力的平衡条件包括:A. 力的合成B. 力的分解C. 力的平衡D. 力的传递答案:C4. 弯矩图的绘制过程中,下列哪项是错误的?A. 弯矩图是力的函数B. 弯矩图是位移的函数C. 弯矩图是力和位移的函数D. 弯矩图是力的积分答案:B二、填空题(每题5分,共20分)1. 结构力学中,______是结构在荷载作用下不发生变形的必要条件。
答案:几何不变性2. 刚体在受力后,若其形状和大小不发生变化,则该刚体处于______状态。
答案:平衡3. 梁的剪力图与弯矩图之间的关系是______。
答案:剪力图是弯矩图的导数4. 根据虚功原理,结构在外力作用下,其虚位移所做的虚功等于______。
答案:外力对虚位移所做的虚功三、简答题(每题10分,共30分)1. 简述什么是超静定结构,并举例说明。
答案:超静定结构是指在静力平衡方程之外还需要其他条件才能确定结构内力和位移的结构。
例如,连续梁和框架结构。
2. 描述一下结构力学中力法的基本步骤。
答案:结构力学中力法的基本步骤包括:建立原结构的刚度方程,引入虚支座,计算虚支座反力,计算弯矩图,最后计算支座反力。
3. 解释什么是结构的稳定性,并举例说明。
答案:结构的稳定性是指结构在外力作用下不发生失稳破坏的能力。
例如,细长杆在受到压缩力作用时,若其长度远大于其截面尺寸,则容易发生屈曲失稳。
四、计算题(每题15分,共30分)1. 给定一个悬臂梁,其长度为L,一端固定,另一端承受一个集中力P,求该梁的弯矩图和剪力图。
答案:略(需要根据具体数值进行计算)2. 给定一个简支梁,其长度为L,两端支撑,梁上分布有均匀荷载q,求该梁的弯矩图和剪力图。
混凝土结构设计课后思考题答案___沈蒲生_第4版
Zz1521.1 结构设计的基本内容及步骤有哪些?试举例说明。
根据结构的概念设计,确定结构材料,结构体系,布置和施工方法;结构分析与设计(其中包括计算简图,内力,变形分析及配筋计算等),结构的构造设计;绘制结构的施工图(其中包括结构布置图,构件末班和配筋图等)1.2 钢筋混凝土梁板结构有几种形式?他们是怎么样划分的?由单向板组成的梁板结构称为单向板梁板结构,由双向板组成的梁板结构称为双向板梁板结构。
当L2/L1<=2时,按双向板设计,2<L2/L1<3时,宜按双向板设计,若按单向板设计时,沿长边方向应配置不少于短边方向25%的受力钢筋。
1.3荷载在整体式单向板结构的板、次梁和主梁中式如何传递的,为什么?按弹性理论和塑性理论计算式两者的计算简图有何区别?单向板以次梁为支座,次梁以主梁为支座,主梁以柱和墙体为支座,作用于结构上的荷载首先由单向板传递给次梁,再由次梁传递给主梁,最后由主梁传递给柱和墙体。
为了减少整体式单向板梁板结构中的跨度,应设置次梁,为了减少次梁的跨度,应设置主梁,为了减小主梁的跨度,应设置柱或墙体;两者计算简图在结构计算跨度处不同。
1.4整体式梁板结构中,欲求结构跨内和支座截面最危险内力时,如何确定活荷载最不利位置?(1)欲求结构某跨跨内截面最大正弯矩时,除恒荷载作用外,应在该跨布置活荷载,然后向两侧隔跨布置活荷载(2)欲求某跨跨内截面最大负弯矩时,除恒荷载作用外,应在该跨不布置活荷载,然后向两侧隔跨布置活荷载(3)欲求结构某支座截面最大负弯矩时,除恒荷载作用外,应在该支座相邻两跨布置活荷载,然后向两侧隔跨布置活荷载(4)欲求结构边支座截面最大剪力,除恒荷载作用外,其活荷载布置与该跨跨内截面最大正弯矩是活荷载布置相同。
欲求结构中间跨支座截面最大剪力时,其活荷载布置和求该支座截面最大负弯矩时活荷载布置相同。
1-5 结构各截面的最大内力值的连线或点的轨迹,即为结构内力包络图。
经典板理论及其变分原理
经典板理论是一种用于研究在弯曲加载下的薄板结构行为的力学理论。
它基于弹性理论和平面假设,适用于较小弯曲和较小挠度的情况。
变分原理是经典板理论的基础,用于推导出平衡方程和位移方程。
以下是关于经典板理论及其变分原理的简要介绍:
1. 经典板理论的基本假设:
-平面应力假设:假定薄板中应力沿板厚方向变化忽略不计。
-线性弹性假设:假定材料服从胡克定律,即应力与应变成线性关系。
-小挠度假设:假定板的挠度较小,即小于板的厚度。
2. 变分原理:
经典板理论中的变分原理是通过最小化势能或最小作用量来推导出平衡方程和位移方程。
其基本思想是,真实的应变状态应该使得系统的势能达到极小值。
3. 平衡方程:
通过应用变分原理,可以得到平衡方程,即拉普拉斯方程。
对于均匀静态情况下的板结构,平衡方程可
以表示为拉普拉斯-贝尔特拉米方程。
4. 位移方程:
利用变分原理,经典板理论可以推导出位移方程,即柯西-里奇扎蒂方程。
该方程描述了板的弯曲变形和剪切变形之间的关系。
5. 应力和应变分布:
经典板理论可以由位移方程得到应力和应变的解析表达式。
这些表达式可用于计算板的内部应力和变形。
经典板理论在设计和分析薄板结构时具有重要的应用价值。
然而,它也有一定的局限性,对于大挠度和非线性材料行为的情况不适用。
因此,在实际应用中需要综合考虑结构的特点和条件,并可能需要使用更高级的理论或数值模拟方法进行分析。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第4章 板的稳定 4.1 板的稳定微分方程 板按照其厚度可分为厚板、薄板和薄膜三种。将板的厚度t与板幅面的最小宽度b相比,如果t/b>1/5~1/8时,板被称为厚板;当1/80~1/100<t/b<1/5~1/8,这个范围内的板称为薄板;而当t/b<1/80~1/100时,板被称为薄膜。薄膜没有抗弯刚度,完全靠薄膜张力来支承横向荷载作用。薄板不仅具有抗弯刚度,还可能存在薄膜张力。 薄板在横向荷载作用下,如果产生的挠度w ≤t/5,则属于薄板的小挠度弯曲问题。此时,薄板中的物理方程和内力表达式与弹性力学中平面应力问题的物理方程和内力表达式相同。薄板中于薄板上下表面等距离的面称为中面。当薄板在中面内承受平行于中面的荷载而失稳时,也可以根据静力平衡准则来确定板的临界荷载。下面根据小挠度理论给出薄板的弹性稳定微分方程 图4-1
2222244224442)2(ywNyxwNxwNywyxwxwDyxyx
(4.1)
式中:)1(1223EtD,为单位宽度板的抗弯刚度。式(4.1)是一个以挠度为未知量的常系数线性四阶偏微分方程。 板的边界条件表达式:
1)简支边:挠度w=0,弯距xM=0,即 0)(2222ywxwD,由于板的边界各点
挠度均为零,则其曲率022yw,故022xw。 2)固定边:挠度w=0,斜率0xw。 3)自由边:弯距xM=0,即 0)(2222ywxwD;剪力0xQ;扭矩xyM=0,均匀分布的扭矩xyM等效于均匀分布的剪力yMxy,xQ与yMxy可合并为0)2(2333yxwx
w。
4.2 受压简支板的弹性失稳 4.2.1 单向均匀受压简支板的弹性失稳 如图4-2所示四边简支矩形板,板的中面上作用有xxPN、xyN=0、yN=0;因此,板的弹性稳定微分方程式(4.1)为
图4-2 0)2(224422444xwPywyxwxwDx (4.2)
根据板的边界条件,当 0x 和 ax时,w=0、022xw、022yw;
0y 和 by时,w=0、022xw、022yw.
符合这些边界条件的板的挠曲面可用二重三角级数表示为 11sinsinmnmnbyna
xmAw
(4.3)
式中:m和n分别是板失稳时,在x和y方向的半波数,mnAnm;,3,2,1,,3,2,1 为各项的待定常数。 对w微分两次和四次后代入偏微分方程,得
11222444224224440sinsin2mnxmnbynaxmamD
Pbnbanmam
Aw
(4.4)
由于axmsin和bymsin均不为零,mnA也不为零,否则板仍然为平面平衡状态,所以 0222244422422444amDPbnbanma
mx
解得 2222mbanambbDPx (4.5) 式中:)1(1223EtD. 只有当n=1时,式(4.5)有最小值,所以有意义的临界荷载为222
mbaambbDPcrx
(4.6)
或 22bDKPcrx (4.7) 式中:K — 为稳定系数,2)(mbaambK. (4.8) 由0dmdK,可得bam,K=4,则
224bD
Pcrx
要使上式成立,bam/ 必须是整数,如果 bam/不是整数,则 222)(mbaambbDPcrx 或 2
2
bD
KPcrx
计算临界荷载时,m的取值应使K值最小。当n=1时,即在y方向成一个半波的条件下,K随m和a / b成曲线关系,如图4-2所示。 板的长宽比在2/ba时,m=1,板以一个半波的形式失稳;长宽比在2与6之间时,m=2,板以两个半波的形式失稳;长宽比在6与12之间时,m=3,余类推。而当长宽比0.4/ba时,K值已非常接近于最小值4minK。 由式(4.6)可得板的临界应力
222)/()1(12tbEKtPcrxcrx
(4.9)
由式(4.9)可知,单向均匀受压板的临界应力与板的宽厚比的平方成反比,而与板的长度无关。对于单向均匀受压的矩形板,当加载边为简支,而非加载边为各种不同的支承条件时,稳定系数K的最小值如表4-1。 表4-1 序号 1 2 3 4 5 非加载边的支承条件 一边简支 一边自由 一边固定 一边自由 两边简支 一边简支 一边固定 两边固定
稳定系数K 0.425 1.280 4.000 5.420 6.970 只有当2/ba时,同时边界条件由简支变为固定,稳定系数K才会有较大提高。通过上述讨论可知,对于单向均匀受压的狭长板,用增加横向加劲肋来改变ba/,从而提高稳定系数的做法并无明显的效果;如果把加劲肋的间距取得小于2b又很不经济。而对于很宽的薄板,如果采用纵向加劲肋以减少板的宽度b倒是有效的。例如在板的纵
向中心加一条加劲肋时2216bPDcrx。
4.2.2 能量法计算简支板的弹性失稳 四边简支的均匀受压板,计算公式中每一项都有bynaxmsinsin,因此,可以从各项中将其分离出来,这样用平衡法求解很方便;而当板的支承条件不是简支时,三角函数则无法分离,这时就需要用能量法来求解。 1)板的总势能 已知 VU,则弹性应变能
dxdyyxwywxwywxwDUab0022222222221
2 外力势能 abxyyxdxdyywxwPywPxwPV0022221 2)瑞利—李兹法 假设符合板的几何边界条件的挠曲面函数为
11),(mnmnyxfAw
将此式代入总势能公式中,经积分后,根据势能驻值原理建立一组0,,0,01211mnAAA的线性代数方程组。其非零解的条件是方程组的系数行列
式为零,即可得板的稳定方程。 3)伽辽金法 已知板的平衡偏微分方程为0)(wL,需假定符合板的几何与自然边界条件的挠曲面函数,现假定挠曲面函数为 niiiyxAw1),(
可得伽辽金方程组
abnababdxdyyxwLdxdyyxwLdxdyyxwL000020010),()(0),()(0),()(
上面方程组经积分后可得到nAAA,,21的线性方程组,为得到它们的非零解,其系数行列式应为零,则可得稳定方程,由稳定方程可解的临界荷载。 4.3 均匀受剪简支板的弹性失稳 均匀受剪的四边简支板如图4-4所示,在其对角线方向因受压而失稳,失稳时板的波长与另一对角线方向的拉力有关。对于长板,失稳时的半波长度约为板宽的1.25倍。当采用能量法求解剪切临界荷载时,板的挠曲面函数可用二重三角函数表示。但是对于均匀受剪的四边简支板,可以利用均匀受剪四边简支的正方形板来求解临界荷载的近似值。 图4-4 如采用伽辽金法求解时,板的中面力yxxyyxxyppNN,而0,0yxNN。板的平衡偏微分方程为 022)(24422444yxwDpywyxwxwwLxy 设满足几何和自然边界条件的挠曲面函数为 ayaxAayaxAw2sin2sinsinsin21 (4.10)
伽辽金方程组为 0sinsin)(00dxdyayaxwLaa (4.11a) 02sin2sin)(00dxdyayaxwLaa (4.11b)
将式(4.10)的偏微分代入式(4.11)经积分后得到
09322124ADpAa
xy
0169322241AaAD
pxy
板的失稳条件是 0169329322424aDpDpaxyxy