索膜结构是一种新型建筑结构形式, 近年来发展相当迅速, 其材料与制作, 试验与检验方法以及施工方法等都与普通建筑结构有所不同[1]。索膜结构裁剪分析是索膜结构设计过程中的一个关键问题, 如果裁剪分析不当, 将极大地改变索膜结构中原来的应力分布, 甚至出现褶皱。通常索膜结构是具有双向曲率的曲面, 多数是不可展的。最好的近似方法当然是用许多具有无限小宽度的膜条来形成整个结构, 但在实际工程中是不可能做到的[2]。用测地线来进行裁剪分析, 可使膜片的边界不产生过大的弓形。本文根据文献[3][4]提出的方法, 将其进行展开和深入, 推导了最终求解测地线的公式, 解决了索膜结构裁剪设计中的一个关键问题。

在求解测地线的过程中, 百富策略白菜网单元—节点拓扑矩阵可以大大的简化和方便计算, 因此有必要先引入单元—节点拓扑矩阵的概念。

设网络划分, 单元数为m, 节点数为n, 其中元素定义如下:对于i单元, 始端节点为r, 终端节点为s, 单元—节点拓扑矩阵C和D为矩阵m×n, 可描述一个网络关系。设第j单元, 始端节点为i, 终端节点为i+1, 单元坐标差为:

uj=xi+1-xi;vj=yi+1-yi;wj=zi+1-zi;

qj=xiti+1-xi+1ti;rj=yiti+1-yi+1ti;

sj=ziti+1-zi+1ti。

设为矩阵形式:

u=C·x;v=C·y;w=C·z。

q=D·x;r=D·y;s=D·z。

其中, u, v, w, q, r, s均为单元坐标向量 (m×1) ;x, y, z均为节点坐标向量 (n×1) 。

已知曲面上的两点, 求经过这两点的测地线, 该问题可以表述为:给定曲面φ (x, y, z) =0, 求曲面上所有连接已知点A (x1, y1, z1) , B (x2, y2, z2) 的曲线中, 长度最短的曲线。问题归结为求泛函的极小值[2,3]:

l=∫x2x11+y′2+z′2−−−−−−−−−−√dxl=∫x1x21+y′2+z′2dx (1)

其中, y (x) , z (x) 满足φ (x, y, z) =0, 求测地线问题为条件极值问题, 用拉格朗日乘子法可化为求解下述辅助函数的极值问题, 并取参数方程。

x=x (t) ;y=y (t) ;z=z (t) ;λ=λ (t) 得到式 (2) :

l=∫T0[x2t+y2t+z2t−−−−−−−−−−√+λ(t)ϕ]dtl=∫0Τ[xt2+yt2+zt2+λ(t)ϕ]dt (2)

其中, xt=dx/dt;yt=dy/dt;zt=dz/dt, 当x=x1时, t=0;x=x2时, t=T。

若已知曲面方程φ (x, y, z) =0, 则可由上式求得极值函数x=x (t) , y=y (t) , z=z (t) 和λ=λ (t) 。由于我们得到的初始曲面是一些离散点, 得不到初始曲面具体方程, 因此得不到极值函数x=x (t) , y=y (t) , z=z (t) 和λ=λ (t) 。然而可以用分段近似的方法确定测地线。

因为无法得到初始曲面的具体方程, 可以利用分段近似的方法确定测地线, 其具体流程如下[4]:

1) 已知膜曲面的离散点p1, p2, p3…点的坐标, 把这些点组成一个集合K;2) 已知测地线的端点A (x1, y1, z1) , B (x2, y2, z2) , A, B点在边界索上;3) 连接A, B两点, 计算线段AB的长度, 根据AB的长度, 将线段n等分, 得到n个单元, (n-1) 个内插点M1′, M2′, M′n-1, 把这些点组成一个集合为G, G内各点的坐标可通过计算求得, 为已知量;4) 取出第一个单元AM1′, 以其中点为球心, R=10×|AM1′|为半径作球, 球内所包含的K集合

two parallel rigid cylinders or between two rigid spheres[J].Rubber Chem.Technol., 1986, 59 (1) :77-83.

[4] O.H.Yeoh.Relation between crack surface displacements and strain energy release rate in thin rubber sheets[J].Mechanics of Materials, 2002 (4) :315.

[5] A.N.Gent, Y.C.Hwang.Elastic behavior of rubber layer bonded between two rigid spheres[J].Rubber Chem.Technol, 1988 (61) :630.

On nonlinear finite element analysis of radial stiffness of pipe rubber bearing

GAO Jian-hong CHEN Jia-bin

Abstract:

Aiming at the big deformation characteristics of rubber, the paper has the nonlinear finite element analysis of the common pipe rubber bearing in projects, concludes the relationship between the stiffness of the radial loading from the bearing and the sizes of the changes and the bearing and obtains the similar calculation model for the radial stiffness by using engineering regression tools.

Key words:

rubber bearing, radial, nonlinear, finite element

