The Application of Differential Quadrature Method in Structural Dynamic Analysis
-
摘要: 微分求积法(DQM)是1种求解微分方程初(边)值问题的数值方法,通常以较小的计算工作量即可获得较高的数值精度。这种方法应用于工程领域时多用来解决梁、板等结构的静力分析或结构特征值分析等问题,即对边值问题的微分方程的求解。结构动力分析属于初值问题,荷载和结构反应都具有特殊性,直接套用DQM求解边值问题并不能获得问题的解。本文尝试利用微分求积原理建立求解结构动力反应的具体方法。借鉴单元法的思想,将荷载持时划分为若干个时步,在每个时步内对动态荷载和结构反应进行离散,然后用DQM对时步逐个进行求解,得到体系在整个时域内的反应过程。通过对3种不同自振周期的线弹性单自由度体系在不同频率简谐激励下反应的计算,阐释了本文方法的可行性以及高精度、高效率的特点,通过数值试验确定了时步内相对较优的节点数,并为时步长度的选取提供了建议。Abstract: The differential quadrature method (DQM) is a numerical technique of solving initial/boundary value problems of differential equations and capable of obtaining a higher numerical accuracy with a smaller calculation workload. This method is often used to solve the problems of structural static analysis of beams and slabs or eigenvalue analysis when it is applied to the engineering fields, which is to solve the differential equation of the boundary value problems. Dynamic analysis of structures is an initial value problem as well as particular loads and structural response. As a result, applying the DQ method of solving the boundary value problem directly cannot obtain solution of problem. The principle of differential quadrature is applied to establish the specific method of solving structural dynamic response in this paper. By analogizing the idea of unit method, the duration of the load is divided into many time steps and dynamic load and structure are discretized in each time step, then the response of the system can be solved in the whole time domain by employing the DQ method step by step. The feasibility of this method and the characteristics of high precision and high efficiency are expounded by calculating the response of three linear elastic single-freedom-degree system of different natural vibration periods excited by simple harmonic loads of different frequencies. By means of numerical experiment, the optimal meshing scheme is determined and the suggestion for the time step is given.
-
引言
龙门山断裂带位于青藏高原东缘,作为青藏高原向东扩展的边界,其地势险峻,自然灾害频发,一直是国内外研究的热点区域。2008年5月12日和2013年4月20日,在龙门山构造带上发生了汶川MS 8.0和芦山MS 7.0大地震,进一步说明该区域构造活动极为强烈。汶川地震发震断裂为北川-映秀断裂和灌县-江油断裂,分别产生了约240km和72km的地表破裂带。地震前后,各学者对龙门山地区都进行了一定的研究(李传友等,2004;李勇等,2006;周荣军等,2006;Densmore等,2007;刘静等,2008;徐锡伟等,2008),基本认为北川-映秀断裂是逆冲为主兼具走滑的断裂,具有明显分段性;然而区域内人为因素和气候对地形改造迅速,地貌标志保存困难且植被茂密,影响滑动位错和地层年代的定量观测,造成龙门山断裂带运动学特征定量研究的不足。构造活动对区域地貌演化研究具有很大的影响,地貌参数值对几千年至几十万年尺度的构造活动有指示作用(Hack,1973;Kirby等,2012),运用高精度数字高程模型(DEM)与GIS等技术,通过对地貌参数的定量研究,可以获取区域构造活动的信息(何祥丽等,2014;常直杨,2015)。
随着GIS技术的不断发展和DEM数据的获取,国内外学者利用各种地貌参数分析构造活动,这些参数包括了坡度、地形起伏度、河长坡降指数、Hack剖面、流域对称度以及面积高程积分等,定量化的地貌参数可以多维度、全方位地描述与分析构造地貌,并获取构造活动的多层次信息(张会平等,2006a;常直杨,2015)。坡度与地形起伏度是传统描述地貌形态的参数。面积高程积分反映了流域盆地的地貌演化过程,其对构造活动的影响十分敏感。
本研究以北川-映秀断裂及邻区流域为研究对象,区域内河流发育,河流地貌对构造活动的响应敏感,利用ASTER GDEM V2数据,空间分辨率30m,基于ArcMap 10.1水文分析模块对该地区坡度、地形起伏度及面积高程积分进行提取,进而定量研究区域的地貌特征,北川-映秀断裂沿走向运动学特征发生变化,本文依据地貌参数,定量化了新构造活动程度,讨论了北川-映秀断裂走向上及其两侧的参数变化特征、地层岩性和水系的变化,探讨了其与构造活动的关系。
1. 区域概况
龙门山构造带由后山断裂(汶川-茂县)、中央断裂(映秀-北川)、前山断裂(灌县-江油)和山前隐伏断裂等4条逆冲断裂及其所夹持的逆冲楔体组成,其形成过程表现为自后山断裂向东南方向发展的前展式逆冲,受到北川附近的NW向走滑断层及都江堰附近的NW向三江口断裂的切割,龙门山构造带可分为北段、中段和南段。研究区域(图 1)位于龙门山的中北段地区,断裂带沿线分布有新生代沉积物,断裂横穿湔江、沱江、凯江及涪江流域支流安昌河、通口河以及平通河,穿过断裂水系发生大型右旋位错,水系形态受北川-映秀断裂逆冲走滑作用明显。
龙门山构造带表现为逆冲推覆模式,逆冲断层下盘产生新的逆冲活动并不断向前陆方向扩展(林茂炳等,1991),区域内自北西—南东方向地层出露由老变新,中段出露大面积元古代岩浆岩,北段出露大面积志留系板岩、凝灰岩,泥盆系泥岩、灰岩、白云岩,向前陆方向,三叠系砂岩、碳酸盐岩和侏罗系砂岩沿断裂走向呈带状展布。北川-映秀断裂处于龙门山中央断裂中段,其南段展布于彭灌杂岩的东南缘,断层两侧地层岩性差异明显,其西侧主要出露彭灌杂岩、下震旦统火山岩,发育一系列飞来峰构造;其东侧为志留系泥质岩和碳酸盐岩建造,泥盆系碳酸盐岩与碎屑岩互层以及上三叠统的须家河组,西侧地层较东侧地层强硬(姚琪等,2012);北川以北段展布于志留系与泥盆系地层中,表现为变质的下古生界逆冲于未变质的古生界之上(陈社发等,1994)。地层单元上,北川-映秀断裂北西侧与南东侧地层存在明显差异。
2. 数据方法
本研究利用ASTER GDEM V2数据(水平分辨率为30m),流域水系和地貌参数的提取分析均在ArcGIS 10.1平台下完成,地理坐标系统为WGS84。
2.1 坡度与地形起伏度
坡度与地形起伏度是传统的宏观地形分析方法(Scherler等,2014)。坡度指地面的倾斜程度,反应地表物质剥蚀和堆积的发展,能够指示地貌发育阶段及程度。坡度具有一定的临界值,约为30°,超过这个值,则不再随侵蚀速率等因素的增大而增大(Schmidt等,1995)。使用ArcGIS 10.1空间分析中的坡度分析自动生成研究区域坡度图。
地形起伏度为一定区域内最高点和最低点高程之差,它能够反映地面侵蚀程度(Zhang等,2011),间接反映了构造活动的强弱。一般使用公式:
$$ R=H\text{max}-H\text{min} $$ (1) 其中,R表示地形起伏度,Hmax和Hmin分别表示分析窗口内的最大高程值和最小高程值。利用ArcGIS10.1的空间分析模块中的焦点统计工具,采用分析窗口4.28km2(23×23)统计出最大高程值和最小高程值(南希等,2017),再使用栅格计算器计算地形起伏度。
2.2 条带状剖面
条带状剖面可统计一定区域范围内地形高程的最大值、最小值和平均值,用以半定量、定量分析研究区域内山峰、河谷的高程变化以及侵蚀程度(张会平等,2006b;梁明剑等,2014;苏琦等,2016)。相对于线状剖面,条带剖面则可以较为准确地表征地球表面本身具有的面状地形起伏特征。为了比较横向地形变化,利用ArcGIS 10.1沿断裂走向,在3个不同区域沿约130°方向提取5km宽条带状剖面,将高程数据导入Origin 8.1成图。
2.3 面积-高程积分(HI)
面积高程积分是Strahler(1952)提出的定量描述流域内高程的分布、表征未被侵蚀部分体积的1个参数,可以指示区域构造活动差异(Alipoor等,2011;Gao等,2013)。
在1个流域盆地演化的不同阶段,其地貌形态会发生较大的变化。在演化早期阶段,河流侵蚀能力强,地貌以深切河谷为主,流域内高海拔区域较大;随着河流演化进入后期,U形谷发育,流域内低海拔区域的面积逐渐增大。面积-高程积分曲线就是1个表征流域盆地内某一特定高程值之上的面积占全流域面积的比例函数(常直杨等,2014)。
在1个流域盆地内,最高海拔为Hmax,最低海拔为Hmin(一般为河流出口处的海拔值),取H=Hmax - Hmin,h为0—H之间的任意值,a为H—h之间的面积,A为全流域面积(图 2(a))。令x=a/A、y=h/H作图,可得到面积高程积分曲线(图 2(b))。
图 2 面积-高程积分曲线示意图(Singh等,2008)Figure 2. Schematic diagram of area-elevation integral curve (fromh Singh et. al., 2008)图 2(b)中,灰色部分表示流域内未被侵蚀部分的比例,即面积高程积分(HI),是1个表征侵蚀程度的值。其计算公式为:
$$ HI=\int_{0}^{1}{\left(\frac{h}{H} \right)\text{d}\left(\frac{a}{A} \right)}=\int_{0}^{1}{y\text{d}x} $$ (2) 根据 Pike等(1971)的研究,可以简化为:
$$ HI=\frac{H\text{max}-H\text{mean}}{H\text{max}-H\text{min}} $$ (3) 其中,Hmax为流域高程最高值,Hmin为最小值,Hmean为平均值。
由图 2(b)可知,流域积分曲线为上凸型时,其积分值高,流域积分曲线为下凹型时,积分值低。高积分值意味着流域内大部分物质体积未被侵蚀,地形演化时间短,地貌处于幼年期。积分曲线呈上凸型与流域内岩性、构造活动和降雨剥蚀等因素有关。地形演化时间越长,地貌处于老年期,积分曲线呈下凹型。介于凸型及凹型之间,呈凹凸型的则为中年期。依据Strahler(1952)对HI的划分,可分为3级:①幼年期:HI>0.5;②中年期:0.4≤HI≤0.5;③老年期:HI<0.4。面积高程积分对构造活动的响应非常敏感,可很好地判断流域盆地内地貌特征及其构造活动响应。
3. 结果分析
面积高程积分的提取在ArcGIS 10.1水文分析模块进行,为了细致分析流域地貌发育特征对于构造活动的差异响应,结合龙门山中北段流域特征,选取阈值7000提取河网,获得1381个亚流域盆地,分别计算各亚流域的HI值,为了更好地反映HI的空间分布特征,对HI值进行克里金插值分析,获取等值分区图。
3.1 宏观地貌分析
龙门山中北段沿断裂横跨断裂两侧,坡度及地形起伏度表现出明显的西高东低的特征。地貌参数沿断裂走向呈现中段高、北段低的特征,且坡度与地形起伏度分布特征具有一致性。龙门山中段坡度与地形起伏度极大,沿北川-映秀断裂两侧参数值差异明显,北西侧坡度达到80°,大部分地区坡度处在临界值之上,跨过断裂坡度值骤减,普遍在40°左右。地形起伏度有相似的递变趋势,北川-映秀断裂北西侧局部地形起伏达到5000m,跨断裂后骤减。龙门山北段坡度与地形起伏值相比中段减少,跨断裂自西向东两侧参数值仍有明显的骤减,大部分地区坡度值处在临界值30°之下,地形起伏在2000—3000m之间,局部达到4000m(图 3)。
3条条带状剖面反映了沿断裂走向上地形的变化(图 4),图中(a)、(b)、(c)分别代表图 1中AA′、BB′、CC′条带状剖面。龙门山中段平均海拔处于3000m左右,在北川附近平均海拔为2000m左右,再往北段降低至1500m;垂直断裂方向,AA′剖面显示龙门山中段地形起伏相对较大,3条剖面显示地形落差最大处在北川-映秀断裂两侧。
3.2 流域特征分析
龙门山中北段发育多条大型水系,断裂横穿湔江、沱江、凯江及涪江流域支流安昌河、通口河以及平通河,北段水系右旋位错明显(贾营营等,2010)。在流域盆地基础上提取亚流域盆地获取HI值,由图 5可见,HI值呈规律性复杂分布。HI高值(大于0.5)基本出现在断层上盘,流域发育程度低,处于幼年期;前山断裂的下盘HI值(小于0.4)较低,流域发育成熟,处于老年期。龙门山中段HI高值出现在逆冲断裂的上盘且集中分布,在茂汶-安县区段内,HI高值分布靠近前缘且较分散,北段断裂两侧的HI值差异较中段小,沿断裂方向,局部的河谷地区出现HI低值区,灌县-江油断裂上盘出现HI高值区。
4. 讨论
地表地形地貌特征主要受控于构造、气候以及岩性等综合性因素。近年来大量的研究认为,在构造非常活跃的地区,断层作用控制的山体隆升是控制区域地貌生长的主要因素(Kirby等,2003;Zhang等,2006;Kirby等,2012)。而在地貌演化中,前人的研究也指出河流地貌参数与气候变化和基岩的抗侵蚀能力有关(赵洪壮等,2010;邵崇建等,2015)。本研究结合龙门山中北段宏观地貌特征与地貌参数,探讨区域地形地貌主控因素及构造活动差异性。
4.1 地貌参数影响因素分析
气候是地貌演化过程中的重要因素,其直接表现形式之一为降水,降水通过影响水系流量进而影响河流对于河道的侵蚀堆积作用,从而控制区域地形的演化(Hartshorn等,2002)。在年均降水量约400mm/a的祁连山地区,前人对降水量与亚流域盆地HI指数做了相关性分析,认为降水量对地貌指数无明显控制(苏琦等, 2016, 2017)。同时HI值具有空间依赖性(郑光佑,2002),在年降水量更高的龙门山地区,邵崇建等(2015)将集水盆地依据高程划分山地、丘陵和平原3种类型,根据同类型地形下的降水与HI指数对比,得出降水不是影响地貌指数的主要因素。
本文利用前人插值获取的全球高分辨降水数据(Hijmans等,2005),获得龙门山地区年平均降水分布图(图 6)。年平均降水沿东南到西北方向递减,900mm/a降水等值线基本与龙门山断裂带平行,研究区域内北川-映秀断裂两侧年均降水量普遍在850—950mm/a之间,局部存在地形雨造成的降水高值区,但是顺降水等值线方向研究区域HI指数变化极大。龙门山地区降水量受地形控制,形成季节性强降雨带,研究区域属于山地与丘陵地形,受到地形雨影响严重(邵崇建等,2015),在不考虑构造活动的前提下,流域内年降水量越多,内部侵蚀越严重,面积-高程积分越小。由图 6可知,研究区域内北川-映秀断裂两侧降水量普遍在800—900mm/a之间,变化不大,然而研究区域所在的山前丘陵地区与山地地区HI指数差异很大。据此认为区域年降水量对地貌指数并无明显控制作用。
岩性也是控制地貌演化的1个因素。不同岩性的抗侵蚀能力不同,坚硬岩抗侵蚀能力强,风化程度小,易于保留原地地貌,软弱岩抗侵蚀能力弱,风化程度大,易被塑造成新的地貌形态。前人根据岩性及其形成时代对比讨论地形发育情况,认为岩性不是影响HI值的主要因素(苏琦等,2017)。陈彦杰等(2005)以不同流域阈值提取集水盆地地层和活动构造做叠图分析时,发现阈值小于1000m2,HI值受岩性与构造影响;阈值大于2000m2,HI值主要受构造活动影响。在龙门山地区地貌参数研究中,梁欧博等(2018)将岩石依据坚硬程度分为坚硬岩、较硬岩、较软岩、软岩和极软岩5个等级,通过同类型岩性下地貌指数的对比得出岩性并非影响地貌参数的主控因素。
研究区域除龙门山中段出露前寒武系火成岩外,区域内发育志留系碳酸盐岩和砂岩,处于较硬岩与较软岩之间,且跨北川-映秀断裂岩性界限分明。然而在北川-映秀断裂北川和南坝的局部河谷,不同的岩性有相似的HI值。因此岩性与气候并非研究区域内地貌特征的主控因素。
印-亚板块碰撞塑造了青藏高原及其周缘的地质地貌格局,沿龙门山断裂系发生了一系列挤压与右旋变形,形成了其东缘快速变化的地形特征(Burchiel等,2008)。在龙门山地区的地貌学研究中,Kirby等(2003,2008)通过排除岩性、沉积通量及地形雨的影响,认为区域内河道陡峭指数梯度变化反映了构造抬升活动的差异。前人在祁连山及龙门山地区对流域盆地的地貌演化研究中,通过影响因素的相干性分析,也认为HI值的变化反映了构造隆升的差异(梁明剑等, 2013, 2014;李奋生等,2015;苏琦等,2017)。根据上述对岩性及降水量的分析,推测构造活动是区域地形地貌演化的主控因素。
4.2 龙门山中北段构造活动差异性
认为龙门山断裂带第四纪构造活动存在差异是1个普遍的认识。龙门山南段表现为强烈的逆冲活动,南段前陆区发育广阔褶皱构造带,延续至中段龙泉山一带并向北逐渐消亡;中北段构造变形主要发生在山区冲断带,而中段成为构造活动的转换区间,整体上呈现逆冲兼走滑的性质(杨晓平等,1999;马保起等,2005;姜大伟等,2018)。在龙门山地区活动构造研究中,马保起等(2005)根据岷江阶地的变形估算了汶川-茂汶断裂、北川-映秀断裂和灌县-江油断裂晚第四纪平均垂直滑动速率为0.5mm/a、0.3—0.6mm/a和0.2mm/a。前人对整个东缘的断裂地貌进行了研究,在茂汶和映秀两地得到汶川-茂汶断裂和北川-映秀断裂平均垂直滑动速率约为0.84mm/a和0.54mm/a,平均水平滑动速率约为1mm/a,灌县-江油断裂则表现出纯逆冲特征,结合对北段青川断裂走滑运动的研究,得出北段的右旋走滑位错量大于中段,且承担龙门山断裂带大部分右旋走滑分量(马保起等,2005;周荣军等,2006;李勇等,2006;Densmore等,2007;樊春等,2008;徐锡伟等,2008)。
坡度与地形起伏度是地表形态的直观反映,研究区域地形起伏与坡度向盆地内梯度递减,沿断裂走向中段比北段地形起伏大,坡度大,表明了龙门山中段构造抬升活动强烈;北段的地貌参数值均变小且跨断裂没有快速变化,表明了构造抬升活动减弱,这种沿断裂走向与垂直走向上的地貌参数差异与断裂的构造性质及滑动速率值的变化呈现一致性的特点。
龙门山地区大量的低温热年代学研究反映了区域剥蚀速率的差异,根据动态均衡原理,假定地形演化处在1个稳定模式下(England等,1990;Willett等,2002),山地基岩构造抬升速率与剥蚀速率相等。前人在龙门山中段的低温热年代学研究中得到剥蚀速率的空间展布模式,从四川盆地西缘到北川-映秀与汶川-茂汶断裂之间的区域,剥蚀速率从小于0.1km/m.y增长到大于1km/m.y,而后再向高原内部递减(Tian等,2013)。面积-高程积分对区域应力挤压造成的构造隆升反映敏感(郑光佑,2002),龙门山中段逆冲推覆构造形成的构造推覆体表现异常高的HI值,向盆地方向,HI值迅速降低;再向北东方向,断裂上盘HI值相对降低并且更集中在断裂附近,茂汶-北川一带HI值呈过渡状态,HI高值区靠近断裂且向北减弱,前山断裂由于纯逆冲性质,其上盘呈现明显HI高值,北川-南坝一带HI值相对减小,局部的河谷盆地呈现HI低值区。HI高值区表明了构造隆升的强烈,区域内HI值的分布模式也与剥蚀速率的变化呈现一致性。 Tan等(2019)在对龙门山中北段的低温热年代的研究中得到岷山地区剥蚀速率约0.9km/m.y,与龙门山中南段剥蚀速率进行对比,提出了高原东缘最大剥蚀速率带,认为该带是高原向东逆冲活动的边界,结合前人在东缘的构造地貌研究以及上述结果,反映了地貌参数与剥蚀速率和构造抬升变化的一致性,同时龙门山北段构造活动的变化以及所引起的地形地貌因素的响应,表明北段不再是高原东向扩展的边界(Zhang等,2011)。
5. 结论
本文通过ASTER GDEM V2数据系统提取了龙门山中北段地区坡度、地形起伏度、条带状剖面和流域水系及其面积-高程积分等地貌参数,并分析其地质意义,研究结论如下:
(1)构造活动是龙门山中北段地区地貌特征的主控因素,地貌地形特点反映了龙门山断裂带的活动特征。
(2)坡度、地形起伏度和条带状剖面等地貌参数反映了青藏高原东缘龙门山断裂带的逆冲演化过程,地貌参数值呈梯度分布,走向上存在差异;HI值分布图反应构造抬升活动的演化,由龙门山中段的逆冲为主转换为北段逆冲兼走滑的特征,反应龙门山中北段的构造差异性。
致谢: 感谢姜大伟老师在论文写作及修改过程中给予的讨论与意见。 -
表 1 体系的基本特性
Table 1. Basic characteristics of systems
体系编号 自振周期Tn/s 阻尼比ξ 1 0.68 0.05 2 0.25 0.05 3 0.08 0.05 表 2 简谐荷载的信息
Table 2. Information of simple harmonic load
荷载函数/N 周期/s 振幅/m 持时/s sin2πt 1 1 40 sin10πt 0.2 1 40 sin20πt 0.1 1 40 表 3 荷载周期1s的平均相对误差
Table 3. Average relative error with load period of 1s
m 体系编号 1 2 3 位移/% 速度/% 位移/% 速度/% 位移/% 速度/% 2 98.734 98.734 98.734 98.734 98.734 98.734 4 88.218 81.849 99.415 389.850 53.532 1293.500 6 42.402 32.547 15.926 4.728 10.979 9.609 8 6.698 8.020 1.543 7.814 2.042 2.741 10 1.623 2.687 0.907 2.338 1.036 0.349 12 0.121 0.138 0.897 3.341 0.819 4.156 14 0.008 0.009 0.935 5.153 0.716 1.290 16 0 0 0.930 3.899 0.635 0.825 18 0 0 1.005 3.292 0.577 2.28 20 0 0 1.268 4.715 0.527 0.899 表 4 荷载周期0.2s的平均相对误差
Table 4. Average relative error with load period of 0.2s
m 体系编号 1 2 3 位移/% 速度/% 位移/% 速度/% 位移/% 速度/% 2 99.749 99.749 99.749 99.749 99.749 99.749 4 799.760 2460.400 403.870 411.780 108.540 225.480 6 90.884 11.229 23.620 19.707 11.436 5.439 8 6.052 15.762 2.102 2.135 1.068 2.520 10 0.264 0.046 0.116 0.140 0.765 1.155 12 0.009 0.022 0.005 0.004 0.807 3.325 14 3681.600 256.240 251.760 2935.800 972.740 1334.400 16 0 0 0 0 0.093 0.167 18 0 0 0 0 0.010 0.013 20 0 0 0 0 0.001 0.001 表 5 荷载周期为0.1s的平均相对误差
Table 5. Average relative error with load period of 0.1 seconds
m 体系编号 1 2 3 位移/% 速度/% 位移/% 速度/% 位移/% 速度/% 2 99.875 99.875 99.875 99.875 99.875 99.875 4 1620.200 1370.200 578.690 459.540 51.868 52.298 6 175.570 9.444 65.345 7.903 6.80×10282 8.63×10282 8 11.571 9.569 4.413 3.254 1.951 2.655 10 0.498 0.046 0.204 0.044 0.142 0.266 12 0.016 0.013 0.007 0.005 0.006 0.007 14 27683 306.960 1322.400 246.430 234.600 3888.600 16 0 0 0 0 0 0 18 0 0 0 0 0 0 20 0 0 0 0 0 0 -
李鸿晶, 王通, 2011a.结构地震反应分析的逐步微分积分方法.力学学报, 43(2):430-435. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=CAS201303040000629782 李鸿晶, 廖旭, 王通, 2011b.结构地震反应DQ解法的两种数值格式.应用基础与工程科学学报, 19(5):758-766. http://d.old.wanfangdata.com.cn/Periodical/yyjcygckxxb201105008 廖旭, 李鸿晶, 孙广俊, 2013.基于DQ原理的结构弹塑性地震反应分析.工程力学, 30(7):161-166. http://www.cnki.com.cn/Article/CJFDTotal-GCLX201307026.htm 王鑫伟, 1995.微分求积法在结构力学中的应用.力学进展, 25(2):232-240. doi: 10.3321/j.issn:1000-0992.1995.02.010 Bellman R., Casti J., 1971. Differential quadrature and long-term integration. Journal of Mathematical Analysis and Applications, 34 (2):235-238. doi: 10.1016/0022-247X(71)90110-7 Bellman R., Kashef B. G., Casti J., 1972. Differential quadrature:a technique for the rapid solution of nonlinear partial differential equations. Journal of Computational Physics, 10 (1):40-52. doi: 10.1016/0021-9991(72)90089-7 Bert C. W., Jang S. K., Striz A. G., 1988. Two new approximate methods for analyzing free vibration of structural components. AIAA Journal, 26 (5):612-618. doi: 10.2514/3.9941 Bert C. W., Wang X. W., Striz A. G., 1993. Differential quadrature for static and free vibration analyses of anisotropic plates. International Journal of Solids and Structures, 30 (13):1737-1744. doi: 10.1016/0020-7683(93)90230-5 Bert C. W., Malik M., 1996. The differential quadrature method for irregular domains and application to plate vibration. International Journal of Mechanical Sciences, 38 (6):589-606. doi: 10.1016/S0020-7403(96)80003-8 Bert C. W., Malik M., 1997. Differential quadrature:a powerful new technique for analysis of composite structures. Composite Structures, 39 (3-4):179-189. doi: 10.1016/S0263-8223(97)00112-8 Fung T. C., 2001a. Solving initial value problems by differential quadrature method——part 1:first-order equations. International Journal for Numerical Methods in Engineering, 50 (6):1411-1427. doi: 10.1002/(ISSN)1097-0207 Fung T. C., 2001b. Solving initial value problems by differential quadrature method——part 2:second-and higher-order equations. International Journal for Numerical Methods in Engineering, 50 (6):1429-1454. doi: 10.1002/(ISSN)1097-0207 Jang S. K., Bert C. W., Striz A. G., 1989. Application of differential quadrature to static analysis of structural components. International Journal for Numerical Methods in Engineering, 28 (3):561-577. doi: 10.1002/(ISSN)1097-0207 Kang K., Bert C. W., Striz A. G., 1995. Vibration analysis of shear deformable circular arches by the differential quadrature method. Journal of Sound and Vibration, 183 (2):353-360. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=cd5d834cd6a527d56839a5876ae0d87a Kang K., Bert C. W., Striz A. G., 1996. Vibration analysis of horizontally curved beams with warping using DQM. Journal of Structural Engineering, 122 (6):657-662. doi: 10.1061/(ASCE)0733-9445(1996)122:6(657) Liew K. M., Han J. B., Xiao Z. M., 1996a. Differential quadrature method for thick symmetric cross-ply laminates with first-order shear flexibility. International Journal of Solids and Structures, 33 (18):2647-2658. doi: 10.1016/0020-7683(95)00174-3 Liew K. M., Han J. B., Xiao Z. M., et al., 1996b. Differential quadrature method for Mindlin plates on Winkler foundations. International Journal of Mechanical Sciences, 38 (4):405-421. doi: 10.1016/0020-7403(95)00062-3 Liu J., Wang X. W., 2008. An assessment of the differential quadrature time integration scheme for nonlinear dynamic equations. Journal of Sound and Vibration, 314 (1-2):246-253. doi: 10.1016/j.jsv.2008.01.004 Newmark N. M., 1959. A method of computation for structural dynamics. Journal of the Engineering Mechanics Division, 85 (3):67-94. http://cn.bing.com/academic/profile?id=04737dd237295d8ef1e115cb61eed609&encoded=0&v=paper_preview&mkt=zh-cn Sherbourne A. N., Pandey M. D., 1991. Differential quadrature method in the buckling analysis of beams and composite plates. Computers & Structures, 40 (4):903-913. doi: 10.1016-0045-7949(91)90320-L/ Striz A. G., Jang S. K., Bert C. W., 1988. Nonlinear bending analysis of thin circular plates by differential quadrature. Thin-Walled Structures, 6 (1):51-62. doi: 10.1016/0263-8231(88)90025-0 Wang X., Bert C. W., Striz A. G., 1993. Differential quadrature analysis of deflection, buckling, and free vibration of beams and rectangular plates. Computers & Structures, 48 (3):473-479. http://cn.bing.com/academic/profile?id=5a95130e65486d67576bb93101054491&encoded=0&v=paper_preview&mkt=zh-cn Wang X., Striz A. G., Bert C. W., 1994. Buckling and vibration analysis of skew plates by the differential quadrature method. AIAA Journal, 32 (4):886-889. doi: 10.2514/3.12071 Zeng H., Bert C. W., 2001. A differential quadrature analysis of vibration for rectangular stiffened plates. Journal of Sound and Vibration, 241 (2):247-252. doi: 10.1006/jsvi.2000.3295 -