CN105447237B - 拉伸载荷作用下复合材料干涉连接应力分析方法 - Google Patents

拉伸载荷作用下复合材料干涉连接应力分析方法 Download PDF

Info

Publication number
CN105447237B
CN105447237B CN201510783235.7A CN201510783235A CN105447237B CN 105447237 B CN105447237 B CN 105447237B CN 201510783235 A CN201510783235 A CN 201510783235A CN 105447237 B CN105447237 B CN 105447237B
Authority
CN
China
Prior art keywords
stress
composite material
tensile load
coordinate system
coefficient
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.)
Active
Application number
CN201510783235.7A
Other languages
English (en)
Other versions
CN105447237A (zh
Inventor
张开富
武涛
李原
程晖
刘平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201510783235.7A priority Critical patent/CN105447237B/zh
Publication of CN105447237A publication Critical patent/CN105447237A/zh
Application granted granted Critical
Publication of CN105447237B publication Critical patent/CN105447237B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明提出一种拉伸载荷作用下复合材料干涉连接应力分析方法,该方法基于Lekhnitskii复势函数理论,在建立接头孔边变形与应力函数关系的基础上,通过孔边变形的求解和应力函数运算,得到孔边各处的应力分量。该方法可以在已知材料属性、外载荷大小和干涉量的条件下求解接头孔边的应力状态,使用方便、计算准确,可以为结构强度的分析和工艺参数的分析提供有利的理论支持;弥补现有方法在计算时间和使用方便性上的缺陷,为结构的强度分析和损伤预测提供依据,为复合材料结构干涉配合连接的设计提供参考。

Description

拉伸载荷作用下复合材料干涉连接应力分析方法
技术领域
本发明涉及复合材料连接应力分析技术领域,具体为一种拉伸载荷作用下复合材料干涉连接应力分析方法。
背景技术
复合材料因其比强度高、比刚度高、可设计性强等优点,在航空、航天和汽车等工业领域得到了越来越多的应用。复合材料在工程中的使用不可避免的导致连接结构的存在。对于承载较大的结构而言,机械连接是最常用的连接形式。但由于机械连接需要预制连接孔,造成材料不连续和应力集中,机械连接接头往往是结构中的最薄弱的区域,相关统计结果表明60%~80%的结构失效发生在接头部位。为提高复合材料机械连接接头的强度和疲劳寿命,美国麦道公司将已在金属结构中得到应用的干涉配合技术使用在复合材料结构上,试验研究表明干涉配合对复合材料接头同样有效。但由于工艺要求较高且相关工艺参数的选取较为复杂,目前该技术的应用尚处于发展状态。
拉伸载荷是连接结构承受的最主要外载形式之一。对于干涉配合的复合材料接头而言,由于存在干涉残余应力并受到复合材料各向异性的影响,其连接孔周围的应力分布沿角度变化剧烈,应力集中严重,极易成为损伤的源头,引起结构的破坏。应力水平的计算是分析结构强度和寿命的基础,对工艺参数的选择具有重要作用。快速、准确地分析干涉连接在外力作用下的应力状态,可以方便干涉量等工艺参数的选择,对于干涉配合技术的应用具有十分重要的意义。
目前已有学者围绕复合材料连接接头的应力状态进行了研究,主要分为两类:间隙配合接头在拉伸载荷作用下应力状态研究(Aluko O and Whitworth HA.“Analysis ofstress distribution around pin loaded holes in orthotropic plates”.CompositeStructures.2008;86:308-13);干涉连接接头周围干涉残余应力的研究(Kim S,He B,ShimC,and Kim D,“An experimental and numerical study on the interference-fit pininstallation process for cross-ply glass fiber reinforced plastics(GFRP)”,Composites:Part B.2013 54:153-162)。孔周应力状态的分析目前主要采用有限元数值模拟的方式(张俊琪,刘龙权,陈昆昆,汪海.“干涉配合对复合材料机械连接结构承载能力的影响”.上海交通大学学报,2013,47(11):1795-1806)。有限元数值模拟按照研究对象的几何尺寸、材料属性和受力状态建立模型进行分析,尽管具有较高的精度但需要较长的计算时间。而且由于模型与研究对象存在直接的对应关系,当研究对象的尺寸、材料、受力状态发生变化时,需要重复建模和分析过程,因而对于大量不同工艺参数及材料属性的对象,其所耗费的计算成本成倍增加。
由于承载能力强、便于拆卸等优点,飞机结构中大量使用机械连接接头。为了满足寿命、强度、密封性等要求,干涉配合在复合材料机械连接接头中得到应用。对于飞机结构件而言,外载作用下的应力状态是分析结构质量、安全性的关键因素,直接影响飞机产品的安全可靠。目前,针对拉伸载荷作用下复合材料干涉配合接头的孔边应力状态主要采用有限元数值模拟的方式,这种方法存在以下问题:1)计算时间代价很高。为保证精度,有限元模型需要划分一定规模的网格单元,因而其计算过程较长。此外,提高计算精度需要成倍增加计算规模,导致计算时间的进一步增加;2)模型具有一定的特异性。有限元模型针对特定材料、特定尺寸的结构建立,当分析对象变化时,需要重复建模及计算过程,而该过程耗时很长,使得其不适于大量对象的分析。
发明内容
为解决现有技术存在的问题,本发明提出了一种拉伸载荷作用下复合材料干涉连接应力分析方法,该方法基于Lekhnitskii复势函数理论,在建立接头孔边变形与应力函数关系的基础上,通过孔边变形的求解和应力函数运算,得到孔边各处的应力分量。该方法可以在已知材料属性、外载荷大小和干涉量的条件下求解接头孔边的应力状态,使用方便、计算准确,可以为结构强度的分析和工艺参数的分析提供有利的理论支持;弥补现有方法在计算时间和使用方便性上的缺陷,为结构的强度分析和损伤预测提供依据,为复合材料结构干涉配合连接的设计提供参考。
本发明的技术方案为:
所述一种拉伸载荷作用下复合材料干涉连接应力分析方法,其特征在于:复合材料板上某一点在笛卡尔坐标系下的应力分量σxyxy通过下面公式表示,所述笛卡尔坐标系为在复合材料板平面上,建立的以干涉连接孔中心为原点,以拉伸载荷方向为x轴正方向的笛卡尔坐标系;复合材料板上该点在笛卡尔坐标系下的坐标为(x,y);
其中Re表示取实部,μ1、μ2为与复合材料属性相关的复参数,Φ′1(z1)为应力函数Φ1(z1)的导数,Φ′2(z2)为应力函数Φ2(z2)的导数;
μ1、μ2通过公式
得到,其中i为虚数单位,k和n为与复合材料参数相关的系数,k=(E1/E2)1/2,n=[2(k-v12)+E1/G12]1/2,E1、E2为复合材料弹性模量,G12为复合材料剪切模量;
应力函数Φ1(z1)和Φ2(z2)为:
应力函数Φ1(z1)和Φ2(z2)的导数为:
其中A和B是与拉伸载荷以及复合材料参数相关的系数,
P为拉伸载荷,对应为μ1和μ2的共轭,Δ=R-r表示干涉量,R为连接螺栓半径,r为复合材料板干涉连接孔初始直径,D=(μ12)g/E1,g为与复合材料参数相关的系数,g=(1-v12v21)/E2+k/G12,v12、v21为复合材料的泊松系数;
系数U1,U2,V1,V2通过以下公式得到:
其中u0为与拉伸载荷大小相关的位移量,c为与复合材料参数相关的系数,
c=-(10(k-v12+nk+2n)-11n(1-k))/(10(k-v12+nk+2n)-n(1-k))
ζj为与坐标有关的参数,j=1,2:
p1,p2,q1,q2为与材料参数有关的系数,
其中a11、a12、a22为复合材料板的柔度系数。
进一步的优选方案,所述一种拉伸载荷作用下复合材料干涉连接应力分析方法,其特征在于:复合材料板上某一点在极坐标系下的应力分量σrθ通过下面公式表示,
所述极坐标系为在复合材料板平面上,建立的以干涉连接孔中心为原点,以拉伸载荷方向为0°方向建立的极坐标系,θ表示点在极坐标系中的角度。
有益效果
本发明具有的优点如下:
(1)本发明建立了拉伸载荷作用下的复合材料干涉配合接头应力分布的计算模型,可以在Matlab软件中编程求解,具有精度高、速度快、计算代价低、使用方便等优点。
(2)该方法具有一定的通用性,对于不同材料属性、不同干涉量的接头只需修改材料及结构参数,不需重复建模及编程,可以降低分析工作的耗时。
(3)该方法可以实现整个复合材料板平面内任何点位的应力计算,可以为结构设计及强度预测提供支持。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1:结构受力状态及主要几何尺寸;
1表示金属螺栓,2表示复合材料板;
图2:螺栓孔边区域变形位移分析;
图3:拉伸载荷作用下螺栓孔边应力分布结果。
具体实施方式
下面详细描述本发明的实施例,所述实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
本发明的目的在于提出一种拉伸载荷作用下复合材料干涉连接应力分析方法,弥补现有方法在计算时间和使用方便性上的缺陷,为结构的强度分析和损伤预测提供依据,为复合材料结构干涉配合连接的设计提供参考。该方法计算成本低、计算量小并具有较高的精度,可以满足工程应用中的需求。
本发明方法涉及到的主要部件为螺栓(1)、复合材料板(2),其中螺栓半径为R、复合材料板螺栓孔初始直径为r(r<R)、复合材料板厚为h、复合材料宽度为W、螺栓孔边距为E。螺栓通过冷挤压方式安装到复合材料螺栓孔中。
申请人基于上述研究目的,提出本发明的研究思路如下:
(1)根据拉伸外载荷与孔周区域变形的关系,分析拉伸载荷引起的孔周变形位移。由于拉伸载荷受力沿板的长度方向对称,因而仅需分析模型的一半。根据变形特点,可以将孔周按角度分为5个区域:区域Ⅰ为连接件与复合材料板的纯挤压区域,区域Ⅱ为连接件与复合材料板相互挤压并存在滑动的区域,区域Ⅲ为连接件与复合材料板间不存在法向挤压的区域,区域Ⅳ为连接件与复合材料板相互分离并存在滑动的区域,区域Ⅴ为连接件与复合材料板相互分离的区域。忽略螺栓的变形,以孔中心为原点,拉伸载荷方向为0°方向建立极坐标系,用u、v表示沿0°和90°的位移分量,各区域的范围及位移关系如下:
区域 Ⅰ θ=0 u=u0 and v=0
区域 Ⅱ 0<θ<π/2 (u0-u)cosθ=vsinθ
区域 Ⅲ θ=π/2 u=u0/c and v=0
区域 Ⅳ π/2<θ<π (u0-u)cosθ=-vsinθ and σr≤0
区域 Ⅴ θ=π u=u0 and σr≤0
其中σr为某一点的径向应力,θ表示孔上一点在坐标系中的角度(逆时针旋转为正),u0为与外载荷大小相关的位移量,其数值可由公式(1)给出
P为外载荷,k,n,g,c均为与材料参数相关的系数,其计算公式为
k=(E1/E2)1/2,
n=[2(k-v12)+E1/G12]1/2,
g=(1-v12v21)/E2+k/G12,
c=-(10(k-v12+nk+2n)-11n(1-k))/(10(k-v12+nk+2n)-n(1-k))。
计算公式中v12,v21为复合材料的泊松系数,E1,E2为弹性模量,G12为剪切模量。最终可以将位移分量u和v表示为公式(2)的形式
系数U1,U2,V1,V2用于确定应力函数。在(2)式带入区域I、Ⅱ、Ⅲ的位移关系得如下关系式
即U1,U2,V1,V2均可由u0和c的值得到。
(2)对于受到拉伸载荷作用的干涉配合接头而言,其应力状态可以视为干涉应力与拉伸应力的叠加。在板平面上,以孔中心为原点,拉伸载荷方向为x轴正方向建立笛卡尔坐标系。依据Lekhnitskii复势函数理论,根据变形位移与应力状态的关系,其应力函数Φ表达式为
式中Δ=R-r表示干涉量,D=(μ12)g/E1。μ1,μ2为与材料属性相关的复参数,其值由式(5)求得。
ζj(j=1,2)为与坐标有关的参数。对复合材料板上任意一点(x,y)而言,其坐标的复参数表达形式为zj=x+μjy,对该表达式进行保角变化得ζj,如公式(6)
p1,p2,q1,q2为与材料参数有关的系数,其值可通过方程组(7)求解得到。amn(m,n=1,2,6)为复合材料板的柔度系数。
A和B是与外载荷和材料参数相关的系数,其表达式如公式(8),i为虚数单位。
(3)由于应力分量是由应力函数的导数表达的,对(3)式求导并进一步化简得到
(4)在笛卡尔坐标系下的应力分量σxyxy可由应力函数的导数公式(9)求得,表达式如(10)所示。考虑到接头孔的形状特点,应力分量采用极坐标的表达意义更为清晰,便于理解。借助转轴公式可以得到极坐标系下的应力分量σrθ表达式(11)
基于上述研究过程,本发明的方法为复合材料板上某一点在笛卡尔坐标系下的应力分量σxyxy通过下面公式表示,所述笛卡尔坐标系为在复合材料板平面上,建立的以干涉连接孔中心为原点,以拉伸载荷方向为x轴正方向的笛卡尔坐标系;复合材料板上该点在笛卡尔坐标系下的坐标为(x,y);
其中Re表示取实部,μ1、μ2为与复合材料属性相关的复参数,Φ′1(z1)为应力函数Φ1(z1)的导数,Φ′2(z2)为应力函数Φ2(z2)的导数;
μ1、μ2通过公式
得到,其中i为虚数单位,k和n为与复合材料参数相关的系数,k=(E1/E2)1/2,n=[2(k-v12)+E1/G12]1/2,E1、E2为复合材料弹性模量,G12为复合材料剪切模量;
应力函数Φ1(z1)和Φ2(z2)为:
应力函数Φ1(z1)和Φ2(z2)的导数为:
其中A和B是与拉伸载荷以及复合材料参数相关的系数,
P为拉伸载荷,对应为μ1和μ2的共轭,Δ=R-r表示干涉量,R为连接螺栓半径,r为复合材料板干涉连接孔初始直径,D=(μ12)g/E1,g为与复合材料参数相关的系数,g=(1-v12v21)/E2+k/G12,v12、v21为复合材料的泊松系数;
系数U1,U2,V1,V2通过以下公式得到:
其中u0为与拉伸载荷大小相关的位移量,c为与复合材料参数相关的系数,
c=-(10(k-v12+nk+2n)-11n(1-k))/(10(k-v12+nk+2n)-n(1-k))
ζj为与坐标有关的参数,j=1,2:
p1,p2,q1,q2为与材料参数有关的系数,
其中a11、a12、a22为复合材料板的柔度系数。
进一步的,复合材料板上某一点在极坐标系下的应力分量σrθ通过下面公式表示,
所述极坐标系为在复合材料板平面上,建立的以干涉连接孔中心为原点,以拉伸载荷方向为0°方向建立的极坐标系,θ表示点在极坐标系中的角度。
本实施例中,复合材料平板通过钛合金高锁螺栓连接,连接采用干涉配合,结构的材料属性及几何参数如表1、表2、表3所示,所受外载为匀布拉伸载荷P=500Mpa。
表1干涉配合接头几何尺寸(mm)
注:E为螺栓孔中心到复合材料构件一端的最短距离,W为复合材料板的宽度,h为复合材料板的厚度,r为复合材料板预制连接孔的半径,R为螺栓半径,如图1所示。
表2螺栓材料及其性能参数
注:剪切模量G可以通过弹性模量E和泊松比v求得。
表3复合材料板材料及其性能参数
注:E1,E2,E3为复合材料的弹性模量,ν121323为复合材料的泊松比,G12,G13,G23为复合材料的剪切模量。
拉伸载荷作用下的接头孔边应力计算具体实施步骤如下:
(1)将材料弹性模量、剪切模量、泊松比带入系数计算公式得
k=3.8357,g=0.0011,n=6.2475
将系数值、材料属性和外载荷数值带入位移计算公式得到
u0=0.00525618590918473
c=-1.424378995998533
U1=0.0044731727,U2=0.0007830131
V1=0.0060391990,V2=0.0007830131
孔周区域位移分量u,v为
(2)将材料参数带入复参数μ1,μ2的计算公式得
μ1=5.55730240104i
μ2=0.69021012986i
将材料参数和外载荷带入系数A、B的计算公式得
A=14.6944803221
B=-3.0196589997
将材料常数带入系数D,p1,p2,q1,q2的计算公式得
D=0.0000000408025506045i
p1=-0.0002437782
p2=-0.0000062217
q1=(-0.0000345764)i
q2=(-0.0001682581)i
(3)将系数带入应力函数的导数表达式得
(4)笛卡尔坐标系下的应力分量表达式如下:
极坐标系下的应力分量表达式如下
将应力函数导数结果带入应力分量表达式即可获得复合材料板孔周区域应力分量的计算公式,带入所求位置的角度即可得到相应点位的应力状态。有限元模拟结果对比如图3所示,由图可知计算得到的应力分布与有限元数值仿真结果一致,且得到的应力数值较为精确,可以满足工程需要。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。

