用彎矩分量求三次樣條插值函數(shù)的彎矩

劉遠鵬，劉光好，金少華

（1.河北工業(yè)大學(xué) 工程訓(xùn)練中心，天津 300401；2.河北工業(yè)大學(xué) 土木工程學(xué)院，天津 300401；3.河北工業(yè)大學(xué) 理學(xué)院，天津 300401）

用彎矩分量求三次樣條插值函數(shù)的彎矩

劉遠鵬1，劉光好2，金少華3

（1.河北工業(yè)大學(xué) 工程訓(xùn)練中心，天津 300401；2.河北工業(yè)大學(xué) 土木工程學(xué)院，天津 300401；3.河北工業(yè)大學(xué) 理學(xué)院，天津 300401）

根據(jù)參考彎矩圖解法原理，將相鄰節(jié)點函數(shù)差和節(jié)點函數(shù)的導(dǎo)數(shù)變換為側(cè)移分量和轉(zhuǎn)角分量．利用各區(qū)間轉(zhuǎn)角分量、側(cè)移分量與函數(shù)的二階導(dǎo)數(shù)形成彎矩-位移方程，以及連續(xù)梁支座相鄰兩跨轉(zhuǎn)角方程，結(jié)合連續(xù)梁彎矩圖的性質(zhì)分析樣條插值函數(shù)，做參考彎矩圖；由此得到三次樣條插值函數(shù)的支座彎矩．通過算例說明未知量參考彎矩圖、荷載參考彎矩圖、轉(zhuǎn)角分量圖和側(cè)移分量圖的使用方法．

特征點；轉(zhuǎn)角分量；側(cè)移分量；參考彎矩圖；三次樣條插值函數(shù)

三次樣條插值函數(shù)有廣泛應(yīng)用[1-3]，且有各種表達方法[4-5]，計算方法也多種．文獻 [4]則采用構(gòu)造三角陣的逆陣方法簡化節(jié)點等距問題．通?？梢詫⑷螛訔l插值函數(shù)比擬為連續(xù)梁的撓曲線；其二階導(dǎo)數(shù)視為連續(xù)梁的彎矩．因此三次樣條函數(shù)的二階導(dǎo)數(shù)分布就相當于無荷載作用有支座沉陷的抗彎剛度EI=1（本文采用無量綱）的連續(xù)梁彎矩圖．用力學(xué)計算這種連續(xù)梁方法很多，可以利用廣義圖乘法及其彎矩分量[4-7]分析連續(xù)梁的原理．利用連續(xù)梁彎矩圖和單跨梁彎矩圖（分別是結(jié)點位移彎矩圖和桿端位移彎矩圖．在不致混淆時的也簡稱為彎矩圖）的性質(zhì)；求彎矩圖從而獲得結(jié)點彎矩，從而產(chǎn)生插值函數(shù)．整體計算過程簡單，力學(xué)意義清楚．其計算量也得到優(yōu)化．

1 彎矩分量與連續(xù)梁結(jié)點轉(zhuǎn)角方程

首先分析節(jié)點區(qū)間內(nèi)的插值函數(shù)（這相當于單跨梁的撓曲線）．如下常見表達式：

根據(jù)結(jié)構(gòu)力學(xué)常用方式，對三次樣條函數(shù)的連續(xù)梁稱弦轉(zhuǎn)角為：

左右兩端的轉(zhuǎn)角分量[5-6]

統(tǒng)稱這3個量為梁彎矩的位移分量，即彎矩分量．稱梁內(nèi)1/3點和2/3點分別為左特征點和右特征點[5-6]；設(shè)為L點和R點．利用式 (1)計算特征點的彎矩，與位移分量比較：

這2個方程是梁的彎矩-位移方程，是梁的普遍形式的簡化[5]．

因連續(xù)梁相鄰兩跨的桿端轉(zhuǎn)角分量表達式不同，但是桿端連接有相同轉(zhuǎn)角，則：

這個方程為連續(xù)梁的結(jié)點轉(zhuǎn)角方程[7]．

根據(jù)樣條函數(shù)的二階導(dǎo)數(shù)，相鄰兩跨的特征點彎矩Mk1+1/3、Mk1+2/3與桿端的彎矩Mk1、Mk有如下的關(guān)系：

由此可知，彎矩分量有2種，它們是特征點的彎矩分解．若梁的弦轉(zhuǎn)角和桿端轉(zhuǎn)角為正方向，那么彎矩分量可用圖1方法表示．

連續(xù)梁彎矩圖的性質(zhì)：

1）區(qū)間內(nèi)2個特征點的側(cè)移分量數(shù)值相同，分別在桿的上下兩側(cè)，與弦轉(zhuǎn)角方向一致．如圖1a）．

圖1 特征點上的彎矩關(guān)系Fig.1 Relationship ofmomenton characteristic point

