梁单元问题

hnpency

各位老师、朋友，我想请教一个问题，我在用ANSYS做一个H单桩的分析，H型钢桩截面为100*100*6*8(mm)，桩总长为5m。桩在土里，我在桩头作用一个水平力，直至桩破坏。其中我将土等效为弹簧，这就变成一个二维的问题，因为我是做的破坏性分析，所以弹性分析的梁单元我不选用，比较过之后，我在犹豫是取beam188还是取beam23单元？我看过两者之间的比较，beam188它有现成的H截面供选择，但是它首先是个三维的，其次它是铁木辛柯梁，即考虑了剪切变形，它的计算精度在很大程度上将取决于单元划分的疏密程度，也就是单元划分的越细，它的精度越高；beam23是个二维的，但是它没有一个现成的H型截面供选择，它的实常数需要计算，假设是绕弱轴，参照ANSYS帮助文件，可以列出如下计算实常数的方程组：

A=0.0625*(A(-50)+A(50))+0.28935185*(A(-30)+A(30))+0.2962963*A(0)

I2=h^2*[0.015625*(A(-50)+A(50))+0.0260417*(A(-30)+A(30))]

I4=h^4*[0.0039063*(A(-50)+A(50))+0.00234375*(A(-30)+A(30))]
由于是对称截面，所以A(－50)＝A(50)，A(30)＝A(－30)
其中I2=Iy=133.48cm^4 , I4如何计算？
所以beam23的实常数我挺头疼的。


请各位老师、朋友指教，小弟不甚感激！