基于渐进变分法的蛋盒型结构等效刚度分析及优化

工程科学学报，第39卷，第9期：1386-1395,2017年9月 Chinese Journal of Engineering,Vol.39,No.9:1386-1395,September 2017 D0l:10.13374/j.issn2095-9389.2017.09.012;htp:/journals..usth.edu.cn 基于渐进变分法的蛋盒型结构等效刚度分析及优化 臧 勇，王远，秦勤区，管奔 北京科技大学机械工程学院，北京100083 ☒通信作者，E-mail:qinqin@me.usth.cdu.cn 摘要对蛋盒型结构的等效刚度特性进行了分析并实现了蛋盒型结构参数的优化设计.首先以蛋盒型结构的单胞为研究 对象，基于渐进变分法，得到了蛋盒型结构等效刚度特性的数值计算方法.随后用该方法计算蛋盒型结构不同参数情况下的 等效刚度特性，并以结构参数为自变量，等效刚度特性为因变量进行拟合.最后应用拟合公式在限定泊松比或等效刚度的情 况下，分别以最大化结构的屈曲载荷和最大化单位质量吸能能力的优化为例，对蛋盒型结构参数进行了量纲为一的优化设 计.计算结果表明：蛋盒型结构拉伸刚度降低，弯曲刚度升高：蛋盒型结构的刚度特性与结构参数之间呈现非线性的特点，结 构表现出负泊松比的特性：在给定优化目标和限定条件时应用拟合公式可以快速实现蛋盒型结构参数的主动优化设计. 关键词蛋盒型结构：渐进变分法：等效刚度：结构优化 分类号TB383 Equivalent stiffness and optimized design of egg-box structure based on variational as- ymptotic method ZANG Yong,WANG Yuan,QIN Qin,GUAN Ben School of Mechanical Engineering,University of Science and Technology Beijing,Beijing 100083,China Corresponding author,E-mail:gingin@me.ustb.edu.cn ABSTRACT The equivalent stiffness of the egg-box structure was investigated to achieve structural optimization of the structure pa- rameters.A unit cell of the egg-box structure was studied using the variational asymptotic method to obtain a numerical calculation method for equivalent stiffness of the whole egg-box structure.This numerical calculation method was then used to analyze the equiva- lent stiffness of the egg-box structure with different structural parameters.A formula for equivalent stiffness characteristics was then fit- ted,and the structural parameters were treated as independent variables.Finally,the fitted formula was used to optimize the dimen- sionless structural parameters to obtain maximum buckling load and maximum energy absorbed per unit weight when Poisson's ratio or equivalent stiffness was restricted.Results show that the tensile stiffness of the egg-box structure is reduced,but the bending stiffness is increased.In addition,there is a nonlinear relationship between stiffness and structural parameters,and the structure exhibits a negative Poisson's ratio.Furthermore,active parameters optimization design of the egg-box structure can be easily achieved with the fitted formula. KEY WORDS egg-box structure;variational asymptotic method;equivalent stiffness;structural optimization 蛋盒型结构是一种具有较高的吸能质量比，较好 是金属薄板、聚合物板或复合材料板等，其结构经过简 的抗弯曲扭转性能，较好的高温稳定性等特点的轻质 单的冷冲压或热冲压即可得到，已经在汽车、船舶行业 结构).同时该种结构易制造易回收，其原材料可以 得到应用 收稿日期：2017-02-22 基金项目：中央高校基本科研业务费资助项目(FF-TP-16-010A3):工信部2016年智能制造综合标准化与新模式应用项目

工程科学学报,第 39 卷,第 9 期:1386鄄鄄1395,2017 年 9 月 Chinese Journal of Engineering, Vol. 39, No. 9: 1386鄄鄄1395, September 2017 DOI: 10. 13374 / j. issn2095鄄鄄9389. 2017. 09. 012; http: / / journals. ustb. edu. cn 基于渐进变分法的蛋盒型结构等效刚度分析及优化 臧 勇, 王 远, 秦 勤苣 , 管 奔 北京科技大学机械工程学院, 北京 100083 苣通信作者, E鄄mail: qinqin@ me. ustb. edu. cn 摘 要 对蛋盒型结构的等效刚度特性进行了分析并实现了蛋盒型结构参数的优化设计. 首先以蛋盒型结构的单胞为研究 对象,基于渐进变分法,得到了蛋盒型结构等效刚度特性的数值计算方法. 随后用该方法计算蛋盒型结构不同参数情况下的 等效刚度特性,并以结构参数为自变量,等效刚度特性为因变量进行拟合. 最后应用拟合公式在限定泊松比或等效刚度的情 况下,分别以最大化结构的屈曲载荷和最大化单位质量吸能能力的优化为例,对蛋盒型结构参数进行了量纲为一的优化设 计. 计算结果表明:蛋盒型结构拉伸刚度降低,弯曲刚度升高;蛋盒型结构的刚度特性与结构参数之间呈现非线性的特点,结 构表现出负泊松比的特性;在给定优化目标和限定条件时应用拟合公式可以快速实现蛋盒型结构参数的主动优化设计. 关键词 蛋盒型结构; 渐进变分法; 等效刚度; 结构优化 分类号 TB383 Equivalent stiffness and optimized design of egg鄄box structure based on variational as鄄 ymptotic method ZANG Yong, WANG Yuan, QIN Qin 苣 , GUAN Ben School of Mechanical Engineering, University of Science and Technology Beijing, Beijing 100083, China 苣Corresponding author, E鄄mail: qinqin@ me. ustb. edu. cn ABSTRACT The equivalent stiffness of the egg鄄box structure was investigated to achieve structural optimization of the structure pa鄄 rameters. A unit cell of the egg鄄box structure was studied using the variational asymptotic method to obtain a numerical calculation method for equivalent stiffness of the whole egg鄄box structure. This numerical calculation method was then used to analyze the equiva鄄 lent stiffness of the egg鄄box structure with different structural parameters. A formula for equivalent stiffness characteristics was then fit鄄 ted, and the structural parameters were treated as independent variables. Finally, the fitted formula was used to optimize the dimen鄄 sionless structural parameters to obtain maximum buckling load and maximum energy absorbed per unit weight when Poisson爷s ratio or equivalent stiffness was restricted. Results show that the tensile stiffness of the egg鄄box structure is reduced, but the bending stiffness is increased. In addition, there is a nonlinear relationship between stiffness and structural parameters, and the structure exhibits a negative Poisson爷s ratio. Furthermore, active parameters optimization design of the egg鄄box structure can be easily achieved with the fitted formula. KEY WORDS egg鄄box structure; variational asymptotic method; equivalent stiffness; structural optimization 收稿日期: 2017鄄鄄02鄄鄄22 基金项目: 中央高校基本科研业务费资助项目(FRF鄄鄄TP鄄鄄16鄄鄄010A3);工信部 2016 年智能制造综合标准化与新模式应用项目 蛋盒型结构是一种具有较高的吸能质量比,较好 的抗弯曲扭转性能,较好的高温稳定性等特点的轻质 结构[1] . 同时该种结构易制造易回收,其原材料可以 是金属薄板、聚合物板或复合材料板等,其结构经过简 单的冷冲压或热冲压即可得到,已经在汽车、船舶行业 得到应用

