版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
结构力学基础概念:虚功原理:虚功原理的历史与发展1结构力学与虚功原理的概述在结构力学领域,虚功原理是一个核心概念,它提供了一种分析结构在各种载荷作用下行为的有效方法。虚功原理基于能量守恒的原理,允许工程师通过计算结构在虚拟位移下的能量变化来预测其真实行为。这一原理在解决复杂结构问题时尤其有用,因为它可以简化分析过程,避免直接求解复杂的力平衡方程。1.1虚功原理的数学表述虚功原理可以表述为:在任何平衡状态下,外力对虚位移做的虚功等于内力对同一虚位移做的虚功。数学上,这可以表示为:δ其中,δW是外力对虚位移做的虚功,δ1.1.1虚功原理的应用虚功原理在工程中有着广泛的应用,特别是在结构分析和优化设计中。例如,它可以用于:结构稳定性分析:通过计算结构在不同载荷下的虚功,可以评估结构的稳定性。结构优化:在设计阶段,虚功原理可以帮助工程师找到在给定约束下结构性能最优的设计方案。振动分析:虚功原理也是分析结构振动特性的基础,通过它,可以求解结构的固有频率和振型。1.2虚功原理在工程中的应用实例1.2.1实例1:梁的弯曲分析考虑一个简支梁,两端固定,受到垂直于梁轴线的集中力F的作用。我们可以通过虚功原理来分析梁的弯曲行为。假设梁在某一虚拟位移δy下弯曲,那么外力F对虚位移δδ同时,梁的内力(主要是弯矩M)对同一虚位移δyδ其中,L是梁的长度,x是梁的坐标轴。通过设置δW1.2.2实例2:结构优化设计在结构优化设计中,虚功原理可以用于最小化结构的重量或成本,同时满足特定的性能要求。例如,设计一个桥梁,我们需要找到在给定载荷和安全系数下的最优截面尺寸。假设桥梁的截面尺寸为A,载荷为P,安全系数为S。我们可以通过计算不同截面尺寸下的虚功,找到满足安全要求的最小重量设计。δδ通过调整A,使δW=δ1.2.3实例3:振动分析虚功原理在振动分析中也扮演着重要角色。例如,分析一个单自由度系统的振动特性,我们可以使用虚功原理来求解系统的固有频率。假设系统受到一个虚拟的位移δx,那么系统的动能变化δT和势能变化δδ其中,m是系统的质量,k是系统的刚度,δx通过设置δT1.3结论虚功原理是结构力学中一个强大的工具,它不仅简化了结构分析的过程,还为结构优化和振动分析提供了理论基础。通过理解和应用虚功原理,工程师可以更有效地解决工程中的结构问题。2虚功原理的历史2.1早期的力学理论与虚功概念的萌芽虚功原理的起源可以追溯到17世纪,当时力学理论正处于萌芽阶段。这一时期,许多伟大的科学家和数学家,如伽利略、牛顿和莱布尼茨,都在探索力与运动之间的关系。虚功概念最初是在研究静力学问题时提出的,它提供了一种分析系统在平衡状态下的方法。虚功原理的基本思想是,当一个系统处于平衡状态时,所有作用在系统上的力对任意虚位移所做的功之和为零。2.1.1伽利略的贡献伽利略在17世纪初对落体运动的研究中,虽然没有直接提出虚功原理,但他对力和运动关系的探索为后来的虚功原理奠定了基础。伽利略通过实验观察和数学分析,提出了自由落体运动的规律,即物体下落的距离与时间的平方成正比,这为理解力的作用提供了新的视角。2.1.2牛顿的三大定律牛顿的三大运动定律进一步深化了对力和运动的理解。牛顿第一定律,即惯性定律,指出在没有外力作用时,物体会保持静止或匀速直线运动的状态。牛顿第二定律,F=ma,描述了力与加速度之间的关系,即物体的加速度与作用在它上面的净力成正比,与它的质量成反比。牛顿第三定律,作用与反作用定律,说明了相互作用的两个物体之间力的大小相等、方向相反。这些定律为虚功原理的数学表述提供了理论框架。2.2达朗贝尔与虚功原理的数学表述虚功原理的数学表述主要归功于法国数学家和物理学家达朗贝尔(JeanleRondd’Alembert)。在18世纪中叶,达朗贝尔在研究力学问题时,引入了虚位移的概念,并将其与力的功联系起来,从而形成了虚功原理的初步数学形式。2.2.1虚位移与虚功虚位移是指在约束条件下,系统中各质点可能发生的无限小位移,这些位移并不实际发生,而是假设的。虚功则是指作用在系统上的力对虚位移所做的功。达朗贝尔认识到,如果一个系统处于平衡状态,那么所有作用在系统上的力对任意虚位移所做的虚功之和必须为零。这一原理可以表示为:i其中,Fi是作用在第i个质点上的力,δ2.2.2达朗贝尔原理达朗贝尔进一步发展了这一思想,提出了达朗贝尔原理,即在任何瞬时,一个系统的动力学状态可以视为一系列平衡状态的叠加,其中每个平衡状态对应于一个虚位移。这一原理为分析复杂动力学问题提供了一种有效的方法。2.3虚功原理在19世纪的发展与完善进入19世纪,虚功原理得到了进一步的发展和完善。这一时期,许多科学家和工程师开始将虚功原理应用于解决实际工程问题,如桥梁设计、机械结构分析等。虚功原理的应用范围也从静力学扩展到了动力学和振动学领域。2.3.1拉格朗日的贡献意大利数学家拉格朗日(Joseph-LouisLagrange)在18世纪末提出了拉格朗日力学,这是一种基于虚功原理的分析力学方法。拉格朗日引入了广义坐标和广义力的概念,使得虚功原理可以应用于更广泛的系统,包括那些具有多个自由度的复杂系统。拉格朗日方程是虚功原理在动力学问题中的具体应用,它描述了系统在任意虚位移下的动力学平衡条件。2.3.2虚功原理在工程中的应用19世纪,随着工业革命的推进,虚功原理在工程设计和分析中的应用日益广泛。工程师们开始使用虚功原理来计算结构的稳定性、强度和刚度,以及预测结构在不同载荷下的响应。虚功原理的应用不仅限于静态分析,还扩展到了动态分析,如振动分析和冲击分析。2.3.3虚功原理的现代形式到了20世纪,虚功原理的数学表述得到了进一步的简化和标准化。现代虚功原理通常以能量的形式表述,即在平衡状态下,系统内部能量的变化等于外部力做功的变化。这一表述不仅适用于线性系统,也适用于非线性系统,为现代结构力学和工程分析提供了强大的工具。2.4结论虚功原理从早期的力学理论萌芽,经过达朗贝尔和拉格朗日等科学家的数学表述和理论发展,到19世纪在工程中的广泛应用,再到20世纪的现代形式,其历史和发展反映了力学理论与工程实践的紧密联系。虚功原理不仅是一种理论工具,也是解决实际工程问题的有效方法,它在结构力学、机械工程、土木工程等领域有着广泛的应用。请注意,上述内容中没有包含代码示例,因为虚功原理的讨论主要基于理论和数学公式,而不是具体的编程实现。然而,如果在工程分析软件中实现虚功原理的计算,通常会涉及到数值方法和线性代数,这些领域的代码示例可以使用Python、MATLAB等编程语言来编写,但具体实现超出了本教程的范围。3虚功原理的基本概念3.1虚位移的定义与性质在结构力学中,虚位移(VirtualDisplacement)是一个重要的概念,它指的是在约束条件下,系统中各质点可能发生的、与时间无关的、无限小的位移。虚位移满足系统的所有约束条件,但并不一定是实际发生的位移。虚位移的性质包括:无限小性:虚位移是无限小的,这意味着它在数学上可以被视为一个微分量。与时间无关:虚位移不随时间变化,它是一个静态的概念。满足约束条件:虚位移必须满足系统的所有约束条件,例如,如果一个质点被固定,那么它的虚位移为零。3.2虚功的计算方法虚功(VirtualWork)是力在虚位移上所做的功。如果一个力F作用在系统上,而系统发生虚位移u,那么力F在虚位移u上所做的虚功W可以表示为:W在连续介质中,虚功可以表示为应力σ和应变ε的内积:W其中,V是系统的体积,:表示双点积,即对应分量的乘积之和。3.2.1示例假设有一个简单的梁,受到一个垂直向下的力F,梁的虚位移为ux,其中xW如果梁的长度为L,虚位移ux在x=0W3.3虚功原理的数学表达式虚功原理(PrincipleofVirtualWork)是结构力学中的一个基本原理,它表明在平衡状态下,所有外力和内力在任意虚位移上所做的虚功之和为零。对于一个静力平衡的系统,虚功原理可以表示为:δ其中,δW外力是外力在虚位移上所做的虚功,3.3.1示例考虑一个受力的梁,其外力为Fx,内力(例如,弯矩)为Mx,虚位移为δδ在平衡状态下,虚功原理要求:δ即:0这个方程可以用来求解梁的弯矩分布Mx,或者验证一个给定的力分布Fx和位移分布以上内容详细介绍了虚功原理的基本概念,包括虚位移的定义与性质、虚功的计算方法以及虚功原理的数学表达式。通过这些概念和公式,我们可以更好地理解和应用虚功原理在结构力学中的作用。4虚功原理的应用4.1静态结构的平衡分析虚功原理在静态结构分析中是一个强大的工具,它基于能量守恒的概念,用于确定结构在给定载荷下的平衡状态。虚功原理的核心是,如果一个结构处于平衡状态,那么所有作用在结构上的外力对任意虚位移所做的虚功之和等于零。4.1.1原理描述考虑一个受外力作用的结构,如果该结构处于平衡状态,那么对于任意虚位移,外力所做的虚功与内力(由变形引起的)所做的虚功相等。数学上,这可以表示为:δ其中,δW外是外力对虚位移所做的虚功,4.1.2示例假设有一个简支梁,长度为L,受到均匀分布载荷q的作用。我们想要使用虚功原理来验证该梁在给定载荷下的平衡状态。计算外力的虚功:假设梁的任意虚位移为δyδ计算内力的虚功:梁的内力主要是弯矩Mx和剪力VM内力的虚功可以通过梁的挠度方程和虚位移来计算。梁的挠度方程yxδ其中,EI验证平衡状态:如果δW4.2动态结构的振动分析虚功原理在动态结构分析中同样重要,特别是在振动分析中。它可以帮助我们建立结构的动力学方程,从而分析结构的振动特性。4.2.1原理描述在动态分析中,虚功原理考虑了动能和势能的变化。如果一个结构在振动过程中处于平衡状态,那么在任意时刻,外力的虚功、内力的虚功、动能的虚变和势能的虚变之和等于零。δ其中,δT是动能的虚变,δ4.2.2示例考虑一个单自由度系统,质量为m,弹簧刚度为k,受到外力Ft计算外力的虚功:假设系统的虚位移为δxδ计算内力的虚功:内力主要是弹簧力,其虚功为:δ计算动能的虚变:系统的动能为12mv2δ计算势能的虚变:系统的势能为12δ建立动力学方程:将上述虚功和虚变代入虚功原理的方程中,可以得到系统的动力学方程:m4.3能量法求解结构问题能量法是基于虚功原理的一种求解结构问题的方法,它通过最小化结构的总势能或总能量来确定结构的平衡状态或最优状态。4.3.1原理描述能量法的核心是能量原理,即在给定的约束条件下,结构的平衡状态对应于总势能或总能量的最小值。对于静态结构,总势能可以表示为:Π其中,V是结构的内能(势能),W外4.3.2示例假设我们有一个悬臂梁,长度为L,受到端部集中载荷P的作用。我们想要使用能量法来确定梁的挠度。建立总势能表达式:悬臂梁的内能(势能)可以表示为:V外力所做的功为:W因此,总势能为:Π求解挠度:为了求解挠度yx,我们需要找到使总势能Π最小的函数yx。这可以通过变分法来实现,即求解δ解这个方程可以得到梁的挠度方程。以上示例展示了虚功原理在静态结构平衡分析、动态结构振动分析以及能量法求解结构问题中的应用。通过这些方法,我们可以更深入地理解结构的行为,并有效地解决工程中的实际问题。5虚功原理的现代发展5.1计算机辅助工程中的虚功原理应用在计算机辅助工程(CAE)领域,虚功原理被广泛应用于结构分析和优化设计中。这一原理允许工程师在不考虑实际位移路径的情况下,通过计算虚拟位移下的功来分析结构的平衡状态,从而简化了复杂结构的分析过程。在CAE软件中,虚功原理通常与有限元方法(FEM)结合使用,以解决静态和动态问题。5.1.1示例:使用虚功原理进行结构优化假设我们有一个简单的梁结构,需要通过虚功原理来优化其设计,以减少材料使用量同时保持结构的稳定性。在CAE软件中,我们可以通过定义虚拟位移和计算虚拟功来实现这一目标。数据样例结构参数:梁的长度为10米,宽度为0.2米,高度为0.1米。材料属性:弹性模量E=200GPa,泊松比ν=0.3。载荷条件:在梁的中点施加垂直向下的力F=10kN。计算过程定义虚拟位移:在梁的中点定义一个垂直向上的虚拟位移δy。计算虚拟功:通过计算虚拟位移下外力所做的功和材料内部应力所做的功,来确定结构的平衡状态。优化设计:调整梁的尺寸,直到虚拟功为零,表明结构在给定载荷下达到平衡状态。5.1.2代码示例以下是一个使用Python和NumPy库来计算梁结构在虚拟位移下的虚拟功的简化示例:importnumpyasnp
#结构参数
length=10.0#梁的长度
width=0.2#梁的宽度
height=0.1#梁的高度
#材料属性
E=200e9#弹性模量
nu=0.3#泊松比
#载荷条件
F=10e3#施加的力
delta_y=0.01#虚拟位移
#计算截面惯性矩
I=(width*height**3)/12
#计算弯矩
M=F*length/4
#计算虚拟功
virtual_work=(M*delta_y)/(E*I)
print("虚拟功:",virtual_work)代码解释计算截面惯性矩:使用梁的宽度和高度计算截面惯性矩I。计算弯矩:基于梁的长度和施加的力F,计算梁中点的弯矩M。计算虚拟功:通过弯矩M、虚拟位移δy、弹性模量E和截面惯性矩I,计算虚拟功。5.2现代结构分析软件中的虚功算法现代结构分析软件,如ANSYS、ABAQUS和NASTRAN,都集成了虚功原理的算法,用于解决复杂的结构力学问题。这些软件通过有限元方法将结构离散化,然后在每个节点上应用虚功原理,以计算结构的响应。5.2.1示例:ABAQUS中的虚功分析在ABAQUS中,虚功原理可以用于进行非线性分析,包括大变形、接触和材料非线性。通过定义虚拟位移和计算虚拟功,软件能够自动调整求解策略,以确保分析的准确性。数据样例模型:一个包含接触面的复杂结构模型。材料:非线性材料属性,如塑性材料。载荷:非线性载荷,如温度变化引起的热应力。分析过程定义模型:在ABAQUS中创建结构模型,包括几何、材料和载荷条件。设置分析类型:选择非线性分析,并启用虚功原理。求解:运行分析,ABAQUS将自动应用虚功原理来解决非线性问题。5.3虚功原理在非线性力学中的扩展虚功原理不仅适用于线性问题,其在非线性力学中的应用也得到了显著的扩展。在处理大变形、接触和材料非线性等问题时,虚功原理提供了一种有效的分析方法,能够准确地预测结构的行为。5.3.1示例:大变形下的虚功分析考虑一个在大变形条件下工作的结构,如橡胶密封件。在这种情况下,虚功原理能够考虑材料的非线性响应,从而提供更准确的分析结果。数据样例材料:橡胶材料,具有非线性弹性行为。载荷:施加在结构上的力,导致显著的变形。分析过程定义材料模型:在分析软件中输入橡胶材料的非线性弹性模型。应用虚功原理:通过计算虚拟位移下的功,分析结构在大变形条件下的平衡状态。结果评估:评估结构的变形和应力分布,确保其在工作条件下能够正常运行。通过上述示例,我们可以看到虚功原理在现代工程分析中的重要性和实用性,特别是在处理复杂和非线性问题时。6案例研究与实践6.1桥梁结构的虚功分析在结构力学中,虚功原理被广泛应用于桥梁结构的分析中,尤其在评估桥梁的稳定性和动态响应时。虚功原理提供了一种有效的方法来计算结构在虚拟位移下的能量变化,从而帮助工程师理解结构在不同载荷条件下的行为。6.1.1桥梁结构分析中的虚功原理应用假设我们有一座简支梁桥,需要分析其在特定载荷下的响应。我们可以通过虚功原理来计算梁在载荷作用下的变形和内力。虚功原理的基本思想是,结构在虚拟位移下的外力虚功等于其在相同虚拟位移下的内力虚功。示例:简支梁桥的虚功分析考虑一个简支梁,长度为L,在中点受到集中力P的作用。我们可以通过虚功原理来计算梁的中点挠度δ。定义虚拟位移:假设梁在中点的虚拟位移为δv计算外力虚功:外力虚功Wext为外力P与虚拟位移δ计算内力虚功:内力虚功Wint需要通过梁的弯矩分布和对应的虚拟转角来计算。对于简支梁,弯矩分布可以通过梁的截面刚度和载荷分布来确定。假设梁的截面刚度为EI,则梁在中点的弯矩为M=PL4。虚拟转角应用虚功原理:根据虚功原理,Wext=W6.1.2Python代码示例importsympyassp
#定义变量
P,L,EI,x,delta_v=sp.symbols('PLEIxdelta_v')
#定义弯矩函数
M=P*L/4
#定义虚拟转角函数
theta_v=delta_v/(L/2)
#计算内力虚功
W_int=egrate(M*theta_v,(x,0,L))
#根据虚功原理求解中点挠度
delta=sp.solve(P*delta_v-W_int,delta_v)
print("中点挠度为:",delta[0])此代码示例使用SymPy库来计算简支梁在中点集中力作用下的中点挠度。通过定义弯矩和虚拟转角的函数,然后计算内力虚功,最后应用虚功原理求解挠度。6.2高层建筑的动态响应分析虚功原理在高层建筑的动态响应分析中也扮演着重要角色。通过虚功原理,工程师可以评估建筑在风载、地震等动态载荷下的响应,确保建筑的安全性和稳定性。6.2.1高层建筑动态响应分析中的虚功原理在动态分析中,虚功原理可以用于计算结构在虚拟位移下的动能和势能变化。这有助于工程师理解结构在动态载荷作用下的振动模式和频率。示例:高层建筑的风载响应分析考虑一座高层建筑,需要分析其在风载作用下的动态响应。我们可以通过虚功原理来计算建筑在风载下的振动。定义虚拟位移:假设建筑在风载作用下的虚拟位移为uv计算外力虚功:外力虚功Wext为风载与虚拟位移uv的乘积,即We计算内力虚功:内力虚功Wi应用虚功原理:根据虚功原理,Wext=Wint,即6.2.2Python代码示例importnumpyasnp
fromegrateimportquad
#定义风载分布函数
defwind_load(x):
return0.5*1.225*20*(10-x)**2#假设风速随高度线性减少
#定义虚拟位移函数
defvirtual_displacement(x):
returnnp.sin(np.pi*x/100)#假设虚拟位移为正弦函数
#计算外力虚功
W_ext,_=quad(lambdax:wind_load(x)*virtual_displacement(x),0,100)
#假设内力分布已知,此处仅示例外力虚功计算
print("外力虚功为:",W_ext)此代码示例使用SciPy库的quad函数来计算高层建筑在风载作用下的外力虚功。通过定义风载分布和虚拟位移的函数,然后计算外力虚功。6.3复合材料结构的虚功应用复合材料因其独特的性能在现代工程中得到广泛应用。虚功原理在复合材料结构的分析中,帮助工程师理解材料的各向异性对结构性能的影响。6.3.1复合材料结构分析中的虚功原理在复合材料结构中,虚功原理可以用于计算结构在虚拟位移下的能量变化,考虑到材料的各向异性。这有助于工程师优化复合材料结构的设计,提高其效率和性能。示例:复合材料梁的虚功分析考虑一根复合材料梁,需要分析其在特定载荷下的响应。我们可以通过虚功原理来计算梁的变形和内力,考虑到材料的各向异性。定义虚拟位移:假设梁在载荷作用下的虚拟位移为uv计算外力虚功:外力虚功Wext为外力与虚拟位移uv的乘积,即We计算内力虚功:内力虚功Wi应用虚功原理:根据虚功原理,Wext=Wint,即6.3.2Python代码示例importnumpyasnp
fromegrateimportquad
#定义载荷分布函数
defload(x):
return1000*(1-x/10)#假设载荷沿长度线性减少
#定义虚拟位移函数
defvirtual_displacement(x):
returnnp.sin(np.pi*x/10)#假设虚拟位移为正弦函数
#计算外力虚功
W_ext,_=quad(lambdax:load(x)*virtual_displacement(x),0,10)
#假设内力分布已知,此处仅示例外力虚功计算
print("外力虚功为:",W_ext)此代码示例使用SciPy库的quad函数来计算复合材料梁在外力作用下的外力虚功。通过定义载荷分布和虚拟位移的函数,然后计算外力虚功,考虑到复合材料的特性。通过这些案例研究,我们可以看到虚功原理在不同结构分析中的应用,以及如何通过Python编程来实现这些分析。虚功原理不仅提供了一种理论框架,还为工程师提供了一种实用的工具,用于理解和优化结构设计。#结论与展望
##虚功原理在结构力学中的重要性总结
在结构力学领域,虚功原理作为分析结构平衡状态和稳定性的一种强大工具,其重要性不言而喻。它不仅提供了结构在各种载荷作用下行为的深刻理解,而且在工程设计和分析中扮演着关键角色。虚功原理的核心在于,它允许我们通过考虑虚拟位移和虚拟力来评估结构的真实状态,而无需直接求解复杂的力平衡方程。
###虚功原理的应用
-**结构分析**:虚功原理可以用于确定结构的内力、位移和应力分布,特别是在处理非线性问题时,其灵活性和通用性尤为突出。
-**能量方法**:在能量方法中,虚功原理是求解结构最小势能状态的基础,从而帮助工程师设计出既经济又安全的结构。
-**有限元分析**:现代工程分析软件广泛采用虚功原理作为其背后的数学模型,通过离散化结构,使用有限元方法来求解复杂的工程问题。
###虚功原理的数学表达
虚功原理的基本数学表达式为:
$$
\deltaW=\delta\int\mathbf{f}\c
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2025年江门货运资格证500道题库
- 单车位租赁合同范例
- 婚礼跟妆合同范例
- 2025年新疆货运车从业考试题
- 显微镜购买合同范例
- 2025年宜春年货运从业资格证考试从业从业资格资格题库及答案
- 天府新区航空旅游职业学院《环境设计专题》2023-2024学年第一学期期末试卷
- 《12 图文并茂-精确设置图片尺寸》教学实录-2023-2024学年清华版(2012)信息技术三年级下册
- 2025年山东货物运输从业资格考试答题软件
- 2025年凉山州驾驶资格证模拟考试
- 低空经济产业园项目可行性研究报告
- 中国神话故事绘本仓颉造字
- MOOC 心理健康与创新能力-电子科技大学 中国大学慕课答案
- 中华传统造型的艺术之美-中国美术史专题精讲智慧树知到期末考试答案章节答案2024年山东工艺美术学院
- 黄蒿界矿井及选煤厂建设项目环境影响报告书
- 2023-2024学年高一下学期家长会 课件
- 感动中国人物张桂梅心得体会(30篇)
- 知识点总结(知识清单)-2023-2024学年人教PEP版英语六年级上册
- 社会医学课件第2章医学模式-2024鲜版
- 德勤测评能力测试题及答案
- 《囚歌》教学课件
评论
0/150
提交评论