radioss接触非线性分析
结构分析介绍radioss-intro

suoxz@free.fr 28th Oct 2003
4
高速撞击模拟
λ λ
高速撞击 : 正面, 后面, 侧面…
可修复性验证 λ 行人碰撞伤残评估 λ 驾驶员伤残评估 λ 安全装置(气囊, 安全带)有效性检测 λ 汽车翻滚模拟 ...
40%覆盖面的正面撞击 汽车与卡车的撞击
静态汽车受侧面运动物的撞击 后部撞击模拟
心血管中的生物力学
当用于流固偶合问题时,比如生物力学领域里血液在 心血管中的流动问题,RADIOSS软件也显示出了强大 的生命力,可逼真地模拟血液在可变形壁管中的流动, 以此来分析涉及心脏阀门、血管扩张 、动静脉所承受 的压力。
汽车正面撞击时 头部的运动方式
RADIOSS模拟人工心脏阀门中的血液流动过程Biblioteka 汽车撞击实验的独立测试机构
美国 (NHTSA 正面56 km/h和侧面62 km/h,IIHS 正面64 km/h), 欧共体(EuroNCAP 正面56 km/h和侧面62 km/h) : 根据人体的伤残程度,将撞击结果按所获星数分为5 个等级,其中5星为顶级。
suoxz@free.fr 28th Oct 2003
suoxz@free.fr 28th Oct 2003
10
RADIOSS HUMOS : 欧洲第一个驾驶位置的全柔性仿真人体模型
建立步骤 - 获取人体几何参数;
- 建立三维CAD模型; - 对各个器官划分有限元网格; - 确定器官的材料特性; - 联结不同器官、建立整体模型; - 有限元模型的实验验证; - 建立人体模型数据库(不同位置和姿态的男性、女性模型)
人体几何参数来自于医学切片
最初的本意是将这两种模型用于创伤生物力学,分析人体胸部、下肢、盆骨在外力打 击下的伤残程度。目前它们也被用于生物医学领域,让人造假肢的力学性能尽可能地 接近人类肢体。
radioss在线性动力学分析中的应用

Radioss提供了丰富的可视化工具,能够直 观地展示线性动力学分析结果,方便用户 进行结果分析和优化。
限制与挑战
01
02
03
适用范围
Radioss主要适用于线性 动力学分析,对于非线性 动力学问题可能存在局限 性。
参数设置
在进行线性动力学分析时, 需要合理设置材料参数、 边界条件等,否则可能导 致分析结果不准确。
展望
进一步优化算法
随着技术的进步,未来可以进一步优化Radioss的算法,提高其计算 效率和准确性。
扩展到非线性动力学分析
目前,Radioss主要应用于线性动力学分析。未来可以尝试将其扩展 到非线性动力学分析领域,以适应更广泛的问题。
与其他软件的集成
为了更好地满足用户需求,未来可以将Radioss与其他工程软件进行 集成,实现数据共享和交互。
02
线性动力学主要关注系统在平衡状态附近的微小振动,即线性
振动。
线性动力学适用于分析各种工程结构的振动特性,如桥梁、建
03
筑和机械系统等。
线性动力学模型建立
根据实际系统的物理 特性和运动规律,建 立线性动力学模型。
模型参数需要根据实 际测量或实验数据进 行标定。
模型通常由微分方程 或差分方程表示,描 述系统各质点的运动 规律。
04 Radioss在线性动力学分 析中的优势与限制
优势
高效性
精确性
Radioss具有高效的计算能力,能够快速进 行线性动力学分析,缩短了仿真时间,提 高了工作效率。
Radioss采用了先进的算法和数值方法,能 够提供高精度的线性动力学分析结果,为 设计决策提供可靠依据。
可扩展性
可视化功能
Radioss支持多种材料模型和边界条件,能 够根据实际需求进行定制化设置,满足不 同领域的线性动力学分析需求。
耦合与接触非线性在铁道罐车车体结构分析中的应用

用来装 运 流体 。设 计 安全 、可 靠 的 罐 车 车体 是 保 证 运 输安全 的首要前 提 。本 次 分 析 的罐 车 车 体 主要 由罐 体 装配 和底 架装 配等 组 成 。罐 体 由封 头 、简 体 、上 鞍 等
零部 件组 焊而成 ,采 用 直筒 结 构 。封 头 为标 准椭 圆形
维普资讯
第 2 第6 7卷 期 2 0 年 l 月 0 7 2 文 章 编 号 :10 —74 (07 6—02 —0 08 82 20 )0 03 3
铁 道 机 车 车 辆
RA IW AY D C0MaⅡ VE & CAR I l I
Vo . 7 No. 12 6
摘 要 利用 A S S 限元计算 软件 ,采用接触非 线性与耦 合 自由度两 种方法 对亚硫 酸罐 车车体进 行有 限元 NY 有
结构分析 ,并就两种 方法 的计算结 果作了详细的对 比。通 过对 比分析表 明,两 种分析 方法 的计 算结果很 接 近 , 最大误差为 1 . %,基本能满足一般 的工程需要 。两 种方法在铁道罐车 的结构分析 中都是切实可行的 。 47 7 关键词 有 限元分析 ;接 触非 线性 ;耦合 ;罐车
维普资讯
铁 道 机 车 车 辆
第 2 卷 7
比 ,拉 格朗 日算法 不 易 引 起病 态 条 件 ,对 接 触 刚 度 的 灵敏 度较小 ,然 而 ,在 有 些分 析 中 ,增 强 的拉 格 朗 日
加 载过 程 中材 料 的变 形 ,点 接触 可能会 发 展为 面接 触 ,
这 都体 现 出接 触 问题 的 高 度 复 杂 性 。在 计 算 过 程 中 , 需 要通 过材 料参 数 、载 荷 、边 界条 件 等 因 素 ,进行 综 合 分析后 才 能判 断或 确定 当前 的接 触状 态 。 ()接 触 条件 往往 具 有 非 线 性 特 征 。接 触 条 件包 2 括 :①接 触 物体不 可 相 互 侵 入 ;② 接 触 力 的法 向力 分 量 只能是 压力 ;③ 切 向接触 的摩 擦 条件 。
第六章 radioss接触(接触)建模

