在材料力学教材中剪切系数如矩形截面常取6/5是基于经典梁理论推导得出的用于计算剪切变形或切应力分布。然而在有限元软件如Abaqus中默认的剪切系数可能有所不同因为软件通常采用更通用的理论如Timoshenko梁理论或根据截面形状自动计算例如对于矩形截面Abaqus默认的剪切系数可能不同与教材中的值存在差异。从Timoshenko在1921年提出问题到Cowper在1966年给出严格解答再到后续学者在各种复杂情况下的推广剪切系数的研究历程折射出力学发展的一个缩影从工程直觉到严格理论从简单情况到复杂应用从单一学科到交叉融合。一篇穿越60年的学术探索带你读懂剪切系数背后的力学本质。这个表一定正确吗材料力学写在前文用有限元软件建梁单元模型时你可能注意过一个不起眼的参数——剪切系数或者叫剪切修正因子。大多数人直接用软件默认值很少去想这个系数从何而来、为什么需要它、取多少才算合理。这个看似微不足道的参数连接着梁弯曲理论中最核心的问题如何用一维梁模型去描述三维弹性体的真实变形行为从20世纪初起,当时的工程师们已经熟练运用欧拉-伯努利梁理论来设计各种结构。这个理论假设梁的横截面在弯曲后仍保持平面且垂直于中性轴也就是著名的平截面假设。对于细长梁这个假设足够精确。但随着工程实践的发展人们发现了一个问题梁比较短粗时理论预测的变形总是比实际测量值小而且频率越高、梁越短粗偏差就越明显。问题的根源在于剪切变形被忽略了。1921年Stephen Timoshenko发表了一篇开创性的论文提出了后来以他名字命名的梁理论。这个理论的核心创新在于引入了两个独立的位移变量挠度w和截面转角φ从而将剪切变形纳入考虑范围。然而引入剪切变形带来了一个新问题如何定义等效的剪切应变梁横截面上的剪应力分布是不均匀的。以矩形截面为例剪应力沿截面高度呈抛物线分布中性轴处最大上下边缘为零。这意味着截面上不同点的剪切变形是不同的。但在梁理论中我们需要用一个统一的截面剪切应变来描述整个截面的剪切行为这就需要一个修正系数将实际的非均匀剪应力分布等效为一个均匀分布。这个修正系数就是剪切系数κ。1966年G.R. Cowper发表了一篇被广泛引用的论文《The Shear Coefficient in Timoshenkos Beam Theory》从三维弹性理论出发给出了剪切系数的严格推导和精确公式。这篇论文不仅解决了理论上的争议更为工程应用提供了可靠的依据。本文将系统梳理剪切系数的理论发展、计算方法与工程应用。剪切系数的起源与定义1.1 问题的物理本质要理解剪切系数的由来需要深入理解梁弯曲时横截面上的应力分布规律。考虑一根承受横向荷载的梁任意横截面上既有弯矩M也有剪力V。弯矩产生正应力σ剪力产生剪应力τ。根据材料力学的基本分析对于矩形截面梁剪应力沿截面高度的分布遵循抛物线规律其中A是截面面积h是截面高度y是从中性轴量起的距离。可以看到中性轴处y0剪应力最大为平均剪应力的1.5倍而在上下边缘y±h/2剪应力为零。这种非均匀分布带来一个根本性的问题在Timoshenko梁理论中我们假设整个截面有一个统一的剪切变形γ这个γ等于截面转角φ对x的导数减去挠度w对x的导数。但是如果截面上各点的实际剪应力因而剪应变都不一样我们该如何定义这个统一的γ答案是我们需要某种等效处理。具体来说我们希望找到一个等效剪应力τ使得当截面上所有点都承受这个均匀剪应力时产生的应变能与实际非均匀剪应力分布产生的应变能相等。或者从另一个角度我们希望找到一个等效剪切刚度使得在相同剪力作用下梁的变形与实际三维弹性体的变形一致。这就是剪切系数κ的物理本质它是一个将非均匀剪应力分布等效为均匀分布的修正因子。数学上如果我们定义平均剪应力为τ那么等效剪应力可以表示为τκτ或者更常见的定义是等效剪切应变为γκ其中G是剪切模量。1.2 Timoshenko的原始工作1921年Timoshenko在《Philosophical Magazine》上发表了题为On the correction for shear of the differential equation for transverse vibrations of prismatic bars的论文。在这篇论文中他首次系统地将剪切变形和转动惯量纳入梁的振动方程。Timoshenko的基本思路是梁的变形由两部分组成——弯曲变形和剪切变形。弯曲变形使截面产生转动剪切变形使截面产生相对于弯曲转动方向的附加转动。因此总转角θ可以分解为两部分由于弯曲产生的转角∂w/∂x以及由于剪切产生的转角γ。在Timoshenko的原始表述中他引入了两个独立的变量挠度w(x,t)和截面转角φ(x,t)其中φ代表弯曲产生的转角而剪切应变为γ ∂w/∂x - φ。在建立运动方程时Timoshenko意识到由于剪应力分布的非均匀性直接使用V/GA来计算剪切应变是不准确的。他引入了一个修正系数k后来常用κ表示将剪切刚度修正为kGA。对于矩形截面他建议k 2/3对于圆形截面k 3/4对于薄壁工字钢截面k ≈ A_web/A腹板面积与总面积之比。这些数值是怎么来的Timoshenko并没有给出严格的推导。对于矩形截面k 2/3可以理解为最大剪应力是平均剪应力的1.5倍而等效剪应力应该介于两者之间取倒数正好是2/3。对于圆形截面k 3/4也有类似的直观解释。但这些只是经验性的估计缺乏理论基础。1922年Timoshenko发表了第二篇相关论文进一步完善了理论框架但仍然没有解决剪切系数的精确定义问题。在接下来的几十年里不同学者提出了不同的定义和计算方法数值结果也存在差异这在一定程度上影响了Timoshenko梁理论的应用。1.3 定义方式的多样性剪切系数的定义方式之所以多样根源在于等效可以有不同的准则。应变能等效定义这是最常用的定义方式。假设截面上的剪应力分布为τ(y,z)则实际剪应变能为如果用等效均匀剪应力τ_eq来表示则应变能为令两者相等得到而剪切系数定义为κ τ_avg/τ_eq (V/A)/τ_eq。这种定义方式从能量角度出发物理意义明确也是Cowper采用的基本思路。变形一致性定义另一种思路是要求等效剪切应变产生的截面变形与实际变形在某种平均意义上一致。具体来说可以定义截面的平均剪切位移然后令等效模型给出相同的位移。这种定义方式得到的剪切系数可能与应变能定义不同。频率匹配定义对于动力学问题可以要求等效模型的振动频率与三维弹性理论预测的频率一致。Mindlin在1951年的工作中采用了这种思路发现剪切系数是频率的函数高频时需要修正。应力平均值定义最简单的定义是直接用最大剪应力与平均剪应力的比值。对于矩形截面最大剪应力是平均剪应力的1.5倍所以κ 1/1.5 2/3。这正是Timoshenko最初采用的方法但这种方法忽略了应力分布的形状效应。不同的定义方式得到的结果可能不同这是早期文献中剪切系数数值存在差异的主要原因。直到Cowper在1966年从三维弹性理论出发给出了一个统一且严格的推导这个问题才得到较为满意的解决。Cowper的突破性贡献2.1 从三维弹性理论出发Cowper在1966年发表的论文其核心贡献在于从三维弹性力学的基本方程出发推导出了剪切系数的一般公式。与之前学者采用的各种近似方法不同Cowper的方法是严格的不依赖于任何关于应力分布的先验假设。Cowper的出发点是三维弹性理论的平衡方程和本构关系。考虑一根等截面直梁承受横向荷载q(x)。在三维弹性理论中梁内任意一点的位移可以表示为u(x,y,z)、v(x,y,z)、w(x,y,z)分别代表沿x、y、z方向的位移分量。这里x沿梁轴线方向y和z是截面内的坐标。Cowper假设位移场具有如下形式其中w(x)是中性轴的挠度φ(x)是截面的弯曲转角f(y,z)是一个待定的截面形状函数。这个位移假设的关键在于引入了形状函数f(y,z)它描述了剪应力分布的非均匀性。通过选择合适的f(y,z)可以使位移场满足平衡方程和边界条件。2.2 推导过程的关键步骤Cowper的推导过程相当精妙这里概述其主要步骤确定形状函数f(y,z)形状函数f(y,z)需要满足截面边界上的条件。Cowper指出f(y,z)应该满足以下方程在截面边界上f的梯度应该满足其中n是边界的外法线方向。这个条件保证了截面边界上没有剪应力分量垂直于边界。实际上f(y,z)就是扭转问题中熟知的翘曲函数Warping Function。对于简单截面f(y,z)可以解析求解对于复杂截面需要数值方法。计算剪应力分布根据位移场和弹性本构关系可以得到剪应力分量这里G是剪切模量。可以看到剪应力由两部分组成一部分与(dw/dx - φ)相关代表均匀剪切另一部分与df/dy或df/dz相关代表非均匀分布的修正。计算应变能剪切的应变能为将剪应力表达式代入并积分可以得到应变能与位移变量的关系。定义剪切系数令应变能用等效剪切应变表示与第三步的结果比较可以得到剪切系数的表达式。Cowper最终给出的公式为其中I是截面惯性矩S是截面对中性轴的静矩b是截面宽度。对于非矩形截面这个公式需要适当修正。2.3 Cowper公式的物理意义Cowper公式的核心在于积分项这个积分反映了剪应力分布的非均匀程度。对于均匀剪应力分布S与y成正比积分结果为A·I²/A I²因此κ 1。对于非均匀分布积分结果大于I²因此κ 1。可以这样理解剪切系数κ反映了截面上剪应力分布的集中程度。剪应力分布越不均匀越集中在中性轴附近κ越小分布越均匀κ越接近1。这也解释了为什么薄壁截面的剪切系数较小——剪应力主要集中在腹板上。Cowper公式还揭示了泊松比ν对剪切系数的影响。对于矩形截面Cowper给出的精确解为当ν 0.3时κ 0.850当ν 0时κ 0.833。可以看到泊松比的影响并不大但对于精确计算是需要考虑的。2.4 Cowper工作的历史意义Cowper的论文发表后迅速成为该领域的经典文献被引用数千次。其重要意义在于统一了定义方式从基本力学原理出发给出了一个明确的定义消除了之前不同定义带来的混淆。提供了精确解对于常见截面形状给出了精确或高精度的数值结果可以直接用于工程计算。揭示了物理本质通过将剪切系数与截面几何特性联系起来帮助人们理解了剪切系数的物理本质——它是截面剪应力分布非均匀性的度量。建立了计算框架对于任意形状的截面提供了一个统一的计算框架可以数值求解。当然Cowper的工作也有其局限性。例如它假设材料是均匀各向同性的对于复合材料或功能梯度材料需要修正它只适用于静力或低频动力问题对于高频振动需要考虑频率依赖性。这些局限性的克服成为后续研究的重要方向。各种截面剪切系数详解3.1 矩形截面矩形截面是最简单也是最常见的截面形式。Cowper给出的精确公式为对于宽度为b、高度为h的矩形截面截面积为A bh惯性矩为I bh³/12。静矩S(y)为根据材料力学剪应力分布为将S(y)代入Cowper公式经过积分运算可以得到上述精确公式。下表给出了不同泊松比对应的剪切系数泊松比 ν剪切系数 κ与传统值2/3的偏差0.00.83325%0.20.84527%0.30.85028%0.40.85528%0.50.85729%可以看到Cowper的结果比传统值2/3大约25-30%。这意味着如果我们使用传统值会高估剪切刚度低估剪切变形。对于短粗梁这可能导致显著的误差。为什么传统值会偏低原因是Timoshenko最初采用的最大剪应力与平均剪应力之比的方法实际上给出的是剪应力分布的峰值因子而不是真正的等效系数。应变能等效原则需要考虑整个应力分布的形状而不仅仅是峰值。3.2 圆形截面对于半径为R的实心圆形截面Cowper给出的精确公式为当ν 0.3时κ 0.886当ν 0时κ 0.857。与传统值3/4即0.75相比Cowper的结果大约18%。圆形截面的推导比矩形截面复杂因为需要使用极坐标。截面上任意一点(r, θ)处的静矩为剪应力分布为将S(r)代入Cowper公式并积分可以得到上述精确公式。圆形截面的剪切系数比矩形截面略大这反映了圆形截面剪应力分布相对更均匀的特点。在圆形截面上剪应力从中心到边缘的变化比矩形截面平缓因此等效系数更接近1。3.3 空心圆形截面空心圆形截面圆管在工程中应用广泛如钢管、管道等。设外半径为R_o内半径为R_i定义半径比α R_i/R_o。Cowper给出的公式为对于薄壁圆管α接近1剪切系数趋近于当ν 0.3时κ ≈ 0.87。有趣的是薄壁圆管的剪切系数与壁厚关系不大主要取决于泊松比。对于厚壁圆管剪切系数随α减小而增大。当α 0实心圆时公式退化为实心圆的结果当α → 1薄壁时趋近于薄壁公式。3.4 椭圆截面对于半长轴为a、半短轴为b的椭圆截面Cowper给出的公式为当a b圆形时公式退化为圆形截面的结果。当a b扁平椭圆时剪切系数趋近于矩形截面的结果。椭圆截面的剪切系数介于圆形和矩形之间这符合直觉椭圆的几何特性介于两者之间。3.5 半圆形截面半圆形截面较少见但在某些特殊结构中可能出现。Cowper给出的数值结果为当ν 0时κ 0.833与矩形相同 当ν 0.3时κ 0.850半圆形截面的剪切系数与矩形截面相近这可能与半圆和矩形的高宽比类似有关。3.6 工字钢与薄壁截面工字钢、H型钢等薄壁截面是钢结构中最常用的截面形式。这类截面的剪切系数计算较为复杂因为剪应力主要集中在腹板上翼缘的贡献很小。对于标准的工字钢截面剪切系数可以近似为其中A_web是腹板面积A是总截面积。这个近似基于以下观察大部分剪力由腹板承担而翼缘对剪切刚度的贡献很小。更精确的计算需要考虑翼缘的剪切贡献。Cowper的方法可以用于任意薄壁截面但需要数值积分。下表给出了几种常见型钢的剪切系数参考值截面类型剪切系数 κ工字钢窄翼缘0.30 - 0.40工字钢宽翼缘0.35 - 0.45H型钢0.40 - 0.50槽钢0.30 - 0.40角钢0.60 - 0.70T型钢0.35 - 0.45可以看到薄壁开口截面的剪切系数明显小于实心截面这是因为剪应力高度集中在腹板或特定部位。对于闭口薄壁截面如箱形截面、空心矩形剪切系数通常比开口截面大因为剪应力可以沿整个截面周长分布。不同学者结果的对比与讨论4.1 主要学者及其贡献自Timoshenko以来许多学者对剪切系数进行了研究。下表总结了主要学者的贡献学者年份核心贡献矩形截面κTimoshenko1921提出概念给出经验值0.667Goens1931实验测量频率分析~0.85Mindlin1951高频振动频率依赖性变量Cowper1966三维弹性理论严格推导0.850Goodman Sutherland1951变分方法~0.85Roark1965工程手册实用公式0.833Kaneko1975综述对比分析0.85Hutchinson2001高精度数值解0.8501除了Timoshenko的早期经验值大多数学者的结果都在0.85左右相互吻合良好。这表明对于矩形截面剪切系数的问题已经基本解决。4.2 结果差异的原因分析为什么不同学者的结果会存在差异主要原因有定义方式不同应变能等效、变形一致性、频率匹配等不同定义可能给出略有不同的结果。Cowper采用应变能等效这是最广泛接受的定义。假设条件不同一些学者采用了不同的假设如忽略泊松效应、假设简化的应力分布等。这些假设可能导致结果偏差。应用场景不同对于静力问题Cowper的结果是精确的但对于高频振动剪切系数可能与频率相关Mindlin的工作就考虑了这一点。计算精度不同早期学者可能使用近似方法或手工计算精度有限现代数值方法可以给出高精度结果。4.3 如何选择合适的剪切系数在实际工程中如何选择剪切系数以下是一些建议对于静力分析使用Cowper公式计算或采用本文提供的参考值。对于常见截面Cowper的结果是可靠的。对于低频动力分析同样可以使用Cowper的结果。对于结构的前几阶模态剪切系数的频率依赖性可以忽略。对于高频动力分析需要考虑剪切系数的频率依赖性。Mindlin的工作表明当振动频率接近剪切波在截面内的传播频率时等效剪切刚度会发生变化。这种情况在超声波、冲击等问题中可能出现。对于复杂截面使用数值方法如有限元计算Cowper公式中的积分。大多数有限元软件都提供了计算剪切系数的功能。对于复合材料需要考虑材料的各向异性。简单的修正方法是将Cowper公式中的剪切模量G替换为等效剪切模量但更精确的方法需要专门推导。现代发展与工程应用5.1 有限元软件中的处理现代有限元软件在处理Timoshenko梁单元时都需要指定剪切系数。以下是几种常用软件的处理方式ANSYSBEAM188和BEAM189单元默认使用剪切系数。对于常见截面软件会自动计算用户也可以通过SECDATA命令手动指定。对于矩形截面ANSYS默认κ 5/6 ≈ 0.833与Cowper的结果略有差异。ABAQUSB21、B22、B31、B32等梁单元使用剪切系数。对于矩形截面默认κ 5/6用户可以通过*Transverse Shear Stiffness选项自定义。SAP2000在框架单元中考虑剪切变形用户可以指定剪切面积A_s κA。对于矩形截面建议使用κ 0.85。OpenSeeselasticBeamColumn和forceBeamColumn单元可以指定剪切刚度。用户需要计算κGA并输入。建议对于精确分析建议使用Cowper公式计算的剪切系数而不是软件默认值。不同软件的默认值可能不同了解其来源很重要。5.2 复合材料梁复合材料如碳纤维增强塑料、玻璃纤维增强塑料在航空航天、汽车等领域应用广泛。复合材料梁的剪切系数计算比各向同性材料复杂原因在于材料各向异性复合材料在不同方向具有不同的弹性模量纵向模量E₁和横向模量E₂可能相差一个数量级。层间剪切复合材料的层间剪切强度较低层间剪切变形可能占总变形的较大比例。截面非均匀复合材料梁可能采用夹层结构或铺层设计截面材料分布不均匀。对于复合材料梁剪切系数的计算需要考虑材料的各向异性。一种常用的方法是使用等效剪切模量然后应用Cowper公式。但这种方法忽略了材料分布的非均匀效应可能不够精确。更精确的方法是直接从三维弹性理论出发考虑材料的各向异性。5.3 功能梯度材料功能梯度材料Functionally Graded Materials, FGM是一种新型材料其材料属性沿某个方向连续变化。例如陶瓷-金属功能梯度材料一侧富陶瓷另一侧富金属中间区域成分梯度变化。对于功能梯度材料梁剪切系数的计算面临新的挑战材料属性是位置的函数G G(y,z)。这意味着剪应力分布不仅取决于截面几何还取决于材料梯度。近年来一些学者发展了功能梯度材料梁剪切系数的计算方法。基本思路是将Cowper的方法推广到非均匀材料情况通过数值积分计算剪切系数。5.4 高频振动问题Mindlin在1951年的经典论文中指出剪切系数可能与振动频率相关。当振动频率较高时截面内的惯性效应不能忽略等效剪切刚度会发生变化。Mindlin的分析表明剪切系数可以表示为其中κ₀是静力剪切系数ω是振动频率ω_s是与剪切波相关的特征频率f是一个无量纲函数。对于低频ω ω_sf ≈ 1剪切系数退化为静力值。对于高频f 1等效剪切刚度降低。高频剪切系数的研究主要应用于超声波、声学、冲击动力学等领域。对于常规结构动力学分析频率效应通常可以忽略。5.5 实验验证研究理论预测的剪切系数是否可靠实验验证是关键。早期的实验研究主要通过测量梁的振动频率来反推剪切系数。Goens1931使用这种方法得到了与Cowper理论相近的结果。现代实验方法包括频率测量法测量不同长度梁的振动频率通过拟合分析确定剪切系数。静态变形测量测量短粗梁在横向荷载下的挠度与理论预测比较。应变测量使用应变片测量梁表面的应变分布推断剪切变形。数字图像相关DIC现代光学测量技术可以全场测量位移和应变为剪切系数的验证提供了新的手段。总体而言实验结果与Cowper理论吻合良好验证了理论的可靠性。结论与建议6.1 核心要点总结经过六十多年的发展Timoshenko梁理论中的剪切系数问题已经基本解决。剪切系数的本质它是将非均匀剪应力分布等效为均匀分布的修正因子反映了截面上剪应力分布的非均匀程度。Cowper的贡献从三维弹性理论出发给出了剪切系数的严格推导和精确公式统一了定义方式提供了可靠的计算方法。数值规律实心截面矩形、圆形的剪切系数在0.85左右薄壁开口截面工字钢、槽钢的剪切系数较小约0.3-0.5剪切系数随泊松比略有变化。工程应用有限元软件中需要指定剪切系数建议使用Cowper公式计算而不是简单采用默认值。6.2 实践指南如何为不同截面选择剪切系数矩形截面使用κ 10(1ν)/(1211ν)当ν 0.3时κ 0.850。圆形截面使用κ 6(1ν)/(76ν)当ν 0.3时κ 0.886。薄壁工字钢近似使用κ A_web/A或参考型钢手册。复杂截面使用数值方法计算Cowper公式中的积分。软件设置注意事项了解软件默认值及其来源对于精确分析手动输入计算值注意单位制的一致性验证计算结果的合理性常见误区与纠正误区剪切系数总是等于2/3或5/6。纠正这只是特定截面的近似值不同截面差异很大。误区剪切系数对结果影响很小。纠正对于短粗梁或高阶模态影响可能很显著。误区软件默认值总是正确的。纠正默认值可能基于简化假设不一定适用于所有情况。6.3 待研究的问题尽管剪切系数的基本问题已经解决但仍有一些方向值得进一步研究复杂截面的精确解对于任意形状截面特别是薄壁开口截面如何快速准确地计算剪切系数非均匀材料复合材料、功能梯度材料的剪切系数计算方法仍需完善。非线性材料弹塑性、粘弹性材料的剪切系数如何定义和计算多场耦合在热-力耦合、流-固耦合等问题中剪切系数是否需要修正智能材料压电、形状记忆合金等智能材料梁的剪切系数如何处理这些问题不仅是学术研究的课题也具有重要的工程应用价值。参考文献[1] Timoshenko, S.P. (1921). On the correction for shear of the differential equation for transverse vibrations of prismatic bars.Philosophical Magazine, 41(245), 744-746.[2] Timoshenko, S.P. (1922). On the transverse vibrations of bars of uniform cross-section.Philosophical Magazine, 43(253), 125-131.[3] Cowper, G.R. (1966). The shear coefficient in Timoshenkos beam theory.Journal of Applied Mechanics, 33(2), 335-340.[4] Mindlin, R.D. (1951). Influence of rotatory inertia and shear on flexural motions of isotropic, elastic plates.Journal of Applied Mechanics, 18(1), 31-38.[5] Goens, E. (1931). Über die Bestimmung des Elastizitätsmoduls von Stäben mit Hilfe von Biegeschwingungen.Annalen der Physik, 402(6), 717-736.[6] Goodman, L.E., Sutherland, J.G. (1951). Discussion of Natural frequencies of continuous beams.Journal of Applied Mechanics, 18, 217-218.[7] Roark, R.J. (1965).Formulas for Stress and Strain(4th ed.). McGraw-Hill.[8] Kaneko, T. (1975). On Timoshenkos correction for shear in vibrating beams.Journal of Physics D: Applied Physics, 8(16), 1927.[9] Hutchinson, J.R. (2001). Shear coefficients for Timoshenko beam theory.Journal of Applied Mechanics, 68(1), 87-92.[10] Stephen, N.G. (1978). On the variation of Timoshenkos shear coefficient with frequency.Journal of Applied Mechanics, 45(4), 825-828.[11] Dong, S.B., Alpdogan, C., Taciroglu, E. (2010). Much ado about shear correction factors in Timoshenko beam theory.International Journal of Solids and Structures, 47(13), 1651-1665.[12] Reddy, J.N. (2006).An Introduction to the Finite Element Method(3rd ed.). McGraw-Hill.完更多精彩关注建源学堂【往期精彩】# 性能分析【JY】基于性能的抗震设计浅析一【JY】基于性能的抗震设计浅析二【JY】浅析消能附加阻尼比【JY】近断层结构设计策略分析与讨论【JY】浅析各动力求解算法及其算法数值阻尼(人工阻尼)理念【JY|体系】结构概念设计之(结构体系概念)【JY|理念】结构概念设计之(设计理念进展)【JY|减震】结构概念之(消能减震黏滞阻尼器)【JY|隔震】结构概念设计之隔震概念设计# 概念机理【JY】砌体的精细化有限元模拟【JY】从一块薄板计算说起算例对比【JY】Abaqus 三维应力单元解析、选择与应用指南【JY】基于Ramberg-Osgood本构模型的双线性计算分析【JY】结构动力学初步-单质点结构的瞬态动力学分析【JY】从一根悬臂梁说起【JY】反应谱的详解与介绍【JY】结构瑞利阻尼与经济订货模型【JY】主成分分析与振型分解【JY】浅谈结构多点激励之概念机理上【JY】浅谈结构多点激励之分析方法下【JY】板壳单元的分析详解【JY】橡胶支座的简述和其力学性能计算# 软件讨论【JY】模态分析关键点笔记【JY】浅析时程分析中的阻尼设置【JY】减隔震元件计算表格分享【JY】复合材料分析利器—内聚力单元【JY】SDOF计算教学软件开发应用分享【JY】Abaqus6.14-4如何关联fortran【JY】如何利用python来编写GUI【JY】如何解决MATLAB GUI编程软件移植运行问题
【JY】Timoshenko梁理论中的剪切系数：从经典理论到现代应用
在材料力学教材中剪切系数如矩形截面常取6/5是基于经典梁理论推导得出的用于计算剪切变形或切应力分布。然而在有限元软件如Abaqus中默认的剪切系数可能有所不同因为软件通常采用更通用的理论如Timoshenko梁理论或根据截面形状自动计算例如对于矩形截面Abaqus默认的剪切系数可能不同与教材中的值存在差异。从Timoshenko在1921年提出问题到Cowper在1966年给出严格解答再到后续学者在各种复杂情况下的推广剪切系数的研究历程折射出力学发展的一个缩影从工程直觉到严格理论从简单情况到复杂应用从单一学科到交叉融合。一篇穿越60年的学术探索带你读懂剪切系数背后的力学本质。这个表一定正确吗材料力学写在前文用有限元软件建梁单元模型时你可能注意过一个不起眼的参数——剪切系数或者叫剪切修正因子。大多数人直接用软件默认值很少去想这个系数从何而来、为什么需要它、取多少才算合理。这个看似微不足道的参数连接着梁弯曲理论中最核心的问题如何用一维梁模型去描述三维弹性体的真实变形行为从20世纪初起,当时的工程师们已经熟练运用欧拉-伯努利梁理论来设计各种结构。这个理论假设梁的横截面在弯曲后仍保持平面且垂直于中性轴也就是著名的平截面假设。对于细长梁这个假设足够精确。但随着工程实践的发展人们发现了一个问题梁比较短粗时理论预测的变形总是比实际测量值小而且频率越高、梁越短粗偏差就越明显。问题的根源在于剪切变形被忽略了。1921年Stephen Timoshenko发表了一篇开创性的论文提出了后来以他名字命名的梁理论。这个理论的核心创新在于引入了两个独立的位移变量挠度w和截面转角φ从而将剪切变形纳入考虑范围。然而引入剪切变形带来了一个新问题如何定义等效的剪切应变梁横截面上的剪应力分布是不均匀的。以矩形截面为例剪应力沿截面高度呈抛物线分布中性轴处最大上下边缘为零。这意味着截面上不同点的剪切变形是不同的。但在梁理论中我们需要用一个统一的截面剪切应变来描述整个截面的剪切行为这就需要一个修正系数将实际的非均匀剪应力分布等效为一个均匀分布。这个修正系数就是剪切系数κ。1966年G.R. Cowper发表了一篇被广泛引用的论文《The Shear Coefficient in Timoshenkos Beam Theory》从三维弹性理论出发给出了剪切系数的严格推导和精确公式。这篇论文不仅解决了理论上的争议更为工程应用提供了可靠的依据。本文将系统梳理剪切系数的理论发展、计算方法与工程应用。剪切系数的起源与定义1.1 问题的物理本质要理解剪切系数的由来需要深入理解梁弯曲时横截面上的应力分布规律。考虑一根承受横向荷载的梁任意横截面上既有弯矩M也有剪力V。弯矩产生正应力σ剪力产生剪应力τ。根据材料力学的基本分析对于矩形截面梁剪应力沿截面高度的分布遵循抛物线规律其中A是截面面积h是截面高度y是从中性轴量起的距离。可以看到中性轴处y0剪应力最大为平均剪应力的1.5倍而在上下边缘y±h/2剪应力为零。这种非均匀分布带来一个根本性的问题在Timoshenko梁理论中我们假设整个截面有一个统一的剪切变形γ这个γ等于截面转角φ对x的导数减去挠度w对x的导数。但是如果截面上各点的实际剪应力因而剪应变都不一样我们该如何定义这个统一的γ答案是我们需要某种等效处理。具体来说我们希望找到一个等效剪应力τ使得当截面上所有点都承受这个均匀剪应力时产生的应变能与实际非均匀剪应力分布产生的应变能相等。或者从另一个角度我们希望找到一个等效剪切刚度使得在相同剪力作用下梁的变形与实际三维弹性体的变形一致。这就是剪切系数κ的物理本质它是一个将非均匀剪应力分布等效为均匀分布的修正因子。数学上如果我们定义平均剪应力为τ那么等效剪应力可以表示为τκτ或者更常见的定义是等效剪切应变为γκ其中G是剪切模量。1.2 Timoshenko的原始工作1921年Timoshenko在《Philosophical Magazine》上发表了题为On the correction for shear of the differential equation for transverse vibrations of prismatic bars的论文。在这篇论文中他首次系统地将剪切变形和转动惯量纳入梁的振动方程。Timoshenko的基本思路是梁的变形由两部分组成——弯曲变形和剪切变形。弯曲变形使截面产生转动剪切变形使截面产生相对于弯曲转动方向的附加转动。因此总转角θ可以分解为两部分由于弯曲产生的转角∂w/∂x以及由于剪切产生的转角γ。在Timoshenko的原始表述中他引入了两个独立的变量挠度w(x,t)和截面转角φ(x,t)其中φ代表弯曲产生的转角而剪切应变为γ ∂w/∂x - φ。在建立运动方程时Timoshenko意识到由于剪应力分布的非均匀性直接使用V/GA来计算剪切应变是不准确的。他引入了一个修正系数k后来常用κ表示将剪切刚度修正为kGA。对于矩形截面他建议k 2/3对于圆形截面k 3/4对于薄壁工字钢截面k ≈ A_web/A腹板面积与总面积之比。这些数值是怎么来的Timoshenko并没有给出严格的推导。对于矩形截面k 2/3可以理解为最大剪应力是平均剪应力的1.5倍而等效剪应力应该介于两者之间取倒数正好是2/3。对于圆形截面k 3/4也有类似的直观解释。但这些只是经验性的估计缺乏理论基础。1922年Timoshenko发表了第二篇相关论文进一步完善了理论框架但仍然没有解决剪切系数的精确定义问题。在接下来的几十年里不同学者提出了不同的定义和计算方法数值结果也存在差异这在一定程度上影响了Timoshenko梁理论的应用。1.3 定义方式的多样性剪切系数的定义方式之所以多样根源在于等效可以有不同的准则。应变能等效定义这是最常用的定义方式。假设截面上的剪应力分布为τ(y,z)则实际剪应变能为如果用等效均匀剪应力τ_eq来表示则应变能为令两者相等得到而剪切系数定义为κ τ_avg/τ_eq (V/A)/τ_eq。这种定义方式从能量角度出发物理意义明确也是Cowper采用的基本思路。变形一致性定义另一种思路是要求等效剪切应变产生的截面变形与实际变形在某种平均意义上一致。具体来说可以定义截面的平均剪切位移然后令等效模型给出相同的位移。这种定义方式得到的剪切系数可能与应变能定义不同。频率匹配定义对于动力学问题可以要求等效模型的振动频率与三维弹性理论预测的频率一致。Mindlin在1951年的工作中采用了这种思路发现剪切系数是频率的函数高频时需要修正。应力平均值定义最简单的定义是直接用最大剪应力与平均剪应力的比值。对于矩形截面最大剪应力是平均剪应力的1.5倍所以κ 1/1.5 2/3。这正是Timoshenko最初采用的方法但这种方法忽略了应力分布的形状效应。不同的定义方式得到的结果可能不同这是早期文献中剪切系数数值存在差异的主要原因。直到Cowper在1966年从三维弹性理论出发给出了一个统一且严格的推导这个问题才得到较为满意的解决。Cowper的突破性贡献2.1 从三维弹性理论出发Cowper在1966年发表的论文其核心贡献在于从三维弹性力学的基本方程出发推导出了剪切系数的一般公式。与之前学者采用的各种近似方法不同Cowper的方法是严格的不依赖于任何关于应力分布的先验假设。Cowper的出发点是三维弹性理论的平衡方程和本构关系。考虑一根等截面直梁承受横向荷载q(x)。在三维弹性理论中梁内任意一点的位移可以表示为u(x,y,z)、v(x,y,z)、w(x,y,z)分别代表沿x、y、z方向的位移分量。这里x沿梁轴线方向y和z是截面内的坐标。Cowper假设位移场具有如下形式其中w(x)是中性轴的挠度φ(x)是截面的弯曲转角f(y,z)是一个待定的截面形状函数。这个位移假设的关键在于引入了形状函数f(y,z)它描述了剪应力分布的非均匀性。通过选择合适的f(y,z)可以使位移场满足平衡方程和边界条件。2.2 推导过程的关键步骤Cowper的推导过程相当精妙这里概述其主要步骤确定形状函数f(y,z)形状函数f(y,z)需要满足截面边界上的条件。Cowper指出f(y,z)应该满足以下方程在截面边界上f的梯度应该满足其中n是边界的外法线方向。这个条件保证了截面边界上没有剪应力分量垂直于边界。实际上f(y,z)就是扭转问题中熟知的翘曲函数Warping Function。对于简单截面f(y,z)可以解析求解对于复杂截面需要数值方法。计算剪应力分布根据位移场和弹性本构关系可以得到剪应力分量这里G是剪切模量。可以看到剪应力由两部分组成一部分与(dw/dx - φ)相关代表均匀剪切另一部分与df/dy或df/dz相关代表非均匀分布的修正。计算应变能剪切的应变能为将剪应力表达式代入并积分可以得到应变能与位移变量的关系。定义剪切系数令应变能用等效剪切应变表示与第三步的结果比较可以得到剪切系数的表达式。Cowper最终给出的公式为其中I是截面惯性矩S是截面对中性轴的静矩b是截面宽度。对于非矩形截面这个公式需要适当修正。2.3 Cowper公式的物理意义Cowper公式的核心在于积分项这个积分反映了剪应力分布的非均匀程度。对于均匀剪应力分布S与y成正比积分结果为A·I²/A I²因此κ 1。对于非均匀分布积分结果大于I²因此κ 1。可以这样理解剪切系数κ反映了截面上剪应力分布的集中程度。剪应力分布越不均匀越集中在中性轴附近κ越小分布越均匀κ越接近1。这也解释了为什么薄壁截面的剪切系数较小——剪应力主要集中在腹板上。Cowper公式还揭示了泊松比ν对剪切系数的影响。对于矩形截面Cowper给出的精确解为当ν 0.3时κ 0.850当ν 0时κ 0.833。可以看到泊松比的影响并不大但对于精确计算是需要考虑的。2.4 Cowper工作的历史意义Cowper的论文发表后迅速成为该领域的经典文献被引用数千次。其重要意义在于统一了定义方式从基本力学原理出发给出了一个明确的定义消除了之前不同定义带来的混淆。提供了精确解对于常见截面形状给出了精确或高精度的数值结果可以直接用于工程计算。揭示了物理本质通过将剪切系数与截面几何特性联系起来帮助人们理解了剪切系数的物理本质——它是截面剪应力分布非均匀性的度量。建立了计算框架对于任意形状的截面提供了一个统一的计算框架可以数值求解。当然Cowper的工作也有其局限性。例如它假设材料是均匀各向同性的对于复合材料或功能梯度材料需要修正它只适用于静力或低频动力问题对于高频振动需要考虑频率依赖性。这些局限性的克服成为后续研究的重要方向。各种截面剪切系数详解3.1 矩形截面矩形截面是最简单也是最常见的截面形式。Cowper给出的精确公式为对于宽度为b、高度为h的矩形截面截面积为A bh惯性矩为I bh³/12。静矩S(y)为根据材料力学剪应力分布为将S(y)代入Cowper公式经过积分运算可以得到上述精确公式。下表给出了不同泊松比对应的剪切系数泊松比 ν剪切系数 κ与传统值2/3的偏差0.00.83325%0.20.84527%0.30.85028%0.40.85528%0.50.85729%可以看到Cowper的结果比传统值2/3大约25-30%。这意味着如果我们使用传统值会高估剪切刚度低估剪切变形。对于短粗梁这可能导致显著的误差。为什么传统值会偏低原因是Timoshenko最初采用的最大剪应力与平均剪应力之比的方法实际上给出的是剪应力分布的峰值因子而不是真正的等效系数。应变能等效原则需要考虑整个应力分布的形状而不仅仅是峰值。3.2 圆形截面对于半径为R的实心圆形截面Cowper给出的精确公式为当ν 0.3时κ 0.886当ν 0时κ 0.857。与传统值3/4即0.75相比Cowper的结果大约18%。圆形截面的推导比矩形截面复杂因为需要使用极坐标。截面上任意一点(r, θ)处的静矩为剪应力分布为将S(r)代入Cowper公式并积分可以得到上述精确公式。圆形截面的剪切系数比矩形截面略大这反映了圆形截面剪应力分布相对更均匀的特点。在圆形截面上剪应力从中心到边缘的变化比矩形截面平缓因此等效系数更接近1。3.3 空心圆形截面空心圆形截面圆管在工程中应用广泛如钢管、管道等。设外半径为R_o内半径为R_i定义半径比α R_i/R_o。Cowper给出的公式为对于薄壁圆管α接近1剪切系数趋近于当ν 0.3时κ ≈ 0.87。有趣的是薄壁圆管的剪切系数与壁厚关系不大主要取决于泊松比。对于厚壁圆管剪切系数随α减小而增大。当α 0实心圆时公式退化为实心圆的结果当α → 1薄壁时趋近于薄壁公式。3.4 椭圆截面对于半长轴为a、半短轴为b的椭圆截面Cowper给出的公式为当a b圆形时公式退化为圆形截面的结果。当a b扁平椭圆时剪切系数趋近于矩形截面的结果。椭圆截面的剪切系数介于圆形和矩形之间这符合直觉椭圆的几何特性介于两者之间。3.5 半圆形截面半圆形截面较少见但在某些特殊结构中可能出现。Cowper给出的数值结果为当ν 0时κ 0.833与矩形相同 当ν 0.3时κ 0.850半圆形截面的剪切系数与矩形截面相近这可能与半圆和矩形的高宽比类似有关。3.6 工字钢与薄壁截面工字钢、H型钢等薄壁截面是钢结构中最常用的截面形式。这类截面的剪切系数计算较为复杂因为剪应力主要集中在腹板上翼缘的贡献很小。对于标准的工字钢截面剪切系数可以近似为其中A_web是腹板面积A是总截面积。这个近似基于以下观察大部分剪力由腹板承担而翼缘对剪切刚度的贡献很小。更精确的计算需要考虑翼缘的剪切贡献。Cowper的方法可以用于任意薄壁截面但需要数值积分。下表给出了几种常见型钢的剪切系数参考值截面类型剪切系数 κ工字钢窄翼缘0.30 - 0.40工字钢宽翼缘0.35 - 0.45H型钢0.40 - 0.50槽钢0.30 - 0.40角钢0.60 - 0.70T型钢0.35 - 0.45可以看到薄壁开口截面的剪切系数明显小于实心截面这是因为剪应力高度集中在腹板或特定部位。对于闭口薄壁截面如箱形截面、空心矩形剪切系数通常比开口截面大因为剪应力可以沿整个截面周长分布。不同学者结果的对比与讨论4.1 主要学者及其贡献自Timoshenko以来许多学者对剪切系数进行了研究。下表总结了主要学者的贡献学者年份核心贡献矩形截面κTimoshenko1921提出概念给出经验值0.667Goens1931实验测量频率分析~0.85Mindlin1951高频振动频率依赖性变量Cowper1966三维弹性理论严格推导0.850Goodman Sutherland1951变分方法~0.85Roark1965工程手册实用公式0.833Kaneko1975综述对比分析0.85Hutchinson2001高精度数值解0.8501除了Timoshenko的早期经验值大多数学者的结果都在0.85左右相互吻合良好。这表明对于矩形截面剪切系数的问题已经基本解决。4.2 结果差异的原因分析为什么不同学者的结果会存在差异主要原因有定义方式不同应变能等效、变形一致性、频率匹配等不同定义可能给出略有不同的结果。Cowper采用应变能等效这是最广泛接受的定义。假设条件不同一些学者采用了不同的假设如忽略泊松效应、假设简化的应力分布等。这些假设可能导致结果偏差。应用场景不同对于静力问题Cowper的结果是精确的但对于高频振动剪切系数可能与频率相关Mindlin的工作就考虑了这一点。计算精度不同早期学者可能使用近似方法或手工计算精度有限现代数值方法可以给出高精度结果。4.3 如何选择合适的剪切系数在实际工程中如何选择剪切系数以下是一些建议对于静力分析使用Cowper公式计算或采用本文提供的参考值。对于常见截面Cowper的结果是可靠的。对于低频动力分析同样可以使用Cowper的结果。对于结构的前几阶模态剪切系数的频率依赖性可以忽略。对于高频动力分析需要考虑剪切系数的频率依赖性。Mindlin的工作表明当振动频率接近剪切波在截面内的传播频率时等效剪切刚度会发生变化。这种情况在超声波、冲击等问题中可能出现。对于复杂截面使用数值方法如有限元计算Cowper公式中的积分。大多数有限元软件都提供了计算剪切系数的功能。对于复合材料需要考虑材料的各向异性。简单的修正方法是将Cowper公式中的剪切模量G替换为等效剪切模量但更精确的方法需要专门推导。现代发展与工程应用5.1 有限元软件中的处理现代有限元软件在处理Timoshenko梁单元时都需要指定剪切系数。以下是几种常用软件的处理方式ANSYSBEAM188和BEAM189单元默认使用剪切系数。对于常见截面软件会自动计算用户也可以通过SECDATA命令手动指定。对于矩形截面ANSYS默认κ 5/6 ≈ 0.833与Cowper的结果略有差异。ABAQUSB21、B22、B31、B32等梁单元使用剪切系数。对于矩形截面默认κ 5/6用户可以通过*Transverse Shear Stiffness选项自定义。SAP2000在框架单元中考虑剪切变形用户可以指定剪切面积A_s κA。对于矩形截面建议使用κ 0.85。OpenSeeselasticBeamColumn和forceBeamColumn单元可以指定剪切刚度。用户需要计算κGA并输入。建议对于精确分析建议使用Cowper公式计算的剪切系数而不是软件默认值。不同软件的默认值可能不同了解其来源很重要。5.2 复合材料梁复合材料如碳纤维增强塑料、玻璃纤维增强塑料在航空航天、汽车等领域应用广泛。复合材料梁的剪切系数计算比各向同性材料复杂原因在于材料各向异性复合材料在不同方向具有不同的弹性模量纵向模量E₁和横向模量E₂可能相差一个数量级。层间剪切复合材料的层间剪切强度较低层间剪切变形可能占总变形的较大比例。截面非均匀复合材料梁可能采用夹层结构或铺层设计截面材料分布不均匀。对于复合材料梁剪切系数的计算需要考虑材料的各向异性。一种常用的方法是使用等效剪切模量然后应用Cowper公式。但这种方法忽略了材料分布的非均匀效应可能不够精确。更精确的方法是直接从三维弹性理论出发考虑材料的各向异性。5.3 功能梯度材料功能梯度材料Functionally Graded Materials, FGM是一种新型材料其材料属性沿某个方向连续变化。例如陶瓷-金属功能梯度材料一侧富陶瓷另一侧富金属中间区域成分梯度变化。对于功能梯度材料梁剪切系数的计算面临新的挑战材料属性是位置的函数G G(y,z)。这意味着剪应力分布不仅取决于截面几何还取决于材料梯度。近年来一些学者发展了功能梯度材料梁剪切系数的计算方法。基本思路是将Cowper的方法推广到非均匀材料情况通过数值积分计算剪切系数。5.4 高频振动问题Mindlin在1951年的经典论文中指出剪切系数可能与振动频率相关。当振动频率较高时截面内的惯性效应不能忽略等效剪切刚度会发生变化。Mindlin的分析表明剪切系数可以表示为其中κ₀是静力剪切系数ω是振动频率ω_s是与剪切波相关的特征频率f是一个无量纲函数。对于低频ω ω_sf ≈ 1剪切系数退化为静力值。对于高频f 1等效剪切刚度降低。高频剪切系数的研究主要应用于超声波、声学、冲击动力学等领域。对于常规结构动力学分析频率效应通常可以忽略。5.5 实验验证研究理论预测的剪切系数是否可靠实验验证是关键。早期的实验研究主要通过测量梁的振动频率来反推剪切系数。Goens1931使用这种方法得到了与Cowper理论相近的结果。现代实验方法包括频率测量法测量不同长度梁的振动频率通过拟合分析确定剪切系数。静态变形测量测量短粗梁在横向荷载下的挠度与理论预测比较。应变测量使用应变片测量梁表面的应变分布推断剪切变形。数字图像相关DIC现代光学测量技术可以全场测量位移和应变为剪切系数的验证提供了新的手段。总体而言实验结果与Cowper理论吻合良好验证了理论的可靠性。结论与建议6.1 核心要点总结经过六十多年的发展Timoshenko梁理论中的剪切系数问题已经基本解决。剪切系数的本质它是将非均匀剪应力分布等效为均匀分布的修正因子反映了截面上剪应力分布的非均匀程度。Cowper的贡献从三维弹性理论出发给出了剪切系数的严格推导和精确公式统一了定义方式提供了可靠的计算方法。数值规律实心截面矩形、圆形的剪切系数在0.85左右薄壁开口截面工字钢、槽钢的剪切系数较小约0.3-0.5剪切系数随泊松比略有变化。工程应用有限元软件中需要指定剪切系数建议使用Cowper公式计算而不是简单采用默认值。6.2 实践指南如何为不同截面选择剪切系数矩形截面使用κ 10(1ν)/(1211ν)当ν 0.3时κ 0.850。圆形截面使用κ 6(1ν)/(76ν)当ν 0.3时κ 0.886。薄壁工字钢近似使用κ A_web/A或参考型钢手册。复杂截面使用数值方法计算Cowper公式中的积分。软件设置注意事项了解软件默认值及其来源对于精确分析手动输入计算值注意单位制的一致性验证计算结果的合理性常见误区与纠正误区剪切系数总是等于2/3或5/6。纠正这只是特定截面的近似值不同截面差异很大。误区剪切系数对结果影响很小。纠正对于短粗梁或高阶模态影响可能很显著。误区软件默认值总是正确的。纠正默认值可能基于简化假设不一定适用于所有情况。6.3 待研究的问题尽管剪切系数的基本问题已经解决但仍有一些方向值得进一步研究复杂截面的精确解对于任意形状截面特别是薄壁开口截面如何快速准确地计算剪切系数非均匀材料复合材料、功能梯度材料的剪切系数计算方法仍需完善。非线性材料弹塑性、粘弹性材料的剪切系数如何定义和计算多场耦合在热-力耦合、流-固耦合等问题中剪切系数是否需要修正智能材料压电、形状记忆合金等智能材料梁的剪切系数如何处理这些问题不仅是学术研究的课题也具有重要的工程应用价值。参考文献[1] Timoshenko, S.P. (1921). On the correction for shear of the differential equation for transverse vibrations of prismatic bars.Philosophical Magazine, 41(245), 744-746.[2] Timoshenko, S.P. (1922). On the transverse vibrations of bars of uniform cross-section.Philosophical Magazine, 43(253), 125-131.[3] Cowper, G.R. (1966). The shear coefficient in Timoshenkos beam theory.Journal of Applied Mechanics, 33(2), 335-340.[4] Mindlin, R.D. (1951). Influence of rotatory inertia and shear on flexural motions of isotropic, elastic plates.Journal of Applied Mechanics, 18(1), 31-38.[5] Goens, E. (1931). Über die Bestimmung des Elastizitätsmoduls von Stäben mit Hilfe von Biegeschwingungen.Annalen der Physik, 402(6), 717-736.[6] Goodman, L.E., Sutherland, J.G. (1951). Discussion of Natural frequencies of continuous beams.Journal of Applied Mechanics, 18, 217-218.[7] Roark, R.J. (1965).Formulas for Stress and Strain(4th ed.). McGraw-Hill.[8] Kaneko, T. (1975). On Timoshenkos correction for shear in vibrating beams.Journal of Physics D: Applied Physics, 8(16), 1927.[9] Hutchinson, J.R. (2001). Shear coefficients for Timoshenko beam theory.Journal of Applied Mechanics, 68(1), 87-92.[10] Stephen, N.G. (1978). On the variation of Timoshenkos shear coefficient with frequency.Journal of Applied Mechanics, 45(4), 825-828.[11] Dong, S.B., Alpdogan, C., Taciroglu, E. (2010). Much ado about shear correction factors in Timoshenko beam theory.International Journal of Solids and Structures, 47(13), 1651-1665.[12] Reddy, J.N. (2006).An Introduction to the Finite Element Method(3rd ed.). McGraw-Hill.完更多精彩关注建源学堂【往期精彩】# 性能分析【JY】基于性能的抗震设计浅析一【JY】基于性能的抗震设计浅析二【JY】浅析消能附加阻尼比【JY】近断层结构设计策略分析与讨论【JY】浅析各动力求解算法及其算法数值阻尼(人工阻尼)理念【JY|体系】结构概念设计之(结构体系概念)【JY|理念】结构概念设计之(设计理念进展)【JY|减震】结构概念之(消能减震黏滞阻尼器)【JY|隔震】结构概念设计之隔震概念设计# 概念机理【JY】砌体的精细化有限元模拟【JY】从一块薄板计算说起算例对比【JY】Abaqus 三维应力单元解析、选择与应用指南【JY】基于Ramberg-Osgood本构模型的双线性计算分析【JY】结构动力学初步-单质点结构的瞬态动力学分析【JY】从一根悬臂梁说起【JY】反应谱的详解与介绍【JY】结构瑞利阻尼与经济订货模型【JY】主成分分析与振型分解【JY】浅谈结构多点激励之概念机理上【JY】浅谈结构多点激励之分析方法下【JY】板壳单元的分析详解【JY】橡胶支座的简述和其力学性能计算# 软件讨论【JY】模态分析关键点笔记【JY】浅析时程分析中的阻尼设置【JY】减隔震元件计算表格分享【JY】复合材料分析利器—内聚力单元【JY】SDOF计算教学软件开发应用分享【JY】Abaqus6.14-4如何关联fortran【JY】如何利用python来编写GUI【JY】如何解决MATLAB GUI编程软件移植运行问题
