Copyright © 2020 生辰八字网 版权所有
免责声明：本站部分资源来源于网络，如有侵权，请与联系我们管理员联系，我们将会尽快处理，谢谢
我也是新手这些资料你看看吧！都是在网上找的，看完应该能明白个大概！初始地应力平衡应该满足：（1） 竖向位移为零或者很小，一般至少10的负6次方，甚至更小 （2） 竖向应力s22要等于土体深度*密度*g 关于地应力的平衡方法，综合了版上的一些意见，结合了自己的想法，对于初始地应力的施加，得到了e-6的效果，方法比较简单，与大家分享！1.先施加重力荷载的作用，可以在cae中实现；2.在inp文件中的output request中写上*el print s，这样就会将施加重力荷载后的应力输出到*.dat文件中了；3.在*.dat文件中，将单元应力的序号及单元的应力拷出，例如 ELEMENT T FOOT- S11 S22 S33 S12 NOTE 1 1 -1.2598E+05 -1434. -3.1852E+04 892.7 2 1 -1.2249E+05 -6287. -3.2194E+04 1223. 3 1 -1.1795E+05 -497.7 -2.9611E+04 1664. 4 1 -1.1210E+05 -7240. -2.9834E+04 1992. 5 1 -1.0485E+05 579.0 -2.6068E+04 2600. 6 1 -9.5803E+04 -8272. -2.6019E+04 3031. 7 1 -8.4709E+04 1915. -2.0698E+04 4083. 8 1 -7.0634E+04 -9746. -2.0095E+04 4339. 9 1 -5.1088E+04 5401. -1.1422E+04 8519. 10 1 -2.4353E+04 -1.1150E+04 -8876. 1.2126E+04 11 1 -1.2847E+05 268.1 -3.2050E+04 738.1 12 1 -1.2786E+05 -9868. -3.4433E+04 629.1 13 1 -1.2938E+05 -4224. -3.3402E+04 502.5 14 1 -1.3039E+05 -3458. -3.3461E+04 165.9 单独存为一个*.dat文件，4.用excel打开该文件，将其中的1所在的列去掉，在每个单元号前面加上其instance. ，即单元编号变为： instance名称.序号 ；注意不同的instance和part要都按照其所在的单元从小到大编号，而不是按照他们在整体单元编号来编号！5.接下来就在excel把该文件另存为*.csv格式的文件（即带有逗号分隔符的格式），6.最后在inp文件的step之前写上*initial conditions,type=stress,input=文件名.csv即可！这种方法不需要用python，比较简单，希望能对大家有用！先说为什么要施加地应力：1、我们所建立的几何模型一般和工程实际情况或尺寸相对应、相一致，比如边坡几何模型和实际边坡尺寸一致，但我们可以夸张一点想像，实际边坡应是由一个更大一点或更高一点的不受重力的初始边坡在n年前突然受重力和类似目前的边界条件作用下逐渐形成了今天的尺寸大小，n年前受重力和类似目前的边界条件作用之前边坡的尺寸大小，我们不得而知，如果能准确知晓，我们就可以建立一个那时的几何模型，再施加重力和边界条件进行计算，变形后形状和现状边坡形状一致，其内力也就是初始应力场或地应力，就不用专门去施加地应力了，但问题是我们不能知晓边坡受力前的形状尺寸，我们现在的几何模型就是边坡现在的实际尺寸，受力后将会变成一个更小的或与现状不一致的边坡，这不符合我们模拟现状边坡的目的。
如果我们知道现状边坡的内力，将其提取出来作为几何模型的内力，再和外力（重力）平衡，则我们建立的模型才能算和实际模型一致。 真实地知道现状边坡的内力是很难的，我们采取的办法是，用我们所建立的几何模型施加和实际模型一致的重力和边界条件进行计算，得到变形后或变得更小或与现状边坡不完全一致的边坡内力近似的作为现状边坡的内力，并重新将其施加于与现状边坡一致的几何模型，再施加重力（当然边界条件也应基本一致）以平衡，这样才算建立了与现状模型基本一致的模型，其下的计算才成为可能。
这就是所谓“地应力平衡”的含义、目的、作用。2. 地应力平衡中的外力和内力的问题，地应力平衡中，显然，重力是外力，应力场是内力，仅有外力重力，没有内力是不可能的，同样，仅有内力（专指初始应力场）而不受重力也是不可能的，否则，整个体系的力不会平衡。
这就是为什么我们将提取出的内力施加于几何模型后必须再施加重力的原因。为的是内力和外力平衡。）
q0 F3 q6 H1 O#3. 地应力场的方向问题，有网友在论坛里问，既然重力是向下，为与重力平衡，那应力场的方向是不是向上呢，这同样是我开始接触abaqus的疑问，相信很初学者也有这样的疑问，我的理解是内力是没有向上、向下或者向其它方向的概念的，内力只有拉力或压力或剪力之分，其方向也按是拉是压是顺时针或逆时针而分，内力往往都是成对出现，如地应力场中的应力以压应力为主，取一个微元，则压应力同时出现在向下和向上，你能说地应力就是向上，与重力反向吗？ 不怕各位笑话，以上几点在高手看来是很简单的问题，却是我经历了漫长而艰辛的摸索才得到的，今天也写给初学者，不要再走我的老弯路了。 aba中初始地应力场平衡一般在表面水平的情况下仅仅和密度相关，密度一样的话平衡的结果很好，别的参数改变之后经过计算，差别很小。
表面不水平的情况则最好通过文件导入初应力的情况进行平衡。一般来讲，表面不平的时候有很多因素造成误差很大。
地应力平衡方法：第一步：建立模型，材料，分析步（GEOSTATIC） 第二步：施加荷载，LOAD，选择施加重力GRAVITY，在你想施加重力的方向输入数值9.8 第三步：在命令行中输入mdb.models['模型名字'].setValues(noPartsInputFile=ON) （请严格按照 。
法1：通过开始菜单的Abaqus Command运行： a.此方法为启动ABQUS CAE的同时运行脚本，也可以用abaqus cae startup=Beam.py达到同样的效果，与此类似abaqus viewer script=Beam.py 或者 abaqus viewer startup=Beam.py可以在启动ABAQUS Viewer的同时运行脚本。 b.与上述方法相对应，当脚本文件中包含后自动前后处理的相关命令，那么就没有启动ABAQUS CAE/Viewer必要，否则徒增加计算成本，相应的命令为abaqus cae noGUI=beam.py或者abaqus viewr no GUI=beam.py，但此方法的不足之处是脚本在执行的过程中无法与外界进行交互。 法2.在启动屏幕时运行脚本
法3. 从File菜单运行，即File菜单下Run Script的子菜单 法4.从命令行接口运行脚本 在命令行中输入命令：execfile('file_name') 例如：execfile('Beam.py')
abaqus产生几类文件： 1. model_database_name.cae 模型信息、分析任务等。
2. model_databse_name.jnl 日志文件：包括用于复制已存储模型数据库的abaqus/cae命令 *.cae和*.jnl构成支持CAE的两个重要文件，要保证CAE下打开一个项目，这两个文件必须同时同在； 3. job_name.inp 输入文件。由abaqus Command支持计算起始文件，它也可由CAE打开； 4. job_name.dat 数据文件：文本输出信息，记录分析、数据检查、参数检查等信息。
ABAQUS/Explicit的分析结果不会写入这个文件 5. job_name.sta 状态文件：包括分析过程信息 6. job_name.msg 是计算过程的详细记录，分析计算中的平衡迭代次数，计算时间，警告信息，等等可由此文件获得。用STEP模块定义 7. job_name.res 重启动文件，用STEP模块定义 8. job_name.odb 输出数据库文件，即结果文件，需要由visuliazition打开 9. job_name.fil 也为结果文件，可被其他应用程序读入的分析结果表示格式。
ABAQUS/Standard记录分析结果。ABAQUS/Explicit的分析结果要写入此文件中则需要转换，convert=select 或 convert=all 10.abaqus.rpy 记录一次操作中几乎所有的ABAQUS/CAE命令 11.job_name.lck 阻止并发写入输出数据库，关闭输出数据库则自行删除 12.model_database_name.rec 包含用于恢复内存中数据库的ABAQUS/CAE命令 13.job_name.ods 场输出变量的临时操作运算结果，自动删除 14.job_name.ipm 内部过程信息文件：启动ABAQUS/CAE分析时开始写入，记录了从ABAQUS/STANDARD或ABAQUS/Explicit到ABAQUS/CAE的过程日志 15.job_name.log 日志文件：包含了ABAQUS执行过程的起始时间等 16.job_name.abq ABAQUS/Explicit模块才有的状态文件，记录分析、继续和恢复命令。
为restart所需的文件。 17.job_name.mdl 模型文件：在ABAQUS/Standard和ABAQUS/Explict中运行数据检查后产生的文件，在analysis和continue指令下被读入并重写，为restart所需的文件。
18.job_name.pac 打包文件：包含了模型信息，仅用于ABAQUS/Explicit，该文件在执行analysis、datacheck命令时写入，执行analysis、continue、recover指令时读入，restart时需要的文件。 19.job_name.prt 零件信息文件：包含了零件与装配信息。
restart时需要 20.job_name.sel 结果选择文件：用于ABAQUS/Explicit，执行analysi、continue、recover指令时写入并有xonvert=select指令读入，为restart所需要的文件 21.job_name.stt 状态外文件：数据检查时写入的文件，在ABAQUS/Sandard中可在analysis、continue指令下读并写入，在ABAQUS/Explict中可在anslysi、continue指令下读入，为restart所需要的文件。 22.job_name.psf 脚本文件：用户定义parametric study时需要创建的文件 23.job_name.psr 参数化分析要求的输出结果，为文本格式 24.job_name.part 参数更改后重写的参数形式表示的inp文件 25.job_name.pes 参数更改后重写的inp文件 ABAQUS指令介绍 一般： *HEADING 定义分析的标题 节点定义： *NCOPY 使用平移、旋转、镜像的方法来产生新的节点集 *NFILL 在两组节点集中产生完整的节点。
节点距离可以是相等，或是成比例级数 *NGEN 在一条直线或曲线中产生节点集 *NODE 定义节点的坐标 *NSET 将某些节点集聚在一起并命名，之后在应用时便可直接使用这个节点集来定义性质 单元定义： *ELCOPY 产生新的单元 *ELEMENT 定义单元 *ELGEN 当以*ELEMENT定义完一个单元时，便可依此来产生新的单元 *ELSET 给予一单元或一单元集名称 元素性质定义 *RIGID SURFACE 在接触问题中定义刚性面 *BEAM SECTION 定义梁面元素 *SHELL SECTION 定义壳元素断面 *SOLID SECTION 定义固体元素 接触问题 *CONTACT PAIR 定义可能互相接触的一对面 *FRICTION 定义摩擦模型 材料性质 *MATERIAL 定义材料性质 *DAMPING 在动态问题中，用来定义阻尼系数 *DENSITY 在模态分析或瞬态分析时，定义材料比重 *ELASTIC 定义线弹性性质，对等向性材料与非等向性材料均可 *PLASTIC 使用Mises或Hill强屈服面来定义弹性材料，要先定义*ELASTIC *EXPANSION 定义热膨胀系数，可以是等向性与非等向性 约束条件 *BOUNDARY 用来描述某些节点固定位移（不能移动）与固定角度（不能转动） *EQUATION 用来约束多个点线性的关系（移动或转动） *MPC 用来约束多个点的约束 历程输入 *STEP 定义一分析步骤的起始 *END STEP 定义一分析步骤的结束 *INITIAL 用来定义分析的初始条件，可以是初始应力、应变、速度等 CONDITION *RESTART 用来控制分析结果（restart file *.res）的存取 *USER 使用子程序 SUBROUTINE 过程定义 *DYNAMIC 使用直接积分法来做动态响应分析 *FREQENCY 模态分析，计算固有频率及模态分析 *MODAL 使用模态叠加来做动态响应分析 DYNAMIC *STATIC 静态分析 *STEADY STATE 动态响应的 DYNAMICS 加载定义 力控制 *CLOAD 施加集中力或集中力矩OP=NEW：去除原本施力状态OP=MOD： 在原本施力状态下多加上其它力或修正的力（要加。
进入interaction模块
①指定裂纹special/creak/assign seam，选中示意图3中的黄色线，done！ ②生成裂纹crack 1,special/crack/create,name:crack 1,type: contour integral. 当提示选择裂纹前端时，选则示意图的红圈区域，当提示裂纹尖端区域时选择红圈的圆心，用向量q表示裂纹扩展方向（示意图3绿色箭头）。用同样的方法建立crack 2（示意图3中的蓝色区域）。
special/crack/edit，对两个裂纹进行应力奇异的设置