ADINA二次开发培训(全面)
ADINA中文培训——adina training3

时间推进格式
多种方案
隐式程序
亚得科技-ADINA中国
显式程序
显式求解
在空间离散单元与隐式相同 时间推进采用显式中心差分格式 只能求解时间相关问题 有条件稳定:即时间步长要小于临界时间步长
亚得科技-ADINA中国
十 后处理
1 云图 修改云图颜色 等值线图 Smooth图 2 时间步控制 3 矢量图 4 接触反力和支反力图 5 同时显示多个结果图 6 列表显示位移 7 节点或时程曲线 8 画合力曲线(combine) 9 生成动画 10 输出图片 11 频域-时域曲线转换
Force N 1e-02 N 1e+04 N 1e+10 N
Stress Pa
Energy Joule
ρ (steel)
7.83e+03 7.83e-03 7.83e-03 7.83e-03 7.83e-06 7.83e+00 7.83e+00 7.83e-03 7.83e-03 7.83e-09 7.32e-04 1.52e+01
第二天
亚得科技-ADINA中国
材料模式本构
线弹 非线弹 塑性 热本构 蠕变 橡胶/泡沫 粘弹 垫片材料 势流体 岩土 混凝土 …
亚得科技-ADINA中国
单元算法
Truss Beam 2D Solid 3D solid Plate/shell Sping Pipe Fluid … MITC4 Shell/Layers
亚得科技-ADINA中国
模型DEBUG 模型
如果模型没有正常求解就退出,参件Jobname.out和Jobname.msg文件中 的信息;前者记录求解过程的所有内容; 如果程序正常求解完毕,不要认为EVERYTHING IS OK,用后处理将所 有结果审查一遍; 如果经过多次检查, 始终有一个错误(问题)存在,设计一个简单模型专门 测试这个问题. 或者, 将大模型中的部分删除,剩余的东西要表达这个问 题,然后继续测试; 如果实在找不出问题所在: 1 让ADINA工程师解决 2 是程序的BUG,转告ADINA中国工程师.
32 adina培训资料(版本2)element types

[I]
[A] [B] [C] [D] [E]
[F] [G] [H]
亚得科技-ADINA中国
单元类型
ADINA-杆单元或桁架单元(Truss)
说明:在其它单元的类型中,关于 [B]:几何非线性 [D]:质量特性 [E]:初始应变 [H]:蠕变时间 几项的定义与Truss相同,因此不再 重复解释。 [I] [A] [B] [C] [D] [E] [F] [G] [H]
说明: 1. Isobeam单元只能定义矩形截面梁; 2. 适用于小位移、大位移、小应变; 3. 网格划分需要指定辅助节点确定R、S、T坐标系; 4. 在线性分析中,哈密顿梁的刚度矩阵为闭式的,而isobeam采用数值积分的方法计算, 因此在直线梁的线性分析,2-node哈密顿梁单元比2-node isobeam单元更有效; 5. Isobeam单元主要用于模拟曲线梁、壳的刚性支撑、大变形的梁、轴对称荷载作用下的 轴对称壳 ;
亚得科技-ADINA中国
单元类型
[A]:单元子类型 包括: 轴对称,平面应力,平面应变 ,3-D平面应力,通用平面应变; 3-D平面应力一般用于从2D网格通过Sweep或 Revolution方式生成3D网格时的初始网格; [B]:单元选项 选项包括: Porous Media(多孔介质) User-Coded(用户定义) Simple gasket(简单垫片材料) General Gasket()通用垫片材料 只有当模型中有单元选项设置为多孔介质时 ,才能在Model/Materials/Porous Property中设置其参数; [C]:3D Plane Strain单元的参考点 可以是Point或Node;
亚得科技-ADINA中国
单元类型
[A]:选择积分算法和积分点 一般为缺省设置。当进行大应变分 析时,如汽车碰撞和板成形模拟中 ,厚度t方向的积分点个数为3~5个 ; [B]:层数 可分层,并且在 Meshing/Element/Element Layer中 给每层定义不同的材料性质多层材 料如仿弹玻璃等;或给每层定义增 强纤维的体积分数等如纤维增强复 合材料。详见Meshing中单元管理一 节内容; [C]:单元失效准则 定义其各种失效形式,包括复合材 料失效的Tsai-Hill,Tsai-Wu,TsaiHashin等;详见理论手册2.7; [D]:开关,是否计算中性面上的等 效力和力矩;
71 adina培训资料(版本2)Post-processing

