下面我们给出计算平面梁单元坐标转换矩阵的程序,因为坐标转换矩阵在单元刚度矩阵转 换,等效结点力转换和计算结点内力时都要用到,因此我们定义一个函数,方便调用。 1..functionT=TransformMatrix(ie) 2..%计算单元的坐标转换矩阵(局部坐标-整体坐标) 3..%输入参数 4..%ie结点号注释 5..%返回值 6..%T从局部...
假设初始状态是:厚度Z0,溶剂浓度是C0,为了分析方便,把Z0均匀的分成了N层,每层的厚度是dZ,溶剂的初始浓度是C0 我是用用最简单的差分算法,这就会涉及到vector的长度缩短问题。 dZ应该是固定的,那么随着Z的变小,N也在变小。那么表达溶剂浓度的vector也在边小。有没有一个方程能够表达这个变化? 比如,初始时,厚度...