对冲击仿真使用最多的接触类型
接触类型 从节点 主面 Gap
7
点
面
是
刚度
使用
非线性罚函数 点面接触,使用任意应用场合
11
线
线
是
非线性罚函数 线线接触,使用任意应用场合
24
点或面
面
实体,否 线性定刚度罚函 适用于点面/面面接触,允许
壳,是
数
初始网格穿透/交叉
• 自动定义的点面接触面使 用变化刚度的罚函数法
Master segment (shell or solid Face)
K
Penetration into the gap
例子,点面接触
定义主面使用 Mast_id 和 定义从节点使用 Slav_ID
Master segment Part # 2
Slave nodes Part #1
#---1----|----2----|----3----|
TYPE2 TYPE3 TYPE5 TYPE6 TYPE7 T YPE10 T YPE11 T YPE14 T YPE15 T YPE19 T YPE20 T YPE23 T YPE24
: tied node to surface : surface to surface : node to surface : Special rigid body to rigid body : general node to surface : Tied node to surface after impact : edge to edge : MADYMO ellipsoidal to node : MADYMO ellipsoidal to surface : TYPE7 + TYPE11 : general contact : special self-impact contact for airbag : general node and surface to surface
Radioss在碰撞接触分析方面的应用

Time Step Control Time Step Demo with an Example Break Best Practices Hands on Three Point Bending Lunch Hands on Cantilever Beam with Bolt Pretension Hands on Full Car Break Hands on Debug
Altair Proprietary and Confidential Information
Copyright © 2008 Altair Engineering, Inc. All rights reserved.
Chapter 1: RADIOSS Introduction
Application Fields Modeling A Physical Problem Formulations Time Integration Explicit and Implicit Method
ij f (ij )
ij t t ij t ijt
•Assemble Fint , Fhrg
xn
1 2
xn
1 2
xnt
xn1
xn
xn
1 2
t
Copyright © 2008 Altair Engineering, Inc. All rights reserved.
Implicit Newmark Flow Chart
Explicit Flow Chart
ttt
Fext Fint
Fcont vi Fi mi
Time integration
Altair Proprietary and Confidential Information
hyperworks接触分析

好问题,学习了
报告 道具 TOP 8#
/viewthread.php?tid=932569&highlight=contactsurf 2011-5-7
hypermesh中定义接触contact的问题 - CA01:HyperMesh&HyperView--前后... 页码,6/8
Powered by Discuz! 7.0.0 Licensed
© 2001-2009 Comsenz Inc.
中国仿真互动 ( 沪ICP备07510919号) | 联系我们 | Archiver | WAP
GMT+8, 2011-5-7 18:36, Processed in 0.112385 second(s), 9 queries, Gzip enabled.
误
CONTACT interface
1 has Slave
node 900 that lies exactly on the
Master surface with OPENGAP option. Impossible to define pushout
direction.
请教下我少做了哪个步骤..谢
谢.
/viewthread.php?tid=932569&highlight=contactsurf 2011-5-7
hypermesh中定义接触contact的问题 - CA01:HyperMesh&HyperView--前后... 页码,2/8
hypermesh中定义接触contact的问题 - CA01:HyperMesh&HyperView--前后... 页码,4/8
ansys非线性接触分析中的接触行为

ansys非线性接触分析中的接触行为ansys非线性接触分析中接触行为接触是状态改变非线性,经典ANSYS版本中共提供了7种接触行为,每一种都有其特点及相应的应用范围,在选用的时候应该谨慎。
(1)标准接触行为(standard)该接触行为包括了法向接触闭合和分开行为,在该接触模式中既考虑粘着摩擦同时也考虑了滑动摩擦。
如图上,AB与BC本来是分开的,中间通过B点连接,当在A点施加力F,AB慢慢贴近BC,最终靠在一起。
但F撤销后,AB在恢复力的作用下慢慢回复到初始分开状态。
标准接触行为包括了分开状态→闭合状态→分开状态。
当AB与BC靠在一起时,既存在正压力,同时还有沿BC圆弧切线方向的摩擦力。
(2)粗糙接触行为(rough)该接触行为包括了法向接触闭合和分开行为,但滑动行为在此是不会发生的。
原因是所有参与接触的表面都被假定为非常粗糙,以致于可以认为摩擦力无穷大而不能够产生相对滑动。
在这种接触行为中,接触的两个物体或部件之间,除了存在正压力外,还有切向摩擦力,但是接触部分之间不可以产生相对滑动。
(3)绑定接触行为(bonded)是指一旦接触关系建立,那么目标面及接触面就被假定为粘结在一起(不可以分开)。
(4)绑定接触行为(始终)(bonded(always))任何初始时在许可接触容差范围内探测到的接触点或者是那些即将进入接触的点在后续的分析中将被绑定在一起。
这种接触行为的典型应用,如在组装分析中将两种不同网络的组件“加”在一起。
线性静态分析也可以用该种接触行为来解决,虽然由于有接触单元的存在,分析中将会提示为非线性分析,但往往只要一步迭代就完成了。
(5)绑定接触行为(初始接触)(bonded(initial))绑定仅发生在初始状态下就接触的面上,初始状态下没有接触的部分将继续保持分开。
典型的例子是通过焊接连接在一起的两个物体,焊接部分始终保持连接,没有焊接的部分保持分离状态。
(6)不分开型(no separation)一旦接触关系建立,目标面及接触面便被约束在一起了,但还是允许接触面之间有滑动。
MSC系列软件介绍

基于有限元技术的 CAE 软件,无论在数量、规模上,还是应用范围上都处于主要 地位,作为 20 世纪中期兴起的技术手段,有限元技术随着计算机技术的迅猛发展,得 到了飞速发展和广泛的应用。
目前,CAE 软件主要应用于汽车、航空、航天、电子、土木工程、通用机械、兵器、 核能、医疗器械、铁道、石油和化工等行业。在国外,汽车行业是有限元软件的主要应 用行业,其所涉及的专业领域相当广泛,并且应用历史长、应用成熟度高。
2.1 MSC 公司介绍
MSC.Software Corporation 创建于 19863 年,总部设在美国洛杉矶,是享誉全球最 大的工程校验、有限元分析和计算机仿真预测应用软件(CAE,computer aided computer) 供 应 商 , 其 最 早 最 著 名 的 产 品 是 大 型 通 用 结 构 有 限 元 分 析 软 件 MSC.Nastran 。 MSC.Nastran 软件始终作为美国联邦航空管理局(FAA)飞行器适航证领取的唯一验证软 件。
-3-
所有这些都是可以实现的。 VPD 建立在企业现有的投资基础上,包括设计和分析/仿真软件、计算机硬件、样
Abaqus_质量缩放及其结果准确性的评估

