分区样条等参元方法分析加肋轴对称组合壳Ξ.PDFVIP

第15卷第1期 计 算 力 学 学 报 V ol. 15 N o. 1 1998年2月 CH IN ESE JOU RNAL O F COM PU TA T IONAL M ECHAN ICS February 1998 分区样条等参元方法分析加肋轴对称组合壳 黄加强　　郭日修 (海军工程学院一系, 武汉, 430033) 摘　要　本文提出分区样条等参元方法分析加肋轴对称组合壳在静水外压作用下的应力及稳定性。 该方法根据母线几何形状、壳板厚度变化及边界约束状况将加肋轴对称组合壳划分为若干个区域, 在不同区域采用不同步长的样条等参元进行分析, 因而便于对复杂的组合壳进行分析, 并可大大减 少计算所需自由度数, 减少计算工作量, 同时计算精度很高, 并且易于进行边界约束处理。应用该方 法对加肋锥柱、加肋锥环柱组合壳进行了计算, 其结果与其它方法计算值及实验值吻合良好。 关键词　样条等参元; 加肋轴对称组合壳; 应力分析; 稳定性 分类号　 241 O 1　引　　言 在工程结构中, 常采用加肋轴对称组合壳, 它由加肋的锥壳、柱壳、环壳、球壳等组合而成。 这种加肋轴对称组合壳在静水外压作用下的强度和稳定性问题受到力学工作者的高度关注。 由于加肋轴对称组合壳结构形式较为复杂, 因此其强度和稳定性用解析方法求解非常困 难[ 1—3 ] ; 而采用一般有限元方法求解, 则很难保证单元间几何形状及应力的连续性[ 4—7 ] , 计算精 度受到影响。 本文提出分区样条等参元方法分析加肋轴、对称组合壳强度和稳定性。该方法根据壳体子 午线几何形状、壳体厚度变化及边界约束状况等将加肋轴对称组合壳划分为若干个区域, 不同 区域采用不同步长的样条等参元进行分析。由于单元几何形状及位移函数均采用三次B 样条 函数逼近, 单元几何形状及位移不仅涉及到该单元结点坐标及结点位移, 而且还涉及到相邻单 元结点坐标及结点位移, 因而能很好地处理相邻单元在连接点应力和位移的连续性; 同时, 由 于在区域内采用一种新型的B 样条函数对边界局部样条函数进行修正, 因此相邻区域刚度 阵、载荷阵可以直接组装, 并且很容易进行边界约束处理; 另 外, 结构计算所需自由度数也较一般有限元法大为减少, 从而 大大减少计算工作量。 2　计算公式 2 1　壳体单元 在某区域上, 给定一个等距结点划分: : s s s …… s s s ( 1) 1 2 3 N - 1 N N + 1 　　结点间距为 h = si+ 1 - si 图1　区域结点划分及其延拓 收稿日期: 　黄加强: 男, 1968年生, 博士, 工程师 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 1期　　　　　　　　黄加强等: 分区样条等参元方法分析加肋轴对称组合壳　　　　　　 59 然后, 将该划分延拓成: ～ : s s s s …… s s s s (2) 0 1 2 3