简介:微分求积法(DQM)是1种求解微分方程初(边)值问题的数值方法,通常以较小的计算工作量即可获得较高的数值精度。这种方法应用于工程领域时多用来解决梁、板等结构的静力分析或结构特征值分析等问题,即对边值问题的微分方程的求解。结构动力分析属于初值问题,荷载和结构反应都具有特殊性,直接套用DQM求解边值问题并不能获得问题的解。本文尝试利用微分求积原理建立求解结构动力反应的具体方法。借鉴单元法的思想,将荷载持时划分为若干个时步,在每个时步内对动态荷载和结构反应进行离散,然后用DQM对时步逐个进行求解,得到体系在整个时域内的反应过程。通过对3种不同自振周期的线弹性单自由度体系在不同频率简谐激励下反应的计算,阐释了本文方法的可行性以及高精度、高效率的特点,通过数值试验确定了时步内相对较优的节点数,并为时步长度的选取提供了建议。
简介:
简介:有效应力〈σij〉的精确表达式,特别是导致有孔隙流体材料的弹性应变的应力〈P〉(有效应力)是基于假设之上的,仅对Hook定律有效,即〈σij〉=σij-αPδij和〈P〉=Pc-Pp,这里α=1-(K/Ks),Pc和Pp是围压和孔隙压力,K和Ks分别是干燥岩石(排水状态下)和岩石基质(岩石固体部分)的体积模量。Gccrtsma(1957年)和Skempton(1960年)在实验的基础上首次提出关于〈P〉的方程。该表达式虽然不直接依赖于孔隙度,但是当有效应力〈P〉等于围压时孔隙消失,因此K=Ks。如果使用〈σij〉方程中的有效应力定理,那么可根据没有孔隙压力的固体弹性模量确定一个具有孔隙压力的多孔固体的应变。有效应力表达式非常准确地描述了砂岩和花岗岩样品在围压和孔隙压力达到2.5kb(250MPa)时的应变。结果表明,该有效应力定理不适用于非弹性过程(如断裂)。