四边支承矩形薄板自振频率计算

四边支承矩形薄板自振频率计算1. 基本假定及振动微分方程弹性板是假定其厚度远小于其他两尺寸的板，且材料假设为各向同性。

板的振动理论是以以下几个假定为基础的：1）板中原来在中面法线上的各点，在板弯曲变形后仍在中面的法线上。

这个假设称为直法线假设，表示横向剪切变形忽略不计。

2）板的挠度比板厚小很多，板弯曲时中面不产生变形，即中面为中性面。

3）板的横向正应力与其他两个方向正应力相比较，可以忽略不计。

在此基础上，若假定板的挠度不从平面位置算起，而从平衡位置算起，对板内平行六面体进行微元分析，由平衡条件、变形协调条件和物理方程得板的弯曲平衡方程式，然后分析板在振动过程中的动力平衡，可得板的自由振动微分方程[1]：022********=∂∂+∂∂+∂∂+∂∂twm y x w D y w D x w D （1） 等式中)1(1223ν-=Eh D ，式中： m 为板的单位面积的质量；D 为板的弯曲刚度,E ,ν分别为板的弹性模量和泊松比，h 为板的厚度。

微分方程(1)的解答形式为薄板上每一点),(y x 的挠度),()sin cos (1y x W t B t A w m m m m m m ωω+=∑∞=。

被表示成无数多个简谐振动下的挠度相叠加，而每一个简谐振动的圆频率是m ω。

另一方面，薄板在每一瞬时t 的挠度，则表示成为无数多种振形下的挠度相叠加，而每一种振形下的挠度是由振形函数),(y x W m 表示的，为求出各种振形下的振形函数m W ，以及与之相应的圆频率m ω，我们取),()sin cos (y x W t B t A w ωω+=代入方程（1）消除因子)sin cos (t B t A ωω+得到振形微分方程：0222244444=-∂∂+∂∂+∂∂W m yx WD y W D x W D ω (2) 2. 边界条件振形函数需要满足各边界条件，板的边界一般有固支边，简支边，自由边三种情况，这里以x=0的边为例，其相应的边界条件为：固定边：沿固定边的位移和转角为0，即0)(0==x W ，0)(0=∂∂=x xW； 简支边：沿简支边的位移和弯矩为0，即0)(0==x W ，0)(022=∂∂=x xW；自由边：沿自由边的弯矩和剪力为0，即0)(02222=∂∂+∂∂=x y W x W ν，0))2((02333=∂∂∂-+∂∂=x yx Wx W ν 对于四边支承板有如下6中不同边界条件：（a ） （b ）（c ） （d ）（e ） （f ）一般而言，假定合适的位移函数，利用边界条件可以求解上述微分方程。

四边简支矩形板计算项目名称_____________日期_____________设计者_____________校对者_____________一、构件编号: LB-1二、示意图三、依据规范《建筑结构荷载规范》 GB50009-2001《混凝土结构设计规范》 GB50010-2010四、计算信息1.几何参数计算跨度: Lx = 11000 mm; Ly = 7500 mm板厚: h = 400 mm2.材料信息混凝土等级: C30 fc=mm2 ft=mm2 ftk=mm2Ec=×104N/mm2钢筋种类: HRB400 fy = 360 N/mm2Es = ×105 N/mm2最小配筋率: ρ= %纵向受拉钢筋合力点至近边距离: as = 55mm保护层厚度: c = 40mm3.荷载信息(均布荷载)永久荷载分项系数: γG =可变荷载分项系数: γQ =准永久值系数: ψq =永久荷载标准值: ｑgk = m2可变荷载标准值: ｑqk = m24.计算方法:弹性板5.边界条件(上端/下端/左端/右端):简支/简支/简支/简支6.设计参数结构重要性系数: γo =泊松比:μ =五、计算参数:1.计算板的跨度: Lo = 7500 mm2.计算板的有效高度: ho = h-as=400-55=345 mm六、配筋计算(lx/ly=11000/7500=< 所以按双向板计算):向底板钢筋1) 确定X向板底弯矩Mx = 表中系数(γG*ｑgk+γQ*ｑqk)*Lo2= +***+**= kN*m2) 确定计算系数αs = γo*Mx/(α1*fc*b*ho*ho)= *×106/**1000*345*345)=3) 计算相对受压区高度ξ = 1-sqrt(1-2*αs) = 1-sqrt(1-2* =4) 计算受拉钢筋面积As = α1*fc*b*ho*ξ/fy = **1000*345*360= 354mm25) 验算最小配筋率ρ = As/(b*h) = 354/(1000*400) = %ρ<ρmin = % 不满足最小配筋要求所以取面积为As = ρmin*b*h = %*1000*400 = 800 mm2采取方案⌲12@140, 实配面积807 mm2向底板钢筋1) 确定Y向板底弯矩My = 表中系数(γG*ｑgk+γQ*ｑqk)*Lo2= +***+**= kN*m2) 确定计算系数αs = γo*My/(α1*fc*b*ho*ho)= *×106/**1000*345*345)=3) 计算相对受压区高度ξ = 1-sqrt(1-2*αs) = 1-sqrt(1-2* =4) 计算受拉钢筋面积As = α1*fc*b*ho*ξ/fy = **1000*345*360= 638mm25) 验算最小配筋率ρ = As/(b*h) = 638/(1000*400) = %ρ<ρmin = % 不满足最小配筋要求所以取面积为As = ρmin*b*h = %*1000*400 = 800 mm2采取方案⌲12@100, 实配面积1131 mm2七、跨中挠度计算：Mk -------- 按荷载效应的标准组合计算的弯矩值Mq -------- 按荷载效应的准永久组合计算的弯矩值1.计算荷载效应Mk = Mgk + Mqk= +**+* = kN*mMq = Mgk+ψq*Mqk= +**+** = kN*m2.计算受弯构件的短期刚度 Bs1) 计算按荷载荷载效应的两种组合作用下，构件纵向受拉钢筋应力σsk = Mk/*ho*As) 混规(7.1.4-3)= ×106/*345*1131) = N/mmσsq = Mq/*ho*As) 混规(7.1.4-3)= ×106/*345*1131) = N/mm2) 计算按有效受拉混凝土截面面积计算的纵向受拉钢筋配筋率矩形截面积: Ate = *b*h = *1000*400= 200000mm2ρte = As/Ate 混规(7.1.2-4)= 1131/200000 = %3) 计算裂缝间纵向受拉钢筋应变不均匀系数ψψk = 混规(7.1.2-2)= =因为ψ不能小于最小值，所以取ψk =ψq = 混规(7.1.2-2)= =因为ψ不能小于最小值，所以取ψq =4) 计算钢筋弹性模量与混凝土模量的比值αEαE = Es/Ec = ×105/×104 =5) 计算受压翼缘面积与腹板有效面积的比值γf矩形截面，γf=06) 计算纵向受拉钢筋配筋率ρρ = As/(b*ho)= 1131/(1000*345) = %7) 计算受弯构件的短期刚度 BsBsk = Es*As*ho2/[ψk++6*αE*ρ/(1+ γf')](混规(7.2.3-1)) = ×105*1131*3452/[*++6**%/(1+*]= ×104 kN*m2Bsq = Es*As*ho2/[ψq++6*αE*ρ/(1+ γf')](混规(7.2.3-1)) = ×105*1131*3452/[*++6**%/(1+*]= ×104 kN*m23.计算受弯构件的长期刚度B1) 确定考虑荷载长期效应组合对挠度影响增大影响系数θ当ρ'=0时，θ= 混规(7.2.5)2) 计算受弯构件的长期刚度 BBk = Mk/(Mq*(θ-1)+Mk)*Bs (混规(7.2.2-1))= *+*×104= ×104 kN*m2Bq = Bsq/θ (混规(7.2.2-2))= ×104/= ×104 kN*m2B = min(Bk,Bq)= min,=4.计算受弯构件挠度f max = f*(q gk+q qk)*Lo4/B= *+*×104=5.验算挠度挠度限值fo=Lo/250=7500/250=fmax=≤fo=，满足规范要求!八、裂缝宽度验算:1.跨中X方向裂缝1) 计算荷载效应Mx = 表中系数(ｑgk+ψｑqk)*Lo2= +**+**= kN*m2) 带肋钢筋,所以取值v i=3) 因为C > 65，所以取C = 654) 计算按荷载效应的准永久组合作用下，构件纵向受拉钢筋应力σsq=Mq/*ho*As) 混规(7.1.4-3)=×106/*345*807)=mm5) 计算按有效受拉混凝土截面面积计算的纵向受拉钢筋配筋率矩形截面积，Ate=*b*h=*1000*400=200000 mm2ρte=As/Ate 混规(7.1.2-4)=807/200000 =因为ρte= < ,所以让ρte=6) 计算裂缝间纵向受拉钢筋应变不均匀系数ψψ= 混规(7.1.2-2)= =7) 计算单位面积钢筋根数nn=1000/dist = 1000/140=78) 计算受拉区纵向钢筋的等效直径d eqd eq= (∑n i*d i2)/(∑n i*v i*d i)=7*12*12/(7**12)=129) 计算最大裂缝宽度ωmax=αcr*ψ*σsq/Es**C+*Deq/ρte) (混规(7.1.2-1)=**×105**40+*12/= ≤ , 满足规范要求2.跨中Y方向裂缝1) 计算荷载效应My = 表中系数(ｑgk+ψｑqk)*Lo2= +**+**= kN*m2) 带肋钢筋,所以取值v i=3) 因为C > 65，所以取C = 654) 计算按荷载效应的准永久组合作用下，构件纵向受拉钢筋应力σsq=Mq/*ho*As) 混规(7.1.4-3)=×106/*345*1131)=mm5) 计算按有效受拉混凝土截面面积计算的纵向受拉钢筋配筋率矩形截面积，Ate=*b*h=*1000*400=200000 mm2ρte=As/Ate 混规(7.1.2-4)=1131/200000 =因为ρte= < ,所以让ρte=6) 计算裂缝间纵向受拉钢筋应变不均匀系数ψψ= 混规(7.1.2-2)= =7) 计算单位面积钢筋根数nn=1000/dist = 1000/100=108) 计算受拉区纵向钢筋的等效直径d eqd eq= (∑n i*d i2)/(∑n i*v i*d i)=10*12*12/(10**12)=129) 计算最大裂缝宽度ωmax=αcr*ψ*σsq/Es**C+*Deq/ρte) (混规(7.1.2-1)=**×105**40+*12/= ≤ , 满足规范要求。

均匀热环境下四边固支矩形PCB薄板的自由振动高军;黄再兴【摘要】表面贴装形式中PCB板可简化为四边固支矩形薄板.基于刚性板的小挠度理论,推导了热载下四边固支矩形PCB薄板的自由振动微分方程.从微分方程中得出,热载下的PCB薄板等效于面内受均布张力的薄板,进而通过结构力学方法将热载下四边固支薄板振动问题转换为受面内均布张力固支薄板振动问题.利用虚位移理论,得出了温度沿厚度均匀线性变化的热载下四边固支矩形PCB薄板固有频率和自由振动的挠度值的计算方法.讨论了热载下温度、薄板的几何尺寸对矩形PCB薄板自由振动固有频率的影响.结论可为矩形PCB薄板在热载下的振动分析以及固有频率计算提供方法上的参考.【期刊名称】《振动与冲击》【年(卷),期】2014(033)012【总页数】5页(P75-79)【关键词】PCB矩形薄板;热环境;四边固支;微分方程;固有频率【作者】高军;黄再兴【作者单位】南京航空航天大学机械结构力学及控制国家重点实验室,南京210016;南京航空航天大学机械结构力学及控制国家重点实验室,南京210016【正文语种】中文【中图分类】O343表面贴装技术(SMT)以其成本低、集成度高、电子组件重量轻、易于自动化等优点广泛应用于微电子电路[1-2]。

影响表面贴装电子产品可靠性的主要环境因素是热和振动冲击,特别是在环境振动和热载荷的复杂环境下，两类载荷共同影响贴装形式元器件的内力情况，导致振动产生的动态应力和热疲劳应力相互叠加引起封装的失效，从而影响整个封装形式可靠性与寿命。

同时，这两种载荷相互间产生耦合，并非仅仅只表现为两种载荷作用的简单叠加。

目前，已发现大型的工作站随工作温度升高到一定程度会产生共振，从而会影响其正常工作。

明显地，这是由热环境温度的变化导致封装结构固有频率改变带来的问题。

该问题涉及封装结构固有频率与环境温度的相互耦合，但目前还缺乏定量的研究。

已有学者分别对振动和热环境下表面贴装形式电子元器件的结构和可靠性进行了一些研究[3-6]。

四边简支矩形板计算项目名称_____________日期_____________设计者_____________校对者_____________一、构件编号: LB-1二、示意图三、依据规范《建筑结构荷载规范》 GB50009-2001《混凝土结构设计规范》 GB50010-2010四、计算信息1.几何参数计算跨度: Lx = 11000 mm; Ly = 7500 mm板厚: h = 400 mm2.材料信息混凝土等级: C30 fc=14.3N/mm2 ft=1.43N/mm2 ftk=2.01N/mm2Ec=3.00×104N/mm2钢筋种类: HRB400 fy = 360 N/mm2Es = 2.0×105 N/mm2最小配筋率: ρ= 0.200%纵向受拉钢筋合力点至近边距离: as = 55mm保护层厚度: c = 40mm3.荷载信息(均布荷载)永久荷载分项系数: γG = 1.200可变荷载分项系数: γQ = 1.400准永久值系数: ψq = 1.000永久荷载标准值: ｑgk = 15.000kN/m2可变荷载标准值: ｑqk = 0.000kN/m24.计算方法:弹性板5.边界条件(上端/下端/左端/右端):简支/简支/简支/简支6.设计参数结构重要性系数: γo = 1.00泊松比:μ = 0.200五、计算参数:1.计算板的跨度: Lo = 7500 mm2.计算板的有效高度: ho = h-as=400-55=345 mm六、配筋计算(lx/ly=11000/7500=1.467<2.000 所以按双向板计算):1.X向底板钢筋1) 确定X向板底弯矩Mx = 表中系数(γG*ｑgk+γQ*ｑqk)*Lo2= (0.0287+0.0707*0.200)*(1.200*15.000+1.400*0.000)*7.52 = 43.374 kN*m2) 确定计算系数αs = γo*Mx/(α1*fc*b*ho*ho)= 1.00*43.374×106/(1.00*14.3*1000*345*345)= 0.0253) 计算相对受压区高度ξ = 1-sqrt(1-2*αs) = 1-sqrt(1-2*0.025) = 0.0264) 计算受拉钢筋面积As = α1*fc*b*ho*ξ/fy = 1.000*14.3*1000*345*0.026/360= 354mm25) 验算最小配筋率ρ = As/(b*h) = 354/(1000*400) = 0.088%ρ<ρmin = 0.200% 不满足最小配筋要求所以取面积为As = ρmin*b*h = 0.200%*1000*400 = 800 mm2采取方案⌲12@140, 实配面积807 mm22.Y向底板钢筋1) 确定Y向板底弯矩My = 表中系数(γG*ｑgk+γQ*ｑqk)*Lo2= (0.0707+0.0287*0.200)*(1.200*15.000+1.400*0.000)*7.52 = 77.430 kN*m2) 确定计算系数αs = γo*My/(α1*fc*b*ho*ho)= 1.00*77.430×106/(1.00*14.3*1000*345*345)= 0.0453) 计算相对受压区高度ξ = 1-sqrt(1-2*αs) = 1-sqrt(1-2*0.045) = 0.0474) 计算受拉钢筋面积As = α1*fc*b*ho*ξ/fy = 1.000*14.3*1000*345*0.047/360= 638mm25) 验算最小配筋率ρ = As/(b*h) = 638/(1000*400) = 0.160%ρ<ρmin = 0.200% 不满足最小配筋要求所以取面积为As = ρmin*b*h = 0.200%*1000*400 = 800 mm2采取方案⌲12@100, 实配面积1131 mm2七、跨中挠度计算：Mk -------- 按荷载效应的标准组合计算的弯矩值Mq -------- 按荷载效应的准永久组合计算的弯矩值1.计算荷载效应Mk = Mgk + Mqk= (0.0707+0.0287*0.200)*(15.000+0.000)*7.52 = 64.525 kN*mMq = Mgk+ψq*Mqk= (0.0707+0.0287*0.200)*(15.000+1.0*0.000)*7.52 = 64.525 kN*m2.计算受弯构件的短期刚度 Bs1) 计算按荷载荷载效应的两种组合作用下，构件纵向受拉钢筋应力σsk = Mk/(0.87*ho*As) 混规(7.1.4-3)= 64.525×106/(0.87*345*1131) = 190.077 N/mmσsq = Mq/(0.87*ho*As) 混规(7.1.4-3)= 64.525×106/(0.87*345*1131) = 190.077 N/mm2) 计算按有效受拉混凝土截面面积计算的纵向受拉钢筋配筋率矩形截面积: Ate = 0.5*b*h = 0.5*1000*400= 200000mm2ρte = As/Ate 混规(7.1.2-4)= 1131/200000 = 0.566%3) 计算裂缝间纵向受拉钢筋应变不均匀系数ψψk = 1.1-0.65*ftk/(ρte*σsk) 混规(7.1.2-2)= 1.1-0.65*2.01/(0.566%*190.077) = -0.115因为ψ不能小于最小值0.2，所以取ψk = 0.2ψq = 1.1-0.65*ftk/(ρte*σsq) 混规(7.1.2-2)= 1.1-0.65*2.01/(0.566%*190.077) = -0.115因为ψ不能小于最小值0.2，所以取ψq = 0.24) 计算钢筋弹性模量与混凝土模量的比值αEαE = Es/Ec = 2.0×105/3.00×104 = 6.6675) 计算受压翼缘面积与腹板有效面积的比值γf矩形截面，γf=06) 计算纵向受拉钢筋配筋率ρρ = As/(b*ho)= 1131/(1000*345) = 0.328%7) 计算受弯构件的短期刚度 BsBsk = Es*As*ho2/[1.15ψk+0.2+6*αE*ρ/(1+ 3.5γf')](混规(7.2.3-1)) = 2.0×105*1131*3452/[1.15*-0.115+0.2+6*6.667*0.328%/(1+3.5*0.0)]= 4.798×104 kN*m2Bsq = Es*As*ho2/[1.15ψq+0.2+6*αE*ρ/(1+ 3.5γf')](混规(7.2.3-1)) = 2.0×105*1131*3452/[1.15*-0.115+0.2+6*6.667*0.328%/(1+3.5*0.0)]= 4.798×104 kN*m23.计算受弯构件的长期刚度B1) 确定考虑荷载长期效应组合对挠度影响增大影响系数θ当ρ'=0时，θ=2.0 混规(7.2.5)2) 计算受弯构件的长期刚度 BBk = Mk/(Mq*(θ-1)+Mk)*Bs (混规(7.2.2-1))= 64.525/(64.525*(2.0-1)+64.525)*4.798×104= 2.399×104 kN*m2Bq = Bsq/θ (混规(7.2.2-2))= 4.798×104/2.0= 2.399×104 kN*m2B = min(Bk,Bq)= min(23990.371,23990.371)= 23990.3714.计算受弯构件挠度f max = f*(q gk+q qk)*Lo4/B= 0.00752*(15.000+0.000)*7.54/2.399×104= 14.879mm5.验算挠度挠度限值fo=Lo/250=7500/250=30.000mmfmax=14.879mm≤fo=30.000mm，满足规范要求!八、裂缝宽度验算:1.跨中X方向裂缝1) 计算荷载效应Mx = 表中系数(ｑgk+ψｑqk)*Lo2= (0.0287+0.0707*0.200)*(15.000+1.00*0.000)*7.52= 36.145 kN*m2) 带肋钢筋,所以取值v i=1.03) 因为C > 65，所以取C = 654) 计算按荷载效应的准永久组合作用下，构件纵向受拉钢筋应力σsq=Mq/(0.87*ho*As) 混规(7.1.4-3)=36.145×106/(0.87*345*807)=149.222N/mm5) 计算按有效受拉混凝土截面面积计算的纵向受拉钢筋配筋率矩形截面积，Ate=0.5*b*h=0.5*1000*400=200000 mm2ρte=As/Ate 混规(7.1.2-4)=807/200000 = 0.0040因为ρte=0.0040 < 0.01,所以让ρte=0.016) 计算裂缝间纵向受拉钢筋应变不均匀系数ψψ=1.1-0.65*ftk/(ρte*σsq) 混规(7.1.2-2)=1.1-0.65*2.010/(0.0100*149.222)=0.2247) 计算单位面积钢筋根数nn=1000/dist = 1000/140=78) 计算受拉区纵向钢筋的等效直径d eqd eq= (∑n i*d i2)/(∑n i*v i*d i)=7*12*12/(7*1.0*12)=129) 计算最大裂缝宽度ωmax=αcr*ψ*σsq/Es*(1.9*C+0.08*Deq/ρte) (混规(7.1.2-1) =1.9*0.224*149.222/2.0×105*(1.9*40+0.08*12/0.0100)=0.0547mm ≤ 0.20, 满足规范要求2.跨中Y方向裂缝1) 计算荷载效应My = 表中系数(ｑgk+ψｑqk)*Lo2= (0.0707+0.0287*0.200)*(15.000+1.00*0.000)*7.52= 64.525 kN*m2) 带肋钢筋,所以取值v i=1.03) 因为C > 65，所以取C = 654) 计算按荷载效应的准永久组合作用下，构件纵向受拉钢筋应力σsq=Mq/(0.87*ho*As) 混规(7.1.4-3)=64.525×106/(0.87*345*1131)=190.077N/mm5) 计算按有效受拉混凝土截面面积计算的纵向受拉钢筋配筋率矩形截面积，Ate=0.5*b*h=0.5*1000*400=200000 mm2ρte=As/Ate 混规(7.1.2-4)=1131/200000 = 0.0057因为ρte=0.0057 < 0.01,所以让ρte=0.016) 计算裂缝间纵向受拉钢筋应变不均匀系数ψψ=1.1-0.65*ftk/(ρte*σsq) 混规(7.1.2-2)=1.1-0.65*2.010/(0.0100*190.077)=0.4137) 计算单位面积钢筋根数nn=1000/dist = 1000/100=108) 计算受拉区纵向钢筋的等效直径d eqd eq= (∑n i*d i2)/(∑n i*v i*d i)=10*12*12/(10*1.0*12)=129) 计算最大裂缝宽度ωmax=αcr*ψ*σsq/Es*(1.9*C+0.08*Deq/ρte) (混规(7.1.2-1) =1.9*0.413*190.077/2.0×105*(1.9*40+0.08*12/0.0100)=0.1282mm ≤ 0.20, 满足规范要求。

四边支承矩形薄板自振频率计算
常为华
【期刊名称】《山西建筑》
【年(卷),期】2012(038)005
【摘要】从弹性薄板的基本振动微分方程出发,利用可以满足各边边界条件的振型函数,根据能量法推导出四边支承六种不同边界条件的矩形薄板的最低自振频率计算公式,为楼板结构的自振频率计算和舒适性设计提供了指导。

【总页数】3页(P63-65)
【作者】常为华
【作者单位】北京市建筑设计研究院,北京100000
【正文语种】中文
【中图分类】TU311.3
【相关文献】
1.具有焊接残余应力的矩形薄板固有频率计算方法研究 [J], 高永毅;唐果;万文
2.四角点支承四边自由矩形板自振分析新方法 [J], 许琪楼
3.四边固定支承矩形薄板振动分析的有限积分变换法 [J], 钟阳;张永山
4.四角点支承四边自由矩形薄板屈曲问题的新解析解 [J], 杨雨诗;安东琦;倪卓凡;李锐
5.四边任意支承条件下弹性矩形薄板弯曲问题的解析解 [J], 钟阳;张永山
因版权原因，仅展示原文概要，查看原文内容请购买。

四边支承矩形薄板自振频率计算四边支承矩形薄板的自振频率是指薄板在四个边界被支承的情况下，能够在固有模态下以多少频率振动。

这在很多工程和物理问题中都非常重要，因为它涉及到材料和结构的固有特性。

以下将详细介绍如何计算四边支承矩形薄板的自振频率。

首先，我们需要了解薄板的振动方程。

对于四边支承矩形薄板来说，其振动方程为二维拉普拉斯方程：∇^2u+k^2u=0其中，u是振幅，∇^2是二维拉普拉斯算子，k是波数，k=2πf/c，f为频率，c为波速。

接下来，我们需要根据边界条件来确定薄板的固有频率，边界条件一般可以是位移边界条件、速度边界条件或应力边界条件。

在四边支承的情况下，我们常常使用位移边界条件。

对于四边支承的矩形薄板，位移边界条件可以表示为：u(0,y)=u(a,y)=0u(x,0)=u(x,b)=0其中，(0,y)和(a,y)表示薄板的两个平行边界，(x,0)和(x,b)表示薄板的两个垂直边界。

这些边界条件表示，在边界上薄板的位移为零，即薄板被四边支撑。

这些边界条件可以用来解二维拉普拉斯方程。

接下来，在振动方程中代入位移边界条件，我们可以得到一个特征值问题。

通过求解特征值问题，我们可以得到薄板的固有频率和对应的振型。

具体来说，我们需要通过使用分离变量法，将二维拉普拉斯方程转化为两个一维波动方程。

然后，我们可以根据一维波动方程的边值条件来解特征值问题。

解特征值问题的方法有很多种，常见的包括解析解法和数值解法。

解析解法适用于一些简单的情况，如正方形或矩形薄板。

对于复杂的几何形状或边界条件，数值解法（如有限元法或边界元法）可能更合适。

在使用数值解法时，我们需要将薄板分割成小的单元，并在每个单元上使用适当的数学模型和数值方法。

然后，我们可以通过迭代计算来获得薄板的固有频率。

在实际计算中，我们还需要确定薄板的材料参数，如杨氏模量、泊松比和密度。

这些材料参数可以通过实验测试获得，或者根据已有的文献和标准进行估算。

矩形薄板四边固支时的自由振动分析张强【摘要】利用Matlab软件编写的程序对矩形薄板在四边固支情况时进行了有限元分析,求得矩形薄板在四边固支时自由振动的固有频率,并阐述了编写程序的基本思路,利用该算法求解的基频与利用能量法求解的基频相当吻合,同时利用该程序求解了第二阶到第四阶频率,将其与ANSYS求解的结果进行对比,其对比结果也相当吻合.说明利用该程序求解薄板自由振动时的固有频率是可行的,且可快速求解出矩形薄板在四边固支时自由振动的频率,具有一定的工程应用价值.【期刊名称】《铁道建筑》【年(卷),期】2010(000)008【总页数】3页(P142-144)【关键词】矩形薄板;Matlab软件;自由振动;固有频率【作者】张强【作者单位】西南交通大学土木工程学院,成都,610031【正文语种】中文【中图分类】TU33+9薄板结构在土木工程中应用较为广泛，土木工程结构由于振动特性造成的破裂、损害、倾覆和坍塌等破坏事故已为实践所证实，正确分析与评估土木工程结构在动力作用下的安全性，在土木工程结构的设计中起着十分重要的作用。

因此，为了满足工程上各种不同的需要，对薄板振动性能的研究是非常必要的。

在一定的横向荷载作用下处于平衡位置的薄板，受到干扰力的作用而偏离这一位置，当干扰力被除去后，薄板在该平衡位置附近作微幅振动。

目前对矩形薄板横向振动的固有频率计算方法主要有解析法与数值法，其中数值法有差分法，Reyleigh-Ritz法等。

Matlab软件是一款应用较为广泛的数学软件，在土木工程、机械工程等领域都有所应用。

Matlab的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用Matlab来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，Matlab软件具有基本的数值运算、矩阵运算、绘制图形等多种功能，特别是Matlab具有十分强大的矩阵运算功能。

由于在弹性矩形薄板的有限元分析中会涉及到大量的大型矩阵，所以利用Matlab软件进行弹性矩形薄板的有限元分析较为方便快捷。

固有频率和振型分析的平板算例Workbench版人们在制造动力机械、建造桥梁等工程实践中遇到大量灾害性振动问题，由此产生的噪声、疲劳等问题，吸引众多力学家和工程师致力于工程振动问题的研究，发展了近似分析方法、实验方法和有限元等方法。

自从20世纪20年代起，振动逐渐成为机械工程师、结构工程师必须了解的知识，也是高等工程教育的重要内容之一。

上一篇采用ANSYS经典版，对四边铰支平板的固有频率和振型进行了分析。

本篇将采用ANSYS Workbench版，对同一问题进行算例验证。

一、问题描述：某一个矩形薄钢板，板的长度a=1 m，宽度b=1 m，厚度h=4 mm。

材料密度ρ=7850 kg/m3，弹性模量E=200 GPa，泊松比μ=0.3。

假设矩形薄板的四边铰支，计算该薄钢板的固有频率和振型。

二、问题分析：弹性薄板是指厚度比平面尺寸小很多的弹性体，它可提供抗弯刚度。

在板中，与两表面等距离的平面成为中面。

对板弯曲振动的分析基于下述Kirchhoff假设：（1）微振动时，板的挠度远小于厚度，从而中面挠曲线为中性面，中面内无应变。

（2）垂直于平面的法线在板弯曲后仍为直线，且垂直于挠曲线后的中面；该假设等价于忽略横向剪切变形。

（3）板弯曲变形时，板的厚度变化可忽略不计。

（4）板的惯性主要由平动的质量提供，忽略由于弯曲而产生的转动惯量。

根据以上Kirchhoff假设，薄板固有频率的解析解为解析解参考文献：《机械振动基础》，胡海岩，pp118-121。

三、计算结果：四、操作步骤：视频录制：蒋豹。

来源：好学ANSYS(ID:ansys-good）转载分享。

四边弹性约束FGM矩形板面内自由振动的DQM求解蒲育;赵海英;滕兆春【摘要】假设矩形板为正交各向异性，材料的物性沿矩形板的宽度方向按幂律连续分布，基于二维线弹性理论，建立了四边弹性约束功能梯度材料（Functionally Graded Material，FGM）矩形板面内自由振动的控制偏微分方程。

控制方程为复杂耦合的变系数偏微分方程，采用微分求积法（Differential Quadrature Method，DQM）数值研究了四边弹性约束FGM矩形板面内自由振动的无量纲频率特性。

通过设置弹性刚度系数为0或∞，梯度指数为0，问题退化为各种典型边界下矩形板的面内自由振动，与已有的各向同性矩形板自振频率结果进行比较，结果表明分析求解方法行之有效。

最后考虑了FGM矩形板边界条件、长宽比、梯度指数及刚度系数对自振频率的影响。

%The material of rectangular plates was assumed to be orthotropic,and material properties change continuously along the width of a rectangular plate according to power law distributions.Based on the two-dimension theory of linear elasticity,the governing partial differential equations for in-plane free vibration of FGM rectangular plates with 4 elastically restrained edges were derived.The partial differential equations were complicated and coupled with variable ing the differential quadrature method,dimensionless frequency characteristics of in-plane free vibration of FGM rectangular plates with 4 elastically restrained edges were investigated.All the typical boundaries for in-plane vibration of isotropic rectangular plates were obtained by setting stiffnesses of restraining springs to be either zero or infinite and taking material gradient index as zero.Then,the results withDQM were compared with those published in literature for isotropic rectangular plates,it was shwon that the proposed DQM iseffective.Finally,The influences of boundary conditions,geometrical parameters,material gradient index and stiffness coefficients on the natural frequencies of FGM rectangular plates were analyzed.【期刊名称】《振动与冲击》【年(卷),期】2016(035)017【总页数】8页(P58-65)【关键词】FGM矩形板;面内自由振动;弹性约束边界;无量纲频率;微分求积法【作者】蒲育;赵海英;滕兆春【作者单位】兰州工业学院土木工程学院，兰州 730050;兰州工业学院土木工程学院，兰州 730050;兰州理工大学理学院工程力学系，兰州 730050【正文语种】中文【中图分类】O343;O327平板的振动表现为三种波形：弯曲波，纵波(p波)和剪切波(s波)。

四边固支矩形薄板固有振动的理论计算和有限元分析四边固支矩形薄板是一种典型的结构，其固有振动特性的计算对于结构的稳定性以及对外载荷的响应有着重要的影响。

本文将从理论计算和有限元分析两个方面来探讨四边固支矩形薄板的固有振动特性。

一、理论计算在理论计算中，四边固支矩形薄板的固有振动频率可以通过以下公式进行计算：f_n = (C_n^2 + D_n^2)^0.5 / (2πt)^0.5 * (EH^3/12ρ(1-μ^2))，其中，f_n为第n阶固有频率；C_n和D_n分别为第n阶水平和竖直模态振型的振幅比；t为薄板厚度；E为材料的弹性模量；H为矩形薄板的一侧长度；ρ为材料的密度；μ为材料的泊松比。

根据上述公式，我们可以对四边固支矩形薄板进行理论计算，得出其固有振动频率，并根据振动模型分析结构的稳定性以及响应能力。

二、有限元分析在有限元分析中，我们可以通过建立合适的有限元模型，利用求解振型特征值和振型模态来得出四边固支矩形薄板的固有振动特性。

有限元分析的主要步骤包括：1.建立有限元模型：根据实际结构情况，选择合适的有限元支撑和单元类型，对结构进行离散化网格化处理，建立结构有限元模型。

2.确定边界条件：对于固支矩形薄板，边界条件为四边界固定支撑。

3.求解特征值和振型：对于固有振动频率，我们可以通过求解振型特征值和振型模态来得出。

4.分析特征值和振型：得出固有振动频率，我们可以进一步分析与理论计算结果的一致性，同时还可以分析振型特征值与振型模态，进一步了解结构的稳定性和响应能力。

通过有限元分析，我们可以更加精确地了解四边固支矩形薄板的固有振动特性，为结构设计和应用提供更加实际的参考依据。

总之，四边固支矩形薄板的固有振动特性对于结构稳定性和响应能力有着重要的影响。

通过理论计算和有限元分析两个方面的探讨，我们可以更好地理解并应用这一结构特性。

为了更加深入地了解四边固支矩形薄板的固有振动特性，我们可以从以下几个方面进行数据的收集和分析：1. 材料弹性模量与密度：材料的弹性模量和密度直接影响到四边固支矩形薄板的固有振动频率。

自振频率计算公式自振频率计算公式是用来计算系统的自振频率的一种公式，也被称为振荡频率计算公式。

它使用物理学中的振荡理论来计算振荡频率，可以用来对系统进行控制。

这个公式也可以用于分析振荡系统中的不同参数，以及预测振荡系统的振荡特性。

自振频率计算公式的基本原理是：量化的质量是一个系统的特性，它是由系统中的外力和内部阻尼决定的。

量化的质量是一个系统的能量在一个振荡周期内的改变，可以用来计算系统的自振频率。

自振频率计算公式可以设计为：自振频率 =k/m其中，k是系统的外力系数，m是系统中量化的质量。

借助该公式，可以根据外力系数和质量，计算出系统的自振频率。

例如，如果系统中的外力系数是k=20，量化的质量是m=10，则根据这个公式，系统的自振频率将会是：自振频率 =20/10 = 2 Hz因此，通过自振频率计算公式，可以结合系统的外力系数和量化的质量，得出系统的自振频率。

除此之外，自振频率计算公式还可以用于分析振荡系统中的振荡特性。

自振频率的变化可以反映出系统中的不同参数，这些参数又能够影响系统的振荡特性。

例如，如果外力系数和质量都发生变化，则振荡频率也会随之发生变化，从而影响振荡系统的性能。

自振频率计算公式主要用于控制和分析振荡系统，可以有效地预测振荡系统的振荡特性，让工程师能够更加准确地控制系统的性能。

因此，它在振荡系统的控制和分析中扮演着重要的角色，有助于保持系统性能的稳定和完美。

自振频率计算公式在物理学和工程科学领域都有广泛的应用。

目前，它被广泛用于各种振荡系统的控制和分析，如电子元件、电路、机械结构和精密仪器等。

此外，它还可以用来分析复杂的控制系统中的不同参数，从而更好地了解和控制振荡系统的振荡特性。

因此，自振频率计算公式在物理学和工程科学领域具有重要的意义，可以有效地控制和分析各类振荡系统，从而更好地维护系统的稳定性和完美性。

自振频率计算公式例题解析在物理学中，自振频率是指一个物体在没有外力作用下，以自然频率进行振动的频率。

这个概念在工程学和物理学中都有着重要的应用，因此了解如何计算自振频率是非常重要的。

本文将通过例题的解析，帮助读者更好地理解自振频率的计算公式和应用。

自振频率的计算公式如下：\[f = \frac{1}{2\pi}\sqrt{\frac{k}{m}}\]其中，f代表自振频率，k代表弹簧的弹性系数，m代表物体的质量。

这个公式告诉我们，自振频率与弹簧的弹性系数和物体的质量有关，而与振幅和阻尼无关。

现在，让我们通过一个例题来解析自振频率的计算过程。

例题，一个质量为2kg的物体悬挂在一个弹簧上，当物体受到外力拉伸弹簧10cm后，弹簧的弹性系数为200N/m。

求这个系统的自振频率。

解析：首先，我们可以利用胡克定律来计算弹簧的弹性系数。

胡克定律表示弹簧的弹性系数与弹簧的弹性形变成正比，即F=kx，其中F为弹簧的弹力，k为弹性系数，x为弹性形变。

根据题目给出的信息，我们可以得到：\[k = \frac{F}{x} = \frac{200N}{0.1m} = 2000N/m\]接下来，我们可以利用自振频率的计算公式来计算系统的自振频率。

根据公式：\[f = \frac{1}{2\pi}\sqrt{\frac{k}{m}}\]代入已知的数值，我们可以得到：\[f = \frac{1}{2\pi}\sqrt{\frac{2000N/m}{2kg}}\]\[f = \frac{1}{2\pi}\sqrt{1000}\]\[f = \frac{1}{2\pi}\times 31.62\]\[f ≈ 5Hz\]因此，这个系统的自振频率约为5Hz。

通过这个例题的解析，我们可以看到自振频率的计算过程并不复杂。

只需要利用弹簧的弹性系数和物体的质量，就可以轻松地计算出系统的自振频率。

这个公式在工程学和物理学中都有着广泛的应用，可以帮助工程师和科学家们更好地设计和研究振动系统。