内的离散点为P1′, P2′, P3′…, 这些点组成一个集合为P′, P′⊆K;5) 利用P′中的离散点, 进行曲面拟合, 设近似曲面方程为:z=h0+h1x+h2y+h3xy+h4x2+h5y2, 用最小二乘法得到曲面方程的系数h0, h1, h5…, 设ϕ1′ (x, y, z) =z-z (x, y) ;6) 将M1′的x, y坐标代入ϕ1′ (x, y, z) =z-z (x, y) 得到膜曲面上点M1 (x, y, z) , 即为测地线上的一个初始点;7) 对第2个～第n个单元重复步骤4) ～步骤6) , 可得到测地线上的一系列初始点A, M1, M2, …, Mn-1, B, 如图1所示;8) 取第i个单元MiMi+1, 用4) , 5) 同样的方法, 以MiMi+1的中点为球心, R=5×|MiMi+1|为直径作圆, 所包含的离散点进行拟合得到的单元MiMi+1附近的曲面方程ϕj (x, y, z) =z-z (x, y) ;9) 对于第i个测地线单元, 单元内任一点的坐标为:

⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪x=xi+1−xili0+xiti+1−xi+1tili0=ajxt+bjxy=yi+1−yili0+yiti+1−yi+1tili0=ajyt+bjyz=zi+1−zili0+ziti+1−zi+1tili0=ajzt+bjzlj0=ti+1−ti{x=xi+1-xil0i+xiti+1-xi+1til0i=axjt+bxjy=yi+1-yil0i+yiti+1-yi+1til0i=ayjt+byjz=zi+1-zil0i+ziti+1-zi+1til0i=azjt+bzjl0j=ti+1-ti

(3)

对每一个单元设λ为常数, 不随t变化。

根据以上9个步骤得到的每一单元的曲面方程和单元内任一点坐标的参数表达式, 进而就能推导求解测地线的具体公式。因此, 式 (1) 即测地线的总长可以表示为:

l=∑(lj1+λjlj2)l=∑(l1j+λjl2j) (4)

把经过拟合所得到的曲面方程ϕj (x, y, z) =hj00j+hj11jx+hj22jy+hj33jxy+hj44jx2+hj55jy2-z (j=1, 2, 3, …, n-1) 和式 (3) 代入式 (4) , 通过化简可以得到:

lj1=∫ti+1tix2t+y2t+z2t−−−−−−−−−−√dt=u2j+v2j+w2j−−−−−−−−−−√l1j=∫titi+1xt2+yt2+zt2dt=uj2+vj2+wj2。

lj2=hj0lj0+hj1(Bjuj+qj)+hj2(Bjvj+rj)+hj3lj0(Ajujvj+Bjujrj+Bjvjqj+qjrj)+hj4lj0(Aju2j+2Bjujqj+q2j)+hj5lj0(Ajv2j+2Bjvjrj+r2j)−Bjwj−sjl2j=h0jl0j+h1j(Bjuj+qj)+h2j(Bjvj+rj)+h3jl0j(Ajujvj+Bjujrj+Bjvjqj+qjrj)+h4jl0j(Ajuj2+2Bjujqj+qj2)+h5jl0j(Ajvj2+2Bjvjrj+rj2)-Bjwj-sj。

其中, Aj=13(t2i+1−ti+1ti+t2i)Aj=13(ti+12-ti+1ti+ti2);Bj=12(ti+1+ti)Bj=12(ti+1+ti)。

把式 (4) 写成向量形式:

l=kl1+λTl2。

或表示为:

l=kl1+lT22Τλ (5)

其中, k=(1‚1‚⋯‚1)k=(1‚1‚⋯‚1)n个1 ;l1, l2均为n×1的向量;λ为n×n阶对角矩阵。

对目标函数式 (5) 求变分得[5]:

⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪∂l∂xδx=k∂l1∂xδx+λT∂l2∂xδx∂l∂yδy=k∂l1∂yδy+λT∂l2∂yδy∂l∂zδz=k∂l1∂zδz+λT∂l2∂yδz∂l∂λδλ=lT2δλ{∂l∂xδx=k∂l1∂xδx+λΤ∂l2∂xδx∂l∂yδy=k∂l1∂yδy+λΤ∂l2∂yδy∂l∂zδz=k∂l1∂zδz+λΤ∂l2∂yδz∂l∂λδλ=l2Τδλ

(6)

令上式各泛函等于0, 即:

⎧⎩⎨⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪⎪∂l∂x=k∂l1∂x+λT∂l2∂x=0∂l∂y=k∂l1∂y+λT∂l2∂y=0∂l∂z=k∂l1∂z+λT∂l2∂y=0∂l∂λ=lT2=0{∂l∂x=k∂l1∂x+λΤ∂l2∂x=0∂l∂y=k∂l1∂y+λΤ∂l2∂y=0∂l∂z=k∂l1∂z+λΤ∂l2∂y=0∂l∂λ=l2Τ=0

(7)

上述方程组为一非线性方程组, 采用牛顿—拉弗逊法和最速下降法联合求解[2], 可得到x, y, z, 即测地线AB上各个分段点M1, M2, …, Mn-1的坐标。

得到测地线AB的各分段点后, 用线段将其连接起来, 作为裁剪设计中的裁剪线。

本文百富策略白菜网有限单元的思想, 采用古典变分求极值的方法, 给出了求测地线的具体步骤和详细流程, 推导了确定测地线的相关公式。此方法可以任意改变测地线两端点的位置以及测地线的方向, 充分利用了材料。并且根据本文的方法, 作者编制了相应的计算机程序, 通过多个算例的计算, 把所得的结果和EASY软件计算出的结果作比较, 发现二者计算结果的差异均在千分之一以内, 具有相当高的精度。采用测地线作为裁剪线, 子曲面展开图形边界曲率较小, 裁剪片幅度大致相等, 从而达到拼接简单和省料的目的。