Modeling and Simulation 建模与仿真, 2023, 12(2), 1041-1047 Published Online March 2023 in Hans. https:///journal/mos https:///10.12677/mos.2023.122098Abaqus 质量缩放及其结果准确性的评估赵益兴上海工程技术大学机械与汽车工程学院,上海收稿日期:2023年2月7日;录用日期:2023年3月8日;发布日期:2023年3月15日摘要 Abaqus 质量缩放的目的是为了减少或者节约计算时间。
主要表现在,在显式动力学的求解过程中,实际中的一秒或者零点几秒就需要花费计算机更长的时间进行运算,无法满足实际工程或研究的需要。
有时可以适当扩大研究对象的质量,来达到减少计算机运行计算的时间。
通过简单例子来揭示质量缩放意义和对结果正确性进行评估。
最后对质量缩放的范围和相关设置进行介绍。
关键词Abaqus ,质量缩放,时间,评估Abaqus Mass Scaling and Evaluation of the Accuracy of Its ResultsYixing ZhaoSchool of Mechanical and Automotive Engineering, Shanghai University of Engineering Science, Shanghai Received: Feb. 7th , 2023; accepted: Mar. 8th , 2023; published: Mar. 15th , 2023AbstractThe purpose of Abaqus mass scaling is to reduce or save calculation time. The main performance is that in the process of solving the explicit dynamics, it takes longer time for the computer to op-erate in one second or a few seconds in practice, which cannot meet the needs of actual engineer-ing or research. Sometimes the quality of the research object can be appropriately expanded to reduce the time of computer operation and calculation. This paper presents the significance of quality scaling and evaluates the correctness of the results through simple examples. The range of mass scaling and related settings is also introduced.赵益兴KeywordsAbaqus, Mass Scaling, Time, EvaluationCopyright © 2023 by author(s) and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License (CC BY 4.0)./licenses/by/4.0/1. 引言在实际工程应用中,常常会遇到接触、碰撞等高度非线性情形,而工业仿真软件Abaqus中的显式动力学的优势正是解决此类问题。
Radioss分析类型汇总视频讲解

Radioss分析类型视频讲解Hypermesh作为一款高效的前处理软件,其强大而高效的网格划分技术正被越来越多的工程师所青睐。
Radioss 是著名的多学科分析软件包,具有显式、隐式求解器,拥有拉格朗日、欧拉、ALE、SPH等算法,及有限元、有限体积、边界元等数值处理技术;可用于解决结构的线性问题,如强度、刚度、模态、振动噪声平顺性(NVH)等问题,以及几何非线性(大位移、大转动和大应变)、材料非线性和接触非线性等非线性问题;也可以用于计算动载荷、静载荷下的固体结构、流体、流固耦合等问题,以及瞬态响应的被动安全/结构适撞性、爆炸/冲击、跌落、陆地/水上迫降,高频振动、超高速碰撞等非线性动力冲击问题,提供无与伦比的先进技术以帮助用户提高产品性能。
Radioss 的线性分析功能极其强大,是一个先进、快速、精确和稳健的有限元求解器。
通过最新的单元格式,稀疏矩阵技术,模态综合法和Lanczos求解技术,Radioss 能够快速而精确地求解大规模的有限元模型。
通过无缝集成在HyperWorks软件中,Radioss 为工程师提供了梦寐以求的灵活、快速和先进的功能,成为有限元线性分析的首选工具,包括如下求解能力:线性静力分析、模态分析、屈曲分析、频响分析、随机响应分析、疲劳分析、热力分析、复合材料分析等。
本视频正是为了各位看客能快速入门使用Radioss进行计算分析录制的,每一种分析类型都选取了典型模型进行了演示,讲解细致透彻,视频画质清晰,语音响亮,特别适合入门学习以及进阶提升。
本讲解包含如下内容:Altair官方有声培训视频,时长一小时五十八分钟分钟。
大小:125Mb视频专用播放器webex player安装文件。
视频详细介绍了:1、Radioss静力计算中采用gap单元定义接触的应用;2、Radioss静力计算中采用接触对contact的应用;3、Radioss控制卡片control card中控制输出文件的选项;4、Radioss计算模态分析(自由模态和约束模态);5、Radioss下计算线性屈曲分析;6、Radioss下惯性释放分析;7、Radioss下直接频率响应分析设置;视频链接地址:/item.htm?id=176********有疑问请咨询QQ:542436072,枫雪天一接触模拟方式——GAP单元接触模拟方式——GAP单元非线性迭代卡片NLPARM设置GAP单元模拟接触结果云图接触模拟方法——接触对设置接触模拟方法——接触对设置接触模拟结果对比分析控制卡片CONTROL CARD选择性输出设置模态计算模型——电机及其机座模态计算结果显示补充说明——Radioss下单位匹配线性屈曲计算惯性释放计算——简化卫星模型直接频率响应计算模型——悬臂梁直接频率响应计算理论讲解。
RADIOSS冲击运动常用接触及设置介绍

