FLAC,FLAC3D基础与应用-陈育民PPT课件

1,FLAC/FLAC3D基础与应用,陈育民河海大学土木与交通学院2015年6月2日,河海大学研究生课程岩土数值分析,关于教师,2002年本科毕业于中国矿业大学2007年博士毕业于河海大学2013.7-2014.7美国普渡大学访问学者2015.22015.5日本东京大学访问学者研究方向：土动力学与岩土地震工程土木工程防灾减灾岩土工程数值分析联系方式：ymchenhhu,3,关于教材,4,关于课程,2005-11-29河海土木院研究生会组织2006-10-13同济大学土木工程学院2006-10-26河海大学金水节2007-04-15东南大学交通学院2007-07-18同济大学土木工程学院2007-11-03河海大学岩土所组织FLAC学术沙龙2007-11-29河南工业大学2008-11-15河海大学河海金水节培训2010-11-10河海大学校庆报告2011-06-18河海大学举办ITASCA技术与应用专题（南京）研讨会2011-10-16河南理工大学2011-11-03南京工业大学交通学院2011-11-24河海大学土木与交通学院研究生会2011-06-18_ITASCA技术与应用专题（南京）研讨会2012-08-31_解放军理工大学FLAC讲座,课程目的,什么是FLAC？为什么要用FLAC？FLAC能做什么？FLAC为何这么流行？怎么学FLAC？,6,课堂“作业”,研究生课程岩土数值分析上课学生调查,7,课堂“作业”,8,什么是FLAC？,FastLagrangianAnalysisofContinua,9,为什么要用FLAC？,检索期刊：岩土工程学报岩土力学岩石力学与工程学报关键词：,UpdatedonJune2,2015,10,为什么要用FLAC？,检索期刊：JournalofGeotechnicalandGeoenvironmentalEngineeringGeotechniqueCanadiangeotechnicalJournalSoilsandfoundations关键词：,UpdatedonJune2,2015,11,为什么要用FLAC？,检索期刊：GEO*SOILS*关键词：,UpdatedonJune2,2015,12,FLAC能做什么？,岩土工程中的绝大多数问题土力学、岩石力学、防灾减灾、隧道、地下空间等采矿工程中的大部分问题水工结构中的部分问题结构工程,国际通用的岩土工程专业分析程序,13,FLAC为何这么流行？,CharlesFairhurst美国工程院、瑞典皇家工程院院士，国际岩石力学学科和岩石力学学会创始人之一，历任国际岩石力学学会主席和副主席，国际岩石力学学会MULLER奖、美国岩石力学学会终生成就奖获得者。PeterCundall美国工程院、英国皇家工程院院士，国际资深计算岩石力学学家。,14,课程安排,第一讲：基本介绍、静力分析、前后处理第二讲：接触面、FISH语言、流固耦合分析第三讲：动力分析、自定义本构、结构单元第四讲：FLAC（2D）基本介绍与应用实例讨论,课程要求,了解FLAC的基本概念了解FLA的适用范围和缺陷了解深入学习FLAC的方法会用FLAC分析简单的路堤填筑问题（考题）,题外话“虚拟”与“现实”,现实的复杂与虚拟的简化,是否一定要“像”,我们也做过很像的东西,认清土体数值分析的影响因素,复杂的网格往往会把核心问题覆盖掉，经常会造成计算结果无法解释。,虚拟是简单的世界,现实中的桩,虚拟的桩,虚拟是简单的世界,冰碛土体结构模拟结果,台阶坡面上的砾石产出状态,小结,数值分析要对实际工程进行大量的、细心的简化，从效率和结果两方面保证数模分析的成果结果的判断，需要深入扎实的理论功底及“丰富”的工程经验数值分析的作用从“锦上添花”到“雪中送炭”，在于各位的努力，指日可待,24,第一讲,FLAC3D基本介绍、静力分析、前后处理,25,软件介绍,FastLagrangianAnalysisofContinua美国Itasca咨询公司开发2D程序(1986)1990年代初引入中国有限差分法(FDM)3D版本：DOS版2.02.13.03.14.05.02D版本：DOS版4.05.06.07.0,26,FLAC,大应变、小应变计算模式.丰富的本构模型、提供自定义的本构模型功能接触面可以模拟不同材料的接触流固耦合实现土体的固结与渗流拥有各种功能的结构单元类型，模拟土与结构的相互作用强大的动力分析功能.流变分析，拥有粘弹性模型和粘塑性模型热力学分析.,Shearstrainratecontours,27,FLAC3D,与FLAC类似，是FLAC的三维版本与FLAC拥有相同的优点,upstream,downstream,28,基本特点,内置材料模型连续介质非线性，大应变模拟显式解题方案，为不稳定物理过程提供稳定解界面或滑动面用来模拟可产生滑动或分离的离散面，从而模拟断层，节理或摩擦边界内置材料模型丰富：零模型,三个弹性模型(各向同性，横观各向同性和正交各向异性),八个朔性模型(德鲁克-布拉格,摩尔-库伦,应变硬化/软化，单一节理，双线性应变硬化/软化单一节理,双屈服，修正剑桥粘土，霍克-布朗),29,可选模块,可选模块包括:热力学,热-力学耦合,热-流体-力学耦合包括热传导和对流;粘弹,粘朔性(蠕变)材料模型;动力学分析,并可以模拟静边界和自由域使用C+定义自己的模型,核废料储存中的热力学研究问题,30,FLAC/FLAC3D基本原理,FLAC/FLAC3D利用有限差分，显示方案，动态松弛方法模拟连续体的非线性力学行为：即使对准静态问题，程序仍然求解完整的动力学方程。这种方法的好处在于可以为物理非稳定过程例如塌方提供稳定解；在“松弛”方法中，使用阻尼来吸收动能以模拟系统的“静态”反应。这种方法可以用比其它方案如解矩阵法更为真实有效地模拟塌方问题。,31,Lagrangian法,源自流体力学中的拉格朗日法跟踪流体质点的运动状态跟踪固体力学中结点，按时步用Lagrangian法研究网格节点的运动节点和单元随材料移动，边界和接触面与单元的边缘一致固体力学大变形理论,法国数学家、物理学家拉格朗日,32,混合离散技术,+,/2,=,每个,为常应力/应变:,体积应变由整个四边形算出,.应变偏量则有两个三角形,和,分别算出,(混合离散过程),解题过程中网格坐标按照“拉格朗日方式更新”(网格随材料移动),且为显式(一个时步内局部变化不会影响邻域),33,混合离散技术,FLAC3D混和离散,+,/2,=,34,FLAC3D混和离散,结构域离散为可由四面体单元组合形成的五面体或六面体等单元;以为基本单元(常应力、常应变);体应变的计算：;偏应变的计算：.,35,动态松弛,动态松弛法在动态松弛法中，网格点根据牛顿运动定律运动.网格点的速度与该点的不平衡力呈正比.这种求解方法所决定的一系列位移将把系统带入平衡状态，或表明破坏模式.在动态松弛法中有两个因素很重要:时步的选择阻尼效应,36,显式算法,显式解与隐式解的比较,显式,逐时推进,隐式,静态,1.无需进行反复迭代来实现非线性本构关系.2.类似问题求解时间呈N3/2规律增长3.物理非稳定性不会引起数值不稳定性.4.因为无需储存矩阵，用较小内存即可模拟大尺度问题.5.对大位移、大应变问题同样适合，无需额外的计算.,1.需进行反复迭代来实现非线性本构关系2.类似问题求解时间呈N2甚至N3规律增长.3.难以模拟物理非稳定性问题.4.需存储刚度矩阵，需克服相关的带宽问题，需要的内存较大.5.对大位移、大应变问题需进行大量的计算.,37,NewFeaturesinFLACVersion6.0,使用IntelFortrancompiler拥有更快的计算速度自动网格重画功能，解决bad-geometry问题.新的模拟颗粒土材料的硬化模型更新的通用网格生成工具,38,NewFeaturesinFLAC3DVersion3.1,多处理器的并行计算功能新结构单元类型“EmbeddedLiner”提供两个方向的接触作用，可以很好地模拟挡土墙对四面体单元采用新的混合离散方法“NodalMixedDiscretization”提供塑性问题更精确的解答64位程序包含命令手册、FISH手册和应用实例的帮助,39,NewFeaturesinFLAC3DVersion4.0,模拟颗粒状材料的硬化模型自动网格重画功能，解决bad-geometry问题.改进的interface更快的渗流计算更新的动力计算功能,40,Lagrangian格式动量平衡方程,F(t),m,牛顿运动定律,对于连续体,在静力平衡条件下，加速度项为0，方程变为平衡方程,41,自由落体的模拟,G=mg,S=1/2gt2=20m,命令流：configdyngenzonbrisize111inixmul0.1ym0.1zm0.1modelelaspropbulk3e8shear1e8inidens1000setgrav00-10solveage2,42,自由落体的模拟(movie),43,FLAC3D中模型术语,节点,gridpoint：节点zone：单元boundary：边界,44,FLAC3D的求解过程,45,FLAC3D中的本构模型,开挖模型null3个弹性模型各向同性弹性横观各向同性弹性正交各向同性弹性8个弹塑性模型Drucker-Prager模型、Morh-Coulomb模型、应变硬化/软化模型、遍布节理模型、双线性应变硬化/软化遍布节理模型、修正剑桥模型和胡克布朗模型,46,FLAC3D中的本构模型,47,一个最简单的例子,genzonbrisize333;建立网格(前处理)modelelas;材料参数propbulk3e6shear1e6inidens2000;初始条件fixzranz-.1.1;边界条件fixxranx-.1.1fixxranx2.93.1fixyrany-.1.1fixyrany2.93.1setgrav00-10solve;求解appnstr-10e4ranz3x12y12solveploconzd;后处理切片功能,RUNFLAC3D,48,分析问题的过程,建立网格,初始条件,边界条件,初始应力平衡,外荷载,求解,前处理,后处理,49,FLAC3D的文件格式,保存文件(*.sav)含有所有状态变量和用户定义条件的二进制文件数据文件(*.dat)数据文件由用户创建的一种ASC格式的文件，它包括一系列的用于描述所分析问题的FLAC3D命令FISH文件(*.fis)FISH程序文件FLAC3D文件(*.flac3d)FLAC3D的网格信息文件历史记录文件(*.his)记录输入输出历史值的文件图形文件图形文件(各种标准格式)电影文件(*.dcx)AVI或PCX图像文件，这些图像文件可以当作电影放映,50,初始应力的生成,为什么要单独列出？分析过程中出现的很多问题都与初始应力是否合理有关手册中的例子五花八门是所有后续分析的基础！生成方法弹性求解更改强度参数的弹塑性求解设置初始应力的弹塑性求解存在水压力的初始应力生成水下建筑的初始应力生成,51,弹性求解,genzonbrisize112melaspropbulk3e7shear1e7fixzranz0fixxranx0fixxranx1fixyrany0fixyrany1inidens2000setgrav00-10solve,Step=162sz=-40e3sx=-21.54e3,52,更改强度参数的弹塑性求解,genzonbrisize112modelmohrpropbulk3e7shear1e7c1e10f15ten1e10fixzranz0fixxranx0fixxranx1fixyrany0fixyrany1inidens2000setgrav00-10solvepropbulk3e7shear1e7c10e3f15ten0solve,Step=163sz=-40e3sx=-21.54e3,Or:solveelastic,53,设置初始应力的弹塑性求解,genzonbrisize112modelmohrpropbulk3e7shear1e7c10e3f15ten0fixzranz0fixxranx0fixxranx1fixyrany0fixyrany1inidens2000iniszz-40e3grad0020e3ranz02inisyy-20e3grad0010e3ranz02inisxx-20e3grad0010e3ranz02setgrav00-10solve,Step=0sz=-40e3sx=-20e3,54,存在水压力的初始应力生成(1),genzonebricksize112modelmohrpropbulk3e7shear1e7coh10e3fri15ten0fixzranz0fixxranx0fixxranx1fixyrany0fixyrany1inidens2000ranz01inidens1500ranz12iniszz-35e3grad0020e3ranz01inisyy-22.5e3grad0015e3ranz01inisxx-22.5e3grad0015e3ranz01iniszz-30e3grad0015e3ranz12inisyy-15e3grad007.5e3ranz12inisxx-15e3grad007.5e3ranz12inipp10e3grad00-10e3ranz01setgrav00-10solve,Step=0sz=-35e3sx=-22.5e3,rd=rsnsrf,55,存在水压力的初始应力生成(2),configfluidgenzonbrisize112modelelaspropbu3e7sh1e7inidens1500modelfl_isoinifdens=1000fmod0proppor0.5setgrav00-10waterdens1500watertableface001,011,111,101inipp10e3grad00-10e3ranz01iniszz-30e3grad0015e3ranz12.setfluidoffsolve,Step=142sz=-40e3sx=-24e3,单元数较少产生的误差,56,水下建筑的初始应力生成,genzonbrisize112modelmpropbulk3e7shear1e7c10e10f15ten1e10fixzranz0fixxranx0fixxranx1fixyrany0fixyrany1inidens2000ranz02iniszz-50e3grad0020e3ranz02inisyy-40e3grad0015e3ranz02inisxx-40e3grad0015e3ranz02inipp30e3grad00-10e3ranz02appnstress-10e3ranz2setgrav00-10solve,Step=0sz=-50e3sx=-40e3,57,前后处理,58,前后处理,基本前后处理命令操作菜单操作dd100GPa(岩石)反分析方法：通过断层中岩石的变形与原岩的变形,单桩承载力分析,软土地基bulk1.6878E6shear3.6167E5coh15E3fric12dens1.73E3桩体bulk5e9shear3.75e9dens2.5e3,0.5m,8m,10m,20m,计算过程,施加桩顶荷载,计算结果,影响因素,7.50E+03,fric,g,k,3,coh/0.7,0.7fric,10g,10k,2,3.00E+04,20,100e,100e,1,coh,fric,ks,kn,水平因素,计算工况设计,1,2,3,3,(9),3,1,2,3,(8),2,3,1,3,(7),2,1,3,2,(6),1,3,2,2,(5),3,2,1,2,(4),3,3,3,1,(3),2,2,2,1,(2),1,1,1,1,(1),coh,fric,ks,kn,水平因素工况,计算结果,Ks取1,Ks取2,Ks取3,96.9%,1,2,3,3,(9),22.4%,3,1,2,3,(8),54.1%,2,3,1,3,(7),96.9%,2,1,3,2,(6),7.1%,1,3,2,2,(5),37.8%,3,2,1,2,(4),96.9%,3,3,3,1,(3),42.9%,2,2,2,1,(2),37.8%,1,1,1,1,(1),Ra差异度,coh,fric,ks,kn,水平因素,最优方案,2,2,2,2,11,1,1,2,2,10,1,1,2,2,最优方案,17.3%,6.8%,72.8%,11.9%,极差,52.4%,52.7%,96.9%,57.8%,k3,64.6%,59.2%,24.1%,47.3%,k2,47.3%,52.4%,43.2%,59.2%,k1,coh,fric,ks,kn,水平因素,合理步骤,单桩分析,简单网格,接触面参数,多次试算,理论、实测,加密网格,接触参数,理想结果,群桩分析,Pile结构单元,单元参数,理想结果,Pile结构单元,单元参数,理想结果,挡土墙的接触面设置,对于未知问题的分析思路3个独立的接触面3个同ID的接触面2个独立的接触面2个独立的接触面并进行底部merge,挡墙,土体,1,2,3,Wall,Soil,挡土墙的接触面设置(2),Z=5.16cm,Z=1.14cm,Z=不能初始平衡,Z=0.35cm,？,挡土墙的接触面设置(3),X=1.95cm,X=1.91cm,X=1.47cm,？,挡土墙的接触面设置(3),前两种方法的差别实质3个ID的独立接触面在相同位置产生互不影响的两个节点共同ID的接触面在相同位置自动设置为1个节点最终的结论需要您自己去判断！,3interfaces,2IDs,1interface,1ID,FISH语言,FISH语言简介,软件自带的编程语言是否一定要学？视情况而定，需要时查询FISH变量即可语法简单xxxend_xxx注意事项与FLAC本身的关键字冲突保留字不可缩写变量可不定义，因此注意检查程序printfishtable,extra等命令使用,一个最简单的FISH程序,defabcabc=1+2*3abcd=1.0/2.0endabcprintfish,数据格式,函数与变量,都可以在FISH函数中进行赋值，赋值操作与常规的编程语言类似，按照运算符的优先级先后顺序来执行。函数和变量的赋值遵守数据类型的规则，即整型的计算结果为整型，浮点型的计算结果为浮点型，因此读者在进行除法运算、开方运算时都需要将数据类型设置为浮点型，数字尽量使用小数点以保证运算正确。变量和函数名的命名规则是不能以数字开头，不能含有中文，并且不能包含如下的字符。.,*/+-=#();变量和函数名不能与FLAC3D、FISH的保留字相冲突，不要采用过于简单的单词，比如a，hist等，这些都与保留字相冲突。即使程序中存在与保留字相冲突的变量，FLAC3D也不会提供任何提示，所以提醒读者在编制FISH程序时尽量使用较长的、复杂的变量和函数名。,函数与变量,对变量进行赋值时，不能使用当前函数的函数名放在“=”的右边，比如采用下面的定义abcd=abc+1.0在FISH程序执行时会提出错误，因为这样会形成递归调用，这种调用方式在FISH程序中是不允许的。变量和函数的作用是全局的，在命令中的任何地方修改变量的值都会立即生效，因此在实际应用中尽量避免不同的函数中含有相同的变量，因为这样可能会造成赋值错误，并难以检查。在FLAC3D中可以用如下的命令来引用FISH函数和变量PRINT用于查看函数和变量的数值；HISTORY命令可以对函数和变量的数值进行记录；SET命令用于变量的赋值。,主要语句,选择语句CASEOF表达式默认语句CASEn1表达式的值为n1时的语句CASEn2表达式的值为n2时的语句ENDCASE,主要语句,条件语句IF条件表达式THENELSEENDIF,FISH中条件运算符没有“并”、“或”、“否”这样的符号,ifaa1.0ifaa2.0执行语句endifendif,表达“1K+4/3G)湿法：耦合体系的短期行为使用排水的K,c,若SETfluidoff,Biot_mod(fmod)真实,力学过程的特征时间流体扩散过程的特征时间,完全耦合分析方法,时间比例短期行为(不排水)ts(分析时间)tc施加扰动的属性流体扰动：渗流可不与力学过程耦合力学扰动：耦合等级取决于流固刚度比流固刚度比,流-固耦合计算,CONFIGfluid;M(Kf);K(渗透系数)真实，则FLAC3D默认耦合计算pevsevp预估流/力特征时间耦合计算前先达到一个平衡状态SETfluidonmechoff;SETfluidoffmechon;STEPSETmechforce;SETmechsubstepnauto;SETfluidsubstepm(=1)STEP:渗流步足够小,流固耦合的计算方法,手动调整的STEP求解SETfluidonmechoffSTEPSETfluidoffmechonSTEP主从进程的SOLVE求解SETmechforceSETmechsubstepnauto(从进程)SETmechsubstepm(主进程)SOLVEage自动STEP求解STEP,渗流问题(CONFIGfluid)分析步骤,时间比例(ts,tc)稳态不排水状态相当扰动类型力学扰动孔压扰动流固刚度比Rk是否1完全耦合模式时间比例相当；力学扰动,心墙土坝的渗流(1),newconfigfluidsetfluidoffgenzonbrickp000-10size20110genzonbrickp0500p11500p2510p3905p41510p5915p61105p71115size1015groupsoilgroupdamranx57z-50groupdamranid201aid211aid221aid231aid241agroupdamranid202aid212aid222aid232aid242amepropbu3e7sh1e7inipp0grad00-10e3ranz0-10inidens2000modelfl_isoproppor0.5perm1e-10inifden1000ften-1e10inisat0.0ranz05modelfl_nullrangrodam;inipp0rangrodamfixzranz-10fixxranx0fixxranx20fixysetgrav10solvesaveelastic.sav,网格模型,初始孔压,心墙土坝的渗流(2),restelastic.savinixd0yd0zd0 xv0yv0zv0appnstress-40e3grad0010e3ranz04x09solvesavepressure.sav,竖向应力,沉降,心墙土坝的渗流(3),restpressure.savsetfluidonmechoffinifmod2e3ften0.0rangrosoilinixd0yd0zd0 xv0yv0zv0apppp40e3grad00-10e3ranz04x09apppp0ranz0 x1520histid=10zoneppid215solve,Ratio=1,Ratio=1E-5,荷载引起的地基土体的超孔隙水压力,计算文件,;-applyloadslowly-deframpramp=min(1.0,float(step)/200.0)endapplynstress=-40e3histramprangex-.13.1z9.910.1;-fluidflowmodel-modelfl_isoinifmod2e9;-porepressurefixedatzeroatthesurface-fixpp0rangez9.910.1;-settings-setfloff;-test-step75011因为本例中没有设置初始应力，这里只进行了750步的求解。,计算结果,对主从进程法的讨论,SETmechforce设置一个不平衡力的大小，达到这个不平衡力系统认为暂时达到平衡状态；SETmechsubstepnauto设置力学进程为从进程，在主进程每执行一步中必须执行n步，当系统达到平衡时也可以少于n步；SETfluidsubstepm设置流体进程为主进程。,讨论1：对收敛准则进行对比分析不平衡力（force）：1E3、5E3、1E4、5E4不平衡力比（ratio）：1E-4、1E-3、1E-2、1E-1,采用1E-3的收敛准则既可以满足计算流固耦合过程中的计算精度要求，同时又具有较高的计算效率,对主从进程法的讨论(2),讨论2：子步数的影响,setmechsub100fluidsub10setmechsub10fluidsub10setmechsub1fluidsub1,设置合理的子步数也很重要。设置过大，则会导致计算时间大大增加，过小又会造成计算结果的误差。,真空预压的简单模拟,孔压边界条件tstc长期分析(排水)Rk1骨架很软孔压扰动进行biot_mod调整,砂层,软土层,粘土层,PVD,2m,8m,10m,Datafile:,数值分析过程(movie),139,课程安排,第一讲：基本介绍、静力分析、前后处理第二讲：接触面、FISH语言、流固耦合分析第三讲：动力分析、自定义本构、结构单元第四讲：FLAC（2D）基本介绍与应用实例讨论,140,FLAC3D非线性动力分析,非常复杂！,SaidbyProf.PeterCundall,141,为什么要用FLAC做动力分析？,FLAC可以模拟体系（土，岩石，结构，流体）受到的外部动力荷载（比如地震）或内部动力荷载（比如基础振动、爆炸）。可以计算塑性引起的永久变形以及孔隙水压力的消散。土动力学中常用的等效线性方法无法直接处理上述问题。,142,动力模拟的3个重要问题,动力荷载与边界条件材料响应与阻尼土体液化,143,动力荷载,动力输入的类型加速度时程速度时程应力(压力)时程力时程APPLYINTERIOR(内部)TABLEFISH,144,Quiet边界,静态(quiet,粘性)边界LysmerandKuhlemeyer(1969)模型边界法向和切向设置独立的阻尼器性能对于法向p波和s波能很好的吸收对于倾斜入射的波和Rayleigh波也有所吸收，但存在反射人工边界仍应当足够远,145,Quiet边界应用,内部振动(如隧道中的列车振动问题)动力荷载直接施加在节点上使用Quiet边界减小人工边界上的反射不需要FF边界外部荷载的底部边界软土地基上的地震荷载不适合用加速度或速度边界条件使用应力条件t=-2Csrvs地震底部输入的侧向边界扭曲了入射波,quiet,quiet,quiet,146,Free-field边界,Cundalletal.(1980)自由场网格与主体网格的耦合粘性阻尼器，自由场网格的不平衡力施加到主体网格边界上设置条件底部水平，重力方向为z向侧面垂直，法向分别为x,y向其他边界条件在APPLYff之前,相当于一个阻尼器,147,Free-field边界,APPLYff将边界上单元的属性、条件和变量全部转移ff单元上；设置以后主体网格上的改动将不会被FF边界所响应可存在任意的本构模型以及流体耦合(仅竖向)FF边界进行小变形计算，主体网格可大变形，FF边界上的变形要相对较小存在attach的边界将不能设置FF边界边界上的Interface将不能连续动力边界设置需在FF边界设置之前,148,Free-field边界与动力荷载,模型底部边界fix施加速度或加速度荷载刚性边界Free施加应力时程荷载柔性边界对于软弱的地基不适合施加速度（加速度荷载），而应当施加应力荷载,Notethatthereisafactorof2becausetheinputenergydividesintoadownward-&upward-propagatingwave.,149,2.材料响应与阻尼,连续的非线性，表观模量随着应变的增大而降低对所有循环应变等级均存在滞回特性，因此导致随着循环应变的增加阻尼比增大。阻尼是率相关的。对于复杂波形的各个成分都产生阻尼。剪切应变会产生的体积应变，相应的，随着剪应变循环次数的增加体积应变逐渐积累。,150,材料响应,土体在循环荷载作用下呈现出模量衰减和能量消散的特点，那么如何用非线性数值方法对其进行模拟呢？,Nonlinearcharacteristicsofsoils(MartinandSeed,1979),151,试验得到的阻尼比、割线模量随循环剪应变的曲线,r,m,a,l,i,z,e,d,S,h,e,a,r,M,o,d,u,l,u,s,G,/,G,m,a,x,Mid-RangeSandCurve,(Seed&Idriss,1970),SandFillInland:,Friction=32,hr=0.47,Go=440,SandFillunderRockDike:,Friction=30,hr=0.43,Go=440,152,等效线性方法,等效线性方法是岩土地震工程中模拟波的传播的最常用的方法。假定土体是粘弹性体，参照实验室得到的切线模量及阻尼比与剪应变幅值的关系曲线，对地震中每一单元的阻尼和模量重新赋值。,Iterationtowardstrain-compatibleshearmodulusanddampingratio(afterKramer,1996),153,等效线性方法的特点,使用振动荷载的平均水平来估算每个单元的线性属性，并在振动过程中保持不变。在弱震阶段，单元会变得阻尼过大而刚度太小；在强震阶段，单元将会变得阻尼太小而刚度太大。对于不同部位不同运动水平的特性存在空间变异性。不能计算永久变形。等效线性方法模型在加荷与卸荷时模量相同，不能计算土体在周期荷载作用下发生的剩余应变或位移。塑形屈服模拟不合理。在塑性流动阶段，普遍认为应变增量张量是应力张量的函数，称之为“流动法则”。然而，等效线性方法使用的塑性理论认为应变张量（而不是应变增量张量）是应力张量的函数。因此，塑性屈服的模拟不合理。大应变时误差大。等效线性方法所用割线模量在小应变时与非线性的切线模量很相近，但在大应变时二者相差很大，偏于不安全。本构模型单一。等效线性方法本身的材料本构模型包括了应力应变的椭圆形方程，这种预设的方程形式减少了使用者的选择性，但却失去了选择其它形状的适用性。方法中使用迭代程序虽然部分考虑了不同的试验曲线形状，但是由于预先设定了模型形式，所以不能反映与频率无关的滞回圈。另外，模形是率无关的，因此不能考虑率相关性。,154,完全非线性分析方法,FLAC3D采用完全非线性分析方法，基于显式差分方法，使用由周围区域真实密度得出的网格节点集中质量，求解全部运动方程。,155,完全非线性分析方法的特点,可以遵循任何指定的非线性本构模型。如果模型本身能够反映土体在动力作用下的滞回特性，则程序不需要另外提供阻尼参数。如果采用Rayleigh阻尼或局部（local）阻尼，则在动力计算中阻尼参数将保持不变。采用非线性的材料定律，不同频率的波之间可以自然地出现干涉和混合，而等效线性方法做不到这一点。由于采用了弹塑性模型，因此程序可以自动计算永久变形。采用合理的塑性方程，使得塑性应变增量与应力相联系。可以方便地进行不同本构模型的比较。可以同时模拟压缩波和剪切波的传播及两者耦合作用时对材料的影响。在强震作用下，这种耦合作用的影响很重要，比如在摩擦型材料中，法向应力可能会动态地减小从而降低土体的抗剪强度。,156,使用弹塑性模型,附加考虑的因素：阻尼，对于屈服面以下应力的循环体积应变积累，是循环周数与幅值的函数模量衰减，基于平均应变水平的表格,最简单的弹塑性模型往往在描述累计塑性应变方面具有很好的效果，但是对于加速度放大系数的估算上效果不好。,157,弹塑性模型,简单的理想弹塑性本构模型仅仅在发生屈服时才会出现滞回特性,注意：即使这样粗糙的模型也能够作出连续的阻尼比和模量衰减曲线。,在屈服条件下会产生体积改变，但通常都是剪胀。,158,FLAC3D中的阻尼比,弹塑性本构模型使用瑞利（粘性）阻尼弹塑性本构模型使用滞后（HD）阻尼.复杂本构模型拥有连续的屈服应力应变关系和对应的加卸载响应（Wang，UBCSand）,159,瑞利阻尼,瑞利阻尼最初应用于结构和弹性体的动力计算中，以减弱系统的自然振动模式的振幅。在计算时，假设动力方程中的阻尼矩阵C与刚度矩阵K和质量矩阵M有关：,瑞利阻尼中的质量分量相当于连接每个节点和地面的阻尼器，而刚度分量则相当于连接单元之间的阻尼器。虽然两个阻尼器本身是与频率有关的，但是通过选取合适的系数，可以在有限的频率范围内近似获得频率无关的响应。,160,瑞利阻尼,frequency,ratioofdampingtocritical,混合,仅有刚度分量,仅有质量分量,采用叠加的方法得到的阻尼比在较大的频率范围内保持定值（3:1）,Combinedcurvereachesminimumat:,161,中心频率的选择,假设弹性模型计算土石坝不同材料的功率谱曲线(assumingelasticmaterial),frequency,frequency,frequency,frequency,frequency,162,阻尼比的选择,根据65%最大应变来选择阻尼比和模量衰减参参数,163,瑞利阻尼的缺点,必须指定中心频率，有时须在一些充满矛盾的数据当中进行选择（场地响应或地震平均频率）刚度比例项会导致时间步减小,164,滞后阻尼,FLACVersion5.0andFLAC3DVersion3.0提供了滞后阻尼功能，该项阻尼独立于材料模型之外，实质是包含了一个与应变幅值相关的切线剪切模量乘子,如果割线模量是由衰减曲线来确定的，那么切线模量可以表示为：,FromSeed&Idriss(1970),Go=small-strainshearmodulus,165,滞后阻尼,表观应变是偏应变，根据先前反转的点进行累计得到的。使用堆栈的数据结构可以描述应变反转现象。,因此，在“迷你滞回圈”中的能量损失与主滞回圈的计算是一样的。,elasticmodelwithhysteresisdamping,特定的曲线可以使用滞后阻尼来描述，并在本构模型计算中根据应变的大小获得不同的切线模量。,166,滞后阻尼的特点,可以直接采用动力试验中的模量衰减曲线；相对于瑞利阻尼而言，滞后阻尼不影响动力计算的时间步；可以应用于任意的材料模型，且可以与其它阻尼格式同时使用。,一个缺点：模量衰减曲线一致，但阻尼比曲线存在差别,“Good”fittoSeed&IdrisdataforG/Gmax(sigmoidal3-parameterfunction)noteinconsistentdampingresult.,G/Gmax,D-%ofcritical,167,滞后阻尼,低循环应变下得到的阻尼比要小于试验结果，这会导致低级的噪声，尤其在高频情况下。可以在中心频率上增加一个小量的Rayleigh阻尼(0.2%刚度比例)，这样也不会降低时步；若初始应力不为0，剪应力-剪应变曲线可能不匹配。因此在生成初始应力时就要调用Hyst阻尼；Hyst阻尼不仅会增加能量损失，还会导致在大循环应变下的平均剪切模量的降低，在输入波的基频接近共振频率的时候，由于可能会导致动力反应幅值的增大；Hyst阻尼之前要做一次弹性无阻尼求解，以获得发生循环应变的最大水平，若循环应变过大导致剪切模量过多的降低，那么用Hyst阻尼是有问题的；即使应变较小，使用屈服模型也会增大应变，因此若有广泛屈服的现象，则使用屈服模型，不用Hyst阻尼,168,动孔压的生成液化,干沙剪应变循环加载试验初始加载阶段，沙土通常先压实再膨胀。卸载时，沙土遵循与加载相似的路径，但在零应变时，有些残余体积应变存在。取决于初始孔隙率，这可能代表纯粹的压实假定孔隙中充满水对于常体积测试，有效应力降低，孔隙水压保持不变对于常荷载测试，(例如，盒子上法向荷载固定)，孔隙水压增加，有效应力减小有效应力为零时发生液化,169,动孔压的生成液化,因此孔隙水压增加不是液化的基本原因由于颗粒间(重组以后)的低接触力导致有效应力的减小描述液化的模型高级模型：BSHP(边界面低塑性本构模型,Wangetal.1990)简单模型：MC+体积应变增量模型Finn模型：Byrne模型：,170,UBCTOT:Post-liquefaction,Assignpost-liquefactionproperties,t,g,t,岩土介质爆炸液化模型,动孔压模式,应力/应变孔压模型,有效应力/瞬态孔压模型,内时/能量孔压模型,基于Byrne有效应力弹塑性模型，提出一个能够考虑主应力轴旋转、饱和砂土含有少量气体、饱和砂土液化后的应变软化和应力重分布特性的弹塑性模型。,模型应用大坝斜坡场地爆炸液化,图1.有效应力等值线,图2.大坝震后液化区预测,172,PL-Finn,0：零有效应力状态1：非零有效应力状态,0有效应力,非0有效应力,173,地震波的调整,基线校正对于地震分析的加速度时程，其积分得到的速度和位移应归0美国地质调查研究所BasicStrong-MotionAccelerogramProcessingSoftware(BAP)对网格施加一个固定速度从而使残余的位移变为0动力荷载的频率与单元尺寸的双向调整高频的输入要求单元尺寸很小一定的单元尺寸对应输入的最大频率一般进行滤波处理滤掉低能量的高频FFT.FISOriginSeismoSignal,174,地震波的调整,El-Centro波,FFT,修正后的时程,修正后FFT,5Hz,175,PL-Finn模型的应用,阪神地震的码头分析,沉箱顶部水平残余位移最大达5m，平均为3.5m，残余沉降为12m，海侧倾斜角35o。,seaSiltreplacedsandstonelandSandcaissonclay,FLAC3Dgrid8050zones10386gps,可液化的砂土,176,分析参数,基本力学参数,流体参数,动力参数,PL-Finn模型参数,177,输入加速度,南北向地震波最大值=0.6g竖直向地震波最大值=0.2g,178,接触面与自由边界条件,沉箱与填石之间的接触面,模型周围自由场边界条件,179,计算结果,位移,震后水平位移等值线,沉箱顶点处的水平位移、沉降时程曲线,3.44m,1.99m,180,计算结果（续）,超孔压比,震后的超孔压比云图,超孔压比时程曲线,A:置换砂B:海砂（底部）C:海砂（中部）,181,自定义本构模型,182,自定义本构模型的基本方法,必要性试验总结的本构模型特定条件下的本构模型交叉学科的本构模型二次开发环境自定义本构模型的功能自定义本构模型的基本方法,183,二次开发环境,FLAC3D采用面向对象的语言标准C+编写本构模型都是以动态连接库文件(.DLL文件)的形式提供VC+2005或更高版本的开发环境优点自定义的本构模型和软件自带的本构模型的执行效率处在同一个水平自定义本构模型(.DLL文件)适用于高版本的FLAC(2D)、3DEC、UDEC等其他Itasca软件中,184,自定义本构模型的功能,主要功能：对给出的应变增量得到新的应力辅助功能：模型名称、版本读写操作模型文件的编写基类(classConstitutiveMod