版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡介
第30卷第6期陜西科技大學(xué)學(xué)報(bào)Vol.30No.62012年12月JournalofShaanxiUniversityofScience&TechnologyDec.2012*文章編號(hào):1000-5811(2012)06-0093-04壓力容器的有限元分析及優(yōu)化設(shè)計(jì)張彩麗(陜西科技大學(xué)機(jī)電工程學(xué)院,陜西西安710021)摘要:針對傳統(tǒng)設(shè)計(jì)中壓力容器的材料浪費(fèi)問題,提出基于有限元分析的壓力容器的優(yōu)化設(shè)計(jì)方法.該方法首先利用有限元法對壓力容器進(jìn)行分析并提取分析結(jié)果中的相關(guān)參數(shù),然后利用優(yōu)化設(shè)計(jì)方法進(jìn)行定量計(jì)算,最終得到既滿足性能指標(biāo)又滿足設(shè)計(jì)指標(biāo)的設(shè)計(jì)參數(shù).實(shí)際結(jié)果表明,基于有限元分析的優(yōu)化設(shè)計(jì)方法,在工程實(shí)際應(yīng)用中可有效發(fā)揮作用.關(guān)鍵詞:優(yōu)化設(shè)計(jì);有限元分析;壓力容器中圖法分類號(hào):TH122文獻(xiàn)標(biāo)識(shí)碼:AFiniteelementanalysisandoptimizationdesignofthepressurevesselZHANGCai-li(CollegeofMechanicalandElectricalEngineering,ShaanxiUniversityofScience&Technology,Xi′an710021,China)Abstract:Accordingtopressurevesselmaterialwasteprobleminthetraditionaldesign,thispaperputsforwardthepressurevesseloptimizationdesignmethodbasedonfiniteelementa-nalysis.Inthismethod,firstlythefiniteelementanalysisisappliedtocarryoutstresscalcu-lation,andtherelatedresultsparametersareextractedtobeusefollowingcalculation.Thenthequantitativecalculationisconductedbyusingoptimizationdesignmethod,andthebestdesignparameterswhichmeetperformanceindexesaregot.Atlast,thismethodwasusedfortheoptimizationdesignofpressurevessels.Practicalresultsprovedthevalidityandtheprac-ticabilityofthismethodinthepressurevesselsdesign.Keywords:optimizationdesign;finiteelementanalysis;pressurevessel0引言壓力容器是化工、煉油、機(jī)械中廣泛使用的承壓設(shè)備,在傳統(tǒng)設(shè)計(jì)中,考慮到壓力容器安全問題的重要性,壓力容器的設(shè)計(jì)往往偏于保守,設(shè)計(jì)的壓力容器既笨重又浪費(fèi)材料,制造成本明顯高,尤其是隨著壓力容器越來越趨于大型化設(shè)計(jì),材料的浪費(fèi)現(xiàn)象愈加嚴(yán)重[1].為此,設(shè)計(jì)出既滿足性能要求又節(jié)約材料的壓力容器就成為生產(chǎn)制造企業(yè)追求的目標(biāo).由于壓力容器的實(shí)際結(jié)構(gòu)一般都比較復(fù)雜,對其進(jìn)行解析求解較困難,同時(shí)要求設(shè)計(jì)人員應(yīng)具有扎實(shí)的理論基礎(chǔ),故在壓力容器設(shè)計(jì)的過程中,最有效、最實(shí)用的方法就是數(shù)值分析的方法.文中結(jié)合有限元分析的特點(diǎn)及傳統(tǒng)優(yōu)化設(shè)計(jì)的不足,利用大型通用有限元分析軟件ANSYS對壓力容器進(jìn)行有限元分析,然后提取有限元分析結(jié)果進(jìn)行優(yōu)化設(shè)計(jì),既滿足了壓力容器的性能要求,又達(dá)到了節(jié)能降耗、降低成本和價(jià)格,提升企業(yè)競爭力的目的[2,3].*收稿日期:2012-10-08基金項(xiàng)目:陜西省教育廳科學(xué)研究計(jì)劃項(xiàng)目(12JK0690)作者簡介:張彩麗(1973-),女,陜西渭南人,副教授,研究方向:現(xiàn)代設(shè)計(jì)理論及方法陜西科技大學(xué)學(xué)報(bào)第30卷1優(yōu)化設(shè)計(jì)基本原理優(yōu)化設(shè)計(jì)是近年來發(fā)展起來的一門新的學(xué)科,它在解決復(fù)雜問題時(shí),能定量地從眾多的設(shè)計(jì)方案中找到盡可能完美的或最適宜的設(shè)計(jì)方案,故在工程實(shí)際中的應(yīng)用越來越廣泛.優(yōu)化設(shè)計(jì)是數(shù)學(xué)規(guī)劃和計(jì)算機(jī)技術(shù)相結(jié)合的產(chǎn)物,是一種將設(shè)計(jì)變量表示為產(chǎn)品性能指標(biāo)、結(jié)構(gòu)指標(biāo)或運(yùn)動(dòng)參數(shù)指標(biāo)的函數(shù),稱為目標(biāo)函數(shù);然后在產(chǎn)品規(guī)定的性態(tài)、幾何和運(yùn)動(dòng)等其它條件的限制范圍內(nèi),稱為約束條件;尋找一個(gè)或多個(gè)目標(biāo)函數(shù)最大或最小的設(shè)計(jì)變量組合的數(shù)學(xué)方法.進(jìn)行優(yōu)化設(shè)計(jì)時(shí),首先要把實(shí)際設(shè)計(jì)問題轉(zhuǎn)化為優(yōu)化設(shè)計(jì)的數(shù)學(xué)模型.在明確設(shè)計(jì)變量、約束條件、目標(biāo)函數(shù)之后,優(yōu)化設(shè)計(jì)問題數(shù)學(xué)模型的一般形式為:求設(shè)計(jì)變量x=[x1x2…xn]T使目標(biāo)函數(shù)f(x)的值最?。恚椋睿妫ǎ┣覞M足約束條件gj(x)≤0,j=1,2,…,mhk(x)≤0,k=1,2,…,phl(x)=0,l=1,2,…,q式中:n-設(shè)計(jì)變量的個(gè)數(shù);m-性能約束條件的個(gè)數(shù);p-幾何約束條件的個(gè)數(shù);q-設(shè)計(jì)變量之間的約束條件個(gè)數(shù).2有限元法進(jìn)行優(yōu)化設(shè)計(jì)的基本過程利用大型通用有限元分析軟件ANSYS進(jìn)行優(yōu)化設(shè)計(jì),可按照以下4個(gè)步驟進(jìn)行:(1)建立參數(shù)化的有限元分析文件有限元分析文件的建立在整個(gè)優(yōu)化設(shè)計(jì)中具有重要的意義,分析文件建立的正確與否會(huì)影響最終的優(yōu)化設(shè)計(jì)結(jié)果.該文件的建立可采用兩種方法:一種是直接編輯法,另一種是基于ANSYS的交互式方法.建立參數(shù)化的有限元分析文件包括以下內(nèi)容:單元類型的選擇、實(shí)常數(shù)的輸入、材料特性參數(shù)的選擇、實(shí)體模型的建立、對實(shí)體模型的網(wǎng)格劃分即有限元模型的建立、分析類型的選擇、約束條件及載荷的確定、求解以及對分析結(jié)果中相關(guān)數(shù)據(jù)的提?。ǎ玻└鶕?jù)求解問題,在ANSYS數(shù)據(jù)庫中建立與分析文件中的變量相對應(yīng)的設(shè)計(jì)參數(shù).(3)執(zhí)行優(yōu)化計(jì)算執(zhí)行優(yōu)化設(shè)計(jì)計(jì)算時(shí),首先進(jìn)入ANSYS優(yōu)化設(shè)計(jì)的模塊,指定已經(jīng)建立的分析文件;然后聲明優(yōu)化設(shè)計(jì)變量及其取值范圍、狀態(tài)變量及其取值范圍,并選擇目標(biāo)函數(shù),即確定有限元優(yōu)化設(shè)計(jì)的數(shù)學(xué)模型;接著選擇優(yōu)化設(shè)計(jì)工具或優(yōu)化設(shè)計(jì)方法、指定優(yōu)化循環(huán)的控制方式;最后進(jìn)行優(yōu)化求解.(4)查看、選取并檢驗(yàn)優(yōu)化設(shè)計(jì)結(jié)果通過有限元優(yōu)化設(shè)計(jì)會(huì)得到一系列可行的和不可行的設(shè)計(jì)方案,設(shè)計(jì)者需要從這些方案中選出最好的設(shè)計(jì)方案,同時(shí)檢驗(yàn)優(yōu)化設(shè)計(jì)結(jié)果的合理性.3壓力容器的有限元優(yōu)化設(shè)計(jì)3.1問題描述圖1所示為一用20鋼制造的壓力容器,根據(jù)工廠生產(chǎn)的要求,該壓力容器的最大內(nèi)壓為pmax=15MPa,選用鋼板的厚度為h=3mm,壓力容器的幾何尺寸R和H滿足下列關(guān)系:H-R≥30mm,材料的屈服極限為245MPa,彈性模量E=206GPa,泊松比μ=0.3.要求:在鋼板厚度不變的情況下,確定壓力容器具有最大體積時(shí)所對應(yīng)的幾何尺寸R和H的大?。畧D1壓力容器結(jié)構(gòu)圖3.2壓力容器優(yōu)化設(shè)計(jì)數(shù)學(xué)模型的建立根據(jù)上述問題描述,選取幾何尺寸R和H作為設(shè)計(jì)變量,以壓力容器的最大體積作為目標(biāo)函數(shù).結(jié)合壓力容器結(jié)構(gòu)特點(diǎn),可得目標(biāo)函數(shù)的表達(dá)式如下:Vmax=3R3+2πR2(H-R)該壓力容器結(jié)構(gòu)及載荷滿足對稱性要求,為提高后續(xù)有限元優(yōu)化設(shè)計(jì)計(jì)算的效率,取壓力容器的四分之一進(jìn)行有限元優(yōu)化設(shè)計(jì).同時(shí),利用AN-·49·第6期張彩麗:壓力容器的有限元分析及優(yōu)化設(shè)計(jì)SYS軟件進(jìn)行有限元優(yōu)化設(shè)計(jì)時(shí),一般是求解目標(biāo)函數(shù)的最小值,與上述目標(biāo)函數(shù)求解壓力容器的最大體積不相符,需要把求解目標(biāo)函數(shù)的最大值轉(zhuǎn)化為求解目標(biāo)函數(shù)的最小值.根據(jù)優(yōu)化設(shè)計(jì)理論,有兩種轉(zhuǎn)化方式:一種是對原目標(biāo)函數(shù)取倒數(shù),另一種是對原目標(biāo)函數(shù)取負(fù)值.文中對原目標(biāo)函數(shù)取負(fù)值,把求解目標(biāo)函數(shù)的最大值轉(zhuǎn)化為求解目標(biāo)函數(shù)的最小值.得到壓力容器優(yōu)化設(shè)計(jì)目標(biāo)函數(shù)的表達(dá)式為:V=-3R3-2R2(H-R)=-6R2(3H-R)=V(R,H)設(shè)定設(shè)計(jì)變量的取值范圍,考慮性能約束條件及給定條件,最終建立的壓力容器優(yōu)化設(shè)計(jì)的數(shù)學(xué)模型如下:minV(X),X=[RH]Ts.t.g1(X)=R>0g2(X)=H>0g3(X)=σmax-245≥0g4(X)=H-R≥303.3壓力容器有限元模型的建立及其分析首先,考慮到后續(xù)進(jìn)行優(yōu)化設(shè)計(jì)時(shí)設(shè)計(jì)變量隨優(yōu)化結(jié)果的變化而處于動(dòng)態(tài)變化中,對壓力容器半徑R、容器高度相關(guān)尺寸H、容器的體積V在AN-SYS軟件中建立參數(shù),初步?。遥剑常埃龋剑福埃浯?,針對壓力容器的薄壁結(jié)構(gòu)及結(jié)構(gòu)和受力具有軸對稱的特點(diǎn),按殼體結(jié)構(gòu)建立參數(shù)化的實(shí)體模型[4].第三,在ANSYS軟件中先建立1/4圓柱面,再建立1/8球面,兩者經(jīng)布爾運(yùn)算得到壓力容器的1/8實(shí)體模型[5],如圖2所示.圖2壓力容器的參數(shù)化模型針對建立的壓力容器實(shí)體模型,選用單元類型shell63,以厚度h=3為實(shí)常數(shù),輸入材料特性參數(shù):彈性模量E=206GPa,泊松比μ=0.3,對其進(jìn)行網(wǎng)格劃分,建立壓力容器的有限元分析模型[6],如圖3所示.圖3壓力容器的有限元模型對上述有限元模型,在其邊界線上施加對稱的位移約束,并對壓力容器的內(nèi)表面施加壓力載荷,如圖4所示,進(jìn)行有限元求解.求解出壓力容器沿x、y、z方向及總的位移云圖,如圖5所示,壓力容器沿x、y、z方向及總的應(yīng)力云圖,如圖6所示.最后,根據(jù)有限元分析結(jié)果,建立單元列表,提取壓力容器上的最大應(yīng)力,以便后續(xù)優(yōu)化設(shè)計(jì)中使用.圖4壓力容器的加載圖從求解結(jié)果發(fā)現(xiàn):該壓力容器的最大變形量在頂部中間,最大變形量為0.039mm;x、y方向的最大應(yīng)力為163.8MPa,z方向的最大應(yīng)力為141.7MPa,總的最大應(yīng)力在筒體中間,最大應(yīng)力為222MPa,小于材料的許用應(yīng)力.3.4壓力容器的有限元優(yōu)化設(shè)計(jì)根據(jù)上述有限元模型建立、分析、求解及有限元應(yīng)力分析結(jié)果的提取過程,建立壓力容器有限元優(yōu)化設(shè)計(jì)的分析文件,并指定分析文件.依據(jù)R和·59·陜西科技大學(xué)學(xué)報(bào)第30卷H的初始數(shù)值,初步選擇設(shè)計(jì)變量R和H的取值范圍,且先設(shè)置參數(shù)R的取值范圍,然后再設(shè)置參數(shù)H的取值范圍.設(shè)置狀態(tài)變量即最大應(yīng)力的取值范圍,選擇體積V為目標(biāo)函數(shù),取一階優(yōu)化設(shè)計(jì)方法,設(shè)定循環(huán)控制方式,對其進(jìn)行優(yōu)化分析.圖5壓力容器的位移云圖圖6壓力容器的應(yīng)力云圖分析第一次優(yōu)化設(shè)計(jì)的結(jié)果發(fā)現(xiàn),設(shè)計(jì)變量R和H的取值范圍比較接近上限值,且最大應(yīng)力比許用應(yīng)力小,為使壓力容器的體積最大,還需進(jìn)一步優(yōu)化.故依據(jù)第一次優(yōu)化設(shè)計(jì)結(jié)果,調(diào)整各個(gè)設(shè)計(jì)變量的原始取值范圍,重新進(jìn)行優(yōu)化設(shè)計(jì).結(jié)果表明,通過再次優(yōu)化設(shè)計(jì)可進(jìn)一步改善優(yōu)化設(shè)計(jì)結(jié)果.經(jīng)過多次尋優(yōu)設(shè)計(jì),使設(shè)計(jì)變量的值不斷靠近最優(yōu)解,最終經(jīng)過29次迭代計(jì)算,得到最佳的設(shè)計(jì)結(jié)果為:R=55.896mm,H=99.549mm時(shí),壓力容器滿足強(qiáng)度要求,目標(biāo)函數(shù)最小,最小值為V=-3.96928×105mm3,經(jīng)轉(zhuǎn)換得到壓力容器的最大體積為Vmax=4×V=15.8768×105mm3.結(jié)合工程實(shí)際需要,對優(yōu)化后壓力容器的結(jié)構(gòu)尺寸進(jìn)行圓整,可?。遥剑担担恚恚龋剑保埃埃恚?,此時(shí)壓力容器的最大體積為Vmax=15.5144×105mm3.圖7迭代次數(shù)與設(shè)計(jì)變量之間的關(guān)系曲線圖8迭代次數(shù)與目標(biāo)函數(shù)之間的關(guān)系曲線從有限元優(yōu)化設(shè)計(jì)的結(jié)果,可獲知整個(gè)優(yōu)化設(shè)計(jì)的過程中設(shè)計(jì)變量、目標(biāo)函數(shù)隨迭代序列的變化曲線,如圖7、圖8所示.由圖可看出,剛開始迭代計(jì)算時(shí),設(shè)計(jì)變量變化范圍較大,對應(yīng)的目標(biāo)函數(shù)值也有較大變化,隨著迭代次數(shù)的增加,設(shè)計(jì)變量的變化趨于平穩(wěn),目標(biāo)函數(shù)也隨之趨于平穩(wěn),迭代進(jìn)行到第29次,得到壓力容器體積的最小值,即優(yōu)化設(shè)計(jì)問題的最優(yōu)解.4結(jié)論通過對壓力容器的有限元分析及優(yōu)化設(shè)計(jì)可得如下結(jié)論:(1)在設(shè)計(jì)變量的取值范圍不易確定的情況下,可根據(jù)初步的優(yōu)化設(shè)計(jì)結(jié)果進(jìn)行估算,然后逐步縮小設(shè)計(jì)變量的取值范圍,進(jìn)行多次優(yōu)化計(jì)算,可進(jìn)一步提高優(yōu)化設(shè)計(jì)的精度,使設(shè)計(jì)方案更符合實(shí)際的需要.(下轉(zhuǎn)第109頁)·69·第6期李頎等:多功能手持終端開發(fā)板的設(shè)計(jì)6.[2]王海玲,胡琨元,朱云龍,等.基于ARM9+Linux的RFID智能終端設(shè)計(jì)[J].微計(jì)算機(jī)信息,2008,24(11):19-21.[3]虞建平.智能手機(jī)中藍(lán)牙系統(tǒng)的研究和實(shí)現(xiàn)[D].長春:吉林大學(xué),2005.[4]馬忠梅.ARM嵌入式處理器結(jié)構(gòu)與應(yīng)用基礎(chǔ)[M].北京:北京航空航天大學(xué)出版社,2002.[5]王金莉,蘇宛新.基于PXA270的嵌入式系統(tǒng)設(shè)計(jì)[J].微計(jì)算機(jī)信息,2008,24(4):11-13.[6]鄭寧漢.基于PXA255的嵌入式實(shí)驗(yàn)裝置的設(shè)計(jì)與實(shí)現(xiàn)[J].計(jì)算機(jī)工程與設(shè)計(jì),2007,12(8):1882-1884.[7]楊建軍,高明煜,黃繼業(yè).基于PXA310的多功能數(shù)碼相框硬件平臺(tái)開發(fā)[J].杭州電子科技大學(xué)學(xué)報(bào),2009,12(1):9-12.[8]NationalSemiconductor.《PowerManagementUnitforAdvanceApplicationProcessors》內(nèi)部資料,2008,2.[9]邵常勇,陳滌,董國鋒.基于ARM的嵌入式系統(tǒng)設(shè)計(jì)方法研究[J].信息技術(shù)與信息化,2006,24(2):檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯檯96-98.(上接第96頁)(2)優(yōu)化設(shè)計(jì)變量的初始值選擇不同,會(huì)影響設(shè)計(jì)變量、目標(biāo)函數(shù)隨迭代序列的變化曲線,但不影響最終的優(yōu)化設(shè)計(jì)結(jié)果.(3)由最終設(shè)計(jì)序列中設(shè)計(jì)變量的變化關(guān)系,可發(fā)現(xiàn)在優(yōu)化設(shè)計(jì)中壓力容器半徑R的變化對體積的變化明顯,與目標(biāo)函數(shù)表達(dá)式中的關(guān)系相一致.(4)在有限元優(yōu)化設(shè)計(jì)中,設(shè)計(jì)變量、約束條件容差選擇不同,對最終的設(shè)計(jì)結(jié)果有一定的影響.(5)各設(shè)計(jì)變量、目標(biāo)函數(shù)隨迭代次數(shù)的增加均向最優(yōu)解逼近,說明了有限元分析法在優(yōu)化設(shè)計(jì)中的應(yīng)用價(jià)值.參考文獻(xiàn)[1]梁基照.壓力容器優(yōu)化設(shè)計(jì)[M].北京:機(jī)械工業(yè)出版社,2010.[2]李斌斌.基于有限元法的壓力容器壁厚優(yōu)化[J].輕工機(jī)械,2009,27(6):38-40.[3]任國棟,孔春元,任秀玲.基于ansys的壓力容器多變量優(yōu)化設(shè)計(jì)[J].機(jī)械工程與自動(dòng)化,2010,39(1):68-70.[4]劉伯玉,丁傳安.薄壁壓力容器的有限元建模研究[J].現(xiàn)代制造技術(shù)與裝備,2009,45(3):32-35.[5]張亞新,石傳美.基于ANSYS的壓力容器壁厚優(yōu)化設(shè)計(jì)[J].機(jī)械與電子,2009,27(8):57-60.[6]Peng-feiLIU.Optimaldesignofpressurevess
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲(chǔ)空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 重慶2025年重慶市涪陵區(qū)教育事業(yè)單位面向應(yīng)屆公費(fèi)師范生招聘75人筆試歷年參考題庫附帶答案詳解
- 部分學(xué)校上學(xué)期期中考試八年級(jí)語文試卷(PDF版無答案)
- 漯河2024年河南漯河市審計(jì)局事業(yè)單位引進(jìn)高層次人才2人筆試歷年參考題庫附帶答案詳解
- 河南2024年河南信陽師范大學(xué)招聘碩士研究生42人筆試歷年參考題庫附帶答案詳解
- 昭通2025年云南昭通綏江縣住建局招聘編外聘用人員筆試歷年參考題庫附帶答案詳解
- 2025年中國發(fā)動(dòng)機(jī)橫梁總成市場調(diào)查研究報(bào)告
- 2025年重氮項(xiàng)目可行性研究報(bào)告
- 2025至2031年中國鏡子玻璃行業(yè)投資前景及策略咨詢研究報(bào)告
- 2025至2031年中國計(jì)量儀器管理系統(tǒng)行業(yè)投資前景及策略咨詢研究報(bào)告
- 2025年真空浸漆烘干機(jī)項(xiàng)目可行性研究報(bào)告
- 2024至2030年中國中水回用行業(yè)發(fā)展監(jiān)測及投資戰(zhàn)略規(guī)劃報(bào)告
- NB/T 11430-2023煤礦TBM掘進(jìn)施工工藝要求
- 全國職業(yè)院校技能大賽培訓(xùn)課件
- 供貨送貨服務(wù)承諾書
- 福建省醫(yī)院大全
- GB/T 16659-2024煤中汞的測定方法
- 《Unit 10 You're supposed to shake hands》單元檢測題及答案
- TSN 解決方案白皮書
- 完整2024年開工第一課課件
- 貨運(yùn)車輛駕駛員安全培訓(xùn)內(nèi)容資料完整
- 風(fēng)神汽車4S店安全生產(chǎn)培訓(xùn)課件
評(píng)論
0/150
提交評(píng)論