RADIOSS冲击运动接触建模接触建模的典型应用定义接触面是为了模拟不同部件接触时,产生的相互作用。
用于模拟接触的罚函数法罚函数法将从节点和主面之间的作用等效为弹簧,通过罚函数方程产生罚函数力。
在节点与投影点之间创建弹簧单元,根据主面和从节点的材料和属性自动计算弹簧的刚度。
接触类型节点—面面—面自接触节点—面接触概述在RADIOSS 中有大量的接触类型可供选择FLUID / FSI applicationsTYPE1 TYPE9 TYPE12 TYPE18: tied fluid to structure: ALE Lagrange with void opening and free surface contact: Fluid to fluid contact: Fluid structure contactManufacturing contactsTYPE8 TYPE21: drawbeads model: node to rigid surfaceCrash structure applicationsTYPE2TYPE3TYPE5TYPE6TYPE7TYPE10TYPE11TYPE14TYPE15TYPE19TYPE20TYPE23TYPE24: tied node to surface: surface to surface: node to surface: Special rigid body to rigid body: general node to surface: Tied node to surface after impact: edge to edge: MADYMO ellipsoidal to node: MADYMO ellipsoidal to surface: TYPE7 +TYPE11: general contact: special self-impact contact for airbag: general node and surface to surfaceSpecial Crash structure applications/LAGMUL/TYPE16/LAGMUL/TYPE17: Node to quadratic brick: quadratic brick to quadratic brick/INTER/HERTZ/TYPE17 : quadratic brick to quadratic brick with hertz pressure contact…对冲击仿真使用最多的接触类型接触类型从节点主面Gap刚度使用7点面是非线性罚函数点面接触,使用任意应用场合11线线是非线性罚函数线线接触,使用任意应用场合24点或面面实体,否壳,是线性定刚度罚函数适用于点面/面面接触,允许初始网格穿透/交叉•自动定义的点面接触面使用变化刚度的罚函数法•主面和从节点之间的空间定义为Gap•接触Gap可以是恒定的也可以是变化•阻力用于保证从节点位于接触Gap范围之外•接触时间步长能影响计算时间和稳定Slave nodesMaster segment(shell or solid Face)Contact gapContact segmentShell thicknessResistiveForcePenetrationinto the gapKgap例子,点面接触定义主面使用Mast_id 和定义从节点使用Slav_IDSlave nodes Part #1Master segmentPart #2#---1----|----2----|----3----|/INTER/TYPE7/2Type 7 node to surface #Slav_id Mast_id Istf1020................/GRNOD/PART/10Slave nodes of part 11/SURF/PART/EXT/20Master surface of part 22创建两个对称的接触例子,面面接触+#---1----|----2----|----3----|/INTER/TYPE7/2Type 7 node to surface#Slav_id Mast_id Istf 1020...... ........../GRNOD/PART/10Slave nodes of part11/SURF/PART/EXT/20Master surface of part22#---1----|----2----|----3----|/INTER/TYPE7/2Type 7 node to surface#Slav_id Mast_id Istf3040...... ........../GRNOD/PART/30Slave nodes of part22/SURF/PART/EXT/40Master surface of part11Slave nodesPart#2MastersegmentPart#1Master segmentPart #2 Slave nodesPart#1例子, 自接触SELF CONTACT对于自接触的选择为:part自身的从节点和的主面。
金属接触非线性引起的无源互调效应的数值分析

