固体力学非线性课程报告.docVIP

固体力学非线性课程报告 结构非线性 载荷与变形曲线之间呈非线性的关系，也即显示了非线性结构的基本特征—变化的结构钢性。引起结构非线性变化的原因有很多，主要是三种类型： 1）状态变化（包括接触）； 许多普通的结构表现出一种与状态有关的非线性行为，如电缆的松弛与紧绷，轴套的接触与非接触，这些系统的结构刚度都是由于状态的改变而突然变化。状态的改变可能和载荷有关，也可能与某种外部原因有关。其中接触问题很典型，它的特点是：属于边界非线性问题，边界条件不在是定解条件，而是待求结果，两接触体之间的接触面积与压力随外载荷的变化而变化，并且与接触体的刚性有关。其中接触模式可以分为三类：点对点，点对面，面对面。目前对这类问题的求解主要是引用经典理论，借助有限元数值方法进行研究。 2）几何非线性 几何非线性就是位移与应变呈非线性（微分)关系,这是因为在几何非线性条件下，物体或结构会发生较大的位移或转动（也叫有限变形），此时不同时刻的物体具有差别不能忽略的不同构型，其单元刚度矩阵会发生改变或者是单刚阵在向总刚阵变换时发生改变，它的位移的导数的二次项不能忽略。 对于大变形问题，连续介质力学理论对经历大变形后的变形，依据不同的参考对象（初始构型与现时构型）有严格的定义，在此基础上定义的应变有两种：格林应变张量（以初始构型为参考），它是客观的、不随刚体转动而变换的对称张量；另一种为阿尔曼斯应变张量（以现时构型为参考）在大变形分析中经常用到它们的增量形式，均有线性部分和非线性部分构成。 3）材料非线性 应力与应变之间呈现非线性的关系。材料的非线性行为异常丰富，有非线性弹性行为，塑性行为，粘弹性行为。非线性弹性行为：当材料由于应力达到某种临界值而出现应力与应变间的非线性变化关系；塑性行为：当物体由于载荷超过某个临界值（弹性极限）而产生的永久变形，即有不可恢复的应变产生；粘弹性行为（包括松弛和蠕变）：具有介于弹性和黏性液体之间的材料性质，其形变不仅和当时作用力的大小有关，而且和温度的改变、力的作用时间以及加载历史都有关，同时应变速率、应力速率对应变也有影响。除此而外还包括非线性行为的耦合。蠕变：应力不变，应变随时间增加而增长；松弛：应变不变，应力随时间增加而减小；弹性后效：加载或卸载时，弹性应变滞后于应力的现象。 几何非线性与线性问题的区别 a)应力应变关系是线性的，但计算位移应变关系时位移的高阶导数项的影响不能忽略，因而位移应变的关系呈非线性（微分）关系。 b)对于有限变形问题，如大位移-大应变，其应力-应变关系也是非线性的。 c)几何非线性问题的平衡方程组建立在结构变形后的位移形态上，而这个位移形态在求解过程中总是变化的。 d)随着有限位形的变化，材料的本构方程亦发生变化，采用不同的参考位移将得到不同的本构方程。 非线性分析要点 1）非线性结构的行为需要一系列带校正的线性近似方程来求解，就是将载荷分成一系列的载荷增量，在一个载荷增量求解完成后继续下一个之前，程序调整刚度矩阵以反映非线性变化，但是纯粹增量式求解会导致结果失去平衡。 2）使用牛顿-拉普森平衡迭代可以克服这种难度。就是在某个容限内，它使每一个载荷增量的末端解都达到平衡收敛，在每次求解前估算出残差矢量，之后程序使用非平衡载荷进行线性求解，并核查收敛性。 非线性求解的组织级别 3）非线性求解被分成三个操作级别：载荷步、子步、平衡迭代。 载荷步由用户自定义，在每一个载荷步，通过逐步加载可以控制程序执行多次求解，在每一个子步内，程序进行一系列的平衡迭代以获得问题的收敛解。 悬臂梁模型理论与数值分析 题目描述 一个悬臂梁，在自由端受竖直向下的集中载荷P，求梁自由端的最大位移与转角。 在材料力学中基于小变形假设可以知道梁端的最大挠度和最大转角，分别为 假设梁长度为L=1m,横截面为圆形，半径为0.04m,受载荷为P=100N， E=2E11Pa，v=0.3，将各量带入上面公式可得到理论值:梁端的最大位移为 0.0082935cm，最大转角为1.244E-4rad。 用ANSYS建模分析 划分网格，施加约束和载荷后进行求解，不打开大位移模式进行静力分析，得到如下结果： 通过计算结果列表可知，ANSYS分析得到的梁端的最大位移为0.0083776cm,最大转角为1.246E-4rad，由此可知与理论值相比误差很小符合材料力学小变形的假设。但是在不断增大载荷的过程中（假设以10倍的顺序增大），计算得到的最大位移、最大转角也以10倍的顺序增大。 当载荷增大时，求得的最大位移为8.2935m，而梁总长才1m，所以此计算结果明显与实际情形不符，它已经发生了较大的位移，需使用新的计算准则来求解。考虑大转动的影响，打开大位移模式，重新计算，即可得到正