abaqus帮助文档翻译 2.1.11 一摞积木在通用接触下的倒塌分析
Abaqus基本操作中文教程
Abaqus基本操作中文教程目录1 Abaqus软件基本操作 (3)1.1 常用的快捷键 (3)1.2 单位的一致性 (3)1.3 分析流程九步走 (3)1.3.1 几何建模(Part) (4)1.3.2 属性设置(Property) (5)1.3.3 建立装配体(Assembly) (6)1.3.4 定义分析步(Step) (7)1.3.5 相互作用(Interaction) (8)1.3.6 载荷边界(Load) (10)1.3.7 划分网格(Mesh) (11)1.3.8 作业(Job) (15)1.3.9 可视化(Visualization) (16)1 Abaqus软件基本操作1.1 常用的快捷键旋转模型—Ctrl+Alt+鼠标左键平移模型—Ctrl+Alt+鼠标中键1.2 单位的一致性CAE软件其实是数值计算软件,没有单位的概念,常用的国际单位制如下表1所示,建议采用SI (mm)进行建模。
例如,模型的材料为钢材,采用国际单位制SI (m)时,弹性模量为2.06e11N/m2,重力加速度9.800 m/s2,密度为7850 kg/m3,应力Pa;采用国际单位制SI (mm)时,弹性模量为2.06e5N/mm2,重力加速度9800 mm/s2,密度为7850e-12 T/mm3,应力MPa。
1.3 分析流程九步走几何建模(Part)→属性设置(Property)→建立装配体(Assembly)→定义分析步(Step)→相互作用(Interaction)→载荷边界(Load)→划分网格(Mesh)→作业(Job)→可视化(Visualization)1.3.1 几何建模(Part ) 关键步骤的介绍: 部件(Part )导入Pro/E 等CAD 软件建好的模型后,另存成iges 、sat 、step 等格式;然后导入Abaqus 可以直接用,实体模型的导入通常采用sat 格式文件导入。
部件(Part )创建简单的部件建议直接在abaqus 中完成创建,复杂的可以借助Pro/E 或者Solidworks 等专业软件进行建模,然后导入。
ABAQUS帮助范例中文索引
帮助文档ABAQUS Example Problems Menual1.静态应力/位移分析1.1.静态与准静态应力分析1.1.1.螺栓结合型管法兰连接的轴对称分析1.1.2.薄壁机械肘在平面弯曲与内部压力下的弹塑性失效1.1.3.线弹性管线在平面弯曲下的参数研究1.1.4.橡胶海绵在圆形凸模下的变形分析1.1.5.混泥土板的失效1.1.6.有接缝的石坡稳定性研究1.1.7.锯齿状梁在循环载荷下的响应1.1.8.静水力学流体单元:空气弹簧模型1.1.9.管连接中的壳-固体子模型与壳-固体耦合的建立1.1.10.无应力单元的再激活1.1.11.黏弹性轴衬的动载响应1.1.12.厚板的凹入响应1.1.13.叠层复合板的损害和失效1.1.14.汽车密封套分析1.1.15.通风道接缝密封的压力渗透分析1.1.16.震动缓冲器的橡胶/海绵成分的自接触分析1.1.17.橡胶垫圈的橡胶/海绵成分的自接触分析1.1.18.堆叠金属片装配中的子模型分析1.1.19.螺纹连接的轴对称分析1.1.20.周期热-机械载荷下的汽缸盖的直接循环分析1.1.21.材料(沙产品)在油井中的侵蚀分析1.1.22.压力容器盖的子模型应力分析1.1.23.模拟游艇船体中复合涂覆层的应用1.2.屈曲与失效分析1.2.1.圆拱的完全弯曲分析1.2.2. 层压复合壳中带圆孔圆柱形面的屈曲分析1.2.3.点焊圆柱的屈曲分析1.2.4. K型结构的弹塑性分析1.2.5. 不稳定问题:压缩载荷下的加强板分析1.2.6.缺陷敏感柱型壳的屈曲分析1.3. 成形分析1.3.1. 圆柱形坯料墩粗:利用网格对网格方案配置与自适应网格的准静态分析1.3.2. 矩形方盒的超塑性成型1.3.3. 球形凸模的薄板拉伸1.3.4. 圆柱杯的深拉伸1.3.5. 考虑摩擦热产生的圆柱形棒材的挤压成形分析1.3.6. 厚板轧制成形分析1.3.7. 圆柱杯的轴对称成形分析1.3.8. 杯/槽成形分析1.3.9. 正弦曲线形凹模锻造1.3.10. 多重复合凹模锻造1.3.11. 平滑辊成形中的瞬态与稳态分析1.3.12. 型钢扎制成形分析1.3.13. 环扎成形分析1.3.14. 轴对称挤压成形中的瞬态与稳态分析1.3.15. 两步成形仿真1.3.16. 圆柱形坯料墩粗:热-位移耦合与隔热分析1.3.17. 金属板热成形中的不稳定静态问题分析1.4. 破裂与损伤1.4.1. 平板局部破裂分析:弹性线弹簧模拟1.4.2. 线弹性无限半空间中的锥裂纹围道积分1.4.3. 带局部轴向裂纹有限长度圆筒的弹塑性线弹簧模拟1.4.4. 三点弯曲试件的裂纹扩展1.4.5. 压力下刚性表面的松解工艺分析1.4.6. 钝角槽光纤金属绝缘板的失效分析1.5. 输入分析1.5.1. 二维拉伸弯曲的回弹1.5.2. 方形盒的深拉伸2. 动态应力/位移分析2.1. 动态应力分析2.1.1. 有局部非弹性失效结构的非线性动态分析2.1.2. Detroit Edison管滑轮试验2.1.3. 刚性抛丸对板的变化及影响2.1.4. 侵蚀抛丸对板的变化及影响2.1.5. 网球拍与球2.1.6. 变厚度燃料水槽壳的受压分析2.1.7. 汽车悬架模拟2.1.8. 管塞爆炸2.1.9. 常规接触的膝垫效应2.1.10. 常规接触的压接成形2.1.11. 常规接触的堆块失稳2.1.12. 带海绵效应限幅器的木桶坠落2.1.13. 铜杆的斜碰2.1.14. 带挡板木桶中的流体晃动2.1.15. 混泥土重力坝的地震分析2.1.16. 准静态或动载下薄壁铝挤压成形中的逐步损坏分析2.2. 基于模态的动态分析2.2.1.利用子结构和循环对称的旋转风扇分析2.2.2.Indian Point反应堆给水线的线性分析2.2.3. 三维框架建筑的感应波谱2.2.4. 使用平行Lanczos本征求解器结构的特征值分析2.2.5. 制动噪声分析2.2.6. 使用剩余模态的天线结构的动态分析2.2.7. 白车身模型的恒定动态分析2.3.联合仿真分析2.3.1. 充气门密封的关闭模拟。
ABAQUS帮助文档
初始损伤对应于材料开始退化,当应力或应变满足于定义的初始临界损伤准则,则此时退化开始。
Abaqus 的Damage for traction separation laws 中包括:Quade Damage、Maxe Damage、Quads Damage、Maxs Damage、Maxpe Damage、Maxps Damage 六种初始损伤准则,其中前四种用于一般复合材料分层模拟,后两种主要是在扩展有限元法模拟不连续体(比如crack 问题)问题时使用。
前四种对应于界面单元的含义如下:Maxe Damage 最大名义应变准则:Maxs Damage 最大名义应力准则:Quads Damage 二次名义应变准则:Quade Damage 二次名义应力准则最大主应力和最大主应变没有特定的联系,不同材料适用不同准则就像强度理论有最大应力理论和最大应变理论一样~ABAQUS帮助文档10.7.1 Modeling discontinuities as an enriched feature using the extended finite element method 看看里面有没有你想要的Defining damage evolution based on energy dissipated during the damage process根据损伤过程中消耗的能量定义损伤演变You can specify the fracture energy per unit area, , to be dissipated during the damage process directly.您可以指定每单位面积的断裂能量,在损坏过程中直接消散。
Instantaneous failure will occur if is specified as 0.瞬间失效将发生However, this choice is not recommended and should be used with care because it causes a sudden drop in the stress at the material point that can lead to dynamic instabilities.但是,不推荐这种选择,应谨慎使用,因为它会导致材料点的应力突然下降,从而导致动态不稳定。
(word完整版)Abaqus基本操作中文教程(2021年整理精品文档)
(word完整版)Abaqus基本操作中文教程(word完整版)Abaqus基本操作中文教程编辑整理:尊敬的读者朋友们:这里是精品文档编辑中心,本文档内容是由我和我的同事精心编辑整理后发布的,发布之前我们对文中内容进行仔细校对,但是难免会有疏漏的地方,但是任然希望((word完整版)Abaqus基本操作中文教程)的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈,这将是我们进步的源泉,前进的动力。
本文可编辑可修改,如果觉得对您有帮助请收藏以便随时查阅,最后祝您生活愉快业绩进步,以下为(word完整版)Abaqus基本操作中文教程的全部内容。
(word完整版)Abaqus基本操作中文教程Abaqus基本操作中文教程(word完整版)Abaqus基本操作中文教程目录1 Abaqus软件基本操作 (4)1.1 常用的快捷键 (4)1。
2 单位的一致性 (4)1。
3 分析流程九步走 (5)1。
3。
1 几何建模(Part) (5)1.3。
2 属性设置(Property) (7)1。
3。
3 建立装配体(Assembly) (7)1.3.4 定义分析步(Step) (9)1。
3.5 相互作用 (Interaction) (10)1。
3.6 载荷边界(Load) (13)1。
3。
7 划分网格 (Mesh) (14)1.3。
8 作业(Job) (18)1.3.9 可视化(Visualization) (19)(word完整版)Abaqus基本操作中文教程1 Abaqus软件基本操作1.1 常用的快捷键旋转模型— Ctrl+Alt+鼠标左键平移模型 - Ctrl+Alt+鼠标中键缩放模型 - Ctrl+Alt+鼠标右键1.2 单位的一致性CAE软件其实是数值计算软件,没有单位的概念,常用的国际单位制如下表1所示,建议采用SI (mm)进行建模。
国际单位SI (m)SI (mm)制长度m mm力N N质量kg t时间s sPa应力MPa (N/mm2)(N/m2)质量密度kg/m3t/mm3加速度m/s2mm/s2例如,模型的材料为钢材,采用国际单位制SI (m)时,弹性模量为2。
abaqus接触分析的常见问题
CAE(计算机辅助工程)是一门复杂的工程科学,涉及仿真技术、软件、产品设计和力学等众多领域。
世界上几大CAE公司各自以其独到的技术占领着相应的市场。
ABAQUS有限元分析软件拥有世界上最大的非线性力学用户群,是国际上公认的最先进的大型通用非线性有限元分析软件之一。
它广泛应用于机械制造、石油化工、航空航天、汽车交通、土木工程、国防军工、水利水电、生物医学、电子工程、能源、地矿、造船以及日用家电等工业和科学研究领域.ABAQUS在技术、品质和可靠性等方面具有卓越的声誉,可以对工程中各种复杂的线性和非线性问题进行分析计算。
《ABAQUS有限元分析常见问题解答》以问答的形式,详细介绍了使用ABAQUS建模分析过程中的各种常见问题,并以实例的形式教给读者如何分析问题、查找错误原因和尝试解决办法,帮助读者提高解决问题的能力。
《ABAQUS有限元分析常见问题解答》一书由机械工业出版社出版。
16.1。
1点对面离散与面对面离散【常见问题16-1】在ABAQUS/Standard分析中定义接触时,可以选择点对面离散方法(node-to-surf ace—dis—cre-tization)和面对面离散方法(surface-to—surfacediscretization),二者有何差别?『解答』在点对面离散方法中,从面(slavesurface)上的每个节点与该节点在主面(maste rsurface)上的投影点建立接触关系,每个接触条件都包含一个从面节点和它的投影点附近的一组主面节点.使用点对面离散方法时,从面节点不会穿透(penetrate)主面,但是主面节点可以穿透从面。
面对面离散方法会为整个从面(而不是单个节点)建立接触条件,在接触分析过程中同时考虑主面和从面的形状变化.可能在某些节点上出现穿透现象,但是穿透的程度不会很严重。
在如图16—l和图16—2所示的实例中,比较了两种情况。
1)从面网格比主面网格细:点对面离散(图16—1a)和面对面离散(图16-2a)的分析结果都很好,没有发生穿透,从面和主面都发生了正常的变形。
ABAQUS帮助文档翻译
节选-ABAQUS帮助文档翻译 reference to: user manual 18.62008-10-10 12:5918.6 理解自适应网格(adaptive meshing)自适应网格可以通过移动独立的材料网格(allowing the mesh to move independently of the material),让你在整个分析过程中即使发生大变形,也能保持高质量的网格。
通常自适应网格只移动节点,网格的拓扑并不改变。
注意:通常自适应网格多用在Dynamic (动态分析),Explicit and Dynamic(显示动态分析), Temp-disp, Explicit 中。
定义模型中某个区域采用自适应网格的设置:other-->Adaptive Mesh Domain 自适应网格的选项控制设置:Other--〉Adaptive Mesh Controls 通常,在每一个step中只能有一个自适应网格区域。
21.2.1 ABAQUS/Standard defines contact between two bodies in terms of two surfaces that may interact; these surfaces are called a “contact pair.”ABAQUS/Standard defines “self-contact,” which is available only in two-dimensional analysis, in terms of a single surface. [if gte vml 1]><![endif][if !vml][endif]Figure 21.2.1–1 Contact and interaction discretization. 从the first surface (the “slave” surface)的节点向the second surface (the “master” surface)做垂线,寻找最近的垂线的垂足,The interaction is then discretized between the point on the master surface and the slave node. Strict master-slave contact 在这种关系下,主面的节点可以穿入从面(副面),但副面不可以穿入主面。
abaqus使用手册[终稿]
Abaqus Example Problems Manual Introduction1 Static Stress/Displacement Analyses2 Dynamic S tress/Displacement Analyses3 Tire and V ehicle Analyses4 Mechanism Analyses5 Heat Transfer and Thermal-Stress Analyses6 Fluid Dynamics and Fluid-Structure Interaction7 Electrical Analyses8 Mass Diffusion Analyses9 Acoustic and Shock Analyses10 Soils Analyses11 Abaqus/Aqua Analyses12 Design Sensitivity Analyses13 Postprocessing of Abaqus Results Files Product Index ABAQUS例题手册介绍1静态应力/位移分析2动态应力/位移分析3轮胎和车辆分析4机理分析5传热和热应力分析6流体动力学和流体结构相互作用7电气分析8质量扩散分析9声和冲击分析10土壤分析11 ABAQUS /水族分析12设计灵敏度分析13后处理结果文件产品索引2 Dynamic Stress/Displacement Analyses2.1 Dynamic stress analyses2.1.1 Nonlinear dynamic analysis of a structure with local inelastic collapse2.1.2 Detroit Edison pipe whip experiment2.1.3 Rigid projectile impacting eroding plate2.1.4 Eroding projectile impacting eroding plate2.1.5 Tennis racket and ball2.1.6 Pressurized fuel tank with variable shell thickness 2.1.7 Modeling of an automobile suspension2.1.8 Explosive pipe closure2.1.9 Knee bolster impact with general contact2.1.10 Crimp forming with general contact2.1.11 Collapse of a stack of blocks with general contact 2.1.12 Cask drop with foam impact limiter2.1.13 Oblique impact of a copper rod2.1.14 Water sloshing in a baffled tank2.1.15 Seismic analysis of a concrete gravity dam2.1.16 Progressive failure analysis of thin-wall aluminum extrusion under quasi-static and dynamic loads2.1.17 Impact analysis of a pawl-ratchet device2.1.18 High-velocity impact of a ceramic target 2动态应力/位移分析2.1动态应力分析2.1.1非线性动力分析与当地的非弹性坍塌的结构2.1.2底特律爱迪生管鞭实验2.1.3刚性弹丸撞击侵蚀板2.1.4冲刷弹丸冲击侵蚀板2.1.5网球球拍和球2.1.6加压燃料箱具有可变壳厚度2.1.7建模的汽车悬架2.1.8爆管封2.1.9膝垫与一般的接触碰撞2.1.10压与一般的接触形成2.1.11折叠堆栈与通用接触块的2.1.12木桶降与泡沫冲击限制器2.1.13斜的影响铜棒的2.1.14水晃荡在挡板罐2.1.15抗震分析混凝土重力坝2.1.16渐进失效分析准静态和动态载荷作用下薄壁铝型材挤压2.1.17的影响分析一个棘爪棘轮装置2.1.18高速冲击陶瓷靶2.2 Mode-based dynamic analyses2.2.1 Analysis of a rotating fan using substructures and cyclic symmetry 2.2.2 Linear analysis of the Indian Point reactor feedwater line2.2.3 Response spectra of a three-dimensional frame building2.2.4 Brake squeal analysis2.2.5 Dynamic analysis of antenna structure utilizing residual modes 2.2.6 Steady-state dynamic analysis of a vehicle body-in-white model 2.3 Eulerian analyses2.3.1 Rivet forming2.3.2 Impact of a water-filled bottle using Eulerian-Lagrangian contact 2.4 Co-simulation analyses2.4.1 Closure of an air-filled door seal2.4.2 Dynamic impact of a scooter with a bump2.2模式为基础的动态分析2.2.1分析用子结构和循环对称旋转的风扇2.2.2线性分析印度点堆给水线2.2.3响应谱三维框架建设2.2.4制动尖叫分析2.2.5动态分析天线结构的利用残余模式2.2.6稳态动态分析汽车车身的白色模型2.3欧拉分析2.3.1铆钉形成2.3.2影响采用欧拉- 拉格朗日接触的充满水的瓶子2.4协同仿真分析2.4.1封闭的充气门封2.4.2动态影响与凸起的摩托车。
Abaqus帮助文档整理汇总
feature)组成,每一个部分至少有一个基
base feature),特征体可以是所创建的实体,如挤压体、
.首先建立“部件”
1)根据实际模型的尺寸决定部件的近似尺寸,进入绘图区。绘图
edit菜
sketcher options选项里调整。
(比如奇异)。 接触刚度的值决
当默认罚刚度设置用于罚函数
拉格朗日乘子默认不使用。如果用于罚函数
1000倍时,则默
-过
1000倍时,默认拉格朗日乘
:设置主面名2 v* c. b: S8 s) l
:设置允许违反接触条件的最大点数。这个条件由perrmx和
:使standard自动计算过盈容差和分离压力
以防止接触中的振荡。该参数不能与maxchp、perrmx和uerrmx
onset:设置其=immediate(默认)则在接触发生时在增量步
=delayed则延迟摩擦的应用。 G) P# q/ q7
:设置其=yes则强迫接触约束为拉格朗日乘子
=no则不使用拉格朗日乘子法。对于高刚度问题不推荐no,因为
3)分配截面特性给各特征体,把截面特性分配给部件的某一区域
.建立刚体
1)部件包括可变形体、不连续介质刚体和分析刚体三种类型,在
一旦建立后就不能更改其类型。采
在绘制轴对称部件的外形轮廓时不能超过其对
2)刚体是不能够施加质量、惯性轴等特性的,建立刚体后必须给
reference point)。在加载模块里对参考点施
solid element)只有平动自由度,没有转动自由度,所
ABAQUS将边界条件传递给其后的每一个分析步。对
ABAQUS命令参考
設定某些值的開啟與關閉 顯示某些值 設定觀看角度,亦可直接從滑鼠點選 增加,移除或修改視窗 放大或縮小視窗
*EL PRINT *ENERGY PRINT *MODAL PRINT *MONITOR *NODE PRINT *PRINT
定義哪些元素的應力,應變等變數要輸出 輸出彈性應變能,動能或塑性能等 輸出模態分析中的大小 觀察某點某一自由度,可用於初步判斷分析正確否,輸出 至 *.sta 輸出節點位移反力等 輸出 CONTACT:用於複雜接觸問題中,可用來觀察接觸或分 離 FREQUENCY:輸出的頻率
ABAQUS輸入檔指令簡錄 指令 一般 *HEADING 節點定義 *NCOPY *NFILL *NGEN *NODE *NSET 使用平移,旋轉,鏡射的方法來產生新的節點群 在兩組節點群中產生完整的節點。節點距離可以是相等, 或是成等比級數 在一條直線或是曲線中產生節點群 定義節點的座標 將某些節點群聚一起並給予命名,之後在應用時便可直接 使用這節點群來定義其性質 定義分析的標題 說明
定義線性彈性性質,對於等向性材料與非等向性材料均可 使用Miaes或是Hill降服曲面來定義彈塑性材料,要先定義 *ELASTIC 定義熱膨脹係數,可以是等向性與非等向性
用來描述某些節點固定位移(不能移動)與固定角度(不 能轉動) 用來約束多個點線性的關係(分析步驟的結束
輸出*.fil *EL FILE *ENERGY FILE *NODE FILE 輸出至.FIL中,可以在post中觀看 似*ENERGY PRINT 似*NODE PRINT
ABAQUS後處理指令整理 *ANIMATE *SET,BC DISPLAY= ON 用來產生動畫 在執行*DRAW時,顯示邊界條件
材料性質 *MATERIAL *DAMPING *DENSITY 定義材料性質 在動態問題中,用來定義阻尼係數 在模態分析或暫態分析時,定義材料比重
abaqus关键字翻译
abaqus关键字翻译接触相关的keyword:*Contact:开始定义通用的接触(该选项表明通用接触定义的开始。
每个step只能用一次,通用接触定义的不同方面可以通过下面的一些选项指定。
)产品:explicit可选参数:OP:设置OP=MOD(默认),更改已存的通用接触定义。
设置OP=NEW删除以前定义的接触并定义新的。
----------------------------*Contact Clearance:定义接触间隙属性(该选项用来创建接触间隙属性定义。
接触间隙属性将通过*Contact clearance assignment选项控制任何接触交互。
)产品:explicit必须参数:Name:定义属性名可选参数:Adjust:设置adjust=yes(默认),是通过调整节点坐标而无需创建约束来解决间隙问题。
adjust=yes只能用在第一个step定义间隙。
设置adjust=no则存储接触偏移以使间隙能被满足而不需调整节点坐标。
Clearance:设置该参数等于一个数值是为整个从节点集定义初始间隙或等于节点分布的名字。
对于实体单元表面上的从节点,间隙值必须是非负的。
默认是Search above:设置该参数等于表面上的距离加上指定的间隙值将作为搜索从节点的距离。
对于实体单元,默认距离是与从节点关联的单元尺寸的1/10。
对结构单元,默认是从节点相关的厚度。
Search below:设置该参数等于表面下的距离设置该参数等于表面上的距离加上指定的间隙值将作为搜索从节点的距离。
对于实体单元,默认距离是与从节点关联的单元尺寸的1/10。
对结构单元,默认是从节点相关的厚度。
------------------------*Contact clearance assignment:在一般接触区域的表面间施加接触间隙(该选项用来在接触面间定义初始接触间隙,并控制初始接触过盈如何解决。
)-------------*Contact controls:为接触指定额外的控制(该选项用来为接触模型提供额外的控制选项。
Abaqus User Subroutines Reference Guide 用户材料子程序帮助文档
1.1.41 UMATUser subroutine to define a material's mechanical behavior.Product: Abaqus/StandardWarning: The use of this subroutine generally requires considerable expertise. Y ou arecautioned that the implementation of any realistic constitutive model requires extensivedevelopment and testing. Initial testing on a single-element model with prescribedtraction loading is strongly recommended.References“User-defined mechanical material behavior,” Section 26.7.1 of the Abaqus Analysis User's Guide“User-defined thermal material behavior,” Section 26.7.2 of the Abaqus Analysis User's Guide*USER MA TERIAL“S D V I N I,” Section 4.1.11 of the Abaqus V erification Guide“U M A T and U H Y P E R,” Section 4.1.21 of the Abaqus V erification GuideOv erv iewUser subroutine U M A T:can be used to define the mechanical constitutive behavior of a material;will be called at all material calculation points of elements for which the material definition includes auser-defined material behavior;can be used with any procedure that includes mechanical behavior;can use solution-dependent state variables;must update the stresses and solution-dependent state variables to their values at the end of theincrement for which it is called;must provide the material Jacobian matrix, , for the mechanical constitutive model;can be used in conjunction with user subroutine U S D F L D to redefine any field variables before they are passed in; andis described further in “User-defined mechanical material behavior,” Section 26.7.1 of the AbaqusAnalysis User's Guide.Storage of stress and strain componentsIn the stress and strain arrays and in the matrices D D S D D E, D D S D D T, and D R P L D E, direct components are stored first, followed by shear components. There are N D I direct and N S H R engineering shear components. The order of the components is defined in “Conventions,” Section 1.2.2 of the Abaqus Analysis User's Guide. Since the number of active stress and strain components varies between element types, the routine must be coded toprovide for all element types with which it will be used.Defining local orientationsIf a local orientation (“Orientations,” Section 2.2.5 of the Abaqus Analysis User's Guide) is used at the same point as user subroutine U M A T, the stress and strain components will be in the local orientation; and, in the case of finite-strain analysis, the basis system in which stress and strain components are stored rotates with the material.StabilityY ou should ensure that the integration scheme coded in this routine is stable—no direct provision is made to include a stability limit in the time stepping scheme based on the calculations in U M A T.Convergence rateD D S D DE and—for coupled temperature-displacement and coupled thermal-electrical-structural analyses—D D S D D T, D R P L D E, and D R P L D T must be defined accurately if rapid convergence of the overall Newton scheme is to be achieved. In most cases the accuracy of this definition is the most important factor governing the convergence rate. Since nonsymmetric equation solution is as much as four times as expensive as the corresponding symmetric system, if the constitutive Jacobian (D D S D D E) is only slightly nonsymmetric (for example, a frictional material with a small friction angle), it may be less expensive computationally to use a symmetric approximation and accept a slower convergence rate.An incorrect definition of the material Jacobian affects only the convergence rate; the results (if obtained) are unaffected.Special considerations for various element typesThere are several special considerations that need to be noted.A v ailability of deformation gradientThe deformation gradient is available for solid (continuum) elements, membranes, and finite-strain shells(S3/S3R, S4, S4R, SAXs, and SAXAs). It is not available for beams or small-strain shells. It is stored as a 3× 3 matrix with component equivalence D F G R D0(I,J). For fully integrated first-order isoparametric elements (4-node quadrilaterals in two dimensions and 8-node hexahedra in three dimensions) the selectively reduced integration technique is used (also known as the technique). Thus, a modified deformation gradientis passed into user subroutine U M A T. For more details, see “Solid isoparametric quadrilaterals and hexahedra,”Section 3.2.4 of the Abaqus Theory Guide.Beams and shells that calculate transv erse shear energyIf user subroutine U M A T is used to describe the material of beams or shells that calculate transverse shear energy, you must specify the transverse shear stiffness as part of the beam or shell section definition to define the transverse shear behavior. See “Shell section behavior,” Section 29.6.4 of the Abaqus Analysis User's Guide, and “Choosing a beam element,” Section 29.3.3 of the Abaqus Analysis User's Guide, for informationon specifying this stiffness.Open-section beam elementsWhen user subroutine U M A T is used to describe the material response of beams with open sections (for example, an I-section), the torsional stiffness is obtained aswhere J is the torsional constant, A is the section area, k is a shear factor, and is the user-specified transverse shear stiffness (see “Transverse shear stiffness definition” in “Choosing a beam element,” Section29.3.3 of the Abaqus Analysis User's Guide).E lements w ith hourglassing modesIf this capability is used to describe the material of elements with hourglassing modes, you must define the hourglass stiffness factor for hourglass control based on the total stiffness approach as part of the element section definition. The hourglass stiffness factor is not required for enhanced hourglass control, but you can define a scaling factor for the stiffness associated with the drill degree of freedom (rotation about the surface normal). See “Section controls,” Section 27.1.4 of the Abaqus Analysis User's Guide, for information on specifying the stiffness factor.Pipe-soil interaction elementsThe constitutive behavior of the pipe-soil interaction elements (see “Pipe-soil interaction elements,” Section 32.12.1 of the Abaqus Analysis User's Guide) is defined by the force per unit length caused by relative displacement between two edges of the element. The relative-displacements are available as “strains” (S T R A N and D S T R A N). The corresponding forces per unit length must be defined in the S T R E S S array. The Jacobian matrix defines the variation of force per unit length with respect to relative displacement.For two-dimensional elements two in-plane components of “stress” and “strain” exist (N T E N S=N D I=2, andN S H R=0). For three-dimensional elements three components of “stress” and “strain” exist (N T E N S=N D I=3, and N S H R=0).Large volume changes with geometric nonlinearityIf the material model allows large volume changes and geometric nonlinearity is considered, the exact definition of the consistent Jacobian should be used to ensure rapid convergence. These conditions are most commonly encountered when considering either large elastic strains or pressure-dependent plasticity. In the former case, total-form constitutive equations relating the Cauchy stress to the deformation gradient are commonly used; in the latter case, rate-form constitutive laws are generally used.For total-form constitutive laws, the exact consistent Jacobian is defined through the variation in Kirchhoff stress:Here, J is the determinant of the deformation gradient, is the Cauchy stress, is the virtual rate of deformation, and is the virtual spin tensor, defined asFor rate-form constitutive laws, the exact consistent Jacobian is given byUse with incompressible elastic materialsFor user-defined incompressible elastic materials, user subroutine U H Y P E R should be used rather than user subroutine U M A T. In U M A T incompressible materials must be modeled via a penalty method; that is, you must ensure that a finite bulk modulus is used. The bulk modulus should be large enough to model incompressibility sufficiently but small enough to avoid loss of precision. As a general guideline, the bulk modulus should be about – times the shear modulus. The tangent bulk modulus can be calculated fromIf a hybrid element is used with user subroutine U M A T, Abaqus/Standard will replace the pressure stress calculated from your definition of S T R E S S with that derived from the Lagrange multiplier and will modify the Jacobian appropriately.For incompressible pressure-sensitive materials the element choice is particularly important when using user subroutine U M A T. In particular, first-order wedge elements should be avoided. For these elements the technique is not used to alter the deformation gradient that is passed into user subroutine U M A T, which increases the risk of volumetric locking.Increments for which only the Jacobian can be definedAbaqus/Standard passes zero strain increments into user subroutine U M A T to start the first increment of all the steps and all increments of steps for which you have suppressed extrapolation (see “Defining an analysis,”Section 6.1.2 of the Abaqus Analysis User's Guide). In this case you can define only the Jacobian (D D S D D E).Utility routinesSeveral utility routines may help in coding user subroutine U M A T. Their functions include determining stress invariants for a stress tensor and calculating principal values and directions for stress or strain tensors. These utility routines are discussed in detail in “Obtaining stress invariants, principal stress/strain values and directions, and rotating tensors in an Abaqus/Standard analysis,” Section 2.1.11.U ser subroutine interfaceS U B R O U T I N E U M A T(S T R E S S,S T A T E V,D D S D D E,S S E,S P D,S C D,1R P L,D D S D D T,D R P L D E,D R P L D T,2S T R A N,D S T R A N,T I M E,D T I M E,T E M P,D T E M P,P R E D E F,D P R E D,C M N A M E,3N D I,N S H R,N T E N S,N S T A T V,P R O P S,N P R O P S,C O O R D S,D R O T,P N E W D T,4C E L E N T,D F G R D0,D F G R D1,N O E L,N P T,L A Y E R,K S P T,K S T E P,K I N C)CI N C L U D E'A B A_P A R A M.I N C'C H A R A C T E R*80C M N A M ED I ME N S I O N S T R E S S(N T E N S),S T A T E V(N S T A T V),1D D S D D E(N T E N S,N T E N S),D D S D D T(N T E N S),D R P L D E(N T E N S),2S T R A N(N T E N S),D S T R A N(N T E N S),T I M E(2),P R E D E F(1),D P R E D(1),3P R O P S(N P R O P S),C O O R D S(3),D R O T(3,3),D F G R D0(3,3),D F G R D1(3,3)user coding to define D D S D D E,S T R E S S,S T A T E V,S S E,S P D,S C Dand, if necessary,R P L,D D S D D T,D R P L D E,D R P L D T,P N E W D TR E T U R NE N DV ariables to be definedIn all situationsD D S D D E(N TE N S,N T E N S)Jacobian matrix of the constitutive model, , where are the stress increments and are the strain increments. D D S D D E(I,J) defines the change in the I th stress component at the end of the time increment caused by an infinitesimal perturbation of the J th component of the strain increment array.Unless you invoke the unsymmetric equation solution capability for the user-defined material,Abaqus/Standard will use only the symmetric part of D D S D D E. The symmetric part of the matrix iscalculated by taking one half the sum of the matrix and its transpose.S T R E S S(N T E N S)This array is passed in as the stress tensor at the beginning of the increment and must be updated in this routine to be the stress tensor at the end of the increment. If you specified initial stresses (“Initial conditions in Abaqus/Standard and Abaqus/Explicit,” Section 34.2.1 of the Abaqus Analysis User's Guide), this array will contain the initial stresses at the start of the analysis. The size of this array depends on the value of N T E N S as defined below. In finite-strain problems the stress tensor has already been rotated to account for rigid body motion in the increment before U M A T is called, so that only the corotational part of the stress integration should be done in U M A T. The measure of stress used is “true” (Cauchy) stress.S T A T E V(N S T A T V)An array containing the solution-dependent state variables. These are passed in as the values at thebeginning of the increment unless they are updated in user subroutines U S D F L D or U E X P A N, in which case the updated values are passed in. In all cases S T A T E V must be returned as the values at the end of the increment. The size of the array is defined as described in “Allocating space” in “User subroutines:overview,” Section 18.1.1 of the Abaqus Analysis User's Guide.In finite-strain problems any vector-valued or tensor-valued state variables must be rotated to account for rigid body motion of the material, in addition to any update in the values associated with constitutivebehavior. The rotation increment matrix, D R O T, is provided for this purpose.S S E,S P D,S C DSpecific elastic strain energy, plastic dissipation, and “creep” dissipation, respectively. These are passed in as the values at the start of the increment and should be updated to the corresponding specific energy values at the end of the increment. They have no effect on the solution, except that they are used forenergy output.Only in a fully coupled thermal-stress or a coupled thermal-electrical-structural analysisR P LV olumetric heat generation per unit time at the end of the increment caused by mechanical working of the material.D D S D D T(N TE N S)V ariation of the stress increments with respect to the temperature.D R P L D E(N TE N S)V ariation of R P L with respect to the strain increments.D R P L D TV ariation of R P L with respect to the temperature.Only in a geostatic stress procedure or a coupled pore fluid diffusion/stress analysis for pore pressure cohesive elementsR P LR P L is used to indicate whether or not a cohesive element is open to the tangential flow of pore fluid. Set R P L equal to 0 if there is no tangential flow; otherwise, assign a nonzero value to R P L if an element is open.Once opened, a cohesive element will remain open to the fluid flow.V ariable that can be updatedP N E W D TRatio of suggested new time increment to the time increment being used (D T I M E, see discussion later in this section). This variable allows you to provide input to the automatic time incrementation algorithms in Abaqus/Standard (if automatic time incrementation is chosen). For a quasi-static procedure the automatic time stepping that Abaqus/Standard uses, which is based on techniques for integrating standard creep laws (see “Quasi-static analysis,” Section 6.2.5 of the Abaqus Analysis User's Guide), cannot becontrolled from within the U M A T subroutine.P N E W D T is set to a large value before each call to U M A T.If P N E W D T is redefined to be less than 1.0, Abaqus/Standard must abandon the time increment andattempt it again with a smaller time increment. The suggested new time increment provided to theautomatic time integration algorithms is P N E W D T × D T I M E, where the P N E W D T used is the minimum value for all calls to user subroutines that allow redefinition of P N E W D T for this iteration.If P N E W D T is given a value that is greater than 1.0 for all calls to user subroutines for this iteration and the increment converges in this iteration, Abaqus/Standard may increase the time increment. The suggestednew time increment provided to the automatic time integration algorithms is P N E W D T × D T I M E, where the P N E W D T used is the minimum value for all calls to user subroutines for this iteration.If automatic time incrementation is not selected in the analysis procedure, values of P N E W D T that aregreater than 1.0 will be ignored and values of P N E W D T that are less than 1.0 will cause the job to terminate. V ariables passed in for informationS T R A N(N T E N S)An array containing the total strains at the beginning of the increment. If thermal expansion is included in the same material definition, the strains passed into U M A T are the mechanical strains only (that is, thethermal strains computed based upon the thermal expansion coefficient have been subtracted from the total strains). These strains are available for output as the “elastic” strains.In finite-strain problems the strain components have been rotated to account for rigid body motion in the increment before U M A T is called and are approximations to logarithmic strain.D S T R A N(N TE N S)Array of strain increments. If thermal expansion is included in the same material definition, these are the mechanical strain increments (the total strain increments minus the thermal strain increments).T I M E(1)V alue of step time at the beginning of the current increment or frequency.T I M E(2)V alue of total time at the beginning of the current increment.D T I M ETime increment.T E M PTemperature at the start of the increment.D TE M PIncrement of temperature.P R E D E FArray of interpolated values of predefined field variables at this point at the start of the increment, based on the values read in at the nodes.D P RE DArray of increments of predefined field variables.C M N A M EUser-defined material name, left justified. Some internal material models are given names starting with the “ABQ_” character string. To avoid conflict, you should not use “ABQ_” as the leading string for C M N A M E. N D INumber of direct stress components at this point.N S H RNumber of engineering shear stress components at this point.N T E N SSize of the stress or strain component array (N D I + N S H R).N S T A T VNumber of solution-dependent state variables that are associated with this material type (defined asdescribed in “Allocating space” in “User subroutines: overview,” Section 18.1.1 of the Abaqus Analysis User's Guide).P R O P S(N P R O P S)User-specified array of material constants associated with this user material.N P R O P SUser-defined number of material constants associated with this user material.C O O RD SAn array containing the coordinates of this point. These are the current coordinates if geometricnonlinearity is accounted for during the step (see “Defining an analysis,” Section 6.1.2 of the Abaqus Analysis User's Guide); otherwise, the array contains the original coordinates of the point.D R O T(3,3)Rotation increment matrix. This matrix represents the increment of rigid body rotation of the basis system in which the components of stress (S T R E S S) and strain (S T R A N) are stored. It is provided so that vector-or tensor-valued state variables can be rotated appropriately in this subroutine: stress and straincomponents are already rotated by this amount before U M A T is called. This matrix is passed in as a unit matrix for small-displacement analysis and for large-displacement analysis if the basis system for thematerial point rotates with the material (as in a shell element or when a local orientation is used).C E L E N TCharacteristic element length, which is a typical length of a line across an element for a first-order element;it is half of the same typical length for a second-order element. For beams and trusses it is a characteristic length along the element axis. For membranes and shells it is a characteristic length in the referencesurface. For axisymmetric elements it is a characteristic length in the plane only. For cohesiveelements it is equal to the constitutive thickness.D F G R D0(3,3)Array containing the deformation gradient at the beginning of the increment. If a local orientation is defined at the material point, the deformation gradient components are expressed in the local coordinate system defined by the orientation at the beginning of the increment. For a discussion regarding the availability of the deformation gradient for various element types, see “Availability of deformation gradient.”D F G R D1(3,3)Array containing the deformation gradient at the end of the increment. If a local orientation is defined at the material point, the deformation gradient components are expressed in the local coordinate system defined by the orientation. This array is set to the identity matrix if nonlinear geometric effects are not included in the step definition associated with this increment. For a discussion regarding the availability of thedeformation gradient for various element types, see “Availability of deformation gradient.”N O E LElement number.N P TIntegration point number.L A Y E RLayer number (for composite shells and layered solids).K S P TSection point number within the current layer.K S T E PStep number.K I N CIncrement number.Example: Using more than one user-defined mechanical material modelTo use more than one user-defined mechanical material model, the variable C M N A M E can be tested for different material names inside user subroutine U M A T as illustrated below:I F(C M N A M E(1:4).E Q.'M A T1')T H E NC A L L U M A T_M A T1(argument_list)E L S E I F(C M N A M E(1:4).E Q.'M A T2')T H E NC A L L U M A T_M A T2(argument_list)E N D I FU M A T_M A T1 and U M A T_M A T2 are the actual user material subroutines containing the constitutive material models for each material M A T1 and M A T2, respectively. Subroutine U M A T merely acts as a directory here. The argument list may be the same as that used in subroutine U M A T.Example: Simple linear viscoelastic materialAs a simple example of the coding of user subroutine U M A T, consider the linear, viscoelastic model shown in Figure 1.1.41–1. Although this is not a very useful model for real materials, it serves to illustrate how to code the routine.Figure 1.1.41–1 Simple linear viscoelastic model.The behavior of the one-dimensional model shown in the figure iswhere and are the time rates of change of stress and strain. This can be generalized for small straining of an isotropic solid asandwhereand , , , , and are material constants ( and are the Lamé constants).A simple, stable integration operator for this equation is the central difference operator:where f is some function, is its value at the beginning of the increment, is the change in the function over the increment, and is the time increment.Applying this to the rate constitutive equations above givesandso that the Jacobian matrix has the termsandThe total change in specific energy in an increment for this material iswhile the change in specific elastic strain energy iswhere D is the elasticity matrix:No state variables are needed for this material, so the allocation of space for them is not necessary. In a morerealistic case a set of parallel models of this type might be used, and the stress components in each model might be stored as state variables.For our simple case a user material definition can be used to read in the five constants in the order , , , , and so thatThe routine can then be coded as follows:S U B R O U T I N E U M A T(S T R E S S,S T A T E V,D D S D D E,S S E,S P D,S C D,1R P L,D D S D D T,D R P L D E,D R P L D T,2S T R A N,D S T R A N,T I M E,D T I M E,T E M P,D T E M P,P R E D E F,D P R E D,C M N A M E,3N D I,N S H R,N T E N S,N S T A T V,P R O P S,N P R O P S,C O O R D S,D R O T,P N E W D T,4C E L E N T,D F G R D0,D F G R D1,N O E L,N P T,L A Y E R,K S P T,K S T E P,K I N C)CI N C L U D E'A B A_P A R A M.I N C'CC H A R A C T E R*80C M N A M ED I ME N S I O N S T R E S S(N T E N S),S T A T E V(N S T A T V),1D D S D D E(N T E N S,N T E N S),2D D S D D T(N T E N S),D R P L D E(N T E N S),3S T R A N(N T E N S),D S T R A N(N T E N S),T I M E(2),P R E D E F(1),D P R E D(1),4P R O P S(N P R O P S),C O O R D S(3),D R O T(3,3),D F G R D0(3,3),D F G R D1(3,3)D I ME N S I O N D S T R E S(6),D(3,3)CC E V A L U A T E N E W S T R E S S T E N S O RCE V=0.D E V=0.D O K1=1,N D IE V=E V+S T R A N(K1)D E V=D E V+D S T R A N(K1)E N D D OCT E R M1=.5*D T I M E+P R O P S(5)T E R M1I=1./T E R M1T E R M2=(.5*D T I M E*P R O P S(1)+P R O P S(3))*T E R M1I*D E VT E R M3=(D T I M E*P R O P S(2)+2.*P R O P S(4))*T E R M1ICD O K1=1,N D ID S T RE S(K1)=T E R M2+T E R M3*D S T R A N(K1)1+D T I M E*T E R M1I*(P R O P S(1)*E V2+2.*P R O P S(2)*S T R A N(K1)-S T R E S S(K1))S T R E S S(K1)=S T R E S S(K1)+D S T R E S(K1)E N D D OCT E R M2=(.5*D T I M E*P R O P S(2)+P R O P S(4))*T E R M1II1=N D ID O K1=1,N S H RI1=I1+1D S T RE S(I1)=T E R M2*D S T R A N(I1)+1D T I M E*T E R M1I*(P R O P S(2)*S T R A N(I1)-S T R E S S(I1)) S T R E S S(I1)=S T R E S S(I1)+D S T R E S(I1)E N D D OCC C R E A T E N E W J A C O B I A NCT E R M2=(D T I M E*(.5*P R O P S(1)+P R O P S(2))+P R O P S(3)+12.*P R O P S(4))*T E R M1IT E R M3=(.5*D T I M E*P R O P S(1)+P R O P S(3))*T E R M1ID O K1=1,N TE N SD O K2=1,N TE N SD D S D D E(K2,K1)=0.E N D D OE N D D OCD O K1=1,N D ID D S D D E(K1,K1)=TE R M2E N D D OCD O K1=2,N D IN2=K1–1D O K2=1,N2D D S D D E(K2,K1)=TE R M3D D S D D E(K1,K2)=TE R M3E N D D OE N D D OT E R M2=(.5*D T I M E*P R O P S(2)+P R O P S(4))*T E R M1II1=N D ID O K1=1,N S H RI1=I1+1D D S D D E(I1,I1)=TE R M2E N D D OCC T O T A L C H A N G E I N S P E C I F I C E N E R G YCT D E=0.D O K1=1,N TE N ST D E=T D E+(S T R E S S(K1)-.5*D S T R E S(K1))*D S T R A N(K1)E N D D OCC C H A N G E I N S P E C I F I C E L A S T I C S T R A I N E N E R G YCT E R M1=P R O P S(1)+2.*P R O P S(2)D O K1=1,N D ID(K1,K1)=T E R M1E N D D OD O K1=2,N D IN2=K1-1D O K2=1,N2D(K1,K2)=P R O P S(1)D(K2,K1)=P R O P S(1)E N D D OE N D D OD E E=0.D O K1=1,N D IT E R M1=0.T E R M2=0.D O K2=1,N D IT E R M1=T E R M1+D(K1,K2)*S T R A N(K2)T E R M2=T E R M2+D(K1,K2)*D S T R A N(K2)E N D D OD E E=D E E+(T E R M1+.5*T E R M2)*D S T R A N(K1)E N D D OI1=N D ID O K1=1,N S H RI1=I1+1D E E=D E E+P R O P S(2)*(S T R A N(I1)+.5*D S T R A N(I1))*D S T R A N(I1)E N D D OS S E=S S E+D E ES C D=S C D+T D E– D E ER E T U R NE N D。
ABAQUS简易培训教材(中文)
指定为静态分析过程 载荷定义,11:节点号,2:自由度 -1200.0:载荷大小
输出数据
历程数据以*end step 选项结束
在输入文件中使用集名引用属性:
*ELEMENT, TYPE=B21, ELSET=BEAMS 1, 1, 3 *BEAM SECTION, SECTION=RECT, ELSET=BEAMS, MATERIAL=MAT1 50.0, 5.0 *MATERIAL, NAME=MAT1 *ELASTIC 2.0E5, 0.3
*BEAM SECTION 为单元集 BEAMS 和材料集 MAT1 建立联系。在*BEAM SECTION 选项中,横截面为长 方形(RECT),宽度为 50.0,高度为 5.0。
在*MATERIAL 选项块中,材料名为 MAT1,弹性模量为 2.0E5,泊松比为 0.3。
边界条件:
*BOUNDARY
564, 102, 103
数据行
572, 103, 104
·
节点号(相对于梁 b21 单元)
·
单元号
每个数据块要么属于模型数据,要么属于历程数据,模型数据必然置于历程数据之前。而在模型数据和历程数
据内部,数据块的顺序和位置是任意的,除了一些特例,如:*HEADING 必须置于输入文件的第一行,*ELASTIC、 *DENSITY 和*PLASTIC 是*MATERIAL 的子选项,则他们必须直接跟在*MATERIAL 后等 典型例题
我们将通过 ABAQUS/CAE 完成右图的建模及分析过程。
首先我们创建几何体 一、创建基本特征:
1、首先运行 ABAQUS/CAE,在出现的对话框内 选择 Create Model Database。
2、从 Module 列表中选择 Part,进入 Part 模块
ABAQUS使用手册(中文版)
ABAQUS使用手册(中文版)abaqus入门使用手册abaqus简介:abaqus是一套先进的通用有限元程序系统,这套软件的目的是对固体和结构的力学问题进行数值计算分析,而我们将其用于材料的计算机模拟及其前后处理,主要得益于abaqus给我们的abaqus/standard及abaqus/explicit通用分析模块。
abaqus存有众多的分析模块,我们采用的模块主要就是abaqus/cae及viewer,前者用作建模及适当的前处置,后者用作对结果展开分析及处置。
下面将对这两个模块的采用融合本人的体会搞一些具体内容的表明:一.abaqus/caecae模块用于分析对象的建模,特性及约束条件的给定,网格的划分以及数据传输等等,其核心由七个步骤组成,下面将对这七个步骤作出说明:1.part步(1)part→creatmodelingspace:①3d代表三维②2d代表二维③aa xisymmetric代表轴对称,这三个选项的选取视乎所演示对象的结构而的定。
type:①deformable为一般选项,适合于绝大多数的模拟对象。
②discreterigid和analyticalrigid用于多个物体组合时,与我们所研究的对象相关的物体上。
abaqus假设这些与所研究的对象相关的物体均为刚体,对于其中较简单的刚体,如球体而言,选择前者即可。
若刚体形状较复杂,或者不是规则的几何图形,那么就选择后者。
需要说明的是,由于后者所建立的模型是离散的,所以只能是近似的,不可能和实际物体一样,因此误差较大。
shape中存有四个选项,其排序规则就是按照维数而的定的,可以根据我们的演示对象确认。
1type:①extrusion用作创建通常情况的三维模型②revolution创建旋转体模型③sweep用作创建形状任一的模型。
approximatesize:在此栏中设定作图区的大致尺寸,其单位与我们选定的单位一致。
设置完,页面continue步入作图区。
ABAQUS帮助里关键字(keywords)翻译
ABAQUS帮助里关键字(keywords)翻译ABAQUS帮助里关键字(keywords)翻译总规则1、关键字必须以*号开头,且关键字前无空格2、**为注释行,它可以出现在中的任何地方3、当关键字后带有时,关键词后必须采用逗号隔开4、参数间都采用逗号隔开5、关键词可以采用简写的方式,只要程序能识别就可以了6、不需使用隔行符,如果参数比较多,一行放不下,可以另起一行,只要在上一行的末尾加逗号便可以----------------------------------------------------------------------------------------------------------------------------------------- *AMPLITUDE:幅值这个选项允许任意的载荷、和其它指定的数值在一个分析步中随时间的变化(或者在ABAQUS/Standard分析中随着的变化)。
必需的参数:NAME:幅值曲线的名字可选参数:DEFINITION:设置definition=Tabular(默认)给出表格形式的幅值-时间(或幅值-频率)定义。
设置DEFINITION=EQUALLY SPACED/PERIODIC/MODULATED/DECAY/SMOOTHSTEP/SOLUTION DEPENDENT或BUBBLE来定义其他形式的幅值曲线。
INPUT:设置该参数等于替换输入文件名字。
TIME:设置TIME=STEP TIME(默认)则表示分析步时间或频率。
TIME=TOTAL TIME表示总时间。
VALUE:设置VALUE=RELATIVE(默认),定义相对幅值。
VALUE=ABSOLUTE表示绝对幅值,此时,行中载荷选项内的值将被省略,而且当温度是指定给已定义了温度TEMPERATURE=GRADIENTS(默认)梁上或壳上的,不能使用ABSOLUTE。
对于DEFINITION=TABULAR的可选参数:SMOOTH:设置该参数等于DEFINITION=TABULAR的数据行第一行1、时间或频率2、第一点的幅值(绝对或相对)3、时间或频率4、第二点的幅值(绝对或相对) 等等基本形式:*Amplitude,name=Amp-10.,0.,0.2,1.5,0.4,2.,1.,1.*BEAM SECTION:当需要数值积分时定义梁截面*BOND:定义绑定和绑定属性*BOUNDARY:定义边界条件用来在节点定义边界条件或在子模型分析中指定被驱动的节点。
Abaqus帮助文档整理汇总
Abaqus 使用日记Abaqus标准版共有“部件(part)”、“材料特性(propoterty)”、“装配(assemble)”、“计算步骤(step)”、“交互(interaction)”、“加载(load)”、“单元划分(mesh)”、“计算(job)”、“后处理(visualization)”、“草图(sketch)”十大模块组成。
建模方法:一个模型(model)通常由一个或几个部件(part)组成,“部件”又由一个或几个特征体(feature)组成,每一个部分至少有一个基本特征体(base feature),特征体可以是所创建的实体,如挤压体、切割挤压体、数据点、参考点、数据轴,数据平面,装配体的装配约束、装配体的实例等等。
1.首先建立“部件”(1)根据实际模型的尺寸决定部件的近似尺寸,进入绘图区。
绘图区根据所输入的近似尺寸决定网格的间距,间距大小可以在edit菜单sketcher options选项里调整。
(2)在绘图区分别建立部件中的各个特征体,建立特征体的方法主要有挤压、旋转、平扫三种。
同一个模型中两个不同的部件可以有同名的特征体组成,也就是说不同部件中可以有同名的特征体,同名特征体可以相同也可以不同。
部件的特征体包括用各种方法建立的基本特征体、数据点(datum point)、数据轴(datum axis)、数据平面(datum plane)等等。
(3)编辑部件可以用部件管理器进行部件复制,重命名,删除等,部件中的特征体可以是直接建立的特征体,还可以间接手段建立,如首先建立一个数据点特征体,通过数据点建立数据轴特征体,然后建立数据平面特征体,再由此基础上建立某一特征体,最先建立的数据点特征体就是父特征体,依次往下分别为子特征体,删除或隐藏父特征体其下级所有子特征体都将被删除或隐藏。
××××特征体被删除后将不能够恢复,一个部件如果只包含一个特征体,删除特征体时部件也同时被删除×××××2.建立材料特性(1)输入材料特性参数弹性模量、泊松比等(2)建立截面(section)特性,如均质的、各项同性、平面应力平面应变等等,截面特性管理器依赖于材料参数管理器(3)分配截面特性给各特征体,把截面特性分配给部件的某一区域就表示该区域已经和该截面特性相关联3.建立刚体(1)部件包括可变形体、不连续介质刚体和分析刚体三种类型,在创建部件时需要指定部件的类型,一旦建立后就不能更改其类型。
abaqus接触分析的常见问题【范本模板】
CAE(计算机辅助工程)是一门复杂的工程科学,涉及仿真技术、软件、产品设计和力学等众多领域。
世界上几大CAE公司各自以其独到的技术占领着相应的市场。
ABAQUS有限元分析软件拥有世界上最大的非线性力学用户群,是国际上公认的最先进的大型通用非线性有限元分析软件之一。
它广泛应用于机械制造、石油化工、航空航天、汽车交通、土木工程、国防军工、水利水电、生物医学、电子工程、能源、地矿、造船以及日用家电等工业和科学研究领域。
ABAQUS在技术、品质和可靠性等方面具有卓越的声誉,可以对工程中各种复杂的线性和非线性问题进行分析计算.《ABAQUS有限元分析常见问题解答》以问答的形式,详细介绍了使用ABAQUS建模分析过程中的各种常见问题,并以实例的形式教给读者如何分析问题、查找错误原因和尝试解决办法,帮助读者提高解决问题的能力。
《ABAQUS有限元分析常见问题解答》一书由机械工业出版社出版。
16。
1.1点对面离散与面对面离散【常见问题16-1】在ABAQUS/Standard分析中定义接触时,可以选择点对面离散方法(node-to-sur face-dis—cre-tization)和面对面离散方法(surface-to-surfacediscretization),二者有何差别?『解答』在点对面离散方法中,从面(slavesurface)上的每个节点与该节点在主面(mastersu rface)上的投影点建立接触关系,每个接触条件都包含一个从面节点和它的投影点附近的一组主面节点。
使用点对面离散方法时,从面节点不会穿透(penetrate)主面,但是主面节点可以穿透从面。
面对面离散方法会为整个从面(而不是单个节点)建立接触条件,在接触分析过程中同时考虑主面和从面的形状变化.可能在某些节点上出现穿透现象,但是穿透的程度不会很严重。
在如图16-l和图16-2所示的实例中,比较了两种情况。
1)从面网格比主面网格细:点对面离散(图16—1a)和面对面离散(图16-2a)的分析结果都很好,没有发生穿透,从面和主面都发生了正常的变形。
ABAQUS帮助文档翻译
节选-ABAQUS帮助文档翻译 reference to: user manual 18.62008-10-10 12:5918.6 理解自适应网格(adaptive meshing)自适应网格可以通过移动独立的材料网格(allowing the mesh to move independently of the material),让你在整个分析过程中即使发生大变形,也能保持高质量的网格。
通常自适应网格只移动节点,网格的拓扑并不改变。
注意:通常自适应网格多用在Dynamic (动态分析),Explicit and Dynamic(显示动态分析), Temp-disp, Explicit 中。
定义模型中某个区域采用自适应网格的设置:other-->Adaptive Mesh Domain 自适应网格的选项控制设置:Other--〉Adaptive Mesh Controls 通常,在每一个step中只能有一个自适应网格区域。
21.2.1 ABAQUS/Standard defines contact between two bodies in terms of two surfaces that may interact; these surfaces are called a “contact pair.”ABAQUS/Standard defines “self-contact,” which is available only in two-dimensional analysis, in terms of a single surface. [if gte vml 1]><![endif][if !vml][endif]Figure 21.2.1–1 Contact and interaction discretization. 从the first surface (the “slave” surface)的节点向the second surface (the “master” surface)做垂线,寻找最近的垂线的垂足,The interaction is then discretized between the point on the master surface and the slave node. Strict master-slave contact 在这种关系下,主面的节点可以穿入从面(副面),但副面不可以穿入主面。
ABAQUS帮助-用户分析手册目录---在帮助文件中查找相关信息
ABAQUS Analysis User’s Manual目录第1章介绍1.1 介绍1.1.1 介绍:概要1.2 ABAQUS构造和约定1.2.1 Input构造规则1.2.2 约定1.3 定义一个ABAQUS模型1.3.1 在ABAQUS中定义一个模型1.4 参数模型1.4.1 参数输入第2章空间模型2.1 定义节点2.1.1 节点定义2.1.2 外形参数变量2.1.3 节点厚度2.1.4 节点的法线定义2.1.5 坐标系统的转换2.2 定义单元2.2.1 单元定义2.2.2 单元建立2.2.3 定义加筋2.2.4 定义钢筋作为一个单元属性2.2.5 方向2.3 定义表面2.3.1 表面:概述2.3.2 定义基于单元的表面2.3.3 定义基于节点的表面2.3.4 定义解析刚体表面2.3.5 对表面进行操作2.4 定义刚体2.4.1 刚体定义2.5 定义积分输出项2.5.1 积分输出项的定义2.6 定义不做结构材料的质量2.6.1 不做结构材料的质量定义2.7 定义分布2.7.1 分布的定义2.8 定义显示体2.8.1 显示体的定义2.9 定义一个装配2.9.1 定义一个装配2.10 定义矩阵2.10.1 定义矩阵第3章执行程序3.1 执行程序:概述3.1.1 执行ABAQUS程序:概述3.2 执行程序3.2.1 用于获得信息的执行程序3.2.2 用于ABAQUS/Standard和ABAQUS/Explicit的执行程序3.2.3 用于ABAQUS/CAE的执行程序3.2.4 用于ABAQUS/Viewer的执行程序3.2.5 用于Python的执行程序3.2.6 用于参数研究的执行程序3.2.7 用于ABAQUS HTML文件的执行程序3.2.8 用于许可证有效性的执行程序3.2.9 用于结果文件(.fil)的ASCII转化的执行程序3.2.10 用于连接结果文件(.fil)的执行程序3.2.11 用于查询关键词/问题数据库的执行程序3.2.12 用于获取例子input文件的执行程序3.2.13 用于用户自定义执行和子程序的执行程序3.2.14 用于input文件和输出数据库升级效用的执行程序3.2.15 用于生成输出数据报告的执行程序3.2.16 用于重启动分析连接输出数据库(.odb)的执行程序3.2.17 用于结合子结构输出的执行程序3.2.18 用于网络输出数据库文件连接器的执行程序3.2.20 用于将NASTRAN大批数据文件转化为ABAQUS中input文件的执行程序3.2.21 用于将PAM-CRASH输入文件转化为部分ABAQUS中input文件的执行程序3.2.22 用于将ABAQUS输出数据库文件转为NASTRAN Output2结果文件的执行程序3.2.23 用于和ZAERO交换ABAQUS数据的执行程序3.2.24 加密和解密ABAQUS输入数据的执行程序3.2.25 用于job执行控制的执行程序3.3 环境文件设置3.3.1 使用ABAQUS环境文件设置3.4 管理内存和硬盘资源3.4.1 在ABAQUS中管理内存和硬盘资源3.5 文件扩展定义3.5.1 通过ABAQUS使用文件扩展定义3.6 FORTRAN单位数3.6.1 通过ABAQUS使用的FORTRAN单位数第4章输出4.1 输出4.1.1 输出4.1.2 数据和结果文件的输出4.1.3 输出数据库的输出4.2 输出变量4.2.1 ABAQUS/Standard输出变量符4.2.2 ABAQUS/Explicit输出变量符4.3 后处理器4.3.1 后处理器第5章文件输出格式5.1 访问结果文件5.1.1 访问结果文件:概述5.1.2 结果文件输出格式5.1.3 访问结果文件信息5.1.4 用于访问结果文件的增效程序第6章分析程序6.1 介绍6.1.1 程序:概述6.1.2 一般的和线性扰动的程序6.1.3 多重荷载情况分析6.1.4 直接线性方程求解6.1.5 迭代线性方程求解6.2 静态应力/位移分析6.2.1 静态应力分析程序:概述6.2.2 静态应力分析6.2.3 特征值崩溃预测6.2.4 不稳定的崩塌和崩溃后分析6.2.5 准静态分析6.2.6 直接循环分析(已译)6.3 动态应力/位移分析6.3.1 动态分析程序:概述6.3.2 使用直接积分的隐式动态分析6.3.3 显示动态分析6.3.4 直接求解的稳定状态动态分析6.3.5 自然频率的提取6.3.6 复杂特征值的提取6.3.7 瞬时模态动态分析6.3.8 基于范数的稳定状态动态分析6.3.9 基于子空间的稳定状态动态分析6.3.10 响应谱分析6.3.11 随机响应分析6.4 稳定状态的运输分析6.4.1 稳定状态的运输分析6.5 热传播和温度-应力分析6.5.1 热传播分析程序:概述6.5.2 非耦合的热传播分析6.5.3 连续耦合的温度-应力分析6.5.4 全耦合的温度-应力分析6.5.5 绝热分析6.6 电分析6.6.1 电分析程序:概述6.6.2 耦合温度-电分析6.6.3 压电分析6.7 耦合多孔流体流动和应力分析6.7.1 耦合多孔流体扩散和应力分析(已译)6.7.2 地应力状态(已译)6.8 质量扩散分析6.8.1 质量扩散分析6.9 声学和振动分析6.9.1 声学、振动和耦合声(波)-结构分析6.10 ABAQUS/Aqua分析6.10.1 ABAQUS/Aqua分析6.11 退火6.11.1 退火程序第7章分析求解和控制7.1 求解非线性问题7.1.1 求解非线性问题7.1.2 接触迭代7.2 分析的收敛控制7.2.1 收敛和时间积分准则:概述7.2.2 普遍使用的控制参数7.2.3 非线性问题的收敛准则7.2.4 瞬态问题中的时间积分精度第8章分析技术:介绍8.1 介绍8.1.1 分析技术:概述第9章连续分析的技术9.1 重启动一个分析9.1.1 重启动一个分析9.2 输入和传递结果9.2.1 在ABAQUS分析中传递结果:概述9.2.2 在ABAQUS/Explicit和ABAQUS/Standard中传递结果9.2.3 将ABAQUS/Standard分析中的结果传递给另一个第10章模型提取10.1 子结构10.1.1 使用子结构10.1.2 定义子结构10.2 子模型10.2.1 子模型10.3 对称模型的生成,结果传递,循环对称模型的分析10.3.1 对称模型的生成10.3.2 将一个对称网格或一个部分三维网格的结果传递到完全三维网格10.3.3 分析存在循环对称的模型10.4 梁横截面网格划分10.4.1 梁横截面网格划分第11章特定目标的技术11.1 惯量解除11.1.1 惯量解除11.2 网格修改或置换11.2.1 单元和接触对的移除和重新激活(弹塑性理论)11.3 几何不完整11.3.1 在模型中引入一个几何不完整11.4 断裂力学11.4.1 断裂力学:概述(已译)11.4.2 围道积分评价11.4.3 裂缝扩展分析11.5 静水力学的流动模型11.5.1 模拟充满流体的空穴11.6 基于表面的流动模型11.6.1 基于表面的流体空穴:概述11.6.2 定义流体空穴11.6.3 定义流体的交换11.6.4 定义充气机11.7 质量数标度11.7.1 质量数标度11.8 稳定状态的探测11.8.1 稳定状态的探测11.9 平行执行11.9.1 ABAQUS中的平行执行11.9.2 ABAQUS/Standard中的平行执行11.9.3 ABAQUS/Explicit中的平行执行第12章自适应技术12.1 自适应技术:概述12.1.1 自适应技术12.2 ALE自适应网格划分12.2.1 ALE自适应网格划分:概述12.2.2 在ABAQUS/Explicit中定义ALE自适应网格划分区域12.2.3 ABAQUS/Explicit中的ALE自适应网格划分和重新绘图12.2.4 ABAQUS/Explicit中Eulerian自适应网格划分区域的模型技术12.2.5 ABAQUS/Explicit中ALE自适应网格划分的输出和诊断12.2.6 在ABAQUS/Standard中定义ALE自适应网格划分区域12.2.7 ABAQUS/Standard中的ALE自适应网格划分和重新绘图12.3 自适应重新网格划分12.3.1 自适应重新网格划分:概述12.3.2 误差指示器12.3.3 基于求解的网格划分尺寸12.4 网格划分置换后的连续分析12.4.1 网格-网格求解映射第13章扩展ABAQUS分析的功能13.1 联合仿真13.1.1 联合仿真:概述13.1.2 为联合仿真准备一个ABAQUS分析13.1.3 使用MpCCI联合仿真13.1.4 含有MADYMO的联合仿真13.2 用户子程序和增效程序13.2.1 用户子程序:概述13.2.2 可用的用户子程序13.2.3 可用的增效程序第14章设计敏感度分析14.1 设计敏感度分析14.1.1 设计敏感度分析第15章参数的研究15.1 脚本参数的研究15.1.1 脚本参数的研究15.2 参数的研究:命令15.2.1 为参数研究结合参数样本15.2.2 在参数研究中约束联合的参数值15.2.3 为参数研究定义参数15.2.4 执行参数研究设计分析15.2.5 聚集参数研究的结果15.2.6 为一个参数研究生成分析任务数据15.2.7 指定参数研究结果的来源15.2.8 创建一个参数研究15.2.9 报告参数研究的结果15.2.10 参数研究的样本参数第16章材料:介绍16.1 介绍16.1.1 材料库:概述(已译)16.1.2 材料数据的定义16.1.3 材料的结合行为16.2 一般属性16.2.1 密度第17章弹性力学性质(可以看看)17.1 概述17.1.1 弹性行为:概述(已译)17.2 线弹性17.2.1 线弹性行为17.2.2 无压缩或无拉伸17.2.3 平面应力各向正交异性失效测量17.3 多孔弹性17.3.1 多孔材料的弹性行为17.4 亚弹性17.4.1 亚弹性行为17.5 超弹性17.5.1 橡胶类材料的超弹性行为17.5.2 泡沫胶的超弹性行为17.6 Mullins效果17.6.1 橡胶类材料的Mullins效果17.6.2 泡沫胶的能量消散17.7 粘弹性17.7.1 时域粘弹性17.7.2 频域粘弹性17.8 滞后作用17.8.1 弹性体(人造橡胶)的滞后作用17.9 状态方程17.9.1 状态方程第18章非弹性力学性质18.1 概述18.1.1 非弹性行为(已译)18.2 金属塑性18.2.1 经典金属塑性(已译)18.2.2 承受循环荷载下的金属的模型(已译)18.2.3 率相关屈服18.2.4 率相关塑性:蠕变和膨胀18.2.5 退火和融化18.2.6 各向异性的屈服/蠕变18.2.7 Johnson-Cook塑性18.2.8 动态失效模型18.2.9 多孔金属塑性18.2.10 灰铸铁塑性18.2.11 两层粘塑性18.2.12 ORNL(Oak Ridge National Laboratory)本构模型18.2.13 变形塑性18.3 其他塑性模型18.3.1 扩展的Drucker-Prager模型18.3.2 修正的Drucker-Prager/CAP模型18.3.3 Mohr-Coulomb塑性18.3.4 临界状态(粘土)塑性模型18.3.5 可压碎的泡沫塑性模型18.4 有接缝的材料18.4.1 有接缝的材料模型18.5 混凝土18.5.1 混凝土涂抹开裂18.5.2 混凝土开裂模型18.5.3 混凝土塑性损伤(已译)第19章累积损伤和失效19.1 累积损伤和失效:概述19.1.1 累积损伤和失效(已译)19.2 延性金属的损伤和失效19.2.1 延性金属的损伤和失效:概述19.2.2 延性金属的损伤开始19.2.3 延性金属的损伤演化和单元移除19.3 加筋复合物的损伤和失效19.3.1 加筋复合物的损伤和失效:概述(已译)19.3.2 加筋复合物的损伤开始(已译)19.3.3 加筋复合物的损伤演化和单元移除(已译)第20章其他材料性质20.1 力学属性20.1.1 材料阻尼20.1.2 热膨胀20.2 热传播属性20.2.1 温度属性:概述20.2.2 传导性20.2.3 比热20.2.4 潜伏热20.3 声属性20.3.1 声媒介20.4 静水力学的流体属性20.4.1 静水力学的流体模型20.5 质量扩散属性20.5.1 扩散能力20.5.2 溶解性20.6 电属性20.6.1 电导率20.6.2 压电行为20.7 多孔流体流动属性20.7.1 多孔流体流动属性(已译)20.7.2 渗透性(已译)20.7.3 多孔体积模量(已译)20.7.4 吸附作用20.7.5 膨胀凝胶体20.7.6 湿度膨胀20.8 用户定义材料20.8.1 用户定义的材料力学行为20.8.2 用户定义的材料温度行为第21章单元:介绍21.1 介绍21.1.1 单元库:概述21.1.2 选择单元的维度21.1.3 对一个分析类型选择合适的单元21.1.4 截面控制21.1.5 根据单元-单元原理分配单元属性第22章连续单元22.1 多用途的连续单元22.1.1 固体(连续)单元22.1.2 一维固体(链接)单元库22.1.3 二维固体单元库22.1.4 三维固体单元库22.1.5 圆柱固体单元库22.1.6 轴对称固体单元库22.1.7 含有非线性、不均匀变形的轴对称固体单元22.2 无限单元22.2.1 无限单元22.2.2 无限单元库22.3 翘曲单元22.3.1 翘曲单元22.3.1 翘曲单元库第23章结构单元23.1 膜单元23.1.1 膜单元23.1.2 一般的膜单元库23.1.3 圆柱膜单元库23.1.4 轴对称膜单元库23.2 杆单元23.2.1 杆单元23.2.2 杆单元库23.3 梁单元23.3.1 梁模型:概述23.3.2 选择一个梁横截面23.3.3 选择一个量单元23.3.4 梁单元横截面的方向23.3.5 梁截面的行为23.3.6 在分析中使用一个梁截面积分来定义截面行为23.3.7 使用一个一般的梁截面来定义截面行为23.3.8 梁单元库23.3.9 梁横截面库23.4 框架单元23.4.1 框架单元23.4.2 框架截面属性23.4.3 框架单元库23.5 弯头单元23.5.1 有变形横截面的管和管弯头:弯头单元23.5.2 弯头单元库23.6 壳单元23.6.1 壳单元:概述23.6.2 选择一个壳单元23.6.3 定义传统壳单元的初始尺寸23.6.4 壳截面行为23.6.5 在分析中使用一个壳截面积分来定义截面行为23.6.6 使用一个一般壳截面来定义截面行为23.6.7 三维传统的壳单元库23.6.8 连续壳单元库23.6.9 轴对称壳单元库23.6.10 含有非线性、非轴对称变形的轴对称壳单元第24章惯性单元、刚体单元和电容单元24.1 点质量单元24.1.1 点质量24.1.2 质量单元库24.2 旋转惯量单元24.2.1 旋转惯性24.2.2 旋转惯性单元库24.3 刚体单元24.3.1 刚体单元24.3.2 刚体单元库24.4 电容单元24.4.1 点电容24.4.2 电容单元库第25章连接器单元25.1 连接器单元25.1.1 连接器:概述25.1.2 连接器单元25.1.3 连接器驱动25.1.4 连接器单元库25.1.5 连接类型库25.2 连接器单元行为25.2.1 连接器的行为25.2.2 连接器的弹性行为25.2.3 连接器的阻尼行为25.2.4 用于耦合行为的连接器功能25.2.5 连接器的摩擦行为25.2.6 连接器的塑性行为25.2.7 连接器的损伤行为25.2.8 连接器的停止和锁定25.2.9 连接器的失效模型第26章特定目标的单元26.1 弹簧单元26.1.1 弹簧26.1.2 弹簧单元库26.2 阻尼器单元26.2.1 阻尼器26.2.2 阻尼器单元库26.3 柔性接头单元26.3.1 柔性接头单元26.3.2 柔性接头单元库26.4 分布耦合单元26.4.1 分布耦合单元26.4.2 分布耦合单元库26.5 粘结单元26.5.1 粘结单元:概述26.5.2 选择一个粘结单元26.5.3 含有粘结单元的模型26.5.4 定义粘结单元的初始尺寸26.5.5 使用连续方法定义粘性单元的本构响应26.5.6 使用牵引-分离描述定义粘性单元的本构响应26.5.7 在粘结单元的缺口内定义流体的本构响应26.5.8 两维的粘结单元库26.5.9 三维的粘结单元库26.5.10 轴对称粘结单元库26.6 垫圈单元26.6.1 垫圈单元:概述26.6.2 选择一个垫圈单元26.6.3 在一个模型中包含垫圈单元26.6.4 定义垫圈单元的初始尺寸26.6.5 使用一个材料模型定义垫圈行为26.6.6 直接使用一个垫圈行为模型定义垫圈行为26.6.7 两维垫圈单元库26.6.8 三维垫圈单元库26.6.9 轴对称垫圈单元库26.7 表面单元26.7.1 表面单元26.7.2 一般的表面单元库26.7.3 圆柱表面单元库26.7.4 轴对称表面单元库26.8 静水力学的流体单元26.8.1 静水力学的流体单元26.8.2 静水力学的流体单元库26.8.3 流体链接单元26.8.4 静水力学的流体连接库26.9 管座单元26.9.1 管座单元26.9.2 管座单元库26.10 线弹簧单元26.10.1 模拟壳中部分贯通裂缝的线弹簧单元26.10.2 线弹簧单元库26.11 弹-塑性接头26.11.1 弹-塑性接头26.11.2 弹-塑性接头单元库26.12 拉链单元26.12.1 拉链26.12.2 拉链单元库26.13 桩-土单元26.13.1 桩-土相互作用单元26.13.2 桩-土相互作用单元库26.14 声学的界面单元26.14.1 声学界面单元26.14.2 声学界面单元库26.15 用户自定义的单元26.15.1 用户自定义的单元26.15.2 用户自定义的单元库第27章施加条件27.1 概述27.1.1 施加条件:概述27.1.2 幅值曲线27.2 初始条件27.2.1 初始条件27.3 边界条件27.3.1 边界条件27.4 荷载27.4.1 施加荷载:概述27.4.2 集中荷载27.4.3 分布荷载27.4.4 温度荷载27.4.5 声荷载27.4.6 多孔流体流动27.5 指定装配荷载27.5.1 指定装配荷载27.6 预先确定场27.6.1 预先确定场第28章约束28.1 概述28.1.1 运动约束:概述28.2 多点约束28.2.1 线性约束方程28.2.2 一般多点约束28.2.3 运动的耦合约束28.3 基于表面的约束28.3.1 网格打结约束28.3.2 耦合约束28.3.3 壳-固体耦合28.3.4 不依赖网格的扣件28.4 植入单元28.4.1 植入单元28.5 单元释放终点28.5.1 单元释放终点28.6 过约束检查28.6.1 过约束检查第29章定义接触的相互作用29.1 概述29.1.1 接触相互作用分析:概述29.2 在ABAQUS/Standard中定义接触29.2.1 在ABAQUS/Standard中定义接触对29.2.2 ABAQUS/Standard接触对的接触公式29.2.3 ABAQUS/Standard接触对的强制约束方法29.2.4 在ABAQUS/Standard中模拟接触的干涉配合29.2.5 在ABAQUS/Standard接触对中调整初始面的位置和指定初始清空29.2.6 移除/重新激活ABAQUS/Standard接触对29.2.7 在ABAQUS/Standard中定义打结接触29.2.8 延伸主面和滑移线29.2.9 如果子结构存在的接触模拟29.2.10 如果不均匀-对称单元存在的接触模拟29.2.11 在ABAQUS/Standard中模拟接触时普遍存在的困难29.2.12 在ABAQUS/Standard中调整接触控制29.3 在ABAQUS/Explicit中定义一般接触29.3.1 定义一般接触的相互作用29.3.2 一般接触的表面属性29.3.3 一般接触的接触属性29.3.4 一般接触的接触公式29.3.5 一般接触的初始过闭合的分解和指定初始清空29.3.6 一般接触的接触控制29.4 在ABAQUS/Explicit中定义接触对29.4.1 在ABAQUS/Explicit中定义接触对29.4.2 ABAQUS/Explicit接触对的表面属性29.4.3 ABAQUS/Explicit接触对的接触属性29.4.4 ABAQUS/Explicit接触对的接触公式29.4.5 在ABAQUS/Explicit接触对中调整初始面的位置和指定初始清空29.4.6 在ABAQUS/Explicit中使用接触对算法模拟接触时普遍存在的困难第30章接触属性模型30.1 接触力学属性30.1.1 接触力学属性:概述30.1.2 接触的压力-过闭合关系30.1.3 接触阻尼30.1.4 接触阻塞30.1.5 摩擦行为30.1.6 用户自定义的界面本构行为30.1.7 贯入式压力荷载30.1.8 松解面的相互作用30.1.9 易碎的结合30.2 温度接触属性30.2.1 温度接触属性30.3 电接触属性30.3.1 电接触属性30.4 多孔流体接触属性30.4.1 多孔流体接触属性第31章ABAQUS/Standard中的接触单元31.1 模拟接触的单元31.1.1 模拟接触的单元31.2 缺口接触单元31.2.1 缺口接触单元31.2.2 缺口单元库31.3 管-管接触单元31.3.1 管-管接触单元31.3.2 管-管接触单元库31.4 滑动线接触单元31.4.1 滑动线接触单元31.4.2 轴对称滑动线接触单元库31.5 刚体表面接触单元31.5.1 刚体表面接触单元31.5.2 轴对称刚体表面接触单元库第32章在ABAQUS/Standard中定义空腔辐射32.1 定义空腔辐射32.1.1 空腔辐射。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2.1.11 Collapse of a stack of blocks with general contactProduct: Abaqus/ExplicitThis example illustrates the use of the general contact capability in a simulation involving a large number of contacting bodies. The general contact algorithm allows very simple definitions of contact with very few restrictions on the types of surfaces involved (see “Defining general contact interactions in Abaqus/Explicit,” Section 35.4.1 of the Abaqus Analysis User's Manual).Problem descriptionThe model simulates the collapse of a stack of blocks. The undeformed configuration of the model is shown in Figure 2.1.11–1. There are 35 blocks, and each block is 12.7 × 12.7 × 76.2 mm (0.5 × 0.5 × 3 inches) in size. The blocks are stacked on a rigid floor. The stack is subjected to gravity loading. It is assumed that a key block near the bottom of the stack has been removed just before the start of the analysis, initiating the collapse.Each block is modeled with a single C3D8R element. The use of a coarse mesh highlights the edge-to-edge contact capability of the general contact algorithm, because the majority of the block-to-block interactions do not result in penetrations of nodes into faces.Two different cases are analyzed. In the first analysis the blocks are rigid. In the second analysis the blocks are deformable. In the latter case, the material of the block is assumed to be linear elastic with a Young's modulus of 12.135 GPa (1.76 × 106 Psi), a Poisson's ratio of 0.3, and a density of 577.098 kg/m3 (5.4 × 10–5 lb s2/in4). Only the density is relevant for the analysis assuming rigid blocks. In addition, ENHANCED hourglass control is used for the deformable analysis. The rigid floor is modeled as a discrete rigid surface using a single R3D4 element.This model involves a large number of contacting bodies. The general contact capability greatly simplifies the contact definition, since each of the 595 possible block-to-block pairings does not need to be specified individually. The general contact inclusions option to automatically define an all-inclusive surface is used and is the simplest way to define contact in the model. Coulomb friction with a friction coefficient of 0.15 is assumed between the individual blocks and between the blocks and floor. The general contact property assignment is used to assign this nondefault contact property.By default, the general contact algorithm in Abaqus/Explicit accounts foredge-to-edge contact of perimeter edges on structural elements. Geometric featureedges of a model can also be considered for edge-to-edge contact by the general contact algorithm; including the geometric feature edges is crucial in this analysis. A cutoff feature angle of 20° is specified for the feature angle criterion of the surface property to indicate that all edges with feature angles greater than 20° should be considered for edge-to-edge contact. The feature angle is the angle formed between the normals of the two facets connected to an edge.The magnitude of the gravity loading is increased by a factor of 10 to facilitate demonstration of the edge-to-edge contact capability with a short analysis time. The analysis is performed for a period of 0.15 seconds. For the analysis with rigid blocks there is no deformable element available in the model to control the stable time increment. A fixed time increment of 1 × 10–6 seconds is specified for this purpose, which is similar to the time increment used by the analysis with deformable blocks. The time increment chosen for the analysis with rigid blocks will affect the penalty stiffness used by the contact algorithm since the penalty stiffness is inversely proportional to the time increment squared.Results and discussionResults are shown for the rigid body case. Results for the deformable case are very similar to the rigid model results.Figure 2.1.11–2 shows the displaced shape of the block assembly after 0.0375 seconds. The stack of blocks has started to collapse under gravity loading. Figure2.1.11–3 shows a close-up view of the collapsing blocks after 0.1125 seconds. This figure clearly shows that the geometric feature edges of individual blocks contact each other during collapse. Figure 2.1.11–4 shows the final configuration of the blocks. The stack has collapsed completely on the rigid surface.Input filesblocks_rigid_gcont.inpInput file for the rigid body analysis.blocks_rigid_assembly.inpExternal file referenced by the rigid body analysis.blocks_deform_gcont.inpInput file for the deformable analysis.blocks_deform_assembly.inpExternal file referenced by the deformable analysis.FiguresFigure 2.1.11–1 Initial configuration of the stack of blocks.Figure 2.1.11–2 Displaced shape after 0.0375 s.Figure 2.1.11–3 Close-up view of the collapsing blocks after 0.1125 s.Figure 2.1.11–4 Final configuration of the model.2.1.11 一摞积木在通用接触下的倒塌分析Product:Abaqus/Explicit这个例子说明通用接触的使用在涉及大量接触物体倒塌的模拟中应用。