材料非线性分析
钢筋混凝土板的非线性分析
钢筋混凝土板的非线性分析钢筋混凝土板的非线性分析钢筋混凝土板是一种常用的结构构件，在建筑和桥梁中广泛应用。
由于其在使用过程中会受到各种荷载的作用，因此需要对其进行非线性分析，以确保其安全可靠。
非线性分析是指在分析过程中考虑材料和结构的非线性特性，包括材料的本构关系、几何非线性和接触非线性等因素。
在钢筋混凝土板的非线性分析中，需要考虑以下几个方面。
1. 材料的本构关系钢筋混凝土板的材料包括混凝土和钢筋两部分，它们的本构关系是非线性的。
混凝土的本构关系可以采用双曲正切模型或Drucker-Prager 模型等进行描述，而钢筋的本构关系则可以采用弹塑性模型或Ramberg-Osgood模型等进行描述。
在进行非线性分析时，需要考虑这些材料的本构关系对结构的影响。
2. 几何非线性钢筋混凝土板在受到荷载作用后会发生变形，这种变形会导致结构的几何非线性。
几何非线性包括平面内的弯曲变形和平面外的扭转变形等。
在进行非线性分析时，需要考虑这些几何非线性因素对结构的影响。
3. 接触非线性钢筋混凝土板在使用过程中会受到多种荷载的作用，其中包括接触荷载。
接触非线性是指结构中两个或多个体之间的接触面会发生变形，从而影响结构的力学性能。
在进行非线性分析时，需要考虑接触非线性对结构的影响。
以上三个方面是钢筋混凝土板非线性分析的关键因素，下面将对其进行详细介绍。
1. 材料的本构关系混凝土的本构关系可以用双曲正切模型或Drucker-Prager模型等进行描述。
其中，双曲正切模型是一种常用的混凝土本构模型，其本构方程如下：σ = f(ε) = σc + α(ε-εc) + β(ε-εc)/(1+(ε-εc)/γ)其中，σ为混凝土的应力，ε为混凝土的应变，σc和εc分别为混凝土的极限应力和极限应变，α、β和γ为模型参数。
该模型可以较好地描述混凝土的非线性本构关系。
钢筋的本构关系可以采用弹塑性模型或Ramberg-Osgood模型等进行描述。
材料非线性
( ) P( ) f 0 0
其中： 表示载荷变化的量。 dP d d f 0 KT f0 0 d d d d 1 K T ( ) f 0 d 切线矩阵
1 1 m1 m KT ( m ) f0m KT ( m )f m
一、材料弹塑性行为的描述
弹塑性材料进入塑性的特征是当载荷卸去后 存在不可恢复的永久变形，因而在涉及卸载的情 况下，应力和应变之间不再存在一一对应的关系， 这是区别于非线性弹性的基本属性。
11
单调加载 对于大多数材料存在屈服应力，应力低于屈服 应力时，材料为弹性，而当应力超出屈服应力时， 材料进入 弹塑性状态。 当应力达到屈服应力后，应力不再增加，而材料 变形可以继续增加—理想弹塑性材料。
第六章 材料非线性问题的有限元法
1
第一节
引言
线弹性力学基本方程的特点： 几何方程的位移和应变的关系是线性的； 物理方程的应力和应变的关系是线性的； 建立于变形前的平衡方程也是线性的。 几何非线性问题 结构的变形使体系的受力状态发生显著变化， 以致于不能用变形前的平衡方程分析，且位移和应 变的关系不是线性的。
K ( ) f 0
增量法 载荷分为若干步： f 0 , f1 , f 2 , f 3 位移分成若干步： 0 , 1 , 2 , 3 每两步之间增长量为增量。 增量解法的一般做法是： 假设第m步的载荷 f m 和位移 m ； 让载荷增加 f m1 ( f m f ) ，再求解 m1( m )。 如果每一步的增量 f 足够小，解的收敛性 可以得到保证
9
NmR-N方法求解非线性方程组时，收敛速度 较慢，特别是对于结构分析时载荷趋近极限载荷或突 然变软的情况下，收敛速度会很慢。为了加速收敛， 可以采用一些方法，比较常用和有效的是Aitken法。 该方法每隔一次迭代进行一次加速。
光学材料中的非线性光学特性分析
光学材料中的非线性光学特性分析光学材料是指能够对光进行控制、调节以及产生新的光学效应的材料。
非线性光学特性是光学材料中一种重要的现象，其研究在光通信、激光技术、光信息处理等领域具有广泛的应用价值。
本文将对光学材料中的非线性光学特性进行分析，探讨其机理以及应用前景。
1. 非线性光学特性简介非线性光学特性是指当光与光学材料相互作用时，产生的光学效应与入射光强度不呈线性关系的现象。
与线性光学特性不同，非线性光学特性由于其强度依赖关系的非线性性质，使得光学材料在应用中具有更加丰富的功能和效果。
常见的非线性光学效应包括二次谐波发生、和频与差频发生、自聚焦、自相位调制等。
2. 非线性光学效应的机理非线性光学效应的产生是由于光照射到光学材料中的原子或分子后，其能级结构发生变化并引发非线性相互作用。
比如，二次谐波发生是由于材料的非线性极化率产生了非线性响应，将入射的光分解为频率为二倍的新光。
自聚焦效应是由于材料的光折射率与光强度的关系非线性，使得光束在传播过程中自动聚焦。
3. 光学材料中的非线性光学特性研究方法为了研究和应用光学材料中的非线性光学特性，科学家们发展了多种实验方法。
其中，著名的方法包括Z-scan技术、功率扭曲、相位匹配等。
Z-scan技术可测量材料的非线性吸收和折射率，并通过测量传播动力学过程来分析非线性效应。
功率扭曲实验通过改变光束强度来研究材料的非线性响应。
相位匹配为材料中的非线性效应提供了最佳的相位条件，以增强非线性光学效应。
4. 非线性光学特性在光通信中的应用非线性光学特性在光通信中具有重要的应用价值。
比如，光纤通信中信号调制和光时钟的生成都离不开非线性光学效应。
非线性光学特性还可用于光通信中的光放大器、光开关和光限幅器等器件的设计和制造。
利用非线性光学特性，还可以实现光通信中的非线性光调制和光波混频等功能。
5. 非线性光学特性在激光技术中的应用非线性光学特性在激光技术中有着广泛的应用。
非线性结构的变形与稳定性分析
非线性结构的变形与稳定性分析随着科技的进步和工程领域的发展，越来越多的非线性结构被广泛应用于各种工程项目中。
非线性结构的变形与稳定性分析成为了一个重要的研究领域。
本文将从非线性结构的变形分析和稳定性分析两个方面进行探讨。
一、非线性结构的变形分析非线性结构的变形分析是指在施加荷载作用下，结构的变形情况以及在变形过程中的力学特性如何变化的研究。
非线性结构的变形分析需要考虑以下几个因素：1. 材料非线性材料的非线性是非线性结构变形的主要原因之一。
传统的线弹性理论无法准确描述结构在大变形情况下的行为。
因此，非线性材料力学性质的研究和建模非常重要。
2. 几何非线性几何非线性是指在变形过程中，结构的形状和尺寸发生变化，相邻杆件之间的夹角和边长发生变化。
几何非线性的存在使得结构的变形情况更为复杂。
3. 边界条件非线性边界条件的非线性是指结构的边界条件随着变形而变化。
例如，施加在结构上的约束力随着变形而变化，从而影响结构的变形情况。
4. 辅助载荷非线性辅助载荷的非线性是指在结构变形过程中，施加在结构上的辅助力随着变形而变化。
这些辅助载荷可能来自于支撑结构的杆件或者其他零部件。
二、非线性结构的稳定性分析非线性结构的稳定性分析是指在施加荷载作用下，结构是否能够保持平衡和稳定的研究。
稳定性分析是保证结构安全性和可靠性的重要手段，需要考虑以下几个因素：1. 局部稳定性局部稳定性是指结构中的局部部分在承受荷载时是否会发生失稳。
局部失稳可能导致结构的整体性能下降，甚至引起局部的崩塌或破坏。
2. 全局稳定性全局稳定性是指整个结构在承受荷载时是否能够保持平衡和稳定。
全局失稳可能导致结构整体的倾覆、折断等严重后果。
3. 塑性转变塑性转变是非线性结构在承受荷载过程中由弹性状态向塑性状态的转变过程。
塑性转变对于结构的稳定性具有重要影响，需要进行充分的分析和设计。
4. 承载能力分析承载能力分析是指在稳定性分析的基础上，对结构的最大承载能力进行评估和计算。
工程力学中的非线性分析方法有哪些？
工程力学中的非线性分析方法有哪些？在工程力学领域，非线性问题的研究至关重要。
与线性问题相比，非线性问题更加复杂，需要采用专门的分析方法来准确描述和解决。
下面我们就来探讨一下工程力学中常见的非线性分析方法。
首先要提到的是有限元法。
这是一种非常强大且广泛应用的数值分析方法。
在处理非线性问题时，它能够有效地模拟材料的非线性行为，比如塑性、蠕变等。
通过将复杂的结构离散为有限个单元，并对每个单元进行分析，最终得到整个结构的响应。
对于几何非线性问题，如大变形、大转动等，有限元法能够通过更新坐标和刚度矩阵来准确捕捉结构的变化。
而对于材料非线性，如弹塑性问题，通过定义合适的本构关系，可以精确地模拟材料在不同应力状态下的行为。
再来看看边界元法。
它是另一种有效的数值方法，特别适用于处理无限域或半无限域问题。
在非线性分析中，边界元法可以结合迭代算法来求解非线性边界条件或非线性材料特性。
与有限元法相比，边界元法通常只需要对边界进行离散，从而降低了问题的维数，减少了计算量。
但在处理复杂的非线性问题时，其数学推导和编程实现可能会相对复杂。
还有一种方法是摄动法。
这是一种基于微扰理论的分析方法。
对于弱非线性问题，通过将非线性项视为对线性问题的小扰动，将问题的解表示为一个级数形式。
通过求解这个级数的各项，可以逐步逼近非线性问题的精确解。
摄动法在处理一些简单的非线性问题时非常有效，但对于强非线性问题，其精度可能会受到限制。
接下来是增量法。
在处理非线性问题时，将加载过程或变形过程分成一系列的小增量。
在每个增量步内，将问题近似为线性问题进行求解，然后逐步累加得到最终的结果。
这种方法适用于各种非线性问题，尤其是在考虑加载历史和路径相关性的情况下。
非线性有限差分法也是常用的手段之一。
它直接对控制方程进行离散，通过差分近似来表示导数项。
在处理非线性问题时，可以采用迭代的方式求解离散后的方程组。
这种方法对于简单的几何形状和边界条件的问题较为适用，但对于复杂的结构可能会面临网格划分和精度控制的挑战。
材料力学的非线性行为分析
材料力学的非线性行为分析材料力学是研究材料在外力作用下的变形和破坏行为的科学，非线性行为是指材料在受力作用时呈现出的非线性特性，即力与应变不成比例关系。
在许多工程和科学领域中，对材料力学的非线性行为进行准确和全面的分析具有重要意义。
本文将着重讨论非线性行为的基本概念、常见的非线性模型以及分析方法。
一、非线性行为的基本概念在材料力学中，强度、刚度、屈服点等参数通常被用来描述材料的特性。
然而，当外力增大到一定程度时，材料的性质将不再呈现线性关系，这时就出现了非线性行为。
非线性行为主要包括弹性-塑性行为、接触-分离行为以及材料的损伤和断裂等。
二、非线性模型的选择1. 弹塑性模型弹塑性模型是描述材料弹性和塑性变形的常用模型。
其中，最经典的是von Mises屈服准则，常用于金属的塑性变形分析。
2. 黏弹性模型黏弹性模型主要用于描述粘弹性材料的非线性行为，包括粘性和弹性两个部分。
常见的黏弹性模型有Kelvin模型和Maxwell模型。
3. 损伤模型损伤模型用于描述材料在加载过程中的损伤积累和破坏行为。
常用的损伤模型有弹塑性损伤模型、粘弹性损伤模型以及断裂力学模型等。
三、非线性行为的分析方法1. 实验测试实验测试是分析材料非线性行为最直接的方法之一。
通过应力-应变测试、拉伸试验等，可以获得材料在不同应力下的应变，进而建立非线性模型。
2. 数值计算数值计算是通过数学方法对材料力学进行模拟和计算的重要手段。
常用的数值计算方法有有限元法、边界元法、网格法等。
通过设定材料的非线性模型及边界条件，可以得到材料的应力分布和变形情况。
非线性分析的结果可用于工程设计、材料选用以及破坏预测等方面。
但是在进行非线性分析时，需要注意模型的参数选择、模型的适用性以及计算误差等因素。
总之，非线性行为是材料力学中重要的研究内容，对于理解材料的变形和破坏行为具有重要意义。
通过选择合适的非线性模型和分析方法，我们可以准确地描述和预测材料的非线性行为，为工程实践和科学研究提供有力支持。
建筑知识：建筑材料的非线性分析与优化
建筑知识：建筑材料的非线性分析与优化建筑工程的质量和稳定性是保证安全和可持续发展的重要保障，而建筑材料的质量直接关系到建筑工程的稳定性和耐用性。
在实际的建设过程中，建筑材料的非线性分析与优化是保证建筑工程质量、提高建筑材料性能的关键技术。
一、建筑材料的非线性分析建筑材料的非线性分析是指当材料承受一定的载荷时，其力学性能发生变化的现象。
材料的非线性分析是不可避免的，在设计中必须考虑到非线性效应对设计的影响，并进行相应的修正和优化。
1.轴向受压的混凝土材料的非线性分析在实际的工程应用中，混凝土出现了“骨架曲线”的特性，在不同的载荷下，它的应变硬化率也不同。
这种情况下，使用线性弹性理论来分析混凝土不能完全符合实际情况。
对于轴向受压的混凝土材料，采用理论模型可以更好地描述非线性物理现象。
通过混凝土骨架的微观分析，建立了各向同性的弹塑性理论模型，这种模型被广泛地应用于混凝土结构设计中。
2.钢筋混凝土的非线性分析在钢筋混凝土中，钢筋和混凝土输送负载的方式不同，因此在载荷作用下，这两种材料的形变和应力响应不同。
另外，在钢筋混凝土中，混凝土的应力-应变关系是非线性的，随着加荷的增加，弹性模量和抗拉强度都会增加。
对于钢筋混凝土，采用非线性有限元方法建立的数值模型可以更精确地描述其非线性特征。
该方法可以模拟出混凝土的非线性应力-应变特性和裂缝的产生和扩展情况，并根据实际材料性能进行相应的修正。
二、建筑材料的优化设计材料优化设计是保证建筑工程质量的基础工作。
优化设计的目的是在满足强度和刚度等基本要求的前提下，通过材料性能的优化实现结构的轻量化和高效化。
1.硅酸盐水泥混凝土的优化设计硅酸盐水泥混凝土作为一种新型的材料，它具有良好的力学性能和化学稳定性。
通过研究混凝土中的微纤维增强体系，可以增强混凝土的耐劈裂性和韧性，提高混凝土的力学性能。
另外，在混凝土中加入微粉、飞灰等物质，可以防止混凝土龟裂、提高混凝土的抗渗透性和耐久性。
材料非线性分析的例子
材料非线性分析的例子例子材料分析非线性光学材料非线性回归分析篇一：ANSYS结构非线性分析相应步骤及命令流ANSYS结构非线性分析相应步骤及命令流屈服准则概念：1.理想弹性材料物体发生弹性变形时，应力与应变完全成线性关系，并可假定它从弹性变形过渡到塑性变形是突然的。
2.理想塑性材料（又称全塑性材料）材料发生塑性变形时不产生硬化的材料，这种材料在进入塑性状态之后，应力不再增加，也即在中性载荷时即可连续产生塑性变形。
3.弹塑性材料在研究材料塑性变形时，需要考虑塑性变形之前的弹性变形的材料这里可分两种情况：Ⅰ.理想弹塑性材料在塑性变形时，需要考虑塑性变形之前的弹性变形，而不考虑硬化的材料，也即材料进入塑性状态后，应力不再增加可连续产生塑性变形。
Ⅱ.弹塑性硬化材料在塑性变形时，既要考虑塑性变形之前的弹性变形，又要考虑加工硬化的材料，这种材料在进入塑性状态后，如应力保持不变，则不能进一步变形。
只有在应力不断增加，也即在加载条件下才能连续产生塑性变形。
4.刚塑性材料在研究塑性变形时不考虑塑性变形之前的弹性变形。
这又可分两种情况：Ⅰ.理想刚塑性材料在研究塑性变形时，既不考虑弹性变形，又不考虑变形过程中的加工硬化的材料。
Ⅱ.刚塑性硬化材料在研究塑性变形时，不考虑塑性变形之前的弹性变形，但需要考虑变形过程中的加工硬化材料。
屈服准则的条件：1.受力物体内质点处于单向应力状态时，只要单向应力大到材料的屈服点时，则该质点开始由弹性状态进入塑性状态，即处于屈服。
2.受力物体内质点处于多向应力状态时，必须同时考虑所有的应力分量。
在一定的变形条件（变形温度、变形速度等）下，只有当各应力分量之间符合一定关系时，质点才开始进入塑性状态，这种关系称为屈服准则，也称塑性条件。
它是描述受力物体中不同应力状态下的质点进入塑性状态并使塑性变形继续进行所必须遵守的力学条件，这种力学条件一般可表示为f（σij）＝C又称为屈服函数，式中C是与材料性质有关而与应力状态无关的常数，可通过试验求得。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
(
)
(4)
6. 当最终应力状态不在屈服面上时使用下面手动回归方法将应力移动到屈服面上。
FC e aT D aC + h C e σ= σ C − δλC D aC D
δλC =
(5)
பைடு நூலகம்
屈服面的形状在各子增量的结束点使用硬化准则进行修正。 卸载时假设为弹性。
判定屈服与否的屈服条件(yield criteria) 计算塑性变形的流动法则(flow rule) 描述塑性变形时屈服面的变化的硬化法则(hardening rule)
We Analyze and Design the Future
2
midas Civil
(1)
Analysis for Civil Structure
图2.8.27 隐式后退欧拉方法
显式方法中的硬化和塑性流动的方向的计算基准位置为‘ 交叉点’，即弹性应力增量与屈 服面的交点(图2.8.25的A)。隐式方法的计算基准位置为最终应力点(图2.8.27的B)。
显式方法计算相对简单，直接对应力进行积分，即不必在高斯点(Gauss Point)重复计 算。但是显式方法有下列缺点： 8
We Analyze and Design the Future
10
midas Civil
Analysis for Civil Structure
隐式后退欧拉方法 在隐式方法中使用下面公式计算最终应力。
σ= σ B − d λ Dea C C
(6)
式中的下标的意义参见图2.8.27。
dF =
其中，h是塑性硬化系数，因此可按下面公式计算应力的变化率。
e DeaaT DeT σ d = − T D dε a Dea + h
(7)
在完全牛顿-拉普森迭代过程中使用一致性(consistent)本构矩阵时，因为其二阶收敛特 性可以让计算更快的收敛。
隐式方法不需要子增分或手动回归方法也能获得充分正确的解，并无条件稳定。但是对 于一般的屈服准则需要在高斯点反复计算。隐式方法可以构成具有关联性的切线刚度矩 阵。当使用Newton-Raphson迭代计算方法时，在高斯点进行迭代计算的效率也是较高 的。
显式前进欧拉方法 1. 首先计算应变增量
midas Civil
∂f dε p = λ ∂σ
Smooth a
midas Civil
σa
Plastic potential
g( = σ ) f( = σ ) const.
σ , dε p
σd
d Corner
dε p
图2.8.23 连续流动法则和奇异点
1
We Analyze and Design the Future
∆σ e
A X B
σ
midas Civil
X : 前一阶段的最终应力状态 A : 应力增量和屈服面的交点 B : 考虑弹性应力增量的测试应力状态
(a) 交点A的位置
∆σ e
A X D
B −∆λCe C
σ
C : 修正后的应力状态 D : 使用手动回归方法后的应力状态
(b) 在 A位置沿切线方向向 C点移动后修正到 D点 图2.8.25 显式前进欧拉方法
屈服条件 定义弹性响应(elastic response)边界的屈服函数(或加载函数)F如下(参见图2.8.23)。
F (σ , ε p ,κ = ) σ e (σ , ε p ) − κ (ε p ) ≤ 0
(2)
其中， σ : 当前的应力 σ e : 等效(equivalent)或有效(effective)应力
8-6-3 应力积分
应力积分可以使用下面两种方法。
包含子增量的显式前进欧拉方法(Explicit forward Euler algorithm with subincrementation) (图2.8.25和图2.8.26) 隐式后退欧拉方法(Implicit backward Euler algorithm) (图2.8.27)
8-6-1 塑性理论
静力塑性应变的假定如下：
本构响应(constitution response)与变形的速度无关。 弹性响应(elastic response)不受塑性变形的影响。 总应变如下：
= εe +ε p ε
其中， ε : 总应变 ε e : 弹性应变 ε p : 塑性应变 另外公式中将使用如下的基本概念：
1
We Analyze and Design the Future
第8章 | 非线性分析
B C A E
D
A, B, C, D : 各子增量中的应力修正后的应力状态 E : 使用手动回归方法后的应力状态
图2.8.26 子增量
B
C X
σx
σc
X : 前一阶段的最终应力状态 B : 考虑了弹性应力增分的测试应力状态 C : 未知的最终应力状态
−1
(
)
We Analyze and Design the Future
6
midas Civil
∂F ∂F ∂F dσ + p d ε p + d κ = aT D e d ε − aT D e a + h d λ = 0 ∂σ ∂κ ∂ε
T
(
)
(6)
Analysis for Civil Structure
因为在式(6)中不知道C点的值，所以需要通过迭代计算求解，用向量r表示当前应力与 后退欧拉应力间的差如下。
r =σ C − σ B − d λ DeaC
(
)
(7)
midas Civil
迭代计算过程就是将r减少为0，最终的应力应满足屈服条件。使用下面公式计算新的残 余量。
De a +λ rn = ro + σ
T
(10)
1
We Analyze and Design the Future
第8章 | 非线性分析
其中， ε p : 有效塑性应变
可按下式计算。 其中， λ
= λ
Fo − aT ro aT D e a + h
(11)
8-6-4 塑性材料本构模型
在程序中提供下面四种塑性材料本构模型。
∂F λ = d ε p d= d λa ∂σ
(4)
在如图2.8.23的角点或平坦的面上不能确定塑性流动的方向(不具有唯一性)，即存在奇 异点(singular point)，分析过程中程序对这些点应进行特别处理。
硬化法则(Hardening rule) 硬化法则决定材料屈服时随塑性应变变化的屈服面的变化规律。
硬化法则根据定义有效塑性应变的方法分为应变硬化(strain hardening)和 加 工 硬化 (work hardening)。应变硬化基于塑性非压缩性(plastic incompressibility)假定，适合于 不受静水压影响的材料本构。因此基于加工硬化的硬化法则的适用性更广泛。 另外如图2.8.24所示，根据屈服面的变化方式硬化法则又分为各向同性硬化(isotropic hardening)、随动硬化(kinematic hardening)、混合硬化(mixed hardening)。
We Analyze and Design the Future
4
midas Civil
此可将上面公式(3)表现为如下形式。
Analysis for Civil Structure
σ2
Translation and Expansion Initial Yield Surface
F(σ ) = κ 2
: σ 的变化量 其中， σ : d λ 的变化量 λ
(8)
可得， 将上式设为0求解 σ
Dea = σ −ro − λ
(9)
使用台劳级数表示屈服函数如下。
FCn = FC o +
∂F ∂F p = FCo + aT + σ ε Cσ − hλ = 0 ∂σ ∂ε p
F(α − α ) = κ 12 > κ 2
σ1
Translation only
σ −α) = F( κ2
图2.8.24 随动硬化和混合型硬化
midas Civil
1
We Analyze and Design the Future
第8章 | 非线性分析
8-6-2 本构矩阵(Constitutive matrix)
形成标准塑性本构矩阵的过程如下。
应力由应变的变化率向量中的弹性部分决定。即，
dσ = De d ε − d ε p = De ( d ε − d λ a )
(
)
(5)
e 其中，D 是弹性本构矩阵。
应力始终在屈服面上，所以满足下面的一致性条件(consistence condition)。
RaaT R T R − T d= σ dε a Ra + h
(8)
e ∂a 其中， R = I + d λD ∂σ
−1 e = I + d λ De A De 。 D
4. 然后计算交叉点应力。测试弹性应力增量可分为容许应力增分和非容许应力增量，
1
We Analyze and Design the Future