3. 一些后处理操作需点选节点(或单元),选择节点前要先点击图标 点符号;
打开节
亚得科技-ADINA中国
后处理
1 云图 修改云图颜色 等值线图 Smooth图 2 时间步控制 3 同时显示多个结果图 4 矢量图 5 接触反力和支反力图 6 显示相对位移 ( relative displacement ) 7 节点或时程曲线 8 列表显示结果results 9 极值列表显示 10 合力计算(combination)
[A]:显示结果的范围 [B]:轴向的定义 [C]:轴力符号 [D]:剪力方向 [E]:弯矩方向 [F]:结果显示符号尺寸 [G]:颜色选择 [G] [F] [A]
后处理功能-修改Beam, Truss单元表
[B]
[C]
[D]
[E]
亚得科技-ADINA中国
后处理
后处理功能-切片显示
切片显示 显示变形网格 变形放大显示 显示未变形 网格
[B]
[A] [C]
点击"切片显示"图标,出现右侧窗口; [A]:选择Cutting Surface,如果从切片 状态回到正常状态,选择None; [B]:切片平面方向和位置 [C]:如果要用一系列等间隔平面切模型 ,则输入相邻两个平面距离值; [D]:控制切开模型的显示方式; [D]
亚得科技-ADINA中国
后处理
后处理功能-切片显示
切片显示 显示变形网格 变形放大显示 显示未变形 网格
[A]
点击"切片显示"图标,出现右侧窗口; [A]:选择Iso Surface,则可以画三维 模型的等值面图;如果从切片状态回到 正常状态,选择None;
ADINA自定义材料开发原理与实际操作

3 二次开发实际操作 ................................................................................................. 8
3.1 3.2 3.3 3.4 自定义材料开发综述 ...................................................................................................... 8 开发环境 .......................................................................................................................... 8 开发步骤 .......................................................................................................................... 9 注意事项 .......................................................................................................................... 9
1
ADINA 中二次开发
图 2-1 小增量步示意图
图 2-1 为 ADINA 结构理论手册中的示意图,Δ eij 即是程序中用到的 DEPS。按 照上述说明,两种增量步之间可以理解为以下关系:
= t INTER D E P= S STRAIN INTER
ADINA培训课程

u
亚得科技-ADINA中国
PDF 文件使用 "pdfFactory Pro" 试用版本创建
Newton-Raphson方法求解非线性问题
Newton-Raphson 法迭代求解使用下列方程: [KT]{∆u} = {Fa} - {Fnr} 这里:
亚得科技-ADINA中国
PDF 文件使用 "pdfFactory Pro" 试用版本创建
时间函数 • 任何载荷都必须由时间函数定 义其随时间的变化; • 缺省的时间函数是随时间没有 变化; • 任何时间步的载荷增量由时间 函数和时间步长共同确定;
Load
1.0
time 例子:荷载为指定速度 ,速度为10mm/s,时间 步长恒定为0.1,共计算 150步,则最终的压下 距离为150mm。
增加删除图标 定义快捷键
所谓Line Draw-Down图标的控制
定义宏命令
• ADINA的国
PDF 文件使用 "pdfFactory Pro" 试用版本创建
• Graph图形显示区的调整Edit
GDI与OpenGL两种方式 设置AUI内存 设置背景颜色,窗口比例 设置图形环境参数
εy
ε
亚得科技-ADINA中国
PDF 文件使用 "pdfFactory Pro" 试用版本创建
塑性多线性模型
都使用多个线段的应力-应变曲线来模拟随动强化效应。使用von Mises 屈服准则, 包括各向同性和随动硬化。
输入弹性模量和应力-应变数据点定义材料参数 : •每条应力-应变曲线必须 用同一组应变值; •曲线的第一个点必须 与弹性模量一致; •每一段的斜度不能超过弹性模量(不允许负斜度); •对于超过输入曲线末端的应变值,假设为理想塑性材 料。
ADINA软件中邓肯—张模型的二次开发与应用