臧勇等：基于渐进变分法的蛋盒型结构等效刚度分析及优化 ·1387· 目前对蛋盒型结构进行应用的主要方面为碰撞吸 单胞结构上具有形貌复杂的特点，并且宏观结构尺寸 能装置的使用，因此对其的研究主要在失效原理和吸 远大于单胞尺寸，如图1所示. 能能力两个方面.Deshpande与Fleck[]运用塑性铰理 论得到了蛋盒型结构的失效机制图，描述了在不同载 荷情况下的失效形式.Zpan等)研究了蛋盒型结构 在落锤试验过程中的失效形式，并对比了蛋盒型结构 与金属泡沫的吸能能力.Akisanya与Fleck4研究了蛋 盒型结构中单个结构在压缩与剪切载荷共同作用时的 失效形式.Chang等s-)研究了纤维材料制造的蛋盒 10 mm 型结构的准静态压缩过程的变形过程和吸能能力，以 图1蛋盒型结构的形貌[3) 及泡沫填充的蛋盒型结构在准静态压缩过程的吸能能 Fig.I Photograph of aluminum egg-box structure 力.Nowpada等s-)对不同结构参数的铝制蛋盒型结 根据图1，图2所示，蛋盒型结构形貌可近似为正 构在不同边界条件下进行了准静态压缩试验，得到了 弦形状，其公式为z=Hsin(2mx/p)sin(2ry/p)/2,结 结构在压缩过程中的最大载荷能力及结构失效情况 构参数主要包括正弦函数的高度H,周期间距P,假设 等.Sashikumar等[o运用有限元的方法研究了蛋盒型 蛋盒型结构各位置厚度均为：，为便于后续计算需对其 结构在不同边界条件时的吸能能力.已有的研究均表 单胞结构的函数进行量纲为一化 明蛋盒型结构在破坏过程中具有较强的吸能能力，在 1.1蛋盒型结构的应变能 某些情况下可替代泡沫结构或蜂窝结构 蛋盒型结构的单胞尺寸远小于蛋盒型结构的宏观 目前的研究对蛋盒型结构的宏观刚度特性问题研 尺寸，单胞在x方向和y方向周期排布并且在两个方 究较少，实际上，相对于蜂窝、瓦楞、泡沫铝等结构，蛋 向具有相同的间距，可将其假设为一种符合Kirchhoff 盒型结构在保证吸能能力的前提下还具有更高的比刚 假设的各向异性薄板，并且不会产生拉伸弯曲耦合刚 度，更强的侧向载荷承载能力等良好的力学特性，同时 度，根据各向异性薄板结构能量的表达方式，应变能可 有研究发现蛋盒型结构还具有负泊松比等特异性的属 表达为下式， 性]，因此蛋盒型结构适合于替代其他轻质结构作为 Ja 承载构件进行使用.等效刚度特性是一种用来表征轻 [Au 0 A 0 0 01e. 质结构承载能力的参数，可以高效计算得到宏观结构 0 An 0 0 0 0 在不同载荷情况下的位移响应[2-]，为结构的应用奠 An 0 A 0 0 0 定基础.蛋盒型结构由于其单胞形貌复杂，用于瓦楞 2 dS. 结构等效刚度特性的解析方法[4-]对蛋盒型结构并 0 0 D Da Kg 不直接适用，利用传统有限元仿真方法对其宏观结构 2K 0 0 0 0 D2 0 2K到 进行分析时需要大量的网格来表征单胞结构，效率低 0 0 0 De 0 , 精度差.渐进变分法是一种计算周期性材料等效刚度 (1) 特性的通用方法[)，其具有严格的数学公式推导，计 式中：E,∈，E,代表拉伸应变；K,K,K,代表弯曲应 算精度高，可用于蛋盒型结构等效刚度特性的计算. 变：A代表拉伸刚度，A和A,代表拉（压）力与中面拉 因此本文拟利用渐进变分法研究蛋盒型结构等效刚度 伸（压缩）应变之间的刚度系数，A,代表结构在单向受 特性与结构参数之间的关系，并实现结构参数的正向 拉或受压时，产生的纵向力与横向正应变之间的刚度 主动优化设计. 系数，A2代表剪切力与中面剪应变之间的刚度系数；D 基于以上分析，本文首先以蛋盒型结构的单胞结构 代表弯曲刚度，D,和D:代表弯矩与曲率之间的刚度 为研究对象，运用渐进变分法和有限单元法[侧建立了蛋 系数，D定义与A类似，D,代表扭转与扭曲率之间的 盒型结构等效刚度特性的数值计算模型，并使用传统有 刚度系数四，因为两个方向周期相同，所以A=A, 限元方法对计算模型进行了对比验证：然后基于该方法 D1=D3,S为薄板的面积范围. 研究了不同结构参数时的等效刚度特性，以等效刚度特 蛋盒型结构的中面位置矢量可以表达为： 性为因变量，结构参数为自变量进行拟合：最后研究了蛋 r(x,y)=xe +ye2 +zes (2) 盒型结构基于数值计算模型结果的最优参数设计方法. 式中：z为x和y的函数，e,e2,C3分别代表笛卡尔坐标 1蛋盒型结构等效刚度特性的计算 系下x,y,2方向的基矢。令x1=x,x2=y,x=z,任意形 状壳体的拉伸应变中。和弯曲应变p。可以写成山：有 蛋盒型结构在宏观上呈现出各向异性的特点，在 关的公式[]

臧 勇等: 基于渐进变分法的蛋盒型结构等效刚度分析及优化 目前对蛋盒型结构进行应用的主要方面为碰撞吸 能装置的使用,因此对其的研究主要在失效原理和吸 能能力两个方面. Deshpande 与 Fleck [2] 运用塑性铰理 论得到了蛋盒型结构的失效机制图,描述了在不同载 荷情况下的失效形式. Zupan 等[3] 研究了蛋盒型结构 在落锤试验过程中的失效形式,并对比了蛋盒型结构 与金属泡沫的吸能能力. Akisanya 与 Fleck [4]研究了蛋 盒型结构中单个结构在压缩与剪切载荷共同作用时的 失效形式. Chang 等[5鄄鄄7] 研究了纤维材料制造的蛋盒 型结构的准静态压缩过程的变形过程和吸能能力,以 及泡沫填充的蛋盒型结构在准静态压缩过程的吸能能 力. Nowpada 等[8鄄鄄9]对不同结构参数的铝制蛋盒型结 构在不同边界条件下进行了准静态压缩试验,得到了 结构在压缩过程中的最大载荷能力及结构失效情况 等. Sashikumar 等[10]运用有限元的方法研究了蛋盒型 结构在不同边界条件时的吸能能力. 已有的研究均表 明蛋盒型结构在破坏过程中具有较强的吸能能力,在 某些情况下可替代泡沫结构或蜂窝结构. 目前的研究对蛋盒型结构的宏观刚度特性问题研 究较少,实际上,相对于蜂窝、瓦楞、泡沫铝等结构,蛋 盒型结构在保证吸能能力的前提下还具有更高的比刚 度,更强的侧向载荷承载能力等良好的力学特性,同时 有研究发现蛋盒型结构还具有负泊松比等特异性的属 性[11] ,因此蛋盒型结构适合于替代其他轻质结构作为 承载构件进行使用. 等效刚度特性是一种用来表征轻 质结构承载能力的参数,可以高效计算得到宏观结构 在不同载荷情况下的位移响应[12鄄鄄13] ,为结构的应用奠 定基础. 蛋盒型结构由于其单胞形貌复杂,用于瓦楞 结构等效刚度特性的解析方法[14鄄鄄16] 对蛋盒型结构并 不直接适用,利用传统有限元仿真方法对其宏观结构 进行分析时需要大量的网格来表征单胞结构,效率低 精度差. 渐进变分法是一种计算周期性材料等效刚度 特性的通用方法[17] ,其具有严格的数学公式推导,计 算精度高,可用于蛋盒型结构等效刚度特性的计算. 因此本文拟利用渐进变分法研究蛋盒型结构等效刚度 特性与结构参数之间的关系,并实现结构参数的正向 主动优化设计. 基于以上分析,本文首先以蛋盒型结构的单胞结构 为研究对象,运用渐进变分法和有限单元法[18]建立了蛋 盒型结构等效刚度特性的数值计算模型,并使用传统有 限元方法对计算模型进行了对比验证;然后基于该方法 研究了不同结构参数时的等效刚度特性,以等效刚度特 性为因变量,结构参数为自变量进行拟合;最后研究了蛋 盒型结构基于数值计算模型结果的最优参数设计方法. 1 蛋盒型结构等效刚度特性的计算 蛋盒型结构在宏观上呈现出各向异性的特点,在 单胞结构上具有形貌复杂的特点,并且宏观结构尺寸 远大于单胞尺寸,如图 1 所示. 图 1 蛋盒型结构的形貌[3] Fig. 1 Photograph of aluminum egg鄄box structure 根据图 1,图 2 所示,蛋盒型结构形貌可近似为正 弦形状,其公式为 z = Hsin (2仔x / p) sin (2仔y / p) / 2,结 构参数主要包括正弦函数的高度 H,周期间距 p,假设 蛋盒型结构各位置厚度均为 t,为便于后续计算需对其 单胞结构的函数进行量纲为一化. 1郾 1 蛋盒型结构的应变能 蛋盒型结构的单胞尺寸远小于蛋盒型结构的宏观 尺寸,单胞在 x 方向和 y 方向周期排布并且在两个方 向具有相同的间距,可将其假设为一种符合 Kirchhoff 假设的各向异性薄板,并且不会产生拉伸弯曲耦合刚 度,根据各向异性薄板结构能量的表达方式,应变能可 表达为下式, J = 1 2 蓦s 缀x 2缀xy 缀y 资x 2资xy 资 ì î í ï ï ï ï ï ï ï ï ü þ ý ï ï ï ï ï ï ï ï y T A11 0 A13 0 0 0 0 A22 0 0 0 0 A13 0 A33 0 0 0 0 0 0 D11 0 D13 0 0 0 0 D22 0 0 0 0 D13 0 D é ë ê ê ê ê ê ê ê ê ù û ú ú ú ú ú ú ú ú 33 缀x 2缀xy 缀y 资x 2资xy 资 ì î í ï ï ï ï ï ï ï ï ü þ ý ï ï ï ï ï ï ï ï y dS. (1) 式中:缀x,缀xy,缀y 代表拉伸应变;资x,资xy,资y 代表弯曲应 变;A 代表拉伸刚度,A11和 A13代表拉(压)力与中面拉 伸(压缩)应变之间的刚度系数,A13代表结构在单向受 拉或受压时,产生的纵向力与横向正应变之间的刚度 系数,A22代表剪切力与中面剪应变之间的刚度系数;D 代表弯曲刚度,D11和 D13代表弯矩与曲率之间的刚度 系数,D13定义与 A13类似,D22代表扭转与扭曲率之间的 刚度系数[22] ,因为两个方向周期相同,所以 A11 = A33 , D11 = D33 ,S 为薄板的面积范围. 蛋盒型结构的中面位置矢量可以表达为: r(x,y) = x^e1 + y^e2 + z ^e3 . (2) 式中:z 为 x 和 y 的函数,^e1 ,^e2 ,^e3 分别代表笛卡尔坐标 系下 x,y,z 方向的基矢. 令 x1 = x,x2 = y,x3 = z,任意形 状壳体的拉伸应变 准琢茁和弯曲应变 籽琢茁可以写成 ui 有 关的公式[17] . ·1387·

·1388· 工程科学学报，第39卷，第9期 图2蛋盒型的形貌及其单胞坐标量纲为一化.(a)结构形貌：(b)A-A截面：(c)B-B截面 Fig.2 Geometry of egg-box structure:(a)entire geometry;(b)section A-A (c):section B-B Ni dui, 2中p=。axg ((a4eP+a0s)小+ Ea（(ape)2+ata"p.Pw)）ds. 24(1+u)T1-u nu+9(en%+e0b2). 2 ax (5) (3) 1.2基于渐进变分法的蛋盒型单胞结构能量表达 式中：4：表示结构在全局坐标系下的位移分量，「代 以单胞为研究对象，建立一个量纲为一的局部坐 表克里斯托弗尔符号，「。代表曲面的切向量分量，n:代 标系X。,X。的范围为-1/2到1/2，根据一般周期结构 表曲面的法向量分量，e.代表Levi-Civita张量，0代表 材料的渐进变分法的表述16)，可以将双方向周期排 浆法向矢景的底转角9=26（之-2)，以= 布板壳结构的位移表述为以下方式： 2（ 山1(X1,X2,x1,x2,x)=(x1,x2,x3)+ abm,a吧为曲面的度量张量，bn为曲面形貌的二次形 p吨，(X1,X2,x,x2x3)-x(X1,X2)D.1 式，a为‖a‖的行列式.上述公式及后续公式中a, 山2（X,X2,x1,x2,x3）=2(x1,x2,x3）+ B,y取值为1,2，i取值为1,2,3. p2(X1,X2,x,x2,)-（X,X)2 则蛋盒型结构的应变能密度为[6： 山3（X,X2,x1,2,x3）=3（x1,x2,x3）+ b=2+a*a*0ns）+ Et p3(X1,X2,x1,2,x3）. (6) 其中的：代表全局坐标系下蛋盒型结构中面的位移， ain(ap产+agp） (4) 山：代表局部坐标系下单胞结构的位移，.。代表对 X。求导.则结合式(1)、式(5)和式(6)，单胞的应变 式中，E为基础材料的杨氏模量，4为材料的泊松比 能可写为公式(7)所示的形式 蛋盒型结构单胞上的应变能可以表达为： 1=R)a, (7) 式中：K为材料的刚度，由式(5)展开可得：R=

工程科学学报,第 39 卷,第 9 期 图 2 蛋盒型的形貌及其单胞坐标量纲为一化. (a)结构形貌;(b)A鄄鄄A 截面;(c)B鄄鄄B 截面 Fig. 2 Geometry of egg鄄box structure: (a) entire geometry; (b) section A鄄鄄A (c); section B鄄鄄B 2准琢茁 = r i 琢 鄣ui 鄣x茁 + r i 茁 鄣ui 鄣x琢 , 2籽琢茁 = 鄣 鄣x ( 茁 ni 鄣ui 鄣x ) 琢 + 鄣 鄣x ( 琢 ni 鄣ui 鄣x ) 茁 - 2祝 酌 琢茁 ni 鄣ui 鄣x酌 + 兹(e酌琢 b 酌 茁 + e酌茁 b 酌 琢 ). (3) 式中:ui 表示结构在全局坐标系下的位移分量,祝 酌 琢茁代 表克里斯托弗尔符号,r i 琢 代表曲面的切向量分量,ni 代 表曲面的法向量分量,e酌琢代表 Levi鄄鄄Civita 张量,兹 代表 绕法向矢量的旋转角 兹 = 1 2 ( a r i 1 鄣ui 鄣x2 - r i 2 鄣ui 鄣x ) 1 ,b 酌 琢 = a 琢茁 b茁酌 ,a 琢茁为曲面的度量张量,b茁酌为曲面形貌的二次形 式,a 为椰a 琢茁椰的行列式. 上述公式及后续公式中 琢, 茁,酌 取值为 1,2,i 取值为 1,2,3. 则蛋盒型结构的应变能密度为[16] : 椎 = Et 2(1 + 滋 ( ) u 1 - 滋 (a 琢茁准琢茁 ) 2 + a 琢茁 a 酌啄准琢酌准茁啄 ) + Et 3 24(1 + 滋 ( ) 滋 1 - 滋 (a 琢茁 籽琢茁 ) 2 + a 琢茁 a 酌啄 籽琢酌 籽茁啄 ). (4) 式中,E 为基础材料的杨氏模量,滋 为材料的泊松比. 蛋盒型结构单胞上的应变能可以表达为: J = 蓦S 椎 adS = 蓦 ( S Et a 2(1 + 滋 ( ) 滋 1 - 滋 (a 琢茁准琢茁 ) 2 + a 琢茁 a 酌啄准琢酌准茁啄 ) + Et 3 a 24(1 + 滋 ( ) 滋 1 - 滋 (a 琢茁 籽琢茁 ) 2 + a 琢茁 a 酌啄 籽琢酌 籽茁啄 ) ) dS. (5) 1郾 2 基于渐进变分法的蛋盒型单胞结构能量表达 以单胞为研究对象,建立一个量纲为一的局部坐 标系 X琢 ,X琢 的范围为 - 1 / 2 到 1 / 2,根据一般周期结构 材料的渐进变分法的表述[16鄄鄄17] ,可以将双方向周期排 布板壳结构的位移表述为以下方式: u1 (X1 ,X2 ,x1 ,x2 ,x3 ) = v1 (x1 ,x2 ,x3 ) + p鬃1 (X1 ,X2 ,x1 ,x2 ,x3 ) - x3 (X1 ,X2 )v3,1 , u2 (X1 ,X2 ,x1 ,x2 ,x3 ) = v2 (x1 ,x2 ,x3 ) + p鬃2 (X1 ,X2 ,x1 ,x2 ,x3 ) - x3 (X1 ,X2 )v3,2 , u3 (X1 ,X2 ,x1 ,x2 ,x3 ) = v3 (x1 ,x2 ,x3 ) + p鬃3 (X1 ,X2 ,x1 ,x2 ,x3 ). (6) 其中的 vi 代表全局坐标系下蛋盒型结构中面的位移, 鬃i 代表局部坐标系下单胞结构的位移,v3,琢代表 v3 对 X琢 求导. 则结合式(1)、式(5)和式(6),单胞的应变 能可写为公式(7)所示的形式. J = 1 2 乙 1 / 2 -1 / 2 乙 1 / 2 -1 / 2 (R TKR)dX1 dX2 . (7) 式中:K 为 材 料 的 刚 度, 由 式 ( 5 ) 展 开 可 得; R = ·1388·

臧勇等：基于渐进变分法的蛋盒型结构等效刚度分析及优化 ·1389· [2p2pp2]';将R展开为e和 构中面的等效应变，各个值分别为： 亚有关的函数，其中∈代表蛋盒型结构中面的应变，与 ∈，=112ew=1.2+%21,∈，=2.2， 式(1)中相同，亚代表单胞结构的位移，计算∈和亚之 K.=-.1,2K=-3,2,K,=-. (11) 间的关系，即可得到等效刚度特性 100x3 0 07 R=Te+TwΨ (8) 0 10 0 0 其中： 00 1 0 0 3 业=[少中中] (9) r.= (12) I12I3 L I 代表蛋盒型结构任意位置在局部坐标系下的位移 MM,M, M. Ms o e=[∈：2e,∈，K,2KwK,].(10) N N2 Na N Ns 应变e与公式(1)中的应变意义相同，代表蛋盒型结 ax 0 910X1 d a aX, ax aX2 +92aX1 a 0 3X2 920X2 a2 a2 (13) h成+l成+l h品+品+l%品 、02 成+成+ MnaX +MaX 82 +MpaX aX: Ma aX, +MeaX 2 +MsaX,aX: 2 M沉+M元 +MoaX,aX: a2 a2 Nn+元+Nnx 武+Ne成+ a2 式(12)和式(13)中的I,M,N代表与∈和亚相乘的 将式(7)写为式(1)的形式： 项，详细表达式由公式(6)和(3)整理展开后得到，带 人公式(7)和(8)可得到蛋盒型结构应变能的详细表 达式 1.3蛋盒型结构等效刚度特性有限元列式求解 2e(太+k.)ead（a8) 式(8)极其复杂，即使结合周期性边界条件，通过 因此可以得到刚度特性： 解析的方法很难得到∈和亚之间的关系，因此运用最 「A01 =石Km+K (19) 小能量原理和有限单元法得到∈和亚关系的数值解. LO DI 对亚进行离散可得： 其中的5。=-KwKe 亚=Cξ. (14) 蛋盒型结构等效刚度特性的求解过程如图3 式中，C为形函数，专代表离散后各点的亚值.为保证 所示 业连续，选用Hermite插值方法建立9节点Cl型单 1.4方法验证 元，将式(8)，式(14)带入式(7)可得： 选用中面形貌函数为x?=0.005sin(2πx,/0.03) J=2传'Kn5+25Ke+eKe. (15) sin(2mx2/0.03),t=0.0005m,E=71GPa,u=0.33的 铝制蛋盒型结构研究对象.根据图3所示蛋盒型结构 其中： 等效刚度特性的求解过程，首先对其形貌函数进行量 =arokrc)a 纲为一化，根据式(5)可得到其单胞结构应变能的表 e)raxds.. 达式，结合式(6)可将单胞结构的应变能运用渐进变 (16) 分法表述.单胞结构等效刚度特性的计算过程运用了 -ra 有限单元法，将单胞结构离散为400个单元，形函数的 表达形式利用Hermite矩形单元插值函数o],分别计 结合周期性结构的边界条件及最小能量原理[]， 算每个单元的刚度矩阵并得到总刚矩阵表达的应变 对式(15)进行求导可得到： 能，即式(15).该矩阵并不能直接求解，运用最小能量 Kws =-Kv.E. (17) 原理对矩阵进行约束，可以求得矩阵关系的等式，此时

