非保守集中力作用下飽和多孔懸臂梁的非線性彎曲

楊 驍, 周冬華

(上海大學(xué)土木工程系,上海 200072)

非保守集中力作用下飽和多孔懸臂梁的非線性彎曲

楊 驍, 周冬華

(上海大學(xué)土木工程系,上海 200072)

基于孔隙流體僅沿梁軸向運(yùn)動(dòng)的微觀不可壓飽和多孔彈性梁大撓度彎曲數(shù)學(xué)模型,利用 Galerkin截?cái)喾?研究固定端不可滲透、自由度可滲透的飽和多孔懸臂彈性梁在自由端處承受突加非保守集中力作用下的擬靜態(tài)非線性彎曲問題,給出了梁彎曲時(shí)撓度、彎矩等隨時(shí)間的響應(yīng)以及沿梁軸線的分布.數(shù)值結(jié)果表明:當(dāng)載荷較小時(shí),非保守集中力、保守集中力以及線性小撓度理論的結(jié)果相差很小;當(dāng)載荷較大時(shí),非線性大撓度理論的結(jié)果小于相應(yīng)線性小撓度理論的結(jié)果,非保守力的結(jié)果大于相應(yīng)保守力的結(jié)果,且這種差異隨著載荷的增大而增大.同時(shí),在集中載荷突加于梁上時(shí),多孔彈性梁骨架最初不變形,但隨著時(shí)間的增加,梁的撓度逐漸增大,并最終趨于穩(wěn)態(tài)值,此時(shí)多孔梁骨架承擔(dān)全部的外載荷.

多孔介質(zhì)理論;飽和多孔彈性梁;大撓度;非保守力;Galerkin截?cái)喾?/p>

Abstract:Assumption made in themathematicalmodel for large deflection of microscopic incompressible saturated poroelastic beam is that fluid in porescan only flow in the direction of a deformed axis.Based on thismodel,nonlinear quasi-static bending of a saturated poroelastic cantilever beam with fixed end impermeable and free end permeable,with a suddenly applied nonconservative concentrated transverse load at its free end,is investigated with the Galerkin truncation method.The curves of deflections and bending moments of the beam skeleton are shown w ith respect to time and the beam axis.Numerical results show that,for small loads,there is little difference between results for the nonconservative and conservative concentrated load and between large nonlinear and small linear deflection models.For large loads,resultsof large nonlinear deflection model are smaller than those of small linear deflection model.Results for the nonconservative loads are larger than those for conservative loads.The difference increases with load.Further,when a concentrated load is suddenly applied,deformation of the beam skeleton does not instantly occur.Deflection of the poroelastic beam increases and finally converges to a stationary state in which the poroelastic beam skeleton takes the entire external load.

Key words:theory of porousmedia;saturated poroelastic beam;large deflection;nonconservative load;Galerkin truncation method

飽和多孔介質(zhì)是指連通孔隙中充滿液體的多孔連通介質(zhì),其結(jié)構(gòu)力學(xué)行為的研究在流-固耦合相互作用理論中具有重要的學(xué)術(shù)價(jià)值和廣泛的應(yīng)用背景.這方面的研究早期起源于生物力學(xué)中的若干問題[1-2].目前,相關(guān)研究成果已應(yīng)用于各種工程領(lǐng)域中,如建筑行業(yè)和傳熱傳質(zhì)工程等.基于 Biot理論,Schanz等[3]研究了飽和多孔柱的波動(dòng)問題,得到了時(shí)間區(qū)域上的精確解;Rémond等[4]研究了飽和多孔空心圓柱在軸向循環(huán)載荷作用下的穩(wěn)態(tài)響應(yīng);Leclaire等[5]建立了一個(gè)新的多孔彈性板模型,研究了孔隙度、滲透率等對(duì)多孔板共振頻率的影響.基于Nunziato-Cowin模型,B?rsan[6]建立了多孔彈性殼的非線性模型,證明了線性化模型初邊值問題解的唯一性.