ADINA软件中邓肯—张模型的二次开发与应用
陈志坚;江涛;陈松
【期刊名称】《水电能源科学》
【年(卷),期】2007(25)2
【摘要】根据AD INA提供的用户子程序,用FORTRAN编写了邓肯—张E-v本构模型,给出了模型开发的方法和步骤,并用常规三轴压缩试验对模型进行了检验。
检验结果表明,开发的模型正确合理、计算精度较理想,可应用于实际工程计算以解决复杂的非线性土工问题。
【总页数】3页(P65-67)
【关键词】ADINA;邓肯-张模型;二次开发
【作者】陈志坚;江涛;陈松
【作者单位】河海大学土木工程学院
【正文语种】中文
【中图分类】TU411;TP311.5
【相关文献】
1.邓肯-张E-B模型的ANSYS二次开发及应用 [J], 孙明权;陈姣姣;刘运红
2.Duncan-Chang E-υ模型在ADINA软件中的开发与应用 [J], 陈志波;朱俊高;陈浩锋
3.邓肯-张模型开发及其在面板坝计算中的应用 [J], 宋志宇;李斌;董莉莉;梁春雨
4.邓肯-张EB模型参数求解的二次优化法 [J], 陈立宏
5.一种改进邓肯张模型及其在土石坝数值模拟中的应用 [J], 邵东琛
因版权原因,仅展示原文概要,查看原文内容请购买。
ADINA中的用户二次开发特征介绍

易于扩充功能的分析平台(ADINA二次开发环境)任何一个商业化软件,都不可能完全满足所有的工程问题以及所有用户的需求,尤其是研究院所、高校等一些处在科学研究最前沿的单位。
为此ADINA提供丰富完善的二次开发功能,允许用户方便地自定义各种用户功能,常见如本构算法、材料破坏准则、接触摩擦形式;断裂力学判据和开裂扩展规律、用户边界条件/荷载等等。
材料本构用户自定义常见的邓肯-张E-B模型,基于Harding的动本构模型;双屈服面动本构模型;各种冻土本构模型等等;每种材料开发时可以直接耦合温度自由度,并且能够同时定义材料(单元)失效条件;接触状态用户自定义土体的接触摩擦状态比机械零件之间的摩擦状态(机械行业往往采用双线性的库仑摩擦模型)复杂得多,ADINA中除了双线性的库仑摩擦模型外,还提供多种摩擦模型,得到变参数摩擦状态。
除此之外,ADINA的用户接口提供,定义摩擦系数随任意接触状态参数如接触压力、相对滑动速度、甚至温度的变化规律。
ADINA理论手册中介绍了常见的13种形式的接触摩擦模型,供用户进行自定义并使用,如下:C = 1 means constant friction coefficientC = 2 means Friction Law #1C = 3 means Friction Law #2C = 4 means different static and dyanamic friction coefficientsC = 5 means friction coefficient varies with sliding velocityC = 6 means anisotropic friction lawC = 7 means linearly varying coefficient of frictionC (linear with respect to normal force)C = 8 means a time varying coefficient of frictionC = 9 means coefficient of friction varies with coordinatesC = 12 means Friction Law #1C = 13 means Friction Law #2材料破坏准则用户自定义对于任何一种材料(包括用户开发材料),用户可以提供自定义的破坏准则。
ADINA中文培训parasolid

