双线性弹塑性模型(三)

本节用Python来实现基于随动硬化模型的当前应力计算。

[算例] 一根各向同性杆,一端固定,另一端施加轴向力做拉伸试验,荷载分级来加。某一时刻应力 ,塑性应变 , .当应变增量 ,计算应力和塑性应变。

import math

#定义符号函数
def Sgn(x):
    if x > 0 :
        return 1
    elif x < 0 :
        return -1
    else :
        return 0

def KinematicHard1D(MP,deps, stressN,alphaN, epsN):

    E  = MP[0#弹性模量
    H  = MP[1#塑性模量
    Y0 = MP[2#初始屈服应力
    
    stress_tr = stressN + E*deps
    eta_tr = stress_tr - alphaN
    f_tr = math.fabs(eta_tr ) - Y0
    if f_tr < 0:
        stress = stress_tr
        alpha = alphaN  #alpha不变
        ep = epsN       #塑性应变不增加,为0
        flag = 0     # 处于弹性状态的标志
    else:
        dep = f_tr / ( E + H )
        stress = stress_tr - Sgn(eta_tr) * E * dep
        alpha = alphaN + Sgn(eta_tr) * H * dep
        ep = epsN + dep
        flag = 1    # 处于塑性状态的标志

    return stress, alpha, ep, flag

E = 200E3; H = 25E3; sYield = 250
Et = E*H / (E+H)
mp = [E, H, sYield]
#上一步的应力,塑性应变以及alpha值
nS = 200
nep = 1E-4
nA = 2.5

delta_eps = -0.003

Snew, Anew, epnew, flag = KinematicHard1D(mp,delta_eps,nS,nA,nep)

print(Snew, Anew, epnew, flag)

手算结果在前一篇

双线性弹塑性模型(二)


扫码_搜索联合传播样式-白色版.png
PYTHON非线性材料

双线性弹塑性模型(三)的评论0条

    暂无评论

    双线性弹塑性模型(三)的相关案例教程

    本节内容为多杆结构的弹塑性有限元计算。 对于弹塑性材料, ,其中 含多个杆单元的结构,需要分别判断每个单元的弹塑性状态,确定是 或者 参与计算。 [算例] 如图所示,两个并联的杆,一段固定,另一端另一端施加轴向力P。荷载逐渐增加到 ,然后逐渐卸载至0。杆一。 杆二。 import matplotlib.pyplot as plt import numpy as np import numpy as
    1 本构理论 1.1 率形式 本构方程为: 单轴拉伸的应力应变的硬化曲线如下: 根据单轴试验得到硬化部分的曲线: 当仅考虑随动硬化时,屈服面的中心在移动,而屈服面的大小不发生改变,即为常数。 屈服条件为: 增加了背应力来表示屈服面中心移动即随动硬化的效果。式中相对应力的表达式为: 流动法则为: 背应力的演化法则为: 式中: 上标k代表随动部分(kinematic),表示以下随动硬化曲线的梯度: 1
    四阶张量双点乘: subroutine aaaa_dot_dot_bbbb(n,a,b,c) include 'ABA_PARAM.INC' dimension a(n,n,n,n), b(n,n,n,n), c(n,n,n,n) do i = 1,n do j = 1,n do k = 1,n do l = 1,n c(i,j,k,l) = 0 do m1 = 1,n do m2 = 1,n c
    一、ABAQUS自带Tresca本构与VUMAT对比 二、Tresca本构模型介绍 以下, 粗体符号表示向量或矩阵,上标“T”表示向量转置。 当屈服函数f(σ)的值为零时,材料屈服。这里σ是应力张量(为列矩阵)。如采用相关联的流动法则,则无穷小的塑性应变增量为 式中,a是屈服函数的梯度,dλ是一个非负塑性乘子(non-negative plastic multiplier)。 如果在分析增量步j
    摘要:串口是通信中最常用的通信方式,可能写串口的驱动,能写几十种方法, 查询方式,中断方式,DMA方式,定时器方式。可能也其中几种方式的组合形式,经典的用法是:发送用查询方式,接收用中断方式,或者DMA+空闲中断。本篇不讲串口是啥,现在还在讲串口是啥,估计会被喷。今天来聊一聊串口常用的几种方式,最简单的方法就不说了。 一、经典方法 查询方式 可靠性很高,要考虑下个数据包覆盖上一个数据包的问题,小数
    本科/工程师
    影响力
    粉丝
    内容
    获赞
    收藏
      项目客服
      培训客服
      0 3