版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、机构的运动分析是已知机构原动件的运动规律,对机构某点或某构件进行位移(角位移)、速度(角速度)和加速度(角加速度)分析。这些分析无论是对了解现有机械的运动性能及设计新的机械都是十分必要的。运动分析是完善机构综合的重要步骤之一,通过运动分析可以计算构件的惯性力、了解机械的受力情况和研究机械的动力性能。因此,机构的运动分析是对机构进行受力分析的基础和必要的前提。机构的力分析包括两部分,一部分是考虑摩擦的受力分析;另一部分是动态静力分析。前者要考虑机构中各构件的相对运动关系,后者要计算出机构在各个位置的速度和加速度计算惯性力。这两部分都要求在机构运动过程中各运动副中的总反力和平衡力(或平衡力矩),为
2、进一步计算各构件的强度、刚度及结构尺寸提供依据。无论是运动分析还是受力分析其具体解法都有图解法和解析法两种,图解法的特点是直观、易懂,但不精确;解析法的特点是将机构放在直角坐标系下,将已知和未知的运动量之间的关系用数学式子表达出来,然后求解。随着计算机的普及和发展,解析法已逐渐推广,并用于生产实际中。在本章学习过程中首先讲解图解法,以求对问题的理解,再着重讲解解析法。3.1 速度瞬心用图解法分析机构的速度,有速度瞬心法和矢量方程图解法等。对有些机构应用速度瞬心法求机构中某点的速度或某构件的角速度是十分简便的。3.1.1 速度瞬心.PijAvA1A2 B(A1A2)(B1B2)vB1B2 .Pi
3、j互作平行平面运动的两构件,在任一瞬时其相对速度为零,绝对速度相等的瞬时重合点称为该两构件的速度瞬心,简称瞬心。若该点上的绝对速度为零,则该点的瞬心称为绝对瞬心;若该点上绝对速度非零,则该点的瞬心称为相对瞬心。一般用符号Pij(或Pji)表示构件i和构件j的瞬心。如图3.1所示,1、2构件互作平行平面运动,在该瞬时1、2构件上A、B各点的相对速度是绕P12这一瞬时重合点运动的,P12即为1、2构件的瞬心。若1、2构件都在运动,则此时的P12点是相对瞬心;若1、2构件中有一个构件固定,则此时的P12即为绝对瞬心。所以判断作平行平面运动的两构件某一瞬时重合点是否是绝对瞬心,主要看其中某构件是否与机
4、架即固定件组成瞬心。3.1.2 机构中瞬心的数目根据瞬心定义和表示方法,可见Pij亦是Pji,与构件i、j排列的次序无关。若机构中有N个构件(包括机架在内),每两个构件存在一个瞬心,则机构中总的瞬心数K的求解是一个组合问题。机构中总的瞬心数为: (3.1)3.1.3 机构中瞬心位置的确定1、直接成副两构件的瞬心位置的确定1)两构件由转动副联接 由转动副相联的两构件,其铰接中心点即为瞬心点。P12P12(b)(a)2121AA如图3.2 (a)、(b)所示,1、2构件在A点铰接,根据瞬心定义vA1A2=0,故此A点即为瞬心P12。(a)、(b)图中的P12分别为绝对瞬心和相对瞬心。2)两构件由移
5、动副联接 由移动副相联的两构件,其瞬心点在垂直于导路的无穷远处。图3.3 由转动副组成的瞬心P12 (b)(a)21B2vB1B2B1vB1B2P12如图3.3 (a)、(b)所示,1、2构件在B点重合,且组成移动副,其相对速度vB1B2方向均沿着导路,可以看作vB1B2是绕垂直于导路无穷远处的一点转动。因此,P12瞬心在垂直导路的无穷远处。3)两构件由高副联接 由高副相联的两构件,其瞬心在过接触点的公法线上。所示,1、2构件在C点接触,且组成高副。若1、2构件之间的运动是无滑动的纯滚动,在接触点C处相对速度为零,则接触点C即为瞬心点P12,如(a)所示;若1、2构件的运动在接触点C处是连滚带
6、滑运动,则瞬心在公法线n-n上,在n-n线上哪一点,应在具体机构上去找,如(b)所示。P121vC1C2=0P121C2(b)(a)CvC1C2vC1C2¹02nn2、三心定理三心定理是解决不直接成副的两构件瞬心位置的确定问题。即:三个互作平行平面运动的构件共有三个瞬心,且这三个瞬心必在一条直线上。P23vK2P1312BAKvK1w2w13现证明如下:如图3.5所示,设构件1、2、3彼此间互作平行平面运动,总的瞬心数为3,其中P13、P23分别处于转动副A、B处。P12的位置根据三心定理应在P13、P23两点所在的连线上,即AB线上。下面利用反证法证明,若P12不在AB线上不成立,
7、则定理正确。1)若设P12在K点,如图3.5所示。由于1、2构件分别绕A、B两点转动,在图中可见,若,由于方向不一致,则该两速度不等。所以P12必定不在K点。2)若,要使两速度相等,则只有方向一致才成立;若使方向一致,K点就必须落在AB的连线上。所以得证P12必在P12 P23两点的连线上。下面举例说明三心定理的应用。例3.1 图3.6为一平面四杆机构,确定机构图示位置的全部瞬心。P12P14P34P24P23P133421四杆机构的全部瞬心 P14、P34P13 P12、P23 P12、P14P24 P23、P34解:机构的全部瞬心P14、P12、P23、P34可由直接成副的两构件的瞬心求法
8、标出。P24根据三心定理观察2、1、4和2、3、4二组构件,可见P24在连线上,也在连线上,这两线的交点即是P24;同理可求P13,P13是与两线的交点。这样就将该机构的六个瞬心全部求出。机构中构件数目较少时用上述方法求解机构的全部瞬心较简单,若构件数目较多时,则不易求解。下面介绍一种称为“瞬心多边形”法来求解机构的全部瞬心。其内容为:多边形的各顶点代表机构中的各构件,用数字表示;多边形的各顶点之间的连线分别代表两构件组成的瞬心,已知瞬心用实线画出,未知或要求的瞬心用虚线画出;三个顶点连线所形成的三角形即表示三瞬心共线,两个三角形公共边即表示两瞬心线的交点。图3.7为一曲柄滑块机构,确定图示机
9、构的全部瞬心。解:机构的瞬心数目P12P14P34®P24P23P133421图3.7 曲柄滑块机构的瞬心(a)(b) P14 4 3 2 1 P12 P23 P13 P24 P34 、可直接标在图上,根据“瞬心多边形”(b)所示。可见是与两连线的交点,其中P34瞬心是由34构件组成移动副的瞬心,该瞬心在垂直于导路的无穷远处,所以其方位线可在垂直于导路的方位上平移。P24是与两连线交点,在的(a)图标出即可。速度瞬心在速度分析中的应用利用速度瞬心对一些简单的平面机构进行速度分析既直观又方便。例3.3 已知:、尺寸和原动件1的角速度w1。比例尺为。求:构件1、3的传动比i13及构件3角
10、速度w3。14P12P34P23P132P23w14P1431243图3.8 求P13瞬心解:已知1构件的运动,求3构件的运动,应将机构中P13瞬心求出。利用“瞬心多边形”画出P13的位置如图示。1P12nP23P132P23w14n3图3.9求P12瞬心已知:、尺寸和原动件1的角速度w1,比例尺为。求:从动件2的速度。解:已知1构件的运动,求构件2的速度,应将机构中的瞬心求出。根据三心定理如图所示。如上所述,速度瞬心只能用来求解机构某点、某构件的速度和角速度,若要求解机构中的加速度,则需用其它方法。 3.2 机构的运动分析3.2.1 简介矢量方程图解法对机构进行运动分析14AEBD2w14C
11、3 p¢bcemac”c²¢ pebcmv(b)(c)(a)图3.10铰链四杆机构矢量方程图解法所依据的基本原理是理论力学中所介绍的刚体平面运动中的基点法和点的复合运动法。运用这两基本原理时,对于不同的构件、不同的点的运动求解时可能要反复利用多次,而且列矢量方程时必须标明各点的字母和各构件的数字。下面举例说明矢量方程图解法的应用。1、同一构件上两点间的速度和加速度求法例3.5 在图3.10(a)所示的铰链四杆机构中,已知:机构的位置、各构件长度及曲柄1的角速度为常数。求:连杆2的角速度,角加速度a2及其上点C和E的速度和加速度;构件3的角速度及角加速度a3。解:首
12、先选取长度比例尺ml,画出机构位置图。1)速度求解:(3.2)式中含二个未知量,可通过画矢量封闭图求解。选定速度比例尺,取p点(绝对速度为零的点)。代表矢量,方向:垂直于AB,图长:。过b点作垂直于BC的线,代表的方向线,过p点作垂直于CD的线代表的方向线,上述两方向线的交点即为c点。代表,代表,如图3-10(b)所示。方向如图。同理,E点也可根据基点法列出如下方程:式中有三个未知量不可解,故也可通过2构件中的EC两点列方程:)式中也有三个未知量不可解。但将式(3.3)式(3.4)联立可画图求解。)在由(3.2)式所画的速度图上,和已在图中画出,过b点做垂直于EB的线,过c点做垂直于EC的线,
13、两线的交点即是e点,即表示E点速度大小。则23构件的角速度分别为:对照图3-10 (a)、(b)可看出,在速度多边形中代表各相对速度的向量、和分别垂直于机构图中的BC、EC和BE。因此,且两三角形顶角字母b、c、e和B、C、E的顺序相同,均为顺时针方向,将速度图中的称为结构图中的影像。由上可见当已知一构件上两点的运动时,要求该构件上其它任一点的运动便可利用影像关系求解,这一原理称为影像原理。可以证明在同一构件上已知两点的加速度,要求该构件上任一点的加速度时,也有同样的加速度三角形与结构三角形相似的情况,也可以用影像关系求解。速度和加速度影像原理:1在同一构件上,若已知该构件上两点的速度和加速度
14、,求该构件上其它任一点的运动时可用影像关系求解;2速度和加速度图形上的字母绕行顺序应与结构图中字母绕行顺序一致。2)加速度求解,利用基点法。式中含二个未知量,可通过画矢量多边形求解。选定加速度比例尺,取点(绝对加速度为零的点)。根据式(3.6)画加速度矢量多边形,如图3.10(c)所示,代表矢量,方向:由B指向A,过点画线段方向由C指向B,;过点作垂直于BC的方向线;过点画线段,方向由C指向D,过点作垂直于CD的方向线;与两条方向线的交点即为点,这时C点绝对加速度可用表示。将连线,得相对加速度。根据影像原理可求出2构件上E点的加速度(c)所示。2、两构件瞬时重合点之间的速度和加速度求法例3.6
15、 如图3.11所示的导杆机构,已知机构的位置各构件长度及曲柄1的等角速度。求:导杆3的角速度w3和角加速度a3。解:根据长度比例尺ml画出机构位置图。1) 速度求解14CAw14B32(b)(c)(a) pb2b3mvpb2rkmab3”b3图3.11 导杆机构分析B点,1与2构件是铰链联接,故;2与3构件是移动副联接,。根据点的复合运动,将B3点作为动点,构件2为动系,方程如下:式中有两个未知量,可画图求解。选速度比例尺,取p点,根据矢量方程(3.7),先画,方向;过b2点作BC的方向线,代表vB3B2方位;过p点作BC的方向线,代表vB3方位,与两方向线的交点即为点,向量即代表。2) 加速
16、度求解,根据点的复合运动式中有两个未知量,可画图求解。选加速度比例尺,取点,根据矢量方程(3.8),先画,方向;过b2点画代表,过k点作方位线,方向平行BC,代表的方位线;过点作代表,过b3”作线,代表方位线。、两线交点为点。哥氏加速度的求解:大小:,由于在纸面内,的方向垂直纸面,所以q=90o,即;方向:将顺着的方向转。以上简介了矢量方程图解法的过程,从求解过程看,此种方法只能求解机构的某一位置上的速度和加速度,若要求解某一运动循环中各个任意位置上的运动,此法则显得慢而繁杂,而且不精确。特别是在计算机普及的今天,解析法则显现出其强大的优势。3.2.2 解析法对机构进行运动分析x.BAw1j1
17、CDl3l4l2l1图3.12封闭向量多边形投影法j2j3C¢y解析法一般是将机构放在直角坐标系下,建立机构在任一位置的位置方程,然后将位置方程对时间求导,即可得机构的速度和加速度方程,然后将所推导的方程编入程序计算。解析方法有向量法复数法封闭向量多边形投影法和拆杆组法等。下面主要向大家介绍封闭向量多边形投影法和拆杆组法。1封闭向量多边形投影法在如图3.12所示的铰链四杆机构中,已知各杆长分别为、,原动件1的转角为及等角速度。求:连杆2和摇杆3的角位移与;角速度与及角加速度a2与a3。解:1)位置分析将铰链四杆机构ABCD看作一封闭向量多边形,建立如图所示直角坐标系。将各杆长度看作是
18、向量,各向量与X轴正向夹角为、,、分别表示各构件的向量,该机构的封闭向量方程式为:将(3.9)式向XY轴投影得:在式(3.10)中、是未知数,消去后得:其中:解(3.11)式得:式(3.12)根号前的符号可按所给机构的装配方案来选择。“-”号适用于图示ABCD机构位置;“+”号适用于图示机构位置。构件2的角位移可通过(3.10)式求得:注:、在程序中求解后要注意其角度所在象限的判断。2)速度分析将(3.10)式对时间求导得:在(3.14)式中只有和未知,为了求应消去得:同理,在(3.14)式中消去得:3)、加速度分析将(3.14)式对时间求导得:在(3.17)式中只有a2和a3未知,为了求a2
19、应消去a3得:同理,为了求a3,在(3.17)式中消去a2得: 将上述求得的j2、j3、w2、w3、a2和a3编程上机计算,以j1为循环变量,对于不同的j1就会得到一组数值。当j1=0o360o变化时,可算出上述六个量的对应值,即可求出机构的运动曲线来,便于分析和比较。当机构的构件多时,为编制通用程序,可利用拆杆组法进行求解。2拆杆组法一般可将平面机构看成是由级机构和若干个自由度为零的基本杆组所组成,将级机构和各种基本杆组的运动方程列出,并分别编写成独立的子程序,在对一个机构进行运动分析时,仅需调用相应的子程序即可。下面我们主要来分析级杆组的方程式的推导,即数学模型的建立。1)级机构OxAyA
20、x图3.13级机构aiBjiliwiAyrBrA如图3.13所示,已知:A点的坐标(xA,yA ),AB杆的杆长li,及li与X轴正向夹角ji,角速度wi,角加速度ai。求:构件B点的速度和加速度。a、 位置分析:在直角坐标系中,B点的位置矢量为:投影方程为: b、速度分析:将(3.21)式对时间求导,得B点速度方程:c、加速度分析:2)式对时间求导,得B点加速度方程:根据(3.21)、(3.22)、(3.23)式,若A为固定点,则、均为零:若A为动点,为求出B点的运动,必须先给出A点的运动参数。表3.1 几种常见的级杆组级杆组2.RRP级杆组级杆组2)、RRR级杆组xyBOCDljli图3.
21、14 RRR级杆组jijjrBrDrC如图3.14所示,RRR级杆组是由两个构件和三个转动副所组成的级杆组。建立如图所示坐标系。已知:BC、CD杆长分别为li、lj,B、D两点的坐标及运动参数。求:C点的位置及运动参数。a、位置分析:C点的矢量为将向量向X、Y轴上投影得:消去jj,求ji,上式经整理得: 其中:有两组解,当BCD三运动副顺时针排列时取“+”,反之取“-”。将代入(3.25)式得:b、速度分析:对(3.25)式求导,得:将上式移项整理得:在(3.31)式中只有、未知,联立得:b、 加速度分析:对(3.30)式求导,得:将上式移项整理得:在(3.36)式中,只有ai、aj未知,联立
22、得: 3)、RRP级杆组RxyBOCDljli图3.15 RRP级杆组jijjsrRrBrC如图3.15所示的RRP级杆组。建立如图所示的坐标系。其中两构件的杆长分别为li、lj,构件li的角位置为,lj杆垂直于滑块导路,滑块导路与X轴正向夹角为jj,滑块上D点相对于参考点R的位移量用S表示。已知:B点和参考点R的位置、jj及运动参数。求:C点的位置及运动参数。 a、位置分析: 运动副C点的矢量方程为:将上式在x、y轴上投影:将上式移项整理得:在(3.42)式中,与s是未知量,联立得:其中:注意,上式中,则s有两个共轭复根,表明此RRP级杆组在机构中不能装配。所以在计算s值之前,应检验的值;当
23、机构中顺时针排列时根号前取“”,反之取“”。将(3.44)式求得的s代入(3.41)式得: b、速度分析:将(3.41)式对时间求导,得:将上式移项整理得:在上式中,与是未知量,联立得: c、加速度分析:将(3-47)式对时间求导,得:将上式移项整理得:在上式中,与是未知量,联立得:4)、RPR级杆组xlirBrDjjyOCDlk图3.16 RPR级杆组sEBljrC如图3-16所示的RPR级杆组,它是由两个构件,两个转动副和一个移动副组成。已知:各构件的长度li、lj、lK,B点和D点的位置及运动参数。求:C点的运动参数,构件j的位置角jj,角速度wj,角加速度aj和E点的运动参数。a、位置
24、分析:C点的矢量方程为:将上述矢量方程在x、y轴上投影得:将上式移项整理得:在上式中,s与jj是未知量,联立得:当BCD顺时针排列时,根号前取“”,反之取“”;为保证正确装配,上式中根号内必须大于零。将(3.61)代入(3.59)得:(3.62) (3.63)E点的矢量方程为: (3.64)将矢量式向X、Y轴投影得: b、速度分析将(3.59)式对时间求导并整理得: (3.66)上式中和wj为未知量,联立求得: (3.67) (3.68)C点速度,对(3.58)求导得: (3.69)E点速度,对(3.65)求导得: (3.70) c、加速度分析将(3-66)式对时间求导并整理得:(3.71)在
25、上式中,和aj为未知量,联立求得:(3.72)(3.73)其中,将(3.69)式对时间求导得C点加速度为:(3.74)将(3.70)式对时间求导得E点加速度为:(3.75)将上述各常见杆组所推导的方程,即数学模型编成各相应杆组的子程序,然后在所要求解的机构中调用即可。进行运动仿真建模的过程在对机构进行运动解析建模设计时,用各种工程设计软件都可以对机构进行运动分析。利用ADAMS高级工程软件建模设计既方便又直观。下面简介一下利用ADAMS软件进行机构参数化建模的方法。所谓参数化建模,是将机构放在直角坐标系下,将机构初始位置中各点的x、y坐标用各杆件长度和角度的方程式来表达,并将各表达式写在参数表
26、中,以便在主界面中建立各点。然后连接各点来创建构件,并在各构件间创建运动副,给机构加上驱动后,机构就可进行仿真运动。在ADAMS的后处理器中,可很方便地看到机构中各点的位移、速度和加速度曲线,及各构件的角位移、角速度和角加速度曲线。并通过曲线中各点值的大小对机构进行对比分析。若要认为何处需要改动,即可将参数表调出,改动构件的长短或原动件驱动量的大小。将改变参数后的机构再次仿真运行后,就可以看到改变后曲线数值的大小。用ADAMS软件对机构进行运动分析是很方便的。 下面以摆动导杆机构为例说明用ADAMS软件建模的过程。HYXP1P4P3P6P5P7P10P9P8yP2AB曲柄滑块1导杆滑块2连杆图
27、3.17 摆动导杆机构OCDna已知:机构中各杆的长度,原动件曲柄的转速n,方向如图所示。求:滑块2的位移、速度和加速度曲线。1、 建立参数化的数学模型将机构放在直角坐标系下,在机构的左极限位置建模。机构中各个点都用Pi(i=1、2、10)表示,为了描述滑块1、滑块2的大小,在A点和C点分别用几个点来表示其尺寸。设滑块1用圆柱体来表示,其长度lP4P5=50;滑块2用长方体表示,其x和y方向的尺寸分别是60和40。如图所示: 机构中各点的坐标为:2、建立参数表,将上述各个点与机构的已知运动参数写进参数表中。3、创建构件,并在构件与构件之间加运动副,在原动件上加驱动。然后进行运动仿真。4、机构运
28、动仿真后,在后处理器中调出运动曲线。以上是用ADAMS软件运动仿真建模的过程,利用建好的模型还可以进一步进行动力学分析。在进行受力分析时,在已知条件中要加入各个构件的质量和转动惯量,及外力。同样要将这些参数写入参数表,以便进行受力分析时用。3.3 机构的力分析3.3.1机构的动态静力分析机构的动态静力分析主要是分析机构在刚体运动范围内受力的情况,也是为将来需要分析机构在高速运转下弹性体受力分析打基础。机构的动态静力分析主要是求解各机构联接处的低副内的支反力和原动件上的平衡力(或平衡力矩)。每个低副中的反力都有两个未知的要素,如转动副中的反力有方向和大小未知,移动副中的反力有大小和作用点未知。若
29、一个机构中有PL个低副,其反力的未知数有2PL,加上一个未知的平衡力(或平衡力矩),故机构中未知数的个数为2PL1,而每个活动构件受力平衡时,可以列出三个平衡方程,即:,若机构中有n个活动构件,则可列出3n个力平衡方程式。为了使机构受力可解,其未知数的个数应等于所列的力平衡方程数,即:(3.76)此式即为自由度F=1的平面机构的自由度计算公式,即:(3.77)说明自由度为1的平面机构受力是静定可解的。本节用解析法中的矩阵法讨论这种自由度为1的机构的受力分析。xyOC(xc ,yc)R23xA(xA ,yA)图3.18 构件受力平衡图2R23yP2xR12yR12xP2yB(xB,yB)M2Mb
30、在讨论机构受力分析前,先做如下规定:(1)力矩:逆时针转向为正;顺时针转向为负。(2)未知铰链点的力向X、Y方向上的投影量、方向:均为X、Y轴的正向为正;反之为负。(3)在机构中,因为,为了求解方便,统一用表示,所以可用表示。下面举例说明。例8中构件2,该构件分别在A点和B点与1构件和3构件相联。已知:作用在2构件质心C点的力和力矩有、和及各点位置。求:平衡力矩和各运动副反力。解:各运动副中的力如图3.18所示:(3.78)将上式中已知量放在方程右侧,未知量放在等式左侧整理得: 对于机构中每一个构件都可以这样分析,然后联立求解。例3.9如图3.19所示的四杆机构。已知:各铰链点及质心点的坐标值,每个构件上作用于质心点Ci(i=1,2,3)的外力分别为Pix ,Piy(包括重力、惯性力和工作阻力),外力矩分别为Mi,3构件上
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 二零二四年度项目进度对赌协议3篇
- 生猪栏位设计与建设技术转让协议04
- 2024年度东莞城市亮化工程项目合同2篇
- 药品购销合同模板完整版
- 2024年度店面物业:租赁合同与物业服务协议2篇
- 铝合金门窗材料购销合同范本 2篇
- 2024年版豪华住宅区公寓销售协议3篇
- 2024年度环境监理工程师服务合同3篇
- 开展青年员工培训提升
- 人教版九年级化学第二单元我们周围的空气2氧气课时1氧气的性质教学教学课件
- 《我的白鸽》课件
- 国开2024年《中国法律史》平时作业1-3答案
- 大学生职业规划大赛成长赛道参赛作品
- GB 17790-2008家用和类似用途空调器安装规范
- 五年级上册数学课件 -《平行四边形的面积》 人教版(共15张PPT)
- PSPC作业指导书用书
- 安培环路定理及应用PPT课件
- TS3000在航天炉煤化工中应用
- 关于公司办公用房超标的情况报告
- 生产车间环境管理制度
- 120个文言实词小故事全文翻译
评论
0/150
提交评论