Claims (2)

1.一种拉伸载荷作用下复合材料干涉连接应力分析方法,其特征在于:复合材料板上某一点在笛卡尔坐标系下的应力分量σxyxy通过下面公式表示,所述笛卡尔坐标系为在复合材料板平面上,建立的以干涉连接孔中心为原点,以拉伸载荷方向为x轴正方向的笛卡尔坐标系;复合材料板上该点在笛卡尔坐标系下的坐标为(x,y);
其中Re表示取实部,μ1、μ2为与复合材料属性相关的复参数,Φ′1(z1)为应力函数Φ1(z1)的导数,Φ′2(z2)为应力函数Φ2(z2)的导数;
μ1、μ2通过公式
得到,其中i为虚数单位,k和n为与复合材料参数相关的系数,k=(E1/E2)1/2,n=[2(k-v12)+E1/G12]1/2,E1、E2为复合材料弹性模量,G12为复合材料剪切模量;
应力函数Φ1(z1)和Φ2(z2)为:
应力函数Φ1(z1)和Φ2(z2)的导数为:
其中A和B是与拉伸载荷以及复合材料参数相关的系数,
P为拉伸载荷,对应为μ1和μ2的共轭,Δ=R-r表示干涉量,R为连接螺栓半径,r为复合材料板干涉连接孔初始直径,D=(μ12)g/E1,g为与复合材料参数相关的系数,g=(1-v12v21)/E2+k/G12,v12、v21为复合材料的泊松系数;
系数U1,U2,V1,V2通过以下公式得到:
其中u0为与拉伸载荷大小相关的位移量,c为与复合材料参数相关的系数,
c=-(10(k-v12+nk+2n)-11n(1-k))/(10(k-v12+nk+2n)-n(1-k))
ζj为与坐标有关的参数,j=1,2:
zj=x+μjy
p1,p2,q1,q2为与材料参数有关的系数,
其中a11、a12、a22为复合材料板的柔度系数。
2.根据权利要求1所述一种拉伸载荷作用下复合材料干涉连接应力分析方法,其特征在于:复合材料板上某一点在极坐标系下的应力分量σrθ通过下面公式表示,
所述极坐标系为在复合材料板平面上,建立的以干涉连接孔中心为原点,以拉伸载荷方向为0°方向建立的极坐标系,θ表示点在极坐标系中的角度。
CN201510783235.7A 2015-11-16 2015-11-16 拉伸载荷作用下复合材料干涉连接应力分析方法 Active CN105447237B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510783235.7A CN105447237B (zh) 2015-11-16 2015-11-16 拉伸载荷作用下复合材料干涉连接应力分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510783235.7A CN105447237B (zh) 2015-11-16 2015-11-16 拉伸载荷作用下复合材料干涉连接应力分析方法

