• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于應(yīng)力約束的框架結(jié)構(gòu)拓?fù)鋬?yōu)化研究

    2021-03-10 00:37:24何一凡趙磊
    河南科技 2021年28期
    關(guān)鍵詞:拓?fù)鋬?yōu)化框架結(jié)構(gòu)

    何一凡 趙磊

    摘 要:針對(duì)工程框架結(jié)構(gòu)的強(qiáng)度需求,建立了應(yīng)力和體積約束下框架結(jié)構(gòu)柔順度最小化的拓?fù)鋬?yōu)化模型。首先,為解決應(yīng)力優(yōu)化過程中的應(yīng)力奇異和大量局部約束問題,利用qp應(yīng)力松弛技術(shù)和p范數(shù)凝聚函數(shù)法構(gòu)建了應(yīng)力約束的歸一化等效應(yīng)力約束方案。其次,提出了基于凝聚應(yīng)力和變體積約束限措施的修正方案,以克服優(yōu)化過程中應(yīng)力約束嚴(yán)重非線性和最大局部應(yīng)力波動(dòng)等問題。再次,導(dǎo)出了目標(biāo)函數(shù)和應(yīng)力約束的靈敏度公式,并采用移動(dòng)漸近線方法算法進(jìn)行優(yōu)化求解。最后,給出優(yōu)化算例,驗(yàn)證了本文方法的正確性與可行性。

    關(guān)鍵詞:拓?fù)鋬?yōu)化;框架結(jié)構(gòu);應(yīng)力約束;應(yīng)力松弛;凝聚函數(shù)法

    中圖分類號(hào):TB21 文獻(xiàn)標(biāo)識(shí)碼:A 文章編號(hào):1003-5168(2021)28-00-07

    Abstract: For the strength requirements of the engineering frame structures, this article establishes a topology optimization model that minimizes the flexibility of the frame structure under stress and volume constraints. In order to solve the stress singularity and deal with a large number of local constraints during the stress optimization process, a normalized equivalent stress constraint scheme for stress constraints is constructed by using the qp stress relaxation technique and the p-norm aggregation function method. Subsequently, a modified scheme based on the aggregation stress and the variable volume constraint limit measures is proposed to overcome the nonlinear problem of aggregation stress constraints and the maximum local stress fluctuation problem during an optimization process. The sensitivity formulas of the objective function and stress constraints are derived, and the MMA algorithm is adopted to optimize the model. Finally, the optimization example given verifies the correctness and feasibility of the proposed method.

    Keywords: topology optimization;frame structure;stress constraint;stress relaxation;aggregate function method

    拓?fù)鋬?yōu)化作為結(jié)構(gòu)創(chuàng)新設(shè)計(jì)的重要手段,已在機(jī)械、土木及航空航天等工程領(lǐng)域得到了廣泛應(yīng)用。目前,大多數(shù)研究仍關(guān)注連續(xù)體結(jié)構(gòu)剛度相關(guān)的優(yōu)化設(shè)計(jì)[1-4],而強(qiáng)度問題是工程結(jié)構(gòu)設(shè)計(jì)需要考慮的首要條件之一。因此,考慮強(qiáng)度的結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)獲得了越來越多的關(guān)注[5-8]。

    近年來,涉及強(qiáng)度需求的連續(xù)體結(jié)構(gòu)拓?fù)鋬?yōu)化設(shè)計(jì)方法取得了較大進(jìn)展,但由于應(yīng)力約束的局部特性,優(yōu)化模型中存在大量的局部約束導(dǎo)致優(yōu)化求解困難,限制了相關(guān)研究的發(fā)展。為此,許多學(xué)者采用p范數(shù)或K-S(Kreisselmeier-Steinhauser)凝聚函數(shù)法對(duì)應(yīng)力約束進(jìn)行等效化處理[9],以解決應(yīng)力優(yōu)化相關(guān)的大規(guī)模局部約束問題。然而,由于應(yīng)力約束本身具有非線性特征,加之結(jié)構(gòu)中尖端、凹槽等區(qū)域出現(xiàn)的應(yīng)力集中現(xiàn)象和應(yīng)力約束凝聚處理導(dǎo)致了應(yīng)力約束的強(qiáng)非線性問題,使得優(yōu)化結(jié)果對(duì)優(yōu)化模型及優(yōu)化求解器的參數(shù)高度敏感。為了規(guī)避優(yōu)化結(jié)構(gòu)中的應(yīng)力集中現(xiàn)象和解決具有強(qiáng)非線性優(yōu)化模型的求解問題,需要對(duì)凝聚參數(shù)合理取值,并修正優(yōu)化模型,以自適應(yīng)方式抑制設(shè)計(jì)變量大的變化。同時(shí),應(yīng)力水平的精確評(píng)估與基于梯度的優(yōu)化算法的選用對(duì)克服該強(qiáng)非線性問題至關(guān)重要[10-11]。為此,張維聲等采用兩種全局應(yīng)力評(píng)估方法(包括應(yīng)力梯度和局部曲率信息)來控制結(jié)構(gòu)的局部應(yīng)力[12]。XIA等采用一種自適應(yīng)的有限元技術(shù)來避免優(yōu)化過程中出現(xiàn)的人工弱材料,以精確計(jì)算局部應(yīng)力水平[13]。PICELLI等提出一種應(yīng)力平均技術(shù)以提高應(yīng)力求解精確性[14],并采用自適應(yīng)歸一化方法以準(zhǔn)確約束局部應(yīng)力。此外,優(yōu)化過程中低密度區(qū)域會(huì)引起應(yīng)力奇異現(xiàn)象,進(jìn)一步導(dǎo)致應(yīng)力相關(guān)優(yōu)化求解存在較大困難。針對(duì)該問題,許多學(xué)者提出了ε松弛和qp應(yīng)力松弛技術(shù),并成功應(yīng)用到了應(yīng)力相關(guān)的拓?fù)鋬?yōu)化中[11]。但上述研究仍局限于連續(xù)體結(jié)構(gòu)的優(yōu)化設(shè)計(jì),而工程中由梁、桿等組成的框架結(jié)構(gòu)應(yīng)用普遍。為推進(jìn)工程框架結(jié)構(gòu)的實(shí)際應(yīng)用,需要進(jìn)一步發(fā)展考慮應(yīng)力約束的框架結(jié)構(gòu)的拓?fù)鋬?yōu)化方法。

    考慮到工程中框架結(jié)構(gòu)的強(qiáng)度需求,筆者提出了考慮應(yīng)力及變體積限約束的框架結(jié)構(gòu)柔順度最小化的拓?fù)鋬?yōu)化設(shè)計(jì)方法。為規(guī)避優(yōu)化結(jié)構(gòu)中的細(xì)長桿件,首先基于Heaviside映射函數(shù)構(gòu)建了設(shè)計(jì)變量與物理變量的合理映射關(guān)系。其次,采用qp應(yīng)力松弛技術(shù)和p范數(shù)凝聚函數(shù)對(duì)應(yīng)力約束進(jìn)行等效化處理,以克服應(yīng)力奇異問題和大量局部約束問題。再次,基于變體積限約束方法和凝聚函數(shù)的修正方案,建立了考慮應(yīng)力約束的框架結(jié)構(gòu)拓?fù)鋬?yōu)化的等效優(yōu)化模型。最后,通過移動(dòng)漸近線方法(Method of Moving Asymptotes,MMA)算法進(jìn)行了優(yōu)化求解,并通過兩個(gè)數(shù)值優(yōu)化算例說明了優(yōu)化設(shè)計(jì)方法的正確性與可行性。

    1 框架結(jié)構(gòu)的懲罰模型和等效應(yīng)力

    1.1 框架結(jié)構(gòu)的等效應(yīng)力

    選取研究對(duì)象為由空心矩形截面的Euler-Bernoulli梁單元[如圖1(a)所示]組成的框架結(jié)構(gòu)。

    對(duì)框架結(jié)構(gòu)而言,其構(gòu)件截面上應(yīng)力與截面中點(diǎn)的位置坐標(biāo)有關(guān)。將如圖1(b)所示的梁單元端面黑色標(biāo)注點(diǎn)作為該截面的臨危點(diǎn)[15],采用臨危點(diǎn)應(yīng)力校核結(jié)構(gòu)強(qiáng)度設(shè)計(jì)的安全性。若b、h和t分別表示梁單元橫截面的寬度、高度及厚度參數(shù),則單元橫截面積可表示為:

    在如圖1(a)所示的局部坐標(biāo)系(x,y,z)下,梁單元的內(nèi)力可表示為:

    式中:Ui、Ki和Fi分別為單元整體坐標(biāo)系中第i號(hào)單元節(jié)點(diǎn)位移矢量、在局部坐標(biāo)系下第i號(hào)單元的固有剛度矩陣及其單元內(nèi)力;Ti為局部坐標(biāo)系與整體坐標(biāo)系之間的變換矩陣,由單元局部坐標(biāo)軸和結(jié)構(gòu)整體坐標(biāo)軸的方向余弦構(gòu)成。在局部坐標(biāo)系下,第i號(hào)梁單元兩端截面內(nèi)力可以表示為:

    式中:qi1和qi2分別為第i號(hào)梁單元兩端截面的內(nèi)力;r11=[I6×6 06×6],r2=[06×6 I6×6],I6×6和06×6分別為維度為6×6的單位矩陣及零矩陣;Fx,ij表示第i號(hào)梁單元j(j=1,2)截面的軸力;Fy,ij和Fz,ij分別為相應(yīng)截面平行于y和z軸的剪力;Mx,ij、My,ij和Mz,ij分別為相應(yīng)截面受到的繞x軸、y軸和z軸的彎矩。

    梁單元截面上任意一點(diǎn)的應(yīng)力狀態(tài)為:

    其中:

    式中:σF是軸向力Fx引起的正應(yīng)力;σMy、σMz是彎矩My和Mz引起的彎曲應(yīng)力;τMy是由扭矩Mx引起的扭轉(zhuǎn)切應(yīng)力;τFy、τFz是由剪力Fy和Fz引起的彎曲切應(yīng)力。

    第i號(hào)梁單元第j截面上臨危點(diǎn)的Von Mises等效應(yīng)力可表示為:

    1.2 框架結(jié)構(gòu)的懲罰模型

    采用固體各向同性材料懲罰模型(Solid Isotropic Material with Penalization,SIMP)插值模型,第i號(hào)單元的剛度矩陣Ki表示為:

    式中:Ki0為第i號(hào)單元局部坐標(biāo)系中固有單元?jiǎng)偠染仃?p為懲罰參數(shù)。為了解決結(jié)構(gòu)單元灰度問題,采用Heaviside映射技術(shù),形成設(shè)計(jì)變量與物理變量ρi映射的表達(dá)式:

    式中:β為Heaviside映射曲率參數(shù)。借鑒相關(guān)文獻(xiàn)[11]的應(yīng)力奇異處理措施,可將單元體積vi、截面應(yīng)力σ(jvm,t)和應(yīng)力限分別表示為:

    式中:vi0和表示第i號(hào)單元的固有體積和j(j=1,2)截面臨危點(diǎn)固有應(yīng)力;αv、p和q為懲罰參數(shù);為結(jié)構(gòu)固有的最大許用應(yīng)力。

    2 優(yōu)化模型

    考慮到工程框架結(jié)構(gòu)的強(qiáng)度需求,可構(gòu)建應(yīng)力及體積約束的框架結(jié)構(gòu)柔順度最小化的拓?fù)鋬?yōu)化模型,如式(17)所示:

    式中:是拓?fù)湓O(shè)計(jì)變量向量,m為結(jié)構(gòu)中梁單元的數(shù)目;為第l組載荷工況下的自適應(yīng)權(quán)重系數(shù);Ctol為組合的目標(biāo)柔順度;Cl為第l組載荷工況下的柔順度,由式(19)給出;K為整體剛度矩陣;Ul為第l組載荷工況下整體坐標(biāo)系中結(jié)構(gòu)的位移矢量;Fl為第l個(gè)載荷工況;n為載荷工況數(shù);ρi和vi0分別為第i個(gè)物理變量和第i號(hào)單元的初始體積;V(0)為初始迭代時(shí)結(jié)構(gòu)的總體積;為在第l組載荷工況下第i號(hào)單元的j號(hào)截面的Von Mises等效應(yīng)力值;和分別為拓?fù)湓O(shè)計(jì)變量上下限,本文取;V*為目標(biāo)體積。自適應(yīng)權(quán)重系數(shù)的表達(dá)式為:

    式中:Cl(w)表示第w迭代步第l(l=1,2,…,n)組載荷工況下的結(jié)構(gòu)柔順度;Cl為第l組載荷工況下的結(jié)構(gòu)柔順度,其可表示為:

    為了解決大量局部應(yīng)力約束導(dǎo)致的計(jì)算量大的困難,基于qp應(yīng)力松弛技術(shù)和p范數(shù)凝聚函數(shù)法,將優(yōu)化模型式(17)中的大量局部應(yīng)力約束近似等效為下列凝聚應(yīng)力約束[15]。

    式中:可以通過凝聚參數(shù)pn(pn>1)的值使得第l組載荷工況下的結(jié)構(gòu)凝聚應(yīng)力約束近似等效于第l組載荷工況下結(jié)構(gòu)大量局部應(yīng)力約束;表示第l組載荷工況下的結(jié)構(gòu)第i號(hào)單元的j(j=1,2)截面臨危點(diǎn)固有應(yīng)力;為第w迭代步的值。但是,pn取值過大將使凝聚應(yīng)力約束的非線性程度增大,會(huì)導(dǎo)致優(yōu)化求解過程反復(fù)震蕩;當(dāng)pn取值較小時(shí),凝聚函數(shù)無法取代包絡(luò)對(duì)象的最大值。相關(guān)參考文獻(xiàn)[16]引入系數(shù)cp,使得式(20)為緊約束,其中:

    3 靈敏度分析

    3.1 結(jié)構(gòu)柔順度及體積的靈敏度

    由優(yōu)化模型式(17)的結(jié)構(gòu)柔順度表達(dá)式可推出目標(biāo)函數(shù)Ctol關(guān)于物理變量ρk的靈敏度:

    式中:Ul(w)為第l組載荷工況下整體坐標(biāo)系中結(jié)構(gòu)第w步的位移矢量。

    3.2 應(yīng)力約束的靈敏度

    由鏈?zhǔn)角髮?dǎo)法則可得凝聚應(yīng)力約束對(duì)物理變量的靈敏度為:

    基于上述靈敏度分析,結(jié)合MMA算法,對(duì)本文所提出的優(yōu)化模型進(jìn)行優(yōu)化求解,當(dāng)滿足如式(37)和式(38)所示的收斂條件時(shí)則終止迭代,得到優(yōu)化解。為了進(jìn)一步解決嚴(yán)重的非線性求解問題,在MMA迭代求解過程中采用xml對(duì)設(shè)計(jì)變量進(jìn)行控制。

    式中:m的初始值為1,并且其值每50個(gè)外循環(huán)迭代步增加1。

    式中:和分別為第w步和第w-1步設(shè)計(jì)變量的值;md為灰度指標(biāo);ε1、ε2為兩個(gè)較小的經(jīng)驗(yàn)參數(shù)值,本文取ε1=0.01,ε2=0.01。

    4 算例分析

    一高聳空心結(jié)構(gòu)高80 m,最底部外尺寸為4 m×4 m的正方形區(qū)域,最頂部外尺寸為直徑4 m的圓形區(qū)域,采用TOMáS等[17]提出的方式生成該模型。中間空心區(qū)域是底面半徑為3 m、高80 m的圓柱形區(qū)域。初始框架基結(jié)構(gòu)中,梁構(gòu)件的截面均為0.1 m×0.1 m且壁厚0.01 m的空心矩形截面。結(jié)構(gòu)底面固定且受到如圖2所示的兩組方向相反的1 000 kN外載荷作用,其中F1沿y軸的反方向(工況1),F(xiàn)2沿y軸的正方向(工況2)。結(jié)構(gòu)的彈性模量E=210 GPa,泊松比v=0.3,結(jié)構(gòu)的初始最大應(yīng)力為3.974 1 MPa。本算例中取SIMP懲罰參數(shù)為p=3,應(yīng)力約束限的懲罰參數(shù)q=2.5,凝聚參數(shù)pn=8,目標(biāo)體積為0.06,變體積限的移動(dòng)步長為0.02。

    結(jié)構(gòu)的初始基結(jié)構(gòu)設(shè)計(jì)域如圖3所示,由1 656個(gè)梁單元構(gòu)成。采取4種不同的應(yīng)力約束方案。方案一=3.2 MPa,方案二=3.5 MPa,方案三=5 MPa,方案四未考慮應(yīng)力約束。

    如圖2所示,在工況1和工況2的共同作用下,結(jié)構(gòu)在不同應(yīng)力約束限下的優(yōu)化拓?fù)浼捌鋺?yīng)力云圖分別如圖4至圖11所示??芍蚣芙Y(jié)構(gòu)的優(yōu)化拓?fù)潆S應(yīng)力約束限的不同產(chǎn)生細(xì)微變化,而相應(yīng)的應(yīng)力云圖隨應(yīng)力約束限的不同產(chǎn)生較大變化。由圖4至圖11的對(duì)比分析還可以看出:有無考慮應(yīng)力約束的優(yōu)化結(jié)果存在明顯不同。

    結(jié)構(gòu)在不同約束限下的柔順度、最大Von Mises應(yīng)力及結(jié)構(gòu)體積的迭代歷程曲線分別如圖12至圖15所示,可知:優(yōu)化過程穩(wěn)定收斂,且其優(yōu)化結(jié)構(gòu)的柔順度和最大應(yīng)力都隨桿件減少而增大。對(duì)應(yīng)圖4至圖9中優(yōu)化構(gòu)型的最大應(yīng)力分別為2.880 9 MPa、2.6345 MPa及2.529 7 MPa,其最終柔順度分別為10.284 6 N·m、11.231 6 N·m和10.435 1 N·m。圖10和圖11中未考慮應(yīng)力約束的優(yōu)化構(gòu)型的最大應(yīng)力為3.014 2 MPa,其最終柔順度為11.088 7 N·m。由此可知,優(yōu)化結(jié)構(gòu)的最大應(yīng)力也隨不同應(yīng)力約束限而不同,且優(yōu)化結(jié)構(gòu)的最大應(yīng)力隨應(yīng)力約束限的增大而增大,而其柔順度未發(fā)生明顯變化。同時(shí),考慮應(yīng)力約束的優(yōu)化結(jié)構(gòu)的最大應(yīng)力明顯小于未考慮應(yīng)力約束的優(yōu)化結(jié)果??梢?,提出的方法有效控制了優(yōu)化結(jié)構(gòu)的最大應(yīng)力,得到的拓?fù)鋬?yōu)化結(jié)果是滿足應(yīng)力約束且合理的。

    5 結(jié)論

    針對(duì)框架結(jié)構(gòu)的強(qiáng)度需求,基于SIMP插值模型和Heaviside映射函數(shù)建立了應(yīng)力及體積約束下框架結(jié)構(gòu)輕量化拓?fù)鋬?yōu)化模型。采用qp應(yīng)力松弛技術(shù)和p范數(shù)凝聚函數(shù)法構(gòu)建了等效應(yīng)力約束函數(shù),通過凝聚應(yīng)力約束修正方案和變體積約束限方法形成了相應(yīng)的等效優(yōu)化模型。此外,進(jìn)一步解決了嚴(yán)重的非線性模型求解問題,對(duì)每一步設(shè)計(jì)變量的變化范圍進(jìn)行了控制。最后,通過數(shù)值算例來說明所提出優(yōu)化方法的正確性與可行性,并得到如下結(jié)論。

    ①本方法能獲得穩(wěn)定收斂的解,且能得到滿足應(yīng)力約束的清晰優(yōu)化拓?fù)洹?/p>

    ②結(jié)構(gòu)的優(yōu)化構(gòu)型隨應(yīng)力約束限的不同而發(fā)生變化,且其柔順度和最大應(yīng)力都隨桿件的消除而增大。

    ③考慮應(yīng)力約束的優(yōu)化結(jié)構(gòu)的最大應(yīng)力明顯低于未考慮應(yīng)力約束的最大應(yīng)力。提出的方法有效控制了優(yōu)化結(jié)構(gòu)的最大應(yīng)力,得到的拓?fù)鋬?yōu)化結(jié)果滿足應(yīng)力約束且合理。

    參考文獻(xiàn):

    [1]DUYSINX P, BENDSOE M P. Topology optimization of continuum structures with stress constraints[J]. International Journal for Numerical Methods in Engineering, 1998(8):1453-1478.

    [2]ZHU J H, GUO W J , ZHANG W H, et al. Integrated layout and topology optimization design of multi-frame and multi-component fuselage structure systems[J]. Structural and Multidisciplinary Optimization, 2017(1):21-45.

    [3]冷國俊, 張卓, 保宏, 等.考慮重疊過濾及穩(wěn)定性約束的桁架拓?fù)鋬?yōu)化方法[J].工程力學(xué), 2013(2):8-13.

    [4]CUI H, AN H, HUANG H . Truss topology optimization considering local buckling constraints and restrictions on intersection and overlap of bar members[J]. Structural and Multidisciplinary Optimization, 2018(2):575-594.

    [5]CHANGIZI N, KABOO DA NIAN H, JALALPOUR M. Stress-based topology optimization of frame structures under geometric uncertainty[J]. Computer Methods in Applied Mechanics and Engineering, 2017, 315:121-140.

    [6]JEONG S H, CHOI D H, YOON G H . Separable stress interpolation scheme for stress-based topology optimization with multiple homogenous materials[J]. Finite Elements in Analysis & Design, 2014(5):16-31.

    [7]GAO X J, LI YX, MA H T, et al. Improving the overall performance of continuum structures: a topology optimization model considering stiffness, strength and stability[J]. Computer Methods in Applied Mechanics and Engineering, 2020,359:112660.

    [8]LE C, NORATO J, BRUNS T, et al. Stress-based topology optimization for continua[J]. Structural and Multidisciplinary Optimization, 2010(4):605-620.

    [9]王選, 劉宏亮, 龍凱, 等.基于改進(jìn)的雙向漸進(jìn)結(jié)構(gòu)優(yōu)化法的應(yīng)力約束拓?fù)鋬?yōu)化[J].力學(xué)學(xué)報(bào), 2018(2):385-394.

    [10]PARíS J, NAVARRINA F, COLOMINAS I, et al. Topology optimization of continuum structures with local and global stress constraints[J]. Chinese Journal of Theoretical & Applied Mechanics, 2009(4):98-104.

    [11]MOON S, YOON G. A newly developed qp-relaxation method for element connectivity parameterization to achieve stress-based topology optimization for geometrically nonlinear structures [J]. Computer Methods in Applied Mechanics and Engineering, 2013(3):226-241.

    [12]ZHANG W, GUO X, WANG M, et al. Optimal topology design of continuum structures with stress concentration alleviation via level set method [J]. International Journal for Numerical Methods in Engineering, 2013(9):942-959.

    [13]XIA Q, SHI T, LIU S, et al. A level set solution to the stress-based structural shape and topology optimization [J]. Computers & Structures, 2012(1):55-64.

    [14]PICELLI R, TOWNSEND S, BRAMPTON C, et al. Stress-based shape and topology optimization with the level set method [J]. Computer Methods in Applied Mechanics and Engineering, 2018(1):1-23.

    [15]ZUO W, YU J, SAITOU K. Stress sensitivity analysis and optimization of automobile body frame consisting of rectangular tubes[J]. International Journal of Automotive Technology, 2016(5):843-851.

    [16]龍凱, 王選, 吉亮.面向應(yīng)力約束的獨(dú)立連續(xù)映射方法[J].力學(xué)學(xué)報(bào), 2019(2):320-329.

    [17]TOMáS ZEGARD, PAULINO G H . GRAND3— Ground structure based topology optimization for arbitrary 3D domains using MATLAB[J]. Structural & Multidisciplinary Optimization, 2015(6):1161-1184.

    3424500338229

    猜你喜歡
    拓?fù)鋬?yōu)化框架結(jié)構(gòu)
    無黏結(jié)預(yù)應(yīng)力框架結(jié)構(gòu)的拆改加固設(shè)計(jì)
    混凝土框架結(jié)構(gòu)抗震加固方法簡述
    建筑結(jié)構(gòu)設(shè)計(jì)中框架結(jié)構(gòu)設(shè)計(jì)的應(yīng)用
    基于負(fù)載均衡的電力通信流量分析和預(yù)測方法研究
    軌道扣件彈性墊板結(jié)構(gòu)優(yōu)化設(shè)計(jì)
    電推衛(wèi)星方圓過渡結(jié)構(gòu)輕量化設(shè)計(jì)
    基于OptiStruct的襟翼搖臂的拓?fù)鋬?yōu)化
    科技視界(2016年22期)2016-10-18 15:37:54
    基于ANSYS的多層框架結(jié)構(gòu)隔震性能分析
    利用HyperWorks進(jìn)行助力缸支架優(yōu)化設(shè)計(jì)
    某型飛機(jī)操縱搖臂拓?fù)鋬?yōu)化研究
    科技視界(2015年25期)2015-09-01 16:31:59
    成人18禁高潮啪啪吃奶动态图 | 秋霞伦理黄片| 九九在线视频观看精品| 18禁裸乳无遮挡动漫免费视频| 亚洲情色 制服丝袜| 最近2019中文字幕mv第一页| 成年女人在线观看亚洲视频| 寂寞人妻少妇视频99o| 少妇人妻久久综合中文| av电影中文网址| 99久久人妻综合| 十八禁网站网址无遮挡| 国产熟女午夜一区二区三区 | 99re6热这里在线精品视频| 亚洲一区二区三区欧美精品| av一本久久久久| 久久久久精品性色| 欧美一级a爱片免费观看看| 黑丝袜美女国产一区| 欧美日韩精品成人综合77777| 中国国产av一级| 熟女人妻精品中文字幕| 另类亚洲欧美激情| 在线观看免费日韩欧美大片 | 日韩熟女老妇一区二区性免费视频| 黄片无遮挡物在线观看| 精品少妇久久久久久888优播| 一级黄片播放器| 国产淫语在线视频| 亚洲一区二区三区欧美精品| 丰满乱子伦码专区| 最近2019中文字幕mv第一页| 尾随美女入室| 街头女战士在线观看网站| av一本久久久久| 精品人妻熟女av久视频| 亚洲精华国产精华液的使用体验| 91国产中文字幕| 久久午夜综合久久蜜桃| 飞空精品影院首页| 免费av不卡在线播放| 久久久久精品久久久久真实原创| 中文字幕av电影在线播放| 国产精品久久久久久久久免| 人人妻人人澡人人爽人人夜夜| 国产欧美日韩综合在线一区二区| 国产成人a∨麻豆精品| 国产爽快片一区二区三区| 国产伦理片在线播放av一区| 26uuu在线亚洲综合色| 日本av手机在线免费观看| 日本黄色片子视频| 欧美日本中文国产一区发布| 草草在线视频免费看| 国产免费视频播放在线视频| 久久久久久人妻| 女人久久www免费人成看片| 欧美日韩精品成人综合77777| 91在线精品国自产拍蜜月| 国产精品久久久久久精品古装| 欧美激情 高清一区二区三区| 这个男人来自地球电影免费观看 | videosex国产| 成人国语在线视频| 亚洲四区av| 国产成人精品一,二区| 肉色欧美久久久久久久蜜桃| 国产精品人妻久久久影院| 一本一本综合久久| 亚洲图色成人| 国产精品不卡视频一区二区| 观看av在线不卡| 国产成人91sexporn| av国产精品久久久久影院| 欧美日韩亚洲高清精品| 精品国产一区二区三区久久久樱花| 美女国产视频在线观看| 丝袜脚勾引网站| 少妇的逼水好多| 永久免费av网站大全| 九色成人免费人妻av| 2021少妇久久久久久久久久久| 亚洲国产毛片av蜜桃av| 亚洲av欧美aⅴ国产| 国国产精品蜜臀av免费| 亚洲精品第二区| 午夜免费男女啪啪视频观看| 3wmmmm亚洲av在线观看| videossex国产| 新久久久久国产一级毛片| 欧美日韩一区二区视频在线观看视频在线| 18禁裸乳无遮挡动漫免费视频| 免费日韩欧美在线观看| 中国美白少妇内射xxxbb| 中文欧美无线码| 天美传媒精品一区二区| 国产精品蜜桃在线观看| 大陆偷拍与自拍| 日韩中文字幕视频在线看片| 丝袜脚勾引网站| 91精品一卡2卡3卡4卡| av播播在线观看一区| 天堂8中文在线网| 亚洲国产成人一精品久久久| 日本黄大片高清| 精品久久久久久久久av| 蜜桃国产av成人99| 你懂的网址亚洲精品在线观看| xxx大片免费视频| 久久久久精品性色| www.av在线官网国产| 亚洲美女视频黄频| 一本一本综合久久| 精品少妇黑人巨大在线播放| 久久久久久伊人网av| 国内精品宾馆在线| 99久久精品一区二区三区| 2022亚洲国产成人精品| 嫩草影院入口| 久久久欧美国产精品| 欧美日韩在线观看h| 啦啦啦视频在线资源免费观看| 久久久久久久久久久丰满| 国产熟女欧美一区二区| 成人漫画全彩无遮挡| 特大巨黑吊av在线直播| 国产精品久久久久久久电影| 亚洲欧洲日产国产| 国产欧美日韩一区二区三区在线 | 97精品久久久久久久久久精品| 色5月婷婷丁香| 69精品国产乱码久久久| 99九九在线精品视频| 夜夜爽夜夜爽视频| 亚洲天堂av无毛| 精品国产国语对白av| 国产一区有黄有色的免费视频| 母亲3免费完整高清在线观看 | 国产一区二区三区av在线| 91久久精品国产一区二区成人| 又黄又爽又刺激的免费视频.| 母亲3免费完整高清在线观看 | 日韩伦理黄色片| 能在线免费看毛片的网站| 亚洲丝袜综合中文字幕| 精品国产乱码久久久久久小说| 成人无遮挡网站| 黄色欧美视频在线观看| 看非洲黑人一级黄片| 亚洲高清免费不卡视频| kizo精华| 99久久中文字幕三级久久日本| a级毛片黄视频| 日韩在线高清观看一区二区三区| 最黄视频免费看| 一区二区三区乱码不卡18| 久久精品熟女亚洲av麻豆精品| 只有这里有精品99| 少妇猛男粗大的猛烈进出视频| 久久久a久久爽久久v久久| 亚洲精品自拍成人| 亚洲国产最新在线播放| 久久av网站| 亚洲欧美成人综合另类久久久| 夜夜骑夜夜射夜夜干| 考比视频在线观看| 日日撸夜夜添| 午夜激情福利司机影院| 国产精品一区二区在线观看99| 自拍欧美九色日韩亚洲蝌蚪91| 韩国高清视频一区二区三区| 男女免费视频国产| 欧美日韩视频高清一区二区三区二| 欧美成人午夜免费资源| 久久热精品热| 一本久久精品| 亚洲精品美女久久av网站| 亚洲精品乱久久久久久| 国产欧美日韩一区二区三区在线 | 妹子高潮喷水视频| 欧美日韩一区二区视频在线观看视频在线| 午夜免费鲁丝| 一级毛片aaaaaa免费看小| 亚洲色图 男人天堂 中文字幕 | 99热全是精品| 日本-黄色视频高清免费观看| 久久ye,这里只有精品| 99久久人妻综合| 男女啪啪激烈高潮av片| av国产精品久久久久影院| 一级爰片在线观看| 亚洲美女视频黄频| 超碰97精品在线观看| 九九在线视频观看精品| av在线播放精品| 免费少妇av软件| 国产精品成人在线| 免费观看的影片在线观看| 欧美老熟妇乱子伦牲交| 国产成人免费观看mmmm| 国产一级毛片在线| 亚洲精品亚洲一区二区| 日本欧美国产在线视频| 国产欧美亚洲国产| 久久ye,这里只有精品| 久久精品国产鲁丝片午夜精品| 成人毛片a级毛片在线播放| 精品少妇久久久久久888优播| 国产精品麻豆人妻色哟哟久久| 91国产中文字幕| 亚洲精品国产av成人精品| 蜜桃在线观看..| 又大又黄又爽视频免费| 久久99热6这里只有精品| 久久国产精品大桥未久av| 国产视频首页在线观看| videos熟女内射| av播播在线观看一区| 我要看黄色一级片免费的| 久久久久久久久久成人| 日韩av免费高清视频| 18禁观看日本| 国产精品成人在线| 精品人妻熟女av久视频| 肉色欧美久久久久久久蜜桃| 亚洲成人手机| 十八禁网站网址无遮挡| 久久 成人 亚洲| 久久这里有精品视频免费| 视频区图区小说| 狂野欧美白嫩少妇大欣赏| 国产白丝娇喘喷水9色精品| 免费久久久久久久精品成人欧美视频 | 有码 亚洲区| 国产欧美日韩一区二区三区在线 | 在线亚洲精品国产二区图片欧美 | 欧美 日韩 精品 国产| 天天操日日干夜夜撸| 中文字幕av电影在线播放| 精品人妻在线不人妻| 国产一区二区在线观看av| 麻豆乱淫一区二区| 夜夜看夜夜爽夜夜摸| 国产成人av激情在线播放 | 日本wwww免费看| 免费黄色在线免费观看| 制服丝袜香蕉在线| 99久久综合免费| 亚洲无线观看免费| 91久久精品国产一区二区三区| av在线app专区| 国产精品无大码| 一本大道久久a久久精品| 亚洲在久久综合| 国产男人的电影天堂91| 免费黄频网站在线观看国产| 亚洲美女视频黄频| 中文字幕制服av| 亚洲一区二区三区欧美精品| av女优亚洲男人天堂| 日韩熟女老妇一区二区性免费视频| 中文字幕最新亚洲高清| 久久这里有精品视频免费| 99热网站在线观看| 黑人巨大精品欧美一区二区蜜桃 | 国产成人av激情在线播放 | 亚洲精品乱久久久久久| 97在线人人人人妻| 亚洲欧洲国产日韩| 欧美日韩亚洲高清精品| 国产 一区精品| 97精品久久久久久久久久精品| 夜夜看夜夜爽夜夜摸| 亚洲欧美日韩卡通动漫| 国产精品不卡视频一区二区| 激情五月婷婷亚洲| 国产淫语在线视频| 久久久国产一区二区| 精品人妻一区二区三区麻豆| 亚洲av成人精品一二三区| 在线观看www视频免费| 我的老师免费观看完整版| 国产免费一区二区三区四区乱码| 久久鲁丝午夜福利片| 国产视频首页在线观看| 男人添女人高潮全过程视频| xxx大片免费视频| 亚洲国产精品国产精品| 又大又黄又爽视频免费| 在线亚洲精品国产二区图片欧美 | 久久精品国产鲁丝片午夜精品| 男人操女人黄网站| 另类精品久久| 欧美日韩成人在线一区二区| 51国产日韩欧美| 亚洲欧美精品自产自拍| 日韩免费高清中文字幕av| 九草在线视频观看| 国产高清有码在线观看视频| 国产精品99久久久久久久久| 亚洲情色 制服丝袜| 80岁老熟妇乱子伦牲交| 色哟哟·www| 亚洲av成人精品一区久久| 搡女人真爽免费视频火全软件| 999精品在线视频| 国产黄色免费在线视频| 成年人午夜在线观看视频| 精品酒店卫生间| 丰满乱子伦码专区| 国产精品国产av在线观看| 大片电影免费在线观看免费| 夜夜爽夜夜爽视频| 国产国拍精品亚洲av在线观看| 乱码一卡2卡4卡精品| 久久精品久久久久久噜噜老黄| 在线观看www视频免费| 男女边摸边吃奶| 在现免费观看毛片| 一级毛片黄色毛片免费观看视频| 最近手机中文字幕大全| 中文字幕人妻丝袜制服| 亚洲综合精品二区| 夜夜看夜夜爽夜夜摸| 午夜91福利影院| av有码第一页| 日日摸夜夜添夜夜爱| 精品少妇黑人巨大在线播放| 久久精品久久精品一区二区三区| 人妻制服诱惑在线中文字幕| 国产成人免费无遮挡视频| 国产成人freesex在线| 18禁在线播放成人免费| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 午夜福利在线观看免费完整高清在| videossex国产| 国产白丝娇喘喷水9色精品| 久久精品久久精品一区二区三区| 九草在线视频观看| 男男h啪啪无遮挡| 人妻制服诱惑在线中文字幕| 久久国产亚洲av麻豆专区| 亚洲久久久国产精品| 久久久精品94久久精品| 亚洲人成网站在线观看播放| 中文字幕最新亚洲高清| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久久亚洲精品成人影院| 曰老女人黄片| 亚洲精品自拍成人| a级毛色黄片| 久久99精品国语久久久| 成人二区视频| 99热国产这里只有精品6| 嘟嘟电影网在线观看| 成人手机av| 韩国高清视频一区二区三区| 亚洲熟女精品中文字幕| 纯流量卡能插随身wifi吗| 国产熟女午夜一区二区三区 | 搡女人真爽免费视频火全软件| 少妇 在线观看| av.在线天堂| 成人亚洲欧美一区二区av| 人妻一区二区av| 国产国语露脸激情在线看| 黄色视频在线播放观看不卡| 午夜免费男女啪啪视频观看| 最新中文字幕久久久久| 国产 精品1| av有码第一页| 一级二级三级毛片免费看| 欧美一级a爱片免费观看看| 日韩不卡一区二区三区视频在线| 青春草亚洲视频在线观看| 亚洲不卡免费看| 黄片无遮挡物在线观看| 王馨瑶露胸无遮挡在线观看| 午夜影院在线不卡| 亚洲精品久久久久久婷婷小说| 久久久久人妻精品一区果冻| 国产午夜精品一二区理论片| 久久久精品94久久精品| 91精品一卡2卡3卡4卡| 激情五月婷婷亚洲| 午夜免费男女啪啪视频观看| 国产高清国产精品国产三级| 午夜福利,免费看| 成人18禁高潮啪啪吃奶动态图 | 全区人妻精品视频| 老司机亚洲免费影院| 亚洲第一区二区三区不卡| videossex国产| 国产男女超爽视频在线观看| 久久精品国产自在天天线| 99久国产av精品国产电影| 九九在线视频观看精品| 久久国产亚洲av麻豆专区| 五月开心婷婷网| 久久97久久精品| 九草在线视频观看| 好男人视频免费观看在线| 黑丝袜美女国产一区| 丰满少妇做爰视频| 日本免费在线观看一区| 人妻人人澡人人爽人人| 国产免费福利视频在线观看| 人妻少妇偷人精品九色| 亚洲精品亚洲一区二区| 婷婷成人精品国产| 人人妻人人澡人人看| 自线自在国产av| 少妇人妻 视频| 亚洲国产欧美日韩在线播放| 又大又黄又爽视频免费| 精品亚洲乱码少妇综合久久| av免费在线看不卡| 国产日韩欧美在线精品| 多毛熟女@视频| 久久久欧美国产精品| 中国国产av一级| 国产伦理片在线播放av一区| av在线播放精品| 少妇高潮的动态图| 国产精品久久久久久av不卡| 精品国产露脸久久av麻豆| 亚洲欧美清纯卡通| 欧美xxⅹ黑人| 狂野欧美白嫩少妇大欣赏| 2022亚洲国产成人精品| 国产精品久久久久久久久免| 亚洲怡红院男人天堂| 老熟女久久久| 国产精品国产三级国产av玫瑰| 又黄又爽又刺激的免费视频.| 大片电影免费在线观看免费| 久久av网站| 啦啦啦视频在线资源免费观看| 大陆偷拍与自拍| 亚洲精品日本国产第一区| 亚洲精品亚洲一区二区| 18禁动态无遮挡网站| 亚洲怡红院男人天堂| 午夜影院在线不卡| 国产日韩一区二区三区精品不卡 | 国产精品一区二区在线不卡| 少妇高潮的动态图| 国产精品熟女久久久久浪| 美女大奶头黄色视频| 欧美丝袜亚洲另类| www.色视频.com| 大又大粗又爽又黄少妇毛片口| 男女高潮啪啪啪动态图| 国产欧美亚洲国产| 最后的刺客免费高清国语| 亚洲av日韩在线播放| 久久精品国产亚洲网站| 80岁老熟妇乱子伦牲交| 亚洲av成人精品一区久久| 美女国产视频在线观看| 老熟女久久久| 又大又黄又爽视频免费| 纯流量卡能插随身wifi吗| 亚洲色图 男人天堂 中文字幕 | 午夜久久久在线观看| 精品亚洲乱码少妇综合久久| 3wmmmm亚洲av在线观看| 亚洲精品第二区| 男女啪啪激烈高潮av片| 丝袜脚勾引网站| 国产成人精品无人区| 亚洲婷婷狠狠爱综合网| 国产白丝娇喘喷水9色精品| 最近中文字幕高清免费大全6| a级片在线免费高清观看视频| 日本猛色少妇xxxxx猛交久久| 日韩视频在线欧美| 亚洲内射少妇av| 亚洲av成人精品一区久久| 亚洲精品乱久久久久久| 91成人精品电影| 一级毛片 在线播放| 三上悠亚av全集在线观看| 男人操女人黄网站| 欧美丝袜亚洲另类| 成人手机av| 妹子高潮喷水视频| 欧美日韩视频精品一区| 亚洲五月色婷婷综合| 午夜福利网站1000一区二区三区| 久久国内精品自在自线图片| 亚洲第一区二区三区不卡| 国内精品宾馆在线| 亚洲国产精品一区二区三区在线| 色婷婷久久久亚洲欧美| 99视频精品全部免费 在线| 亚洲欧美中文字幕日韩二区| 亚洲色图 男人天堂 中文字幕 | kizo精华| 午夜91福利影院| 91午夜精品亚洲一区二区三区| 亚洲国产欧美日韩在线播放| 五月玫瑰六月丁香| 精品少妇黑人巨大在线播放| 国内精品宾馆在线| 制服诱惑二区| 午夜影院在线不卡| 黄色怎么调成土黄色| 欧美日韩一区二区视频在线观看视频在线| 国产精品国产三级国产专区5o| 精品久久久久久电影网| 少妇的逼水好多| av卡一久久| 男的添女的下面高潮视频| 亚洲国产欧美日韩在线播放| 啦啦啦在线观看免费高清www| 国产午夜精品一二区理论片| 大片免费播放器 马上看| a级毛片在线看网站| 亚洲人成网站在线播| 久久久国产一区二区| 伊人久久精品亚洲午夜| 国产老妇伦熟女老妇高清| 亚洲精品日韩av片在线观看| 精品久久久噜噜| 免费大片18禁| 免费观看在线日韩| 热re99久久精品国产66热6| 九九久久精品国产亚洲av麻豆| av视频免费观看在线观看| 少妇猛男粗大的猛烈进出视频| 男的添女的下面高潮视频| 18禁观看日本| 久久久精品94久久精品| 久久久国产精品麻豆| 人妻少妇偷人精品九色| 亚洲av成人精品一二三区| 80岁老熟妇乱子伦牲交| 国产免费福利视频在线观看| 9色porny在线观看| 肉色欧美久久久久久久蜜桃| 亚洲精品中文字幕在线视频| xxx大片免费视频| 国产在线视频一区二区| 曰老女人黄片| 亚洲欧美成人综合另类久久久| 欧美精品高潮呻吟av久久| 免费大片18禁| 满18在线观看网站| 特大巨黑吊av在线直播| 国产亚洲午夜精品一区二区久久| 色94色欧美一区二区| 亚洲人与动物交配视频| 美女国产视频在线观看| 高清视频免费观看一区二区| 一本一本综合久久| 日韩伦理黄色片| 考比视频在线观看| 一区二区三区乱码不卡18| av播播在线观看一区| 日本欧美视频一区| 91在线精品国自产拍蜜月| av网站免费在线观看视频| 黑人巨大精品欧美一区二区蜜桃 | 国产亚洲一区二区精品| 亚洲一级一片aⅴ在线观看| 91精品三级在线观看| 五月天丁香电影| 99热国产这里只有精品6| 亚洲人成网站在线观看播放| 一本久久精品| 欧美人与性动交α欧美精品济南到 | 色哟哟·www| 欧美最新免费一区二区三区| 久久人人爽人人片av| 一级毛片 在线播放| 国产毛片在线视频| 老司机亚洲免费影院| 又大又黄又爽视频免费| 亚洲精品一二三| 国产精品熟女久久久久浪| 久久久久精品久久久久真实原创| 另类亚洲欧美激情| 免费观看的影片在线观看| 最黄视频免费看| 午夜免费男女啪啪视频观看| 中文字幕精品免费在线观看视频 | 久久 成人 亚洲| 国产黄色免费在线视频| 一级爰片在线观看| 在线观看免费高清a一片| 18在线观看网站| av在线老鸭窝| 啦啦啦视频在线资源免费观看| 久久久久久久久久久丰满| 国产成人av激情在线播放 | 亚洲欧美清纯卡通| 69精品国产乱码久久久| 国产午夜精品一二区理论片| 中文字幕精品免费在线观看视频 | 18禁在线无遮挡免费观看视频| 2022亚洲国产成人精品| 久久国产精品男人的天堂亚洲 | 在线观看人妻少妇| 亚洲,一卡二卡三卡| 亚洲精品亚洲一区二区| av线在线观看网站| 综合色丁香网| 男女边摸边吃奶| 在线观看一区二区三区激情| 国产成人a∨麻豆精品| 91国产中文字幕| 午夜免费观看性视频| 最近手机中文字幕大全| 国产极品粉嫩免费观看在线 |