<?xml version="1.0" encoding="utf-8" ?><rss version="2.0"><channel><title><![CDATA[Niklas4088的博客]]></title><description><![CDATA[]]></description><link>https://blog.youkuaiyun.com/Niklas4088</link><language>zh-cn</language><generator>https://blog.youkuaiyun.com/</generator><copyright><![CDATA[Copyright &copy; Niklas4088]]></copyright><item><title><![CDATA[利用python在abaqus中画Voronoi多面体简单示例]]></title><link>https://blog.youkuaiyun.com/Niklas4088/article/details/135077952</link><guid>https://blog.youkuaiyun.com/Niklas4088/article/details/135077952</guid><author>Niklas4088</author><pubDate>Tue, 19 Dec 2023 15:05:47 +0800</pubDate><description><![CDATA[为Voronoi多面体各面上的顶点序号。或 CAE 新一个part 后删除几何。为Voronoi多面体顶点坐标。从rpy文件读取响应指令。从rpy文件读取响应指令。从rpy文件读取响应指令。]]></description><category></category></item><item><title><![CDATA[状态方程-层裂相关]]></title><link>https://blog.youkuaiyun.com/Niklas4088/article/details/120879103</link><guid>https://blog.youkuaiyun.com/Niklas4088/article/details/120879103</guid><author>Niklas4088</author><pubDate>Sun, 24 Oct 2021 10:14:33 +0800</pubDate><description><![CDATA[层裂相关的问题状态方程与体模量
状态方程与体模量
U−uU-uU−u Hugoniot 线
U=a0+suU=a_0+suU=a0​+su
Euler形式的冲击突跃条件：
[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​+]]></description><category></category></item><item><title><![CDATA[FCC单晶塑性abaqus-vumat]]></title><link>https://blog.youkuaiyun.com/Niklas4088/article/details/112713431</link><guid>https://blog.youkuaiyun.com/Niklas4088/article/details/112713431</guid><author>Niklas4088</author><pubDate>Wed, 19 May 2021 09:14:00 +0800</pubDate><description><![CDATA[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]]></description><category></category></item><item><title><![CDATA[python 画空间多边形]]></title><link>https://blog.youkuaiyun.com/Niklas4088/article/details/115339049</link><guid>https://blog.youkuaiyun.com/Niklas4088/article/details/115339049</guid><author>Niklas4088</author><pubDate>Tue, 30 Mar 2021 23:30:41 +0800</pubDate><description><![CDATA[python 画空间多边形
import mpl_toolkits.mplot3d as a3
import matplotlib.colors as colors
import matplotlib.pyplot as plt
import numpy as np
ax = 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]]></description><category></category></item><item><title><![CDATA[python个人练习]]></title><link>https://blog.youkuaiyun.com/Niklas4088/article/details/115017797</link><guid>https://blog.youkuaiyun.com/Niklas4088/article/details/115017797</guid><author>Niklas4088</author><pubDate>Wed, 24 Mar 2021 08:53:03 +0800</pubDate><description><![CDATA[python 曲线拟合（个人练习）数据读取自定义函数最小二乘法+梯度下降算法
利用最小二乘法+梯度下降算法，用自定义函数拟合离散数据。（注程序是通过jupyter notebook调试的）
数据读取
自定义函数
最小二乘法+梯度下降算法

...]]></description><category></category></item><item><title><![CDATA[Abaqus python 实例]]></title><link>https://blog.youkuaiyun.com/Niklas4088/article/details/114198784</link><guid>https://blog.youkuaiyun.com/Niklas4088/article/details/114198784</guid><author>Niklas4088</author><pubDate>Mon, 01 Mar 2021 11:23:21 +0800</pubDate><description><![CDATA[Abaqus python 实例--个人练习题目：第一步，录制宏，建立模型
学习了一段时间的Abaqus的python脚本，给自己出道题练习一下。先说一下，学习资料，主要是学习的是曹金凤老师的《Python语言在Abaqus中的应用》，如果有一定基础，需要更多实例练习可以参考江丙云老师的《ABAQUS Python二次开发攻略》。
题目：
带孔方板拉伸，逐渐加密孔边网格，孔边最大Mises应力，会相应变化，当最大Mises应力变化不大时，输出相应网格密度（孔边的种子个数）。
第一步，录制宏，建立模型
启动A]]></description><category></category></item><item><title><![CDATA[单晶塑性abaqus-umat 自学—续]]></title><link>https://blog.youkuaiyun.com/Niklas4088/article/details/112747130</link><guid>https://blog.youkuaiyun.com/Niklas4088/article/details/112747130</guid><author>Niklas4088</author><pubDate>Tue, 19 Jan 2021 12:39:38 +0800</pubDate><description><![CDATA[单晶塑性abaqus-umat 自学—续1.刚度矩阵的旋转2. LU分解3. 反向代入法
本文是针对黄老师单晶umat程序中的一些细节问题。都是本人一下没看明白，通过思考和查资料弄明白后，写在这里。
1.刚度矩阵的旋转
C-----  Rotation matrix: ROTD to transform local elastic matrix DLOCAL 
C     to global elastic matrix D
C
      DO J=1,3
         J1=1+J/3
      ]]></description><category></category></item><item><title><![CDATA[各向同性+随动硬化+过应力-vumat-理论推导]]></title><link>https://blog.youkuaiyun.com/Niklas4088/article/details/112504110</link><guid>https://blog.youkuaiyun.com/Niklas4088/article/details/112504110</guid><author>Niklas4088</author><pubDate>Thu, 14 Jan 2021 14:53:55 +0800</pubDate><description><![CDATA[过应力模型:
