人工边界转换方法概要
- 格式:doc
- 大小:492.00 KB
- 文档页数:17
人工智能与智能科学的边界与交叉研究人工智能(Artificial Intelligence,简称AI)和智能科学(Cognitive Science)是两个在近几十年内逐渐兴起的学科领域。
两者在研究内容和方法上有重要的联系和交叉,旨在解决人类智能及其表现的相关问题。
然而,AI和智能科学也在某些方面存在着区别和各自的研究边界。
本文将探讨人工智能与智能科学的边界以及两者之间的交叉研究。
一、人工智能的边界人工智能是一门研究如何使计算机拥有类似人类智能的技术与方法的学科。
它主要关注计算机系统的自主决策、学习能力、自然语言处理、计算机视觉和机器人等方面。
人工智能的研究可以分为强人工智能和弱人工智能。
强人工智能旨在开发具有与人类完全相似的智能,而弱人工智能则集中于特定任务的智能化和模拟。
在人工智能的边界方面,一个基本的界定是:人工智能研究的目标是创造和开发新的智能技术和应用,而不是对智能本质的详细解释或模拟。
这就意味着,人工智能关注的是如何利用计算机技术来完成复杂的任务,而不仅仅是模仿人类智能的理论与逻辑。
二、智能科学的边界智能科学是一门跨学科的科学,旨在研究和理解智能现象及其产生的基本原理。
它包括心理学、神经科学、计算机科学和哲学等多个学科的综合研究。
智能科学的核心问题是解释和理解人的智能、知觉、认知、学习以及决策等方面。
与人工智能不同,智能科学更侧重于理论框架和实证研究,以及对智能现象的深入分析。
在智能科学的边界方面,研究者们努力寻找智能背后的基本机制和原则,旨在解释人类智能活动的本质。
智能科学不仅关注智能的外部表现,也关注智能的内在作用机制。
因此,智能科学的研究范围更广泛,并且包括对生物智能、动物智能以及机器智能的研究。
三、人工智能与智能科学的交叉研究虽然人工智能和智能科学在某些方面有着各自的研究边界,但也存在着广泛的交叉研究。
这种交叉研究旨在将计算机科学和认知科学相结合,以更好地理解人类智能及其计算化实现。
静-动力分析中人工边界转换方法的研究摘要:通过将粘弹性动力人工边界应用于同时考虑静力效应和动力效应的工程算例,阐明了此类问题静-动力分析人工边界转换时保证模型为静力平衡体的必要性。
通过将粘弹性静-动力统一人工边界应用于半无限空间体有限元模型的静力分析中,验证了静力计算中的误差将使模型动力分析的稳态反应出现相近的误差。
在此基础上,系统阐述了适用于同时考虑静力效应和动力效应的工程问题的静-动力分析人工边界转换方法。
关键词:人工边界,静力分析,动力分析,边界转换Abstract:Though the application of dynamic viscous-spring artificial boundary to an engineering case with a consideration of both static and dynamic effect, and the application of the unified viscous-spring boundary for static and dynamic analysis to static analysis of a finite modal of half space, the problems of the applications of viscous-spring artificial boundary to this kind of engineering calculation was pointed out, and its corresponding solving method was proposed. On the base, a systematic switching method of these artificial boundaries was specified.Keywords: artificial boundary, static analysis, dynamic analysis, switching of boundaries1 前言人工边界从广义上可分为静力人工边界和动力人工边界。
无界区域上非齐次抛物型方程的人工边界条件法抛物型方程在科学和工程领域中具有广泛的应用,例如热传导、扩散、波动等问题。
然而,对于无界区域上的抛物型方程,由于边界条件的缺失,往往导致数值求解的困难。
为了克服这一困难,人工边界条件法被提出并得到了广泛的应用。
人工边界条件法是一种将有界区域上的问题转化为无界区域上的问题进行求解的方法。
其核心思想是通过引入人工边界,将无界区域转化为有界区域,从而可以使用传统的数值方法求解。
在非齐次抛物型方程中,人工边界条件法的基本步骤如下:首先,选择合适的人工边界,通常选择与原问题的特征线相切的直线或曲线作为人工边界。
然后,根据边界条件的类型,确定人工边界上的数值近似值。
常见的方法有零边界条件、第一类边界条件和第二类边界条件。
接下来,根据人工边界条件的近似值,将无界区域划分为有界区域和无界区域。
有界区域是指与人工边界相切的部分,通常在这个区域内使用传统的数值方法求解。
无界区域是指人工边界以外的部分,通常采用特殊的数值方法求解,例如人工边界元法、人工边界积分法等。
最后,将有界区域和无界区域的数值解进行整合,得到整个无界区域上的数值解。
通常采用插值或外推的方法进行整合,确保数值解的精确性和稳定性。
与传统的有界区域上的数值方法相比,人工边界条件法具有以下优势:一是可以处理无界区域上的问题,扩展了数值求解的范围;二是可以避免边界条件的缺失带来的困扰,提高了数值求解的精度和稳定性;三是可以减少计算量和存储空间的需求,提高了求解效率。
总之,无界区域上非齐次抛物型方程的人工边界条件法是一种有效的数值求解方法。
通过引入人工边界,将无界区域转化为有界区域,可以克服边界条件缺失的困扰,实现对无界区域上抛物型方程的精确求解。
随着数值方法的不断发展,人工边界条件法在科学和工程领域中将发挥越来越重要的作用。
人工智能赋能政府治理转型的逻辑、边界与路径作者:司林波谭筱波来源:《中共天津市委党校学报》2023年第05期[摘要]人工智能技術变革赋能政府治理转型,为政府治理从数字治理走向智慧治理提供机遇和支撑。
人工智能赋能政府治理转型,遵循情境转换、理念嬗变、技术创新、功能转化、机制完善和价值实现的逻辑机理。
面对人工智能赋能政府治理转型可能带来的挑战和风险,需要从技术应用、决策辅助、数据安全、民主参与和价值实现层面厘清赋能边界,坚持技术祛魅与技术化归相结合,强化智能决策与人类决策相融合,注重数据运作与规范管理相配合,突出政府主导与民主参与相适应,实现多元价值与统一规范相协调,通过人工智能技术的规范应用助力政府智慧治理的实现。
[关键词]数字治理;智慧治理;人工智能;治理转型;治理边界中图分类号:D63文献标识码:A文章编号:1008-410X(2023)05-0054-12一、问题的提出:人工智能时代政府治理转型人工智能技术的变革及其在社会各领域的广泛应用,重塑着个体的价值观念、行为模式、工作结构和交互关系,也影响着政府运作的内外环境,成为政府治理变革的重要变量。
党的十九届四中全会强调要推进数字政府建设,建立健全运用互联网、大数据、人工智能等技术手段进行行政管理的制度规则[1];“十四五”规划明确了加快数字化发展,建设数字中国的重点任务,鼓励支持新技术赋能[2]。
建设数字中国成为数字时代推进中国式现代化的重要引擎[3]。
2022年11月,美国Open AI公司发布的ChatGPT聊天机器人程序,引发了人们对于人工智能发展的讨论。
ChatGPT不同于以往的狭义人工智能,其以大规模预训练和自然语言处理为基石,具备上下文感知性和多领域适应性,能够通过强大的自主学习能力进行语言编辑和创作,实现与人类的深度智能交互。
人工智能技术的变革与发展引发了人们对其应用领域的讨论和探索,尤其是对政府治理变革带来的机遇和价值。
作为数字中国建设的重要课题,数字政府建设及其治理能力的提升正是对新一轮产业革命和国际竞争新焦点等战略态势的准确把握,也是对大数据驱动知识学习、跨媒体协同处理、人机协同增强智能等人工智能技术发展的积极回应[4]。
FLAC3D动力分析中的人工透射边界和地震波施加方法从动力学的角度上看,动力响应是确定惯性(质量效应)和阻尼起着重要作用时质点或质点系动力学特性和响应的技术,它包括自振、冲击、谐振动、随机振动等分支。
动力学最早应用于结构抗震设计,自上世纪50年代逐步借鉴到岩土抗震设计中。
动力发展历程可总结为静力理论,反应谱理论和时程分析理论三个阶段。
我们知道,地震的三要素为振幅、频谱和持时。
静力理论只考虑了地震引起的最大振幅,属于拟静力法;反应谱理论考虑了振幅和频谱,但在设计中仍然把地震惯性力视为静力,只能算准动力法;时程分析理论考虑了振幅、频谱和持时,是严格意义上的动力分析法。
通常时程动力分析选用的地震波来自:(1)根据设计反应谱人工合成的场地波;(2)场地附近地震台记录的实测地震波。
由于实测地震波中掺杂了许多噪声和干扰信号,因此在使用前必须滤波去噪、频谱分析、积分变换和基线修正。
滤波去噪是为了消除噪声和高频波,频谱分析是为了检测地震波持时内所含的频率分量和振幅,积分变换可以转换地震加速度波为速度波或位移波,基线修正则是为了消除非平稳地震波中的弹性位移零线漂移、基线偏移等现象,大崎顺彦在其著作《地震动的谱分析入门》中做了详细而生动的说明,并附出了地震波处理的Fortran源程序。
鉴于FLAC3D软件是岩土领域广泛应用的时程动力分析软件,这里以著名的埃尔森特罗波(El Centro)为输入激励,研究基于FLAC3D软件的地震波处理和计算方法。
网站“http://www. /data.htm”提供了31秒的El Centro加速度波数据。
有兴趣者可按《地震动的谱分析入门》的方法选取了前8秒的地震加速度波(共401个记录),然后补零配成了512个记录的加速度波以采用快速傅里叶变换法,首先采用FLAC3D Fish函数库的filter函数进行滤波去噪,然后采用fft函数进行快速傅里叶变换,得到傅里叶加速度谱和功率谱,接着采用integrate函数积分两次求得速度波和位移波,并计算地震位移零线漂移值。
无界区域问题的高阶近似人工边界条件摘要:在无界区域问题中,边界条件的确定是非常关键的。
传统的边界条件方法存在很多问题,因此人工边界条件的提出成为了一种有效的数值方法。
本文介绍了高阶近似人工边界条件的原理及其在无界区域问题中的应用。
关键词:无界区域问题,边界条件,人工边界条件,高阶近似1. 引言无界区域问题是数学、物理、工程等领域中常见的问题,如声波、电磁波、流体力学等。
在数值求解无界区域问题时,边界条件的确定是非常关键的。
传统的边界条件方法存在很多问题,如难以确定、精度不高等。
因此,人工边界条件的提出成为了一种有效的数值方法。
人工边界条件的基本思想是将无界区域的边界虚拟化为有限的区域,然后再施加适当的边界条件。
这种方法可以有效地避免传统边界条件方法中的一些问题,如边界的不确定性、精度不高等。
在实际应用中,人工边界条件已经得到了广泛的应用。
2. 传统边界条件方法存在的问题在传统的边界条件方法中,通常采用的是数学模型的边界条件。
但是,在无界区域问题中,边界往往是无穷远处,难以确定。
此外,传统边界条件方法的精度也存在一定的问题,难以满足实际应用的需求。
例如,在求解电磁波问题时,传统的边界条件方法是通过将电场或磁场在边界处等于零来确定边界条件。
但是,在实际应用中,往往难以确定边界的位置,导致边界条件的确定存在误差。
3. 人工边界条件的基本思想人工边界条件的基本思想是将无界区域的边界虚拟化为有限的区域,然后再施加适当的边界条件。
这种方法可以有效地避免传统边界条件方法中的一些问题,如边界的不确定性、精度不高等。
人工边界条件的实现需要考虑以下几个方面:(1)虚拟边界的位置和形状的确定;(2)虚拟边界处的数值计算方法;(3)虚拟边界处的边界条件的施加。
一般而言,虚拟边界的位置和形状可以根据物理问题的特点来确定。
例如,在求解电磁波问题时,可以将虚拟边界设置在离散介质的边界处。
虚拟边界处的数值计算方法可以根据具体问题来确定,一般采用数值方法进行计算。
地震波斜入射下土质边坡的稳定性分析
王雪艳;刘超;尹超;李伟华
【期刊名称】《防灾减灾工程学报》
【年(卷),期】2020(40)4
【摘要】地震诱发的边坡失稳是重要的地质灾害类型之一,数值方法是对地震边坡进行分析的一项重要工具。
基于黏弹性人工边界,通过将地震记录转化为等效地震荷载实现地震波在场地中的施加;并采用Mohr-Coulomb准则模拟土体的非线性特征。
根据非线性有限元软件ABAQUS,建立了双面边坡的数值模型,通过分析地震波在不同角度入射下不同坡角边坡的地震响应,讨论了地震波入射角以及边坡坡角对破坏的影响。
结果表明,地震波入射角和边坡形状改变了地震波在坡体内的传播规律,因此影响了边坡的破坏。
随着地震波入射角的增加,靠近波源一侧坡面的滑动位移和滑动区域逐渐增加,而远离波源一侧的滑动位移和滑动区域逐渐减小。
边坡坡度的增加引起坡面危害的逐渐加剧。
本文讨论了不同尺寸边坡在不同地震作用下的破坏规律,从而为地震诱发边坡的失稳破坏提供了理论支撑,具有重要的现实工程意义。
【总页数】7页(P589-595)
【作者】王雪艳;刘超;尹超;李伟华
【作者单位】北京交通大学土木建筑工程学院
【正文语种】中文
【中图分类】TU43;P315.9
【相关文献】
1.地震波斜入射下复杂地质中隧道动力响应分析
2.地震波斜入射下混凝土重力坝的塑性损伤响应分析
3.地震波斜入射下水下地基场地的等效线性化分析方法及应用
4.地震波斜入射下重力坝的动响应分析
5.地震波斜入射下层状TI饱和场地地震反应分析
因版权原因,仅展示原文概要,查看原文内容请购买。
水利与建筑工程学报Joprnai of Watee Resoprees and Architectural Engineering 第3卷第(期2 0 2 1年2月Vol. 19 No. 1Feb. ,202)DOI :19. 3969/C issc. 9672 - 1144.2021 61631基于黏弹性人工边界条件的岩质边坡动力反应分析景鹏旭6尹 超4,门丽君5,许杰夫6杨清逸4(1中国地震应急搜救中心,北京60042; 2.北京交通大学,北京60055;3.山西工程技术学院,山西阳泉045000 ;4.北京清控人居光电研究院有限公司,北京60089)摘要:基于数值分析方法,以二维岩质边坡为例,通过布置在边界处的弹簧和阻尼器构成黏弹性人工 边界条件,探讨边坡尺寸、形状以及地震波的类型对放大效应的影响。
结果表明地震波频率和坡角是地 震响应的决定性因素,坡高对于边坡特征点(坡脚、坡顶)的变化敏感性较低。
结论如下:对于同一地震 波频率,不同坡高对于地震放大效应的影响较小,差别仅在于边坡越高,坡脚距离坡顶的距离越长,这部分的动力响应越明显;对于同一边坡角度,随着入射地震波的频率越高,坡面的放大系数逐渐增加;对于 同一边坡高度,随着边坡的角度变陡,坡面的放大系数逐渐增加。
综上所述,在岩质边坡抗震设计时应 考虑放大效应的影响。
关键词:岩质边坡;黏弹性边界人工边界;放大效应;频率中图分类号:TU443 文献标识码:A 文章编号:372—134(2021)01—030—05Dynamic Response Analysis of Rock Slope Based onViscoelastic ArtificiaO Boundary ConditionJING Pengxu 1, YIN Chav -, MEN Lijun 5, XU Jiefu 1, YANG Qingyl 3(1. National Earthquake Response Support Service , Beijing 100049, China ;2. Beijing Jiaotong Universito , Beijing 140055 , China ;3. Shanxi Instituto r Techalogy , Yaagquaa , Shanxi 045400 , China ;4. Bejinq Tsinghua Holdings Humag Sethemenie Piqhhng Research Instituto Cn. ,Lth , Bejing 100089 , China)Abstrrci : BwS on tUe dumericai vOysis methoU , this pwer tafes a two - 47x 60X01 roch slope as an example , and discusses the inUuences of the size and shape of the slope and the type of seismic wavv on the ampliXcakon ebeO by forming a viscoelastic artificial bouudam condition with springs and dampem arraoaen at the boundarp. The msuVsshow that seismic wavv 001X0 and slope angle are the decisivv factori Ur seismic response , and slope height is less sehitav te Ovges in slope characteristic 00X01:0 ( slope U o W slope top) : The conclusions are as follows. F os the sameseismic wavv Uequency , diOerent slope heightu havv dss impact on the seismic ampliXcakon effect. The only (11X60- e n ee is that the highes the slope , the longe s the distance betwe e n the slope toe and the slope top , the more odvious thedynamic response of this pan. Fos the same slope angle , as the Uequency of the incident sb s mic wavv is highes , the slope amplification fOctos grafualip increases. Foe the same slope heighw as the slope angle becomes steepe, the slopeampliOcakon factor grafuakp increases. In summarp , the impacts of the amplification effect should be cousideren in theseismic design of roch slopes.Key w ord t : rrck slope ; viscoelastio boundary ; artiOciaO boundary ; ampliOcotion effeci : frequency收稿日期:2222-S6-77 修稿日期:22229974基金项目:中国地震应急搜救中心青年基金(SJ1991 -作者简介:景鹏旭(699—)男,山西太原人,硕士,工程师,主要从事岩土工程、地震工程等方面的研究工作。
静-动力分析中人工边界转换方法的研究摘要:通过将粘弹性动力人工边界应用于同时考虑静力效应和动力效应的工程算例,阐明了此类问题静-动力分析人工边界转换时保证模型为静力平衡体的必要性。
通过将粘弹性静-动力统一人工边界应用于半无限空间体有限元模型的静力分析中,验证了静力计算中的误差将使模型动力分析的稳态反应出现相近的误差。
在此基础上,系统阐述了适用于同时考虑静力效应和动力效应的工程问题的静-动力分析人工边界转换方法。
关键词:人工边界,静力分析,动力分析,边界转换Abstract:Though the application of dynamic viscous-spring artificial boundary to an engineering case with a consideration of both static and dynamic effect, and the application of the unified viscous-spring boundary for static and dynamic analysis to static analysis of a finite modal of half space, the problems of the applications of viscous-spring artificial boundary to this kind of engineering calculation was pointed out, and its corresponding solving method was proposed. On the base, a systematic switching method of these artificial boundaries was specified.Keywords: artificial boundary, static analysis, dynamic analysis, switching of boundaries1 前言人工边界从广义上可分为静力人工边界和动力人工边界。
静力人工边界由来已久,通常有固定边界、滚轴边界等。
动力人工边界经过几十年的研究发展,已形成具有全局人工边界和局部人工边界的两大类别,并应用于各自适应的工程计算中[1]。
动力人工边界发展到现在已有透射边界、粘性边界、粘弹性边界等几种类型。
1994年,Deeks 提出粘弹性人工边界[11]。
1998年,刘晶波等人发展了二维的黏弹性人工边界[3],又于2005年将其发展为三维时域黏弹性人工边界[4]。
2006年,刘晶波等人再将二维黏弹性边界发展成一致粘弹性人工边界及其对应的粘弹性边界单元[5],并于2007年推导了三维一致粘弹性人工边界及等效粘弹性边界单元[6]。
目前对静-动力分析的普遍做法是采用静力人工边界和动力人工边界分别对静力问题和动力问题进行计算,将计算结果进行叠加后得到完整的结果[1]。
但由于叠加原理仅在线弹性小变形范围内适用,原则上不能应用于涉及非线性或大变形问题的分析。
目前对涉及非线性或大变形问题的静-动力分析,常用的人工边界转换方法主要有以下几种:(1)静力分析和动力分析都采用滚轴边界或固定边界;(2)静力分析中采用滚轴边界或固定边界,动力分析采用粘弹性边界、透射边界、粘性边界等人工边界;(3)静力分析和动力分析都采用静-动力统一边界,如粘弹性静-动力统一人工边界。
对第(1)种方法,由于固定边界使波动全部反射,已有许多文献证明其具有放大振动效应的作用,目前已经使用得不多。
刘晶波等人基于黏弹性动力人工边界和半无限空间中静力问题的基本解,建立了对动力问题和静力问题均适用的三维黏弹性静-动力统一人工边界,从而上述第(3)种方法得以解决[1]。
然而,在使用人工边界对地下结构进行动力分析时,还存在一些问题。
如第(2)种方法,由于在静-动力分析的人工边界转换时的方法存在问题,致使产生错误的结果。
在第(3)种方法中,将粘弹性静-动力统一人工边界应用于地下结构的静力分析时,其解与准确值存在误差。
本文将就此两问题进行论证和分析,并阐述合理的地下结构静-动力分析人工边界转换方法。
2 静力和动力有限元分析原理2.1 静力分析原理地铁等地下工程初始应力场的确定须先计算未开挖状态下围岩的自重应力场,进而根据施工步骤,采用释放荷载法,计算出衬砌结构和围岩的静应力场[13]。
许多地下结构的自重应力场模型可以假设为半无限空间体,根据经典围岩压力理论和弹性力学理论,半无限空间体中距地表面任一深度h 处的应力状态可定义为[7][8]:h V γσ= (1) h V H λγλσσ== (2)式中,σV 为竖向应力,γ为围岩重度,σH 为横向应力,λ为侧压力系数。
对浅层围岩,可假设其为各向同性介质,侧压力系数可用泊松比表示如下[7]:μμλ-=1 (3) 对有限元计算中所取的有限区域,可以据此确定有限域边界条件。
2.2静力分析边界条件根据上述静力分析方法,在静-动力共同作用问题的计算中,在进行动力分析之前须先确定地下结构模型的静应力场。
地下结构开挖前可将大地假设为半无限空间体,其在重力作用下的静力计算,根据对称性,模型中任一处的水平位移0=H u 。
因此,计算该应力场时,有限区域模型两侧可用水平约束即法向约束,底部可用全约束或仅约束竖直方向,顶面即地面应为自由边界。
2.3 动力分析原理在动力荷载作用下,有限元体系在t +∆t 时刻的运动平衡方程为:t t t t t t t t F u K u C uM ∆+∆+∆+∆+=++ (4) 式中M 为体系的总质量矩阵;C 为体系的总阻尼矩阵;K 为体系的总刚度矩阵; ut t +∆为体系的节点加速度向量; ut t +∆为体系的节点速度向量;u t t +∆为体系的节点位移向量;t t F ∆+为外荷载向量。
体系的总阻尼矩阵采用瑞利阻尼[9]:C M K =+αβ (5)式中α、β为常数,可按两种不同的振动频率下测得的阻尼比ξ加以确定。
计算中常数α、β可由α+βω2i =2ωi ξi 和α+βω2j =2ωj ξj 求得[9]。
则α和β可表示为:j i j i j i j i i j ωωωωωωωξωξα))(()(2-+-=(6)))(()(2j i j i j j i i ωωωωωξωξβ-+-=(7) ωi 、ξi 分别为振型向量φi 对应的自振圆频率和阻尼比。
根据振型分析结果可求得ωi 和ωj ,阻尼比在计算中取ξi =ξj =0.05。
2.4 粘弹性人工边界粘弹性人工边界从用途上可分为动力人工边界和静-动力统一人工边界;从具体实现方法上可分为弹簧-阻尼器边界单元和一致粘弹性边界单元。
下面先介绍动力人工边界。
具体计算中模型边界材料参数由其相邻的围岩介质材料决定。
则当人工边界采用等效的弹簧和阻尼器物理元件来模拟时,其弹簧系数和阻尼系数的计算分别如下[3][4]:法向边界:)1(2μαα+==R ER G K TTT , s T c C ρ= (8) 切向边界:)1(2μαα+==R ER G K NNN , p N c C ρ= (9) 式中K T 、K N 分别为法向与切向的弹簧刚度;R 为波源至人工边界点的距离;c s 和c p 分别为S波和P 波波速;E 和G 分别为介质弹性模量和剪切模量;μ为介质泊松比;ρ为介质质量密度;αT 与αN 分别为切向与法向粘弹性人工边界参数,具体取值情况见表1。
表1 粘弹性动力人工边界中参数α的取值模型类型 方向 α 二维人工边界平面内法向 2.0 平面内切向 1.5 出平面切向 0.5 三维人工边界法向 4.0 切向2.0若人工边界采用一致粘弹性边界单元来模拟,其边界单元等效剪切模量、等效弹性模量和阻尼系数可分别用如下几式计算[5][6]:)1(2~μαα+==R EhR G hG T T (10) ~~~~~~~)1()21)(1()1(2)1()21)(1(νννμανννα--++=--+=R E h R G hE N N (11)])1[()1(2])1[(~Np T s NpTsc c n nE R c c n nGRααμρααρη+-+=+-=(12)其中等效阻尼系数取的是法向和切向的平均值。
式中h 为等效边界单元厚度,~ν为等效泊松比,n 为计算模型维数。
其余符号意义同上。
表2 粘弹性静-动力统一人工边界中参数α的取值模型类型 人工边界位置 方向 α 二维人工边界模型底面法向 α*/3 切向 2/3 模型侧面 法向 α*/4 切向 1/8 三维人工边界模型底面法向 α* 切向 2 模型侧面法向 α* 切向1/2考虑到实施的方便,实际计算中边界材料常采用各向同性材料,这时上述等效剪切模量和等效弹性模量之间存在隐含关系式~E =2(1+~ν)~G 。
考虑到普通有限元材料泊松比应限制在0~0.5范围内,则该等效泊松比可按如下取值[5]。
⎪⎩⎪⎨⎧≥--=其它,,T N T N T N 02/)1/(22/~ααααααν (13) 即计算等效泊松比时应先确定T N αα/大小范围,再确定其计算式。
静-动力统一人工边界是在上述动力人工边界的基础上对人工边界参数α进行调整,其具体取值见表2[1]。
其中参数α*具体计算公式:])(1[])(1)[1(26222*RdR d +++-=μα (10)式中μ为围岩泊松比,d 为位置坐标,R 为荷载作用点到人工边界点的距离。
对于底面人工边界,式中D 取荷载作用点至边界单元的水平距离;对于侧面人工边界,式中D 取荷载作用点至边界单元的垂直距离。
3 静-动力分析中的人工边界转换的几种方法在引言中已经介绍了地下结构静-动力分析中常用的三种人工边界转换方法,下面采用算例对第二种和第三种方法存在的问题进行分析。
3.1算例一某双线铁路隧道,其净空跨度10.2m ,为三心圆断面;衬砌初期支护15cm ,二次衬砌30cm ;隧道埋深45m ,Ⅳ类围岩。
采用ANSYS 大型有限元分析软件对该隧道进行建模和计算,模型区域竖向取100m ,横向取100m 。
隧道及围岩物理力学参数见表4。
静力分析采用的边界条件如图1所示,其中g 为重力加速度。
动力分析采用如图2如示动力模型,其最外一层即为粘弹性动力人工边界,厚度2米。
图中P (t )为某实测的列车振动荷载,其时程曲线见图3。
第(2)方法在静力分析中采用滚轴边界或固定边界(可称之为传统静力边界条件),在动力分析中采用粘弹性动力人工边界,在静力分析和动力分析中的荷载和边界条件见表3。
图1隧道静力分析边界条件 图2 动力分析边界条件表3 计算步骤及条件表表4 围岩物理力学参数表图3 列车荷载时程曲线由于隧道结构是对称的,整个模型有限域以及荷载也是对称的,监测点可以取自其对称轴及其一侧即可;取隧道衬砌的拱顶、拱腰、拱脚以及仰拱底等几个点监测其动力响应情况。