资源描述:
《应用力学学报-弹性梁式 薄板在横向绕流中的大 变形.pdf》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、第26卷第2期应用力学学报Vol.26No.22009年6月CHINESEJOURNALOFAPPLIEDMECHANICSJun.2009文章编号:10004939(2009)02030404*弹性梁式薄板在横向绕流中的大变形郝亚娟白象忠(燕山大学066004秦皇岛)摘要:采用相容拉格朗日欧拉(ULE)法,给出了弹性薄板理想流体横向绕流条件下变形与应力的理论算法,其中:对固体采用拉格朗日法;对流体采用欧拉法;对相互接触面采用拉格朗日法和欧拉法。建立了不间断横向绕流条件下弹性梁式薄板的大弯曲变形的非线性微分方程。求解
2、该方程时,将纵向位移分量和曲率的改变量用挠度表示。通过具体算例分析了各参数对悬臂梁式薄板挠度、纵向位移及应力大小的影响。理论解与数值解进行比较,验证了理论解的可靠性。关键词:流固耦合;相容拉格朗日欧拉(ULE)法;梁式板;挠度;纵向位移;应力中图分类号:O33;TB12文献标识码:A结合这两种方法。其优点在于:在求解流固耦合问1引言题时,可以直接利用流体力学和固体力学中的基本方程。在接触面上应用这两种描述方法,建立接触作为流固耦合问题的研究,多见结构在流场中面的运动学条件和动力学条件,从初始状态下表面的振动和稳定性问题的
3、报道,如文献[13]研究了圆接触条件的导出到运动状态下接触条件的确定,通柱振动问题,文献[46]研究了流体中板的振动问过动力学方程、边界条件就可以使这两种方法相互题。此外,文献[7]考虑横向剪切变形及流固耦合效联系起来。采用此法建立有关的流固耦合方程和关应时,建立了迭层板的非线性动力屈曲控制方程。系式,可方便地求解壳体的内力和变形及流场流动文献[8]介绍了板垂向位移与流体相互耦合的频域状态的变化。分析方法等等。而处于流场中结构的变形问题的报本文应用ULE法,求解了弹性悬臂梁式薄板在道却很少见。对于求解流固耦合问题所遇到的最大理想流
4、体横向绕流条件下的变形和应力,并讨论了流困难,就在于采用统一坐标系及两相界面的协调问速、板的几何尺寸、材料性质对变形和应力的影响。题。两相界面的位形事先未知,开始重合的点将随运动固体点的变形而发生位移。如何协调界面点的2弹性悬臂长薄板大变形问题非线性耦合问题,比起单纯固体或流体非线性问题的基本方程来,具有特殊的难度。由于非线性,叠加原理失效,[9]必须探讨全场求解的途径。如图1所示,厚度为h的弹性梁式薄板(0x为了解决非线性流体弹性力学问题理论解的困b,-5、欧拉(ULE)法,x=0边固定即研究弹性薄壳和流体相互作用时,壳体采用拉格w(x)u(x)=w(x)==0朗日法描述,流体采用欧拉法描述,在相互接触面处x*基金项目:河北省自然科学基金(A2006000190)来稿日期:20080104修回日期:20081007第一作者简介:郝亚娟,女,1977年生,燕山大学理学院,博士生;研究方向!!!耦合场理论。Email:moonhyj@sina.com.cn第2期郝亚娟,等:弹性梁式薄板在横向绕流中的大变形3052x=b边自由*11*11
6、Q1=D,N11=-D,***x2N11(x)=Q1(x)=M11(x)=023*1111Z3其中:u和w分别为薄板位移矢量在x轴和z轴的投2+=-(6)***x2D影;N11、Q1、M11分别为薄板变形后内力中的轴力、[10]曲率改变量的表达式为剪力、弯矩。222wuwuw求解在连续理想不可压缩的横向绕流条件下,11=-2+2-2(7)xxxxx弹性梁式薄板大变形时的挠度和应力状态。设:流体由于薄板的中面相对伸长量为零,有压力为p;流体为定常势流;无穷远处(z∀)流体u1u21w2∀11=+()+(
7、)=0(8)的压力、密度、速度分别为x2x2xup、、V。对于薄板的变式(8)舍掉后,得到xwu1w21w2形:认为中面不可延伸;=-()[1+()](9)hx2x4xwww[11]式(9)代入式(7),舍掉的高次项得到;#5。在上述指xbh2w1w211=-2[1+()](10)定状态下,流体运动的速度x2x图1悬臂长板绕流示意图势满足以下各式。x方向的位移u和曲率改变量11均用挠度w表方程组示,将式(4)和式(10)代入式(6),得22!=0*Dw2w2N11=-(2)[1+()]
8、,22(1)2xxp=p+[V-(!)]22*w1w2Q1=-D{2[1+()]},无限远处的条件xx2x4222w1ww21w2=Vz,(z∀)(2)4+2[2()]+(2)x2