σ=L:εe+η ε˙p\sigma=\mathbb{L}:\varepsilon^e+\eta\,\dot\varepsilon^pσ=L:εe+ηε˙p
同时考虑各向同性硬化和随动硬化
η ε˙p=σ′−α\eta\,\dot\varepsilon^p=\sigma'-\alphaηε˙p=σ′−α

]]></description><category></category></item><item><title><![CDATA[单晶塑性abaqus-umat自学]]></title><link>https://blog.youkuaiyun.com/Niklas4088/article/details/111942686</link><guid>https://blog.youkuaiyun.com/Niklas4088/article/details/111942686</guid><author>Niklas4088</author><pubDate>Thu, 07 Jan 2021 10:51:55 +0800</pubDate><description><![CDATA[单晶塑性abaqus-umat 摘要1. 简介2. 单晶弹塑性本构模型回顾2.1运动学
本文以翻译黄永刚老师的umat文章为主，这个颜色的为个人的理解。
摘要
包含单晶塑性的abaqus用户材料子程序。本文回顾有限元形式的单晶弹塑性和黏弹塑性，包含小变形理论，以及严格的有限变形和有限旋转理论。单晶的非弹性变形由于晶体滑移，假定遵从Schmid准则。提出多种多样的分切应力与剪切变形的自硬化和潜硬化关系，并且包含在子程序内。
1. 简介
  有限元程序ABAQUS已经广泛应用在固体的变形与应力分析。除了自身大量]]></description><category></category></item><item><title><![CDATA[晶体塑形自学1—大变形]]></title><link>https://blog.youkuaiyun.com/Niklas4088/article/details/111770388</link><guid>https://blog.youkuaiyun.com/Niklas4088/article/details/111770388</guid><author>Niklas4088</author><pubDate>Mon, 28 Dec 2020 17:44:14 +0800</pubDate><description><![CDATA[晶体塑形自学1—大变形引言应变度量
引言

如上图所示，变形前后构型发生变化，当变化足够下，就可以忽略初始构型和变形后构型的区别，即小变形假设。
而大变形，更准确的说是有限变形，就是在描述应力和应变时考虑构型变化的影响。
小变形假设下，名义应力 p\bm{p}p 约等于 真应力 σ\bm{\sigma}σ
[p1100p22]≈[σ1100σ22] \left[ \begin{matrix}p_{11}&amp;0\\0&amp; p_{22}    \end{matrix}\right]\approx \]]></description><category></category></item><item><title><![CDATA[ABAQUS材料子程序学习（线性随动硬化塑性）]]></title><link>https://blog.youkuaiyun.com/Niklas4088/article/details/111288775</link><guid>https://blog.youkuaiyun.com/Niklas4088/article/details/111288775</guid><author>Niklas4088</author><pubDate>Fri, 18 Dec 2020 20:38:31 +0800</pubDate><description><![CDATA[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}]]></description><category></category></item><item><title><![CDATA[ABAQUS材料子程序学习（线性各向同性硬化塑性）]]></title><link>https://blog.youkuaiyun.com/Niklas4088/article/details/111044562</link><guid>https://blog.youkuaiyun.com/Niklas4088/article/details/111044562</guid><author>Niklas4088</author><pubDate>Wed, 16 Dec 2020 14:04:09 +0800</pubDate><description><![CDATA[ABAQUS材料子程序学习（线性硬化塑性）前言塑性力学基础umat子程序DeBug计算结果
前言
记录自己学习abaqus软件umat子程序的t过程，本文主要参考了《非线性本构关系在ABAQUS中的实现》第四章和技术邻的视频课程“非线性各向同性强化弹塑性umat子程序教程”
塑性力学基础
应变分为弹性应变和塑性应变：Δε\Delta\bf{\varepsilon }   \quadΔε
umat子程序
      SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SC]]></description><category></category></item><item><title><![CDATA[ABAQUS材料子程序学习（20年12月2日）]]></title><link>https://blog.youkuaiyun.com/Niklas4088/article/details/110497548</link><guid>https://blog.youkuaiyun.com/Niklas4088/article/details/110497548</guid><author>Niklas4088</author><pubDate>Wed, 02 Dec 2020 20:24:48 +0800</pubDate><description><![CDATA[@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,]]></description><category></category></item><item><title><![CDATA[ABAQUS材料子程序学习（20年12月1日）]]></title><link>https://blog.youkuaiyun.com/Niklas4088/article/details/110452489</link><guid>https://blog.youkuaiyun.com/Niklas4088/article/details/110452489</guid><author>Niklas4088</author><pubDate>Tue, 01 Dec 2020 21:21:34 +0800</pubDate><description><![CDATA[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,
   ]]></description><category></category></item></channel></rss>