提亚峰 季宝锋 黄 诚 王世勋 李雄魁
(北京宇航系统工程研究所,北京 100076)
纤维缠绕是制造格栅(geodesic)增强复合材料结构(正交格栅和等边三角形格栅结构壁板)的理想方法,它可使结构的质量很轻,而强度刚度却非常高[1]。
复合材料网格结构在国内应用于工程实践时间尚短,但自动化缠绕成型趋势已十分明显。自动化缠绕成型可提高产品质量、提高生产效率,便于批量生产。
圆柱段网格加筋结构,除纵筋沿母线、环筋沿圆环外,其余筋的中心线均为圆柱螺旋线。截圆锥段网格加筋结构,除纵筋沿母线、环筋为圆环外,其余筋的中心线均为圆锥螺旋线。蒙皮的缠绕除90°沿圆环外,其余缠绕角度均沿圆柱圆锥侧面上的短程线缠绕。因此应对圆柱圆锥侧面上的螺旋线和短程线的几何性质进行深入系统的研究,得出工程上可行的方法以应用于圆柱段和截圆锥段网格加筋结构的几何建模、结构承力计算、自动缠绕时缠绕编程,为复合材料网格结构自动化制造创造条件。本文主要介绍了锥柱面螺旋线短程线与网格缠绕成型技术。
圆柱段复合材料网格结构,除纵筋和环筋外,其余筋的中心线均为圆柱螺旋线,其轨迹方程为:
(1)
式中,α为螺旋角。因螺旋角α为定值,当θ=2π为时,z为定值,即圆柱螺旋线是等螺旋角等螺距螺旋线。
圆柱螺旋线方程(1)的切线方程为:
(2)
设圆柱螺旋线的切线与z轴的夹角为ω,则:
(3)
由于圆柱的母线与z轴平行,圆柱螺旋线的切线与z轴的夹角就是圆柱螺旋线的切线与圆柱母线的夹角。圆柱螺旋线的切线与母线的夹角ω与其螺旋角α相等且为定值。
截圆锥段网格加筋结构,除纵筋和环筋外,其余筋的中心线均为圆锥螺旋线。
将螺旋角设为定值,得到等螺旋角圆锥螺旋线轨迹方程[2]:
(4)
式中,β为半锥角。ρ0为坐标原点到截圆锥体小端端面的圆锥侧面上的母线长度。因方程中含有指数函数,等螺旋角的圆锥螺旋筋槽对加工设备要求较高。
将螺距设为定值,得到等螺距圆锥螺旋线轨迹方程[3]:
(5)
螺距h=2πacosβ
曲线切线与锥面母线之间夹角的余弦:
式中,aθ为坐标原点到截面锥体小端端面的圆锥侧面上的母线长度。等螺距圆锥螺旋筋槽对加工设备要求较低。模具等速旋转,刀具等速沿圆锥侧面母线运动,即可加工出等螺距圆锥螺旋线。
等螺旋角和等螺距螺旋线在圆锥侧面上是两种曲线,在圆柱侧面上是同一种圆柱螺旋线。
为便于圆柱段和圆锥段复合材料网格结构硬模成型装置的设计生产和纤维缠绕,需要求出圆柱和锥段螺旋线的弧长和曲率。
圆柱螺旋线的弧长为:
(6)
等螺旋角圆锥螺旋线的弧长为:
(7)
(8)
等螺距圆锥螺旋线的弧长为:
(9)
(10)
arsh(θ2sinβ)-arsh(θ1sinβ))
(11)
圆柱螺旋线:
(12)
曲率[2]:
(13)
即圆柱螺旋线曲率仅与螺旋角和圆柱半径有关,α=0°时为母线,曲率为0,对应复合材料网格结构纵筋;α=90°时曲率最大,对应环筋。
等螺旋角圆锥螺旋线
(14)
(15)
等螺距圆锥螺旋线
(16)
其曲率为:
(17)
筋条的缠绕,由于模具上加工了筋槽,纤维丝束可在存在张力的情况下,缠绕到筋槽内,从而成型各种螺旋角的螺旋筋条。在表面上缠绕纤维最稳妥的方法是沿着两点之间的最短路线进行缠绕[1],即蒙皮的缠绕,纤维丝束在张力作用下应沿短程线缠绕。
圆柱段蒙皮,除0°和90°外,纤维沿圆柱侧面上的短程线进行缠绕,如图1所示。
图1 圆柱面上的短程线
设圆柱面方程为:x2+z2=1,(x>0)。并设在圆柱面上连接A(x0,y0,z0),B(x1,y1,z1)两点的曲线方程为:
(18)
(19)
其欧拉方程组:
(20)
由于H不显含y,欧拉方程组化为:
(21)
求其首次积分,并用弧长s作参数,得:
(22)
欧拉方程组化为:
(23)
较之青壮年而言,中老年群体生理、心理均处于特殊阶段,情绪波动容易受外界刺激、劳累、疾病等因素影响,因此脑卒中后极易发生抑郁症状。调查显示,抑郁不仅会导致患者产生“焦虑、悲观”等消极情绪,加强身体负性应激及影响治疗配合度,另外还会进一步加重疾病症状及异常体征,致使患者出现更为显著的神经及心理功能异常[3]。由此可见,上述两者可相互促进、互为因果,因此临床治疗中单纯使用药物治疗很难消除负性情绪对疾病治疗的影响,因此临床多辅助针对性护理工作。
(24)
根据ds的定义有:
(ds)2=(dx)2+(dy)2+(dz)2
(25)
(26)
(27)
(28)
于是:
(29)
再对欧拉方程组第一式积分得:
y=C1s+C3
(30)
综合式(28)、(29)、(30):
(31)
由式(1),圆柱侧面上的短程线是其上的螺旋线。
由式(22),用弧长s作参数,弧长s与图1所示xOz面的极坐标θ及r间的关系为s=cscα·r·θ,为推导方便,已设x2+z2=1,(x>0),即r=1,因此s=cscα·θ,y=cotα·θ。由式(23)第一式,积分常数c1=cosα。由式(28),特定的取积分常数c2为0,事实上当圆柱螺旋线的起点是z轴时,c2=0。
与式(1)相比,这是圆柱半径r=1,以xOz面为圆柱底面向y轴延伸的圆柱螺旋线,而式(1)是圆柱半径为r以xOy面为圆柱底面向z轴延伸的圆柱螺旋线。即式(1)是式(31)的起点在坐标轴的特例。
圆锥段蒙皮,除0°和90°外,纤维沿圆锥侧面上的短程线进行缠绕,如图2所示。设圆锥面方程为:x2+z2=y2tan2θ,(x>0),并设在圆锥面上连接A(x0,y0,z0),B(x1,y1,z1)两点的曲线方程为:
图2 圆锥面上的短程线
作辅助函数:
(32)
(33)
其欧拉方程组:
(34)
其欧拉方程组化为:
(35)
不易分离变量,难以求出短程线方程的解析解。
设有直角三角形纸片的一个锐角为α,将此纸片卷在一圆柱侧面上,使角α的直角边与圆柱的母线重合,角α的斜边则在圆柱面上盘旋上升,这样在圆柱面上就形成一条空间曲线,这条空间曲线就是圆柱螺旋线。上述过程的逆过程,就是将圆柱螺旋线展开,如图3所示。圆柱面展开其上的圆柱螺旋线是一条直线。能展开的平面,可利用两点之间线段最短来求解。还可以此验证等螺旋角及等螺距圆锥螺旋线是否是圆锥面上的短程线。
图3 圆柱螺旋线的展开
另外三角网格的称谓对于圆柱结构是指展开后环筋、纵筋和螺旋筋构成三角形,实际结构三个节点间不会象平面三角形一样稳定。筋条中心线与受力方向夹角越大、曲率越大,越不利于承受对应载荷,即纵筋轴压承载效率高,环筋外压承载效率高,单独轴向直筋承载效率高于螺旋筋。
圆锥侧面展开图如图4所示。
图4 圆锥侧面展开图
假设曲线的极坐标方程为:
r=r(Θ)
(36)
利用直角坐标与极坐标的关系;x=rcosΘ,
y=rsinΘ则方程组:x=r(Θ)·cosΘ及y=r(Θ)·sinΘ是曲线的参数方程,其中参数为极角Θ,曲线的切线斜率是:
(37)
根据几何关系可得出公式:
α=γ-Θ
(38)
(39)
将y′的表达式(37)代入(39):
(40)
(41)
取y=0为初始点(图中B点),此时Θ=0,代入式(41)则:
c=-lnx0
(42)
将式(42)代入式(41)得:
(43)
因此该曲线极坐标方程为:r(Θ)=x0ecotα·Θ,是形如r=aekθ的对数螺旋线。
等螺距圆锥螺旋线在展开平面上的曲线方程为:r=x0Θ[2],是形如r=aθ的阿基米德螺线。
阿基米德螺线上A1,A2,A3,A4间的距离相等,如图5所示。
图5 阿基米德螺线
等螺旋角及等螺距圆锥螺旋线不是圆锥面上的短程线,因为圆锥侧面展开之后这两条曲线分别是对数螺线和阿基米德螺线,不是直线。
由于圆锥侧面可展开为扇形,等螺旋角及等螺距圆锥螺旋线在圆锥侧面展开之后分别是对数螺线和阿基米德螺线,对于脱模困难而网格尺寸精度要求不高的双锥等结构,可将锥段阳模展开成平面,平面浇注后,再贴合到芯模上形成锥面螺旋筋槽,应注意阿基米德螺线的加工易于对数螺线。
为便于锥段复合材料网格结构阳模成型装置的设计生产和纤维缠绕,需要求出对数螺线和阿基米德螺线的弧长和曲率。
对数螺线r(Θ)=x0ecotα·Θ(-∞<Θ<+∞)从矢径Θ=Θ1到矢径Θ=Θ2的弧长。
(44)
阿基米德螺线r=x0Θ(-∞<Θ<+∞)从矢径Θ=Θ1到矢径Θ=Θ2的弧长。
(45)
对比展开前和展开后可知,Θ=sinβ·θ,
(46)
对数螺线的曲率:
(47)
(48)
(49)
曲率[4]:
(50)
(51)
阿基米德螺线的曲率:
(52)
(53)
(54)
(55)
(56)
利用两点之间线段最短来求解圆锥侧面上的短程线,如图6所示。设圆锥高为h,圆锥面方程:
x2+y2=(z-h)2tan2β,r=htanβ,A(r,0,0),B(x0,0,z0),0 图6 圆锥面短程线的展开 (57) (58) (59) 从而在xOz上线段AB的方程为: (60) 在线段AB上任取一点N′(x′,z′),则当展开面重又绕回到锥面上时,其对应点N的坐标设为(x,y,z),求N点坐标,如图7所示。 |MN|=|MN′|,∠AMN′=φ=θ·sinβ [与前述(46)式Θ=sinβ·θ结论相同] 图7 圆锥面短程线 (61) -cos(φ+β)) (62) (63) ∵N′在AB上, (64) (65) (66) 再将式(58) 、(59)中xB,zB带入,即得到最终解。方程非常烦琐,不适合用于筋条的缠绕。但蒙皮的缠绕,纤维丝束在张力作用下是沿此方程缠绕的,即圆锥面蒙皮的缠绕其缠绕角度只能是近似按设计规定进行。在半锥角较小的情况下,此偏差影响不大[5]。缠绕时要注意对纤维束起点位置的调整。 另外,锥面短程螺旋线对于平面浇注软模成型的复合材料网格结构极具应用价值。在浇注阳模的模具上按图4所示加工与扇形径向线夹角为α的直线筋槽即可。这将大大降低加工难度,提高加工效率,降低成本。为此计算直线与对数螺线的曲率差,由公式(43)r(Θ)=x0ecotα·Θ,结合工程应用中的小直径锥段,取α=30°,x0=1244.7。根据公式(50),θ=0°时,k得最大值为0.0004017,即曲率很小,可用直线代替对数螺线,不影响结构承载。对于大直径锥段,曲率更小,更接近直线。 通过对圆柱、圆锥侧面上的螺旋线和短程线的探讨,解决了圆柱、圆锥段复合材料网格结构筋条及蒙皮的缠绕问题;得出了圆柱和锥段螺旋线的弧长和曲率计算方法,用于圆柱段和圆锥段复合材料网格结构成型装置的设计生产和纤维缠绕;得出了对数螺线和阿基米德螺线的弧长和曲率计算方法,用于锥段复合材料网格结构软模成型装置的设计生产和纤维缠绕。为复合材料网格结构自动化缠绕生产制造创造了条件。同时针对工程需要,对不同的网格结构选择不同的加工成型方式,可大大降低加工难度,实现低成本快速批量生产。 [1] 牛春匀.实用飞机复合材料结构设计与制造[M].北京:航空工业出版社,2010:176. [2] 王连祥,方德植.等.数学手册[M].北京:人民教育出版社,1979:409,411-412. [3] 提亚峰,张铎,董波.C/E复合材料网格缠绕结构一体化设计[J].宇航材料工艺,2010,40(6):33-37. [4] 同济大学数学教研室主编.高等数学(第三版,上册)[M].北京:高等教育出版社,1988:217. [5] 提亚峰.复合材料网格结构参数调整及优化计算[J].宇航系统工程,2013,3(3):26-34.3 结论