《损伤力学基础》课程作业二(崔玮 0820020163).docVIP

《》课程： 班级 建筑工程系2班 姓名 学号 0820020163 教师 2009-12-7  目 录 1 单轴受拉随机损伤本构关系 1 2 损伤变量及本构关系的确定 3 2.1 断裂应变为均匀分布 3 2.2 断裂应变为三角形分布 4 2.3 断裂应变为抛物线分布 6 3 结论 7 作业要求： 假设微弹簧断裂应变分别服从一维均匀分布、三角形分布、抛物线分布，给出三种情况下损伤变量D与应变的关系、应力—应变关系 1 单轴受拉随机损伤本构关系 受拉损伤属于损伤的一种基本模式，对于混凝土单轴受拉应力状态，可用图1所示的细观随机断裂模型建立其随机损伤本构关系，模型中每一个微弹簧用于模拟混凝土细观受拉损伤单元。在外力作用下，微弹簧拉应力增加、伸长变形，当微弹簧达到其极限断裂应变时，该单元退出工作。由于每个微弹簧具有不同的破坏应变，弹簧退出工作时的过程可以视作随机发生的，所以各弹簧的断裂应变也应是随机物理量。 图1 受拉单元与细观随机断裂模型 图2 微弹簧应力应变关系 上述模型包含的假定主要有： (1) 在单轴受拉应力状态下，混凝土单元体中损伤的产生和演化可以用图1中所示并联弹簧模型来描述； Δ为随机变量，服从某一概率分布。 根据上述模型，定义受拉损伤变量为： (1) 式中：At为细观单元断裂而导致材料退出工作的面积；A为受拉试件横截面积。 在任一拉应变状态，由于试件均匀受拉，故存在： (2) 式中，Ai为微弹簧的面积；N为微弹簧总数量；为微弹簧的极限破坏应变，为服从某一概率分布的随机变量；为Heaviside函数： (3) 当单元中弹簧个数N趋近于无穷大，并联弹簧可以看作一连续体，根据随机积分的定义，若上述极限存在，则有： (4) 式中： 为一维随机场，反应受拉微弹簧断裂应变沿x方向的随机分布性质。 利用受拉单元静力平衡原理，可以给出单轴受拉时的随机损伤本构关系： (5) 若给定，由于断裂应变为一随机变量，且损伤的发展与断裂应变有关，所以公式5给出的仍是服从某一概率密度的随机变量，利用概率论的运算法则，可以给出均值关系和给定条件下的均方差表达式，关键是先要求出损伤变量的均值与方差。 现令，则由Heaviside函数定义可知服从0-1分布，再结合概率论的相关知识可得到的均值函数为： (6) 的自相关函数为： (7) ① 若两两弹簧间的断裂是相互独立的，则其自相关函数为： (8) ② 若两两弹簧间的断裂具有相关性，则其自相关函数为： (9) 式中：为自相关系数。 损伤变量的均值函数（即期望值）为： (10) 的方差函数为： (11) 式中：。 ① 若两两弹簧间的断裂是相互独立的，则其方差： (12) 当然也可从另外一个角度来考察这个问题： (13) ② 若两两弹簧间的断裂具有相关性，则结合上述公式7-11可得其方差为： (14) 上述为细观随机断裂模型的基本概念，下面将利用这些原理，针对微弹簧断裂应变服从不同概率分布的情况进行论述。 2 损伤变量及本构关系的确定 2.1 断裂应变为均匀分布 (1) 假设断裂应变服从一维均匀分布，且微弹簧间的断裂相互独立时，其概率密度函数，式中εb为断裂应变，如图3所示。则损伤变量的均值函数为： (15) 如前所述，由于各微弹簧的断裂是相互独立的，损伤变量的方差=0。所以此时的损伤可以看作是确定的，用均值表示为： (16) 损伤与应变的关系如图4所示，当应变取为断裂应变εb时，微弹簧断裂退出工作，对应的损伤值为1。 由公式5可得此种情况下对应的应力为： (17) σ-ε关系εb/2时，应力值达到最大为Eεb/4。 图4 均匀分布时D-ε关系σ-ε关系(2) 若考虑弹簧断裂间的相互影响，即两两弹簧间断裂应变不再独立，假设此时的自相关系数： (18) 结合公式79 可得 (19) 代入公式14可得损伤变量的方差为： (20) 此时方差不为0，从而考虑了随机性的影响。所以引入损伤变量后得到的应力是服从某一概率密度的随机变量，只需确定参数后就可描述完整损伤变量的分布。 2.2 断裂应变为三角形分布 (1) 假设断裂应变服从三角形分