臧 勇等: 基于渐进变分法的蛋盒型结构等效刚度分析及优化 [准 0 11 2准 0 12 准 0 22 籽 0 11 2籽 0 12 籽 0 22 ] T ;将 R 展开为 缀 和 追 有关的函数,其中 缀代表蛋盒型结构中面的应变,与 式(1)中相同,追 代表单胞结构的位移,计算 缀和 追 之 间的关系,即可得到等效刚度特性. R = 祝缀缀 + 祝追追. (8) 其中: 追 = [鬃1 鬃2 鬃3 ] T . (9) 代表蛋盒型结构任意位置在局部坐标系下的位移. 缀 = [缀x 2缀xy 缀y 资x 2资xy 资y] T . (10) 应变 缀 与公式(1) 中的应变意义相同,代表蛋盒型结 构中面的等效应变,各个值分别为: 缀x = 淄1,1 ,2缀xy = 淄1,2 + 淄2,1 ,缀y = 淄2,2 , 资x = - 淄3,11 ,2资xy = - 淄3,12 ,资y = - 淄3,22 . (11) 祝缀 = 1 0 0 x3 0 0 0 1 0 0 x3 0 0 0 1 0 0 x3 I1 I2 I3 I4 I5 I6 M1 M2 M3 M4 M5 M6 N1 N2 N3 N4 N5 N é ë ê ê ê ê ê ê ê ê ù û ú ú ú ú ú ú ú ú 6 . (12) 祝追 = 鄣 鄣X1 0 渍1 鄣 鄣X1 鄣 鄣X2 鄣 鄣X1 渍1 鄣 鄣X2 + 渍2 鄣 鄣X1 0 鄣 鄣X2 渍2 鄣 鄣X2 I71 鄣 鄣X1 + I72 鄣 鄣X2 + I73 鄣 2 鄣X 2 1 I81 鄣 鄣X1 + I82 鄣 鄣X2 + I83 鄣 2 鄣X 2 1 I91 鄣 鄣X1 + I92 鄣 鄣X2 + I93 鄣 2 鄣X 2 1 M71 鄣 鄣X1 + M72 鄣 鄣X2 + M73 鄣 2 鄣X1 鄣X2 M81 鄣 鄣X1 + M82 鄣 鄣X2 + M83 鄣 2 鄣X1 鄣X2 M91 鄣 鄣X1 + M92 鄣 鄣X2 + M93 鄣 2 鄣X1 鄣X2 N71 鄣 鄣X1 + N72 鄣 鄣X2 + N73 鄣 2 鄣X 2 2 N81 鄣 鄣X1 + N82 鄣 鄣X2 + N83 鄣 2 鄣X 2 2 N91 鄣 鄣X1 + N92 鄣 鄣X2 + N93 鄣 2 鄣X é ë ê ê ê ê ê ê ê ê ê ê ê ê ê ê ê ê ù û ú ú ú ú ú ú ú ú ú ú ú ú ú ú ú 2 ú 2 . (13) 式(12)和式(13) 中的 I,M,N 代表与 缀 和 追 相乘的 项,详细表达式由公式(6)和(3)整理展开后得到,带 入公式(7)和(8)可得到蛋盒型结构应变能的详细表 达式. 1郾 3 蛋盒型结构等效刚度特性有限元列式求解 式(8)极其复杂,即使结合周期性边界条件,通过 解析的方法很难得到 缀和 追 之间的关系,因此运用最 小能量原理和有限单元法得到 缀 和 追 关系的数值解. 对 追 进行离散可得: 追 = C孜. (14) 式中,C 为形函数,孜 代表离散后各点的 追 值. 为保证 追 连续,选用 Hermite 插值方法建立 9 节点 C1 型单 元,将式(8),式(14)带入式(7)可得: J = 1 2 (孜 TK追追 孜 + 2孜 TK追着 缀 + 缀 TK着着 缀). (15) 其中: K追追 = 乙 1 / 2 -1 / 2 乙 1 / 2 -1 / 2 ((祝追 C) TK(祝追 C))dX1 dX2 , K追着 = 乙 1 / 2 -1 / 2 乙 1 / 2 -1 / 2 ((祝追 C) TK祝缀)dX1 dX2 , K着着 = 乙 1 / 2 -1 / 2 乙 1 / 2 -1 / 2 (祝 T 缀 K祝缀)dX1 dX2 . (16) 结合周期性结构的边界条件及最小能量原理[19] , 对式(15)进行求导可得到: K追追孜 = - K追着 缀. (17) 将式(7)写为式(1)的形式: J = 1 2 乙 1 / 2 -1 / 2 乙 1 / 2 -1 / 2 缀 T A 0 0 [ ] D æ è ç ö ø 缀÷dX1 dX2 = 1 2 乙 1 / 2 -1 / 2 乙 1 / 2 -1 / 2 (缀 T (孜 T 0K追着 + K着着 )缀)dX1 dX2 . (18) 因此可以得到刚度特性: A 0 0 [ ] D = 孜 T 0K追着 + K着着 . (19) 其中的孜0 = - K - 1 追追K追着 . 蛋盒型结构等效刚度特性的求解过程如图 3 所示. 1郾 4 方法验证 选用中面形貌函数为 x3 = 0郾 005sin (2仔x1 / 0郾 03) sin (2仔x2 / 0郾 03),t = 0郾 0005 m,E = 71 GPa,滋 = 0郾 33 的 铝制蛋盒型结构研究对象. 根据图 3 所示蛋盒型结构 等效刚度特性的求解过程,首先对其形貌函数进行量 纲为一化,根据式(5) 可得到其单胞结构应变能的表 达式,结合式(6) 可将单胞结构的应变能运用渐进变 分法表述. 单胞结构等效刚度特性的计算过程运用了 有限单元法,将单胞结构离散为 400 个单元,形函数的 表达形式利用 Hermite 矩形单元插值函数[20] ,分别计 算每个单元的刚度矩阵并得到总刚矩阵表达的应变 能,即式(15). 该矩阵并不能直接求解,运用最小能量 原理对矩阵进行约束,可以求得矩阵关系的等式,此时 ·1389·

