Estimating the spatial distribution of coseismic slip is an ill-posed inverse problem, and solutions may be extremely oscillatory due to measurement errors without any constraints on the coseismic slip distribution. I...Estimating the spatial distribution of coseismic slip is an ill-posed inverse problem, and solutions may be extremely oscillatory due to measurement errors without any constraints on the coseismic slip distribution. In order to obtain stable solution for coseismic slip inversion, regularization method with smoothness-constrained was imposed. Trade-off parameter in regularized inversion, which balances the minimization of the data misfit and model roughness, should be a critical procedure to achieve both resolution and stability. Then, the active constraint balancing approach is adopted, in which the trade-off parameter is regarded as a spatial variable at each model parameter and automatically determined via the model resolution matrix and the spread function. Numerical experiments for a synthetical model indicate that regularized inversion using active constraint balancing approach can provides stable inversion results and have low sensitivity to the knowledge of the exact character of the Gaussian noise. Regularized inversion combined with active constraint balancing approach is conducted on the 2005 Nias earthquake. The released moment based on the estimated coseismic slip distribution is 9.91×1021 N·m, which is equivalent to a moment magnitude of 8.6 and almost identical to the value determined by USGS. The inversion results for synthetic coseismic uniform-slip model and the 2005 earthquake show that smoothness-constrained regularized inversion method combined with active constraint balancing approach is effective, and can be reasonable to reconstruct coseismic slip distribution on fault.展开更多
柔性绳索结构因质量轻、高柔性和可设计性强等特点,在航空航天及机器人等工程领域得到广泛应用.由于几何大变形与材料非线性等因素影响,利用非线性有限元等高精度方法分析这类柔性结构的动力学特性时,存在模型复杂、计算量大和难以高效...柔性绳索结构因质量轻、高柔性和可设计性强等特点,在航空航天及机器人等工程领域得到广泛应用.由于几何大变形与材料非线性等因素影响,利用非线性有限元等高精度方法分析这类柔性结构的动力学特性时,存在模型复杂、计算量大和难以高效仿真等问题.文章改进了起源于计算机图形学的XPBD(extended position based dynamics)算法,设计了新的约束函数与迭代方法,使其适用于柔性绳索结构的动力学仿真.主要贡献包括:引入旋转向量以更准确描述结构姿态,基于Cosserat弹性杆理论设计了更精确的约束函数,采用约束能量作为迭代收敛判据以提高仿真精度,提出分组求解约束的方法提高迭代计算效率.还以索-杆复合结构为例验证了改进算法的有效性:与ADAMS软件相比,改进算法具有更好的稳定性,且计算效率大致相当;与原XPBD算法相比,改进算法显著提高了求解精度.改进算法引入了旋转向量,使约束函数的设计更灵活,计算框架也更利于并行化,具有在复杂柔性结构动力学仿真中进一步应用的优势和潜力.展开更多
基金Projects(41604111,41541036) supported by the National Natural Science Foundation of China
文摘Estimating the spatial distribution of coseismic slip is an ill-posed inverse problem, and solutions may be extremely oscillatory due to measurement errors without any constraints on the coseismic slip distribution. In order to obtain stable solution for coseismic slip inversion, regularization method with smoothness-constrained was imposed. Trade-off parameter in regularized inversion, which balances the minimization of the data misfit and model roughness, should be a critical procedure to achieve both resolution and stability. Then, the active constraint balancing approach is adopted, in which the trade-off parameter is regarded as a spatial variable at each model parameter and automatically determined via the model resolution matrix and the spread function. Numerical experiments for a synthetical model indicate that regularized inversion using active constraint balancing approach can provides stable inversion results and have low sensitivity to the knowledge of the exact character of the Gaussian noise. Regularized inversion combined with active constraint balancing approach is conducted on the 2005 Nias earthquake. The released moment based on the estimated coseismic slip distribution is 9.91×1021 N·m, which is equivalent to a moment magnitude of 8.6 and almost identical to the value determined by USGS. The inversion results for synthetic coseismic uniform-slip model and the 2005 earthquake show that smoothness-constrained regularized inversion method combined with active constraint balancing approach is effective, and can be reasonable to reconstruct coseismic slip distribution on fault.
文摘柔性绳索结构因质量轻、高柔性和可设计性强等特点,在航空航天及机器人等工程领域得到广泛应用.由于几何大变形与材料非线性等因素影响,利用非线性有限元等高精度方法分析这类柔性结构的动力学特性时,存在模型复杂、计算量大和难以高效仿真等问题.文章改进了起源于计算机图形学的XPBD(extended position based dynamics)算法,设计了新的约束函数与迭代方法,使其适用于柔性绳索结构的动力学仿真.主要贡献包括:引入旋转向量以更准确描述结构姿态,基于Cosserat弹性杆理论设计了更精确的约束函数,采用约束能量作为迭代收敛判据以提高仿真精度,提出分组求解约束的方法提高迭代计算效率.还以索-杆复合结构为例验证了改进算法的有效性:与ADAMS软件相比,改进算法具有更好的稳定性,且计算效率大致相当;与原XPBD算法相比,改进算法显著提高了求解精度.改进算法引入了旋转向量,使约束函数的设计更灵活,计算框架也更利于并行化,具有在复杂柔性结构动力学仿真中进一步应用的优势和潜力.