Publications (2)

Publication Number Publication Date
CN105447237A CN105447237A (zh) 2016-03-30
CN105447237B true CN105447237B (zh) 2019-01-08

Family

ID=55557408

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510783235.7A Active CN105447237B (zh) 2015-11-16 2015-11-16 拉伸载荷作用下复合材料干涉连接应力分析方法

Country Status (1)

Country Link
CN (1) CN105447237B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110175399B (zh) * 2019-05-27 2022-10-21 四川领航空天智能科技有限公司 一种考虑衬套膨胀效应的衬套螺栓干涉量计算方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102622472A (zh) * 2012-02-27 2012-08-01 西北工业大学 分析复合材料机械连接钉载和应力分布及应力水平的方法
CN102682170A (zh) * 2012-05-11 2012-09-19 中国航空工业集团公司西安飞机设计研究所 一种复合材料湿装配连接处材料性能的处理方法
CN103559390A (zh) * 2013-10-22 2014-02-05 北京航空航天大学 一种基于平均失效指数的复合材料π形胶接连接结构拉伸强度预测方法
CN104794299A (zh) * 2015-04-29 2015-07-22 西北工业大学 一种复合材料干涉配合接头应力分布计算方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120259593A1 (en) * 2011-04-07 2012-10-11 El-Zein Mohamad S Method for the prediction of fatigue life for welded structures

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102622472A (zh) * 2012-02-27 2012-08-01 西北工业大学 分析复合材料机械连接钉载和应力分布及应力水平的方法
CN102682170A (zh) * 2012-05-11 2012-09-19 中国航空工业集团公司西安飞机设计研究所 一种复合材料湿装配连接处材料性能的处理方法
CN103559390A (zh) * 2013-10-22 2014-02-05 北京航空航天大学 一种基于平均失效指数的复合材料π形胶接连接结构拉伸强度预测方法
CN104794299A (zh) * 2015-04-29 2015-07-22 西北工业大学 一种复合材料干涉配合接头应力分布计算方法