·1390· 工程科学学报，第39卷，第9期 蛋盒结构基本参数 应变能的计算 等效刚度特性有限元求解 周期边 蛋盒 蛋盒 单胞 界条件 型 计算 结构 离散 等 单胞 应变 刚度 结构 能渐 量纲 应变 进变 分法 胞结构离散化 总刚矩阵的计算 求 属性 一化 表述 最小能 量原理 图3 蛋盒型结构等效刚度特性的求解过程 Fig.3 Process used to solve egg-box equivalent stiffness 该等式的解不唯一，即式(17)，结合周期性边界条件 特性的影响规律，也可对计算结果进行拟合以实现蛋 对总刚矩阵进行处理】，即可求得唯一解，即公式 盒型结构参数的优化设计 (19). 2.1不同结构参数情况下的刚度特性 假设蛋盒型结构在x方向和y方向均具有10个 假设基础材料为杨氏模量E=71GPa,泊松比μ= 周期，三维模型建立详细的蛋盒型结构细节有限元模 0.33,厚度为t的铝薄板，改变H和p,得到不同参数时 型，等效模型直接将得到的等效刚度特性赋予壳单元， 的等效刚度特性，并除以对厚度1平板的刚度特性，得 结构四面固定支撑，受到P=0.01MPa的压强，建立1/ 到蛋盒型结构刚度特性量纲为一的变化规律，定义量 4模型并对比等效模型与三维模型对称边上的位移如 纲为一的等效拉伸刚度为A,量纲为一的等效弯曲刚 图4所示. 度为D.常见的蛋盒型结构均由冲压得到，所以H≤P, 假定H/1在0~50范围内变化，p/L在50~150范围内 一三维模型 一等效模型 变化，当H/t=0时即为初始材料平板，选取其中的不 同的结构参数进行计算，得到其等效刚度特性的变化 1/4模型 规律. -2 载荷P=0.0IMPa 图5~7为蛋盒型结构等效拉伸刚度特性与结构 参数之间的关系，蛋盒型结构任意截面均为正弦曲线， 平面内应变的过程需先将曲线拉直，这就使得蛋盒结 构面内拉伸刚度A较普通平板降低，当结构横向受拉 时纵向产生负的应变，因此A:产生负值，并且该值是 与蛋盒结构的结构参数有关的.当H/≤5时，A基本 50 100 150200 250 300 不受p/L变化的影响.当H/1≥5时，A,随H/L的增大 x方向坐标/mm 而降低，并且随p/1的增大而增大；A,均为负值且与 图4等效模型与三维模型同等载荷下位移响应 H/L和p/L呈现非线性的特点：A,随H/L的增大而降低 Fig.4 Displacement of equivalent model and 3D model under same 但基本不受p/1变化的影响. load 图8~10为蛋盒型结构等效弯曲刚度特性与结构 计算过程蛋盒型结构三维模型需要250000个壳 参数之间的关系，蛋盒型结构各个截面均为正弦曲线， 单元，而等效模型仅需要900个壳单元，等效模型计算 蛋盒型结构的相对高度较普通平板有很大的提升，同 效率大大提高，两种模型的位移趋势相同，且计算误差 时结构符合Kirchhoff假设，因此蛋盒型结构的弯曲刚 在0.5%以内.对比三维模型与等效模型的响应可知 度有所提升.蛋盒型结构板的弯曲刚度与H呈现出高 该计算方法适用于蛋盒型结构等效刚度特性的计算. 度非线性的关系，D,随H/L的增大而增大，随p/L的增 同时三维模型最大位移仅为同等厚度基础材料板材模 大而减小；D.与H/1和p/L呈现出非线性的关系；D 型响应的1/2，有效提升了板材的刚度 随H/1的增大而增大，随p/1的增大而增大，相较于 D,和D变化较大 2蛋盒型结构等效刚度特性分析与优化 根据计算得到的蛋盒型结构在不同参数情况下的 蛋盒型结构板参数主要包括H,P和1，计算不同 等效刚度特性比值，以结构参数为自变量进行拟合，假 结构参数时的等效刚度特性可得到结构参数等效刚度 设x,=pL,xn=H/h.刚度特性的比值A和D的拟合公

