静-动力人工边界转换，粘弹性人工边界

Hi@田野

在地下结构静-动力分析中，以静力分析结果中的应力场为初始条件进行动力分析时，为保证模型在动力计算初始时刻为静力平衡状态，可在输入静应力场的同时输入引起此应力场的荷载和约束，转换边界时，需要去掉约束条件，同时施加响应的约束反力。计算步骤：（1）采用静力约束条件对模型进行静力分析，计算结果应包括应力场和约束条件（2）将约束条件去掉，施加粘弹性边界（3）在原静力约束被去掉的地方施加约束反力，对模型施加静力场和静力荷载（4）对模型施加动力荷载，进行动力计算
根据以上方法，我建立了简单的10*10*10m的模型
1.模型1，底部固定，其余侧面约束法向方向的位移，施加重力荷载，运行并输出应力场（S11-S23）。图1为运行后模型的应力场和位移云图。


2.模型2，底部固定，其余侧面约束法向方向的位移，施加重力荷载，修改关键词，加入应力场，运行并输出约束反力（RF1-RF3）。图2为运行后模型的应力场和位移云图。

3.模型3，底部和侧面的约束全部去掉，修改关键词，加入应力场和约束反力，施加重力荷载，运行。图3为运行后模型的应力场和位移云图。

在模型3，把固定边界去掉在相同位置增加了约束反力，但是为什么运行后的位移（-3）与模型1存在边界条件的位移（-18）相差这么大？

Hi@田野

dk贾队长 发表于 2020-5-15 10:28
额，抱歉，实在不好意思，这个论坛没有提示我有新的回复，我就没注意到你的回复信息，我的粘弹性边界是参 ...

地应力平衡不用这么麻烦，直接地应力平衡，把土留下其他结构杀死就可以平衡，加微信吧ty414468780

Hi@田野

静-动力边界转换：先在静力模型中设置两个分析步，静力通用（静力平衡，模型加重力）和动力隐式（动力计算），无边界约束，生成inp文件之后，在inp文件的相应位置把加粘弹性边界加入，导入模型，在荷载中，把等效节点力全部右移只存在step2中（动力分析步），然后生成inp，在inp文件中的**step之前把初始应力场加进去，在**loads之后通过*cload加入约束反力，再次把inp导入，运行就可以。

冒冒在进化

你好，可以加一个qq讨论一下静-动力边界这个转换的一些问题吗？我做了你和小木虫上差不多的算例。但是我移除固定边界+反力和应力场就去算地应力平衡，算不过。我的qq 862145668

dk贾队长

       楼主，我想请教一下，我的模型part很多，在对导出的应力文件修改的时候 修改节点编号前的part名太麻烦了。我能不能在模型计算前在命令行里面输入mdb.models['Model-1'].setValues(noPartsInputFile=ON) 这样的话就没有不同part的节点编号。
      因为看其他人说的输入那句py语句的时候，这样的话模型节点编号是从全局编写的，不存在节点编号重复的问题，就不需要标part名了？。
      我的想法是，1.在先把模型的粘弹性边界设置好，然后×掉，分析步先只设置静力平衡分析步，第一次在静力边界条件下地应力平衡后，导出应力文件，修改好应力文件.csv，打开第一次算的inp文件，在step前导入保存好，再在abaqus软件里面的job模块，计算新生成的inp文件，再导出RF文件，修改后，在新生成的inp再导入*Cload保存好。然后再CAE中设置一个静力平衡分析步，和一个动力分析步。把粘弹性边界√上。固定约束边界delete掉。在模型底部加上地震加速度时程。然后生成inp文件，再把修改好的应力文件和反力文件导入到inp文件里，保存好后进行计算。
老哥，帮我瞅瞅我计算思路有什么问题没？还有一个问题就是导出反力的时候 是导出全部节点的反力吗？因为我也有看到我土体里面的衬砌part也有反力

Hi@田野

dk贾队长 发表于 2020-4-28 16:24
楼主，我想请教一下，我的模型part很多，在对导出的应力文件修改的时候 修改节点编号前的part名太麻 ...

你的粘弹性边界是在inp文件中输入粘弹性边界的文件吗？我是在计算模型的inp文件中写几句话把粘弹性边界的inp文件写进去，然后在软件中导入Inp,导入粘弹性边界的同时还有地震波时程曲线转换的节点应力，特别多有几千个，如果x掉的话，要好长时间，所以我一般是没有粘弹性边界先进行静力，然后再加粘弹性以及地应力以及约束反力进行动力计算。我觉得修改part不麻烦啊，你下载一个破解版的ueditor，有列设置，一秒就能修改。到处节点的反力可以只是土体的也可以是全部的，你可以做个简单的模型验证一下。

dk贾队长

Hi@田野 发表于 2020-4-29 23:58
你的粘弹性边界是在inp文件中输入粘弹性边界的文件吗？我是在计算模型的inp文件中写几句话把粘弹性边界的 ...

额，抱歉，实在不好意思，这个论坛没有提示我有新的回复，我就没注意到你的回复信息，我的粘弹性边界是参考的技术邻冷月老师的用的CAE操作对模型表面的点设置set，然后用接地弹簧设置的弹簧刚度和阻尼刚度。这样的话确实麻烦一点，我不会你说的写入代码导入inp的方法。现在我是在运算第一次地应力平衡的模型之前输入了那行代码，计算出来的就只有一个PART了。然后在处理应力，反力文件的时候，不需要修改PART名，直接就是节点编号的。现在我这边地应力平衡还是有问题，我导入应力.csv文件计算不了， 我就采用的导入odb的方式进行地应力平衡。但是中途会死掉。另外我能不能加一下你的联系方式，想请教一下这个边界转化问题和等效节点荷载的问题。

shoupan

学习了:lol:lol:lol

永不放弃的慧

请问约束反力和应力具体是加到inp文件的那个位置啊？拜托看到回复一下