弹性力学切应力大作业

1、 弹性力学1题目请分别用极坐标系下的半逆解法和有限元方法解答图中平面薄板在周边剪切荷载作用下的应力集中问题；并比较应力分布及应力集中系数的解析解和数值解2弹性力学解答【解】求出两个主应力，即原来的问题变为矩形薄板在左右两边受均布拉力q而在上下两边受均布压力q,如图所示:应力分量Ox=g6=q在孔边,边界条件是T=0代入坐标变换式，得到边界上的边界条件（6）p=r=qcos2qp（rw）P-r=_gsin2cpa）b）（c）d）（CFp）pr=0（CT/xp）p=r=0由边界条件式（a）、（b）、（c）、（d）可见，用半逆解法时，可假设为P的某一函数乘以心2q：而匸:为P的另一函数乘以而1d1护

2、pdpp2dp?d（1。环甲dpIpap因此可假设=/(p)cos2qpe)将(e)式代入相容方程得dY(p)|2d7(p)9d7(p)+9L(p)cospdp3p2dp2p3dp=0删去因子cos2(P以后，求解这个常微分方程，得/(p)=Ap4+Bp3+C+-其中A,B,C,D为待定常数，代入式(e)，得应力函数中WBp()()由应力函数得应力分量的表达式z4C6D、0P=-cos2qp(2B+),6=cos2cp(12Ap2+2B+卑)，32C6Dtp(p=sin2cp(6Ap+2B-P4)将上述式代入应力边界条件由式(a)得由式(b)得4C6D2B+f+f=:64R2+2B-2C6P(

3、g)由式(c)得由式(d)得R2用4C6D八2B+=Q,rr“2c2C6D6AR2+2B-r(h)(i)(j)联立求解式(g)(j)，并命二卩()，得C=qr2将各式系数值代入应力分量的表达式得r2L26=qCOS2cp(l-7)(1_3)o(p=-qcos2q)(1+3-y),P,屛ip=Tpp=-qsin2q)(l)(1+3)。沿着孔边P二厂，环向正应力是6=-4gcos2(p它的几个重要数值如下表所示：表1不同角丿度处的应力值030456090Otpq込0Tq沿着y轴，：皿，环向正应力是5)=g(l+3),它的几个重要数值如下表所示。表2距孔处不同距径处的应力值Pr1r3r斗F1.188

4、L037fL012(?1.005可以看出应力在孔边达到均匀拉力的4倍，但随着远离孔边而急剧趋近于q。3基于PDE求解结果有限元法是一种将连续体离散为有限大小单元体的集合,用以求解连续体力学问题的数值方法,有限元法具有求解精度高、通用性强等特点。3.1数值分析模型计算模型为长宽均为2mm,中间开一半径为0.2mm的圆孔的薄板，薄板的杨氏模量为1*103Mpa,泊松比为0.3。建立如下的几何模型，拉压应力均为1。2)划分网格3)求解结果X方向应力集中系数将理论解与数值解进行对比从上图可以看出，虽然理论解与数值解之间存在一定的误差，主要是由于PDE建模时的诸多因素如网格密度、约束条件简化等引起的，但

5、总的来说沿孔边的轴向正应力的理论解和PDE计算的数值解拟合的很好。3.2其它计算模型改变计算模型，薄板的尺寸不变，改变中心圆孔半径，分别为0.5mm和0.7mm,材料参数不变。1）孔径为0.5的计算模型2）孔径为0.7的计算模型2.理论解与数值解对比X方向应力集中系数从上图可以看出，在改变小孔半径之后，理论解与数值解存在明显差距，主要是因为孔口尺寸与薄板尺寸接近，弹性力学小孔应力计算方法已不适合该情况。同时对比pde计算结果可以发现，小孔的尺寸的大小显著影响孔口应力集中系数，直径越大，应力集中越明显。4误差分析产生误差的原因有：1.计算原理产生的误差有限元计算本身就是一种简化近似计算，在计算过

6、程中由于一些条件的简化，以及一些假设，使得计算结果与理论值有差别。模型的简化产生的误差本例实际是模拟一种小孔处应力集中现象，而在建模时，为了能够清楚看到洞口边应力应变变化，所以小孔选取较大，而整个模型本该延伸无限远处，所以这些造成模型简化误差。网格划分误差有限元软件计算的重要步骤就是划分网格，一般而言网格划分越密结果精度越高，但实际中考虑到经济问题，在一些我们不需要过多信息的地方我们选取粗网格，这也在一定程度上产生误差。结果选取误差由于关键部位网格划分较密，在选取结果时，只能选取在较近处的值，在比较积分点处的值由于网格差异积分点也会有一定的出入。6结论本次作业通过对小孔处应力集中现象的弹性力学有限元不同方法的模拟计算能得出如下结论：1.弹性力学通过几何方程，物理方程，力学平衡方程和边界条件能够很好地求解问题的理论解。弹塑性力学中,将带圆孔平板受双向对称拉压力的问题,近似为平面应力问题。通过弹性力学分析得到,孔边切向应力Oe在9=90。时最大，为4q。当0=90。时，Oe随