亚得科技-ADINA中国
Part two 布尔运算
Boolean Operator:
• Intersect
亚得科技-ADINA中国
Part two 布尔运算
Boolean Operator:
• Subtract
• 相减后可能会生成多个体
亚得科技-ADINA中国
Part two 布尔运算
亚得科技-ADINA中国
Part three 修改体
Boolean Modifier:
• Blend
• Multiple Edges Blended by a Constant Radius • A Single Edge Blended by Two Radius
亚得科技-ADINA中国
Part three 修改体
定义体:
• CLINDER:
亚得科技-ADINA中国
Part one 创建体(Body)
定义体:
• CLINDER:
亚得科技-ADINA中国
Part one 创建体(Body)
定义体:
• SPHERE:
• 指定坐标或点为球心 • 指定半径或球面上点
亚得科技-ADINA中国
定义体:
Part one 创建体(Body)
Boolean Modifier:
• Blend Multiple Faces Blended by a
Constant Radius
亚得科技-ADINA中国
Part three 修改体
Boolean Modifier:
• Chamfer
亚得科技-ADINA中国
Part three 修改体
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
FOCUSED ON EXCELLENCE
支持的二次开发工作
材料本构开发(2D:ovl30u*.f; 3D: ovl40u*.f,分别10个源文件) 摩擦算法开发(ovlusr.f) 单元算法开发 (ovl160u.f) 单元失效准则开发(ovlusr.f) 结构断裂力学和裂纹扩展规律开发(2D:ovl30u.f;3D: ovl40u.f) 压电材料/固结材料开发(2D:ovl30u.f;3D: ovl40u.f) 热本构开发(at00u.f) 流体本构开发(ovlf0u.f) 边界、载荷条件开发(ovl170u.f)
FOCUSED ON EXCELLENCE
开发详细过程
• 材料二次开发,与模型中采用的单元类型相关, 因此2D和3D的应用是分开的; • 前面提到,二次开发需先选择一个已经存在的源 文件,可以是任何一个,因为对于开发来说其提供 的入口条件是相同的; • 现在要自己开发一个最简单的线弹性材料,用于 3D模型,下面是详细的开发过程; • 由于ovl40u_vp1.f中详细介绍了所以参数的含义, 因此选择这个源文件进行修改。
εx应变 εy应变 εz应变 εxy应变 εyz应变 εzx应变
c c C CTI(1)=young's modulus CTI(2)=poisson's ratio CTI(n) =…
C
CTDD(98)=…
C c c c c
常用的其它参量
TIME = 当前时间 DT = 当前时间步长 TEMP1 = 当前温度
子程序说明: 传递的变量很多,完全满足绝大多数开发所需 常用变量:STRAIN(应变) STRESS(应力) TIME(当前时间) DT(当前时间步长) TMP1(当前温度)…
FOCUSED ON EXCELLENCE
开发详细过程
C ...................................................................... C. C . ELEMENTS: C. . Linear Elastic Model for 3-D SOLID ELEMENTS . .
这意味着我们的代码都是嵌入到已经指定的那些源文件中,而不是任意生成 一个新的源文件。
• 开发材料本构时,需要编辑Makefile.adusr文件,修改其 中的参数:
MAT2D_OBJ = ovl30u_pl1.obj (2D模型) MAT3D_OBJ = ovl40u_pl1.obj (3D模型) 这里,ovl30u_pl1.obj和ovl40u_pl1.obj分别指向ovl30u_pl1.f和ovl40u_pl1.f源 文件,原因是可以进行编译的源文件很多,而ADINA每个只能指定一个进入 DLL。
FOCUSED ON EXCELLENCE
材料本构开发入口条件
未知变量
STRESS(1)= σ xx STRESS(2)= σ yy STRESS(3)= σ zz STRESS(4)= σ xy STRESS(5)= σ yz STRESS(6)= σ zx
历史变量
ARRAY(LGTH1) 提供多于60个可以输出的实数型历史变量,常用 于做如屈服应力、塑性应变、应力水平等数据的存储; IARRAY( LGTH2 ) 提供多个全局整型历史变量,往往用于程序控制 或其它用途;
FOCUSED ON EXCELLENCE
开发详细过程
C*I K E Y = 2 C*I INTEGRATION OF ELEMENT STRESSES ( CALCULATE STRESS(6) ) 2 CONTINUE C*I * * * C C YM=CTI(1) PR=CTI(2) g2=YM/(1.0+PR) g=0.5*g2 C davg=(-STRAIN(1)-STRAIN(2)-STRAIN(3))*C1D3 p=-davg*YM/(1.0-2.0*PR) C STRESS(1)=p+g2*(STRAIN(1)+davg) STRESS(2)=p+g2*(STRAIN(2)+davg) STRESS(3)=p+g2*(STRAIN(3)+davg) STRESS(4)=g*STRAIN(4) STRESS(5)=g*STRAIN(5) STRESS(6)=g*STRAIN(6) C RETURN INSERT
FOCUSED ON EXCELLENCE
开发代码是如何进入到求解器的? • 任何程序段是通过以上列出的动态链接库 将用户开发的内容传送给求解器的:
ADINA-Structure模块:adusr.dll->adina.exe ADINA-THermal模块:atusr.dll ->adinat.exe ADINA-CFD模块:afusr.dll ->adinaf.exe ADINA-FSI模块:adfusr.dll ->adfsi.exe ADINA-TMC模块:adtusr.dll ->adtmc.exe
变量声明: 采用ADINA开发实例中的变量声明即可。
FOCUSED ON EXCELLENCE
开发详细过程
C GO TO (1,2,3,4), KEY C*I C*I C*I C*I C*I C*I C*I 1 CONTINUE C*I DO 11 I = 1,LGTH1 11 ARRAY(I) = 0.0 DO 12 I = 1,LGTH2 12 IARRAY(I) = 0 C RETURN 初始化两个数组,用于存放 历史变量(ADINA在后处理时 像位移、应力等变量显示这 些数值),LGTH1和LGTH2是 在AUI中输入的整数,有缺省 值。 执行第一段程序 INITIALIZE COMPONENTS OF REAL AND INTEGER WORKING ARRAYS ( INITIALIZE ARRAY(60) AND IARRAY(2) ) K E Y = 1 将整个程序分为四段
FOCUSED ON EXCELLENCE
材料本构二次开发
FOCUSED ON EXCELLENCE
材料本构开发入口条件
已知参数
c c c c c c strain(1)= strain(2)= strain(3)= strain(4)= strain(5)= strain(6)=
AUI输入参数
FOCUSED ON EXCELLENCE
材料本构开发入口条件
SUBROUTINE CUSER3 [DLLEXPORT] (NG,NEL,IPT,IDEATH,ITE,IUPDTL, 1 STRESS,EPS,STRAIN,DEPS,DEPST,THSTR1,THSTR2, 2 KTR,INTER,SCP,ARRAY,LGTH1,IARRAY,LGTH2, 3 D,ALFA,CTD,ALFAA,CTDD,CTI,TMP1,TMP2, 4 TIME,ETIMV,ETIMV2,DT,PHIST,PRST,RN,PHIST1, 5 DPSP,TGRAD,INTEG,ISUBM,INDNL, 6 DP,NELP,DPJE1D,DPJE2D,AKAPPA,PBAR, 7 NNODE,NODNUM,XYZ,DCA,IOUT,KEY) C C ......................................................................
C ......................................................................
程序中的注释部分: 一般程序中主要使用的中间变量给予说明; 如上面说明YM在程序中代表弹性模型;PR代表波松比;
FOCUSED ON EXCELLENCE
开发详细过程
• 当启动任何求解器时,DLL中的程序自动成 为求解器的一部分
FOCUSED ON EXCELLENCE
开发代码和输入数据是如何匹配的? •开发代码中指定了各种变量,包括所谓已知变量和未知变 量,已知变量的一部分是用户通过界面的相应菜单输入 的,如本构二次开发过程中的材料参数,控制参数等等
例如在材料开发时: 已知变量:应变--程序计算确定数值,出现在等式的右侧 弹性模量、波松比--用户输入参数,程序从dat模型文件识别读入 未知变量:应力--程序计算确定,出现在等式的左侧
FOCUSED ON EXCELLENCE
DLL是如何编译的? • Fortran源代码修改后,需要重新编译,生成一个新的DLL 文件。用这个文件替代求解器当前目录(\x32或者\x64子 目录下)的相应文件就实现了替换。DLL文件是如何编译的 呢?
每个文件都是用一个ADINA指定的文件进行编译: ADINA-Structure模块:Makefile.adusr->adusr.dll->adina.exe ADINA-THermal模块: Makefile.atusr-> atusr.dll ->adinat.exe ADINA-CFD模块: Makefile.afusr-> afusr.dll ->adinaf.exe ADINA-FSI模块: Makefile.adfusr-> adfusr.dll ->adfsi.exe ADINA-TMC模块: Makefile.adtusr-> adtusr.dll ->adtmc.exe
• 编译过程:
..\DF98\bin\dfvars(启动编译环境) nmake /f Makefile.<usr> (编译生成DLL文件) FOCUSED ON EXCELLENCE
DLL是如何编译的? • Makefile文件已经编写了有关路径、指定源代码文件和 OBJ文件等指令,因此除了材料本构可能需要修改外,其它 的开发不需修改,自Байду номын сангаас都是包括在内的。