對(duì)于各向同性可壓飽和多孔彈性梁,Nowinski等[7]、Zhang等[8]研究發(fā)現(xiàn),梁橫向的孔隙流體壓力梯度通常大于其軸向的壓力梯度.于是,一般情形下可忽略孔隙流體的軸向運(yùn)動(dòng).然而,在自然和工程應(yīng)用中存在一類多孔材料,其微觀結(jié)構(gòu)使得流體以軸向運(yùn)動(dòng)為主,可忽略其橫向的運(yùn)動(dòng),如植物的根莖、傳熱管道等.基于 Biot理論,Li等[9-11]建立了分析飽和多孔彈性梁、桿變形的線性和非線性數(shù)學(xué)模型,詳細(xì)研究了飽和多孔彈性梁的靜、動(dòng)態(tài)力學(xué)行為.實(shí)際上,一些情形下可將固相骨架和孔隙流體視為微觀不可壓材料.對(duì)此類問題,楊驍?shù)萚12]基于飽和多孔介質(zhì)理論的三維控制方程,建立了微觀不可壓飽和多孔梁、桿線性小變形的數(shù)學(xué)模型,并利用 Lap lace變換數(shù)值研究了飽和多孔懸臂梁的擬靜態(tài)和動(dòng)態(tài)彎曲.Yang等[13]建立了不可壓飽和多孔彈性梁非線性彎曲的數(shù)學(xué)模型,并分析了飽和多孔懸臂梁在端部保守集中載荷作用下的非線性響應(yīng).基于此模型,楊驍和李麗[14-15]研究了飽和多孔簡支彈性梁在均布載荷作用下的靜、動(dòng)力非線性彎曲.

基于飽和多孔彈性梁非線性彎曲數(shù)學(xué)模型,本工作研究了飽和多孔懸臂梁在始終垂直于梁軸線自由端集中力作用下的非線性擬靜態(tài)彎曲問題.首先,給出了以橫截面轉(zhuǎn)角θS為基本未知量的初邊值問題;其次,利用 Galerkin截?cái)喾ê统Ｎ⒎址匠痰?Runge-Kutta數(shù)值法,得到了不同載荷集度下梁骨架轉(zhuǎn)角、彎矩等隨時(shí)間的響應(yīng).數(shù)值結(jié)果表明,與小撓度線性理論的結(jié)果相比,大撓度非線性理論的結(jié)果偏小,非保守力的結(jié)果大于相應(yīng)保守力的結(jié)果,且大撓度理論結(jié)果趨于其穩(wěn)態(tài)值的時(shí)間明顯小于相應(yīng)的小撓度理論結(jié)果.另外,飽和多孔彈性梁的撓度呈現(xiàn)黏性蠕變特征,且其彎矩仍隨時(shí)間的增加而增大.這些結(jié)果有助于揭示相關(guān)工程結(jié)構(gòu)的力學(xué)行為.

1 數(shù)學(xué)模型

設(shè)長為 L、慣性矩為 I的流體飽和多孔彈性懸臂梁由不相溶的微觀不可壓流相 RF和微觀不可壓多孔彈性固相骨架 RS組成 (見圖1),其孔隙度為 φ,側(cè)表面不透水.

圖1 自由端承受非保守集中載荷的多孔懸臂梁Fig.1 Cantilever poroelastic beam subjected to a nonconservative concentrated load at the free end

由于梁骨架微觀結(jié)構(gòu)的特殊性,假定多孔彈性梁中的孔隙流體僅能沿梁的軸向方向運(yùn)動(dòng),且固定端不可滲透,自由端可滲透.若彈性梁骨架 RS關(guān)于其橫截面橫觀各向同性,t=0時(shí)刻在自由端處承受突加集中載荷 P的作用,且載荷 P始終垂直于梁的軸線.根據(jù)飽和多孔彈性梁非線性彎曲的一般數(shù)學(xué)模型,可得如下無量綱的初邊值問題[13-15]:

式中,uS(s,t)和 wS(s,t)分別為梁骨架軸線的無量綱水平位移和撓度,θS(s,t)為梁骨架橫截面的轉(zhuǎn)角,N(s,t)為無量綱軸力,Mp(s,t)為無量綱的孔隙流體壓力等效力偶,而 s為無量綱的軸線弧長[13].

初邊值問題(1)～(6)可進(jìn)一步簡化為

式中 ,θ1S=θS(1,t).

確定轉(zhuǎn)角θS(s,t)后,N(s,t),Mp(s,t),uS(s,t)和 wS(s,t)可分別由下式確定:

2 Galerk in截?cái)喾ê蛿?shù)值結(jié)果

顯然,求解初邊值問題 (7)的精確解是相當(dāng)困難的.為此,下面利用 Galerkin截?cái)喾ê?Runge-Kutta法數(shù)值求初邊值問題 (7)的近似解.作近似

顯然,θS(s,t)滿足初邊值問題 (7)中的邊界條件.將式 (10)代人式 (7)的方程中,利用 Galerkin截?cái)喾傻?/p>

式中,

利用四階 Runge-Kutta法,可容易數(shù)值求解初值問題 (12).這樣,由式 (8)和 (9)可數(shù)值得到多孔彈性懸臂梁骨架的撓度、水平位移、彎矩及其孔隙流體壓力等效力偶等.

當(dāng)集中力 P較小時(shí),即在小撓度情形下,始終垂直于梁軸線的非保守集中力導(dǎo)致的飽和多孔梁的彎曲,應(yīng)接近于保守的始終豎直向下集中力導(dǎo)致的多孔彈性梁的彎曲.這一性質(zhì)在一定程度上提供了校驗(yàn)近似表達(dá)式 (10)的有效性方法,而豎直集中力作用下飽和多孔彈性梁彎曲的精確解可由 Lap lace變換求得[15].表 1給出了當(dāng)載荷參數(shù) P=0.1時(shí),在不同時(shí)刻 t,多孔懸臂梁自由端處撓度 wS的精確解和近似解,其中“ES”和“AS”分別表示精確解和本工作的近似解.可見,在小撓度范圍內(nèi),近似表達(dá)式 (10)可給出滿意的結(jié)果.

圖2示出了不同載荷 P下,飽和多孔懸臂彈性梁自由端 (即 s=1)處無量綱撓度 wS隨無量綱時(shí)間t的變化曲線.圖中可見,撓度隨時(shí)間 t單調(diào)增大,并且當(dāng)載荷參數(shù) P較小,即當(dāng) P<1時(shí),大撓度和小撓度理論,以及始終豎直向下保守集中力引起的自由端撓度和始終垂直于梁軸線的非保守集中力引起的自由端撓度沒有明顯的差別.然而,隨著載荷 P的增加,相關(guān)結(jié)果相差越來越大,當(dāng) P=5時(shí),多孔懸臂梁自由端處保守力情形下,大撓度理論的穩(wěn)態(tài)撓度只有非保守力情形下穩(wěn)態(tài)撓度的 82.3%,并且此撓度只是小撓度理論穩(wěn)態(tài)撓度的 48.1%.當(dāng)載荷 P較大時(shí),線性小撓度模型的自由端處撓度將大于 1,這在幾何上是不可能的.同時(shí),始終垂直于梁軸線的非保守集中力引起的自由端撓度大于豎直向下保守集中力引起的自由端撓度,并且二者差異隨載荷 P的增大而增大.圖3給出了當(dāng) P=3時(shí),飽和多孔懸臂梁在不同時(shí)刻 t的撓曲線.

表 1 近似解與精確解的比較Table 1 Compar ison between the approximate and analytical solutions

圖2 不同載荷 P下,自由端撓度 wS隨時(shí)間 t的變化Fig.2 Deflections wSvs.the time t at the free end for d ifferent load P

圖3 當(dāng) P=3時(shí),不同時(shí)刻 t的撓曲線 wS分布Fig.3 Deflection curves wSat d ifferent time t when P=3

圖4 不同載荷 P下,固定端處彎矩 mSE隨時(shí)間 t的變化Fig.4 Bend ing m om en ts mSEvs.the tim e t at the f ixed end for d ifferent load P

圖5 當(dāng) P=3時(shí),不同時(shí)刻 t彎矩 mSE沿軸線 s的分布Fig.5 D istr ibution s of the bend ing m om en ts mSEalong the beam axis s at d ifferent time t when P=3

3 結(jié) 束 語