工程科学学报,第 39 卷,第 9 期 图 3 蛋盒型结构等效刚度特性的求解过程 Fig. 3 Process used to solve egg鄄box equivalent stiffness 该等式的解不唯一,即式(17),结合周期性边界条件 对总刚矩阵进行处理[21] ,即可求得唯一解,即公式 (19). 假设蛋盒型结构在 x 方向和 y 方向均具有 10 个 周期,三维模型建立详细的蛋盒型结构细节有限元模 型,等效模型直接将得到的等效刚度特性赋予壳单元, 结构四面固定支撑,受到 P = 0郾 01 MPa 的压强,建立1 / 4 模型并对比等效模型与三维模型对称边上的位移如 图 4 所示. 图 4 等效模型与三维模型同等载荷下位移响应 Fig. 4 Displacement of equivalent model and 3D model under same load 计算过程蛋盒型结构三维模型需要 250000 个壳 单元,而等效模型仅需要 900 个壳单元,等效模型计算 效率大大提高,两种模型的位移趋势相同,且计算误差 在 0郾 5% 以内. 对比三维模型与等效模型的响应可知 该计算方法适用于蛋盒型结构等效刚度特性的计算. 同时三维模型最大位移仅为同等厚度基础材料板材模 型响应的 1 / 2,有效提升了板材的刚度. 2 蛋盒型结构等效刚度特性分析与优化 蛋盒型结构板参数主要包括 H,p 和 t,计算不同 结构参数时的等效刚度特性可得到结构参数等效刚度 特性的影响规律,也可对计算结果进行拟合以实现蛋 盒型结构参数的优化设计. 2郾 1 不同结构参数情况下的刚度特性 假设基础材料为杨氏模量 E = 71 GPa,泊松比 滋 = 0郾 33,厚度为 t 的铝薄板,改变 H 和 p,得到不同参数时 的等效刚度特性,并除以对厚度 t 平板的刚度特性,得 到蛋盒型结构刚度特性量纲为一的变化规律,定义量 纲为一的等效拉伸刚度为 A,量纲为一的等效弯曲刚 度为 D. 常见的蛋盒型结构均由冲压得到,所以 H臆p, 假定 H/ t 在 0 ~ 50 范围内变化,p / t 在 50 ~ 150 范围内 变化,当 H/ t = 0 时即为初始材料平板,选取其中的不 同的结构参数进行计算,得到其等效刚度特性的变化 规律. 图 5 ~ 7 为蛋盒型结构等效拉伸刚度特性与结构 参数之间的关系,蛋盒型结构任意截面均为正弦曲线, 平面内应变的过程需先将曲线拉直,这就使得蛋盒结 构面内拉伸刚度 A 较普通平板降低,当结构横向受拉 时纵向产生负的应变,因此A13产生负值,并且该值是 与蛋盒结构的结构参数有关的. 当 H/ t臆5 时,A 基本 不受 p / t 变化的影响. 当 H/ t逸5 时,A11随 H/ t 的增大 而降低,并且随 p / t 的增大而增大;A13 均为负值且与 H/ t 和 p / t 呈现非线性的特点;A22随 H/ t 的增大而降低 但基本不受 p / t 变化的影响. 图 8 ~ 10 为蛋盒型结构等效弯曲刚度特性与结构 参数之间的关系,蛋盒型结构各个截面均为正弦曲线, 蛋盒型结构的相对高度较普通平板有很大的提升,同 时结构符合 Kirchhoff 假设,因此蛋盒型结构的弯曲刚 度有所提升. 蛋盒型结构板的弯曲刚度与 H 呈现出高 度非线性的关系,D11随 H/ t 的增大而增大,随 p / t 的增 大而减小;D13与 H/ t 和 p / t 呈现出非线性的关系;D22 随 H/ t 的增大而增大,随 p / t 的增大而增大,相较于 D11和 D13变化较大. 根据计算得到的蛋盒型结构在不同参数情况下的 等效刚度特性比值,以结构参数为自变量进行拟合,假 设 xp = p / t,xH = H/ t. 刚度特性的比值 A 和 D 的拟合公 ·1390·