CN102096736A - 一种基于渐近变分法的复合材料层合板仿真及优化方法 - Google Patents
一种基于渐近变分法的复合材料层合板仿真及优化方法 Download PDFInfo
- Publication number
- CN102096736A CN102096736A CN2011100320114A CN201110032011A CN102096736A CN 102096736 A CN102096736 A CN 102096736A CN 2011100320114 A CN2011100320114 A CN 2011100320114A CN 201110032011 A CN201110032011 A CN 201110032011A CN 102096736 A CN102096736 A CN 102096736A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- asymptotic
- analysis
- strain
- plate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于渐近变分法的复合材料层合板仿真及优化方法,属于材料力学分析领域。具体包括以下步骤:基于旋转张量分解概念构建由一维广义应变和翘曲函数表示的三维板能量方程;基于渐近变分法将原三维问题分析严格拆分为非线性二维板分析和沿厚度方向的截面分析;利用板固有的小参数将降维模型的近似能量泛函渐近修正到第二阶,并通过平衡方程将近似能量泛函转换为Reissner模型形式以便于实际应用;由得到的全局响应的渐近修正翘曲函数精确重构三维应力/应变/位移场;采用极大化截面分析得到的弯曲和扭转刚度系数的优化策略对复合层合板进行优化。本发明实用性强,通用性高,可显著提高此类问题的解算速度和效率。
Description
技术领域
本发明属于材料力学性能分析领域,尤其是一种能够有效模拟和准确重构复合材料层合板三维应力/应变/变形场和极大化刚度矩阵的弯曲和扭转项为优化策略、层合复合材料铺设倾角为设计变量的优化设计方法。
背景技术
近20年来,先进复合材料结构因其高强度、高模量、可设计性等优点已广泛应用于航天航空、机械、土木等领域。复合材料结构很多是厚度尺寸比其它两个方向尺寸小得多的平板铺层结构,复合材料铺层设计的灵活性为结构的应力计算、变形分析以及强度预测带来额外的复杂性。基于Kirchhoff假设(薄板)、由三维弹性理论推导出的古典层合理论是最简单的复合板分析理论,但若计入板厚,其精度较三维有限元分析下降很多。各国学者利用板厚相对基准面上的变形很小的特点,将板的厚度坐标从偏微分方程的独立变量中去掉,用精确的二维模型代替三维模型以对古典模型不足之处进行改进,但大部分研究都是基于特定的动力学假设(如沿厚度方向的位移分布假设等),无法较好地反映层合板三维效应和铺层之间的相互作用。
随着近年来计算硬件和软件的发展,可使用三维有限元软件ANSYS或NASTRAN对这类复合层合板进行三维计算和分析,如王跃进等以ABAQUS有限元软件为平台,对三维复合层合板渐近损伤进行非线性分析;Apalak的分析模型中对复合材料层合板采用ANSYS软件中的八节点三维层状体单元进行应力求解。但仍有必要对这种结构进行简化分析,原因在于:首先,许多工程问题目前无法用三维有限元分析处理,如实际的转子叶片由200多层非常薄的复合层构成,若使用三维有限元模型,每一铺层至少需要一个单元进行模拟,该叶片模型很容易超过109个自由度。因此,由4个转子叶片构成的旋翼气弹性分析目前无法在任一台计算机上通过三维有限元分析解决;二是在对整个结构进行有限元分析的起始准备工作,包括材料种类的选择、层合和夹层结构的设计以及层合板铺设方式设计,有限元软件的实用性不大,当在层合板结构的不同层的级别上来对复合材料行为(应力、应变分布)进行细节研究时,有限元软件提供的后处理能力尤其有限。另外,尽管三维有限元在连续介质力学框架内提供很好的准确性,但计算太费时,且耗费大量的计算机资源,难以在规定的设计和分析时间内完成预定目标,因此迫切需要一种高效、快速的专业复合材料层合板分析方法和程序,以缩短设计时间,降低计算成本。本专利技术就是在这样的背景下展开的。
发明内容
针对现有技术存在分析效率低、精度差不足,尤其是沿厚度方向的应力分布无法精确预测之不足,本发明的目的是提供一种计算量小,占用计算机资源少,且效率高的基于渐近变分方法的复合材料层合板三维场仿真及优化方法。
本发明采用的技术方案如下:一种基于渐近变分法的复合材料层合板仿真及优化方法,包括如下步骤:
1)基于旋转张量分解方法将复合材料层合板三维连续方程用一维广义应变和翘曲函数表示;
2)将非线性三维复合层合板模型分解为二维截面分析和沿法线方向的一维非线性分析;采用渐近变分算法实现简化模型的求解过程:通过二维截面分析得到本构关系,即刚度矩阵和截面特性和渐近修正翘曲函数;通过沿结线方向的一维非线性分析得到二维板的全局响应性能;
3)基于推导的重构关系,由计算得到的二维板变形、全局旋转张量分量和渐近修正翘曲函数重构三维变形;由一维广义应变和渐近修正翘曲函数重构三维应变场;使用材料本构关系得到三维应力;
4)以弯扭屈曲荷载最大为目标,以不发生屈曲为约束条件,层合复合材料铺设倾角为设计变量的复合梁优化设计。
相比现有技术,本发明具有如下有益效果:
本发明是一种有效模拟和准确重构复合材料层合板三维应力/应变/变形场和极大化刚度矩阵的弯曲和扭转项为优化策略、层合复合材料铺设倾角为设计变量的优化设计方法,借助渐近变分技术提供一种高效、快速、实用性强、通用性强,能够通过渐近修正翘曲函数建立复合材料层合板简化数值模型和三维重构关系的方法。
本发明首先建立用一维广义应变和翘曲函数表示的三维板能量方程数值模式;然后将渐近变分引入到复合材料层合板简化数值建模问题研究中,将原三维问题分析严格拆分为二维截面分析(等效单层板模型)和沿法线方向的一维非线性分析。本发明通过二维截面分析得到沿法线一维非线性分析所需的本构关系(刚度矩阵和截面特性)和渐近修正翘曲函数;利用一维非线性分析得到的板全局响应性能和翘曲函数精确重构三维场;采用极大化二维截面分析得到的弯曲和扭转刚度系数的优化策略对复合层合板进行优化,以达到通过改变复合材料铺设倾角的方法提高复合板所承受的最大弯扭屈曲荷载的目的。极大化刚度系数的优化策略能有效提高复合梁的屈曲临界荷载,并可节约优化搜索的时间和步骤,为工程技术人员提供了有效的设计手段。本发明实用性强,通用性高,可显著提高此类问题的解算速度和效率。
附图说明
图1为本发明三维场仿真及优化方法流程图;
图2为板变形前后坐标系图;
图3为复合材料层合板柱形弯曲问题模型图;
图4为σ11随厚度坐标的变化关系图;
图5为σ12随厚度坐标的变化关系图;
图6为σ22随厚度坐标的变化关系图;
图7为σ13随厚度坐标的变化关系图;
图8为σ23随厚度坐标的变化关系图;
图9为σ33随厚度坐标的变化关系图;
图10为ID3在极大化刚度系数优化策略下设计变量和目标函数的评估。
具体实施方式
下面结合具体实例和附图对本发明作进一步说明。
1、三维表达式
参见图2,是板变形前后坐标系图。复合层合板基准面上任意一点位置可由其笛卡尔坐标系xi表示,其中,xα是基准面上相互正交的坐标,x3是法向坐标(字母下标i、j、k代表1、2、3;字母下标α、β代表1、2,下同)。引入一组沿xi方向的正交参考坐标向量bi,变形前板上任一点的位置可由固定点O到xα确定点的位置向量描述:
当定义板的中间层为其基准面时
其中尖括号表示沿板厚度的积分,下同。
Bi=Cijbj Cij=Bi·bj (3)
变形板上任一点的位置向量可表示为
其中wi是翘曲分量,在这里视为未知的三维函数求解,以考虑包括古典层合理论未考虑的局部翘曲变形在内的所有变形。翘曲的引入使式(4)有6次冗余,需6个约束来求解方程,可通过适当定义和Bi消除冗余,与式(2)类似,可定义在板的中间层,这样,翘曲函数必须满足3个约束:
<wi(x1,x2,x3)>=0 (5)
将B3设为与板变形面垂直可确定其它两个约束。
基于旋转张量分解理论,由局部小旋转的条件,Jauman-Biot-Cauchy应变分量可表示为
Γij=1/2(Fij+Fji)-δij (6)
其中,Fij是变形梯度张量的混合基分量
假设应变很小,忽略一维广义应变和翘曲的乘积项,可得到一维广义应变和翘曲函数表示的三维应变场:
Γ=[Γ11 2Γ12 2Γ22 2Γ13 2Γ23 Γ33]T (8)
w=[w1 w2 w3]T (9)
ε=[ε11 2ε12 ε22 K11 K12+K21 K22]T (10)
板每单位面积上的应变能可表示为
其中D是6×6阶对称材料参数矩阵。
由于应变很小,可安全地忽略虚拟旋转中翘曲与荷载的乘积项,得到式(5)约束下荷载τiBi,βiBi,φiBi(分别作用在板的顶面、底面和厚度方向)虚拟位移产生的虚功:
W=-τTw|x3=h/2-βTw|x3=-h/2-<φTw> (12)
根据虚功原理,层合板的三维能量问题的完整表述为
δU-δW=0 (13)
总能量的泛函为
∏=U+W (14)
δ∏=0 (15)
在(5)式约束下的总能量中,仅翘曲是变化的。未知的翘曲函数可由对总能量最小化求得。
(14)式仅是原三维板弹性问题的另一表达形式。若直接进行求解,将会遇到与原三维板弹性问题相同的困难。若板由多层复合层构成,该计算过程将变得十分繁琐。本文借助渐近变分法渐近计算三维翘曲函数,可使计算得到简化。为考虑任何横截面形状和各向异性材料,首先应用有限元方法将三维翘曲场离散为一维有限元单元形式
w(x1,x2,x3)=S(x3)V(x1,x2) (16)
其中:S(x3)表示单元形函数,V为横向法线方向翘曲场的节点位移。
将(16)代入(13)中,可得到离散形式的总能量为
其中L=-S+Tτ-S-Tβ-<STφ>为荷载相关项。
新引入的几何形状和材料属性有关的变量包括:
E=<[ΓhS]T D[ΓhS]>,Dhε=<[ΓhS]T DΓε>
Dεε=<Γε T DΓε>,
式(5)的离散形式为
VTHψ=0 (19)
其中H=<STS>和ψ是E的正交化核心矩阵,ψTHψ=1。三维翘曲函数的求解问题转化为在式(19)约束下(17)式最小化问题。
2、降维后的近似能量公式推导
从所周知,弹性体的形态完全由其能量所确定。渐近变分法是一种有效的数学工具,可用二维能量尽可能准确地再现三维能量。对复合层合板结构,可选择h/l(h为板厚,l为板翘曲变形值)和广义二维应变的阶数ε作为渐近变分计算所需的小参数。
应用渐近变分法,首先需根据泛函的不同阶数找到主导项。由于总能量中只有翘曲是变化的,只需找到仅含有翘曲的主导项或翘曲与其它量乘积项(如广义应变和荷载)的主导项。
对式(17)零阶近似后泛函的主导项为
在式(19)约束下(20)式的欧拉-拉格朗日方程可借助拉格朗日乘子Λ通过对变量的常规计算得到
EV+Dhεε=HψΛ (21)
考虑内核矩阵ψ的属性,拉格朗日乘子Λ为
A=ψTDhεε (22)
将式(22)代回式(21),可得
EV=(HψψT-I)Dhεε (23)
由于式(23)等式右边与V的零空间E正交,存在唯一解与零空间E线性独立,可选择任何约束以方便求解:
V=V*+ψλ (24)
其中:V*为线性体系解,λ可由式(19)确定
λ=-ψTHV* (25)
式(25)代入式(24),则式(19)约束下式(17)的最小化解为
将式(26)代入式(20)得到渐近修正到零阶近似的总能量泛函为
这一近似能量与古典分层板理论得出的结果相吻合,但未建立在Kirchhoff动力学假设基础上,且尽管能量相同,但推导出的横向法线方向应变并不为零。注意到零阶近似得到的翘曲是ε阶,根据渐近变分法,为接受零阶近似结果,必须检查下一阶近似是否比这还要高。为得到一阶近似翘曲,对零阶近似翘曲进行如下摄动
V=V0+V1 (28)
将式(28)代入式(7)及(17),可得到总能量泛函一阶近似的主导项
同零阶近似一样,可求解一阶翘曲场
V1=V1αε,α+V1L (30)
代入(29)式,可得到渐近修正到二阶的总能量泛函:
其中:
(32)
3、将近似能量转换为Reissner模型形式
尽管总能量经二次渐近修正,但因包含广义应变量的导数,涉及到更加复杂的边界条件。为得到实用的能量方程,可将这种近似能量转换为Reissner模型形式。
在Reissner模型中有两个增加的横向剪切应变自由度γ=[2γ13 2γ23]T,将其纳入横向法线方向的旋转变量中。古典应变量的Reissner形式可表示为
ε=R-Dαγ,α (33)
其中:
将(33)式代入(31)式,可得到Reissner应变形式表示的修正到第二阶的总能量泛函
4、三维场重构关系推导
经渐近修正的Reissner模型形式的近似能量尽可能接近总能。但这还远不足以准确重现复合层合板的三维位移场,还需提供重构关系以完善降维模型。本文通过得到二维变形和翘曲函数准确重构复合层合板三维场。
对经二阶渐近修正的应变能,可基于渐近修正的严格定义,由式(1),(3),(4)得到一阶渐近修正的三维变形
其中:U3d是三维变形的列阵;u2d是板的二维变形;Cij由式(3)得到的全局旋转张量分量。
由式(7)可得到一阶渐近修正的三维应变场
然后,可使用材料本构关系得到三维应力σij。
5、优化问题和算法
优化的主要目的是寻找复合层合梁最优铺层倾角,以提高弯扭屈曲载荷。在这里考虑在屈曲临界荷载约束下极大化刚度矩阵的弯曲和扭转项。由刚度矩阵的扭转项s44和弯曲项s55作为目标函数
minmize ψ0(θ)=1/(D44+D55) (38)
优化算法选用以梯度为基础的序列二次规划算法,在当前设计点将目标函数做二阶近似,求解后得到新的设计变量,重复这一过程直至迭代收敛。该方法在连续的、凸设计空间被证明是有效的,在非凸空间可能会收敛到局部最优解,可从不同出发点进行优化设计以得到全局收敛解。
为了验证上述方法的可靠性,对一具有20层对称复合层合板([30°/-30°/30°/-30°]5,从底部到顶部编号)柱形弯曲问题进行分析。板长L=10cm(沿x1轴方向),板厚h=2.5cm,长厚比L/h=4。板两端简支,承受正弦变化荷载:
参见图4~9,为重构的沿厚度方向6个应力分量与精确解的对比结果。可看出,所有模型的平面应力分量σ11,σ12,σ22均与精确解吻合较好;对于面外横向应力分量σ13,σ23,本方法解可获得比古典层合理论和一阶剪切变形理论解更精确的结果;对于面外横向正应力σ33,古典层合理论和一阶剪切变形理论解皆无法获得满意的数值计算结果,而本方法解与精确解一致,且可在一般台式机上1s完成所有计算,而三维有限元计算需30s(台式机配置为P41.7G,512M内存,320G硬盘)。
优化过程所需要的初始设计变量和屈曲荷载列于表2,共考虑三种不同的初始铺层倾角:ID1、ID2、ID3。设计变量的优化范围定义为[-90°,90°]。采用极大化刚度系数的优化策略、铺层倾角为设计变量、屈曲荷载最大为优化目标进行优化,优化结果见表3。对比分析表明:ID1、D2、ID3的屈曲荷载经优化后都有明显增加,ID3的屈曲荷载由初始的400N增加到586N,增幅达46.5%;优化过程中使用不同的初始设计会产生不同的优化设计结果,这是因为在非凸问题中使用确定性算法的结果。在这种情况下,优化往往只能找到局部最优解,而不是所希望的整体优化设计。可通过从几个不同的初始点进行优化设计以获取最佳解。
参见图10,为极大化刚度系数的优化策略下ID3的目标函数和设计变量随迭代次数的收敛情况。从图10(a)可看出,经第一次函数调用后目标函数就基本收敛到最优解收敛速度很快。由图10(b)可看到,设计变量θ2由60°变化到20°,θ4由-45°变化到-5°,且基本上都是在第一次函数调用后就收敛到最佳倾角值,大大节约了优化搜索的时间和步骤,同时也表明该优化策略的高效性。
表1为复合材料工程弹性常数表;
表2为复合层合板的初始设计变量和屈曲荷载;
表3为使用极大化刚度系数优化策略进行优化设计结果。
表1
表2
表3
Claims (6)
1.一种基于渐近变分法的复合材料层合板仿真及优化方法,其特征在于,包括如下步骤:
1)基于旋转张量分解方法将复合材料层合板三维连续方程用一维广义应变和翘曲函数表示;
2)将非线性三维复合层合板模型分解为二维截面分析和沿法线方向的一维非线性分析;采用渐近变分算法实现简化模型的求解过程:通过二维截面分析得到本构关系,即刚度矩阵和截面特性和渐近修正翘曲函数;通过沿结线方向的一维非线性分析得到二维板的全局响应性能;
3)基于推导的重构关系,由计算得到的二维板变形、全局旋转张量分量和渐近修正翘曲函数重构三维变形;由一维广义应变和渐近修正翘曲函数重构三维应变场;使用材料本构关系得到三维应力;
4)以弯扭屈曲荷载最大为目标,以不发生屈曲为约束条件,层合复合材料铺设倾角为设计变量的复合梁优化设计。
3.根据权利要求1所述的基于渐近变分法的复合材料层合板仿真及优化方法,其特征在于:
所述基于渐近变分法将原非线性三维分析严格分解为沿厚度方向的一维分析和沿法线方向的二维非线性板分析;利用板固有的小参数将三维总能量泛函渐近扩展为系列二维泛函,∏=∏0+∏1(h/l)+∏2(h2/l2)+O(h2/l2),其中∏0,∏1,∏2分别为零阶、一阶、二阶近似二维能量泛函;由于假设应变很小,能量泛函∏中仅翘曲是变化的;根据变分原理,未知翘曲函数可通过泛函取驻值得到,δ∏=0;未知翘曲函数求解后,将其代回总能量泛函得到二维板分析所需的二维能量泛函;如零阶近似得到的二维板渐近变分模型可用于分析板的全局性能
其中:ε为含一维广义应变量的矩阵,ε=[ε112ε12ε22K11K12+K21K22]T;Dhε,Dεε为与几何形状和材料属性有关的矩阵;
更高阶近似可用来分析横向剪切效应,从而构建更精确板模型;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110032011 CN102096736B (zh) | 2011-01-28 | 2011-01-28 | 一种基于渐近变分法的复合材料层合板仿真及优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110032011 CN102096736B (zh) | 2011-01-28 | 2011-01-28 | 一种基于渐近变分法的复合材料层合板仿真及优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102096736A true CN102096736A (zh) | 2011-06-15 |
CN102096736B CN102096736B (zh) | 2013-08-14 |
Family
ID=44129831
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110032011 Expired - Fee Related CN102096736B (zh) | 2011-01-28 | 2011-01-28 | 一种基于渐近变分法的复合材料层合板仿真及优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102096736B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102841978A (zh) * | 2011-06-20 | 2012-12-26 | 波音公司 | 用于复合材料层合板的弯曲纤维路径的设计 |
CN103034766A (zh) * | 2012-12-15 | 2013-04-10 | 北京航空航天大学 | 一种确定层合板复合材料的铺层角方向与厚度的方法 |
CN103294862A (zh) * | 2013-05-29 | 2013-09-11 | 上海飞机制造有限公司 | 碳纤维增强树脂基复合材料加工变形的仿真模拟方法 |
CN103353916A (zh) * | 2012-09-10 | 2013-10-16 | 中国航空工业集团公司科学技术委员会 | 基于工程的复合材料层合板铺层优化后处理方法 |
CN103955587A (zh) * | 2014-05-13 | 2014-07-30 | 重庆大学 | 一种压电复合材料层合壳压电弹性分析方法 |
CN104091003A (zh) * | 2014-06-30 | 2014-10-08 | 南京航空航天大学 | 一种基础运动时柔性壳结构大变形响应的有限元建模方法 |
CN105631073A (zh) * | 2014-10-31 | 2016-06-01 | 中国航空工业集团公司西安飞机设计研究所 | 一种优化制造复合材料波纹梁的方法 |
CN106682292A (zh) * | 2016-12-14 | 2017-05-17 | 西安交通大学 | 一种降维模拟退火算法的叶根结构优化方法 |
CN107526898A (zh) * | 2017-09-13 | 2017-12-29 | 大连理工大学 | 一种变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计方法 |
CN107808048A (zh) * | 2017-10-26 | 2018-03-16 | 哈尔滨工程大学 | 一种均匀截面梁结构的振动分析方法 |
CN107944177A (zh) * | 2017-12-07 | 2018-04-20 | 重庆大学 | 一种梯形直波纹板生产方法 |
CN109002630A (zh) * | 2018-08-01 | 2018-12-14 | 北京大学 | 一种超弹性材料的快速仿真方法 |
CN109032074A (zh) * | 2018-07-12 | 2018-12-18 | 苏州大学 | 一种手机金属中框机加工翘曲变形重构方法 |
CN109740268A (zh) * | 2019-01-07 | 2019-05-10 | 重庆大学 | 一种建筑用夹芯板生产方法 |
CN115421460A (zh) * | 2022-09-15 | 2022-12-02 | 重庆大学 | 基于计算机数值模拟的铸造残余应力控制优化方法及应用 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060036978A1 (en) * | 2004-08-11 | 2006-02-16 | Fujitsu Limited | Board design aiding apparatus, board design aiding method and board design aiding program |
CN101089861A (zh) * | 2007-06-22 | 2007-12-19 | 北京航空航天大学 | 树脂基复合材料l形层板在成型中树脂流动与层板变形的模拟方法 |
-
2011
- 2011-01-28 CN CN 201110032011 patent/CN102096736B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060036978A1 (en) * | 2004-08-11 | 2006-02-16 | Fujitsu Limited | Board design aiding apparatus, board design aiding method and board design aiding program |
CN101089861A (zh) * | 2007-06-22 | 2007-12-19 | 北京航空航天大学 | 树脂基复合材料l形层板在成型中树脂流动与层板变形的模拟方法 |
Non-Patent Citations (1)
Title |
---|
钟秩峰等: ""用变分渐近法进行复合材料层合板仿真及三维场重构"", 《复合材料学报》, vol. 27, no. 4, 31 August 2010 (2010-08-31), pages 174 - 179 * |
Cited By (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102841978A (zh) * | 2011-06-20 | 2012-12-26 | 波音公司 | 用于复合材料层合板的弯曲纤维路径的设计 |
CN102841978B (zh) * | 2011-06-20 | 2017-01-18 | 波音公司 | 用于复合材料层合板的弯曲纤维路径的设计 |
CN103353916A (zh) * | 2012-09-10 | 2013-10-16 | 中国航空工业集团公司科学技术委员会 | 基于工程的复合材料层合板铺层优化后处理方法 |
CN103034766A (zh) * | 2012-12-15 | 2013-04-10 | 北京航空航天大学 | 一种确定层合板复合材料的铺层角方向与厚度的方法 |
CN103034766B (zh) * | 2012-12-15 | 2016-05-25 | 北京航空航天大学 | 一种确定层合板复合材料的铺层角方向与厚度的方法 |
CN103294862A (zh) * | 2013-05-29 | 2013-09-11 | 上海飞机制造有限公司 | 碳纤维增强树脂基复合材料加工变形的仿真模拟方法 |
CN103955587B (zh) * | 2014-05-13 | 2017-07-07 | 重庆大学 | 一种压电复合材料层合壳压电弹性分析方法 |
CN103955587A (zh) * | 2014-05-13 | 2014-07-30 | 重庆大学 | 一种压电复合材料层合壳压电弹性分析方法 |
CN104091003A (zh) * | 2014-06-30 | 2014-10-08 | 南京航空航天大学 | 一种基础运动时柔性壳结构大变形响应的有限元建模方法 |
CN104091003B (zh) * | 2014-06-30 | 2017-05-10 | 南京航空航天大学 | 一种基础运动时柔性壳结构大变形响应的有限元建模方法 |
CN105631073B (zh) * | 2014-10-31 | 2019-05-21 | 中国航空工业集团公司西安飞机设计研究所 | 一种优化制造复合材料波纹梁的方法 |
CN105631073A (zh) * | 2014-10-31 | 2016-06-01 | 中国航空工业集团公司西安飞机设计研究所 | 一种优化制造复合材料波纹梁的方法 |
CN106682292A (zh) * | 2016-12-14 | 2017-05-17 | 西安交通大学 | 一种降维模拟退火算法的叶根结构优化方法 |
CN106682292B (zh) * | 2016-12-14 | 2019-10-11 | 西安交通大学 | 一种降维模拟退火算法的叶根结构优化方法 |
CN107526898A (zh) * | 2017-09-13 | 2017-12-29 | 大连理工大学 | 一种变刚度复合材料板壳结构精确建模分析与可靠度优化一体化设计方法 |
CN107526898B (zh) * | 2017-09-13 | 2019-12-27 | 大连理工大学 | 变刚度复合材料板壳结构建模分析与可靠度优化设计方法 |
CN107808048A (zh) * | 2017-10-26 | 2018-03-16 | 哈尔滨工程大学 | 一种均匀截面梁结构的振动分析方法 |
CN107808048B (zh) * | 2017-10-26 | 2020-12-25 | 哈尔滨工程大学 | 一种均匀截面梁结构的振动分析方法 |
CN107944177A (zh) * | 2017-12-07 | 2018-04-20 | 重庆大学 | 一种梯形直波纹板生产方法 |
CN109032074A (zh) * | 2018-07-12 | 2018-12-18 | 苏州大学 | 一种手机金属中框机加工翘曲变形重构方法 |
CN109032074B (zh) * | 2018-07-12 | 2020-11-17 | 苏州大学 | 一种手机金属中框机加工翘曲变形重构方法 |
CN109002630A (zh) * | 2018-08-01 | 2018-12-14 | 北京大学 | 一种超弹性材料的快速仿真方法 |
CN109002630B (zh) * | 2018-08-01 | 2023-06-06 | 北京大学 | 一种超弹性材料的快速仿真方法 |
CN109740268A (zh) * | 2019-01-07 | 2019-05-10 | 重庆大学 | 一种建筑用夹芯板生产方法 |
CN115421460A (zh) * | 2022-09-15 | 2022-12-02 | 重庆大学 | 基于计算机数值模拟的铸造残余应力控制优化方法及应用 |
Also Published As
Publication number | Publication date |
---|---|
CN102096736B (zh) | 2013-08-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102096736B (zh) | 一种基于渐近变分法的复合材料层合板仿真及优化方法 | |
Caliri Jr et al. | A review on plate and shell theories for laminated and sandwich structures highlighting the Finite Element Method | |
Herrema et al. | Penalty coupling of non-matching isogeometric Kirchhoff–Love shell patches with application to composite wind turbine blades | |
Madeo et al. | Post-buckling analysis of variable-angle tow composite plates using Koiter's approach and the finite element method | |
Yin et al. | Isogeometric locking-free plate element: a simple first order shear deformation theory for functionally graded plates | |
Yu et al. | Generalized Timoshenko theory of the variational asymptotic beam sectional analysis | |
Mohammadi et al. | An equivalent model for trapezoidal corrugated cores based on homogenization method | |
Ferreira et al. | Buckling analysis of isotropic and laminated plates by radial basis functions according to a higher-order shear deformation theory | |
Ye et al. | Three-dimensional vibration analysis of laminated functionally graded spherical shells with general boundary conditions | |
Yasin et al. | An efficient layerwise finite element for shallow composite and sandwich shells | |
Tornabene et al. | Static analysis of laminated composite curved shells and panels of revolution with a posteriori shear and normal stress recovery using generalized differential quadrature method | |
Zucco et al. | Mixed shell element for static and buckling analysis of variable angle tow composite plates | |
Shahbazi et al. | A mesh-free method using exponential basis functions for laminates modeled by CLPT, FSDT and TSDT–Part I: Formulation | |
Yu et al. | Asymptotically accurate 3-D recovery from Reissner-like composite plate finite elements | |
Nemeth | A treatise on equivalent-plate stiffnesses for stiffened laminated-composite plates and plate-like lattices | |
Yu et al. | Asymptotic generalization of Reissner–Mindlin theory: accurate three-dimensional recovery for composite shells | |
Pagani et al. | Exact solutions for static analysis of laminated, box and sandwich beams by refined layer-wise theory | |
Polit et al. | A multilayered/sandwich triangular finite element applied to linear and non-linear analyses | |
Nguyen-Minh et al. | Free vibration analysis of corrugated panels using homogenization methods and a cell-based smoothed Mindlin plate element (CS-MIN3) | |
Li et al. | Linear statics and free vibration sensitivity analysis of the composite sandwich plates based on a layerwise/solid-element method | |
Shahmohammadi et al. | Stability of laminated composite and sandwich FGM shells using a novel isogeometric finite strip method | |
Patel et al. | Buckling response of laminated composite stiffened plates subjected to partial in-plane edge loading | |
Chakravarty | On the modeling of composite beam cross-sections | |
Ye et al. | A spectral-sampling surface method for the vibration of 2-D laminated curved beams with variable curvatures and general restraints | |
CN104281730A (zh) | 一种大转动变形的板壳结构动响应的有限元分析方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130814 Termination date: 20160128 |
|
EXPY | Termination of patent right or utility model |