金属接触非线性引起的无源互调效应的数值分析毛煜茹;刘莹;谢拥军;程志华;王新波;李韵【摘要】针对时域物理光学方法求解无源互调问题的局限性,提出基于“场路结合”的数值分析方法,建立非线性金属结的等效电路模型并将其引入到全波方法中,通过全波方法与等效电路相结合的方式来分析系统的无源互调干扰.计算实例表明,该方法可以精确地计算由于金属接触非线性引起的无源互调及其作为二次辐射源所激发的电磁场,避免了时域方法在电大尺寸问题上计算的积累误差以及只能对处于远场区的简单微波结构的无源互调效应进行求解的问题.本文提出的场路结合数值分析方法为微波部件无源互调效应的研究提供了新的思路与方法.【期刊名称】《电子学报》【年(卷),期】2015(043)006【总页数】5页(P1174-1178)【关键词】场路结合;等效电路模型;电磁数值分析;无源互调【作者】毛煜茹;刘莹;谢拥军;程志华;王新波;李韵【作者单位】北京航空航天大学电子信息工程学院,北京100191;西安电子科技大学天线与微波技术重点实验室,陕西西安710071;北京航空航天大学电子信息工程学院,北京100191;北京航空航天大学电子信息工程学院,北京100191;中国空间技术研究院西安分院空间微波技术国家重点实验室,陕西西安710100;中国空间技术研究院西安分院空间微波技术国家重点实验室,陕西西安710100【正文语种】中文【中图分类】TN003两个或两个以上的载波信号经过具有非线性响应的无源器件时,会产生无源互调(Passive InterModulation,PIM)现象.无源互调一旦产生,会造成严重后果.数值方法是分析无源互调对通信系统影响的有效手段,但是数值仿真技术的困难限制了无源互调的研究进展.无源互调产生的机理可以分为两类:接触非线性和材料非线性.接触非线性表示具有非线性电流与电压特性的接触,如铆钉连接处的松动、氧化和腐蚀等;材料非线性指具有固有非线性电特性的材料[1~3].微波部件的表面一般会进行电镀处理,粗糙的表面和不合适的镀层厚度将引起非线性,产生无源互调,可以用“趋肤效应”来解释.生产实践证明,金属表面镀层的材料及厚度对PIM的影响非常大,而趋肤深度决定了导电镀层的厚度,当趋肤深度大于镀层厚度时,PIM电平将会增大,只要镀层厚度大于趋肤深度,其PIM特性与镀层金属相近[1].对于非线性材料,由于镀层的面积比较大,电磁场在非线性材料中的分布不均匀,对材料非线性的研究需采用特殊的等效方法,本文暂不作研究.金属接触非线性引起的无源互调的散射机理及其数值计算由Pelosi课题组于1999年提出[4],将时域物理光学方法(Time-Domain Physical Optics,TD-PO)的应用[5]推广于非线性环境,给出金属相互接触时产生无源互调的计算模型,在平面波照射下求解金属结产生的空间散射场,并对该散射场进行傅里叶变换求得各阶无源互调,实现对金属接触产生的无源互调的预测[6,7].之后该课题组使用基因算法对模型中的非线性部分做了优化处理[8].但是,TD-PO方法中的入射波被近似为平面波,仅适用于辐射源位于散射目标远区的情形,无法计算实际复杂微波部件产生的无源互调,并且在电大环境下时域算法所带来的误差积累是不可忽视的.电磁频域方法可以准确地求解物质结构的电磁响应,但该方法不适用于求解非线性问题;而等效电路法是分析具有非线性响应的无源部件的有效手断,白春江[9]对行波管的非线性问题进行了研究,利用等效电路模型推导了行波管的非线性模型,该模型能够准确地计算行波管的互调分量.姚常飞[10]、郝建红[11]将场路结合的分析方法用于电子设备的分析,该方法具有计算过程简单且计算精度高的优势.本文提出了场路结合的数值分析方法用于微波部件无源互调的计算,在金属接触的非线性区域采用等效电路模型来模拟,而激发无源互调的外加场和非线性区域散射的无源互调波的计算均采用准确的频域计算方法.这种场路结合的数值分析方法既使得频域计算方法可以用于非线性问题的求解,同时解决了TD-PO方法计算累积误差较大这一难题.2.1 金属接触非线性引起的无源互调金属相互接触产生的非线性是最常见的PIM根源.当暴露在空气中的两个导体接触时,由于表面粗糙度的影响以及氧化作用会在导体接触表面产生微米量级的粗糙度和纳米量级的介质层.其非线性接触面及电路等效模型如图1所示,金属与金属的连接处形成金属接触①、接触面之间有金属膜氧化物形成“金属-氧化物-金属结”结构②、金属面之间夹有绝缘介质③、微小空气间隙④、大的空气间隙⑤[1].①和②形成电流的主要通道,金属-氧化物-金属连接处②属于半导体行为,是非线性的.当电流流过阻抗Z,产生间隙电压,如图1(b)所示,可能激活任何一个半导体产生无源互调效应,具体如天馈系统的波导连接、反射面天线中的铆钉铆接等.天线和微波接头处的金属结的非线性效应可用式(1)所示的电路模型来等效[4,6], 式中的h0和eT为待定的系数,决定了非线性金属结的材料与非线性程度;hx为沿非线性金属结的磁场强度,ey为垂直于非线性金属结的电场强度.假设两个频率分别为f1和f2的大功率电磁波在金属结处产生的激励如下:式中V是合成信号,A1和A2是两个激励信号的幅度.采用电路分析方法对非线性金属结的等效电路模型式(1)展开成n阶幂级数形式[12,13]:式中,hx是输出信号,C1,C2,C3是依赖于非线性金属结特性的系数,由式(1)中的h0、eT决定.将式(2)代入式(3)整理后,就可得到hx的非线性频域表达式(4),由激励信号、谐波以及互调信号组成.式中θ1=2πf1t,θ2=2πf2t.互调产物的频率组成由方程fIM=mf1+nf2决定,m、n均为整数,定义为互调产物的阶数.至此,可得到在两个大功率辐射场作用下非线性金属结产生的无源互调.由于双载波激励下的金属结引起的PIM是在m,n的各种取值下的各阶干扰,但是在收发系统中,靠近载波的3阶和5阶互调产物最有可能落入工作频带,对系统的危害最严重,因此需要重点关注.对于一些高灵敏度的接收机,在PIM达到-140dBm时,仍然会对接收机产生较大影响,淹没实际有用信号.2.2 无源互调的场路结合计算方法基于加权残数法的电磁频域方法,在复杂结构的电磁响应的计算过程中可以利用数值残差的设定来控制计算结果的准确性,因而被广泛应用于电磁问题的数值求解.本文采用准确性较高的电磁频域方法与等效电路方法相结合的方式对非线性金属结在大功率波照射下产生的无源互调进行分析,场路结合数值分析方法计算无源互调的具体流程如图2所示.非线性金属结在两个大功率波的照射下,在金属结表面会产生切向电场E1(x,y,z)和E2(x,y,z),采用频域数值方法可以准确地求解金属结表面的切向电场.切向电场作用于非线性金属结的等效电路时将会产生电压降V1和V2,由式(5)决定:将V1和V2代入(4)式,可得到各阶无源互调产物的功率电平,如靠近载波处的3阶差频(2f1-f2)频率处产生的无源互调产物如式(6):在实际计算各阶无源互调的功率电平时,由于每一次谐波都会对比它阶次低的谐波产生影响,因此需要对式(1)的幂级数展开式进行高次幂截断,最高次幂的取值应大于需要评估的PIM的最高次幂.非线性金属结在无源互调产物的激励下,会向外辐射电磁波,并在接收机端口处引起干扰电平.无源互调作为激励源所激发的空间二次场的求解亦是采用电磁频域方法,对空间二次场在接收机端口处所激励的电场进行面积分即可获得无源互调对系统产生的干扰电平.这样就可以用准确的电磁频域方法结合等效电路方法评估微波部件的无源互调效应.3.1 金属结无源互调的计算考虑双载波频率情况,建立图3所示的平板验证模型,使用喇叭天线作为馈源,由金属平板形成的单个金属结距离喇叭天线0.3m,金属结部分的等效电路模型采用式(1).假设频率分别为f1=11GHz,f2=12.6GHz的两个电磁波垂直照射到平板上,计算分析该金属结的3阶、5阶互调干扰.使用图2所示的场路结合数值分析方法对图3所示的验证模型进行计算,首先对图3的单个金属结以及喇叭馈源进行电磁建模,采用频域计算方法求得金属结的表面电场,并在非线性区域对表面电场进行积分求得非线性区域的电压降,代入式(4)求得3阶和5阶互调产物的功率电平.再次利用频域计算方法求解3阶和5阶互调产物作为激励源所激发的空间二次辐射场,在接收端口处对电场进行积分得到的无源互调功率结果如表1所示,与测量数据[6]比较可以看出,虽然存在一定的误差,但仍可认为预测结果与测量数据基本符合,初步验证了本文场路结合数值分析方法的准确性.3.2 反射面天线无源互调的计算对处于喇叭馈源近场区域的反射面天线的无源互调产物进行计算.计算模型如图4所示,每个金属结的非线性接触区域的宽度为2mm,取f1=8GHz,f2=7.5GHz,反射面直径为3m,约100个波长的电大尺寸,对馈源以及反射面天线进行有限元网格剖分以便进行数值计算,剖分网格数为50万左右.在近场区的反射面天线的入射波并非平面波,如使用TD-PO算法对该电大尺寸模型进行计算,则会产生不可接受的积累误差,导致计算结果出现错误.采用本文方法对反射面上由于接触非线性引起的无源互调现象的40个金属结在馈源接收端口处产生的互调功率进行计算,金属结的分布及数量与参考文献一致,选取时重点考虑产生互调干扰较大的区域,即在喇叭馈源照射下反射面表面电流较强的区域.在等幅载波条件下,馈源接收端的3阶(2f1-f2)频点和5阶(3f1-2f2)频点处的互调功率随输入功率的变化情况如图5所示.在PIM达到-140dBm时,对于高灵敏度的接收机仍会产生较大影响,例如在卫星通信系统中,要求接收机端口处的PIM值小于-140dBm.从图5可以看出,3阶互调和5阶互调具有相似的变化规律,当输入功率增加时,会导致各个金属结的压降增大,从而导致金属结作为二次辐射源的激励功率上升,向空间辐射的二次场增加.即在接收机端口处所激发的电场会随输入功率的增加而增加,导致被接收机接收到的PIM干扰电平也随之增加,淹没有用信号.采用场路结合数值分析方法计算非线性金属结无源互调产物的误差来源主要是对金属结的等效电路模型进行幂级数展开时,由于高次幂的截断会带来一定的误差;且本文在计算互调电平时只考虑了两个载波所产生的无源互调,忽略了载波与互调产物、互调产物之间产生的二次无源互调的影响,因此在计算无源互调时也会出现一定的误差.基于具有非线性效应的金属结的导电特性,本文提出的场路结合数值分析方法既有效地解决了TD-PO方法在求解电大尺寸微波部件产生的无源互调时只适用于辐射源位于散射目标的散射远区的情形以及计算累积误差较大这一难题,又克服了频域计算方法无法求解非线性物质结构的电磁响应的问题.结合近场条件下反射面天线无源互调计算的数值实例验证了本文方法的有效性,为数值计算无源互调这类非线性问题提供了新的计算方法.本文通过建立基于无源互调分析的金属接触非线性模型,可应用于更为复杂的微波部件,如网状反射面天线的无源互调的计算,研究不同的丝网编织方式和网孔大小对无源互调的影响.通过本文场路结合数值分析方法构建无源互调有效预测的分析手段,并可为无源互调的抑制方案进行数值验证且提出其重点关注区域,建立无源互调预测及验证平台.毛煜茹女,1984年生于山西晋城.现为北京航空航天大学电子科学与技术博士研究生.研究方向为计算电磁学、天线工程等.E-mail:******************刘莹女,1980年生于湖北武汉.现为西安电子科技大学博士研究生.研究方向为电磁兼容、天线与微波技术及电磁场数值算法.谢拥军男,1968年生于陕西西安.现为北京航空航天大学教授、博士生导师.研究方向为计算电磁学、太赫兹技术与天线工程等.【相关文献】[1]罗一锋.舰船通信系统的无源互调研究[J].现代电子学报,2010,23:39-44. Luo Yifeng.Approach to passive intermodulation in ship communication system[J].Modern Electronics Technique,2010,23:39-44.(in Chinese)[2]叶鸣,贺永宁,王新波,崔万照.金属波导连接的无源互调非线性物理机制和计算方法[J].西安交通大学学报,2011,45(2):82-86. Ye Ming,He Yongning,Wang Xinbo,Cui Wanzhao.Nonlinear physical mechanism calculation method of passive intermodulation at metal waveguide connection[J].Journal of Xi’an Jiaotong University,2011,45(2):82-86.(in Chinese)[3]王海宁,梁建刚,王积勤,张晨新.高功率微波条件下的无源互调问题综述[J].微波学报,2005,21(Supplement):1-6. Wang Haining,Liang Jiangang,Wang Jiqin,Zhang Chenxin.Review of passive intermodulation in HPM condition[J].Journal of Microwaves,2005,21(Supplement):1-6.(in Chinese)[4]P Bolli,A Naldini,G Pelosi,C Pochini,S Selleri.A time domain physical optics approach to passive intermodulation scattering problems[A].IEEE International Symposium on Antennas and Propagation Digest[C].Florida:IEEE,1999.2018-2021.[5]E Y Sun,W V T Rusch.Time-domain physical-optics[J].IEEE Transactions on Antennas and Propagation,1994,42(1):9-15.[6]S Selleri,P Bolli,G Pelosi.A time domain physical optics heuristic approach to passiveintermodulation scattering[J].IEEE Transactions on ElectromagneticCompatibility,2001,43(2):203-209.[7]P Bolli,P Pelacchi,G Pelosi,S Selleri.Passive intermodulation on large reflector antennas[J].IEEE Antenna’s and Propagation Magazine,2002,44(5):13-20.[8]S Selleri,P Bolli,G Pelosi.Genetic algorithms for the determination of the nonlinearity model coefficients in passive intermodulation scattering[J].IEEE Transactions on Electromagnetic Compatibility,2004,46(2):309-311.[9]白春江,李建清,胡玉禄,杨中海,李斌.利用等效电路模型计算耦合腔行波管注-波互作用[J].物理学报,2012,61(17):178401-1-178401-6. Bai Chunjiang,Li Jianqing,Hu Yulu,Yang Zhonghai,Li Bin.Calculation of beam-wave interaction of coupled-cavity TWT using equivalent circuit model[J].Acta Physica Sinica,2012,61(17):178401-1-178401-6.(in Chinese)[10]姚常飞,周明,罗运生,许从海,寇亚男,陈以钢.基于肖特基势垒二极管的太赫兹固态倍频源和检测器研制[J].电子学报,2013,41(3):438-443. Yao Changfei,Zhou Ming,Luo Yunsheng,Xu Conghai,Kou Yanan,Chen Yigang.Development of terahertz frequency solid state multiply sources and sensors with Schottky barrier diodes[J].Acta Electronica Sinica,2013,41(3):438-443.(in Chinese)[11]郝建红,孙娜燕,高璞,范杰青.场路结合法分析设备电磁辐射照应[J].电波科学学报,2012,27(2):359-364. Hao Jianhong,Sun Nayan,Gao Pu,Fan Jieqing.Analysis of the electromagnetic radiation effect of the equipment by field-circuit method[J].Chinese Journal of Radio Science,2012,27(2):359-364.(in Chinese)[12]张国艳,廖怀林,黄如,Mansun CHAN,张兴,王阳元.适用于数模混合集成的SOI MOSFET 的失真分析[J].电子学报,2002,30(2):153-156.Zhang Guoyan,Liao Huailin,Huang Ru,Mansun Chan,Zhang Xing,Wang Yangyuan.The distortion analysis of SOI MOSFET for mixed mode integrated circuit application[J].Acta Electronica Sinica,2002,30(2):153-156.(in Chinese)[13]张世全,葛德彪.基于傅里叶级数法的互调产物一般特性分析[J].电波科学学报,2005,20(2):265-268. Zhang Shiquan,Ge Debiao.Analysis of general nature of intermodulation products based on Fourier series method[J].Chinese Journal of Radio Science,2005,20(2):265-268.(in Chinese)。
工程有限元方法接触非线性讲义PPT(30页)

