应用分解刚度法分析金属点阵夹芯板屈曲问题
ABAQUS屈曲分析

1
精品课件
几何非线性
非线性的来源:
➢ 几何非线性 • 大位移、大转动、大变形
➢ 材料非线性 • 非线性弹性、塑性、 • 损伤、失效……
➢ 边界非线性 • 接触、摩擦
2
精品课件
几何非线性
No Image
21
几何非线性的来源:
1 位移增量和应变增量之间的非线性关系 (应变矩阵); 2 针对当前未知体积 V 积分,不满足弹性理论中的小变形假定;
➢ 在这类分析中,载荷-变形的响应表现出负刚度的特点,并且必须释放一定应 变能来维持结构平衡
25
精品课件
静态后屈曲分析
• 为避免稳定力(Stabilizing Forces)效应,可以 在不施加稳定力的前提下,对静态平衡方程进行求 解。
• 在该算法中,荷载的施加是自动调整的 ➢ 同时求解荷载和位移
➢ 需要考虑非线性的位移-应变关系
3 1
1
7
变形后的网格
7
原始网格
5
精品课件
几何非线性
实例2:框架结构的整体失稳分析
➢ 结构的稳定性是工程分析及设计人员经常面对的问题; ➢ 该实例中,在矩形截面框架的角点处施加点荷载,分析其后屈曲行为。
A
A
矩形横截面 线弹性材料 端点铰接
A-A截面
6
精品课件
几何非线性
载达到临界值后,刚度突然大幅降低。
8
精品课件
特征值屈曲分析
欧拉柱的荷载-位移响应
9
精品课件
特征值屈曲分析
欧拉柱的变形
10
精品课件
特征值屈曲分析
特征值屈曲分析
➢ 分析结构刚度矩阵在线性摄动过程中的奇异性
屈曲分析

问题概述:
一个可靠的产品设计,不仅强度要满足设计要求,而且结构要有足够 的刚度来保证产品性能。现代电子产品(其他产品也一样)已经越来越小, 电子元件之间的空隙非常狭小,在刚度不够的情况下往往导致零部件之间 的干涉。例如,手机从高处坠落,有可能会出现摔坏的情况,这可能是因 为外壳变形过大破坏了内部结构。通常,运行一个静态分析就可以得到结 构在载荷作用下的变形。在某些结构,如承受压应力的部件,在压力载荷 到达一定程度以后会发生于静态分析相比大的多的不可思议的变形,这就 是由于结构已经在这一载荷作用下发生了失稳,这时就需要稳定性分析即 屈曲分析。 实际上结构发生失稳也是由于应力刚度矩阵在影响,应力刚度矩阵可 以加强或减弱结构刚度,这与应力是拉应力还是压应力有关。正如前面计 算出的结果一样,拉应力会使结构的横向刚度增强;结构受压时,会导致 结构的刚度减弱,当压力越来越大时,刚度弱化超出了结构固有的刚度, 结构就表现的很脆弱,位移急剧增大,发生屈曲。
半径i的数值就能使 减小。可见,如果不增加截面面积,尽可能的把材料放在离截 面形心较远处以取得较大的I和i值,就能提高临界应力。
改进措施:在相同截面积下,将杆的结构改为空心杆,截面见图1.其余 各项设置与实心杆相同,计算出的BLF值和实心杆的BLF对比,见图2,失稳 临界载荷因子有明显提高,说明上述分析是正确的。
图1 空心杆截面图
图2 实心杆与空心杆BLF值对比
屈曲分析示例
一端固定一端自由的薄壁圆筒屈曲模态振型
屈曲分析示例
细长圆杆失稳分析及改进
细长圆杆如下图,直径15mm,长200mm,一端固定,一端自由,且受 到100N的压力作用,进行失稳分析并改进。
圆杆三维模型
各阶失稳临界载荷因子(BLF) 失稳的屈曲模态振型(10阶)
屈曲分析-雷晋芳