2）轉(zhuǎn)角分量的方向分別與其所表示的轉(zhuǎn)角方向一致．如圖1b）．

3）在每個區(qū)間內(nèi)，彎矩圖是直線，如圖1c）．

4）特征點的彎矩及其分量滿足梁的彎矩-位移方程和連續(xù)梁結(jié)點轉(zhuǎn)角方程．連續(xù)梁相鄰兩跨的轉(zhuǎn)角分，如圖1d）．

只要區(qū)間內(nèi)插值函數(shù)是三次多項式；無論采用那種表達式等，上述性質(zhì)都成立．而分析樣條函數(shù)過程中，對結(jié)點轉(zhuǎn)角的估計相當于對特征點的轉(zhuǎn)角分量的估計．

2 參考彎矩圖分解與彎矩分量圖

2.1 參考彎矩圖分解

計算過程中連續(xù)梁彎矩圖是用已知量和未知量共同表示的，則為區(qū)別稱其為參考彎矩圖（用M表示）．那么據(jù)線性關(guān)系，將參考彎矩圖M分解為[6]：

荷載參考彎矩圖是利用連續(xù)梁的支座位移和一個已知條件（一端的已知彎矩或轉(zhuǎn)角）計算形成；且表達式不含未知量．而未知量參考彎矩圖根據(jù)參考彎矩圖的未知量計算形成的彎矩圖；不含已知條件（位移條件或桿端彎矩條件）．其各彎矩表達式是x的一次齊次式．

荷載參考彎矩圖和未知量參考彎矩圖都是桿內(nèi)無荷載的連續(xù)彎矩圖，則在每個彎矩圖內(nèi)分別具有連續(xù)梁彎矩圖的上述4條性質(zhì)．對于具體的計算過程來說，的值由已知的邊界條件和0到n的xk，yk的值確定；這樣對于來說yk的值全為0，則也全為0，帶入式 (6)、式 (7)后得，，根據(jù)公式后得．

根據(jù)三次樣條插值函數(shù)連續(xù)梁彎矩圖的上述性質(zhì)，通過作參考彎矩圖方法求支座彎矩．將支座彎矩代入公式 (1)就生成區(qū)間內(nèi)三次插值樣條函數(shù)．

2.2 彎矩分量圖與算例

已知函數(shù)值分布如表1所示[8]．

表1 節(jié)點數(shù)值Tab.1 Digiton point

根據(jù)支座位移計算側(cè)移分量：

將它們標在連續(xù)梁上形成側(cè)移分量圖，如圖2a）．

提取圖2a）和圖2b）的位移分量和轉(zhuǎn)角分量，根據(jù)式 (6)計算第1跨的L點彎矩：

作荷載參考彎矩圖．將1跨 L點彎矩標在連續(xù)梁上，如圖3a）．因最左端彎矩等于0，根據(jù)式 (9)、式 10，作第1跨參考彎矩圖直線，產(chǎn)生R特征點和右端的彎矩，如圖3a）．

根據(jù)式 (7)轉(zhuǎn)角分量

將其標在轉(zhuǎn)角分量圖上，如圖2b）．

利用結(jié)點轉(zhuǎn)角方程計算第2跨左端轉(zhuǎn)角分量，如圖2b）．

同理依次分析第2跨和第3跨的桿端彎矩和特征點的彎矩；依次產(chǎn)生各跨參考彎矩圖直線．同時得到轉(zhuǎn)角分量的分布．如圖2b）（稱其為轉(zhuǎn)角分量圖）．因此荷載參考彎矩圖如圖3a）所示．

根據(jù)未知量假設(shè)作未知量參考彎矩圖．

這個圖的A端轉(zhuǎn)角等于0，且無側(cè)移；則各個特征點彎矩等于轉(zhuǎn)角分量．

利用A端彎矩等于未知量x，A端轉(zhuǎn)角分量等于0作第1跨彎矩圖，如圖3b）．

利用特征點未知量參考彎矩與桿長成反比關(guān)系，則

圖2 彎矩分量的分布圖Fig.2 Distribution ofmomentcomponent

圖3 參考彎矩圖求三次樣條插值函數(shù)Fig.3 Reference BMD for cubic spine interpolation function

依次作第2跨、第3跨和第4跨的未知量參考彎矩圖，如圖3b）．

荷載參考彎矩圖與未知量參考彎矩圖相加形成參考彎矩圖．因第4跨的右轉(zhuǎn)角分量和側(cè)移分量是已知量，其和是R特征點彎矩；與參考彎矩圖的特征點彎矩比較表達式建立方程：

因此x=-2.028 6，計算各桿端的彎矩值，用式 (11)計算各桿端的支座彎矩：

最后結(jié)果與其它方法一致．

