在ansys中导入自定义本构模型---邓肯-张模型
- 格式:doc
- 大小:27.50 KB
- 文档页数:5
两种本构模型的对比分析
土体本构模型的研究和探讨由来已久,但考虑到土的复杂特性,当其受到外部荷载的作用之后,其变形具有剪胀性、非线性、各向异性、流变性等特性,因此目前所发展的几百个本构模型在实际应用的过程中普遍存在局限性。
为了进一步推动土的本构理论的发展,本文主要从以下几个方面着手进行研究。
第一,系统介绍了目前岩土工程领域广泛应用的本构模型——邓肯-张模型,该模型能够较为完整地描述土体的主要变形特征,并且模型参数也都具有明确的物理意义,一般通过常规的三轴试验来确定其邓肯-张模型的参数。
因此,邓肯-张模型在岩土工程中受到了广泛的应用。
第二,本文引入卢应发教授所提出的新的本构模型,并且对新的本构模型的基本方程、参数的决定、相关系数以及模型特征进行了详细介绍。
新的本构模型阐明了当前在岩土工程领域中现有的众多类型的本构模型的数学基础和假设及其相关联系,针对具有极其复杂力学特性的岩土材料,新的本构模型为建立统一的岩土理论模型提供了新的尝试。
第三,本文分别用邓肯-张模型和新的本构模型对同一试验算例进行计算,两种本构模型的计算结果均取得了令人满意的效果。
最后,本文在有限元软件ABAQUS的平台中实现了邓肯-张模型和新的本构模型的二次开发,利用本文所开发的邓肯-张模型和新的本构模型对同一试验进行模拟分析,验证两种本构模型的适用性,从而为土体本构模型的研究找到新的方法。
CAD导入ANSYS的方法ANSYS是一款强大的工程仿真软件,支持多种文件格式,包括CAD格式。
DXF(Drawing Exchange Format)是一种CAD文件格式,可以被许多CAD程序使用。
在ANSYS中将DXF文件导入可以通过以下步骤实现。
1.准备要导入的DXF文件。
确保DXF文件已经完成设计和准备好导入到ANSYS中。
如果需要作出任何修改,请在导入之前完成。
2.打开ANSYS软件。
启动ANSYS工程仿真软件,确保已经加载了所有必要的模块和组件。
3.创建新工程。
在ANSYS的主界面上,选择“创建新工程”或类似的选项,创建一个新的工程。
4.导入DXF文件。
在ANSYS工程界面上,选择“文件”->“导入”->“三维CAD模型”(文件路径和选项可能会有所不同,具体根据软件版本而定)。
在导入对话框中,找到DXF文件所在的位置,并选择它。
然后,点击“打开”按钮。
5.设置导入选项。
在导入对话框中,你可以设置导入选项,如坐标系、单位、长度、尺寸等。
根据DXF文件的属性和设计需求,选择适当的选项。
6.等待导入完成。
ANSYS将开始导入DXF文件,并在导入过程中显示进度条。
等待导入过程完成,不要中途中断操作。
7.检查模型。
在导入完成后,ANSYS将显示导入的DXF文件。
你可以使用ANSYS的多种功能和工具,如缩放、旋转、平移、选择等,来查看和检查导入的模型。
如果需要,可以进行进一步的调整和修改。
8. 保存模型。
在导入的模型得到满意的结果后,选择“文件”->“保存”(或类似选项),将模型保存到ANSYS项目文件中(.wbpj格式),以便以后的使用和分析。
总结:DXF格式是通用的CAD文件格式,可以被许多CAD程序使用。
在ANSYS中将DXF文件导入,需要先创建新工程,然后选择导入选项并指定DXF文件,等待导入完成后,使用ANSYS的功能和工具进行进一步的检查和修改,最后保存模型。
导入DXF文件到ANSYS中可以用于许多工程仿真应用,如结构分析、流体动力学、电磁场等。
Adams-Ansys单位制验证为了在adams做模态分析,就要用到ansys制作mnf文件。
两个软件之间存在单位转换的问题,于是,我用adams做了一个曲柄,采用的单位制是mm-kg-s,长度100mm,宽度20mm,厚度10mm,保存成iges格式,导入到ANSYS中。
在ANSYS中,首先创建了Brick 8node 45单元,然后设置弹性模量和密度的时候,要特别注意了:采用第一种:采用mm-t-s单位(大侠讲:这是工程问题,由于工程上常用mm做长度单位,被迫质量采用吨,这样单位才会封闭。
)1、钢的弹性模量是2*e+005 MPa,根据换算关系,弹性模量就应该是2*e+011Pa。
钢的密度是7810kg/m3,转换成7.81*e-9t/mm3。
泊松比是0.3。
2、划分网格,设置外连点。
在导出到adams中之前,单位应采用user defined,所有单位系数都是unit,user,1000,0.001,1,,,,1。
3、导入adams后发现,匹配成功。
采用第二种:采用m-kg-s单位。
这种单位制的转化在ANSYS前处理中,一直采用m-kg-s单位,1、钢的弹性模量是2*e+011MPa,钢的密度是7801kg/m3。
泊松比是0.3。
中间步骤与上述一致,省略了。
直到输出mnf文件的时候,设置单位user defined,/unit,user,0.001,1,1,,,,1。
(之所以这么设置,我的初衷是因为我在adams中建模用的mm,而在ansys中建模却是m,如果不把长度单位缩小,图形肯定不一致)。
相当于下图操作。
结果证明,第一种方法在图形大小和柔性体质量方面是合适的,第二种方法还需要解决,图形和质量都非常大。
存在疑问:第一种方法的质量差别较小,但转动惯量有较大差别,怎么办?。
ANSYS 理论基础一、钢筋混凝土模型1、Solid65单元——模拟混凝土和岩石等抗压能力远大于抗拉能力的非均匀材料开发的单元,可以模拟混凝土中的加强钢筋(或玻璃纤维、型钢等);普通8节点三维等参元,增加针对混凝土材料参数和整体式钢筋模型;基本属性:——可以定义3种不同的加固材料;——混凝土具有开裂、压碎、塑性变形和蠕变的能力;—-加强材料只能受拉压,不能承受剪切力。
三种模型:分离式模型——把混凝土和钢筋作为不同的单元来处理,各自划分单元,或钢筋视为线单元(杆件link-spar8或管件pipe16,20);钢筋和混凝土之间可以插入粘结单元来模拟界面的粘结和滑移;整体式模型——将钢筋分布于整个单元中,假定混凝土和钢筋粘结很好,并把单元视为连续均匀材料;组合式模型—-分层组合式:在横截面上分成许多混凝土层和若干钢筋层,并对截面的应变作出某些假设(如应变沿截面高度为直线);或采用带钢筋膜的等参单元。
2、本构模型线性弹性、非线性弹性、弹塑性等;强度理论——Tresca、V on Mises、Druck —Prager等;3、破坏准则单轴破坏(Hongnested等)、双轴破坏(修正的莫尔库仑等)、三轴破坏(最大剪应力、Druck—Prager等),三参数、五参数模型;混凝土开裂前,采用Druck—Prager屈服面模型模拟塑性行为;开裂失效准则,采用William-Warnke五参数强度模型.4、基本数据输入混凝土:ShrCf-Op—张开裂缝的剪切传递系数,0~1ShrCf—Ol—闭合裂缝的剪切传递系数,0。
9~1UnTensSt—抗拉强度,UnCompSt—单轴抗压强度,(若取-1,则以下不必要)BiCompSt—双轴抗压强度,HydroPrs—静水压力,BiCompSt—静水压力下的双轴抗压强度,UnCompSt-静水压力下的单轴抗压强度,TenCrFac—拉应力衰减因子。
加固材料(材料号、体积率、方向角)二、其他材料模型在Ansys中,可在Help菜单中查阅各种不同单元的特性.例1、矩形截面钢筋混凝土板在中心点处作用-2mm的位移,分析板的受力、变形、开裂(采用整体模型分析法).材料性能如下:1、混凝土弹性模量E=24GPa,泊松比ν=0。
UPFs教程1Pomato157300@中国矿业大学(徐州)力学与建筑工程学院:师访1.简单介绍UPFs(User ProgrammableFeatures,用户可编程接口)是用户在ANSYS提供的Fortran和C 语言源代码(并非真正的源代码,仅仅是ANSYS公司定义好的子程序函数)的基础上,修改其用户可编程子程序和函数,从源代码的层次上扩充ANSYS的功能。
用户需要在相应的Fortran和C语言编译器的支持下,将编译修改后的源代码与ANSYS库相连接形成用户版本的ANSYS可执行文件。
所谓的“用户可编程子程序和函数”,以ANSYS12.0版本为例,在安装时不默认安装,需要选择自定义安装方式,并选择ansys customization files。
安装后这些子程序和函数存在于C:\Program Files\ANSYS Inc\v120\ansys\customize\user目录下,多为*.F,即fortran 文件,用户可对其进行修改,拷贝到相关文件夹并编译连接。
2.利用UPFs可以完成的工作:(1)对ANSYS数据库进行操作。
(2)定义特定的荷载类型。
(3)定义新的单元类型。
(4)在非线性功能中实现对用户实参的运算。
(5)接触问题中实现用户摩擦系数。
(6)定义特定的非线性材料特性———如用户塑性屈服准则、蠕变方程、超弹性特性、失效准则、粘弹性性质等。
(7)用户优化———可以用自己的算法和中断准则替换ANSYS优化过程。
(8)定义用户介入计算过程的方式———在每个ANSYS运行求解、载荷步、子步和平衡迭代的开始和结束处允许用户介入,并允许用户在求解过程中评估计算结果。
3.常用的子程序及函数(1)用户自定义命令user01.F-user10.F用户可以使用它们自定义ANSYS命令在ANSYS中执行usr1(或user2)等命令调用以上user01.F(user02.F)中定义的功能。
这里也可以使用/UCMD命令为user01.F指定一个新名称,再通过该名称执行命令。
ANSYS、FLENT、ALGOR、ROBOT、CAESAR II、STAAD PRO、3D3S的区别ANSYS和FLOTRAN1、两个软件的关系:ANSYS公司2005年收购fluent,如今在ansys12版本中已集成fluent2、两个软件使用方向不一样:ANSYS用于固体力学,FLUENT专用于流体力学3、ANSYS的FLOTRAN流体模块是基于有限元方法,FLUENT则是基于有限体积法4、对于机械方向,除了流体机械专业,其他专业更多的使用有限元,也就是说,使用ANSYS更多一些。
流体机械专业则两者都要使用。
ansys侧重于固体传热和应力应变分析等,在求解流体问题是,没有FLUENT好使,airpak主要用于气流组织的模拟,可以作为FLUent的前处理软件fluent专门做流体分析的,热流之类的机械和建筑方向肯定要用ansys流动传热的问题建议用fluent,纯导热问题用ansys。
因为ansys的热分析模块只能处理纯传热问题,不计算流场。
ansys的强项在于处理固体问题,流体有关的问题不是它的擅长,是fluent的擅长。
所以楼主的问题属于固壁传热问题,原来的ansys就可以较好的解决。
另外,虽然说ansys先后收购了CFX和fluent这两个软件,但是ansys仅是指ansys本身的软件,不包括上面的两个。
个人看法:暂时不会出新的混合了上述三种的所谓的新ansys,因为从算法上讲,ansys用的是有限元算法,而fluent和CFX用的是有限体积法,所以暂时无法整合到一起。
ALGORALGOR是新一代的CAE分析工具,在汽车、电子、航空航天、医学、军事、电力系统、石化、土木工程、微机电系统、日用品生产等诸多领域中均得到了广泛的应用。
ALGOR核心代码起源于1970年开发的SAP程序,它是由美国加州大学伯克利分校的K.J.Bathe、E.L.Wilson和F.E.Peterson等人共同研制。
A N S Y S中混凝土的本构关系-CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN一、关于模型钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种,即分离式、分布式和组合式模型。
考虑钢筋和混凝土之间的粘结和滑移,则采用引入粘结单元的分离式模型;假定混凝土和钢筋粘结很好,不考虑二者之间的滑移,则三种模型都可以;分离式和分布式模型适用于二维和三维结构分析,后者对杆系结构分析比较适用。
裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型,后者目前尚处研究之中,主要应用的是前两种。
离散裂缝模型和分布裂缝模型各有特点,可根据不同的分析目的选择使用。
随着计算速度和网格自动划分的快速实现,离散裂缝模型又有被推广使用的趋势。
就ANSYS而言,她可以考虑分离式模型(solid65+link8,认为混凝土和钢筋粘结很好,如要考虑粘结和滑移,则可引入弹簧单元进行模拟,比较困难!),也可采用分布式模型(带筋的solid65)。
而其裂缝的处理方式则为分布裂缝模型。
二、关于本构关系混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类,其中研究最多的是非线性弹性和弹塑性本构关系,其中不乏实用者。
混凝土破坏准则从单参数到五参数模型达数十个模型,或借用古典强度理论或基于试验结果等,各个破坏准则的表达方式和繁简程度各异,适用范围和计算精度差别也比较大,给使用带来了一定的困难。
就ANSYS而言,其问题比较复杂些。
1 ANSYS混凝土的破坏准则与屈服准则是如何定义的?采用tb,concr,matnum则定义了W-W破坏准则(failure criterion),而非屈服准则(yield criterion)。
W-W破坏准则是用于检查混凝土开裂和压碎用的,而混凝土的塑性可以另外考虑(当然是在开裂和压碎之前)。
理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的,例如在高静水压力下会发生相当的塑性变形,表现为屈服,但没有破坏。
ANSYS中混凝土的计算问题最近做了点计算分析,结合各论坛关于这方面的讨论,就一些问题探讨如下,不当之处敬请指正。
一、关于模型钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种,即分离式、分布式和组合式模型。
考虑钢筋和混凝土之间的粘结和滑移,则采用引入粘结单元的分离式模型;假定混凝土和钢筋粘结很好,不考虑二者之间的滑移,则三种模型都可以;分离式和分布式模型适用于二维和三维结构分析,后者对杆系结构分析比较适用。
裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型,后者目前尚处研究之中,主要应用的是前两种。
离散裂缝模型和分布裂缝模型各有特点,可根据不同的分析目的选择使用。
随着计算速度和网格自动划分的快速实现,离散裂缝模型又有被推广使用的趋势。
就ANSYS而言,她可以考虑分离式模型(solid65+link8,认为混凝土和钢筋粘结很好,如要考虑粘结和滑移,则可引入弹簧单元进行模拟,比较困难!),也可采用分布式模型(带筋的solid65)。
而其裂缝的处理方式则为分布裂缝模型。
二、关于本构关系混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类,其中研究最多的是非线性弹性和弹塑性本构关系,其中不乏实用者。
混凝土破坏准则从单参数到五参数模型达数十个模型,或借用古典强度理论或基于试验结果等,各个破坏准则的表达方式和繁简程度各异,适用范围和计算精度差别也比较大,给使用带来了一定的困难。
就ANSYS而言,其问题比较复杂些。
1 ANSYS混凝土的破坏准则与屈服准则是如何定义的?采用tb,concr,matnum则定义了W-W破坏准则(failure criterion),而非屈服准则(yi eld criterion)。
W-W破坏准则是用于检查混凝土开裂和压碎用的,而混凝土的塑性可以另外考虑(当然是在开裂和压碎之前)。
理论上破坏准则(failure criterion)和屈服准则(y ield criterion)是不同的,例如在高静水压力下会发生相当的塑性变形,表现为屈服,但没有破坏。
ANSYS导出柔性体MNF⽂件⼊ADAMS的详细步骤ANSYS导出柔性体MNF⽂件⼊ADAMS的详细步骤基本思路:在ANSYS中进⾏模态中性⽂件(.mnf)⽂件的输出,然后把输出的.mnf ⽂件输⼊到ADAMS中,进⾏零件更换。
最后在ADAMS中进⾏加载约束,仿真,查看结果。
建模仿真软件:ANSYS10.0,ADAMS2013。
具体步骤:1. ANSYS输出.mnf柔性⽂件1.1 准备⼯作和导⼊模型a、设置ANSYS⼯作⽬录,定义⼯作⽂件名和⼯作标题。
b、在其它三维软件如SolidWorks中创建⽬标模型,并将模型另存为*.x_t格式。
将*.x_t⽂件拷贝⾄上述设置的ANSYS⼯作⽬录中,c、File→Import→PARA...→,弹框左侧框中就会看到刚才⽣成的*.x_t⽂件,选中,OK,导⼊完成。
备注:⽂件名和⼯作⽬录最好没有中⽂字符出现。
d、默认模型是线框模式,PlotCtrls→Style→Solid Model Facets→下拉框中选择Normal Faceting→OK, Plot→Replot,即可看到实体模式了【可不进⾏此操作】。
1.2显⽰关键点和按⽐例放缩模型a、显⽰关键点:PlotCtrls→Numbering…弹出以下界⾯,勾选keypoint相关项查看关键点坐标信息:List→Keypoint→Coordinates Only,弹出以下界⾯。
根据关键点坐标信息可知,模型的长度单位为mb、此步可选许可不选,只要单位统⼀即可。
由于在ADAMS中长度单位⽤的是mm,⽽且本实例中其他参数也是按mm给定,如密度、弹性模量等。
将模型进⾏缩放处理,如下所⽰1.3创建单元和材料a、单元1:Solid(Brick8node185)或者其他3D单元;单元2:StructuralMass(3Dmass21),此单元只⽤于连接点单元。
b、设置材料属性:密度,弹性模量,泊松⽐3个参数,以N,mm,kg,s作单位,EX 为2.1e5,PRXT为0.3,DENS为7.85e-6。
亚克力的lsdyna本构模型亚克力是常用的一种塑料材料,具有很好的抗冲击、透明度和质感,广泛应用于家居、建筑和电子产品等领域。
在工程设计和模拟中,了解亚克力的材料本构模型是非常重要的。
本构模型是用来描述材料行为的数学模型,通常包括弹性模量、屈服强度、塑性应变等参数。
亚克力的本构模型可以使用LS-DYNA软件进行建模和模拟,下面将按照步骤阐述亚克力的LS-DYNA本构模型。
第一步:在LS-DYNA软件中新建模型,选择亚克力作为模拟材料,并设置模型的尺寸和几何结构。
第二步:选择亚克力的本构模型,常用的有线性弹性模型、可退化弹性模型、Von Mises本构模型和Drucker-Prager本构模型等。
其中,Drucker-Prager本构模型可以更真实地描述亚克力的力学行为,建议选择该模型进行模拟。
第三步:设置材料参数,包括杨氏模量、泊松比、屈服强度、塑性应变等。
材料参数的设置需要参照实验数据或文献资料,以确保模型的准确性。
第四步:设置加载条件,包括加载速度、初始速度、加载方向等。
在亚克力的模拟中,通常会考虑冲击加载、剪切加载和拉伸加载等情况。
第五步:运行模拟,观察模拟结果,并进行修正和验证。
在模拟过程中,需要注意模型的收敛性、速度和应力变化等情况,以减少误差。
在以上步骤中,亚克力的本构模型选择和材料参数设置是比较关键的环节。
如果参数设置不正确,会导致模型的准确性和可靠性下降,因此需要进行反复修改和验证。
同时,需要提醒的是,本构模型只是对亚克力材料行为的近似描述,实际情况可能存在误差和偏差。
总之,了解亚克力的LS-DYNA本构模型对于工程设计和模拟是非常重要的。
通过选择合适的本构模型、设置准确的材料参数以及优化加载条件,可以更加真实地反映亚克力材料的力学行为,提高模拟结果的准确性和可信度。
ansys中两种方法给材料添加材料属性1 第一种在划分网格之前指定1.1 main menu/preprocessor/meshing/mesh attributes/default attribs 出现meshing attributes 对话框,在【mat】material number下拉框中选择你需要的材料序号。
单击ok1.2 然后划分网格,则此次划分的网格的材料属性为选择的材料序号的属性。
2 第二种在划分网格之后指定2.1 先划分好网格2.2 点击select/entities/ 第一项选择areas ,第二项选择by num/pick,然后点击ok ,弹出面积选择框,选定面积,点击ok,完成面积选择2.3 点击select/entities,第一项选择elements,第二项选择attached to ,第三项选择areas,表示所要选择的单元为已选定面积中的单元,点击ok,选中面中的所有单元。
2.4 点击plot/replot,将只显示已选定的单元和面积。
2.5 点击main menu/preprocessor/material pros/change mat num,在new material number 文本框中输入你需要的材料序号,在elements No. to modefied 输入all 表示所选定的所有单元对应的材料属性转为此材料属性。
ansys多种材料怎样设置材料属性呀,用什么命令?GUI方式楼上正解,或者Proprecessor->Meshing->Mesh Attributes->Picked Volumes 命令为VATT本人喜欢在划分单元前先选好材料、实常数等再划分,命令流如下:type,1mat,1real,11vmesh,all对于其他不同材料,方式相同:ansys中的等效应力是什么物理含义?它与最大应力s1有什么区别,平常讨论应力分布,应该用等效应力还是最大应力s1呢?1)计算等效应力时是否需要输入等效泊松比呢?好像有效泊松比的默认值是0.5。
A N S Y S导入F L A C D全过程标准化工作室编码[XX968T-XX89628-XJ668-XT689N]从C A D导入A N S Y S再到F L A C3D的全过程?很多朋友在学习ANSYS时,发现它是一个大型通用软件,很多模型都可以通过它进行数值模拟。
但就是因为它的通用性,在涉及到具体各个专业领域时就显得不那么精确。
对于岩土工程来说,ANSYS中只有D-P准则,这对于岩土的本构模型是远远不够的。
我们就不得去学习相对专业性强一点的软件,这时FLAC3D就进入了我们的视线。
FLAC3D几乎拥有了岩土的所有的本构模型,并且在模拟岩土工程中的流固耦合有相当强的优势。
当大家兴致勃勃的学习FLAC3D,发现它的前后处理能力比较弱,尤其时建立网格相当复杂。
如果碰到复杂的边界时,大家就只能望而却步了,十分打击学习FLAC3D的积极性,大家可能开始怀念ANSYS的超强前处理能力了。
这时,就不禁要问,能不能把建好的ANSYS网格导入FLAC3D中,然后再在FLAC3D里面进行计算呢?有,当然有。
下面我就分享一下我学习的过程,解密从CAD导入ANSYS再到FLAC3D的全过程。
有点啰嗦了,开始正题。
下面我将从一个隧道的建模开始,当然对于其他的岩土工程也同样适用。
首先从CAD里建立隧道的平面模型,值得注意的时,在CAD中建立的线与线之间交叉必须有断点,否则导入ANSYS中将会出现问题。
下面是隧道的CAD图形:在CAD里面画好隧道模型的平面图后,需要保存成CAD2004的DXF文件格式,不然无法导入ANSYS中。
我使用的CAD导入ANSYS的软件是Dxf2Ansys测试版本,转换前,把保存好的DXF文件拷贝到Dxf2Ansys测试版本文件夹中,然后就可以点击Dxf2Ansys生成导入的ANSYS命令流文件。
下面是导入的过程:点击Writeansysfile就完成生成命令流格式。
然后将命令流复制拷到ANSYS中,在ANSYS中就生成拥有keypoints和lines的模型。
一、Ansys中建模导入aba qus 中虽然很早就看到有人在ANS YS建模,然后得出几何模型数据,再进入到ABA QUS软件进行计算,但是过程相对来说比较繁琐一些。
某些设计院或者研究人员自己也编制了一些导入程序,但都属于内部使用。
前几天看到AB AQUS 6.10版本可以实现ANSYS模型的直接导入,便尝试了一下,发现效果还是很不错的。
基本流程如下:1. 首先在ANSY S模型中输入cdb格式的文件。
(下面是世博会中国的ANSY S模型)我用的输出的方法是:命令行中输入cdwrite,all,filenam e,cdb,然后在文件下便出现我们所需要的文件。
可能还有其他办法,我没有尝试。
2. 打开ABAQU S 6.10文件中的i mport命令,导入模型即可。
3. 在ABAQUS中即可显示出导入后的模型。
看效果还是很不错的。
然后查看信息,包含的还是很全的。
在导入的时候很顺利,没有遇到什么情况,但我对导入后的模型没有没有做详细的检查,大家有兴趣了可以自己试试看。
二、请问各位大侠:请问ABAQU S后处理云图中,其图例的字体大小怎样调整?谢谢。
如果我没有理解错的话。
viewpor t//viewpor t annotat ion options//进入viewp ort annotat ion options对话框,选择Legen d标签,在Text区域选择SetFront按钮进行设置。
inp编码介绍(一).ABAQUS头信息文件段(1-4)1.*PREPRIN T 输出求解过程所要求的信息(在dat文件中)ie:*PREPRIN T, ECHO=YES, HISTORY=YES, MODEL=YES2.*HEADING标题输出文件(出现在POST/VIEW窗口中,且出现在结果输出文件中)ie:*HEADINGSTRESSANALYSI S FOR A PLA TE WITH A HOLE3.*RESTART要求abaqu s/standar d输出其PO ST/view模块所需要的.res文件。
邓肯-张模型的关键点是材料的弹性模量随大小主应力差及小主应力(围压)的变化而变化,用
APDL实现之的基本思路是:给每个单元定义一个材料号,分级施加荷载,在每个荷载步结束
时提取出各单元的大小主应力,据此计算出下个荷载步的弹性模量Et,修改各单元之MP,用
于下一步计算。
以下是一个简单算例,copy出去可直接运行。
!!!常规三轴试验模拟
!!!by taomingxing,NWPU
!!!2003.7.16
FINISH
/CLEAR
/TITLE,Numerical Simulation of three axes testing of soils
/PREP7
*dim,SUy,array,50 !Settlement records
*dim,MaxPs,array,120 !Max history p1-p3
*dim,MaxDs,array,120 !Max history Ds
!*dim,EEt,array,50 !Et of elememt
!!!Duncan-Chang Model
!!!Symbols:c-粘滞力,Fai-内摩擦角,Sf-破坏强度(p1-p3)f,Ds-应力水平,Pa-大气压,P3-
围压
*CREATE,Duncan-Chang !Creat Macro file
*afun,deg !Unit of angle
*set,Pa,1e5
*set,P1,-ArrS3(i) !注意:岩土工程中应力为拉负压正
*set,P3,-ArrS1(i)
*if,P3,LT,0.1*Pa,then
P3=0.1*Pa !围压最小取值
*endif
Sf=2*(c*cos(Fai)+P3*sin(Fai))/(1-sin(Fai)) !Mohr-Coulomb破坏强度(p1-p3)f
Ds=(P1-P3)/Sf !应力水平,
*if,Ds,GT,0.95,then
Ds=0.95 !应力水平最大取值
*endif
!判断加卸荷,如果(P1-P3)小于历史最大值视为卸荷-再加荷过程
*if,MaxPs(i),LT,P1-P3,then
Ei=k*Pa*(P3/Pa)**n
Et=Ei*(1-Rf*Ds)**2 !加荷情况的切线模量
MaxPs(i)=P1-P3 !保存历史最大应力
*elseif,MaxPs(i),GE,P1-P3
Et=Kur*Pa*(P3/Pa)**n !卸荷模量
*endif
mp,ex,i,Et !修改单元i的Et
mp,nuxy,i,Mu
*END
!!!单元类型
et,1,42 !平面四节点单元
KEYOPT,1,3,2 !平面应变
!!!以下定义材料初始模量
mp,ex,1,3.728e7 !砂土的弹性模量
mp,nuxy,1,0.33
mp,dens,1,1800
!!!建立几何模型
blc4,0,0,0.08,0.15 !8cm X 15cm
/PNUM,AREA,1
/REPLOT
!!!网格划分
aesize,all,0.01
mat,1
amesh,all
!!!边界条件
nsel,s,loc,y,0
d,all,Uy !底边界竖向约束
nsel,s,loc,x,0
d,all,Ux !左侧边界水平向约束
nsel,all
/replot
fini
/SOLU
time,0.01 !施加围压
sfl,all,pres,2e5 !200kPa
solve
!!!分级施加荷载,实现非线性计算
!!!荷载增量10kPa,共50级
*DO,ti,1,50
!!!取出计算结果,修改弹性模量
/POST1
*get,SUy(ti),node,29,u,y !Settlement record of time ti
ETABLE,EtabS1,S,1 !取各单元第一主应力
ETABLE,EtabS3,S,3 !取各单元第三主应力
*dim,ArrS1,array,120
*dim,ArrS3,array,120
*do,Num,1,120 !Num为单元编号
*get,ArrS1(Num),elem,Num,etab,EtabS1 !将单元结果存入数组
*get,ArrS3(Num),elem,Num,etab,EtabS3
*enddo
/PREP7
!^^修改砂土单元的Et,单元号1-120
c=0 $Fai=35 $Rf=0.7 $k=400 $n=0.6 $Mu=0.33 $Kur=326.7
*do,i,1,120 !各单元循环计算
*use,Duncan-Chang,c,Fai,Rf,k,n,Mu,Kur !调用Duncan-Chang宏文件
*enddo
!EEt(ti)=ET !保存第120单元之ET
/SOLU
time,ti
sfl,3,pres,2e5+1e4*ti !施加荷载,增量1e4
solve !对ti级荷载情况求解
*ENDDO