反向,使得解的过程不稳定
4 接触问题的有限元方程
• 4.1 接触界面的离散处理
– 接触块
接触块(或线):单元上处于接触面上的面(或边)
二维问题的主被动接触线
– 接触点对
罚函数法
惩罚因子
z* x, y. 2x2 2xy y2 18x 6y x y2
z* 4x 2 y 18 2 x y 0
x
z* 2x 2 y 6 2 x y 0
y
x 12 ; y 12 15 / 11/
y 12
3 接触问题的求解方案
原问题中不包含 接触约束条件的
3. 利用接触面上和上述等式约束所对应的动力学或运动学的不等式 约束条件作为校核条件对解的结果进行检查 ➢ 若不违反接触条件,则完成本步求解并转入下一增量步 ➢ 否则转入1再次进行搜寻和迭代求解,直到每一点的解都满足 校核条件
3 接触问题的求解方案
• 3.2 接触面的定解条件和校核条件
接触状态
定解条件
有限元方法与
• 接触非线性概述 • 接触界面条件 • 接触问题的求解方案 • 接触问题的有限元方程 • 有限元方程的求解方法
1 接触非线性概述
接触非线性的来源:由随时间发生变化的接触状态引起,主要 包括
(1) 接触界面的区域大小和相互位置以及接触状态不仅事先未知,且都随时 间发生变化,需要在求解过程中确定
2 接触界面条件
• 2.3 切向接触条件
– (1) 无摩擦模型
若两个物体的接触面是绝对光滑的,或 者相互间的摩擦可以忽略,则采用无摩 擦模型
非线性-接触分析