具體計算過程中，由于各個分量均有系數(shù)2，可以先將系數(shù)2提出，最后得到的數(shù)值乘以2以獲得各段點彎矩．這樣可以進一步簡化計算過程．

如果左端邊界條件與此例相同，那么所需未知量以及建立的參考彎矩圖的形式可不變；而依據(jù)右端邊界條件建立不同的方程．對于第2種邊界條件和第3種邊界條件也可以用類似的方法求解．

2.3 算例的MATLAB程序

S為邊界條件，com_la為側(cè)移分量，com_l為左端轉(zhuǎn)角分量，com_r為右端轉(zhuǎn)角分量，MF為荷載彎矩（l為左端，r為右端），Mx為未知量彎矩（l為左端，r為右端），M為每個結(jié)點的彎矩．

3 結(jié)論

本文利用連續(xù)梁參考彎矩圖解法形成分析計算三次插值樣條函數(shù)的連續(xù)梁結(jié)點彎矩方法．其中利用參考彎矩圖的兩個層次分解：荷載 未知量圖分解，以及彎矩位移分量圖分解．利用了梁的特征點的彎矩的側(cè)移分量、轉(zhuǎn)角分量，以及彎矩 位移方程和連續(xù)梁結(jié)點轉(zhuǎn)角方程．通過未知量參考彎矩圖、荷載參考彎矩圖、轉(zhuǎn)角分量圖和側(cè)移分量圖分析產(chǎn)生未知量方程，進一步得到三次插值樣條函數(shù)的連續(xù)梁結(jié)點彎矩．計算過程僅需一個未知量，具有明顯的并行計算特點．

[1]陳盛，趙東標，陸永華，等．基于二維輪廓點云的螺紋中徑計算 [J]．光學(xué)精密工程，2015，23（6）：1791-1799．

[2]常景娜，高慧斌，喬冠宇．一種基于傳感器的測量目標運動參數(shù)方法研究 [J]．傳感器與微系統(tǒng)，2015，34（11）：19-22．

[3]梁棟，劉良棟，何英姿．月球精確軟著陸最優(yōu)標稱軌跡在軌制導(dǎo)方法 [J]．中國空間科學(xué)技術(shù)，2011（6）：27-35．

[4]Moon BS．Anexplicitsolution for thecubic spline interpolation for functionofasinglevariable[J]．Applied Mathematicsand Computation，2001，117（2-3）：251-255．

[5]劉光好，趙光鑒．剛架參考彎矩圖解法的側(cè)移問題 [C]//崔京浩．第十一屆全國結(jié)構(gòu)工程學(xué)術(shù)會議論文集（VolⅠ）．長沙，2002：395-398．

[6]劉光好，張揚，王俊杰，等．連續(xù)梁支座反力、彎矩、位移及三支座反力方程 [J]．公路交通科技（應(yīng)用技術(shù)版），2007，31（3）：106-108．

[7]Liu Guanghao，Yang Shaobin，LiuYuanpeng．Three-columniationendmomentequationof framewithsupportsinkageand angle[C]//2011INternationA l Conferenceon Remote Sensing，Environmentand Transportation Engineering（Vol5）．Nanjing,2011：4670-4682．

[8]朱曉臨．數(shù)值分析 [M]．合肥：中國科學(xué)技術(shù)大學(xué)出版社，2010：148-150．

[責(zé)任編輯 楊 屹]

Themomentcomponent formomentof cubic spline interpolation function

LIU Yuanpeng1，LIU Guanghao2，JIN Shaohua3

(1.Engineering Training Center,HebeiUniversity of Technology,Tianjin 300401,China;2.SchoolofCivilEngineering and Transportation,HebeiUniversity of Technology,Tianjin 300401,China;3.Schoolof Sciences,HebeiUniversity of Technology,Tianjin 300401,China)

According to theory of referencebendingmoment,the difference of function atpointbetween neighbor points and derivative of function atnodeare translated into sideway componentand angle component.The sideway component, angle componentand the second derivativeof the function assemblemoment-deflection formula in each section.Themoment-deflection formulaand the cornerequation between adjacentspan on continuousbeam combinedw ith referencebendingmomentdiagram are used for solving Cubic Spline Interpolation function.The example in article illustrates using methodsof the load referencebendingmomentdiagram,theunknown referencebendingmoment,thesideway component and the angle componentdiagram for solving.

characteristic point;angle component;sideway component;reference bendingmomentdiagram;cubic spline interpolation function

TU311.4

A

1007-2373(2016)04-0093-06

10.14081/j.cnki.hgdxb.2016.04.015

2015-11-19

河北省高等學(xué)校科學(xué)技術(shù)研究重點項目（ZD2014051）

劉遠鵬（1984-），男（漢族），實驗師．

：劉光好（1948-），男（漢族），副教授．