import sympy as sp
# 定义变量
b, h, l, t = 50, 100, 15, 2
# 计算各部分的面积
A_flange = b * t
A_web = (h - 2 * t) * t
A_lip = l * t
# 计算形心位置
y_flange = t / 2
y_web = (h - 2 * t) / 2 + t
y_lip = h - t / 2
A_total = 2 * A_flange + A_web + 2 * A_lip
# 形心相对于底部的高度
y_centroid = (2 * A_flange * y_flange + A_web * y_web + 2 * A_lip * y_lip) / A_total
print(f'形心相对于底部的位置&