限制飽和多孔彈性梁中孔隙流體僅沿變形后梁軸線方向擴(kuò)散,本工作研究了固定端不可滲透、自由端可滲透的飽和多孔懸臂彈性梁在自由端非保守集中力作用下的彎曲問題.利用 Galerkin截?cái)喾?將問題轉(zhuǎn)化為以橫截面轉(zhuǎn)角為基本未知量的初值問題,并通過 Runge-Kutta法進(jìn)行數(shù)值計(jì)算,得到了飽和多孔懸臂彈性梁在不同載荷作用下的撓度、彎矩等隨時(shí)間的響應(yīng).通過比較多孔彈性梁在載荷較小時(shí)的研究結(jié)果與線性小撓度精確解,說明了本研究Galerkin截?cái)嘟品椒ǖ挠行?數(shù)值結(jié)果表明:相比于保守力作用下飽和多孔懸臂梁的彎曲,非保守集中力將引起較大的撓度和彎矩;多孔彈性梁的撓度呈現(xiàn)黏性蠕變特征,且其彎矩仍隨時(shí)間的增加而增大.這些結(jié)果有助于揭示相關(guān)工程結(jié)構(gòu)的力學(xué)行為.

[1] LANYON L E.Functional strain as a determinant for bone remodeling[J].Calcified Tissue International,1984,36(S1):S56-S61.

[2] COW IN S C. Bone poroelasticity [J]. Journal of Biomechanics,1999,32(3):217-238.

[3] SCHANZM,CHENG A H D.Transientwave propagation in a one-dimensional poroelastic column [J]. Acta Mechanica,2000,145(1):1-18.

[4] AGNèS R,SALAH N.Transverse isotropic poroelastic osteon model under cyclic loading[J].Mechanics Research Communications,2005,32(6):645-651.

[5] LECLA IRE P, HOROSHENKOV K V. Transverse vibrationsof a thin rectangular porousplate saturated by a fluid[J].Journal of Sound and Vibration,2001,247(1):1-18.

[6] B?RSAN M.On the theory of elastic shellsmade from a material w ith voids[J]. International Journal of Solids and Structures,2006,43(10):3106-3123.

[7] NOW INSKI JL,DAV IS C F.The flexural and torsion of bones viewed as anisotropic poroelastic bodies[J].International Journal of Engineering Science,1972,10(12):1063-1079.

[8] ZHANG D,COW IN S C.Oscillatory bending of a poroelastic beam [J].Journal of the Mechanics and Physicsof Solids,1994,42(10):1575-1599.

[9] L IL P,SCHULGASSER K,CEDERBAUM G.Theory of poroelastic beamsw ith axial diffusion[J].Journal of the Mechanics and Physics of Solids,1995,43(12):2023-2042.

[10] L I L P,SCHULGASSER K,CEDERBAUM G.Large deflection analysis of poroelastic beams [J].International Journal Non-linear Mechanics,1998,33(1):1-14.

[11] CEDERBAUM G,SCHULGASSER K,L IL P.Interesting behavior patternsof poroelastic beams and columns[J].International Journal of Solids and Structures,1998,35(34):4931-4943.

[12] 楊驍,李麗.不可壓飽和多孔彈性梁、桿動(dòng)力響應(yīng)的數(shù)學(xué)模型[J].固體力學(xué)學(xué)報(bào),2006,27(2):159-166.

[13] YANG X,WANG C.A nonlinearmathematicalmodel for large deflection of incompressible saturated poroelastic beams[J].ApplMath Mech:Engl Ed,2007,28(12):1587-1598.

[14] 李麗,楊驍.簡支飽和多孔彈性梁的非線性彎曲 [J].力學(xué)季刊,2007,28(1):86-91.

[15] 楊驍,李麗.簡支飽和多孔彈性梁的非線性動(dòng)力響應(yīng)[J].力學(xué)季刊,2008,29(1):132-136.

(編輯:劉志強(qiáng))

Nonlinear Bend ing of Saturated Poroelastic Cantilever Beam w ith Nonconservative Concentrated L oad

YANG Xiao, ZHOU Dong-hua
(Department of Civil Engineering,ShanghaiUniversity,Shanghai200072,China)

O 357.3

A

1007-2861(2010)03-0221-05

10.3969/j.issn.1007-2861.2010.03.001

2009-06-08

上海市自然科學(xué)基金資助項(xiàng)目 (06ZR14037);國家自然科學(xué)基金資助項(xiàng)目(10872124)

楊 驍 (1965～),男,教授,博士生導(dǎo)師,博士,研究方向?yàn)槎嗫捉橘|(zhì)理論、流固耦合.E-mail:xyang@shu.edu.cn