接触分析接触问题是一种高度非线性行为,需要较大的计算资源,为了进行实为有效的计算,理解问题的特性和建立合理的模型是很重要的。
接触问题存在两个较大的难点:其一,在你求解问题之前,你不知道接触区域,表面之间是接触或分开是未知的,突然变化的,这随载荷、材料、边界条件和其它因素而定;其二,大多的接触问题需要计算摩擦,有几种摩擦和模型供你挑选,它们都是非线性的,摩擦使问题的收敛性变得困难。
一般的接触分类接触问题分为两种基本类型:刚体─柔体的接触,半柔体─柔体的接触,在刚体─柔体的接触问题中,接触面的一个或多个被当作刚体,(与它接触的变形体相比,有大得多的刚度),一般情况下,一种软材料和一种硬材料接触时,问题可以被假定为刚体─柔体的接触,许多金属成形问题归为此类接触,另一类,柔体─柔体的接触,是一种更普遍的类型,在这种情况下,两个接触体都是变形体(有近似的刚度)。
ANSYS接触能力ANSYS支持三种接触方式:点─点,点─面,平面─面,每种接触方式使用的接触单元适用于某类问题。
为了给接触问题建模,首先必须认识到模型中的哪些部分可能会相互接触,如果相互作用的其中之一是一点,模型的对立应组元是一个结点。
如果相互作用的其中之一是一个面,模型的对应组元是单元,例如梁单元,壳单元或实体单元,有限元模型通过指定的接触单元来识别可能的接触匹对,接触单元是覆盖在分析模型接触面之上的一层单元,至于ANSTS 使用的接触单元和使用它们的过程,下面分类详述。
点─点接触单元点─点接触单元主要用于模拟点─点的接触行为,为了使用点─点的接触单元,你需要预先知道接触位置,这类接触问题只能适用于接触面之间有较小相对滑动的情况(即使在几何非线性情况下)如果两个面上的结点一一对应,相对滑动又以忽略不计,两个面挠度(转动)保持小量,那么可以用点─点的接触单元来求解面─面的接触问题,过盈装配问题是一个用点─点的接触单元来模拟面─与的接触问题的典型例子。
《2024年基于ANSYS软件的接触问题分析及在工程中的应用》范文

