P01_Chap_04_几何非线性.docVIP

4-1 概要 当结构发生较大变形或旋转时，几何形状发生较大变化，致使应变-位移关系图呈非线性。如图4.1.1所示，A状态结构受荷后发生了较大位移变成（B）状态时的荷载-位移关系图呈非线性。 图4.1.1 非线性分析概念图 几何非线性分析时，变形前形状（初始形状）与变形后的形状是不一样的。几何非线性分析时的结构刚度（K）变成了位移（u）的函数。 一般几何非线性分析使用的方法有TLF(total Lagrangian formulation)与 ULF(updated Lagrangian formulation)。TLF方法中使用的应力与应变为PK2应力(2nd Piola-Kirchhoff stress)与Green应变(Green-Lagrange strain)，ULF方法中使用的应力与应变为Cauchy应力(Cauchy stress)与线性欧拉应变增量(linear Euler strain increment)。TLF方法中的应变与应力表现的是变形前结构形状，ULF方法表现的是变形后结构形状。例如，线单元分析时的PK2应力表示力与变形前面积的关系，Cauchy应力表示力与变形后面积的关系。 一个结构随时间发生位移的状态如图4.1.2所示。图中表示初始状态（变形前状态），表示当前状态（变形后状态）。表示初始状态的空间坐标，表示经过时间之后的空间坐标。 图4.1.2 结构的初始及当前状态 Green应变使用以下式计算： (4.1.1) 其中， : Green应变(Green-Lagrange strain) : 变形斜率(deformation gradient) 对线单元的结构类型可将式(4.1.1)整理如下： (4.1.2) 式(4.1.2)中，为变形前长度，为变形后长度。当变形微小时，Green应变将回归为线性应变，见式(4.1.3)。 (4.1.3) 为了查看实体单元的PK2应力的意义，下面推导与Cauchy应力的关系。不管是使用Cauchy应力还是PK2应力，能量是一样的，虚功的计算式如下： (4.1.4) 其中, : PK2应力张量 : Cauchy应力张量 : 应变张量变分量， 将Green应变与变形斜率(deformation gradient)变分后如下： (4.1.5) (4.1.6) 将式(4.1.6)代入至式(4.1.5)整理如下： (4.1.7) 将式(4.1.7)代入至(4.1.4)整理得出式(4.1.8)。 (4.1.8) Cauchy应力与PK2的应力关系式如下： ， (4.1.9) 利用PK2应力和Green应变计算在时间时刻的虚功如下；其中为外力做的虚功。 (4.1.10) (4.1.11) 将式(4.1.11)代入至式(4.1.10)整理后可以得出下式： (4.1.12) 其中， : 线性项 : 非线性项 采用泰勒级数展开如下： (4.1.13) 将式(4.1.13)代入至式(4.1.12)的第一项，忽略2次以上的项后，第一项将变成线性的公式(4.1.14)。 (4.1.14) 利用，整理式(4.1.12)的第二项如下： (4.1.15) 将式(4.1.12)右侧的各项表示为式(4.1.16)。 (4.1.16) 其中， : 将式(4.1.12)的右侧项利用式(4.1.16)整理后得出下式(4.1.17)。 (4.1.17) 将式(4.1.14-15,17)代入至式(4.1.12)后得出式(4.1.18)。从式(4.1.18)中可以看出，左侧第一项为线性应变做的虚功，第二项为非线性应变做的虚功。 (4.1.18) 式(4.1.18)可表示为式(4.1.19)，其中刚度矩阵与由形函数决定，详细说明请参照下一节。 (4.1.19) 4-2 桁架单元 桁架(Truss)单元在单元坐标系上具有位移、、，位移利用形函数表示如下： (4.2.1) 桁架单元的应力与应变如下： (4.2.2) 其中, : 单元坐标系 式(4.1.18)中虚功应变（线性）项按下式表示。 (4.2.3) 将式(4.2.3)用虚位移与矩阵的乘积表示如下： (4.2.4) 应变增量的线性项也可使用类似的形态表现。 (4.2.5) 式(4.2.4)与式(4.2.5)的位移-应变关系矩阵如下： (4.2.6) (4.2.7) 式(4.1.18)中虚应变的非线性项表现如下： (4.2.8) 将式(4.2.8)用虚位移与矩阵的乘积表示如下： (4