2013国赛数模C题第二问扭曲.docVIP

我们计算扭曲的方法为：计算一层的中心到某一角与其上面一层的中心到对应的该一角的向量之间的夹角作为角度，类似的两层之间可以得到八个夹角，将该八个角加起来即为总的扭曲度。 计算式为：若第一层中点为第二层中心为，第一层对的某一角的对应点为得到向量第二层对应的对的某一角的对应点为得到向量，将向量、的、平移到同一点有、，则、的扭转角即为、对应的扭转角，设为，有 内为向量乘以顺时针为正，当sin很小的时候可近似等于 ，总的扭转角即为求出剩下七个点（B、C、D、E、F、G、H）对应的、、、、、、，八个角相加即可。 Matlab示意图为： Our method to calculate the the angle of twist the way we the consider the angle which is the sum of eight angles as the criterion to evaluate the torsion resistance and each little angle is the angle between the two vectors ,one vector is from the center of one floor to one specific corner and the other vector is from the center of the floor above to the corresponding corner . The formula is that we set the center of one floor as ,the center of floor above as , is the vector from to a specific corner , is from to ,is the angle between the two vectors.,””means the vector cross product ,the value is positive when clockwise,and sin≈,when is very small and close to 0;we can calculate the total ~ the same way .The total angle is +++++++; Matlab Figure： Computed data： angle 2 3 4 5 6 7 2011 00.010857143 0.00座机电话号码 0.010714385 0.0座机电话号码 0.011893283 2009 0.010212275 0.010853803 0.00座机电话号码 0.010718923 0.0座机电话号码 0.011904419 1996 0.01座机电话号码 0.01座机电话号码 0.01座机电话号码 000.00座机电话号码 1986 0.01座机电话号码 0.01座机电话号码 0.01座机电话号码 0.01座机电话号码 0.01座机电话号码 0.00座机电话号码 angle 8 9 10 11 12 13 2011 0.0座机电话号码 0.0座机电话号码 -0.001557151 -0.00座机电话号码 -0.00座机电话号码 -0.00座机电话号码 2009 0.0座机电话号码 0.0座机电话号码 -0.001515942 -0.00座机电话号码 -0.00座机电话号码 -0 1996 0.010694323 00.01座机电话号码 -0.001687675 -0.001770576 0.00座机电话号码 1986 0.010704724 0.01座机电话号码 0.01座机电话号码 -0.001705758 -0.001718286 0.000248296 我们认为倾斜方向就是13层每层的中心拟合后的斜线与水平面的夹角。通过对中心拟合的计算和观察我们发现这四个年份的倾斜应该分为三段倾斜直线来拟合（这样基本可以保证误差最小），第一段为1-5层，第二段为6-10层，第三段为11-14（包括塔尖），即最终得到3个倾斜的斜率。 We consider the angle between the fitting curve of the centers of 13 floors and the horizontal line as the angle of the tilt .After calculating those values,we find that those f