《基于ANSYS软件的接触问题分析及在工程中的应用》篇一一、引言随着现代工程技术的快速发展,接触问题在各种工程领域中变得越来越重要。
ANSYS软件作为一款强大的工程仿真软件,其在接触问题上的分析和处理能力得到了广泛应用。
本文将介绍基于ANSYS软件的接触问题分析及在工程中的应用。
二、ANSYS软件接触问题分析1. 接触问题基本理论接触问题是一种典型的非线性问题,涉及到两个或多个物体在力、热、电等作用下的相互作用。
在ANSYS中,接触问题主要通过定义接触对、设置接触面参数、定义接触刚度等方式进行模拟。
2. ANSYS软件接触问题处理流程(1)建立模型:在ANSYS中建立涉及接触问题的物理模型。
(2)定义材料属性:设置模型中各部分的材料属性,包括弹性模量、密度、泊松比等。
(3)划分网格:对模型进行网格划分,以便更好地进行后续的数值分析和计算。
(4)定义接触对:根据实际需求,定义接触对,并设置相应的接触面参数。
(5)求解设置:设置求解器、求解参数等。
(6)结果分析:对求解结果进行分析,包括应力分布、位移变化等。
三、ANSYS软件在工程中的应用1. 机械工程领域在机械工程领域,ANSYS软件被广泛应用于分析各种机械零件的接触问题。
例如,齿轮传动中齿轮与齿轮之间的接触问题、轴承中滚动体与内外圈的接触问题等。
通过ANSYS软件的分析,可以有效地预测机械零件的应力分布、疲劳寿命等,为机械产品的设计和优化提供有力支持。
2. 土木工程领域在土木工程领域,ANSYS软件被广泛应用于分析土与结构之间的接触问题。
例如,桥梁、大坝等结构物与地基之间的相互作用、地震作用下建筑结构的动力响应等。
通过ANSYS软件的分析,可以有效地评估结构的稳定性和安全性,为土木工程的设计和施工提供有力支持。
3. 汽车工程领域在汽车工程领域,ANSYS软件被广泛应用于分析汽车零部件的接触问题。
例如,汽车发动机的缸体与缸盖之间的密封问题、汽车轮胎与地面的摩擦问题等。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
HyperWorks SolversThe purpose of this tutorial is to demonstrate how to carry out nonlinear implicit small displacement analysis in OptiStruct, involving elasto-plastic materials, contact and continuing the nonlinear solution sequence from a preceding nonlinear loadcase.Model and Loading DescriptionFigure 1 illustrates the structural model used for this tutorial: Two square solid blocks made of elasto-plastic steel material. The dimensions of the blocks and the material parameters can be obtained in the table below.In the first nonlinear subcase, pressure loading is be applied to the top solid block, the top corners of which are constrained in X and Y directions. The top solid is in contact with the bottom solid, the bottom corners of which are constrained in X, Y and Z directions. The second nonlinear subcase isto simulate the un-loading and is a continuation of the nonlinear solution sequence from the previous loading subcase.Figure 1. Model and Loading DescriptionUnits Length: mm; Time: s; Mass: Mgg; (Force: N; Stress: MPa)Top block72 mm x 72 mmBottom block100 mm x 100 mmThickness of20. mmblocksMaterial Steel, Elasto-plasticInitial density ():7.90e-9 kg/mm3Young's modulus (E):210000 MPaPoisson coefficient (ν):0.3Yield Stress (0):850.0 MPaImposed pressure1000.0 MPa, applied at the center of top blockThe following exercises are included in this tutorial:•Create elasto-plastic material•Define contact between the two blocks•Define nonlinear implicit parameters•Set up NLSTAT analysis for the 1st subcase (loading)•Set up NLSTAT analysis for the 2nd subcase (unloading)•Submit job and view resultExerciseStep 1: Import the modelunch HyperMesh.2. A User Profiles Graphic User Interface (GUI) will appear, select OptiStruct.3.Click OK. This loads the User Profile.4.Click File>Open.Note:If HyperMesh Desktop was launched, use:File > Open >Model.5.Select the file nlstat.hm located in the HyperWorks installation directory under<install_directory>/tutorials/hwsolvers/optistruct/.6.Click Open.Step 2: Create the elasto-plastic materialFirst, the stress vs plastic strain curve for the material needs to be defined.1.In the Tab area, click on the Utility menu. If not active, click View > Utility.2.Click on the TABLE Create tool.3.Click on Create/Edit Table and select the table type as TABLES1.4.Click Next.5.Click the radio button for Create New Table and enter the Name of the table as stress-strain.6.Next, populate the X and Y fields of the table, as shown below.7.Click Apply.A dialog box will appear stating: 'The load collector “stress-strain”with TABLES1 card image iscr e ated’.8.Click OK to close the dialog box.Now, the elasto-plastic material needs to be updated.9.In the toolbar, click the Material icon .10.Toggle the Update subpanel.11.Click on mats and select the material steel.12.Click update/edit.13.Click on MATS1 to define the elastic-plastic material for NLSTAT analysis.14.Click TID and select stress-strain.15.Input the values, as shown below. [TYPSTRN] of 1 signifies specifying stress (Y) vs plasticstrain (X).See material parameters for details.16.Click return twice to go to the main menu.Step 3: Define contact between the two blocksThe contact surfaces for the two blocks need to be defined.1.In the Analysis page, click on entity sets panel.2.For name=, enter top.3.Set card image to SET_ELEM.4.Toggle entity and select props.5.Click on props and select the top solid Solid1.6.Click create.7.Next, for name=, enter bottom.8.Repeat steps 3 through 6, for bottom block select the bottom solid Solid2.9.Click return to go to the main menu.Next, the interface needs to be defined.10.Click on interfaces panel with the toggle set to create.11.For name =, enter solid_contact.12.For type =, select CONTACT, then click create.13.Next, click the radio button for add to add the master and slave surfaces to the interfacedefinition, as shown below.14.Toggle master: to sets and select bottom.15.Click update.16.Toggle slave: to sets and select top.17.Click update.18.Click review to review the interface, as shown in Figure 2.Figure 2. Slave and master surface definition19.Click on card image and then click edit to specify the contact definition, as shown below.20.Click return twice to return to the main menu.Step 4: Define nonlinear implicit parameters1.In the toolbar, click on Load Collectors icon .2.Toggle create and for loadcol name =, enter nlparm.3.Toggle card image =, and select NLPARM.4.Click create/edit and specify the defaults, as shown below.See nonlinear implicit parameters for details.5.Click return twice to go to the main menu.Step 5: Create the 1st nonlinear (loading) subcase1.Click on the loadsteps panel.2.For name =, enter loading.3.Toggle type: to non-linear quasi-static.4.Check the box preceding SPC. An entry field appears to the right of SPC.5.Click on the entry field and select the SPC load collector.6.Check the box preceding LOAD. An entry field appears to the right of LOAD.7.Click on the entry field and select the pressure load collector.8.Check the box preceding NLPARM. An entry field appears to the right of NLPARM.9.Click on the entry field and select the nlparm load collector.10.Click create.Step 6: Create the 2nd nonlinear (unloading) subcase1.For name =, enter unload.2.Uncheck the box preceding LOAD.3.Keep the boxes preceding SPC and NLPARM checked, as shown below.4.Click create.5.Click edit and select CNTNLSUB to continue the solution sequence from the previous loadingsubcase.6.Click return twice to go to the main menu.Step 7: Define Output Control parameters1.From the Analysis page, select control cards.2.Click on GLOBAL_OUTPUT_REQUEST.3.Under CONTF, DISPLACEMENT, STRAIN and STRESS, set Option to Yes.4.Under STRAIN, set TYPE(1) to PLASTIC.5.Click return twice to go to the main menu.Step 8: Submit the job1.From the Analysis page, select OptiStruct.2.Click save as and select the directory where you want to write the file.3.For File Name, enter nlstat_complete.fem.4.Click Save. The file extension .fem is necessary for OptiStruct to recognize it as an input file.5.Set the export options toggle to all.6.Set the run options toggle to analysis.7.Set the memory options toggle to memory default.8.Click OptiStruct. This launches the OptiStruct job.If the job is successful, you should see new results files in the directory where HyperMesh was invoked. The nlstat_complete.out file is where you will find error messages that will help you debug your input deck, if any errors are present.The default files that will be written to your directory are:nlstat_complete.html HTML report of the analysis, giving a summary of theproblem formulation and the analysis results.nlstat_complete.out ASCII based output file of the model check run before thesimulation begins and gives nonlinear iteration history as wellas basic information on the results of the run.nlstat_complete.stat Summary of analysis process, providing CPU information foreach step during the process.nlstat_complete.h3d HyperView compressed binary results file.Step 9: View analysis resultsing HyperView, plot the Displacement, the von Mises stress, plastic strains and contactpressure contours at the end of the 1st (loading) step.Figure 3. Contour of Displacements in blocks subject to LoadingFigure 4. Contour of von Mises Stress in blocks subject to LoadingFigure 5. Contour of Plastic Strains in the blocks subject to LoadingFigure 6. Contour of Contact Pressure in the block interface after the 1st (loading) subcaseNext, change the subcase to the 2nd that is unloading subcase and plot the displacement contour to see the change in displacements in the blocks subject to unloading.Figure 7. Contour of Displacements in blocks subject to Unloading in 2nd subcaseSee Also:OptiStruct Tutorials。