Niklas4088-CSDN博客

原创利用python在abaqus中画Voronoi多面体简单示例

为Voronoi多面体各面上的顶点序号。或 CAE 新一个part 后删除几何。为Voronoi多面体顶点坐标。从rpy文件读取响应指令。从rpy文件读取响应指令。从rpy文件读取响应指令。

2023-12-19 15:05:47 365

原创状态方程-层裂相关

层裂相关的问题状态方程与体模量状态方程与体模量U−uU-uU−u Hugoniot 线U=a0+suU=a_0+suU=a0+suEuler形式的冲击突跃条件：[u]=−ρ+U[V][p]=ρ+Uu−[u]=-\rho^+U[V]\\[p]=\rho^+Uu^-[u]=−ρ+U[V][p]=ρ+Uu−代入：[p]=2ρ+2U2[V][p]=2{\rho^+}^2U^2[V][p]=2ρ+2U2[V]u=−ρ+(a0+su)[V]u=-\rho^+(a_0+su)[V]u=−ρ+(a0+

2021-10-24 10:14:33 366

原创 FCC单晶塑性abaqus-vumat

FCC单晶塑性abaqus-vumat 1.整体思路1.1 弹性各向同性以黄永刚老师的umat文章为基础，构建单晶塑性的vumat。本文是之前“单晶塑性abaqus-umat 自学”的后续，大量引用前文已推导出的结论，用（umat）代指上文。求解过程针对FCC。1.整体思路先从最简单的入手。1.1 弹性各向同性弹性采用各向同性，且分析过程中不考虑体积变化，动态分析中可根据状态方程(EOS)单独给出P-V关系。体变： ΔσV=KΔεv(1.1.1a)\Delta\sigma_V=K\Delta\v

2021-05-19 09:14:00 4621 16

原创 python 画空间多边形

python 画空间多边形import mpl_toolkits.mplot3d as a3import matplotlib.colors as colorsimport matplotlib.pyplot as pltimport numpy as npax = a3.Axes3D(plt.figure())A=[[[0,0,0],[1,1,0],[1,0,1]],[[0,0,0],[1,1,0],[0.5,1,0.5],[0,0.5,0.5],]]for p0 in A: tr

2021-03-30 23:30:41 556

原创 python个人练习

python 曲线拟合（个人练习）数据读取自定义函数最小二乘法+梯度下降算法利用最小二乘法+梯度下降算法，用自定义函数拟合离散数据。（注程序是通过jupyter notebook调试的）数据读取自定义函数最小二乘法+梯度下降算法...

2021-03-24 08:53:03 233

原创 Abaqus python 实例

Abaqus python 实例--个人练习题目：第一步，录制宏，建立模型学习了一段时间的Abaqus的python脚本，给自己出道题练习一下。先说一下，学习资料，主要是学习的是曹金凤老师的《Python语言在Abaqus中的应用》，如果有一定基础，需要更多实例练习可以参考江丙云老师的《ABAQUS Python二次开发攻略》。题目：带孔方板拉伸，逐渐加密孔边网格，孔边最大Mises应力，会相应变化，当最大Mises应力变化不大时，输出相应网格密度（孔边的种子个数）。第一步，录制宏，建立模型启动A

2021-03-01 11:23:21 4988 15

原创单晶塑性abaqus-umat 自学—续

单晶塑性abaqus-umat 自学—续1.刚度矩阵的旋转2. LU分解3. 反向代入法本文是针对黄老师单晶umat程序中的一些细节问题。都是本人一下没看明白，通过思考和查资料弄明白后，写在这里。1.刚度矩阵的旋转C----- Rotation matrix: ROTD to transform local elastic matrix DLOCAL C to global elastic matrix DC DO J=1,3 J1=1+J/3

2021-01-19 12:39:38 3734 21

原创各向同性+随动硬化+过应力-vumat-理论推导

过应力模型:σ=L:εe+η ε˙p\sigma=\mathbb{L}:\varepsilon^e+\eta\,\dot\varepsilon^pσ=L:εe+ηε˙p同时考虑各向同性硬化和随动硬化η ε˙p=σ′−α\eta\,\dot\varepsilon^p=\sigma'-\alphaηε˙p=σ′−α

2021-01-14 14:53:55 3747

原创单晶塑性abaqus-umat自学

单晶塑性abaqus-umat 摘要1. 简介2. 单晶弹塑性本构模型回顾2.1运动学本文以翻译黄永刚老师的umat文章为主，这个颜色的为个人的理解。摘要包含单晶塑性的abaqus用户材料子程序。本文回顾有限元形式的单晶弹塑性和黏弹塑性，包含小变形理论，以及严格的有限变形和有限旋转理论。单晶的非弹性变形由于晶体滑移，假定遵从Schmid准则。提出多种多样的分切应力与剪切变形的自硬化和潜硬化关系，并且包含在子程序内。1. 简介有限元程序ABAQUS已经广泛应用在固体的变形与应力分析。除了自身大量

2021-01-07 10:51:55 9948 29

原创晶体塑形自学1—大变形

晶体塑形自学1—大变形引言应变度量引言如上图所示，变形前后构型发生变化，当变化足够下，就可以忽略初始构型和变形后构型的区别，即小变形假设。而大变形，更准确的说是有限变形，就是在描述应力和应变时考虑构型变化的影响。小变形假设下，名义应力 p\bm{p}p 约等于真应力 σ\bm{\sigma}σ[p1100p22]≈[σ1100σ22] \left[ \begin{matrix}p_{11}&0\\0& p_{22} \end{matrix}\right]\approx \

2020-12-28 17:44:14 2066 1

原创 ABAQUS材料子程序学习（线性随动硬化塑性）

ABAQUS材料子程序学习（线性随动硬化塑性）前言塑性力学增量形式实现前言继续塑性本构，随动硬化塑性。记录自己学习abaqus软件umat子程序的t过程，本文主要参考了《非线性本构关系在ABAQUS中的实现》第四章和技术邻的视频课程“非线性各向同性强化弹塑性umat子程序教程”塑性力学增量形式实现计算过程中，体应变和体应力是弹性关系Δσv=K⋅Δεv(1){Δσ_v}={K}\cdot{Δε_v}\tag{1}Δσv=K⋅Δεv(1)KKK为体积模量，K=E3(1−ν)K=\cfrac{E}

2020-12-18 20:38:31 7297 9

原创 ABAQUS材料子程序学习（线性各向同性硬化塑性）

ABAQUS材料子程序学习（线性硬化塑性）前言塑性力学基础umat子程序DeBug计算结果前言记录自己学习abaqus软件umat子程序的t过程，本文主要参考了《非线性本构关系在ABAQUS中的实现》第四章和技术邻的视频课程“非线性各向同性强化弹塑性umat子程序教程”塑性力学基础应变分为弹性应变和塑性应变：Δε\Delta\bf{\varepsilon } \quadΔεumat子程序 SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SC

2020-12-16 14:04:09 8358 17

原创 ABAQUS材料子程序学习（20年12月2日）

@ABAQUS材料子程序学习（20年12月2日）前言继续记录自己学习过程，本文针对《非线性本构关系在ABAQUS中的实现》第三章“黏弹性“本构的学习，UMAT子程序，有VS code编写。与abaqus帮助文档里的umat例子相似二、umat代码如下（示例）： SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD, 1 RPL,DDSDDT,DRPLDE,DRPLDT, 2 STRAN,DSTRAN,TIME,DTIME,

2020-12-02 20:24:48 2151 2

原创 ABAQUS材料子程序学习（20年12月1日）

ABAQUS材料子程序学习（20年12月1日）前言umat子程序DeBug计算结果前言记录自己学习过程，本文针对《非线性本构关系在ABAQUS中的实现》第二章“动态回弹模量“本构的学习，书中自带程序可能由于转成PDF的原因有一些排版上的混乱，下文的UMAT子程序，有VS code编写，和书中例子基本相同。umat子程序 SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD, 1 RPL,DDSDDT,DRPLDE,DRPLDT,

2020-12-01 21:21:34 5081 3

Niklas4088的博客