稳定系数
图13.查看结果
线性屈曲求得的特征值为5.7,非线性屈曲求得的值为1.89 说明:线性屈曲分析可能会产生非保守结果
操作流程
midas Gen 做屈曲分析的一般流程:
一、线性屈曲——初步评估 二、初始缺陷 三、非线性屈曲
构件计算长度确定——利用屈曲分析获得构件计算长度
GB50017-2003钢结构设计规范 CECS28:90钢管混凝土结构设计与施工规程
① ② ③ ④ ⑤ 自动生成荷载组合, 建立或修改需要转换成非线性荷载工况的荷载组合,如图9 生成非线性荷载工况:主菜单>荷载>建立荷载工况>使用荷载组合建立荷载工况, 如图10 查看在该工况下线弹性分析位移最大的点,做非线性分析的控制节点,如图11(预估) 设定非线性控制数据:主菜单>分析>分析控制>非线性分析控制,如图12 查看荷载-位移曲线:结果>时程>阶段/步骤时程图表,如图13
检查斯图姆序列:勾选该项可检查任何丢失的屈服荷载 系数,若存在,会在信息窗口给出报错提示。
常量、变量: 常量+屈曲荷载系数x变量=屈曲荷载值
图5.屈曲分析控制数据
操作流程
3、定义屈曲分析控制数据
恒载与活载: 结构可能会同时受到恒载与活载 的作用,而得到的屈曲荷载系数将会对 所有荷载进行缩放——不分恒载与活载。 这时候需要将二者区分开来,毕竟在多 次试算过程中,恒载的作用效应是不应 该变化的。这时的操作方法就是:恒载 作为不变荷载,活载作为可变荷载,结 果,屈曲荷载就等于“恒载+修正之后 的活载”
2、非线性分析不收敛,如何调整?
加载步骤数:可适当提高,但不是关键因素 子步骤迭代次数:该值对收敛影响较大,增加迭代次数,一般位移范数会减小 收敛条件:一般为0.001,适当增大,可改善收敛性,但不建议修改 非线性分析时,子步骤内迭代次数达到了设定值,各位移范数仍大于收敛容差, 则程序认为非线性未收敛
弹性力学_第十二章板弯曲(ding

• Plate thickness--the distance between the two plate faces. It is denoted by δ.
• 板的厚度:两个板面之间的距离( δ )
Plate middle plane(中面)—— The plane parallel to the faces of the plate and bisecting the thickness δ is called the middle plane of the plate.
• Coordinate system(坐标系)-x and y are in the middle plane and z axis is perpendicular to the middle plane . The system is a right hand system.
x
y z
19
薄板的小挠度弯曲理论
由于 xz=0, yz=0和 z=0,可见中面的法线在薄板弯曲时
保持不伸缩,并且称为弹性曲面的法线。
26
薄板假设2:应力分量 xz, yz和 z远远小于其余三个应力
分量,因而是次要的,它们所引起的形变可以不计。 注意:这3个次要应力分量本身是维持平衡所必须的,不 能不计。
薄板小挠度问题中的物理方程与薄板平面应力问题的 物理方程相同(但两种问题中应力和形变分量沿厚度方 向的分布是不同的)。
第十二章 薄板弯曲问题。经典解答。
§9.1 有关概念和基本假定 §9.2 弹性曲面的微分方程 §9.3 薄板横截面上的内力 §9.4 边界条件 扭矩的等效剪力 §9.5 四边简支矩形薄板的重三角级数解 §9.6 矩形薄板的单三角级数解 §9.7 矩形薄板的差分解(*) §9.8 圆形薄板的弯曲
应用扰动广义微分求积法的复合材料层合板剪切屈曲分析与优化

l
i
d
i
t
ft
hepe
r
t
urbed GDQ me
t
hod.Thes
t
a
ck
i
ngs
equenc
eo
f
yo
c
ompo
s
i
t
el
ami
na
t
e
si
nc
l
ud
i
ngshe
a
rl
oadswa
sop
t
imi
z
edcomb
i
nedpe
r
t
urba
t
i
onGDQ me
t
hodwi
t
ht
he
d
i
r
e
c
ts
a
rf
o
r
c
e
s.Ther
e
su
l
t
sp
r
ov
i
dear
e
f
e
r
enc
ef
o
rt
heshe
a
rbuck
l
i
ng
r
f
o
rmanc
ede
s
i
fc
ompo
s
i
t
el
ami
na
t
e
s.
pe
gno
:
;
Ke
w
o
r
d
s
s
h
e
a
rb
u
c
k
层板分层屈曲后屈曲问题中的Mindlin模型_杨刚

收稿日期 1998-09-02第一作者:39岁,男,副教授,沈阳建筑工程学院土木工程系,沈阳,110015.层板分层屈曲后屈曲问题中的Mindlin 模型杨 刚 张凤鹏 黄宝宗摘 要 针对复合材料层合板壳结构分层问题,提出二维分层M indlin 模型,给出分层前缘处的位移连续条件及转换关系,并通过算例证明该分析模型的可靠性和正确性.关键词 复合材料;分层;屈曲;后屈曲中图法分类号 O 343.9分层模型在分层研究中一直受到关注,它简化得是否合理将直接影响分析结果的实际应用.目前无论用解析法还是数值法大多采用的是薄膜模型 1~3.它把分层(薄膜)视为无限厚物体的一部分,或将分层固支在一定厚度的母体上,不考虑层板(母体)总体行为的影响.由于 薄膜 近似模型简单,子层边界条件容易处理,因此仍被广泛采用.薄膜模型中分层厚度远小于母体厚度,因此,认为子层屈曲而层板总体不弯曲是可行的,并能得到相对满意的结果.但对于分层厚度与母体厚度可比的层板,在分层屈曲的同时很可能发生整体屈曲,这也是工程中常见的分层屈曲形式,如果再用薄膜模型进行简化,将会产生较大的误差.图1 分层模型一个较好的分析模型是把分层层板考虑成几块板或壳的组合体(见图1).它由三部分组成:一是上面的分层部分,二是底层部分,三是未分层部分.通常分层厚度与层板厚度可比.由于分层边缘位移场不象薄膜模型限制的那么强,故在局部屈曲的同时,可伴随总体屈曲.采用这一模型能得到较合理的分层屈曲和后屈曲解,分析分层的出现对总体屈曲载荷的影响,分析局部屈曲和总体屈曲的相互作用.这一模型在分层前缘处的变形比较复杂,考虑了横向剪切变形影响.而剪切变形理论有很多种,其中使用比较广泛且又比较简单的是M indlin 假设.这种二维分层Mindlin 模型多用于数值法,如有限元法.本文首先简单阐述二维分层Mindlin 模型的基本思想及相关问题,然后给出一些算例,证明这种斜法线假设在一定条件下可以得到较精确的结果.1 基本假设和基本方程考虑一层板由各种铺设方向的单向纤维增强材料构成.板内含有一任意形状的分层,将层1998年10月第14卷第4期沈 阳 建 筑 工 程 学 院 学 报Journal of Sheny ang Arch.and Civ.Eng.I nstOct 1998V ol 14,No 4图2 M indlin 板分层前缘区示意图板分为三部分:分层区 1,底层区 2和未分层部分为 0区(如图1所示).假设每个区域均为非线性变形的Mindlin 板壳,在分层前缘处3块板壳位移和转角连续,这意味着裂尖处变形前中面法线在变形后仍保持直线(如图2所示).这一M indlin 模型与只考虑分层区非线性的薄层(Thin -film)模型比较,包括了整个板壳的非线性影响,因而具有更好的精度.按二维分层M indlin 模型,区域 1, 2和 0是3个M indlin 板壳,在各自中面坐标系下的基本关系见文献 4 .在交界面处(见图2)需满足下列连续条件.假设在A(i)点(3块板中面与交界面的交点)的位移和转角为u (i ),v (i ),w (i ), (i )x , (i )y ,而 (i )z 为横向剪切应变,则(i)=w ,(i )- (i )z ( =x ,y ; i =0,1,2)(1)由连续条件,A (1)和A (2)点的位移和转角与A (0)点有以下关系,(1) = (2) = (0)u (i )=u (0)- (0)x z i =u (0)-w ,(0)x z i + (0)x z z iv(i )=v(0)-w ,(0)y z i + (0)y z z iw (i)=w (2)=w (0) i =1,ch2(2)2 分层前缘位移转换关系本文采用有限元法分析层板屈曲和后屈曲问题,基本理论可见文献 3,4 ,这里主要根据文献 4 中给出的连续条件得到分层前缘处的位移转换关系.设 (0), (1)和 (2)分别代表结点A(0),A (1)和A (2)点的广义位移矢量,即(0)=[u (0),v (0),w (0),w ,(0)x ,w ,(0)y , (0)x z , (0)yz , (1)x z , (1)y z , (2)x z , (2)yz ]T(i)=[u(i ),v(i ),w(i ),w ,(i )x ,w ,(i )y , (i )xz , (i )yz ]T(i =1,2)(3)可见,除分层边缘上的点,节点自由度为11外,板内其它点自由度均为7.根据(2)式,它们之间的转换关系可用矩阵表示为(i )=R i 0(0)(4)这里Rz 10000z 1000000-10100-101001001(5)328 沈 阳 建 筑 工 程 学 院 学 报第14卷R 20=100-z 20z 2000000100-z 20z 200000010********0010-100010000010-10001000000000100001((6)在含有分层的层板有限元分析中,使用这一模型时,必须把 (0)作为A(0),A (1),和A (2)点的广义位移矢量,因此当组集总刚时,必须借助矩阵R 10和R 20对 1和 2中含A(1)和A(2)点的单元进行单刚转换.其转换过程是K 0=T T K (i )T(7)式中,K (i)是按 1和 2区坐标系计算得到的单元刚度矩阵,T 为单元刚度转换矩阵:T =T 11T 22T 33)x(8)T ii 按节点是交界上的点,还是 1或 2上的点依次取T ii =I 7 7 T ii =R 10 T ii =R 20(9)经过这种转换处理后,含有这样点的不同壳上的单元刚度矩阵即可进行叠加,进入总刚.3 算例与数值分析例1 含有贯穿分层的条形层板(如图3所示).受单向压缩,参考载荷p =50 103N/m ,板长宽为100mm 20mm,共16层(单向铺设),每层厚0 1125m m.材料为石墨/环氧树脂,E 1=135 4GPa,E 2=9 6GPa, 12=0 31,G 12=G 13=G 23=4 8GPa.分层在1~2层之间,分层尺寸为a .由于对称性取1/2板分析,得到临界载荷P cr 与分层相对尺寸a/l 之间的关系曲线(见图4),可见本文M indlin 模型解比文献 1 的解更接近实验解.图3 贯穿分层模型 图4 临界载荷与a/l 的关系曲线例2 中间含有圆形分层的双向受压层合板(见图5).参考载荷p =103N/m,共10层(0 )10,板长和宽100mm 100mm,厚度10 0 5mm.各层均为各向同性材料,E =6 5GPa,G =2 5GPa, =0 3.分层在1~2层之间,分层半径为r =30mm,考虑对称性,取1/4板分329第4期杨 刚等:层板分层屈曲后屈曲问题中的M indlin 模型析.按限制总体弯曲(在 0, 2区z =0,面内w =0)和不限制总体弯曲,计算了载荷与分层中点挠度关系曲线(见图6),载荷与能量释放率(45 分层边缘上的能量释放率)关系曲线见图7,能量释放率沿分层边缘分布曲线见图8.图5 圆形分层模型 图6 载荷与分层中点挠度曲线图7 载荷与能量释放率曲线 图8 能量释放率沿分层边缘分布曲线从图6,7中可以看到,用M indlin 模型限制总体弯曲的解与文献 2 的解析解非常吻合;而不考虑总体弯曲影响时,两个解差别有所增加,但由于分层较薄,这种差别仍然很小.例3 中部含有圆形分层的单向受压层合板(见图5).参考载荷p =106N/m,共8层( 45 /0 /90 )s ,板长和宽100mm 100m m,第一层厚0 4mm.其余各层厚0.5143mm,分层在1~2层之间,分层半径为r =15mm,材料为石墨/环氧树脂,E 1=134GPa,E 2=10.2GPa, 12=0.3,G 12=G 13=5.52GPa,G 23=3.43GPa,限制总体弯曲.得到分层中点挠度与受载边位移(应变)之间的关系曲线见图9;能量释放率沿分层边缘的分布曲线见图10.图9 分层中点挠度与受载边应变关系曲线 图10 能量释放率沿分层边缘分布曲线可见,本文M indlin 模型解与文献 5 的三维解很接近.330 沈 阳 建 筑 工 程 学 院 学 报第14卷参 考 文 献1 W oo -M in Kyoung and Chun -Gon Kin.Delamination Buckling and Gr owth of Composite L aminated Plates w ithT r ansverse Shear Deformation.Journal of Composite M aterials,1995,29(15):2047~20682 Evans A G and Hutchinson J W.On the M echanics of Delamination and Spalling in Compressed Films.Int.J.Solids Structures.1984,20(5):455~4663 Srivatsa K S,V idy ashankar B R,Krishna M urty A V and V ijaykumar K.Buckling of Laminated platesContaining Delaminatio puters &structures,1993,48(5):907~9124 Huang B Z,Shenoy V B,Atluri S N.Aquas-i conforming triang ular lamminated composite shell element based ona refined first -order put,Mech.,1994,13(4):295~3145 Jo hn D.Whitcomb.T hree -Dimensional Analysis of a Postbuckling Embedded Delamination.Journal of CompositeM aterials,1989,23(10):862~889Min dlin model in delamination bu ckling an d postbu ckling of lamin ate Yan g G an g (Dept.ofCiv.Eng.,Shenyang Arch.and Civ.Eng.Inst.,Sheny ang ,110015,China)Zhang Fengpeng,Huang Baozhong (Northeast university,Shenyang,110006,China)Received Sep. 2 1998Abstract A M indlin model of 2-D delamination is proposed for delamination on laminatedcomposite plat/shell structure;Continuous condition and conversion relation of foredisplacement of delamination are show n;The model is proven to be reliable and correct by ex amples.Keywords Composite materials;delamination;buckling;postbuckling331第4期杨 刚等:层板分层屈曲后屈曲问题中的M indlin 模型。
12-屈曲分析

Fapp
用载荷控制能达到 Fapp吗?
u
3、非线性特征值屈曲
载荷控制:
Training Manual
Advanced Contact & Fasteners
• 使用 Newton-Raphson 载荷控制的困难是求解不能通过不稳定点。 在不稳定点 (Fcr), 切线刚度矩阵 KT 是奇异的,使用载荷控制, Newton-Raphson 法不收敛。然而, 该类型的分析对描述结构的前 屈曲 行为是有用的。
分叉点 极限点
理想载荷路径 有缺陷结构的载荷路径
理想静态行 为
实际动态响应
前屈曲
后屈曲
u
3、非线性特征值屈曲
Training Manual
Advanced Contact & Fasteners
有几种分析技术用于计算结构的非线性 静力变形响应,这些技术包括: -载荷控制 -位移控制 -弧长法 载荷控制: 如下图所示, 考虑浅拱的快速通过分析,当以增量载荷 (F) 求解该问题时, 求解采用载荷控制来完成。
2、线性特征值屈曲
Training Manual
Advanced Contact & Fasteners
• 然而, 缺陷和非线性行为阻止大多数实际结构达到理想的弹性屈曲 强度,特征值屈曲一般产生非保守 解, 使用时应谨慎。 F
极限载荷 分叉点 理想载荷路径 有缺陷结构的载荷路径
前屈曲
u
尽管特征值屈曲一般产生非保守的结果, 线性屈曲分析仍有两个优点: -相对不费时(快捷)的分析。 -为了提供更真实的结果, 屈曲模态形状可用作非线性屈曲分析的初始几 何缺陷
2、线性特征值屈曲
Training Manual
ANSYS Workbench 17·0有限元分析:第13章-特征值屈曲分析

第13章 特征值屈曲分析
屈曲分析主要用于研究结构在特定载荷下的稳定性以及确定结构失稳的临界载荷,屈曲分析包括线性屈曲分析和非线性屈曲分析。
线性屈曲分析可以考虑固定的预载荷,也可使用惯性释放;非线性屈曲分析包括几何非线性失稳分析、弹塑性失稳分析、非线性后屈★ 了解线性屈曲分析。
13.1 屈曲分析概述
特征值屈曲分析(Eigenvolue Buckling)是以特征值为研究对象的,特征值或线性屈曲分析预测的是理想线性结构的理论屈曲强度(分歧点),特征值方程决定了结构的分歧点。
然而,非理想和非线性行为阻止了许多真实的结构达到它们理论上的弹性屈曲强度。
线性屈曲通常产生非保守的结果,应当谨慎使用。
尽管屈曲分析是非保守的,但是也有许多优点。
屈曲分析比非线性屈曲分析计算省时,并且应当作第一步计算来评估临界载荷(屈曲开始时的载荷)。
通过线性屈曲分析可以预知结构的屈曲模型形状,结构可能发生屈曲的方法可以作为设计中的向导。
13.1.1 关于欧拉屈曲
结构的丧失稳定性称为(结构)屈曲或欧拉屈曲。
L.Euler
从一端固定、另一端自由的受压理想柱出发,给出了压杆的临
界载荷。
所谓理想柱,是指起初完全平直而且承受中心压力的
受压杆,如图13-1所示。
设此柱完全是弹性的,且应力不超过比例极限,若轴向外
载荷P小于它的临界值,则此杆将保持直的状态而只承受轴向图13-1 受压杆。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
应用分解刚度法分析金属点阵夹芯板屈曲问题杨丽红;韩笑;吴林志【摘要】Based on Ressiner sandwich panel theory,buckling of truss-core panel with simply-supported boundary conditions was studied under in-plane compression using split rigidity method.For simplicity,truss-core was equivalent to a continuum and the flexural rigidity of the core and the shear stiffness of the face plate were ignored.The buckling critical loads factors were obtained based on the energy principle.The results obtained by using the split rigidity method were compared with those using Hu Haichang method.Buckling critical loads factors corresponding to three kinds of core materials (alloy steel,cast steel and aluminum alloy)were,respectively,determined.The effects of cell configuration of truss-core,the angle of the core rod and the thickness of the face plate on the buckling critical loads were discussed.%基于Ressiner夹层板理论,采用分解刚度法对四边简支点阵夹芯板承受面内压缩荷载时的屈曲问题进行了分析.将点阵夹芯均匀等效为连续体,假定夹芯只提供抗剪切刚度,面板只提供抗弯刚度,基于能量原理得到了屈曲临界载荷系数.将分解刚度法和胡海昌法得到的结果进行了对比.分别针对合金钢、铸钢和铝合金三种材料夹芯进行了临界载荷系数计算.讨论了点阵胞元构型、夹芯杆倾斜角度和面板厚度对屈曲临界载荷的影响规律.【期刊名称】《强度与环境》【年(卷),期】2017(044)003【总页数】7页(P18-24)【关键词】Ressiner夹层板理论;分解刚度法;点阵夹芯板;屈曲【作者】杨丽红;韩笑;吴林志【作者单位】哈尔滨工程大学航天与建筑工程学院,哈尔滨150001;哈尔滨工程大学航天与建筑工程学院,哈尔滨150001;哈尔滨工程大学航天与建筑工程学院,哈尔滨150001;哈尔滨工业大学复合材料与结构研究所,哈尔滨150001【正文语种】中文【中图分类】O3432001年,以普林斯顿大学Evans教授、哈佛大学Hutchinson教授、剑桥大学Ashby教授等为首提出了点阵夹芯结构概念。
点阵夹芯结构由于具有千变万化的微结构和高孔隙率,可实现结构、热控、储能等多功能一体化。
此外,点阵夹芯结构还具有超轻、高比强度、高比刚度、高能量吸收率等优异的力学性能,因此,在航空航天、机械、铁路、汽车、建筑、能源等领域具有广阔的应用前景。
目前,已经提出多种拓扑构型点阵夹芯结构,如金字塔型、四面体型、三维Kagome型等[1],如图1所示。
对于夹芯结构的整体屈曲问题,国内外学者已经进行了大量的研究。
但是对于由离散型点阵芯子构成的夹芯结构,准确地预测其屈曲临界载荷以指导工程实践,是需要我们重点关注的问题。
娄佳[2]基于最小势能原理,对碳纤维复合材料点阵夹芯梁在面内压缩载荷作用下发生整体屈曲时的临界载荷进行了理论预测。
袁武等[3]通过将点阵夹芯等效为均匀连续体,对均匀温度场下点阵夹层板的热屈曲温度进行了求解. Wu Yuan等[4]采用实验和仿真的方法对含局部缺陷的金字塔点阵夹芯板受面内载荷作用时的屈曲问题进行了研究。
M.G. Leea等[5]采用实验的方法研究了WBK夹芯结构在轴向受压时的屈曲问题,并利用理论和仿真的方法对破坏机理进行了分析。
当前分析线弹性夹层板的理论主要有Reissner型理论、Hoff型理论和прусаков-杜庆华理论。
其中,Reissner型理论忽略了面板的抗弯刚度,并认为夹芯只承受抗剪作用;Hoff型理论则考虑面板的抗弯刚度,夹芯仍认为只承受抗剪作用;прусаков-杜庆华理论在考虑了面板抗弯刚度的同时也考虑了夹芯的横向弹性变形的作用[6]。
对于点阵夹芯板而言,由于夹芯的相对密度一般较小,且面板很薄,因而通常采用Reissner型理论进行分析。
分解刚度法是一种分析复合夹层板的有效计算方法,由美国学者Bijlaard[7]提出。
这一方法工程概念直观,结果精确,受到工程设计人员的欢迎。
胡海昌[8]清楚地定义了分解刚度的含义。
李拓等[9]采用分解刚度法对点阵多孔金属夹芯板振动特性进行了分析。
本文基于Reissner夹层板理论,采用分解刚度法分析点阵夹芯板的屈曲问题。
应用Reissner夹层板理论对点阵夹芯板做基本假设1)芯层可等效为各向同性均匀连续介质;2)面板处于薄膜应力状态,对剪切变形没有贡献,只承受弯曲变形;3)芯层只承受横向剪切变形,对弯曲变形没有贡献;4)夹芯横向不可压缩,即σz=0。
夹芯结构芯层的相对密度是指芯层等效密度与材料密度的比值,图2为金字塔点阵夹芯结构芯层的单胞示意图,其相对密度可表示为[10]式中,Ac为夹芯杆横截面面积,hc为芯层高度,lc为夹芯杆长度,ω为夹芯杆与面板的夹角。
金字塔点阵夹芯的等效剪切模量为式(2)中,E为夹芯杆材料的杨氏模量。
于是芯层的等效剪切刚度为对于图3所示点阵夹芯板等效分析模型,基于一阶剪切变形理论,其位移模式为式(4)中,u、v、w分别为夹芯板在x、y、z方向上的位移,φx、φy分别为夹芯中面法线在xz和yz面内的转角,w0为板中面在z方向的位移。
由线弹性物理方程,面板内应力可表述为式(5)中,E为面板材料的杨氏模量,ν为面板材料泊松比。
u+、v+、u-、v-,分别为上下面板内各点在x方向和y方向上的位移,“+”代表上面板,“-”代表下面板。
芯层内应力为由假设式(2)-式(3),点阵夹芯板的内力可表述为式(7)中,为夹芯板的抗弯刚度,其中t为面板厚度。
矩形点阵夹芯板受面内压缩荷载p作用,夹芯板的总势能为式中,A为夹芯板中面面积。
利用分解刚度法,设C→∞、D不变时板的挠度为wb;D→∞、C不变时板的挠度为ws,此时有考虑式(7),可以将夹芯板总势能分解为式(10)中对于四边简支矩形板,假设屈曲模态为式(12)中, Ab、As为任意常数,m、n为整数,a、b为矩形板边长。
将式(12)代入式(11),根据最小势能原理得式(14)中,为板的边长比。
当n=1时,得屈曲临界载荷式(15)中,k为临界载荷系数,其计算式为式(16)中,。
分析金字塔点阵夹芯板,其结构参数如表1。
面板和夹芯材料均为合金钢,材料参数如表2。
下文如无特殊说明,所研究的点阵夹芯板结构参数均如表1所示。
由分解刚度法和胡海昌法[8]得到的临界载荷系数随边长比的变化曲线如图4所示。
由图4可以看出,在相同条件下,采用分解刚度法得到的临界载荷系数比经典的胡海昌法得到的结果偏高;当边长比较小,取m=1时,两者偏差相对较大;随着边长比的增加,两者的偏差逐渐减小,当m取3时,两者结果基本相同。
分解刚度法概念直观,便于应用能量法进行夹层结构力学问题分析,计算结果精度高,因此采用分解刚度法分析点阵夹芯结构屈曲问题有重要意义。
Guoqi Zhang[11]通过实验和仿真的方法分析金字塔点阵夹芯结构,发现夹芯选用铝合金材料时其吸能效果比选用复合材料更好。
但由于铝合金材料刚度较低,因此含铝合金夹芯的点阵结构更容易出现屈曲问题。
本文针对夹芯材料为合金钢、铸钢和铝合金三种情况的金字塔点阵夹芯结构进行分析,面板仍选用合金钢材料,各种材料的材料参数如表2所示。
由式(16)得各种材料对应的临界载荷系数随边长比的变化曲线如图5所示。
由图5可以看出,材料对屈曲临界载荷系数影响较大。
此外,对于铝合金夹芯结构,不同边长比对应的临界载荷系数变化不大。
点阵夹芯结构具有可设计性,当胞元构型、芯杆倾斜角度、面板厚度t等结构参数发生变化时,夹芯板的力学性能均会随着发生变化,因此有必要进一步分析结构参数对点阵夹芯结构屈曲问题的影响。
4.1 胞元构型的影响本文在芯层相对密度相等的情况下,分别讨论金字塔点阵夹芯、四面体点阵夹芯和3D-Kagome点阵夹芯三种常用胞元构型的夹芯板屈曲问题。
四面体夹芯代表性胞元如图6所示,3D- Kagome夹芯代表性胞元如图7所示。
四面体点阵夹芯的相对密度和等效剪切刚度分别为3D- Kagome夹芯的相对密度和等效剪切刚度分别为令各种胞元构型夹芯板的芯层高度相等,取金字塔点阵夹芯的夹芯杆倾斜角度ω=45°,在相对密度相等的情况下,四面体点阵夹芯的芯杆倾角ω=55°,3D-Kagome点阵夹芯的芯杆倾角ω=60°,胞元尺寸B=1.1068·hc。
取各种胞元构型夹芯板结构参数为:面板厚度t=1mm,芯杆直径d=2mm,芯层高度hc=15mm,,板边长b=300mm,代入式(18)和式(20),得各胞元构型的等效刚度,于是由式(16)可以计算临界载荷系数。
各种胞元构型夹芯板对应的临界载荷系数随边长比的变化曲线如图8所示。
由图8可以看出,在相同的相对密度下,金字塔构型的临界载荷系数高于其他两种构型,随边长比的增加三种胞元构型对应的临界载荷系数差距趋于稳定。
4.2 夹芯杆倾斜角度的影响对于点阵夹芯板,夹芯杆倾斜角度ω直接影响芯层的相对密度和等效剪切刚度C,从而间接影响夹芯板的屈曲性能。
通常情况下,ω增大,芯层的相对密度和等效剪切刚度C也随之增加,进而导致板的屈曲临界载荷增大。
对于金字塔点阵夹芯方板,α=1,m=1,临界载荷系数随ω的变化情况如图9所示。
由图9可以看出,临界载荷系数随夹芯杆的倾斜角度近似按线性规律增加,且各种材料的变化率基本相同。
4.3 面板厚度的影响由于点阵夹芯板的夹芯很软,剪切刚度很小,因此如果增大面板厚度会增大结构的稳定性。
取正方形金字塔点阵夹芯板边长为500mm,芯杆直径d=2mm,芯层厚度t=15mm,α=1,m=1,临界载荷随面板厚度的变化情况如图10所示。