薄壳结构几何非线性全过程分析
GSQ24壳体单元用于板壳结构的几何非线性分析
收稿日期：2001-04-21；修改日期：2002-05-14
元能顺利通过弯曲、膜小片检验，无任何剪切和膜 闭锁及多余零能模式，并且适用于厚、薄板壳结构 的弹塑性分析和组合结构分析。 对于较薄的板壳结构，小挠度理论夸大了挠度 和应力，通常按小挠度理论进行计算和设计的结 果，是偏于安全和浪费的，对实际较薄板壳结构宜 采用大挠度分析。由于壳体单元局部坐标到整体坐 标需转换，本文基于 Von Kaimam 大挠度、中转角 下采用全 Lagrange 描述方法分析板壳结构的大挠 度问题。
1
引言
在壳体结构有限元分析中，尤以平板型壳元最 为简单和直观，因而最早得到发展和广泛地应用。 由于平板型壳元将壳体离散为一些折板，这可避免 复杂的曲面几何描述，并能方便地满足刚体位移及 常量应变的位移模式，其单元刚度矩阵也非常简 洁，而其几何上的近似可通过改善平板弯曲及膜向 拉伸性能来补偿，因而目前人们的兴趣又重新转向 平板壳元[1]。但要指出目前非常理想的平板单元较 少，所遇到的主要困难是：单元共面造成刚度矩阵 奇异、膜向拉伸与平板弯曲两部分的位移插值阶次 如何匹配及相邻单元边界位移不协调等问题。 GSQ24 壳体单元是一种带有平面内旋转自由度的 四结点平板单元，每个结点六个自由度，此种单元 由基于 Mindlin 板理论的平板弯曲单元和基于修正 后的 Reissner 泛函中引入独立转动场变分原理的平 面内膜单元组合而成。 GSQ24 壳体单元一方面避免 了常规壳体单元中每个结点五个自由度在单元共 面时造成刚度矩阵奇异问题；另一方面通过在膜元 位移模式中引入适用于任何四边形单元的不协调 位移函数，使得膜元和弯曲元中位移插值函数阶数 匹配问题也得到妥善解决。另外在应用上易于模拟 复杂的壳面连接和易于与别的带有转动自由度单 元相容。经文献[2,3,4]中的应用表明，GSQ24 壳体单
几何非线性对单索面薄壁斜拉桥的影响分析
Abstract： Taking Wusu cable － stayed bridge practical engineering as an example，through finite element sim － ulation calculation， this paper obtains the calculation results of the internal force and deformation of different parts under the influence of geometric nonlinear factors. Key words： thin － wall cable － stayed bridge； geometric nonlinearity； mechanical property
分析计算指标选取乌苏斜拉桥不同组成部位处 效应值，分别 为 主 梁 纵 向 位 移、主 梁 弯 矩、主 梁 轴 力、 桥塔位移、桥塔弯矩和桥塔轴力值六个指标。分析阶 段则选择成桥后的运营阶段。 3. 2 计算结果
为验算乌苏斜拉桥在运营阶段的几何非线性影 响，计算加载形式采用最不利影响线加载方式。根据 规范要求，加载荷载采用均布车道荷载和集中荷载， 在影响线最不利同号区域进行加载计算。
1 工程概况
乌苏斜拉桥地处东北边陲黑龙江省抚远县黑瞎 子岛，跨越乌苏里江，与俄罗斯毗邻。该桥是《中俄区 域合作规划纲 要 》的 重 点 项 目，其 社 会 意 义 和 战 略 作 用巨大［1］。
乌苏斜拉桥采用独塔单索面钢 － 混结合梁形式， 跨径组合为 140 m + 140 m = 280 m。乌苏斜拉桥采用 塔、梁、墩固结 体 系，主 塔 为 钢 筋 混 凝 土 独 柱 式 结 构， 主梁应用两侧带长挑臂的钢箱结合梁结构形式，斜拉 索采用竖琴形布置方式，斜拉索长度为 20. 575 m 至 147. 651 m，全桥共有 52 根斜拉索。乌苏斜拉桥设计
第三章薄壳结构
第二节 结构分析
一．结构内力
以应力表示
以内力表示 3个薄 膜内 力 Nx , Ny , Nxy = Nyx
5个弯 曲内 力 Mx , M y , Qx , Qy , Mxy = M yx 共计 个内力 8
二．分析方法
1. 解析法
直接以数学方式得到基本方程的解。
2. 半解析法
对平衡方程或几何方程、 物理方程中省略 某些项。
2 2 2 2 k
Et D= 12 1−ν 2
3
(
)
薄膜 内力 ∂Φ Nx = 2 ∂y
2
弯曲 力 内 ∂2w ∂2w Mx = −D 2 +ν 2 ∂x ∂y ∂ w ∂ w My = −D 2 +ν 2 ∂y ∂x ∂2w Mxy = −D(1−ν ) ∂x∂y ∂ 2 Qx = −D ∇ w ∂x ∂ 2 Qy = −D ∇ w ∂y
重200kg/m2。 t R = 60 12500 ≈ 1 208
◆ 移动曲面
筒壳
锯齿形锥壳
山西平遥县棉织厂厂房扩建工程，建于1983年。
折板
◆ 组合曲面
圆柱面切割组合
组合扭面
美国TWA环球航空公司候机楼
美国著名建筑师 沙里宁1961年设计，用4片钢筋砼扁壳组成，形似一只 正要起飞的大鸟。
q
qϕ = q sin ϕ，qr = qconϕ q
rd ϕ 0
r sinϕ0
ϕ0
ϕ
r
Nϕ0
ϕ0 dϕ 0
ϕ0
r
Nϕ 0
当只考虑壳自重荷载时，由 ∑Z = 0 有
∫
ϕ
0
q2πr2 sin ϕ0dϕ0 = −2πr2 (1− cosϕ)q
薄壳结构 (1)ppt课件
1.理解为四边支承的曲板。 2.主要依靠曲面内的双向轴力和顺剪力承重。 3.强度和刚度主要依靠几何形状的合理性，而不是结构
截面尺寸得到。 4.空间整体工作性能良好，内力均匀，结构自重小； 5.强度高、刚度大、材料省、经济合理。 6.曲面多样化，丰富建筑造型。
精选课件PPT
11
特点
构，它的内力还是弯矩。
精选课件PPT
7
• 把平板做成曲板， 曲板的内力就改变 为受压为主，受压 比受弯更能发挥材 料的性能，尤其是 多向受压，处于空 间状态更加有利。
精选课件PPT
8
• 横向受荷传力的梁起“担” 的作用，不能材尽其用， 并非经济的结构形式；以 曲梁承荷传力的拱起“顶” 的作用，能进一步发挥材 力，是较先进的结构形式；
横隔板的型式
框架横隔 精选课件PPT 空腹桁架横隔
32
二 、筒壳的受力特点
是否有横隔板是筒壳和筒拱的区别 筒壳 横向 与拱类似 壳身产生环向压力
纵向 与梁类似 把上部竖向荷载传递给横隔
➢长壳 l1/l2 3
横向拱的作用小；纵向梁传力显著。近似梁作用 按梁理论计算 ➢短壳 l1/l21/2 横向拱的作用明显；纵向梁传力作用很小。近似拱作用 内力主要为薄膜内力，按薄膜理论计算
墨西哥霍奇米洛科餐厅
立面
平面
几何形体
建筑由八瓣鞍壳交叉组成，相交处加厚形成刚度极大 的拱肋，直接支承在八个基础上
建筑平面为30m*30m的正方形，壳厚40mm
壳体的外围八个立面是倾斜的，整个建筑犹如一朵覆 地莲花，造型别致室内采光、通风效果好
精选课件PPT
41
立面图 脚点处理
美国麻省理工学院礼堂
18
第六章薄壳
壳体的厚度比其它尺寸（如长度、曲率 半径等）小得多的壳体称为薄壳。 由壳体厚度中点构成的曲面称为中曲面。 当薄壳受载荷而发生微小的变形时，与 薄板相似，也可以忽略沿壳体厚度方向 的挤压变形和应力，且认为符合直法线 假设，即薄壳中曲面法线上各点在变形 过程中仍保持在变形后中曲面的法线上。
(6.1)
其中，平面变形的子刚阵，弯曲变形的子刚阵， 对角零元素为对应于θzi’的空白部分。其余空白部 分皆为零元素，这反映出，在局部坐标系内，平 面变形与弯曲变形是互不相关的。
单元刚阵的坐标变换与叠加
式(6.1)是在单元的局部坐标系x’y’z’内给出的 矩阵，而组成柱壳的各个平板单元显然不会都 在一个平面内，因而各单元的局部坐标一般是 不一致的。 这样，与同一节点相邻的单元就可能有不同的 局部坐标。 与第一章中刚架分析相似，这种不同坐标系给 出的单元刚阵是不能直接叠加的，必须变换到 统一的坐标系。
概述
工程结构中，为了减轻重量、美观实用，常常 采用一些壳体结构，如柱面及球面屋顶、各种 压力容器、航天及航空器的薄壁结构等。 还有不少机器的部件本身就是一个壳体，如水 轮机的叶片、各种外壳等。 壳体的几何形状及其变形都是很复杂的。为了 便于分析，针对各种壳体情况做了不同的假设， 建立有不同的壳体理论和相应的有限元分析方 法。
如单元沿法线z’方向受有单位面积的横向力pz’， 则分配到4个节点有弯曲的单元节点载荷
Q N p dx' dy'
' e b T e b ' z
其中[Nb]即四节点矩形平板弯曲单元形状函数 矩阵(5.11)式.
由此，单元e的任一节点i将分配有z’方向的集 中力Z i’和绕x’、y’轴的集中力偶M xi’、 M yi’ 。
第5章_薄壳结构
•支承在框架上
斜拱
•直接落地并支承在基础上
Central South University of Forestry & Technology
第5章 钢筋混凝土空间薄壁结构 5.2 圆顶 5.2.2 圆顶的受力特点
1）圆顶的破坏
2）圆顶的薄膜内力
壳面单元体的主要内力 经向应力状态 环向应力状态
圆顶的坐标及薄膜内力
主方向：对应于每一个主曲率的方向称为曲面在o 点的主方向，两个主方向 是相互正交的。
曲率线：在曲面的每点处在主方向上的两根切线，这些切线围成了两组互成 正交的曲线网，称为曲率线
Central South University of Forestry & Technology
第5章 钢筋混凝土空间薄壁结构 5.1 概述 5.1.1 薄壳结构的概念 按高斯曲率分类
第5章 钢筋混凝土空间薄壁结构 5.1 概述 5.1.2 薄壳结构的曲面形式
3、直纹曲面
•锥面: 由一段直线（母线）沿一竖向曲线（导线）移动并始终通过 一定点形成的曲面
•劈锥曲面 由一段直线（母线）沿一直导线和一根竖向曲导线移动并始 终平行于一导平面形成的曲面
劈锥曲面
锥面
Central South University of Forestry & Technology
椭圆抛物面
双曲抛物面
Central South University of Forestry & Technology
第5章 钢筋混凝土空间薄壁结构 5.1 概述 5.1.2 薄壳结构的曲面形式
2、平移曲面
椭圆抛物面
双曲抛物面
Central South University of Forestry & Technology
第5章_薄壳结构
第5章 钢筋混凝土空间薄壁结构 5.1概述
梁式结构 拚架结构 刚架结构 拱式结构
面外:需设支撑体系保证安 全及稳定 平面受力结构
面内:承受屋面板传来的竖 向荷载
平面受力结构体系的特点: 优点:荷载为单向传递，计算分析方便，结构施工吊装方便。 缺点:结构内力较大，材料强度得不到充分发挥，随着结构材料用量
第5章 钢筋混凝土空间薄壁结构
System and Selection of Building Structures 土木工程学院
第5章 钢筋混凝土空间薄壁结构 5.1 概述 5.1.1 薄壳结构的概念 思考:一个鸡蛋能受多大的压力？ F1=?
壳体结构的强度和刚度主要是利用其几何形状 的合理性，而不是以增大其结构截面尺寸取得的， 这是薄壳结构与拱式结构相似之处。
现浇混凝土壳体 预制单元、高空装配成整体壳体 地面现浇壳体或预制单元装配后整体提升
装配整体叠合壳 采用柔模喷涂成壳 预应力混凝土结构
第5章 钢筋混凝土空间薄壁结构 5.1 概述 5.1.4 薄壳结构的施工方法
蜂巢芯 薄壁箱体
Central South University of Forestry & Technology
第5章 钢筋混凝土空间薄壁结构 5.2 圆顶 5.2.0 概述
1.屋顶为球面薄壳，薄壳曲面
由1/8球面构成，是由三个与 水平面夹角相等且通过球心的
大圆从球面上切割的
2.平面形状为48m*41.5m的曲 边三角形 3.壳面荷载通过薄壳的三个边 传至支座。
Central South University of Forestry & Technology
薄壳结构1
0
K k1 k2
0
0
零高斯曲率 正高斯曲率 负高斯曲率
零高斯曲率
正高斯曲率
负高斯曲率
二 、薄壳结构的曲面形式
➢旋转曲面
由一条平面曲线绕着该平面内某一指定的直线旋转一周所 形成的曲面
二 、薄壳结构的曲面形式
➢平移曲面
由一条竖向曲线做母线沿着另一条竖向曲线（导线）平行 移动所形成的曲面
顺剪力使拱方向的支座受推力 拉杆 锚于地下的斜 拉杆
➢扭壳的边缘构件 •落地拱单块扭壳屋盖
顺剪力
落地拱扭壳屋盖边框推力的平衡
合力R作用于A、B的基础上
当地基抗侧移能力不足时， 在基础之间设置拉杆
三、双曲抛物面扭壳的工程实例 ➢大连海港转运仓库
平面图
透视图
钢筋混凝土组合型双曲 抛物面扭壳屋盖
边缘构件为人字形拉杆 拱
➢壳体稳定性好 上凸——薄拱（受压）
下凹——索网（受拉）
避免屈曲失稳，提高稳定性。 壳板可以很薄
➢双曲抛物面是直纹曲面，配筋和模板制作简单
➢扭壳——从双曲抛物面中沿直纹曲面切取的一部分
扭壳单块作屋面或组合
一 、扭壳的结构组成和型式
➢组成：壳板和边缘构件 ➢形式
双倾单块扭壳
单倾单块扭壳
组合型扭壳
屋顶四边采光， 排水方便。
边缘构件为人字形拉杆 拱，壳体及边拱围现浇 刚筋混凝土结构
扭壳的几何图形
曲面的切割与组合
双曲扁壳 柱壳
两柱壳相贯 两双曲抛物面 壳相贯
锯齿形状 （柱状） 劈锥壳
锥形壳
扭壳 锥形壳
锥形壳
鸟瞰图 室外透视
➢美国圣路易航空港候机室
由三组壳体组成 每组有两个圆柱形曲面正交形成 建筑平面为四边形 每组壳面切割成八角形的覆盖平面丰富了建筑造型 两个柱形曲面的交线为十字形交叉拱，加强壳体， 并将荷载传至支座 拱的断面突出壳面使室内形成光滑简洁的曲面 壳体边缘有加劲肋 三组壳体的相交处为采光带
非线性屈曲分析方法在薄膜褶皱研究中的应用与进展
非线性屈曲分析方法在薄膜褶皱研究中的应用与进展引言：现代意义上的膜结构起源于20世纪初。
由于膜结构自重轻，透光率高，抗震性能好等优点，使得膜结构迅速发展，出现了一系列优秀的建筑作品。
1970年，在日本大阪万国博览会上，膜结构第一次集中展示并引起广泛的关注和兴趣。
1995年以后，薄膜结构在我国的应用也日益增多，规模较大的已有130多座[1]杨庆山，姜忆南. 张拉索—膜结构分析与设计[M]. 北京: 科学出版社，2004。
随着薄膜结构的广泛应用，膜材的各项性能也引起了人们的广泛关注。
膜材作为柔性材料，最重要的特性就是它的弯曲刚度特别小，其抗压缩能力很差。
这种结构，在面外荷载作用下所产生的弯矩、剪力需要通过结构的变形转换成面内拉力或压力，当压缩应力超过膜材的抗压能力时，结构上的部分节点就会偏离其原来的平衡位置，出现局部屈曲现象，即产生褶皱。
随着薄膜结构的广泛应用，褶皱带来的不利影响也就见凸显。
褶皱的产生会不仅影响建筑物的美观，更重要的是影响结构的稳定性，同时对结构的动态性能也会产生不利影响。
目前对薄膜结构褶皱研究的方法主要有两种：数值模拟方法和实验分析方法。
实验分析方法受到薄膜自身特性、实验工具、测量手段等的限制，使得目前仅能对部分简单的结构形式采用实验分析的方法进行研究，所得到得实验研究数据不但数量有限，而且只是针对几种非常简单的结构形式。
与此相比，数值方法则灵活的多。
数值分析方法不受实验空间和测量手段等的限制，可以用于计算分析大型复杂的空间结构。
基于多种数值理论的数值分析方法，已经越来越广泛的应用于薄膜结构褶皱的研究。
数值分析方法的发展平面薄膜结构褶皱数值分析方法主要有两种一种是基于薄膜理论采用不可压缩材料模型的数值分析。
该方法包含基于Stein-Hedgepeth理论的迭代薄膜性能（IMP）方法、基于张力场理论的修正变形梯度法、修正弹性张量法、二变量参数（T-VP）法、修正本构矩阵法等，基于薄膜理论的褶皱数值分析方法假定薄膜没有弯曲刚度，不能够承受压缩应力，可以确定褶皱的走向和区域。
薄壳结构调研讲解
薄壳结构班级学号：1101404-25姓名：刘益宁指导老师：彭懿日期：2013.11.20调研建筑：星海音乐厅·悉尼歌剧院·国家大剧院1薄壳结构的定义：壳，是一种曲面构建，主要承受各种作用产生的中面内的力。
薄壳结构就是曲面的薄壁结构，按曲面生成的形式分为筒壳、圆顶薄壳、双曲扁壳和双曲抛物面壳等，材料大都采用钢筋和混凝土。
壳体能充分利用材料强度，同时又能将承重与围护两种功能融合为一。
2薄壳结构的特点：壳体结构一般是由上下两个几何曲面构成的空间薄壁结构。
两个曲面之问的距离即为壳体的厚度(δ)，当δ比壳体其他尺寸(如曲率半径R，跨度等)小得多时，一般要求δ/R≤1/20(鸡蛋壳的δ/R≈1/50)称为薄壳结构。
现代建筑工程中所采用的壳体一般为薄壳结构。
而薄壳结构为双向受力的空间结构，在竖向均布荷载作用下，壳体主要承受曲面内的轴向力(双向法向力)和顺剪力作用，曲面轴力和顺剪力都作用在曲面内，又称为薄膜内力。
而只有在非对称荷载(风，雪等)作用下，壳体才承受较小的弯矩和扭矩。
由于壳体内主要承受以压力为主的薄膜内力，且薄膜内力沿壳体厚度方向均匀分布，所以材料强度能得到充分利用；而且壳体为凸面，处于空间受力状态，各向刚度都较大，因而用薄壳结构能实现以最少之材料构成最坚之结构的理．想。
由于壳体强度高、刚度大、用料省、自重轻，覆盖大面积，无需中柱，而且其造型多变，曲线优美，表现力强，因而深受建筑师们的青睐，故多用于大跨度的建筑物，如展览厅、食堂、剧院、天文馆、厂房、飞机库等。
不过，薄壳结构也有其自身的不足之处，由于体形多为曲线，复杂多变，采用现浇结构时，模板制作难度大，会费模费工，施工难度较大；一般壳体既作承重结构又作屋面，由于壳壁太薄，隔热保温效果不好；并且某些壳体(如球壳、扁壳易产生回声现象，对音响效果要求高的大会堂、体育馆、影剧院等建筑不适宜。
双曲抛物面案例星海音乐厅星海音乐厅位于广州二沙岛，造型奇特的外观，富于现代感，犹如江边欲飞的一只天鹅，与蓝天碧水浑然一体，形成一道瑰丽的风景线。
周边固支深薄球壳的非线性动态分析
果讨论 壳体 的前 三阶振动频率 与壳体各参数之 间的关 系． 结果表 明壳体展 开 角较小 时， 阶振 动的频 率大干 一阶 高 振动的频率． 向载荷对高阶振动频率的影响小于其 对一阶振 动频率的影响． 横 关键词 ：深薄球壳；打靶法 ；非线性
中 图分 类号 ：０ ２ ； ４ ３ ６ ０３ ３ 文 献 标 识 码 ：Ａ
式 中
Ｎ＝ （ 素ｄ ． ） ＝ ＋ ＝ ｒ
Ｎ ｚ＋ ） ｚ （ 素ｄ 一 １ ｚ
ｃ ７
（ ８ ）
Ｍ一 （ 素 ｄ Ｊ １ ） ． ＋
Ｍ一 （ 素 如 ２ 。＋ ） ｊ １
一
（ ９ ）
（） １ ｏ
合材料层合扁球壳的非线性强迫振动问题．
『 （＋ 兰 ￡ 胆） －
（ ｃ ｏ ｌｏ ｃｅ ｃ Ｓ ｈ ｏ ｆＳ ｉｎ ｅ，Ｉ ｚ ｕＵｎｉ． ｏ ｃ ．，ｌ ｚ ｏ ７ ０ ５ ｍｎ ｈｏ ｖ ｆＴｅ ｈ ｍｎ ｈ ｕ ３ ０ ０，Ｃｈｎ ） ｉａ
Ａｂ ｔ ａ ｔ ｓ ｒ ｃ ：Ｔａ ｉｇ ｔ ｅｇ ｏ ｔｙｎ ｎ ｉ ｅ ｒｄ ｆｒ ｔ ｎ ｉｔ ｃ ｏ ｎ ，ａｄｓ ｌｃ ｍｅ ｔｔｐ ｙ ａ ｃ ｃ ｎ ｒ ｌ ｋ ｎ ｈ ｅ ｍｅｒ ｏ ｌ ａ ｅｏ ｍａ ｉ ｎ ｏａ ｃ ｕ ｔ ｉｐａ ｅ ｎ—ｙ ｅｄ ｎ ｍｉ ｏ ｔｏ ｎ ｏ
壳 的动态 问题研究 就更 少．
点法求解线性微分方程组 ， 得到临界载荷的数值解．
王新 志和 王永 刚［研究 扁球 壳在静 载荷 作用 下 的轴 ４ ］
对称非线性 自由振动问题 ， 给出非线性静力边值问
题 的一 次 近 似解 析 解 ． ０ ５年 ， 新 志 等【 用 拟 壳 ２０ 王 ］ 法建立 正 三角形 网格 的三 向扁锥 面单层 网壳 的非 线 性 动力 学微分 方程 ． 在周 边 固定条 件下 ， 用分 离变 量
超空泡运动体圆柱薄壳的非线性动力屈曲分析
超空泡运动体圆柱薄壳的非线性动力屈曲分析王杰方;安海;安伟光【摘要】超空泡运动体的动力屈曲失稳具有隐蔽性、突发性和危险性,因而必须研究清楚运动体的失稳区域边界及失稳振幅.将超空泡运动体模拟成受轴向周期载荷作用的细长圆柱薄壳,给出非线性几何方程、物理方程和平衡方程,建立细长圆柱薄壳带有非线性项的动力屈曲微分方程组;依据非线性项的形式,给出合理的非线性位移表达式,得到具有周期性系数的非线性横向振动微分方程;采用伽辽金交分法和和鲍洛金方法,获得带有周期性系数和非线性项的马奇耶方程;求解非线性马奇耶方程,得到第一、第二阶不稳定区域内的定态振动振幅的解析表达式;绘制超空泡运动体的非线性参数共振曲线,分析航行速度、载荷比例系数、轴向载荷频率和振型对参数共振曲线的影响.以上研究为建立基于参数共振的圆柱薄壳动力失稳的可靠性分析及基于参数共振可靠性的结构动力优化设计的奠定了理论基础.【期刊名称】《力学学报》【年(卷),期】2016(048)001【总页数】11页(P181-191)【关键词】超空泡运动体;圆柱薄壳;非线性;动力屈曲;参数共振【作者】王杰方;安海;安伟光【作者单位】哈尔滨工程大学航天工程系,哈尔滨150001;南昌工程学院土木与建筑工程学院,南昌330099;哈尔滨工程大学航天工程系,哈尔滨150001;哈尔滨工程大学航天工程系,哈尔滨150001【正文语种】中文【中图分类】O343.9对于细长超空泡运动体（长细比约为10～22）而言，仅在头部和尾部与水接触，由于工作环境特殊（气、液两项介质）且航行速度很高，超空泡运动体比常规武器（单一介质）对结构的动力稳定性能更加敏感，特别是轴向周期载荷作用下的参数失稳问题（轴向载荷作用下横向振动的共振）. 只有在特定的情况下，即在横向振动频率、激发系数及外载荷频率等满足一定关系的参数平面内，参数共振才会被激发出来，但参数共振一旦被激发，其振动振幅会迅速地增长. 因此，对于高速甚至超高速运动的细长体而言，参数共振是一种突发的、危险的、隐蔽的值得引起重视的失效模式.国内外对参数振动的研究有：文献 [1-3] 利用弗吕格、桑德尔、拉甫、唐奈尔方程建立了在轴向力作用下的圆柱壳的微分方程，并转化为马奇耶方程分析了圆柱壳的动力稳定性问题. 周承倜等[4]用汉密尔顿原理导出了复合材料叠层圆柱壳的非线性动力稳定性微分方程. 曾潇等[5]对受轴向压缩的有随机初始缺陷的圆柱壳进行了动力稳定性及其可靠性分析. 张善元等[6]对圆柱壳在轴向周期载荷作用下的参数振动进行了研究. 以上的文献都是以圆柱壳作为研究对象的. 实际上，参数振动有广泛的工程背景，例如，赵晶瑞等[7]对一种新型的深海顺应式采油平台（“Spar”平台）的纵摇运动进行了分析，得到了具有三次非线性项的有阻尼马奇耶方程. 徐万海等[8]研究了张力腿、立管等海洋细长结构参数振动的动力不稳定区域. 王俊荣等[9]研究了深水半潜式平台的纵摇和横摇的参数共振，给出了半潜式平台的马奇耶方程，分析了平台失稳的参数条件. 吴学敏等[10]开展了深水顶张式立管参数振动与涡激振动耦合振动分析方法研究. 文献 [11] 研究了参数激励作用下立管的动力稳定性问题. 桑松等[12]用希尔无穷行列式法推导得到了纵摇响应的马奇耶稳定性图谱，并用数值方法分析了平台纵摇马奇耶不稳定运动的发生过程. 除了以上海洋工程结构以外，大跨度斜拉桥的斜拉索参数振动研究也十分热门[13-16].另外，张晓湘[17]分析了钢框架结构承受周期载荷时的动力稳定性. 邱良[18]研究了径向均布周期载荷作用下拱结构的动力稳定性.目前，国内外关于超空泡航行体参数共振研究公开发表的文献很少. 国外，佐治亚理工学院在超空泡航行体的参数共振方面做了一些工作：文献 [19-20]根据汉密尔顿原理建立超空泡航行体结构壳体运动方程，通过有限元离散求解航行体结构动力屈曲问题，采用鲍洛金方法得到超空泡航行体模型的失稳区域. 文献 [21] 采用有限元方法对超空泡航行体普通、加肋和渐缩壳体结构进行静力和动力屈曲研究. 文献[22] 引入自适应空化器概念，采用考虑剪切变形的梁单元建立超空泡航行体有限元模型，对航行体结构进行动力屈曲分析. 国内，施连会等[23]基于弹性壳体的一般理论，利用鲍洛金方法对航行体的动力屈曲问题进行了数值研究. 麻震宇等[24-25]基于更新拉格朗日格式的超空泡航行体结构有限元模型，将结构动力屈曲理论与非线性有限元方法相结合，开展了超空泡航行体双层壳结构动力屈曲研究.宋向华等[26-27]建立了超空泡射弹截锥形结构的动力微分方程，利用鲍洛金方法进行动力稳定性分析，求解出动力不稳定区域边界; 随后，给出射弹结构动力稳定性的安全余量方程，计算了结构的动力稳定性的非概率可靠度. 刘明等[28-29]在超空泡射弹动力稳定性分析的基础上，采用随机因子法求出随机参数射弹结构的动力不稳定区域边界，并对结构的可靠性进行了分析. 王杰方等[30]将超空泡运动体模拟成受动态轴向载荷作用的圆柱薄壳，推导了结构的动力稳定性微分方程，计算了舱段动力屈曲的可靠性指标.从现有的文献来看，基于参数共振的高速航行体圆柱薄壳结构非线性动力失稳问题的研究文献较少. 较大的运动体横向变形，会直接破坏原有的空泡稳定性，而空泡的稳定性决定了高速航行体的运动稳定性. 线性的参数共振理论只适用于小变形的情况，随着横向振幅的增长，非线性因素的影响开始显露出来，成为不可忽略的必要因素. 因此，本文将从非线性的角度开展基于参数振动的失稳研究.文中将线应变和剪应变中计入非线性因素，建立圆柱薄壳横向振动的非线性微分方程; 采用伽辽金变分法和鲍洛金方法，将非线性微分方程转化为具有周期性系数和非线性项的马奇耶方程; 求解非线性马奇耶方程，得到考虑非线性因素细长圆柱薄壳的第一、第二阶不稳定区域内的定态振动的振幅解析表达式，进而对影响圆柱薄壳的非线性参数共振曲线的因素进行分析.将超空泡运动体模拟成受轴向周期载荷作用的细长圆柱薄壳. 对于圆柱薄壳，采用未变形前的曲线坐标系（α,β,z ）来描述有限变形的几何位置，α,β 为壳体中面的主曲率线.在扭转变形的几何关系中计入切向位移的影响，在线应变和剪应变中计入非线性项，细长圆柱薄壳的非线性几何方程为式中，（u, v,w ）为轴向、周向和径向位移，R 为圆柱薄壳的半径为中面内的线应变为中面内的剪切变形为弯曲变形，κ12 为扭转变形.物理方程中，中面力中计入弯曲变形的影响，中面外的力 M 1中计入中面应变的影响，即式中，E为材料弹性模量，h为薄壳厚度，v为泊松比，是中面内的法向力和剪力，是中面外的弯矩和扭矩.平衡方程中，在圆周切线方向（周向）考虑横向剪切力Q2的影响，即式中，X为沿圆柱薄壳轴向的表面力，Y为沿周向的表面力，Z为沿法向的表面力，为横向剪切力.综合以上3个方程，得到细长圆柱薄壳的非线性微分方程组为方程组（4）是关于位移的3个微分方程组，相关微小项已省略，为非线性项，仅与位移w有关，式中式中为了求得方程组（4）的近似解，在轴向周期载荷作用下的圆柱薄壳，设其横向振动的弯曲形式为式中圆柱薄壳两端没有法向位移为轴向和周向的半波数，取整数值，L为圆柱薄壳的长度. 式（10）给出了在均布轴向载荷作用下的圆柱薄壳线性问题的精确解，将其代入中得设方程组（4）中前两式的解的形式为式中满足线性方程组且应该满足边界条件. 将式（14）代入式（4）的前两式（略去切向表面力，），得到系数的表达式为联立几何方程式（1）、物理方程式（2）、式（14）和式（16），得到轴向力和周向力为式中设的表达式为既可以满足边界条件（径向和周向位移为0，轴向不为0），又可以满足方程组（15）. 将式（19）代入式（15）中得对于受轴向周期载荷作用的圆柱薄壳，其轴向的平衡条件为式中为不随时间变化的部分为随时间变化的部分，θ为动载荷的角频率.联立式（17）～式（21），求得内力为横向振动微分方程组（4）的第3式中，沿圆柱薄壳法向的表面力为[31]式中，m是圆柱薄壳中面内单位面积的质量，m＝ρh.综合式（4）中的第3式、式（13）、式（14）、式（16）、式（19）、式（20）、式（22）和式（23），得到圆柱薄壳的非线性横向振动微分方程式，即式中是振型为（n,k）时圆柱薄壳的自由振动固有频率式中g'（n,k）与文献 [30] 中的表达式完全一致，经验证，文献 [30] 中的线性问题得到的自由振动固有频率与非线性问题中的表达式一致，这说明非线性微分方程式（24）中的线性部分与文献 [30]中的线性微分方程式是一致的，也进一步验证了圆柱薄壳非线性横向振动微分方程式（24）的正确性.对圆柱薄壳的非线性横向振动微分方程式（24）进行伽辽金变分，得到非线性横向振动微分方程为引入临界载荷的记号 P∗，令这一记号形如文献 [30] 中的临界载荷另外，引入非线性系数的记号令于是，横向振动微分方程写成下面的形式，即式（31）的微分方程是带有周期性系数和非线性项的二阶齐次微分方程，用于描述细长圆柱薄壳在轴向周期载荷作用下的参数振动问题. 按照鲍洛金方法，将式（31）写成马奇耶方程的形式式中，为在轴向力的定值分量作用下的横向振动的固有频率，为激发系数，即采用鲍洛金方法确定圆柱薄壳在动力不稳定区域内的定态振动振幅.2.1 第一阶不稳定区域内的定态振动振幅为了确定第一阶不稳定区域内的定态振动振幅，取将式（34）代入式（32），得到包含系数和的方程组式中，系数和的表达式为将式（34）中的解代入式（36）中，忽略谐波项，获得式（35）中的两个非线性项和的表达式，并将式（35）变换为式中表示第一共振区内的定态振幅. 显然时，方程组（37）是成立的，相当于受轴压作用的圆柱薄壳没有横向振动的情况.为了获得方程组（37）的非零解，将其视为和的线性齐次方程组，只有当未知数的系数组成的行列式等于零时和才有非零解，即解出第一阶不稳定区域内的定态振幅 A1为马奇耶方程具有以下性质[31]：不稳定区域和稳定区域被周期为T和2T的周期解隔开了，也就是说，周期相同的两个解包围着不稳定区域，周期不同的两个解包围着稳定区域. 因此，第一阶不稳定区域由振幅为和周期为的两个周期解所包围. 2.2 第二阶不稳定区域内的定态振动振幅为了确定第二阶不稳定区域内的定态振幅，取将式（40）代入式（32），得到包含系数 a2和 b0，b2的方程组式中和按式（36）计算，系数为将式（40）代入式（36）和式（42）中，忽略各项谐波，获得式（41）中的 3 个非线性项和并将式（41）变化为式中，为第二共振区附近的定态振幅.为获得方程组（43）非零解，将其视为和的线性齐次方程组，只有当未知数的系数组成的行列式等于零时，和才有非零解（由于故略去式（43）中的微小项[31]），将行列式等式变换为解出第二阶不稳定区域内的定态振幅的解为计算表明，在不稳定区域内，式（45）中的始终为虚数，因此，第二阶不稳定区域内的定态振动振幅为和因此，第二阶不稳定区域由振幅为和周期为T的两个周期解所包围.由式（39）和式（45）可知，在考虑几何非线性因素的圆柱薄壳非线性动力屈曲问题中，动力不稳定区域边界上周期解的定态振动振幅会随着外载荷频率变化而变化，而在线性问题中，认为动力不稳定区域边界上为无限增长的解[31]. 因此，动力屈曲的非线性分析能定量地给出超空泡运动体在不同外载荷频率下的参数振动振幅，进而为运动体尾部的浸水深度和沾湿面积的确定以及水动力分析打基础.超空泡运动体航行深度流场密度对于自然超空泡，空泡内的饱和蒸汽压标准大气压圆柱薄壳舱段的几何参数：半径长度厚度空化器直径材料物理参数[32]：弹性模量材料密度泊松比3.1 超空泡运动体受力分析对于超空泡运动体水平向前运动时的动力屈曲分析，其受力可以简化为轴向的均布载荷，即头部阻力和尾部推力，二者大小相等，作用在头部空化器的阻力[21]为式中为流体密度是空化器的横截面积是空化器的阻力系数，V是航行体的运动速度. 零攻角时，圆盘空化器的阻力系数为式中为空化器锥角，圆盘空化器的锥角为180°为空化数为空泡内压力为环境压力，其表达式为式中，H为航行深度，p为标准大气压.运动体在高速航行过程中，空泡形状和尺寸的不稳定性会导致轴压幅值随时间变化，本文将这一动态的轴向载荷简化为[19]式中，δ 为扰动载荷的比例系数，dn为圆盘空化器的直径.3.2 圆柱薄壳舱段非线性动力屈曲计算结果及分析综合第2节和第3节中的内容可知，动力不稳定区域附近的定态振幅不仅与横向振动的振型有关，与轴向载荷频率有关，还与载荷比例系数δ和航行速度V有关. 因此，下文将分析这些因素对非线性参数共振曲线的影响.3.2.1 给定速度和载荷比例系数时，不同振型下的非线性参数共振曲线图1、图2和图3给出了一定速度和载荷比例系数时，不同振型下，考虑几何非线性因素的第一阶非线性参数共振曲线图，其中虚线表示的是不稳定解. 从式（39）可知，若则以振型应于图1中的 G 点对应于图1中的 C 点. 所以，某一振型对应的两条曲线（一条实线和一条虚线）所夹的横坐标上的区域即为这一振型下的激发区.从以上 3 幅图中可以得出以下结论：（1）当圆柱薄壳所受的外载荷的频率处于激发区以外且在小于激发区一侧（图1中的 BC 段）时，圆柱薄壳没有横向振动.（2）考虑几何非线性的圆柱薄壳舱段，其各阶参数共振曲线都向大于激发区频率的一侧倾斜. 因此，① 外载荷频率从小于激发下界的一侧开始逐渐增大（图1中的路线为BCGF）时，圆柱薄壳的横向振动振幅沿着 BCHD 增大.② 外载荷频率从大于激发上界的一侧开始逐渐减小（路线为FGCB）时：若从 F 到 G 点没有外界干扰或者干扰不足以使壳体的横向振动越过不稳定解（EG 曲线）到达稳定解（DH 曲线）上去，则壳体的横向振动振幅依然为0，直到频率到达 G 点时，壳体才会发生“突变”的、振幅位于 H 点的横向振动.若从 F 到 G 点时，存在干扰使得壳体的横向振动越过不稳定解（EG 曲线）到达稳定解（DH 曲线）上去，横向振动的振幅由零突然增大到 DH 曲线，这使得 GF 段成为潜在的激发区域，导致危险的激发区域扩大.（3）对比图1、图2和图3的纵坐标可知，当 k＝2或 k＝3 时，不稳定区域内的定态振动振幅与壳体的厚度（h＝3 mm）是同一个数量级的，但当 k＝1 时，不稳定区域内的定态动振幅比k＝2 或 k＝3 时高出两个数量级.以振型 i＝1，k＝1为例，对比参数振动的线性和非线性理论，可知：（1）在参数振动的线性理论中，不存在参数共振曲线的“倾斜”问题，激发区域只有图1中的 CG 段，但在考虑几何非线性的参数振动非线性理论中，由于参数共振曲线向大于激发区频率一侧“倾斜”，导致激发区不仅有 CG 段，还包括 GF 段. 因此，考虑几何非线性后，激发区的扩大会导致外载荷安全频率范围的缩小和结构动力稳定可靠性的降低.（2）参数振动的线性理论认为激发区域内振动的振幅是无限大的，这不符合工程实际. 考虑几何非线性，激发区域内的振动振幅实际上为有限振幅. 对于超空泡运动体而言，结构的变形量直接影响其沾湿面积，而沾湿面积决定了空泡的稳定性，这是一个典型的流固耦合问题. 因此，考虑非线性因素，准确地获得激发区域内的参数振动振幅，是进行流固耦合分析的必要前提.3.2.2 给定的振型（i＝2，k＝2）时，速度和载荷比例系数对非线性参数共振曲线的影响图4给出了载荷比例系数δ＝0.4 时，不同航行速度V下的第一阶非线性参数共振曲线图. 图5给出了航行速度 V＝400 m/s 时，不同载荷比例系数δ下的第一阶非线性参数共振曲线图.由图4和图5可以得出以下结论：（1）激发区的范围随着超空泡运动体的航行速度V和载荷比例系数δ增大而增大; 激发区内的定态振动的振幅随着超空泡运动体的航行速度V和载荷比例系数δ增大而增大.（2）结论“（1）”与线性理论的结论相同[30]，事实上，从式（38）和文献 [31] 中的线性临界频率方程式的对比中可知，若令式（38）中的那么，式（38）和线性临界频率方程式是一致的，也就是说，图1中的CG段的激发区域在线性理论和非线性理论中是相同的. 不同之处在于，非线性理论中要将GF段也计入激发区域范围内，使激发区比线性理论中给出的激发区扩大.通过对第二阶非线性参数共振的共振曲线进行分析，也能得出与第一阶相同的结论. 图6给出了第一阶和第二阶非线性共振曲线的对比图（δ＝0.4，V＝400 m/s）. 由图6可知，第二阶非线性共振曲线与第一阶的形状相似，都向较大频率的一侧倾斜; 第二阶不稳定区域内的定态振动的振幅远小于第一阶的定态振动的振幅，这也说明了在圆柱薄壳的横向参数振动中，第一阶不稳定区域不仅激发区远大于第二阶，考虑非线性因素后，第一阶激发区内的振动振幅也远大于第二阶.本文建立了圆柱薄壳的非线性动力屈曲计算模型，并给出了超空泡运动体的圆柱薄壳舱段的分析算例，得到以下结论：（1）考虑几何非线性因素后，参数共振曲线向大于激发区频率的一侧倾斜，使得激发区的范围变大，从而导致外载荷安全频率范围的缩小和结构动力稳定可靠性的降低.（2）考虑非线性因素，能准确地获得激发区域内的参数振动振幅，是进行超空泡运动体流固耦合分析的基础.（3）激发区的范围随着超空泡运动体的航行速度V和载荷比例系数δ增大而增大; 激发区内的定态振动的振幅随着超空泡运动体的航行速度 V和载荷比例系数δ 增大而增大.（4）当 k＝2 或 k＝3 时，不稳定区域内的定态振动振幅与壳体的厚度（h＝3 mm）是同一个数量级的，但当k＝1 时，不稳定区域内的定态振动振幅比k＝2 或 k＝3 时高出两个数量级.【相关文献】1 Lam KY，Ng TY. Dynamic stability of cylindrical shells subjected to conservative periodic axial loads using different shells theories. Journal of Sound and Vibration，1997，207（4）: 497-5202 Ng TY，Lam KY，Liew KM，et al. Dynamic stability analysis of functionally graded cylindrical shells under periodic axial loading. International Journal of Solids and Structures，2001，38: 1295-13093 Darabi M，Darvizeh M，Darvizeh A. Non-linear analysis of dynamic stability for functionally graded cylindrical shells under periodic axial loading. Composite Structures，2008，83: 201-2114 周承倜，王列东. 复合材料叠层圆柱壳的非线性动力稳定性分析.大连大学学报，1993，01: 1-15（Zhou Chengti，Wang Liedong. Nonlinear analysis of dynamic stability for laminated composite cylindrical shells. Journal of Dalian University，1993，01: 1-15（in Chinese））5 曾潇，卓曙君. 有随机初始缺陷的轴压圆柱薄壳的动力稳定性分析. 昆明理工大学学报，1996，21（3）: 59-63（Zeng Xiao，Zhuo Shujun. Dynamic stability analysis of axially compressed cylindrical shell with random initial geometric imperfections. Journal of Kunming University of Science and Technology，1996，21（3）: 59-63（in Chinese））6 张善元，张涛. 圆柱壳的轴向动力屈曲、参数共振与混沌运动. 振动与冲击，2010，29（12）:34-38，66（Zhang Shanyuan，Zhang Tao. Axial dynamic buckling parametric resonance and chaotic motion of a closed cylindrical shell. Journal of Vibration and Shock，2010，29（12）: 34-38，66（in Chinese））7 赵晶瑞，唐友刚，王文杰. 传统 Spar 平台参数激励 Mathieu 不稳定性的研究. 工程力学，2010，27（3）: 222-227（Zhao Jingrui，Tang Yougang，Wang Wenjie. Study on the parametrically excited Mathieu instability of a classic Spar platform. Engineering Mechanics，2010，27（3）: 222-227（in Chinese））8 徐万海，吴应湘，钟兴福等. 海洋细长结构参数激励不稳定区的确定方法. 振动与冲击，2011，30（9）: 79-83（Xu Wanhai，Wu Yingxiang，Zhong Xingfu，et al. Methods for parametric excitation instability analysis of slender flexible cylindrical structures in offshore engineering. Journal of Vibration and Shock，2011，30（9）: 79-83（in Chinese））9 王俊荣，谢彬. 深水半潜式平台Mathieu不稳定问题研究. 工程力学，2012，29（10）: 347-353（Wang Junrong，Xie Bin. Mathieu instability study of a deepwater semi-submersible platform. Engineering Mechanics，2012，29（10）: 347-353（in Chinese））10 吴学敏，黄维平，滕文刚. 深水顶张式立管参数振动与涡激振动耦合振动分析方法研究. 中国海上油气，2014，04: 100-105（Wu Xuemin，Huang Weiping，Teng Wengang. Study on analysis method for coupled vibration of parameter excited vibration and vortex-induce vibration on deep water top-tensed riser. China Offshore Oil and Gas，2014，04: 100-105（in Chinese））11 Zhang J，Tang YG. Mathieu instability analysis of deep-water toptensioned Risers. Journal of Ship Mechanics，2014，09: 1142-115012 桑松，石晓，李长东等. 深海SPAR平台垂荡-纵摇耦合运动Mathieu稳定性分析. 中国海洋平台，2014，04: 34-40（Sang Song，Shi Xiao，Li Changdong，et al. Study on Mathieu stability of heave-pitch coupled motions of SPAR platform in deep sea. China Offshore Platform，2014，04: 34-40（in Chinese））13 李为洲. 斜拉索面内参数振动数值分析及其控制. [硕士论文]. 长沙: 湖南大学，2013（Li Weizhou. Numerical analysis of in-plane parametric vibration and the controlling measures for stay cables. [Master Thesis]. Changsha: Hunan University，2013（in Chinese））14 孙超. 大跨度铁路斜拉桥拉索参数振动研究. [硕士论文]. 成都: 西南交通大学，2014（Sun Chao. Parametric vibration of cables for long-span railway cable-stayed bridges. [Master Thesis]. Chengdu: Southwest Jiaotong University，2014（in Chinese））15 李周. 大跨度悬索桥悬吊体系参数振动研究. [博士论文]. 广州: 华南理工大学，2013（Li Zhou. The research on parametric vibration of suspended system of large-span suspension bridge. [PhD Thesis]. Guangzhou: South China University of Technology，2013（in Chinese））16 杨咏漪，陈克坚. 大跨度铁路斜拉桥斜拉索参数振动分析. 铁道工程学报，2013，10: 60-65（Yang Yongyi，Chen Kejian. Research on parametric oscillation of cables for long span railway cable-stayed Bridge. Journal of Railway Engineering Society，2013，10: 60-65（in Chinese））17 张晓湘. 竖向周期荷载作用下单层钢框架的动力稳定性分析. [硕士论文]. 长沙: 中南大学，2013（Zhang Xiaoxiang. Dynamic stability analysis of single-story steel frames under vertical cyclic load. [Master Thesis]. Changsha: Central South Univerdity，2013（in Chinese））18 邱良. 径向均布周期荷载作用下拱的动力稳定性能研究. [硕士论文]. 长沙: 中南大学，2013（Qiu Liang. Dynamic stability research of arch under radial uniform cyclic loading. [Master Thesis]. Changsha: Central South Univerdity，2013（in Chinese））19 Ruzzene M. Dynamic buckling of periodically stiffened shells: appli- cation to supercavitating vehicles. International Journal of Solids and Structures，2004，41: 1039-105920 Ruzzene M. Non-axisymmetric buckling of stiffened supercavitating shells: static and dynamic analysis. Computers and Structures，2004，82: 257-26921 Ahn SS，Ruzzene M. Optimal design of cylindrical shells for enhanced buckling stability: Application to supercavitation underwater vehicles. Finite Elements in Analysis and Design，2006，42: 967-97622 Choi JY，Ruzzene M. Stability analysis of supercavitating underwater vehicles with adaptive cavitator. International Journal of Mechanical Sciences，2006，48: 1360-137023 施连会，王安稳. 轴向载荷下超空泡航行体动力稳定性的数值研究. 振动与冲击，2011，30（2）: 55-59（Shi Lianhui，Wang Anwen. Numerical study on dynamic stability of supercavitating vehicles subjected to axial loads. Journal of Vibration and Shock，2011，30（2）: 55-59（in Chinese））24 麻震宇. 超空泡航行体结构动力学仿真研究. [博士论文]. 长沙: 国防科学技术大学，2012（Ma Zhenyu. Numerical research on dynamics of structures for supercavitating vehicles. [PhD Thesis]. Changsha: National University of Defense Technology，2012（in Chinese））25 麻震宇，胡凡，陈广南等. 超空泡航行体双层壳结构动力稳定性分析. 国防科技大学学报，2011，33（4）: 43-47（Ma Zhenyu，Hu Fan，Chen Guangnan，et al. Dynamic stability analysisof supercavitating double shells. Journal of National University of Defense Technology，2011，33（4）: 43-47（in Chinese））26 宋向华，安伟光，刘明. 轴向周期载荷下超空泡射弹的动力稳定性分析. 哈尔滨工程大学学报，2012，33（10）: 1238-1243（Song Xi-anghua，An Weiguang，Liu Ming. Dynamic。
一种浅壳的非线性屈曲分析_朱寻焱
壳稳定容许承载力( 荷载取标准值) 应等于网壳稳定
极限承载力除以安全系数 K。弹塑性全过程分析安全
系数 K 取 为 2. 0; 弹 性 全 过 程 分 析 安 全 系 数 K 取 为
4. 2。按 L /300 考虑初始缺陷，网壳的跨度为 27m，初
始缺陷按竖向最大位移为 90mm 考虑。
2. 1 钢材本构关系
W1x = 1525 × 103 mm3
N+
βmx Mx
φxA
γx W1x (
1
－
0.
8
N) N' Ex
= 362 × 103 0. 241 × 15381
欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁欁
无洞室存在时的自由场地表位移幅值明显不同，对地 表地震动有一定的 放 大 作 用，另 外 入 射 波 的 入 射 角 度 和入射频率对地表位移幅值也有显著的影响。
32
低温建筑技术
2012 年第 3 期( 总第 165 期)
一种浅壳的非线性屈曲分析
朱寻焱
( 天津大学建筑工程学院， 天津 300072)
【摘 要】 介绍了一种浅壳的非浅性分析过程，依次进行了模态分析、特征值屈曲分析、引入初始缺陷弹性
极限承载力分析、引入初始缺陷弹塑性极限承载力分析。通过分析可知，浅壳的稳定性满足 JGJ7 － 2010《空间网
从图 9 可以看出结构的最大塑性应变约为 0. 0084，主要分布于腹杆上，结构有较高的材料强度安 全储备。 3 结语
通过对本网壳的计算与分析，可得到如下结论: 对本网壳 分 别 进 行 了 引 入 初 始 缺 陷 几 何 非 线 性 的弹性极限承载分 析，引 入 初 始 缺 陷 同 时 考 虑 几 何 非 线性和材料非线性 的 弹 塑 性 极 限 承 载 分 析，计 算 分 析 结果均满足 JGJ7 － 2010《空间风格结构技术规程》设 计要求，结果汇总于表 3。
旋转几何非线性复合材料薄壁梁的自由振动分析
第３ ２卷 第 １ ４期
动
与
冲
击
Ｖｏ １ ． ３ ２ Ｎｏ ． １ ４ ２ ０ １ ３
ｊ ｏＵＲＮＡＬ ＯＦ ＶＩ ＢＲＡＴ Ｉ ＯＮ ＡＮＤ Ｓ ＨＯＣＫ
旋 转 几 何 非线 性 复 合材 料 薄 壁梁 的 自由振 动分 析
任勇生 ，代其义 ，孙 丙磊 ，张纯金
（ 山东科技大学 机械 电子工程学院 ， 青岛 ２ ６ ６ ５ ９ ０ ）
ｌ ａ ｒ ｇ ｅ ｄ ｅ ｆ ｏ ｍ ａ ｒ ｔ ｉ ｏ ｎ ｗ ａ ｓ ｓ ｔ ｕ ｄ ｉ ｅ ｄ ｈ ｅ ｒ ｅ ．Ｔ ｈ ｅ ｇ ｏ ｖ ｅ ｒ ｎ ｉ ｎ ｇ ｎ ｏ ｎ ｌ ｉ ｎ ｅ ａ ｒ ｅ ｑ ｕ ａ ｔ ｉ ｏ ｎ ｓ ｏ ｆ ｍｏ ｔ ｉ ｏ ｎ ｏｒ ｆ ｔ ｈ ｅ ｒ ｏ ｔ ａ ｔ ｉ ｎ ｇ ｃ ｏ ｍｐ ｏ ｓ ｉ ｔ ｅ ｔ ｈ ｉ ｎ — ｗａ ｉ ｌ ｅ ｄ
ｉ ｎ ｔ ｈｅ ｏ ｒ ｙ ａ ｎ ｄ ｐ ｒ ａ ｃ ｔ ｉ ｃ ｅ ． Ｈｏ ｗｅ ｖ ｅ ｒ ，Ｓ Ｏ ｆ ａ ｒ ｔ ｈ ｅ ｄ ｙ ｎ ａ ｍｉ ｃ ｂ ｅ ｈ ａ ｖ ｉ ｏ ｒ ｓ ｔ ｕｄ ｙ ｏ ｎ ｔ ｈ ｅ ａ ｂ ｏ ｖ ｅ— ｍｅ ｎ ｔ ｉ ｏ ｎ ｅ ｄ ｓ ｔ ｒ ｕ ｃ ｔ ｕ ｒ ｅ ｉ ｓ ｌ ｉ ｍｉ ｔ ｅ ｄ ｏ ｎ ｌ ｙ ｆ ｏ ｒ ｒ ｏ ｔ ａ ｔ ｉ ｎ ｇ ｃ ｏ ｍｐ ｏ ｓ ｉ ｔ ｅ ｂ ｅ ａ ｍｓ ｗｉ ｔ ｈ ｓ ｍａ ｌ ｌ ｅ ｌ ａ ｓ ｔ ｉ ｃ ｄ ｅ ｆ ｏ ｒ ｍａ ｔ ｉ ｏ ｎ．Ｔｈ ｅ ｆ ｒ ｅ ｅ ｖ ｉ ｂ ｒ ａ ｔ ｉ ｏ ｎ ｏ ｆ ａ ｒ ｏ ｔ ａ ｔ ｉ ｎ ｇ ｃ ｏ ｍｐ ｏ ｓ ｉ ｔ ｅ ｔ ｈ ｉ ｎ － ｗａ ｌ ｌ ｅ ｄ ｂ ｅ ａ ｍ ｗｉ ｔ ｈ
(八)薄壳结构
• • 圆顶薄壳的曲面是球形曲面，故也称“球面薄壳”。 圆顶薄壳是最早出现的一种古老型式，它的受力简单。壳板的径向和环向 弯矩极小，可以忽略。所以，壳体主要承受压力，压力沿整个球面扩散均 匀分布，可以充分利用材料的强度。边缘构件为一圆环，对圆顶壳面起箍 的作用，圆顶通过它搁在支承构件上，边缘构件主要承受环向拉力和弯矩。 根据壳板构造，圆顶薄壳可以分为：平滑圆顶、肋形圆顶和多面圆顶三种。
（二）长壳受力特点和内力计算简介
薄壳是空间结构，在荷载作用下产生的内力与普通梁板结构不 同，内力的计算也比普通梁板结构复杂得多。
二、短壳
短壳一般也是由壳板、边梁和横隔构件三部分组成，其跨度与波长的比值l1/ l2<1，通常等于或小于0.5。短壳一般是多跨的。 （一）壳板 壳板的矢高f1不应小于1/8l2。壳板内的应力不大，通常不必计算，可按跨度及 施工条件决定其厚度。对普通跨度（l1＝6～12 m， l2＝18～30 m）的屋盖，当矢 高不小于1/8l2时，厚度可按表9－1选定。 （二）边梁 边梁宜采用矩形截面，其高度一般为（1/15-1/10）l1，而不应小于1/15 l1，宽 度为高度的1/5-2/5。
现代建筑结构欣赏
薄壳结构
嘉庚学院土木工程系 许云燕
第一节
薄壳结构的受力特点与优缺点
• 薄壳结构概念:曲板的内力改变为受压为主，受压比受弯有利得多，尤其 • 薄壳结构的优点：
是多向受压，处于空间工作状态更加有利薄壳必须具备两个条件：一是“曲面 的”，二是“刚性的”。所以，薄壳可以简单理解为四边支承的曲板。薄壳的结 构组成一般包括：曲面的壳板和周边的边缘构件两部分。
3.锥面与锥状面 锥面是一直母线沿一竖向曲导线移动，并始终通过一定点而 形成的曲面，如锥面壳。锥状面是由直母线沿一根直导线和一根竖 向曲导线移动，并始终平行于一导平面而形成的曲面。
均匀轴压圆柱薄壳的稳定性分析
图&
圆柱壳屈曲模态
" 非线性有限元法全过程分析
仍以图 & 所示圆柱壳为计算模型，采用 1’4 单元并划分与前面相同的单元网格，初始压力仍为 78 ’! 6 !*7 9 ，利用 -.-/01 程序进行非线性分析，分析中采用 4=>? 法 + ’ , 。为完成这种计算，需对壳
!"
哈
尔
滨
建
筑
大
学
学
报
第 ## 卷
56789/’8: NE< D-5==9D5- 82DP-94F >E<,3?I -94<53 54K 4,4-94<53 6949>< <-<Q<4> Q<>E,K= 53< <QS-,?<K >, 545-?W< >E< 82DP-94F ,6 >E94 D?-94K39D5- =E<--= 24K<3 2496,3Q 5R95- D,QS3<==9,4) NE< 949>95- 82DP-94F -,5K= S3<K9D5><K 8? >E< >E3<< Q<>E,K= 53< 5SS3,R9Q5><-? >E< =5Q<) .,4-94<53 6949>< <-<Q<4> =>2K? 3<X ;<5-= >E5> >E<3< 9= 5 4<53-? # E,39W,4>5- -<4F>E 94 >E< S,=> # 82DP-94F -,5KXK<6-<D>9,4 =S5D<I D,33<=S,4KX 94F >, >E< 4<53-? # D,4=>54> S,=> # 82DP-94F -,5K) C4 <QS939D5- 6,3Q2-5 >, S3<K9D> >E< 82DP-94F =>3<== 9= K35Y4 63,Q =>5>9=>9D= ,4 <RS<39Q<4>5- K5>5I D,QS539=,4= 8<>Y<<4 >E< 3<=2->= 63,Q >E< 6,3Q2-5 54K >E< ,4<= 63,Q 4,4-94<53 6949>< <-<Q<4> =>2K? -<5K >, >E< D,4D-2=9,4 >E5> >E< <RS<39Q<4>5- 82DP-94F =>3<== 9= 5D>25--? >E< >E<,3<>9D5- S,=> # 82DP-94F =>3<==I 35>E<3 >E54 >E< 949>95- 82DP-94F =>3<==) NE9= 9= 5 3<5=,4X 58-< <RS-545>9,4 >, >E< K9=>23894F 65D> >E5> >E< <RS<39Q<4>5- 82DP-94F -,5K= Y<3< ,6><4 Q2DE 8<-,Y >E< S3<K9D5>9,4= ,6 D-5==9D5- >E<,3?) ;3< =>917: >E94 D?-94K39D5- =E<--M 82DP-94FM 6949>< <-<Q<4>M 4,4 # -94<539>?；<QS939D5- 6,3Q2-5
机械与运载工程学院非线性非线性的来源
16
机械与运载工程学院
牛顿-拉普森求解技术 • 第一次迭代 (i=1) – 假定前面收敛增量步的解u0， P0 ，为已知的。 – 在当前增量步中，将一个小 的增量载荷P施加到结构上。 – Abaqus基于u0处的切线刚度 K0 确定位移修正c1 ；前一增 量步结束时，总载荷PTOTAL和 内力间的关系为：
20
机械与运载工程学院
一般每个分析步(*STEP)需要几个增量步。 每个增量步包含若干迭代步。
21
机械与运载工程学院
几何非线性静力学分析
机械与运载工程学院
问题描述
确定图示斜板跨中的扰度
平面图
正视图
学习要点： （1）应用壳单元；（2）几何非线性
2
机械与运载工程学院
壳单元
壳单元
3-D continuum
surface model
当结构厚度远小于其他方向尺度 （<1/10），并忽略厚度方向应力时， 可用壳单元近似模拟三维连续体。可 较精确地模拟弯曲和面内变形。
K0 c1 P TOTAL I 0 .
17
机械与运载工程学院
– Abaqus更新模型的状态为u1 ，形成K1 并计算 I1 。
– 总载荷PTOTAL与内力 I1的差称为残差， R1: R1= PTOTAL I1.
– 如果 R1 在模型的每个自由度上都非常小 (在容差范围之内)， 结构就是平衡的。 • 默认的容差R1必须小于结构对时间平均力的0.5%。 • 位移修正值c1是否小于增量位移的1%。 – 如果迭代不能得到收敛的解，Abaqus执行另外的迭代，以 找到收敛的解。
7
机械与运载工程学院
非线性
非线性的来源 材料非线性
12钢筋混凝土薄壳理论 (续)
1.054/1.541 混凝土结构力学与设计 (3-0-9)内容提要12钢筋混凝土薄壳结构薄壳{定义－薄壳是指曲面板，其厚度h与其它尺寸和主曲率半径相比很小。
{中间面将壳体平分为两半的面称为中间面。
它指明了该面的形式及任意点的厚度。
{薄壳的分析包括以下步骤：从壳体上截取一微元体，建立该微元体的平衡方程。
满足应变协调，以使得变形后每个微元体与其相邻的微元体保持连续变形。
{应力的合力及应力偶壳理论{基尔霍夫-乐夫理论－壳的一次近似假定：(1) 与壳中间面的最小曲率半径相比，壳的厚度可忽略。
(2) 壳内的应变和位移很小。
(3) 变形前垂直于中间面的直线变形后仍保持直线，并与中间面垂直，直线的长度不变。
(与梁的内维尔假说类似－梁的伯努利-欧拉理论)(4) 与壳中间面垂直的应力可忽略。
假定的推论：－参考面的法线方向在变形后保持直线，并与变形后的参考面垂直。
－假定不包含任何横向剪应变，也即，面的法线与面内的任意直线维持直角关系。
－仅适用于薄壳。
－无法描述局部荷载或连接附近的性能。
（假定(4)在集中横向力附近失效）{ Flüggy-Byrne理论－壳的二次近似假定：－仅采用假定(2)。
－ 该假定也称为尔霍夫-乐夫假定的“高阶近似”。
壳的分类{ 按几何方程分：旋转抛物面旋转双曲面柱面椭圆抛物面双曲抛物面锥面壳的几何分析{ 正交曲线坐标考虑位置向量：123(,)(,)(,)(,)r f f f αβαβαβαβ=++u v w其中，1(,)f αβ,2(,)f αβ和3(,)f αβ是单值连续函数。
面的位置由α和β确定。
α和β称为曲线坐标。
u , v 和w 为笛卡尔坐标系中的单位向量。
正交可表示为：0r r αβ∂∂⋅=∂∂Æ 内积为零 点(,)αβ和点(,)d d ααββ++间的距离为：r r ds d d αβαβ∂∂=+∂∂ ds 的自标量积：2222222r r r r ds ds d d A d B d ds αβαβααββ ∂∂∂∂ ⋅=⋅+⋅=+= ∂∂∂∂其中，A 和B 称为Lame 参数或量测数。