Also Published As

Publication number Publication date
CN105447237A (zh) 2016-03-30

Similar Documents

Publication Publication Date Title
CN104794299B (zh) 一种复合材料干涉配合接头应力分布计算方法
Vescovini et al. A variable-kinematic model for variable stiffness plates: Vibration and buckling analysis
Aydogdu A new shear deformation theory for laminated composite plates
Altenbach Theories for laminated and sandwich plates: A review
Chishti et al. Experimental investigation of damage progression and strength of countersunk composite joints
Cognard et al. Numerical analysis and optimisation of cylindrical adhesive joints under tensile loads
Song et al. Stress distribution modeling for interference-fit area of each individual layer around composite laminates joint
Takeda Micromechanics model for three-dimensional effective elastic properties of composite laminates with ply wrinkles
CN109117504B (zh) 一种双向功能梯度曲壳振动分析方法
Her et al. Adhesively bonded patch repair of composite laminates
Fan et al. Dynamic buckling of cylindrical shells with arbitrary axisymmetric thickness variation under time dependent external pressure
CN105447237B (zh) 拉伸载荷作用下复合材料干涉连接应力分析方法
Jonnalagadda et al. An analytical model for composite tubes with bend–twist coupling
Aabid et al. Stress concentration analysis of a composite patch on a hole in an isotropic plate
Avhad et al. On the deformation of laminated composite and sandwich curved beams
Kang et al. Parametric study on the buckling load after micro-bolt repair of a composite laminate with delamination
Zhu et al. Shear buckling of ship plates with different holes
CN112541284B (zh) 一种纤维增强聚合物基复合材料疲劳寿命的计算方法
Sayyad et al. Bending and free vibration analysis of isotropic plate using refined plate theory
Kaleel et al. Integration of CUF micromechanics framework into NASMAT for multiscale analysis of fiber-reinforced composites
Ke et al. A theoretical model used for determining the stiffness of composite leaf springs with a main spring and an auxiliary spring
Miri et al. Out-of-plane stresses in composite shell panels: layerwise and elasticity solutions
Yazdani Sarvestani et al. Analysis of free edge stresses in composite laminates using higher order theories
Starovoytov et al. Cylindrical bending of an elastic rectangular sandwich plate on a deformable foundation
Zhang et al. Research on vibration transfer characteristics of the bolt connection structure

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant